动力翼伞系统拟坐标形式的多体动力学建模

动力翼伞系统拟坐标形式的多体动力学建模张青斌; 高峰; 郭锐; 丰志伟; 葛健全【期刊名称】《《兵工学报》》【年(卷),期】2019(040)009【总页数】8页(P1935-1942)【关键词】翼伞系统; 多体动力学; 拟坐标; 飞行仿真; 18自由度【作者】张青斌; 高峰; 郭锐; 丰志伟; 葛健全【作者单位】国防科技大学空天科学学院湖南长沙410073; 96796部队吉林白山135200【正文语种】中文【中图分类】V211.30 引言翼伞系统不仅具有折叠体积小、质量轻的优点，而且拥有优良的操稳和滑翔性能，可以在一定初始偏差和风场等干扰条件下实现定点无损着陆，广泛应用于高空投放、远距离渗透和精确着陆等任务中，具有极大的军事和民用价值。

目前国外研发了一系列颇具代表性的空投系统，例如美国Atair Aerospace 公司研制的自主导航空投系统、美国Strong Enterprises公司研制的Screamer空投系统、美国Alabama大学和海军研究院联合研发的Snowflake空投系统，以及加拿大米斯特机动综合系统技术公司研制的Sherpa系统和德国Cassidian公司研制的Paralander空投系统等[1]。

上述空投系统能完成从数千克到数十吨的货物精确空投，落点偏差甚至在几十米以内。

翼伞空投系统已成为世界各军事强国空投技术的发展方向，但是由于理论研究难度大且技术复杂，目前世界上掌握精确空投核心技术并有产品正式列装的国家较少。

翼伞系统飞行动力学是开展翼伞系统设计研制的理论基础。

文献[2-3]建立了多种自由度的翼伞多体动力学模型，分析了展弦比、后缘下拉方式、安装角及系统长度对雀降性能的影响。

Goodrick[4]提出并建立了翼伞和回收物刚性连接的6自由度模型，研究了翼伞系统纵向运动的稳定性和飞行性能。

Zhang等[5]在建立翼伞6自由度模型基础上，通过实验验证了6自由度模型的可行性，用该模型仿真得到了翼伞各项参数变化对翼伞飞行状态的影响。

Slegers[6-7]建立了9自由度多体动力学模型，分析了安装角、伞衣弯曲度及后缘下偏量对翼伞侧向运动的影响。

蒋万松等[8]采用Lagrange乘子法建立了两体8自由度动力学仿真模型，并通过空投试验验证了仿真模型的有效性。

Meirovitch[9]利用分布系统的扩展哈密顿原理，导出混合系统的标准Lagrange方程，随后以拟坐标表示刚体运动，将混合系统的2阶Lagrange方程转化为一组适用于控制的状态方程。

陈建平等[10]通过引入Lagrange乘子，获得了大型翼伞系统封闭的飞行动力学方程组，对某大型翼伞系统飞行动力学过程进行了数值仿真。

杨海鹏等[11]假定伞绳和吊带分别通过铰链与翼伞和载荷相连，考虑了伞绳和吊带的弹性以及伞绳相交点与伞衣的相对运动，通过建模仿真得到了吊带、伞绳交点运动的情况。

但是在目前文献中，还很少有关于载荷连接模式和稳定性的研究成果。

为此，本文借鉴航天器多柔体动力学建模，以某动力的弹载翼伞系统为工程背景(见图1)，考虑伞绳、载荷吊带的连接方式，针对翼伞稳定飞行的运动特性，将伞衣和载荷简化为刚体，伞绳和吊带交汇点处理为质点；分析翼伞系统各体的受力特性，给出了翼伞多项式形式的气动表达式；将各体平动速度、角速度作为18自由度拟坐标，采用Lagrange拟速度形式的建模原理[12]，建立了18自由度的动力学模型。

利用空投试验证明了模型的有效性，仿真分析下拉后缘、突风扰动和不同吊挂模式下的相对运动特性，对比了6自由度单刚体和18自由度多体动力学模型的仿真结果。

图1 某小型翼弹载系统Fig.1 A small wing projectile-borne parafoil system 1 模型描述如图2所示，翼伞系统4体模型主要包括伞衣、伞绳、吊带和载荷部分。

本文将稳定飞行状态下的伞衣、载荷以及2个连接点作为4个刚体，伞衣和载荷分别记为刚体A和刚体B；伞绳和吊带的交汇点简化为质点，记为L和R；伞绳简化为4根绳索LA1、RA2、RA3和LA4，吊带简化为4根绳索LB1、RB2、RB3和LB4. 上述多刚体模型中刚体A和B分别有6个自由度，质点L和质点R分别有3个自由度，总共18个自由度。

该模型可以较完整地描述伞衣与载荷的相对运动，反映了实际的空投工况和载荷吊挂模式。

图2 翼伞系统4体模型Fig.2 4-body model of parafoil system为了便于进行运动描述和受力分析，引入如下4个坐标系：惯性坐标系OIxIyIzI；伞衣坐标系OPxPyPzP；载荷坐标系OSxSySzS；刚体A体坐标系OAxAyAzA，原点为伞衣和伞绳整体的质心；刚体B体坐标系OBxByBzB，原点为载荷和吊带整体的质心；气流坐标系OQxQyQzQ(见图3)。

图3中，ΓP和ΓS分别为伞衣和载荷的安装角，α为伞衣攻角，v∞为来流速度。

刚体A体坐标系相对惯性坐标系的滚转角、俯仰角和偏航角记为φA、θA和ψA，相应的转换矩阵记为刚体B体坐标系相对惯性坐标系的滚转角、俯仰角和偏航角记为φB、θB和ψB，相应的转换矩阵记为伞衣侧滑角记为β，从伞衣坐标系到气流坐标系的转换矩阵记为则从载荷坐标系到刚体B体坐标系的转换矩阵为从伞衣坐标系到刚体A体坐标系的转换矩阵为侧向图如图4所示。

图3 翼伞系统4体模型坐标系示意图Fig.3 Schematic diagram of 4-bodymodel coordinate system of parafoil system图4 侧面示意图Fig.4 Side view of parafoil system2 翼伞系统受力分析当翼伞系统在其稳定状态飞行时，主要受力包括重力、气动力和螺旋桨推力等。

伞衣具有柔性大的变形特性，而且伞绳和载荷的连接吊带均为柔性绳索，很难获得其翼伞系统精确的气动特性。

为此，本文不考虑翼伞下拉后缘过程中的迟滞效应，并忽略各部件的气动耦合作用，采用工程简化模型表示各体的气动力。

2.1 伞衣受力分析关于伞衣气动模型的建立，刘媛等[13]通过有限体积元法求解Navier-Stokes湍流模型的控制方程，并对美国国家航空航天局兰利研究中心的三维翼伞进行了气动性能分析，验证了该方法的可行性；进而利用该方法分析得到了不同迎角来流条件以及伞衣后缘不同程度下偏的气动力数据。

汪龙芳等[14]对伞衣透气性进行了二维和三维数值模拟，得到了透气性和迎角对翼伞升阻力系数变化的影响。

为了便于后续归航控制设计，本文在计算流体力学基础上，采用多项式表示伞衣质心处的气动力和气动力矩[1]，即(1)(2)式中：ρ为大气密度；v为伞衣质心相对于气流的速度大小；SP为伞衣面积；δS 和δa为对称舵偏和反对称舵偏；p、q和r分别为伞衣滚转角速度、俯仰角速度和偏航角速度；b为伞衣展长；c为伞衣弦长；其他参数k、CD0、CL、CDδS、CDδa、CL0、CLα、CLδS、CLδa、CYβ、CYp、CYr、CYδa、Clβ、Clp、Clr、Cl、Clδa、Cm0、Cmα、Cmq、CmδS、Cmδa、Cnβ、Cnp、Cnr、Cnδa均为气动力常系数，由试验得到。

由于伞衣充气张满后，其平均密度与空气密度相当，必须考虑附加质量的影响。

本文基于文献[15-16]，将伞衣附加质量等价地处理为伞衣质心的作用力和力矩，并在伞衣坐标系下表示为FPM=-ωP×(MavP)，(3)MPM=-ωP×(JaωP)，(4)式中：ωP和vP为伞衣角速度和质心相对气流速度；Ma和Ja分别为附加质量相对伞衣质心的质量矩阵和惯量矩阵，即(5)(6)ma,ii和Ia,ii(i=1，2，3)分别为附加质量相对伞衣质心的平动和转动分量。

作用于刚体A质心的合力，可在刚体A体坐标系下表示为(7)式中：为从惯性坐标系到刚体A体坐标系的转换矩阵；GA为伞衣所受重力；FrtA 和FrqA分别为与伞衣相连的各伞绳张力和气动力对刚体A作用的矢量和。

作用于刚体A的合力矩可在刚体A体坐标系下表示为(8)式中：MrtA和MrqA分别为各伞绳张力和气动力对刚体A作用的力矩矢量和。

2.2 载荷受力分析考虑到坐标变换，作用于刚体B质心处的合外力在其体坐标系下表示为(9)式中：Ft为作用于载荷的螺旋桨推力；GB为刚体B所受重力；FrtB和FrqB分别为与载荷相连的各吊带张力总和与气动力矢量和；FqB为载荷所受气动力，忽略螺旋桨对载荷气动力的影响，则载荷所受气动力为(10)式中：vB和Sd分别为载荷质心速度和阻力面积。

作用于刚体B的合力矩在体坐标系下表示为(11)式中：MrtB和MrqB分别为与载荷相连的各吊带张力和气动力对载荷的作用力矩矢量和，在刚体B体坐标系下分别表示为(12)(13)式中：rOBi为载荷体坐标系下吊带端点Bi的位置矢量，i=1，2，3，4；FrtLB1、FrtLB4、FrtRB2、FrtRB3分别为吊带LB1、LB4、RB2、RB3的张力，方向从载荷接触点指向吊带交汇点；FrqLB1、FrqLB4、FrqRB2、FrqRB3分别为吊带LB1、LB4、RB2、RB3受到的气动力在B1、B4、B2、B3上的集中力。

2.3 伞绳和吊带交汇点的受力分析作用于连接点L的合力可表示为FL=GL+FrtL+FrqL,(14)式中：GL为质点L所受重力；FrtL和FrqL分别为与质点L相连的各伞绳(吊带)张力和气动力。

类似地，作用于连接点R的合力在惯性坐标系下可表示为FR=GR+FrtR+FrqR,(15)式中：GR为质点R所受重力；FrtR和FrqR分别为与质点R相连的各伞绳(吊带)张力和气动力。

另外，在具体仿真计算中将作用在绳索上的气动力等效为在两端点的集中力。

3 基于拟坐标形式的翼伞动力学方程针对图2所示的翼伞系统18自由度模型，将伞衣质心相对惯性坐标系的位移和速度分别记为RA=[xA yA zA]T和姿态角和角速度分别记为载荷质心相对惯性坐标系上的位移和速度记为RB=[xB yB zB]T和姿态角和角速度记为θB=[φB θB ψB]T 和连接点R和L在惯性坐标系下的速度记为vL和vR，在惯性坐标系下位置记为RL和RR，伞衣和载荷的质量记为mA和mB，伞绳和吊带连接点的质量记为mL 和mR. 同时将翼伞系统的拟坐标和拟速度分别写为(16)本文采用拟坐标形式Lagrange方程进行动力学建模，经过复杂推导可得到如(17)式的动力学方程：(17)式中：相关状态变换矩阵(18)(19)广义质量矩阵(20)I3×3为3阶单位矩阵，03×3为3阶0矩阵；广义力的相关表达式为(21)(22)(23)(24)(25)(26)(27)分别为ωA、ωB、vA、vB的反对称矩阵。