基于伞载多体动力学的降落伞开伞过程研究_贾华明.pdf

基于伞载多体动力学的降落伞开伞过程研究贾华明，王文强，刘乃彬，林汝领，李少腾，（北京空间机电研究所，北京；中国航天科技集团有限公司航天进入、减速与着陆技术实验室，北京）摘要：针对降落伞开伞过程中柔性大变形、强非线性和强耦合性的特点，采用伞载多体动力学方法，以充气时间为自变量，分别建立拉直、充气和稳降 个阶段的动力学模型。对降落伞开伞过程中物伞系统运动学和动力学进行数值分析，研究降落伞开伞过程，通过某型号降落伞飞行试验验证，表明：基于多体动力学的降落伞开伞过程计算简单方便，能够快速指导降落伞的设计迭代，计算结果可以很好地体现降落伞的开伞特性。关键词：降落伞；多体动力学；充气时间；开伞过程；运动学；动力学中图分类号：.文献标识码：文章编号：（）收稿日期：；修回日期：基金项目：国防科技创新特区项目（）第一作者：贾华明，男，硕士，高级工程师，研究方向为航天器回收系统降落伞设计及高性能柔性材料。：，（，；，）：、，、，：；引言 降落伞作为一种典型的气动减速装置，在航空航天领域应用广泛。随着中国载人航天、火星探测等重大专项工程的持续推进，对降落伞的减速性能提出了更高的要求，同时也给降落伞的理论研究和工程应用带来了新的挑战。目前降落伞理论研究上的不足造成设计分析中某些缺陷无法避免，难点主要集中在降落伞开伞过程研究。降落伞从拉直到充气，伞衣会经历急剧的结构大变形，而且伞衣内外的气流流动状况非常复杂，材料属性和外界环境对降落伞第 卷 第 期 年 月 载 人 航 天 DOI:10.16329/ki.zrht.2023.01.007的开伞过程都有重要影响，存在多种因素剧烈耦合、学科交叉的鲜明特点，这是降落伞开伞过程研究的重点和难点。目前降落伞开伞过程的研究方法有流固耦合力学方法（，）、计算流体力学方法（，）和伞载多体动力学方法（，）。流固耦合力学方法是目前研究的热点，它同时考虑伞衣变形和周围流场的变化，可得到流场、结构、运动等多场参数，但模型建立复杂，计算消耗量巨大，依然存在诸多技术难题，其中解决负体积问题，提高计算效率是当前需要重点关注的问题，一般用于短时的瞬态分析。计算流体力学方法能够得到内外流场结构细节，计算消耗量中等，但不考虑伞衣变形，一般用于稳降状态下的气动性能和静稳定性能研究，其中外形、透气、是否考虑尾流是该方法分析问题的关键点。伞载多体动力学方法模型建立较简单，计算消耗量最小，适用于长时间的运动计算。由于流固耦合方法和流体力学方法建模复杂、计算量极大，尤其是大型伞系统，不能够满足快速设计迭代、性能分析和预测评估的工程需要。本文基于伞载多体动力学，以充气时间为自变量，将轴向动量方程作为参考模型，进行了降落伞开伞过程研究，并将计算结果与飞行试验数据进行了对比，用以验证降落伞开伞特性。典型的降落伞开伞过程 如图 所示，降落伞的开伞过程包括拉直、充气和稳降 个阶段。一般情况下，包装于伞包内的主伞是由引导伞将其拉出并拉直。目前，航天器回收领域降落伞基本采用倒拉法，即先拉伞绳的拉直方法，这是因为在相同条件下，倒拉法的拉直力会小于顺拉法。在倒拉法中首先由引导伞将伞包提起，使其与载荷分离，然后从伞包中先拉出伞绳，再逐步拉出伞衣。从伞系统全长被拉直，进入充气阶段，空气团从伞衣底边进入伞衣内部，引起伞衣迅速充满，在工程上，为了降低开伞载荷，通常采用底边收口方式，这样在充气过程中就增加了收口工作阶段，由火工切割器切断收口绳后，伞衣底边快速张开，伞衣底边进气口迅速增大，直至伞衣完全充满，之后载荷在降落伞的作用下稳定下降。图 降落伞开伞过程 动力学模型.拉直阶段.物伞系统动力学方程 整个伞系统的拉直过程是引导伞和载荷二者相对运动的结果，从理论上看，是 个变质量体的相对运动。图 所示为倒拉法开伞拉直过程的动力学模型，图中 为拉直中的主伞系统微元质量，为载荷连同已拉直主伞系统的质量；为引导伞、伞包连同未拉出的主伞系统的质量；为伞绳拉出阻力，为拉直力，为航迹倾斜角，为伞系拉直距离。图 倒拉法开伞拉直过程 为简化模型，做如下假设：）拉直过程中，载荷和引导伞每一瞬时的运动遵循同一条轨迹，航迹倾斜角 按载荷计算；）物伞系统没有升力，在二维平面内运动；第 期 贾华明，等 基于伞载多体动力学的降落伞开伞过程研究）在拉直过程中，不考虑伞绳伸长；）把引导伞、载荷和拉直中的主伞系统微元质量 作为 个质点处理。拉直过程中运动方程如公式（）所示：（）|（）式中，为载荷的气动阻力，为主伞系统已拉出部分的气动阻力，为引导伞的气动阻力，为主伞包的气动阻力，为拉动的伞系统单位长度的质量，为载荷速度，为引导伞速度，为时间，为重力加速度。拉直阶段的轨迹如公式（）所示：|（）式中，为载荷的水平位移，为载荷的竖直位移。拉直力计算公式如公式（）所示：（）（）为使计算结果接近实际情形，动力学模型中加入了参考风场数据模型，其中风向定义为逆风速矢量（即风的来向）与正北方向的夹角，以顺时针方向为正，并进行三次样条插值，获得任意高度上的风速和风向。.伞系质量变化规律 在拉直过程中，已拉出主伞系统质量 和未拉出主伞系统质量 是随时间而变化的，倒拉法降落伞质量变化规律如公式（）、（）所示。，（）（）|，|（）（）式中，为伞绳总质量，为伞衣总质量，为伞衣底边沿伞绳方向的单位长度质量，为伞绳单位长度质量，为伞衣底边宽度，为伞绳长度，为伞系统全长，为伞衣名义直径。.拉直阶段计算 根据式（）、（）的 个微分方程组以及式（）、（）的伞系质量变化规律，以时间 为计算步长，可以计算出各个时间 内的各种运动学参数、。根据计算得出的各个时间 内的、值，代入公式（）中可求出各个时间 内的伞绳张力，而最大的张力即为拉直力。.充气阶段.伞载组合体动力学方程 伞载组合体充气阶段如图 所示，其中，是降落伞阻力，是降落伞重力，是开伞载荷，是载荷阻力，为载荷重力，为伞载组合体航迹倾斜角。图 降落伞充气过程 为简化模型，做如下假设：）不考虑伞载组合体伸长；）不考虑伞载组合体的升力，在二维平面内运动；）伞载组合体是两质点运动，载荷的质量 集中在载荷重心位置，伞的质量 集中在伞衣底边中心位置。充气阶段，伞的质心相对底边保持位置不变；）伞载组合体在充气阶段轴线呈直线。充气过程中运动学方程如公式（）所示：载 人 航 天第 卷（）（）（）（）（）|（）其中，是载荷质量，是降落伞质量，是伞附加质量，是伞附加质量变化率，是载荷水平速度，是载荷竖直速度，是时间，是重力加速度。开伞载荷算计算如公式（）所示：（）.降落伞单级收口阻力面积 图 是降落伞单级收口阻力面积示意图，（）是收口状态下阻力面积，（）是全张满状态下阻力面积，、为充气阶段的 个过程，是收口过程，其中。图 充气阶段阻力面积变化规律 降落伞单级收口充气阶段分为 个过程：过程是伞开始充气到伞衣呈乌贼状外形，这个过程阻力面积与时间为线性；过程阻力面积维持恒定值（）；过程是结束收口至伞全张满，这时候阻力面积维持在另一个恒定值（），这个过程中阻力面积是时间的 次方，用公式（）表示充气阶段阻力面积变化。（）（），（），（）（）（）|，|（）其中，（）是阻力面积；是充气参数，一般在 之间。公式（）为 过程中阻力面积与时间呈现的线性关系。（）.（）其中，是设定参数，是伞开始充气速度。公式（）为 过程中阻力面积与时间的关系。（）.（）.（）其中，是设定参数，是收口完成时伞载组合体速度。.充气阶段附加质量变化 充气阶段会受到降落伞附加质量的影响，附加质量由降落伞特征和大气参数等决定。公式（）为降落伞附加质量的计算公式。（）.（）其中，是设定参数，是空气密度。降落伞充气阶段附加质量变化率 如公式（）所示：.（）（），.（）（）（）()（）（）（）（），|（）.充气阶段计算 使用公式（）（）的伞衣阻力面积变化和公式（）、（）的附加质量变化，按照公式（）、（）能够计算降落伞单级收口充满阶段所有运动参数和开伞动载。.稳定阶段 在伞衣充满后，物伞系统在气动阻力的作用下继续减速，直至系统达到稳定下降。稳定下降阶段，伞衣已经充满后，降落伞的阻力面积和附加质量不再发生变化，即（），。第 期 贾华明，等 基于伞载多体动力学的降落伞开伞过程研究继续使用公式（）、（）、（）就可以完成稳定下降阶段计算。.开伞过程计算 以时间 为自变量，采用 法，可以进行开伞过程拉直阶段、充气阶段和稳定阶段的物伞系统全部运动学和动力学计算，充气阶段用到的 个设定参数、，可通过已获取的试验数据进行归纳整理得出。某十字型降落伞充气过程动力学模型的 个设定参数分别取为，.，本文建立的模型增加了多种伞型，如环帆伞、带条伞、平面圆伞和盘缝带伞，这大大提高了动力学模型的适用范围。设定参数 和 与伞型密切相关，美国的标准环帆伞.，.；美国的改型环帆伞 .，.；中国的某环帆伞.，.。设定参数 与伞型和透气量密切相关，美国阿波罗环帆伞.；中国某平面圆伞.。验证与分析 某航天器返回舱采用了带条伞环帆伞降落伞方案，飞行试验后通过对数据的分析整理，分别拟合出带条伞和环帆伞 个参数、的取值，将其固化在建立的开伞过程的动力学模型中，进行编程计算，根据多次飞行试验数据对模型中的一些参数进行反复修正优化，计算结果精度，尤其是开伞动载得到了很大提高，更加接近于真实情况。利用参数固化后的动力学模型对某型号降落伞某次飞行试验开伞过程进行计算，降落伞特性参数见表 所示，开伞工况见表 所示。计算结果与试验数据的对比曲线见图 所示，其中图 为弹道轨迹的对比曲线，图 为高度变化对比曲线，图 为速度变化对比曲线，图 为航迹倾斜角变化对比曲线。图 中，将开伞时刻的载荷经纬度与高度置为零，减速伞工作阶段轨迹较陡，试验数据曲线与计算数据曲线一致性较好。主伞工作阶段，轨迹变得较平滑，试验数据与计算数据出现明显偏差，这是因为仿真计算使用的是参考大气数据，与试验场当天真实风场数据偏差较大，主伞工作阶段的载荷速度已经降到很低，风的影响会产生较大的水平位移，从而造成了试验数据与计算数据的偏差。表 降落伞参数 参数名称数值减速伞伞型带条伞减速伞收口时间 减速伞收口阻力面积.减速伞全展开阻力面积.减速伞工作时间 主伞伞型环帆伞主伞收口时间 主伞收口阻力面积 主伞全展开阻力面积 表 计算工况 项目名称工况值载荷开伞海拔高速 回收速度（）航迹倾斜角（）回收质量 图 弹道轨迹对比 图 高度对比 在图、图 的高度和速度的对比曲线中，试验与计算数据的变化趋势一致性很好。计算曲线载 人 航 天第 卷图 速度对比 图 航迹倾斜角对比 较为平滑，而试验数据曲线有一定的峰值变化，这主要是由于大气参数偏差和气象条件的影响。在图 航迹倾斜角对比曲线中，试验数据曲线和计算数据曲线存在一些差异，减速伞及主伞收口阶段，二者的差异不大，在范围内波动，主伞解除收口后，载荷处于稳定下降阶段，计算数据在参考大气数据的情况下，航迹倾斜角保持在，而试验数据在实际风场的情况下，载荷随风飘动，航迹倾斜角在左右波动。开伞动载是降落伞开伞过程中的一个重要性能指标参数，是降落伞结构设计的依据，也是回收系统接口设计的输入。表 分别对减速伞工作阶段和主伞工作阶段的开伞动载进行了试验数据和计算数据的比较，计算偏差均小于，计算结果很好地体现了开伞特性，与试验结果的一致性较好。其中减速伞工作段的开伞载荷偏差要高于主伞工作段，最大为减速伞收口段，偏差率为.。最小为主伞全展开段，偏差率为.。减速伞和主伞分别采用了带条伞和环帆伞 种不同的伞型，所以计算过程中的、也分别表 开伞动载对比 阶段开伞力试验值 计算值 偏差率 减速伞收口段.减速伞全展开段.主伞收口段.主伞全展开段.采用不同的取值。结论 本文采用伞载多体动力学方法对某型号降落伞的开伞过程进行了计算，并将计算数据与试验数据进行了对比分析，得出如下结论：）减速伞工作段，试验数据和计算数据的弹道和航迹倾斜角一致性较好，主伞工作段，尤其是主伞全展开后，载荷速度明显降低，由于参考大气数据与真实风场数据偏差较大，造成试验数据和计算数据产生明显偏差，尤其是航迹倾斜角。）在速度和高度的对比曲线中，试验数据和计算数据的变化趋势一致性较好。）