发明名称 一种机械装备结构运行状态下的频响函数识别方法
摘要 本发明公开了一种机械装备结构运行状态下的频响函数辨识方法,其将工作模态分析方法辨识所得结构的固有频率与阻尼比和实验模态分析方法辨识所得结构的模态振型相结合,合成后得到频响函数,其包括:(1)对装备进行实验模态分析实验,得到结构的模态振型;(2)对装备进行工作模态分析实验,得到结构固有频率和阻尼比;(3)对(1)中得到的模态振型和(2)中得到的固有频率和阻尼比按一定方式进行合成,得到结构工作状态下频响函数。本发明的方法弥补了现有工作模态分析方法不能得到结构频响函数的不足,进一步提高了OMA方法的实用性。
申请公布号 CN106502199A 申请公布日期 2017.03.15
申请号 CN201610914557.5 申请日期 2016.10.20
申请人 华中科技大学 发明人 毛新勇;彭伊丽;李斌;秦潮;和焕斌;荣闹;徐宛丛;冯勤龙;万文杰
分类号 G05B19/18(2006.01)I;G01M7/02(2006.01)I 主分类号 G05B19/18(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 李佑宏
主权项 一种机械装备结构运行状态下的频响函数辨识方法,其特征在于,该方法包括如下步骤:(1)对机械装备结构进行试验模态分析(EMA)实验,相应获得包含所设定的n阶模态参数的频响函数[H(ω)]<sub>0</sub>,其中模态参数包括固有频率ω<sub>0n,r</sub>、阻尼比ξ<sub>0r</sub>和模态阵型{φ}<sub>0r</sub>:<maths num="0001"><math><![CDATA[<mrow><msub><mrow><mo>&lsqb;</mo><mi>H</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>0</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><mrow><msub><mrow><mo>{</mo><mi>&phi;</mi><mo>}</mo></mrow><mrow><mn>0</mn><mi>r</mi></mrow></msub><msubsup><mrow><mo>{</mo><mi>&phi;</mi><mo>}</mo></mrow><mrow><mn>0</mn><mi>r</mi></mrow><mi>T</mi></msubsup></mrow><mrow><mo>-</mo><msup><mi>&omega;</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msub><mi>j&xi;</mi><mrow><mn>0</mn><mi>r</mi></mrow></msub><msub><mi>&omega;</mi><mrow><mn>0</mn><mi>n</mi><mo>,</mo><mi>r</mi></mrow></msub><mi>&omega;</mi><mo>+</mo><msubsup><mi>&omega;</mi><mrow><mn>0</mn><mi>n</mi><mo>,</mo><mi>r</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001134559890000011.GIF" wi="1309" he="175" /></maths>其中:j为虚数单位,ω为自变量,单位为rad/s,{φ}<sub>0r</sub>为敲击实验得到的第r阶模态的模态振型,ω<sub>0n,r</sub>为敲击实验得到的第r阶模态的固有频率,ξ<sub>0r</sub>为敲击实验得到的第r阶模态的阻尼比;(2)对机械装备结构进行工作模态分析实验,相应获得反映各阶模态的振动响应功率谱,该功率谱中具有多个极点,并且各个极点分别对应于一组用于表征机械装备结构动力学特征的模态参数,即第r阶模态的固有频率ω<sub>1n,r</sub>和阻尼比ξ<sub>1r</sub>;(3)将步骤(2)所得固有频率ω<sub>1n,r</sub>、阻尼比ξ<sub>1r</sub>代入步骤(1)的频响函数[H(ω)]<sub>0</sub>中,分别代替ω<sub>0n,r</sub>、ξ<sub>0r</sub>,即可获得合成的频响函数[H(ω)]:<maths num="0002"><math><![CDATA[<mrow><mo>&lsqb;</mo><mi>H</mi><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><mrow><msub><mrow><mo>{</mo><mi>&phi;</mi><mo>}</mo></mrow><mrow><mn>0</mn><mi>r</mi></mrow></msub><msubsup><mrow><mo>{</mo><mi>&phi;</mi><mo>}</mo></mrow><mrow><mn>0</mn><mi>r</mi></mrow><mi>T</mi></msubsup></mrow><mrow><mo>-</mo><msup><mi>&omega;</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msub><mi>j&xi;</mi><mrow><mn>1</mn><mi>r</mi></mrow></msub><msub><mi>&omega;</mi><mrow><mn>1</mn><mi>n</mi><mo>,</mo><mi>r</mi></mrow></msub><mi>&omega;</mi><mo>+</mo><msubsup><mi>&omega;</mi><mrow><mn>1</mn><mi>n</mi><mo>,</mo><mi>r</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001134559890000012.GIF" wi="1342" he="183" /></maths>该频响函数[H(ω)]即为机械装备结构在运行状态下的频响函数。
地址 430074 湖北省武汉市洪山区珞喻路1037号