发明名称 一种回转体立式极转动惯量测量方法
摘要 本发明涉及一种回转体立式极转动惯量测量方法,具体包括以下步骤:(1)采用转动惯量测量台来测量,在该测量台上能够测得回转体的极转动惯量和赤道转动惯量。(2)摆动周期测量:周期通过采样时间序列进行测量,采用光电计时系统;(3)回转体转动惯量的测量:用扭摆法测量转动惯量时将回转体绕指定轴自由摆动,测出摆动周期,然后由摆动周期计算出回转体转动惯量。(4)转动惯量测量误差分析。(5)阻尼对赤道转动惯量测试结果的影响。该发明能有效地针对回转体立式极转动惯量进行测量,使用方便,改善了测试效果,方便根据需要使用。
申请公布号 CN105157918A 申请公布日期 2015.12.16
申请号 CN201510250711.9 申请日期 2015.05.18
申请人 长春理工大学 发明人 李俊烨;赵友;戴正国;王德民;张心明;宋斌;沈军;刘建河
分类号 G01M1/10(2006.01)I 主分类号 G01M1/10(2006.01)I
代理机构 代理人
主权项 一种回转体立式极转动惯量测量方法,其特征在于:具体包括以下步骤:(1)采用转动惯量测量台来测量,在该测量台上能够测得回转体的极转动惯量和赤道转动惯量;该惯量测量台在立式装夹回转体时可以测量回转体的极转动惯量,在横置装夹时可以测量回转体的赤道转动惯量;该惯量测量台为了适应较小回转体的惯量测量,夹具尽可能的采用了轻质材料如铝合金等;该测量台由夹具、丝杠、释放机构、导轨和基座组成;夹具的装夹部位设计成横、纵向卡槽结构,能够实现回转体的横向和纵向的定位装夹,导轨丝杠的设计能够保证回转体的夹紧并使夹具的回转轴线与扭杆的回转轴线重合;回转体赤道转动惯量测量台由基座、支架、导轨和装夹机构组成,测量时将回转体放置在支架上,一端固定在装夹机构内,施力释放后,回转体在装夹机构上的弹簧作用下沿着水平轴线方向往复转动,进而测得回转体的赤道转动惯量;回转体极转动惯量测量台由基座、扭摆机构、支撑和测量平台组成;测量时,将回转体放置在测量平台的V型槽上,并使回转体的质心位置尽量与扭摆机构的回转中心线重合,对测量平台施力释放后,测量平台与回转体一起在扭摆弹簧的作用下沿着水平方向往复摆动,进而测得回转体的极转动惯量;(2)摆动周期测量:周期通过采样时间序列进行测量,采用光电计时系统,用高速单周期c8051F单片机定时器组成计时电路,分辨率为10ns,满足转动惯量测量相对误差小于±0.5%的要求;为了避免环境光干扰,采用交流调制光完成光电采集,在信号处理中,需要对定时器开始定时的时刻进行准确的测量,电路中为了抗干扰,采用了具有一定门限的电平参考和一定延时的展宽整形电路,但此时定时开始时刻必定有一个随机误差,同理在下一个计数脉冲到达时也必然在准确时刻上存在一个随机误差,作为一个完整周期的测量用两个脉冲之间的时间来表征周期量显然不可能达到要求,为了最大限度地降低该误差,根据振动周期基本不变的条件,测量多个周期并取平均数作为周期的测得值;在此时,若振动连续,则每个周期摆端将两次经过反射点,在一个周期内形成两个计数脉冲,可以知道每相邻三个脉冲之间的时间应为一个周期,即相邻两个时间段之和为一个周期;最终的处理方法是:去掉最初始有外力扰动的几个周期,采用每相邻的两个时间段为一个周期数,从有效计数脉冲开始计时,每两个计数脉冲为一次周期数据,并用总的计时时间去除得到的完整周期数进而得到该时刻的周期值;(3)回转体转动惯量的测量:用扭摆法测量转动惯量时将回转体绕指定轴自由摆动,测出摆动周期,然后由摆动周期计算出回转体转动惯量;由力学分析可得摆动物体的摆动方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><mi>&theta;</mi></mrow><msup><mi>dt</mi><mn>2</mn></msup></mfrac><mo>+</mo><mn>2</mn><mi>n</mi><mfrac><mi>d&theta;</mi><mi>dt</mi></mfrac><mo>+</mo><mfrac><mi>k&theta;</mi><mi>J</mi></mfrac><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000021.GIF" wi="528" he="150" /></maths>式中,J为物体的转动惯量;n=C/2J,n为空气阻尼系数;θ为物体的转角(即摆角);k为扭杆的刚度系数;一般情况下,空气阻尼很小,即n<<1,则<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><mfrac><mi>k</mi><mrow><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup></mrow></mfrac><msup><mi>T</mi><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA0000718662830000022.GIF" wi="276" he="148" /></maths>令<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfrac><mi>k</mi><mrow><mn>4</mn><msup><mi>&pi;</mi><mn>2</mn></msup></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000718662830000023.GIF" wi="218" he="142" /></maths>则J=AT<sup>2</sup>   (2)式(2)是一种理想情况,因为回转体必须装在摆动体上进行摆动,因此有J<sub>d</sub>+J<sub>o</sub>=AT<sub>d</sub><sup>2</sup>   (3)式中,J<sub>d</sub>是回转体的转动惯量;J<sub>o</sub>是摆动体本身的转动惯量;A是与摆动体结构有关的常数;T<sub>d</sub>为回转体和摆动体一起摆动时的摆动周期;为求出J<sub>d</sub>必须首先测出摆动体空摆的摆动周期T<sub>o</sub>:J<sub>o</sub>=AT<sub>o</sub><sup>2</sup>   (4)然后测出标准件与摆动物体一起摆动时的摆动周期Ts:J<sub>s</sub>+J<sub>o</sub>=AT<sub>s</sub><sup>2</sup>   (5)由式(3)、式(4)和式(5)式可得:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>J</mi><mi>d</mi></msub><mo>=</mo><msub><mi>J</mi><mi>s</mi></msub><mfrac><mrow><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup><mo>-</mo><msubsup><mi>T</mi><mi>o</mi><mn>2</mn></msubsup></mrow><mrow><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>-</mo><msubsup><mi>T</mi><mi>o</mi><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000024.GIF" wi="418" he="168" /></maths>式中,J<sub>s</sub>是标准样件的转动惯量,预先给出出,因此,可以根据测出的T<sub>o</sub>、T<sub>d</sub>、T<sub>s</sub>以及已知的J<sub>s</sub>计算出J<sub>d</sub>;(4)转动惯量测量误差分析:测量误差由式(3)可得:ΔJ<sub>d</sub>=2AT<sub>d</sub>·ΔT<sub>d</sub>略去J<sub>0</sub>后可得:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&mu;</mi><msub><mi>J</mi><mi>d</mi></msub></msub><mo>=</mo><mfrac><msub><mi>&Delta;J</mi><mi>d</mi></msub><msub><mi>J</mi><mi>d</mi></msub></mfrac><mo>&ap;</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>AT</mi></mrow><mi>d</mi></msub><mo>&CenterDot;</mo><msub><mi>&Delta;T</mi><mi>d</mi></msub></mrow><msup><msub><mi>AT</mi><mi>d</mi></msub><mn>2</mn></msup></mfrac><mo>&ap;</mo><mn>2</mn><mfrac><msub><mi>&Delta;T</mi><mi>d</mi></msub><msub><mi>T</mi><mi>d</mi></msub></mfrac><mo>=</mo><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><msub><mi>T</mi><mi>d</mi></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000025.GIF" wi="910" he="158" /></maths>由(7)可见,转动惯量测量的相对误差,仅为时间测试相对误差的两倍,是很小的;经过推导对以上公式进行了进一步的整理,有:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&mu;</mi><msub><mi>J</mi><mi>d</mi></msub></msub><mo>=</mo><msup><mrow><mo>[</mo><mfrac><msup><msub><mrow><mn>4</mn><mi>T</mi></mrow><mi>d</mi></msub><mn>2</mn></msup><msup><mrow><mo>(</mo><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><msup><msub><mrow><mn>4</mn><mi>T</mi></mrow><mi>s</mi></msub><mn>2</mn></msup><msup><mrow><mo>(</mo><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><mn>4</mn><msup><mrow><mo>(</mo><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><msup><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></msup></mrow><mrow><msup><mrow><mo>(</mo><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup><msup><mrow><mo>(</mo><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>&CenterDot;</mo><mo>|</mo><mi>&Delta;T</mi><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000031.GIF" wi="1390" he="167" /></maths>由此式还可以估算出转动惯量J<sub>0</sub>,J<sub>s</sub>和J<sub>d</sub>之间的关系;计算验证表明,误差满足技术指标1%要求;位置误差:a)回转体纵向轴线位置偏移对极转动惯量精度影响:转轴偏移对极转动惯量的影响,设偏心为e,转轴与回转体形心轴偏离距离为△回转体对形心轴极转动惯量为:J<sub>极(形心)</sub>=J<sub>极(质心)</sub>+Me<sup>2</sup>回转体对偏离的转轴的极转动惯量为:J<sub>极(偏心)</sub>=J<sub>极(质心)</sub>+(e<sup>2</sup>+Δ<sup>2</sup>‑2e·Δ·cosθ)MJ<sub>极(偏心)</sub>=J<sub>极(形心)</sub>+Δ<sup>2</sup>·M‑2e·Δ·cosθ·MΔJ<sub>极</sub><sup>max</sup>=J<sub>极(偏心)</sub><sup>max</sup>‑J<sub>极(形心)</sub>=MΔ<sup>2</sup>+2e·Δ·M极转动惯量的相对误差(最大):<img file="FDA0000718662830000032.GIF" wi="803" he="164" />相对误差大小,与质偏有关;b)转轴偏移对赤道转动惯量的影响:设y轴是通过质心轴线,那么回转体对y轴的转动惯量为(J<sub>y</sub>)由于定位误差造成的以y'轴为转轴,那么J<sub>y'</sub>与J<sub>y</sub>差值就是误差;由平行轴定理可知:J<sub>y'</sub>=J<sub>y</sub>+d<sup>2</sup>Mc)转轴倾斜对极和赤道转动惯量的影响:当倾斜α角时应用转轴公式得:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>J</mi><mi>x</mi></msub><mo>=</mo><msup><mi>cos</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>&CenterDot;</mo><msub><mi>J</mi><msub><mi>x</mi><mn>1</mn></msub></msub><mo>+</mo><msub><mi>J</mi><msub><mi>y</mi><mn>1</mn></msub></msub><msup><mi>sin</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>+</mo><msub><mrow><mn>2</mn><mi>J</mi></mrow><mrow><msub><mi>x</mi><mn>1</mn></msub><msub><mi>y</mi><mn>1</mn></msub></mrow></msub><mi>cos</mi><mi>&alpha;</mi><mo>&CenterDot;</mo><mi>sin</mi><mi>&alpha;</mi></mrow>]]></math><img file="FDA0000718662830000033.GIF" wi="950" he="80" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>J</mi><mi>y</mi></msub><mo>=</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>&CenterDot;</mo><msub><mi>J</mi><msub><mi>x</mi><mn>1</mn></msub></msub><mo>+</mo><msub><mi>J</mi><msub><mi>y</mi><mn>1</mn></msub></msub><msup><mi>cos</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>+</mo><msub><mrow><mn>2</mn><mi>J</mi></mrow><mrow><msub><mi>x</mi><mn>1</mn></msub><msub><mi>y</mi><mn>1</mn></msub></mrow></msub><mi>cos</mi><mi>&alpha;</mi><mo>&CenterDot;</mo><mi>sin</mi><mi>&alpha;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000034.GIF" wi="1078" he="98" /></maths>当α角很小,且J<sub>xy</sub>也很小时,上式可以简化为:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>J</mi><mi>x</mi></msub><mo>=</mo><msup><mi>cos</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>&CenterDot;</mo><msub><mi>J</mi><msub><mi>x</mi><mn>1</mn></msub></msub><mo>+</mo><msub><mi>J</mi><msub><mi>y</mi><mn>1</mn></msub></msub><msup><mi>sin</mi><mn>2</mn></msup><mi>&alpha;</mi></mrow>]]></math><img file="FDA0000718662830000041.GIF" wi="561" he="81" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>J</mi><mi>y</mi></msub><mo>=</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>&CenterDot;</mo><msub><mi>J</mi><msub><mi>x</mi><mn>1</mn></msub></msub><mo>+</mo><msub><mi>J</mi><msub><mi>y</mi><mn>1</mn></msub></msub><msup><mi>cos</mi><mn>2</mn></msup><mi>&alpha;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000042.GIF" wi="690" he="98" /></maths>当<img file="FDA0000718662830000043.GIF" wi="212" he="162" />α≤0.5°时,倾斜角α对极转动惯量J<sub>x</sub>相对误差的影响不超过2×10<sup>‑3</sup>,对赤道转动惯量的影响就更小了;d)转轴倾斜和偏移同时存在时对转动惯量的影响当转轴倾斜和偏移同时存在时,可由上述分析结果叠加得到;(5)阻尼对赤道转动惯量测试结果的影响:在分析阻尼扭转振动之前,先分析无阻尼摆动测转动惯量的过程;a)无阻尼测赤道转动惯量J<sub>d</sub>:测试过程如下:第一步:测出空转台的振动周期T<sub>0</sub>,此时转台上不放任何东西;第二步:测出转台和标准件一起的振动周期T<sub>s</sub>;第三步:测出转台和待测物体一起的振动周期T<sub>d</sub>;测得三个已知数T<sub>0</sub>,T<sub>s</sub>,T<sub>d</sub>,可根据振动频率与转动惯量之间的关系,列出三个方程,求解三个未知量:空转台转动惯量J<sub>o</sub>,弹簧刚度k和待测物体转动惯量J<sub>d</sub>;由以下的方程组:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msup><msub><mi>T</mi><mi>o</mi></msub><mn>2</mn></msup><mo>=</mo><msup><mrow><mn>4</mn><mi>&pi;</mi></mrow><mn>2</mn></msup><mfrac><msub><mi>J</mi><mi>o</mi></msub><mi>k</mi></mfrac></mtd></mtr><mtr><mtd><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>=</mo><msup><mrow><mn>4</mn><mi>&pi;</mi></mrow><mn>2</mn></msup><mfrac><mrow><msub><mi>J</mi><mi>s</mi></msub><mo>+</mo><msub><mi>J</mi><mn>0</mn></msub></mrow><mi>k</mi></mfrac></mtd></mtr><mtr><mtd><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup><mo>=</mo><msup><mrow><mn>4</mn><mi>&pi;</mi></mrow><mn>2</mn></msup><mfrac><mrow><msub><mi>J</mi><mi>d</mi></msub><mo>+</mo><msub><mi>J</mi><mn>0</mn></msub></mrow><mi>k</mi></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000718662830000044.GIF" wi="421" he="431" /></maths>J<sub>s</sub>是标准件转动惯量;可得:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>k</mi><mo>=</mo><mfrac><mrow><msup><mrow><mn>4</mn><mi>&pi;</mi></mrow><mn>2</mn></msup><msub><mi>J</mi><mn>0</mn></msub></mrow><msup><msub><mi>T</mi><mn>0</mn></msub><mn>2</mn></msup></mfrac><mo>,</mo><msub><mi>J</mi><mi>s</mi></msub><mo>=</mo><mfrac><msup><msub><mi>T</mi><mi>o</mi></msub><mn>2</mn></msup><mrow><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup></mrow></mfrac><msub><mi>J</mi><mn>0</mn></msub><mo>,</mo><msub><mi>J</mi><mi>d</mi></msub><mo>=</mo><mfrac><mrow><msup><msub><mi>T</mi><mi>d</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mi>o</mi></msub><mn>2</mn></msup></mrow><mrow><msup><msub><mi>T</mi><mi>s</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>T</mi><mi>o</mi></msub><mn>2</mn></msup></mrow></mfrac><msub><mi>J</mi><mi>s</mi></msub></mrow>]]></math><img file="FDA0000718662830000045.GIF" wi="1032" he="168" /></maths>b)有阻尼时测得的赤道转动惯量:假定摆动阻力矩与摆角速度一次方成正比;这时摆动方程为:<img file="FDA0000718662830000046.GIF" wi="454" he="85" />式中:k——扭簧刚度上式可化为:<img file="FDA0000718662830000047.GIF" wi="342" he="76" />式中:<img file="FDA0000718662830000051.GIF" wi="181" he="128" /><img file="FDA0000718662830000052.GIF" wi="162" he="130" />分别是阻尼系数和无阻尼频率,当n<P时,微分方程的解为:<img file="FDA0000718662830000053.GIF" wi="587" he="97" />阻尼频率:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msup><mover><mi>p</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup><mo>=</mo><msup><mi>p</mi><mn>2</mn></msup><mo>-</mo><msup><mi>n</mi><mn>2</mn></msup></mrow>]]></math><img file="FDA0000718662830000054.GIF" wi="290" he="78" /></maths>阻尼周期:<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><mover><mi>T</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mover><mi>p</mi><mo>&OverBar;</mo></mover></mfrac><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msqrt><msup><mi>p</mi><mn>2</mn></msup><mo>-</mo><msup><mi>n</mi><mn>2</mn></msup></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000718662830000055.GIF" wi="554" he="169" /></maths>再通过测量阻尼振动周期<img file="FDA0000718662830000056.GIF" wi="69" he="83" />来计算出待测物体的转动惯量J;阻尼的实用测量方法:阻尼的振动方程为:<img file="FDA0000718662830000057.GIF" wi="398" he="84" />其一般解为:<img file="FDA0000718662830000058.GIF" wi="502" he="78" />相继两个振幅之比为:<img file="FDA0000718662830000059.GIF" wi="464" he="147" /><img file="FDA00007186628300000510.GIF" wi="62" he="84" />是有阻尼情况下测得的周期,n是阻尼系数;为提高测量精度,可取相隔m个周期后的两个振幅之比;即:<img file="FDA00007186628300000511.GIF" wi="275" he="129" /><img file="FDA00007186628300000512.GIF" wi="474" he="81" /><img file="FDA00007186628300000513.GIF" wi="1291" he="146" />只要测得第i个周期的振幅和第i+m个周期的振幅及周期,便可求出阻尼系数n;例如:<img file="FDA00007186628300000514.GIF" wi="452" he="85" />T<sub>d</sub>=0.8642″则阻尼系数:<maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><mi>n</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mn>1.470</mn><mo>+</mo><mn>1.560</mn></mrow><mrow><mn>10</mn><mo>&times;</mo><mn>0.8642</mn></mrow></mfrac><mo>=</mo><mn>0.0104</mn></mrow>]]></math><img file="FDA00007186628300000515.GIF" wi="599" he="133" /></maths>经10个周期,因空气阻尼振幅减小1°,阻尼已经相当大,但计算结果表明其对周期的影响仍然可以略去,见式(13)。
地址 130022 吉林省长春市卫星路7089号长春理工大学