基于响应面法的结构损伤识别研究.doc

基于响应面法的结构损伤识别研究摘要本文以通过建立悬臂梁的有限元模型获得的数据拟合出输入输出之间关系的响应面模型，来证明响应面法可以有效快速的进行结构损伤识别研究。1. 综合介绍近年来用于结构损伤检测方面的方法，说明各方法的优劣。随后详细的讲解本文要介绍的响应面法的基本内容和设计方法，说明响应面法在结构损伤识别中的有效性和独特优越性。2. 建立一个悬臂梁的有限元模型，简单介绍ANSYS在有限元问题中运用的方便性，并简述ANSYS所采用的分析方法，之后用ANSYS进行有限元分析，得到悬臂梁不同定位和不同刚度下的固有频率值。为拟合响应面做准备。3. 运用DESIGN EXPERT软件，得到输入与输出的二次响应面方程。通过该方程拟合出三维响应面模型，最后进行实际值与预测值的对比，证明响应面法的正确性和可操作性。4. 总结工作成果。关键词 响应面法 结构损伤识别中图分类号：T 文献标识码：A 文章编号：1001-2400(2XXX)0X-0-0Research on Stuructre Damage Identification Based on Response Surface MethodCaokai(School of Electronic Engineering, Xidian Univ., Xian 710071, China)Abstract: In this paper，through establishment of a finite element model for a cantilever beam to obtain the data and use the data to create a response surface between the input and the output,in order to prove the response surface method can efficiently and quickly make the research of structure damage identification.1. Introduce the methods of structure damage identification in recent years, explain pros and cons of the various of the methods.Then make a detail explanation of the basic content and design method of the response surface method,to describe the response surface method is efficient and have unique advantages in structure damage identification.2. Create a finite element model for a cantilever beam, make a brief introduction of ANSYS, to display the convenience of ANSYS in the finite element model problem. Then use ANSYS to get the natural frequency values of the cantilever beam in the different location and stiffness. To prepare fitting the response surface model.3. Through the design expert software, we can obtain the quadratic polynomial of the response surface between the input and the output. We can create the 3D model of the response surface through the quadratic polynomial, make the comparison between the actual and the predict in the end to prove that the response surface method is correct and operable.4. Summing up the result of the work.Key Words: Response surface method Structure damage identification0引言近20年来，随着现代传感技术、无线通讯技术、微处理技术、信号采集与处理、信息融合以及系统建模等技术的日益精确和完善，结构损伤动力检测技术也迅速发展起来，从而在航空航天，土木工程，机械工程等学科都有了较大的发展和应用，为结构的安全情况进行评估和预测。针对于大型桥梁等工程测量大部分方法费时费力，所以提出一个有效的快速的方法使人们所期望的。而本文介绍的响应面法就是这样一个方法。响应面法通过建立一个输入/输出关系的响应面模型，从而实现对结构损伤部位和程度的预测。而本文采取的检测参数为固有频率，通过固有频率的测定，进行反向预测刚度和定位的情况，来确定响应面法准确性。因为不需要有限元的计算，所以响应面的运行是高效和快速的1响应面法的基本理论假定参数或设计点是n维向量，它是待求性能函数的自变量，二者藏在函数关系为。尽管未知的函数可能找不出准确的表达式，但是只要给定了参数值或设计点值，即去定了一个样本点，总可以通过实体的或数值的实验得到相应的性能值，这是对应一个参数值或设计点值的响应值。如果做了足够多的实验，例如m个实验，那么我们就可以利用m个样本点及其m个响应。利用待定系数的方法求出函数的近似函数 (2-1)式中：，待构造的响应面函数。由于性能响应与变量之间的函数关系是未知的，因此事先必须选择函数的形式。选择好的函数会使近似更精确，而且会使适合使用的设计空间更宽广。响应面函数形式的选取要符合以下两个要求：第一个，响应面函数数学表达式应该在能够描述真实函数的前提下尽可能简单。第二个，应在响应面函数中设计尽可能少的待定系数以减少实体实验或数值分析的工作量。因此，实际运用时，通常选取线性或二次多项式的形式。线性或二次多项式表达如下：线性型 (2-2)不含交叉项的二次型 (2-3)含交叉项的二次型 (2-4)式中 常数项待定系数 一次项待定系数 二次想待定系数为了推导统一和简便，令 (2-5) (2-6)将(1-5)(1-6)带入(1-4)中，得到统一的形式 (2-7)式中：待定系数。的个数k根据响应面函数的形式确定，见表2-1。表2-1 函数形式与待定系数个数函数形式待定系数个数k线性型n+1可分离二次型（不含交叉项）2n+1完整二次型(n+1)(n+2)/2为了确定系数，需要做m次(mk)独立试验，每次试验中个变量的趋之不同，得到m个样本点对应的响应值(i=0，m-1)，根据式(1-5)进行换算，得到的数据为 (2-8)将上述m个样本点(j=0，1，m-1)带入式中，得到响应面函数值为 (2-9)因为响应面函数是性能函数的近似函数，所以式(2-9)计算出的结果通常不等于试验都得出的响应值，也就是说，二者存在一个误差，即 (2-10)目前，式(1-10)中的 (i=0，1,，m-1)尚未确定，为此，我们可以通过使极小化的途径，使误差最小，同时也合理的确定了系数(i=0，k-1)，换句话说，为找到最接近所有试验数据点的响应面，利用最小二乘原理使误差的平方和最小，即 (2-11)式(2-11)取极小值的必要条件为 (i=0,，k-1) (2-12)这是k个方程k个未知数的线性方程组，化简并整理，可得 (2-13)写成矩阵形式为 (2-14)其中 若矩阵奇异，则需进行奇异值分解，或采用松弛法，或设计变量归一法，若不奇异，则 (2-15)将式(2-15)得到的代入式(1-7)即得到响应面函数的表达式。以上推导过程也可以采取矩阵形式，令误差平方和去最小，即只需，可求得2 响应面建模方法的统计评价指标上述求得响应面模型能否作为有意义的近似模型，还应对响应面的预测能力进行评估。复相关系数 (2-18)式中：响应值与响应均值差的平方和； 响应值与响应估计值差的平方和； 响应估计值与响应均值差的平方和，其表达式分别为 (2-19) (2-20) (2-21)是一个在之间变化的值，它的值越接近1说明误差的影响越小，即回归方程越准确。若等于1，说明回归方程可以精确的描述y的变化，即观测点全部落在回归方程所确定的曲面上。可以描述响应面的拟合程度，但它有一个缺陷，即其值随回归方程中自变量个数的增加而增加。当所有自变量均在回归方程中时，达到最大，因此不能认为越大回归方程的逼近程度就越好，因为冗余参数的存在也会提高的值。尽管如此，可用于比较具有相同参数个数的不同回归方程的逼近程度，此时越大说明回归效果越好。3模态分析本文要进行ANSYS里面的悬臂梁建模，之后进行模态分析，模态分析一般用于确定结构的振动特性，即结构的固有频率和振型（模态）。 典型的无阻尼模态分析求解的基本方程是经典的特征值问题如下所示： (3-1)式中刚度矩阵 第i阶模态的振型矢量（特征矢量） 第i阶模态的固有频率（特征值） 质量矩阵本文采取的模态提取方法是Block Lanczos（分块兰索斯）法。建模图形如图3.1所示。表3.1不同弹性模量下的五阶固有频率弹性模量(Pa)定位(Hz)(Hz)(Hz)(Hz)(Hz)（Hz）1e1110.278842.67817.354414.15724.44137.0892e1110.373572.88788.209715.75126.21239.9193e1110.435543.02408.606816.70127.52741.4794e1110.480743.13488.858617.30828.53442.7245e1110.515593.23119.047217.72929.30643.7781e1120.352102.84617.460514.73923.81537.8712e1120.431423.04248.012116.27726.29740.4953e1120.472623.13578.435716.97127.75541.8634e1120.498153.19258.772917.04828.67842.9025e1120.515593.23119.047217.72929.30643.7781e1130.435992.22117.905614.92923.00837.8742e1130.481382.63578.577815.86526.09740.5703e1130.499612.90268.832216.59527.71141.9234e1130.509443.09088.965317.20628.67442.9285e1130.515593.23119.047217.72929.30643.7781e1140.4979823.2896.735314.15423.64037.8642e1140.508882.78587.562715.98626.18740.5063e1140.512593.01148.177616.81227.69941.8794e1140.514463.14428.660017.33728.65742.9115e1140.515593.23119.047217.72929.30643.7781e1150.514893.11937.679914.16024.37237.0202e1150.515333.18988.499815.84026.17839.8683e1150.515473.21298.802616.79327.52441.4464e1150.515543.22438.955617.36228.54242.7105e1150.515593.23119.047217.72929.30643.778最终我们可以根据表3.1求得编码后的响应面曲面方程如下： A代表弹性模量 B代表定位下面是对弹性模量和定位的预测，图中斜直线为理想值，可以看到图4.7中，弹性模量很好的识别，而图4.8对定位的预测则稍显离散。图4.7为响应面对弹性模量的预测图4.8为响应面对定位的预测4 综述本文以基于响应面法的结构损伤为课题，进行了深入的分析和研究。以响应面法为核心，介绍了响应面法的不同试验