基于摄动理论的结构损伤识别优化算法.docxVIP

基于摄动理论的结构损伤识别优化算法

在土木工程、航空等领域，大型结构健康监测技术是国内外科技、工程专业领域高度重视和研究的主题。结构健康监测的关键和关键是识别缺陷。根据结构动态特征的变化，识别整体或局部结构损伤是简单快速的，也是识别结构损伤的主要趋势。

结构损伤识别方法可以分为有模型的识别方法和无模型的识别方法。目前国内外对有模型的识别方法开展研究较多,有模型的识别方法使用与模型有关的特征量识别损伤,常用的特征量有固有频率、模态振型及其衍生量等,由于测试模态集不完备、测试自由度不足和测量噪声等影响,所测模态振型的误差较大,导致利用模态振型及其衍生量的损伤识别结果不可靠,而目前固有频率的测试精度却很高,低阻尼结构频率的测试分辨率可以达到0.1%,甚至更高,因此,基于频率变化的结构损伤识别方法可以提高损伤识别的精度。

Cawley和Adams最早发现结构两阶不同频率的变化比仅仅是损伤位置的函数,并依此识别了板的单个损伤位置。1990年Stubbs和Osegueda建立了频率变化量、广义刚度对单元刚度参数的灵敏度矩阵和损伤定位参数之间的线性方程组,由于所得的方程组高度亚定,迭代过程需要反复计算灵敏度矩阵,计算量很大,而且迭代没有终止的可靠判据,很难得出准确的结果。1993年Hasiotis提出用优化方法求解欠定方程组,将求解变成了带有等式约束和不等式约束的二次规划问题,但所引入的权重因子受人为因素的影响较大。2000年FabrizioVestroni和DaniloCapecchi给出了利用结构固有频率变化量识别损伤结果的唯一性证明。

本文根据工程中常见的结构损伤情况在有限元模型中定义了结构单元的损伤识别参数,将摄动理论与振动理论相结合导出结构振动特征值的一、二阶摄动方程以及结构的一、二次损伤识别方程,由此建立了损伤识别参数、模态参数和物理参数之间的函数关系,并给出了损伤识别方程的优化解法。该方法仅利用结构损伤前后若干低阶频率测量值便可确定结构多处损伤位置和损伤程度以及结构老化程度,而且避免了方程求解时需要反复计算灵敏度矩阵与人为因素影响的缺陷。

1等效差分法的基本原则

1.1损伤识别的区别

结构损伤识别的目标主要是定位和定量,因此,首先应定义一个或多个能够给出损伤程度与位置的量化参数,即损伤识别参数。

总结结构常见的各类损伤及其特点如表1所示。

由表1可以看出,对于“可见性”大的结构损伤,一般无需进行识别,可以直接观察出损伤的存在,而对于可见性小的损伤,它们一般有两个共同的特点:(1)这类损伤一般对单元的质量没有影响或影响很小。(2)对单元刚度的影响常是负的,即损伤使单元刚度降低,而且损伤的大小一般与单元刚度的减少成正变关系。

根据上述分析,可以把单元未损伤的刚度阵与损伤后刚度阵表示成如下关系

式中,Kej0表示未损伤单元刚度阵;Kej表示损伤后单元刚度阵;εj是一个小参数,表示单元刚度损失量与原刚度的比值,εj即是定义的结构损伤识别参数,其大小表示损伤程度,脚标j表示损伤位置,即损伤所在的单元。式(1)实际上把损伤单元的刚度阵视为未损伤单元刚度阵的摄动,它建立了结构单元损伤与结构参数之间(某种平均意义)的量化关系。

1.2再识别单元的确定

对于实时监测的大型结构,一次识别中一般只会遇到少量可能损伤位置(偶然事件除外),而且损伤位置常集中在应力较大的区域,因此,可先根据结构的实际情况或经验判断一部分可能损伤区域,给损伤区域的每个单元刚度矩阵一个待定的损伤识别参数εj,而可能损伤区域以外的损伤识别参数直接置零,从而减少识别未知量。为防止漏判,可以再对可能损伤区域以外的单元重新进行一次识别,因为使用的输入量是个绝对量,两次或多次识别不会相互干扰。

若结构共有m个单元可能损伤,即m个单元刚度发生变化,则损伤后结构的总体刚度阵为

式中,K1j是j单元刚度矩阵关于参数εj的一阶变化率,单元形式确定后其为已知量;∑表示按有限单元法组集。

根据摄动有限元方法,可将系统的振动特征值和特征向量按小参数ε展开成Taylor级数,忽略3阶以上小量,可写为

式中,s表示系统特征值和特征向量的阶数,λs1j、us1j、λs2ij、us2ij分别为特征值和特征向量的一阶变化率和二阶变化率,且

由广义特征值问题知

将式(2)~式(4)代入式(5),由等式两边ε同次幂系数相等原则,可得如下递推求解方程组

由此可求出

再根据展开定理可得

同理可得

上述公式给出了结构损伤识别参数、结构物理参数与结构模态参数三者之间的量化关系的正向分析过程。

2损坏识别参数的求解

2.1未损伤情况下的有限元模型

结构损伤识别过程实际上是第1节所述由结构损伤所导致的结构动态特性变化分析过程的逆过程,为利用式(3)求解结构离散模型的损伤识别参数ε,假设未损伤结构的有限元模