主权项 |
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>ω</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>ω</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>ω</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>Y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>ω</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>F</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>ω</mi><mo>)</mo></mrow><mo>+</mo><mi>m</mi><msub><mi>Y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>ω</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>ω</mi><mi>A</mi></msup><mrow><mo>(</mo><msup><mi>μ</mi><mn>0</mn></msup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>{</mo><msubsup><mi>ω</mi><mn>1</mn><mi>A</mi></msubsup><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mo>,</mo><msubsup><mi>ω</mi><mn>2</mn><mi>A</mi></msubsup><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><msubsup><mi>ω</mi><mi>n</mi><mi>A</mi></msubsup><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow><mo>}</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>ω</mi><mi>i</mi><mi>A</mi></msubsup><mrow><mo>(</mo><mi>μ</mi><mo>)</mo></mrow></mrow></math>]]></maths>表示结构的损伤因子为μ时在第i个构件上附加最佳质量后所挑选出来的灵敏度最高的频率,Δω为实测构造频率与理论有限元频率的差值。 |