发明名称 大型结构物无损伤的健康检测方法
摘要 本发明是大型结构物无损伤的健康检测方法,其步骤有:存储步骤;识别定位步骤;损伤程度评估步骤;比较确定步骤。根据这个方法可有效地监测出结构物损伤位置并验证诊伤结果的正确性。应用本检测方法中的结构拉压模态应变能诊伤指标(AMSECR);和/或结构弯曲模态应变能诊伤指标(TMSECR),计算诊伤指标的数值:从而诊断确定结构损伤位置。在诊断确定损伤程度的过程中,根据计算的诊伤指标数值,依据SD<sub>j</sub>≥-1,确定结构损伤程度;进一步确定结构物安全/或不安全。本发明提出适用于大型结构实时监测的方法,用于建立结构物健康监测系统,提高结构物的疲劳寿命和极大环境载荷下的生存能力,以达到提高结构物的安全性、经济性以及使作业人员有安全感。
申请公布号 CN100575924C 申请公布日期 2009.12.30
申请号 CN200510044576.9 申请日期 2005.09.18
申请人 中国海洋大学 发明人 李华军;杨和振;王树青;杨永春;石湘
分类号 G01N19/08(2006.01)I;G06F17/00(2006.01)I 主分类号 G01N19/08(2006.01)I
代理机构 青岛海昊知识产权事务所有限公司 代理人 韩振东
主权项 1、一种基于大型在役结构物动力特性的无损检测方法,其特征在于:所述的方法包括如下步骤:(1)存储步骤:首先,将上述结构物未损伤状态下即在结构设计阶段的环境激励下的结构动力响应数据存储入专用存储器中;其次,将上述结构物的在役结构动力状态下即结构服役阶段的环境激励下的传感器实测的结构动力响应数据存储入专用存储器中;(2)识别定位步骤:利用自然激励技术识别出上述结构物损伤前、后的模态参数,确定结构构件损伤出现的位置;即首先,将模态参数扩阶;其次,应用有限元模型将传感器实测的自由度凝聚;再次,进行结构环境激励下的模态参数识别:应用模态应变能分解法诊断上述结构物的具体损伤位置;(3)损伤程度评估步骤:利用结构拉压模态应变能变化指标和结构弯曲模态应变能变化指标,分析评估上述结构物损伤构件的损伤程度;(4)比较确定步骤:根据上述步骤(3),确定上述结构物构件损伤程度的评估指标数值在:-1≤SD<sub>j</sub>≤0范围内,式中SD<sub>j</sub>表示结构构件的损伤程度;所述的(2)识别定位步骤中应用模态应变能分解法诊断结构物的损伤位置具体为:1)、使用结构拉压模态应变能诊伤指标:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mi>j</mi><mi>a</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>&beta;</mi><mi>j</mi><mi>a</mi></msubsup><mo>-</mo><mover><msup><mi>&beta;</mi><mi>a</mi></msup><mo>&OverBar;</mo></mover></mrow><msub><mi>&sigma;</mi><msup><mi>&beta;</mi><mi>a</mi></msup></msub></mfrac></mrow></math>]]></maths>①式中结构拉压模态应变能变化为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&beta;</mi><mi>j</mi><mi>a</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>m</mi></msub></munderover><mfrac><mrow><mrow><mo>(</mo><msup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mi>T</mi></msup><msubsup><mi>K</mi><mrow><mi>j</mi><mn>0</mn></mrow><mi>a</mi></msubsup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><msup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>a</mi></msubsup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow><msup><msub><mi>&Phi;</mi><mi>i</mi></msub><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>a</mi></msubsup><msub><mi>&Phi;</mi><mi>i</mi></msub></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&Phi;</mi><mi>i</mi></msub><mi>T</mi></msup><msubsup><mi>K</mi><mrow><mi>j</mi><mn>0</mn></mrow><mi>a</mi></msubsup><msub><mi>&Phi;</mi><mi>i</mi></msub><mo>+</mo><msup><msub><mi>&Phi;</mi><mi>i</mi></msub><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>a</mi></msubsup><msub><mi>&Phi;</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>a</mi></msubsup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac></mrow></math>]]></maths>②式中Φ<sub>i</sub>,Φ<sub>i</sub><sup>*</sup>分别指结构损伤前和损伤后的模态振型,其中上标″T″指转置,K<sub>j0</sub>=K<sub>j</sub>/E<sub>i</sub>,其中K<sub>j</sub>是j单元未损伤系统刚度矩阵,K<sub>j0</sub><sup>a</sup>是j单元的拉压刚度矩阵,K<sub>0</sub><sup>a</sup>是通过单元矩阵K<sub>j0</sub><sup>a</sup>组装起来的拉压系统矩阵;β<sub>j</sub><sup>a</sup>还表示<img file="C2005100445760002C3.GIF" wi="150" he="75" />之比,E<sub>j</sub><sup>a</sup>,<img file="C2005100445760002C4.GIF" wi="64" he="72" />各表示结构j单元未损伤和损伤后的结构构件弹性模量,来估计上述结构物的构件损伤位置;和2)、使用结构弯曲模态应变能诊伤指标:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>Z</mi><mi>j</mi><mi>t</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>&beta;</mi><mi>j</mi><mi>t</mi></msubsup><mo>-</mo><mover><msup><mi>&beta;</mi><mi>t</mi></msup><mo>&OverBar;</mo></mover></mrow><msub><mi>&sigma;</mi><msup><mi>&beta;</mi><mi>t</mi></msup></msub></mfrac></mrow></math>]]></maths>③式中结构弯曲模态应变能变化为;<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&beta;</mi><mi>j</mi><mi>t</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>NM</mi></munderover><mfrac><mrow><mrow><mo>(</mo><msup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mi>T</mi></msup><msubsup><mi>K</mi><mrow><mi>j</mi><mn>0</mn></mrow><mi>t</mi></msubsup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><msup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>t</mi></msubsup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow><msup><msub><mi>&Phi;</mi><mi>i</mi></msub><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>t</mi></msubsup><msub><mi>&Phi;</mi><mi>i</mi></msub></mrow><mrow><mrow><mo>(</mo><msup><msub><mi>&Phi;</mi><mi>i</mi></msub><mi>T</mi></msup><msubsup><mi>K</mi><mrow><mi>j</mi><mn>0</mn></mrow><mi>t</mi></msubsup><msub><mi>&Phi;</mi><mi>i</mi></msub><mo>+</mo><msup><msub><mi>&Phi;</mi><mi>i</mi></msub><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>t</mi></msubsup><msub><mi>&Phi;</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup><mi>T</mi></msup><msubsup><mi>K</mi><mn>0</mn><mi>t</mi></msubsup><msubsup><mi>&Phi;</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac></mrow></math>]]></maths>④式中Φ<sub>i</sub>,Φ<sub>i</sub><sup>*</sup>分别指结构损伤前和损伤后的模态振型,其中上标″T″指转置,K<sub>j0</sub>=K<sub>j</sub>/E<sub>j</sub>,其中K<sub>j</sub>是j单元未损伤系统刚度矩阵,K<sub>j0</sub><sup>t</sup>是j单元的弯曲刚度矩阵,K<sub>0</sub><sup>t</sup>是通过单元矩阵K<sub>j0</sub><sup>t</sup>组装起来的弯曲系统矩阵;β<sub>j</sub><sup>t</sup>还表示<img file="C2005100445760003C3.GIF" wi="175" he="77" />之比,E<sub>j</sub><sup>t</sup>,<img file="C2005100445760003C4.GIF" wi="56" he="71" />各表示结构j单元未损伤和损伤后的结构构件弹性模量,来估计上述结构物的构件损伤位置;β<sup>a</sup>、β<sup>t</sup>是对应指标β<sup>a</sup>和β<sup>t</sup>的均值,<img file="C2005100445760003C5.GIF" wi="86" he="56" /><img file="C2005100445760003C6.GIF" wi="50" he="47" />是指标准方差;所述的(3)损伤程度评估步骤中分析评估结构损伤程度是,利用上述结构物的拉压模态应变能变化指标即<maths num="0005"><![CDATA[<math><mrow><msub><mi>SD</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>&alpha;</mi><mi>j</mi><mi>a</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>E</mi><mi>j</mi><msup><mi>a</mi><mo>*</mo></msup></msubsup><mo>-</mo><msubsup><mi>E</mi><mi>j</mi><mi>a</mi></msubsup></mrow><msubsup><mi>E</mi><mi>j</mi><mi>a</mi></msubsup></mfrac><mo>=</mo><mfrac><mn>1</mn><msubsup><mi>&beta;</mi><mi>j</mi><mi>a</mi></msubsup></mfrac><mo>-</mo><mn>1</mn></mrow></math>]]></maths>⑤式中SD<sub>j</sub>≥-1,来估计上述结构物的构件损伤程度;和利用上述结构物的弯曲模态应变能变化指标即<maths num="0006"><![CDATA[<math><mrow><msub><mi>SD</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>&alpha;</mi><mi>j</mi><mi>t</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>E</mi><mi>j</mi><msup><mi>t</mi><mo>*</mo></msup></msubsup><mo>-</mo><msubsup><mi>E</mi><mi>j</mi><mi>t</mi></msubsup></mrow><msubsup><mi>E</mi><mi>j</mi><mi>t</mi></msubsup></mfrac><mo>=</mo><mfrac><mn>1</mn><msubsup><mi>&beta;</mi><mi>j</mi><mi>t</mi></msubsup></mfrac><mo>-</mo><mn>1</mn></mrow></math>]]></maths>⑥式中SD<sub>j</sub>≥-1,来估计上述结构物的构件损伤程度;确定上述结构物的构件损伤程度判断指标:如果SD<sub>j</sub>=0,则表示上述结构物的j单元结构构件没有损伤;如果SD<sub>j</sub>=-1,则表示上述结构物的j单元结构构件的刚度完全丧失;结构构件的损伤程度评估指标数值在:-1≤SD<sub>j</sub>≤0范围内。
地址 266003山东省青岛市市南区鱼山路5号