发明名称 基于附加虚拟质量的框架结构损伤方法
摘要 本发明公开了一种基于附加虚拟质量的框架结构损伤方法,属于质量结构监测技术领域。其特征是,利用一个加速度传感器和模态力锤对框架中各构件分别测量来识别框架结构的损伤。对框架中的构件逐个进行局部动力测试,利用虚拟变形法由测量的加速度响应和激励时程构造加虚拟质量后的频率响应,从而实现在各个构件上附加虚拟质量,提高各构件对损伤的灵敏度。本发明克服了框架结构对局部损伤不敏感和损伤定位的困难,实际操作简单、成本低廉和识别准确等特点,具有很好的应用前景。
申请公布号 CN102539537B 申请公布日期 2013.10.16
申请号 CN201110454236.9 申请日期 2011.12.30
申请人 大连理工大学 发明人 侯吉林;李冬生;喻言;欧进萍
分类号 G01N29/12(2006.01)I;G01N29/44(2006.01)I 主分类号 G01N29/12(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉;李宝元
主权项 1.一种基于附加虚拟质量的框架结构损伤方法,其特征在于,利用一个加速度传感器和一个模态力锤对框架结构中各个构件进行局部动力测试,如果框架由n个构件组成,则需进行n次局部动力测试;在第i次局部动力测试时,将传感器垂直于第i个构件布置在其中间,用模态力锤沿传感器测量方向和位置激励,同时测量加速度响应y<sub>i</sub>(t)和激励f<sub>i</sub>(t),对测量数据进行傅立叶变换分别获得频域的Y<sub>i</sub>(ω)和F<sub>i</sub>(ω);利用理论有限元模型计算附加质量m后各阶频率关于第i个构件损伤的灵敏度,从中估计和挑选出最佳质量m<sub>i</sub>和模态阶数k<sub>i</sub>,使其相对灵敏度达到最大,即m<sub>i</sub>和k<sub>i</sub>满足公式(4);<img file="FDA00003349754700011.GIF" wi="1532" he="118" />式(4)中<img file="FDA00003349754700012.GIF" wi="398" he="188" />为在第i个构件附加质量m后第j阶频率关于第i个构件的相对灵敏度,<img file="FDA00003349754700015.GIF" wi="832" he="189" />为对应的第j阶频率关于第i个构件的灵敏度计算公式,ω<sub>ji</sub>(m)为在第i个构件中间附加质量m所对应理论模型的第j阶频率;将m<sub>i</sub>、Y<sub>i</sub>(ω)和F<sub>i</sub>(ω)代入公式(3)构造出附加质量后的频率响应<img file="FDA00003349754700016.GIF" wi="241" he="70" />由<img file="FDA00003349754700017.GIF" wi="214" he="72" />识别出频率,并将第k<sub>i</sub>阶频率挑选出来,记为<img file="FDA00003349754700018.GIF" wi="91" he="62" /><maths num="0001"><![CDATA[<math><mrow><msubsup><mi>H</mi><mi>i</mi><mi>v</mi></msubsup><mrow><mo>(</mo><mi>&omega;</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>H</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mi>m</mi><msub><mi>H</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>Y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>F</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>+</mo><mi>m</mi><msub><mi>Y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>框架含有n个构件,针对每个构件都在附加一定虚拟质量后构造出一个对其灵敏度最高的频率<img file="FDA00003349754700019.GIF" wi="457" he="76" />那么构造n个频率,每个频率都会对对应的构件具有最高的灵敏度,将这n个构造的频率<img file="FDA000033497547000110.GIF" wi="63" he="63" />作为优化的目标,结合理论有限元模型,由公式(8)和(9)利用基于灵敏度的方法迭代优化各个构件的损伤因子;<img file="FDA00003349754700021.GIF" wi="1643" he="654" />Δω=β<sup>A</sup>Δμ               (9);上式(8)和(9)中参数含义如下:μ={μ<sub>1</sub>,μ<sub>2</sub>,…,μ<sub>n</sub>}<sup>T</sup>为构件损伤因子排列而成的列向量,为待识别的值,<img file="FDA00003349754700023.GIF" wi="456" he="87" />为有限元损伤因子的初始值,Δμ为损伤因子的实际值与初始值的差值,<img file="FDA00003349754700026.GIF" wi="606" he="92" />为利用实测数据构造出的附加虚拟质量m<sub>i</sub>后,关于第i个构件灵敏度最高的频率,<maths num="0002"><![CDATA[<math><mrow><msup><mi>&omega;</mi><mi>A</mi></msup><mrow><mo>(</mo><msup><mi>&mu;</mi><mn>0</mn></msup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>{</mo><msubsup><mi>&omega;</mi><mn>1</mn><mi>A</mi></msubsup><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>,</mo><msubsup><mi>&omega;</mi><mn>2</mn><mi>A</mi></msubsup><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msubsup><mi>&omega;</mi><mi>n</mi><mi>A</mi></msubsup><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>}</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&omega;</mi><mi>i</mi><mi>A</mi></msubsup><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow></mrow></math>]]></maths>表示结构的损伤因子为μ时在第i个构件上附加最佳质量后所挑选出来的灵敏度最高的频率,Δω为实测构造频率与理论有限元频率的差值。
地址 116024 辽宁省大连市凌工路2号