磁气悬吊微重力模拟系统动力学研究*

周梅 张欢 宋晓东

（北京理工大学宇航学院，北京 100081）

引言

随着太空探测和开发的不断深入，空间技术得到了急速发展，在轨操作任务的需求也在不断增加，并且任务特性呈现越发复杂繁重的特点.因此，利用机械臂代替宇航员舱外作业是现阶段和未来探索太空的必然发展趋势［1，2］.空间机械臂是典型的刚柔耦合系统，并且由于臂杆细长、结构刚度低等特点展现明显的柔性特性，在大范围刚性运动的同时伴随小幅度的柔性运动.特别是末端执行器在抓取目标时的残余柔性振动将严重影响空间机械臂的操作稳定性和末端定位精度［3-5］.因此，机械臂投入太空使用之前，需要对其进行地面微重力模拟试验指导在轨服务，复现真实的微重力环境保证在轨准确服务［6-8］.

目前，空间机械臂地面微重力模拟已经发展了吊丝配重法、静平衡法、水浮法、水浮磁悬浮混合法、气浮法等方法［9，10］.吊丝配重法是指通过滑轮组悬挂配重，并调节配重物的质量来补偿空间机械臂的重力.该方法同样可以实现空间机械臂的三维微重力模拟，但是由于吊丝和滑轮之间的摩擦以及吊丝的颤振导致系统的微重力模拟精度较差［11］.静平衡法利用钢丝、弹簧、连杆机构、滑轮等部件，遵循系统的能量守恒定律即系统的重力势能和弹性势能的总和保持不变实现微重力模拟，但是模拟精度比较差［12］.水浮法可以用于三维微重力模拟试验，但是需要对空间机械臂系统进行密封性改造，在水池中通过配重的变化使得空间机械臂在水中的浮力与重力相平衡来实现空间机械臂的微重力模拟［13］.当机械臂进行动力学试验时，水阻力和惯性将严重影响试验的正确性，而且利用液体浮力配平重力不具有实时操作性，可控性差，产生的误差在长时间操作过程中会导致被实验物体发生竖直向上的运动.为了解决水浮法的这些问题，电磁力系统结合水浮法悬浮微重力环境实验模拟方法利用混合磁悬浮技术在线调控微重力状态和稳定实验物体的高度，但是该方法除了需要液体密封外，还需要在机械臂内布置永磁铁和精确的电磁力补偿控制，无疑引入了附加惯量和增加了系统复杂性［14］.气浮法采用空气轴承的喷气反作用力来抵消支撑在光滑气浮台上空间机械臂的重力.气浮法在二维平面的微重力模拟应用广泛，但是难以用于三维空间运动下的空间机械臂微重力模拟［15］.虽然气浮法具有结构简单、承载能力大、精度高等特点，但是气浮装置由于引入系统的附加惯量较大会导致气浮装置和机械臂的动力学耦合特性明显.从微重力模拟精度上看，气浮法最优，但是从与机械臂动力学特性的耦合和三维微重力模拟能力方面看，气浮法存在明显的不足.

针对现有微重力模拟方法的不足，本文基于磁悬浮气足的吸附特性和通气低阻尼特性提出一种新型的微重力模拟装置—磁气悬吊微重力模拟装置，该装置由磁悬浮钢板、磁悬浮气足和吊绳组成，通过磁吸力、磁悬浮气足重力、气浮力和吊索张力的平衡实现微重力模拟.为了研究所提装置的微重力模拟特性，利用多体动力学方法建立系统的动力学模型并进行仿真计算，并与气浮式微重力模拟装置进行了对比，验证所提装置具有与机械臂低耦合动力学特性.

1 磁气悬吊微重力模拟装置设计

1.1 悬吊式磁悬浮气足

悬吊式磁悬浮气足，如图1所示，由进气口、套轴、出气口、磁块、吊环等部件组成，简称为气足.吊环用于连接吊索悬挂机械臂，并提供重力卸载力.气足在磁块与空间磁性钢板之间的磁吸力、气垫作用力和吊索拉力的作用下在空间钢板上随机械臂移动，从而实现机械臂的二维微重力模拟.

图1 磁悬浮气足示意图Fig.1 Schematic diagram of hybrid magnetic air bearing

一种磁气悬吊微重力模拟装置主要包括磁性钢板、磁悬浮气足和吊绳三部分，实现机械臂的微重力模拟的布局如图2所示.特制的马氏体磁性钢板通过支撑或者悬吊等方式置于平台上方.钢板对气足中磁块的向上磁吸力用于重力卸载和约束气足在钢板平面运动.为了降低气足在钢板上运动的摩擦力，系统参考气浮法在气足与钢板之间形成高压气膜，使气足呈悬浮状态.柔性吊绳一端固连在气足的吊环上，另一端捆绑在机械臂的吊点位置.通过调节磁块的磁吸力、高压气体的压强和绳索长度设计绳索的张力，保证在初始状态机械臂在竖直方向平衡.气足处于工作状态时，能沿着钢板的表面做低摩擦被动跟随式移动.

相比于气浮法，磁气悬吊微重力模拟装置与机械臂之间通过柔性索连接，降低了装置与机械臂之间的动力学耦合特性.又由于系统采用悬吊形式，若将吊环替换为滑轮，并在索的另一端连接配重，可将其拓展到三维微重力模拟.相比于吊丝配重法，系统结构简单易实现.另外，改变磁悬浮钢板的尺寸规格、气足数量、吊点布局、悬吊方式等，能够使得磁气悬吊微重力模拟装置在不同情况下高精度和高可靠性地进行空间机械臂的地面微重力模拟实验.

1.2 臂杆-磁气悬吊微重力模拟系统

为了分析磁气悬吊微重力模拟装置对空间机械臂的重力卸载能力和对其动力学特性的影响，本文以柔性臂杆为研究对象，研究臂杆绕固定点在平面内以恒定角速度ω=π/10 rad/s转动的动力学特性.臂杆-磁气悬吊微重力模拟系统，如图2所示，臂杆由两端和中部三个刚性连接段（编号0，1，2）和两个柔性段（编号3，4）组成，采用两点磁气悬吊形式.系统的具体结构和材料参数如表1所列.

表1 系统的主要参数Table 1 Main parameters of the system

图2 臂杆-磁气悬吊微重力模拟系统Fig.2 The system of a manipulator and the hybrid magnetic air suspension

2 系统多体动力学建模

基于多体动力学方法，建立臂杆-磁气悬吊微重力模拟系统的动力学模型，如图3a所示.臂杆的连接段、磁悬浮气足等不需要考虑变形或者变形很小的，但是在空间大范围运动和转动的物体利用基于旋转向量的刚体建模.考虑臂杆的柔性段和吊索在空间的大幅运动和柔性振动，分别利用基于转动向量的几何精确Timoshenko梁单元和Lagrange索单元建模.不考虑钢板与磁悬浮气足之间的磁吸力和高压气膜，利用平面约束将磁悬浮气足约束在钢板表面上.另外，臂杆的刚性段和柔性段连接、吊索与臂杆刚性段的连接以及吊索与吊环的连接采用固定约束建模，臂杆的一端施加转动约束实现平面转动.磁悬浮气足在钢板上运动时还受到磁阻尼力和气浮阻力的作用，在系统的多体动力学方程中以广义外力的形式加入.因此，臂杆-磁气悬吊微重力模拟系统是典型的刚柔耦合多体动力学系统.

鉴于刚体、几何精确梁、旋转副、固定副等建模理论与方法已经完备［16-18］，下面对系统中的刚体、柔性体、约束和外力的建模方法进行简要介绍.

2.1 刚体

如图3c所示，在刚体r的质心位置固连局部坐标系orxryrzr，选择刚体的广义坐标为

图3 系统多体动力学模型Fig.3 Multi body dynamic model of system

刚体r在全局坐标系下的速度为，局部角速度向量为，其中H是传递矩阵

则刚体的动能可表示为

其中质量矩阵Mr为

其中，mr为刚体质量，Jr为刚体在局部坐标系下的惯性矩阵.具体建模方法见参考文献［19］.

2.2 柔性体

如图3c所示，基于转动向量的两节点Timoshenko梁的广义坐标为

为形函数，l为单元长度.对式（8）求导，可得梁单元的动能为

其中ρ，A和J分别为梁单元的密度、截面面积和截面惯性矩阵，矩阵H的形式与式（4）相同.梁单元的应变向量γ和曲率向量κ分别为

其中，上标一撇表示对弧长的求导，矩阵A的形式与式（2）相同.则Timoshenko梁的弹性势能为

其中，CN和CM为线弹性本构关系.索单元的建模方法与梁单元类似，只需将梁单元的动能、弹性力和阻尼力消除转动项即可.具体建模见文献［20］.

2.3 约束

如图3d所示，系统包含固定约束、平面约束和转动约束，分别建立这三类约束的约束方程为

其中，rk=[xkykzk]T和[xkykzk]T分别为k(k=I，J)在刚体约束对应位置或梁/索节点的位置向量和旋转矩阵A分量.索节点与刚体的固定约束只有位置约束.另外，在臂杆的一端施加转动角速度约束

其中，ωz为刚体0的角速度的分量.

2.4 磁气阻尼力

磁悬浮气足内部对称分布四个圆柱形永磁铁，由周恩权［21］和谢晓［22］等针对圆柱形永磁铁磁场建模和实验可得，磁铁在运动过程中将受到与运动速度方向相反的磁阻尼力fmag的作用，其大小与运动速度相关

其中，k可由实验测得，本文中的k为0.18 kg/s.除了磁阻尼外，磁悬浮气足在运动过程中还受到气浮阻尼fair的作用，其方向与气足运动方向相反，大小同样由实验测得

2.5 系统控制方程

利用第一类拉格朗日方程，将系统中的刚体、柔性体整合，可建立描述多体系统的微分代数方程Differential-Algebraic Equations，DAEs）［23］：

其中，T和U为系统的总动能和总势能，q为系统广义坐标向量，Q为系统的广义外力向量，Φ为系统约束方程，λ为拉格朗日乘子向量.完成系统动力学方程建模后，采用向后差分方法高效求解［24］.

3 系统的动力学特性

为了验证磁气悬吊微重力模拟装置对臂杆的动力学特性的影响，利用上一节的建模方法分别对臂杆-气浮微重力模拟系统（如图3b所示）和零重力臂杆系统建立多体动力学模型，其中气浮法中两个气足和气足支撑杆的总质量均为0.422 kg，采用刚体建模，与臂杆刚性段（1，2）的连接采用固定约束建模.在相同的求解框架下，分别对臂杆-磁气悬吊微重力模拟系统、臂杆-气浮微重力模拟系统和零重力臂杆的多体动力学模型进行仿真计算.

臂杆在电机驱动下在平面oxy内转动，对比臂杆刚性段1，2的水平方向位移如图4所示，磁气悬吊微重力模拟法下的位移与零重力几乎吻合，而气浮法出现2倍转动周期的较大振动，且末端位置2处振动更明显.将图4中的位移减去臂杆的刚性运动得到臂杆的柔性振动，如图5所示.磁气悬吊法和零重力下臂杆柔性振动较小，而气浮法柔性振动较大，且出现2倍转动周期的较大振动.

图4 连接段1和2的平面内位移Fig.4 The in-plane displacements of connecting segments 1 and 2

图5 连接段1和2的平面内柔性变形Fig.5 The in-plane flexible deformations of connecting segments 1 and 2

对臂杆柔性振动结果进行FFT变换法得到臂杆频谱图如图6所示，磁气悬吊法和零重力法的臂杆柔性振动频率约为2 Hz，而气浮法下臂杆振动频率约为1 Hz，这是因为气浮法的气足与臂杆采用刚性连接，增加了臂杆的附加惯量，臂杆振动频率变小.因此，相对于气浮法，磁气悬吊装置和臂杆的耦合特性较弱.臂杆的一阶弯曲固有频率和振动频率如表2所列.图7为微重力模拟装置对臂杆的水平作用力，同样可得到上述结论.

表2 臂杆一阶柔性振动频率Table 2 Main parameters of system modeling

图6 连接段1和2的平面内柔性变形的频谱Fig.6 The FFT spectrum of the in-plane flexible deformations of con‑necting segments 1 and 2

图7 微重力模拟装置对机械臂的水平力Fig.7 Horizontal force on connecting segments by microgravity simulator

图8为对臂杆1和2位置处竖直方向位移和竖直方向作用力.磁气悬吊法由于绳索柔性的影响，臂杆的竖直方向位移和竖直卸载力在额定值附近存在微小波动，影响卸载能力，但是本方案可拓展到三维空间微重力模拟运动.

图8 连接段的竖直方向位移和受力Fig.8 Vertical displacement and force of connecting segments

4 结论

针对现有的微重力模拟装置和机械臂的耦合动力学特性明显，提出一种新型的磁气悬吊微重力模拟装置.本文围绕微重力模拟装置对臂杆动力学特性影响展开研究，针对臂杆-磁气悬吊微重力模拟系统、臂杆-气浮微重力模拟系统和零重力臂杆系统动力学建模，微重力模拟装置对臂杆动力学特性的影响进行深入研究.结果表明，相对于气浮微重力模拟装置，磁气悬吊微重力模拟装置对臂杆的动力学特性影响较小，这对实现机械臂在空间运动的操作平稳性和定位精度都有重要意义.本装置还可通过加滑轮的方式实现机械臂在三维空间的微重力模拟.