基于ls-dyna软件的升降水气过程数值模拟.docxVIP

基于ls-dyna软件的升降水气过程数值模拟 0 起落架系统的数值仿真研究是未来市场 涂层系统是一种高效可靠的恢复装置，在中国的载人宇宙飞船、返回卫星的回收和美国的宇宙勘探中得到了广泛应用。充气是降落伞整个工作过程中最为复杂的一个阶段,也是研究人员迫切地想通过理论分析来解释的。目前,对于降落伞充气性能(除最大开伞动载外)的了解或掌握主要是通过一些试验结果的统计分析来获取的,充气性能与开伞条件之间的关系也是通过试验结果统计、拟合出来的。这使得充气时间、伞衣投影面积变化的一些经验公式存在着一定的局限性,只能在某些情况下适用。 随着“嫦娥一号”的成功发射,我国对于火星探测的研究即将进入实质阶段。火星大气层的成分、物理性质与地球的大气层存在着较大的差别,因此,对于降落伞系统在火星大气环境下的工作性能评定,不能像地球上回收航天器所用的降落伞系统那样可以通过空投试验来直接验证,而需要通过理论分析或者数值仿真来进行。目前回收系统工作过程的数值仿真系统所用到降落伞开伞过程的一些参数,如充气时间、伞衣投影面积的变化一般是采用经验公式来处理的。有关研究表明,对于火星这样稀薄的大气环境,这些常用的经验关系式将不再适用。因此,如何通过数值仿真来获取研究降落伞充气过程性能的主要参数,已经被提上了降落伞理论研究和工程实践的日程。 进入21世纪以来,应用LS-DYNA有限元分析软件来分析降落伞的工作过程(充气阶段、稳定阶段等)成为一个热点。国外在这方面做了大量的工作,但是对降落伞充气过程的仿真还不能完全地解决。国内使用LS-DYNA软件对降落伞充气过程仿真研究尚属空白。本文基于LS-DYNA软件对数值模拟降落伞充气过程的方法进行了初步研究。 1 起落架性能表征 降落伞的充气过程是指从伞系统全长地拉直到伞衣第一次充满(伞衣投影直径第一次能达到稳定下降所需的伞衣投影直径,不考虑过度充气)为止的整个工作过程,如图1所示。 降落伞的充气性能主要由临界开伞速度、伞衣充满条件、充气时间、充气距离、伞衣投影面积变化和最大开伞动载等参数来描述和表征。本文基于LS-DYNA软件对降落伞的充气时间和伞衣投影面积变化的数值模拟分析,来研究降落伞的充气过程。 2 基于流固耦合的数值模拟 LS-DYNA软件是著名的显式动力分析程序,具有强大精确的有限元仿真性能,拥有大量不同种类的单元模型、材料模型和算法选择,能够很方便地处理各种高度非线性问题。到目前为止,它是对降落伞充气过程进行仿真的最优工具之一。LS-DYNA 971是一个主要的程序版本。它增加了透气性材料的流固耦合算法,这样就可以通过控制流固耦合过程的各个参数,真实地模拟降落伞的充气过程。 LS-DYNA与一般的CAE辅助分析程序操作过程相似,LS-DYNA的一个完整的显式动力分析过程包括前处理、求解以及后处理3个基本操作环节,如图2所示。 在解决流固耦合问题时,需要一种将Lagrange和Euler的优点结合起来的方法,即Arbitrary Lagrange-Euler,简称ALE方法。LS-DYNA程序中采用关键字*CONSTRAINED_LAGRANGE_IN_SOLID实现固体与流体间的耦合。采用流固耦合算法模拟问题时,往往要对Lagrange算法中的固体结构进行约束,将固体结构的相关参数传递给流体单元。在LS-DYNA中数值模拟降落伞的充气过程时,固体结构的主要约束方法是罚函数耦合约束。对于罚函数约束而言,其原理是耦合系数追踪Lagrange节点(结构,即从物质)和Euler流体物质(主物质)位置间的相对位移d,如图3所示。 3 模拟模型ls-dyna的构建 3.1 模拟模型的构建 3.1.1 结构模型的建立 对降落伞仿真模型作如下假设: 1) 仿真充气过程前模型为轴对称的结构; 2) 模型的初始形状为有一定的进气孔尺寸,截面为梅花形状的模型; 3) 模型只考虑伞绳的连接,对于伞衣径向带、顶孔绳等不作具体处理。 3.1.2 基本参数设置 1 c-9平面圆形伞结构参数 降落伞仿真模型选用美国的经典伞型—C-9平面圆形伞,如图4(a)所示,具体结构参数如表1所示。 完整的伞衣和伞绳几何模型如图4(b)所示。 2 材料的密度和密度 降落伞仿真模型的伞衣材料选用美军标的MIL-C-7020 III材料,假设该材料为各向同性材料,密度为533.77 kg/m3,弹性模量为0.4309 GPa,泊松比为0.14,伞衣厚度为0.000 1 m。伞绳材料的密度为462.00 kg/m3,弹性模量为97.0 GPa。假设降落伞伞绳的横截面为圆形,面积为4.91×10-6m2。 3 有限单元网格划分 LS-DYNA的前处理软件用于对降落伞仿真模型进行有限单元的网格划分。伞衣单元选用薄壳单元,伞绳单元选用离散的梁/索单元,得到降落伞的有