发明名称 一种基于模态信息的结构损伤诊断新方法
摘要 本发明一种基于大型在役结构物动力特性的损伤检测方法利用损伤前后的模态信息的差异进行损伤定位和损伤程度评估,利用结构的单元刚度矩阵和质量矩阵作为诊断项,根据损伤表征向量一次性进行损伤定位和损伤程度评估。本发明不需要质量归一化的振型,仅利用结构损伤前后的低阶模态信息即可,可以同时考虑结构刚度和质量的变化。该方法适用于单损伤,多种损伤工况下,能够精确定位出框架结构的损伤位置和损伤程度,具有一定的实际应用价值。
申请公布号 CN101034053A 申请公布日期 2007.09.12
申请号 CN200710013349.9 申请日期 2007.02.13
申请人 中国海洋大学 发明人 李华军;胡少龙;王树青;王俊荣;张敏
分类号 G01N19/08(2006.01);G01M7/00(2006.01);G01M19/00(2006.01) 主分类号 G01N19/08(2006.01)
代理机构 青岛联智专利商标事务所有限公司 代理人 崔滨生
主权项 1、一种基于模态信息的结构损伤诊断新方法,其特征在于:包括如下步骤:(1)基准模型建立:首先,根据结构物未损伤状态下的结构建立有限元模型;其次,将上述结构物的在役结构动力状态下,存储结构物服役阶段的环境激励下传感器实测的结构动力响应数据,利用模态参数识别技术识别其模态参数;最后通过模型修正方法修正所述在役结构物的有限元模型,得到能够反映未损伤结构实际动力特性的模型,作为损伤诊断的基准模型;(2)损伤检测:利用传感器测量结构物损伤后的结构动力响应数据,利用模态参数识别技术识别其模态参数,从实验模型中提取出Nj阶模态,从基准模型中提取出Ni阶模态,通过基准模型和实验模型的交叉,以及不同阶的模态交叉形成Nm=Ni×Nj个交叉模型交叉模态方程,写为矩阵形式,有:Cα+Eβ=f式中C和E为Nm×Ne阶矩阵;α和β为Ne阶列向量;f为Nm阶列向量,由此式可写为:Gγ=f式中:G=[C E],<math> <mrow> <mi>&gamma;</mi> <mo>=</mo> <mfenced open='{' close='}'> <mtable> <mtr> <mtd> <mi>&alpha;</mi> </mtd> </mtr> <mtr> <mtd> <mi>&beta;</mi> </mtd> </mtr> </mtable> </mfenced> </mrow> </math> 通过最小二乘法解出<math> <mrow> <mover> <mi>&gamma;</mi> <mo>^</mo> </mover> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>G</mi> <mi>T</mi> </msup> <mi>G</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>G</mi> <mi>T</mi> </msup> <mi>f</mi> <mo>,</mo> </mrow> </math> 上标T表示转置,其中<math> <mrow> <msub> <mi>C</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>ij</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>&Phi;</mi> <mi>i</mi> <mi>T</mi> </msubsup> <msub> <mi>K</mi> <mi>n</mi> </msub> <msubsup> <mi>&Phi;</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mo>,</mo> </mrow> </math> <math> <mrow> <msub> <mi>E</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msubsup> <mi>&lambda;</mi> <mi>j</mi> <mo>*</mo> </msubsup> <msub> <mi>D</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <mo>,</mo> </mrow> </math> <math> <mrow> <msub> <mi>D</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>m</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>&Phi;</mi> <mi>i</mi> <mi>T</mi> </msubsup> <msub> <mi>M</mi> <mi>n</mi> </msub> <msubsup> <mi>&Phi;</mi> <mi>j</mi> <mo>*</mo> </msubsup> <mo>,</mo> </mrow> </math> <math> <mrow> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>C</mi> <mi>m</mi> </msub> <mo>+</mo> <msubsup> <mi>&lambda;</mi> <mi>j</mi> <mo>*</mo> </msubsup> <msub> <mi>D</mi> <mi>m</mi> </msub> <mo>,</mo> </mrow> </math> K和M分别表示结构的刚度矩阵和质量矩阵,Kn表示第n个单元的单元刚度矩阵,Mn表示第n个单元的单元质量矩阵,Фi,Фj*分别指结构基准模型和实验模型的模态振型,λi和λj*是基准模型和实验模型的特征值,这样利用基准模型和实测模态信息即可确定表示结构单元刚度和质量发生变化的修正系数αi(i=1,…,Ne)和βi(i=1,…,Ne),其中α为结构单元的刚度修正系数,若-1≤αi<0则表示第i单元发生损伤;β表示结构单元质量的变化,若βi≠0表示第i单元的质量发生了变化;(3)损伤程度评估:所得到的修正系数αi(i=1,…,Ne)和βi(i=1,…,Ne)的值的大小代表了相应单元的损伤程度和质量变化情况。
地址 266100山东省青岛市崂山区松岭路238号