航天器的转动惯量等物理量的变化对于其在空间轨道的运行状态、运动姿态有直接影响,转动惯量的准确测量对于航天器等设计和控制具有重要的意义。单轴气浮平台是一种测量大型不规则物体转动惯量以及角速度等参数的实验装置,用于模拟航天器在空间微重力环境下的单轴运动。文中将在有限元仿真软件和现场实验的基础上,系统地分析气浮平台调平精度、弹簧刚度、系统阻尼大小等因素对转动惯量测量的影响。
1 单轴气浮平台机械结构单轴气浮平台测量角动量的方法为扭摆法,文献中较常见的基于扭摆法的转台系统,其基本结构通常由扭杆、角接触球轴承和工作台组成[1]。扭杆是转台系统的重要部件,利用它的储存能量的性能,实现工作台和待测物周期性摆动。文中所述单轴气浮平台,采用气浮轴承代替角接触球轴承,利用压缩空气将气浮平台悬浮,增加了平台的承载能力,轴承和轴承座之间形成气膜,又可极大程度地减少工作台运动中的摩擦阻力,实现较高精度的角动量测量,同时使得转台系统不再局限于扭摆运动,还可实现小摩擦的连续圆周转动,用于完成角速度的测量。文中的单轴气浮平台采用相对于旋转中心对称的一对弹簧代替扭杆,气浮平台结构如图 1所示。工作台被气浮轴承浮起之后,通过弹簧与平台底座相连,平台扭摆时通过交替地压缩和拉伸弹簧产生扭矩。
基于扭摆法的转动惯量测量方法[2-3],是利用给出气浮台初始角度后,使其在弹簧恢复力作用下自由摆动,记录角度数据、后期处理进而完成转动惯量测量的一种方法。单轴气浮平台在做小角度自由扭摆运动时,弹簧产生的恢复力与其形变成正比,即F=-kx,转换为相对于回转中心的扭矩M1=-Kθ,其中k为弹簧的刚度,K为弹簧等效扭转刚度,弹性限度内为常数,θ为平台转角;空气阻力与扭摆运动中的转动角速度成正比,即
$ J\frac{{{{\rm{d}}^2}\theta }}{{{\rm{d}}{\mathit{t}^2}}} + C\frac{{{\rm{d}}\theta }}{{{\rm{d}}\mathit{t}}} + K\theta = 0 $ |
定义
$ \frac{{{{\rm{d}}^2}\theta }}{{{\rm{d}}{\mathit{t}^2}}} + 2\xi {\omega _0}\frac{{{\rm{d}}\theta }}{{{\rm{d}}\mathit{t}}} + \omega _0^2\theta = 0 $ |
其特征方程为
$ {r^2} + 2\xi {\omega _0}r + \omega _0^2= 0 $ |
特征根为
$ {r_{1, 2}} = - \xi {\omega _0} \pm {\omega _0}\sqrt {{\xi ^2} - 1} $ |
假如ξ<1,即弱阻尼情况时,
$ \theta \left( t \right) = \frac{{{\theta _0}}}{{\sqrt {1 - {\xi ^2}} }}{{\rm{e}}^{\left( { - \xi {\omega _0}t} \right)}}{\rm{cos}}\left( {{\omega _0}t\sqrt {1 - {\xi ^2}} } \right) $ | (1) |
定义有阻尼振动频率
$ J = \frac{K}{{\omega _0^2 }} = \frac{K}{{\left( {2\pi } \right){^2}f_0^2}} = \frac{K}{{\left( {2\pi } \right){^2}f_d^2}}\sqrt {1 - {\xi ^2}} $ | (2) |
由于弹簧等效扭转刚度K不易准确测量,为避免因测量引起的误差,文中采用在气浮平台上附加已知质量为m0、半径为r的圆柱体砝码,砝码质量均匀,其质心距离气浮平台回转中心的距离为L,根据转动惯量计算中的平行轴定理,可知气浮台增加的转动惯量为
$ \Delta J = 2{m_0}{L^2} + {m_0}{r^2} $ |
此时根据式(2)气浮平台的转动惯量变为
$ J + \Delta J = \frac{K}{{\left( {2\pi } \right){^2}f_{d1}^2}}\sqrt {1 - {\xi ^2}} $ | (3) |
式中fd1为增加砝码之后气浮平台的有阻尼振动频率,将式(2)和(3)进行运算,得到
$ J = \frac{{f_{{\rm{d1}}}^2}}{{f_{\rm{d}}^2 - f_{{\rm{d1}}}^2}}\Delta J $ | (4) |
由式(4)可知,气浮平台的转动惯量测量与弹簧的等效扭转刚度无关。实验中采集平台转动角度随时间的变化曲线如图 2所示,利用MATLAB对式(1)进行数据拟合[4-7],进而得到增加砝码前后精确的频率值,以及装载待测物体之后的频率值,即可根据式(4)分别计算出单轴气浮平台的转动惯量以及待测物体的转动惯量。
气浮平台采用4个地脚螺栓与大地相连,通常情况下,难以实现完全水平,而且由于地壳运动,长时间未使用之后气浮平台水平角度也会发生变化。即使在假设气浮轴承供气均匀、台体加工精度极高的情况下,平台的倾斜依然会导致转台的质心偏离转台的旋转中心,其相当于在平台整周旋转中加入了重力引起的不平衡力矩。
3.1 力学分析气浮轴承及转台的总重为G,假设其重心为C,与转台旋转中心的距离为s,在转台平面上建立以旋转轴为z轴、扭摆运动平衡位置为x轴的坐标系。在转台倾斜的情况下,重心C随转台的摆动而摆动。当重心C在最低位置时,转台达到重力矩的平衡状态。
首先分析扭摆运动的平衡点和重力矩平衡位置相同的情况。假设转台倾斜角度如图 3,建立矢量P,方向垂直于旋转轴z,且通过重心C。
当转台旋转一定角度时,重心C沿矢量P方向的分量Gp为
${G_p} = G{\rm{sin}}\;\alpha $ |
式中α为转台回转轴与重力方向的夹角,则由此产生的扭矩为
$ {M_G} = {G_s}{\rm{sin}}\;\alpha {\rm{sin}}\;\theta $ |
则扭摆系统的运动方程为
$J\frac{{{{\rm{d}}^2}\theta }}{{{\rm{d}}{\mathit{t}^2}}} + C\frac{{{\rm{d}}\theta }}{{{\rm{d}}\mathit{t}}} + K\theta = {G_s}{\rm{sin}}\;\alpha {\rm{sin}}\;\theta $ |
扭摆运动的转动角度较小,因此θ≈sin θ,则方程可写为
$ J\frac{{{{\rm{d}}^2}\theta }}{{{\rm{d}}{\mathit{t}^2}}} + C\frac{{{\rm{d}}\theta }}{{{\rm{d}}\mathit{t}}} + N\theta = 0 $ |
式中N=(K-Gssin α)为定值,即相当于增大了弹簧的刚度,则转动惯量的测量仍可根据式(4)完成。
当运动的平衡位置和矢量P不重合时,如图 4所示。
此时扭摆运动的方程变为
$ J\frac{{{{\rm{d}}^2}\theta }}{{{\rm{d}}{\mathit{t}^2}}} + C\frac{{{\rm{d}}\theta }}{{{\rm{d}}\mathit{t}}} + K\theta = \pm Gs{\rm{sin}}\;\alpha {\rm{sin}}\;\left( {\theta - \varphi } \right) $ |
式中φ为扭摆运动平衡位置和重力矩平衡位置向量P间的夹角。此方程无法求出通常意义下的解析解,因此利用MATLAB中的ode45()函数求取方程的数值解并绘图[8-10],只要选取的时间步长足够小,文中选取步长为10 ms,就可以保证数据的精度。通过对系统转动惯量J、阻尼系数ξ、弹簧等效扭转刚度K、转台重量m以及倾斜角度α、平衡位置夹角φ等的赋值,对比观察倾斜角度α以及平衡位置夹角φ变化所引起的数值解图形的变化。MATLAB计算出的数值解曲线如图 5、6所示。
为了便于观察,图 5中仅展示了数值解的第30~50 s,图 5中虚线代表倾斜角度为1°时的数据点,实线代表倾斜角度为10°时的数据点,观察曲线可知,倾斜角度会影响转台摆动周期,随着角度增大,摆动周期随之减小,进而影响转动惯量的测量精度。
图 6中虚线和实线分别代表扭摆运动平衡位置和重力矩平衡位置向量P间的夹角为0.8π和π/3。通过对130~150 s时间段的观察可以看出:随着平衡位置夹角的改变,转台摆动周期随之改变。
在实际转台实验中,转台的调平精度会控制在1°之内,文中按照1°最大偏角极限情况进行分析。根据数值解仿真可知,最大偏角引起的周期改变小于0.001 s,且当弹簧刚度较小、转台及待测物重量较大的情况下,其相对于扭摆周期的比例极小,因此调平问题对于转动惯量测量的影响极小。
3.2 仿真分析由于实验中难以调整气浮平台的倾斜角度,进行不同倾角下的数据采集分析,本节利用有限元仿真,定量地分析了气浮平台的调平情况对于转动惯量测量的影响。
仿真中采用1 N/mm的弹簧组,阻尼设置为固定阻尼,分别设置气浮平台与水平面间的夹角为0°、0.3°、0.6°、1°、2°、3°,待测物为0.767 kg的砝码,其转动惯量计算值为0.173 331 04 kg·m2,进行仿真,数据见表 1。
根据表 1可知,随着气浮平台与地面间夹角的增加,平台的转动惯量会逐渐增大。在气浮平台的实际调平中,平台的倾斜角度会被控制在1°之内,且平台自身转动惯量较大,因此平台的调平对于平台转动惯量的影响很微小,0°~3°平台转动惯量变化仅为0.005 kg·m2。平台倾斜后待测物体的转动惯量也逐渐增大,取0°和1°的数据进行计算,因倾斜所引起的测量误差变化仅为10-5 kg·m2。转动惯量测量与平台倾斜角度关系如图 7所示。
在气浮平台的实际实验中,难以精确测量平台的实际倾斜角度,此处仅进行定性分析。调平前平台倾斜角度约为1.5°左右,调平后倾斜角度小于0.8°,同时由于存在制造误差、操作误差等随机误差的影响,平台和待测物体转动惯量的测量值会有一定波动。待测物体转动惯量理论计算值为3.947 92 kg·m2。如表 2所示,列出了调平前后转动惯量测量值。
调平前转动惯量平均值为3.956 82 kg·m2,调平后转动惯量平均值为3.954 17 kg·m2。根据实际测量数据可知,调平前后待测物体转动惯量的测量绝对误差可控制在0.01 kg·m2以内。因此对于质量较大的待测物体,该气浮平台的测量相对误差更小。
4 弹簧刚度对于转动惯量测量的影响文中所述气浮平台采用成组圆柱螺旋弹簧提供扭摆力矩,随着使用时间的增加,弹簧的刚度系数会发生变化,且由于制造等原因也会存在弹簧刚度系数不均的情况,前文利用公式推导,排除弹簧刚度系数对于转动惯量测量的影响,本节将利用有限元仿真配合实际实验对其进行验证。表 3中列出了不同弹簧刚度下转动惯量的仿真测量值。图 8为不同弹簧刚度系数下的扭摆实验数据。
根据表 3和图 8的数据可知,弹簧刚度的大小会直接影响无阻尼振动频率ω0的大小,但其对于平台转动惯量以及待测物体转动惯量的测量的影响较小。在实际实验中,弹簧刚度难以直接测量,仅更换不同弹簧进行实验。实验数据图像疏密程度类似于图 8(a)、(b),两次更换弹簧后,分别进行了4次转动惯量测量实验,实验数据的平均值分别为3.939 45和3.949 01 kg·m2,数据绝对误差小于0.01 kg·m2。实验结果与仿真分析结果基本一致。因此只要选用在刚度系数k在小变形范围内近似为定值的弹簧,且保证每次实验的初始摆角尽量相同,即可实现转动惯量的较高精度测量。
表 4数据中仅列出几组有代表性的仿真数据,综合各种弹簧刚度不同的情况,砝码转动惯量测量的误差均小于0.001 kg·m2。因此弹簧组的刚度不同所引起的测量误差较小,可不必进行考虑。本节验证了式(4)结论:转动惯量测量与弹簧的等效扭转刚度无关。
气浮平台在扭摆过程中会受到空气阻力以及气浮轴承排出气体的影响,导致平台扭摆的幅度逐渐减小,扭摆速度和角度减小。文中假设系统阻尼大小与扭摆角速度成正比,如图 9为阻尼逐渐增大情况下的气浮平台扭摆角度曲线。从图可看出随着线性阻尼的增大,振幅的衰减速度越来越快,但扭摆周期不变。文中所述的气浮平台阻尼较小,衰减情况和图 9(a)的情况接近,实际实验中为保证转动惯量的测量精度,使用有机玻璃罩住待测物体进行试验,可保证系统的阻尼情况不因待测物体的加载而改变。综上,系统阻尼大小对于转动惯量几乎无影响。
文中介绍了基于单轴气浮平台的转动惯量测量方法,并利用ADAMS有限元分析软件完成仿真,结合具体现场实验得出以下结论:
1) 单轴气浮平台的倾斜角度以及弹簧组刚度的不均对于转动惯量的测量影响很小,测量绝对误差均可控制在0.01 kg·m2以内。
2) 文中介绍的转动惯量测量方法具有很高的测量精度。
3) 文中弹簧刚度和空气阻尼均为线性,而在实际实验中,当扭摆角度过大时,弹簧的刚度系数会呈现非线性性质,文中所述转动惯量测量原理和各种测量影响因素结论,还可在加入非线性因素后进行更进一步的研究。
[1] | 秦远田, 程月华. 基于气浮台的卫星惯量测试方法[J]. 理论与方法, 2013, 32(2): 17-19, 22. (0) |
[2] | 穆继亮. 基于扭摆法的弹体转动惯量测量系统及误差分析[J]. 机械工程与自动化, 2009(1): 103-105. (0) |
[3] | 田留德, 赵建科, 薛勋, 等. 基于单轴气浮台的角动量输出测量方法[J]. 中国测试, 2011, 37(6): 41-44. (0) |
[4] | 刘寅琛, 崔剑, 王亮. 单轴气浮平台的角动量测试系统设计[J]. 单片机与嵌入式系统应用, 2014(4): 36-39. (0) |
[5] | 武颖丽, 李平舟. 扭摆法测量刚体转动惯量的误差分析[J]. 大学物理实验, 2015, 28(3): 101-103. (0) |
[6] | BERNSTEIN D S, MCCLAMROCH N H, BLOCH A. Development of air spindle and triaxial air bearing testbeds for spacecraft dynamics and control experiments[C]//Proceedings of the 2001 American Control Conference. Arlington, VA, USA, 2001:3967-3972. (0) |
[7] | 齐国清, 吕健. 正弦曲线拟合若干问题探讨[J]. 计算机工程与设计, 2008, 29(14): 3677-3680. (0) |
[8] | 唐家德. 基于MATLAB的非线性曲线拟合[J]. 计算机与现代化, 2008(6): 15-19. (0) |
[9] | 赵岩. 扭摆法转动惯量测量中的非线性问题研究[D]. 哈尔滨: 哈尔滨工业大学, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10213-1014080676.htm (0) |
[10] | 林洪文, 马强, 唐文彦, 等. 非线性阻尼扭摆振动模型及在转动惯量测量中应用[J]. 振动与冲击, 2014, 33(10): 13-16. (0) |
[11] | 李广龙, 魏政君, 上官文斌. 非线性试验数据的拟合方法[J]. 新技术新工艺, 2016(8): 18-21. (0) |
[12] | 冯晓伟, 李正生. 高速衰减正弦信号的采样及曲线拟合方法[J]. 兵工自动化, 2011, 30(10): 75-78. DOI:10.3969/j.issn.1006-1576.2011.10.023 (0) |
[13] | HANDEL P. Properties of the IEEE-STD-1057 four-parameter sine wave fit algorithm[J]. IEEE transactions on instrumentation and measurement, 2000, 49(6): 1189-1193. DOI:10.1109/19.893254 (0) |
[14] | 陈向华, 赵国忠. 非线性单摆运动的数值解[J]. 内蒙古科技大学学报, 2007, 26(1): 94-96. (0) |
[15] | 林振华, 董云峰. 基于单轴气浮台摆动特性的调节平衡方法[J]. 科技导报, 2010, 28(2): 46-49. (0) |