发明名称 |
一种基于气浮扭摆台的转动惯量的频域检测方法 |
摘要 |
一种基于气浮扭摆台的转动惯量的频域检测方法,属于转动惯量检测技术领域。本发明解决了现有的扭摆法转动惯量的时域测量方法受到周期测量精度的影响,最终影响转动惯量的测量精度的问题。本发明的技术要点为:阻尼的线性模型条件下,根据转动定律可写出扭摆运动微分方程,解该方程可得扭摆运动转角关于时间的函数,对表达式做近似处理,对其乘积项做幂级数展开忽略其高阶项,然后对表达式进行连续傅里叶变换可得解析式,最后根据峰值频率和阻尼比可计算出被测产品无阻尼自振频率的大小,进而得到产品转动惯量的大小。本发明主要应用于导弹、拦截器、无人机等飞行器的转动惯量测量中,也可应用到其他大型装备的转动惯量测量中。 |
申请公布号 |
CN104155054B |
申请公布日期 |
2016.09.21 |
申请号 |
CN201410412481.7 |
申请日期 |
2014.08.20 |
申请人 |
哈尔滨工业大学 |
发明人 |
张晓琳;张烈山;唐文彦;李加福 |
分类号 |
G01M1/10(2006.01)I |
主分类号 |
G01M1/10(2006.01)I |
代理机构 |
哈尔滨市松花江专利商标事务所 23109 |
代理人 |
杨立超 |
主权项 |
一种基于气浮扭摆台的转动惯量的频域检测方法,其特征在于所述方法包括以下步骤:步骤一、根据虎克定律,扭杆受扭转而产生的恢复力矩M<sub>n</sub>与所转过的角度θ成正比,即M<sub>n</sub>=‑Kθ (1)式中,K为扭杆的刚度系数,考虑空气阻力、轴承机械摩擦产生的阻尼力矩时,气浮扭摆台的合力矩表示为扭杆恢复力矩和阻尼力矩之和,即M=M<sub>n</sub>+M<sub>z</sub> (2)式中,M<sub>n</sub>为扭杆作用下的恢复力矩,M<sub>z</sub>为阻尼作用产生的阻尼力矩,阻尼的线性模型条件下,阻尼力矩用下式表达:M<sub>Z</sub>=c*ω (3)c为阻尼系数,ω为摆动的瞬时角速度;步骤二、根据转动定律写出扭摆运动微分方程,得到扭摆运动转角θ关于时间t的函数表达式:<maths num="0001"><math><![CDATA[<mrow><mi>θ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>θ</mi><mn>0</mn></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>ζω</mi><mi>n</mi></msub><mi>t</mi><mo>)</mo></mrow><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>ω</mi><mi>n</mi></msub><mi>t</mi><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001017397060000011.GIF" wi="1358" he="143" /></maths>式中,ζ为阻尼比,θ<sub>0</sub>为初始摆角,<img file="FDA0001017397060000012.GIF" wi="212" he="127" />ω<sub>n</sub>为无阻尼自振频率,K为扭杆的刚度系数,I为转动惯量;步骤三、对式(4)做近似处理,对其乘积项exp(‑ζω<sub>n</sub>t)做幂级数展开忽略其高阶项,然后对θ(t)进行连续傅里叶变换得解析式:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>F</mi><mrow><mo>(</mo><mi>ω</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>θ</mi><mn>0</mn></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt></mfrac><mo>{</mo><mi>π</mi><mo>[</mo><mi>δ</mi><mrow><mo>(</mo><msub><mi>ω</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mi>ω</mi><mo>)</mo></mrow><mo>+</mo><mi>δ</mi><mrow><mo>(</mo><mo>-</mo><mi>ω</mi><mo>-</mo><msub><mi>ω</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>)</mo></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mfrac><mrow><msup><mi>πζ</mi><mn>2</mn></msup><msup><msub><mi>ω</mi><mi>n</mi></msub><mn>2</mn></msup></mrow><mn>2</mn></mfrac><mo>[</mo><msup><mi>δ</mi><mrow><mo>′</mo><mo>′</mo></mrow></msup><mrow><mo>(</mo><msub><mi>ω</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mi>ω</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>δ</mi><mrow><mo>′</mo><mo>′</mo></mrow></msup><mrow><mo>(</mo><mo>-</mo><mi>ω</mi><mo>-</mo><msub><mi>ω</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>)</mo></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><msub><mi>πζω</mi><mi>n</mi></msub><mo>[</mo><msup><mi>δ</mi><mrow><mo>′</mo><mo>′</mo></mrow></msup><mrow><mo>(</mo><msub><mi>ω</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mi>ω</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>δ</mi><mo>′</mo></msup><mrow><mo>(</mo><mo>-</mo><mi>ω</mi><mo>-</mo><msub><mi>ω</mi><mi>n</mi></msub><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mi>ζ</mi><mn>2</mn></msup></mrow></msqrt><mo>)</mo></mrow><mo>]</mo><mi>i</mi><mo>}</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001017397060000013.GIF" wi="1558" he="381" /></maths>由式(5)可知,在频域内表达扭摆运动时,<img file="FDA0001017397060000014.GIF" wi="333" he="103" />处将出现频谱峰值,根据峰值频率f和阻尼比ζ计算出被测对象无阻尼自振频率的大小,进而得到被测对象的转动惯量;δ为单位冲激函数,δ′为单位冲激函数的一阶导数,δ″为单位冲激函数的二阶导数。 |
地址 |
150001 黑龙江省哈尔滨市南岗区西大直街92号 |