十字形伞开伞充气过程数值仿真研究.pdf
《十字形伞开伞充气过程数值仿真研究.pdf》由会员分享,可在线阅读,更多相关《十字形伞开伞充气过程数值仿真研究.pdf(9页珍藏版)》请在咨信网上搜索。
1、航天返回与遥感 第 44 卷 第 3 期 32 SPACECRAFT RECOVERY&REMOTE SENSING 2023 年 6 月 收稿日期:2023-01-06 基金项目:工业和信息化部重点实验室开放基金(KLAECLS-E-202004)引用格式:谢淮,刘宇,王臻,等.十字形伞开伞充气过程数值仿真研究J.航天返回与遥感,2023,44(3):32-40.XIE Huai,LIU Yu,WANG Zhen,et al.Numerical Simulation Study of Cruciform Parachute Deployment and Inflation ProcessJ.
2、Spacecraft Recovery&Remote Sensing,2023,44(3):32-40.(in Chinese)十字形伞开伞充气过程数值仿真研究 谢淮 刘宇 王臻 昌飞(北京空间机电研究所,北京 100094)摘 要 十字形降落伞是航空航天减速领域常用的伞型,开伞充气过程是其工作过程中较为复杂的一个阶段。为了研究十字形伞开伞充气过程中的动力学特性,文章基于任意拉格朗日-欧拉方法对十字形伞开伞过程进行流固耦合仿真计算,并将计算结果与风洞试验结果进行对比。对比发现:仿真计算得到的十字形伞稳态外形与试验结果一致,且仿真获得的降落伞气动载荷历程曲线与试验结果基本相符,载荷呈现出随着开伞
3、过程逐渐增大,在开伞结束达到最大值后略有降低并逐渐稳定的趋势。仿真结果还表明,在十字形伞开伞过程中,伞衣的最大应力点位于伞臂的中心区域,充气顺序的先后会影响伞臂的应力分布,降落伞稳定后伞衣应力呈对称分布。文章采用的仿真方法能较好地模拟十字形伞开伞充气的动力学过程,得到的伞衣应力分布特点及影响因素可为十字形伞的设计与优化提供参考。关键词 十字形降落伞 计算流固耦合动力学 开伞 充气 中图分类号:V445 文献标志码:A 文章编号:1009-8518(2023)03-0032-09 DOI:10.3969/j.issn.1009-8518.2023.03.004 Numerical Simulat
4、ion Study of Cruciform Parachute Deployment and Inflation Process XIE Huai LIU Yu WANG Zhen CHANG Fei(Beijing Institute of Space Mechanics&Electricity,Beijing 100094,China)Abstract The cruciform parachute is a commonly used parachute design in the field of aerospace deceleration,and the deployment a
5、nd inflation stage is a complex stage of its operation.In order to study the dynamics characteristics of cruciform parachute deployment and inflation process,this paper employs the arbitrary Lagrangian-Eulerian method to carry out the fluid-structure interaction simulation calculations on the deploy
6、ment process of a cruciform parachute.The computed results are compared with wind tunnel test results.The comparison shows that the simulated steady-state shape of the parachute is consistent with the experimental results,and the simulated aerodynamic load history curve of the parachute is largely i
7、n agreement with the test results.The load exhibits a gradual increase during the deployment process,reaching its maximum value after the completion of the deployment,followed by a slight decrease and gradual stabilization.The simulation results also show that,during the cruciform parachute deployme
8、nt process,the maximum stress point of the canopy is located in the central region of the canopy arms,and the inflation sequence affects the stress distribution along the arms.After the parachute stabilizes,the canopy stress exhibits a symmetrical distribution.第 3 期 谢淮 等:十字形伞开伞充气过程数值仿真研究 33 The simu
9、lation method employed in this study effectively simulates the dynamic process of cruciform parachute deployment and inflation,and the characteristics of canopy stress distribution and influencing factors can provide valuable insights for the design and optimization of cruciform parachutes.Keywords
10、cruciform parachute;fluid-structure interaction(FSI);deployment;inflation 0 引言 降落伞是目前航空航天领域应用最广泛的减速手段1-2,根据伞衣结构形式的不同,可以分为平面圆形伞、方形伞、带条伞、环帆伞以及十字形降落伞等。相较于其他伞形,十字形降落伞具有工艺简单、稳定性好、开伞动载小等优点,其独特的结构使得它在高空高速降落和复杂的气象条件下表现出色,成为众多型号任务的理想选择3-4,被广泛用于军事空投、紧急救援、探测器着陆以及科学实验等领域。由于十字形伞具有独特的设计特点和广泛的应用背景,对其力学特性进行深入研究具有重要
11、意义5。在降落伞方案设计中,充气过程中的力学特性6,包括充气时间、气动力变化以及伞衣载荷变化是需要考虑的关键因素7。相较于高成本的试验手段8-10,仿真计算方法是目前研究降落伞开伞过程动力学的有效手段11。流固耦合是降落伞开伞过程仿真的核心问题,因为降落伞伞衣材料为柔性结构,同时具有透气性,这使得降落伞开伞过程流固耦合处理及其困难12。目前应用于降落伞开伞计算的流固耦合仿真方法有CFD-MSD(Computational Fluid Dynamics-Mass Spring Damper)方法13、DSD/SST(Deforming-Spatial-Domain/Stabilized Spac
12、e-Time)方法14及浸入边界法(Immersed Boundary Method,IBM)15等,但是这些方法均具有局限性:CFD-MSD 方法侧重于伞衣外形的计算,无法得到结构应力信息;DSD/SST 方法侧重于外部流场的数值计算,难以对充气阶段进行研究;IBM 侧重于固体的响应,对流场的细节计算不够精确。近些年来,商业软件 LS-DYNA 的流固耦合方法可以克服以上方法的缺点,满足降落伞仿真的需要,因此被广泛应用于亚声速降落伞、超声速降落伞等降落伞仿真领域16-19。在 LS-DYNA 的流固耦合(Fluid-Structure Interaction,FSI)方法中,流场单元通过任意
13、拉格朗日-欧拉(Arbitrary Lagrange-Euler,ALE)方法进行求解,降落伞单元通过拉格朗日法进行求解,最后通过罚函数方法对流场及伞衣之间的相互作用力进行耦合计算。本文基于 ALE/FSI方法对十字形伞开伞过程进行仿真计算,并将计算结果与风洞试验结果进行对比,验证仿真模型的正确性,在此基础上分析十字形伞开伞过程中载荷变化以及伞衣内的应力分布变化的特点,从而能够更深入地理解十字形伞开伞过程动力学特性,为其设计和应用提供更准确的指导。1 研究方法 1.1 ALE 方法 ALE 方法是进行流固耦合计算的重要方法,最初是为了克服拉格朗日法和欧拉法的缺点而被提 出20-22,其基本思想
14、是:求解使用的计算网格既不像欧拉法一样固定不动,也不像拉格朗日法一样跟随质点运动,而是可以相对于坐标系做任意运动。基于 ALE 坐标,一个变量相对于时间 t 的导数可以表示为 LAAd(,)(,)()(,)dttttt=+XXvwX(1)34 航 天 返 回 与 遥 感 2023 年第 44 卷 式中 XL为拉格朗日描述下的质点坐标;XA为ALE描述下的质点坐标;v 和 w 分别为质点运动速度和相应的ALE坐标点运动速度。由式(1)可以看出,当 w=v 时,有 XL=XA,公式转变为拉格朗日形式;当 w=0时,右边第一项为0,公式转变为欧拉形式。基于ALE方法的N-S方程可表示为 ()0()(
15、)tteet+=+=+=+vvwvvwvfvwvf v:(2)式中 为密度;f 为体积力;e为比内能;为柯西应力张量,可表示为 12()3p=+Iv I (3)式中 p为压力;I为单位张量;为流体动力粘度;为应变率张量。通过算子分裂的方法23,式(2)可以分为两步求解:拉格朗日步和输运步。首先通过一个显式有限元方法实施拉格朗日步,在该过程中网格点随流体质点运动,即在方程中令w=v。拉格朗日步的方程为 tet=+=+vfvf v:(4)质量方程也可以通过拉格朗日法求解,但是在此可使用更简便和精确的方法 0J=(5)式中 0为初始密度;J为雅克比矩阵表示的体积应变,可表示为 A,L,detijJ=
16、XX(6)式中 XA,i表示矢量的第i个分量,XL,j表示矢量的第j个分量,由于XA和XL均为三维矢量,故i=1,2,3,j=1,2,3。接下来实施输运步,网格进行运动,质量、动量和能量均因网格运动而产生通量。在本文中,考虑网格运动回初始构型位置,即完整的ALE过程完成后网格保持不变,因此输运步的方程可表示为 ()0t+=vw (7)1.2 结构控制方程 降落伞结构为柔性连续介质,其控制方程为 sssst=+uf (8)式中 s为降落伞材料密度;su为降落伞上质点的速度;s为降落伞材料柯西应力张量;sf为降落伞体积力。1.3 罚函数方法 在流固耦合问题中,流体与固体之间的相互作用是通过模拟两者
17、之间的相互影响来实现的。其中一种常用的方法是罚函数方法。在罚函数法中,为了模拟接触力的相互作用,在接触节点处使用一个类似 第 3 期 谢淮 等:十字形伞开伞充气过程数值仿真研究 35 于弹簧的力模型。在该力模型中,假设存在两个接触节点,引入一个假想的“弹簧”,它对两个节点施加的力大小相等,方向相反。这个“弹簧”的作用类似于一个力的传递媒介,使得节点之间的接触力可以相互传递。当一个节点施加力到弹簧上时,弹簧会将等大、相反方向的力传递给另一节点。这种模拟接触力的处理方法实施简单,可以使得节点之间的接触行为在流固耦合问题的模拟中得到有效的考虑,同时保证了计算过程中动量及能量的守恒性。在本文中用罚函数
18、方法来解决两个问题:伞衣表面的相互接触(自接触)和气体与伞衣表面的接触(流固耦合)。1.3.1 伞衣自接触 当未充满状态下伞衣受来流与伞绳共同作用发生变形时,伞衣单元之间会发生相互接触,本文使用单面接触算法模拟这一现象。以二维情况为例,定义穿透距离s1为一个单元节点(Slave)与面(Master)之间的垂直距离,则 1s=l n (9)式中 l为面上一点到节点的矢径;n为面法向量。由s1的定义可知,当s1=0时,节点与面发生接触,当s10时,节点与面发生穿透。当穿透发生时,定义约束力F1为 11 1Fk s=(10)式中 k1为刚度系数。1.3.2 流固耦合 在流固耦合模型中,流体单元为三维
19、ALE实体单元,伞衣单元为二维拉格朗日壳单元,流体单元与伞衣单元之间同样通过罚函数法耦合来传递作用力。当拉格朗日单元与ALE单元之间相对运动时,将拉格朗日单元节点视为从节点,不考虑网格运动的ALE描述与拉格朗日描述的单元节点坐标相等,考虑网格运动后可进行计算得到主节点,定义穿透位移s2为 2ms()dt=svv (11)式中 vm,vs分别为主节点及从节点运动速度。耦合的作用力定义为 22 2k=Fs(12)式中 k2为耦合刚度系数。1.4 伞衣透气性方程 织物透气量对伞衣充气过程影响很大,而在充气过程中随着伞衣的内外压力差变化,织物透气量也产生变化。伞衣的透气速率与伞衣两侧压差的关系可由经验
20、公式表示为24 21q2qpKvKv=+(13)式中 p为织物两面压差;1K、2K为描述织物透气性的常数,仅由织物决定,一般通过试验测量得到;qv为织物透气速度。2 数值计算模型 2.1 研究对象 根据试验任务,本文仿真计算所使用十字形伞外形尺寸如图1所示。降落伞的材料选用织物透气性参数为K1=35.4和K2=5.29106/m的锦丝绸。36 航 天 返 回 与 遥 感 2023 年第 44 卷 (a)俯视图 (b)侧视图(a)Top view (b)Side view 图 1 十字形伞外形示意 Fig.1 The diagram of cruciform parachute shape 为降
21、低模型复杂度,作以下基本假设:不考虑前体的影响,连接绳一端为伞绳集束点,另一端固定于流场中;伞衣部分只考虑加强带,其余部件忽略;流场为亚声速不可压缩流场,且为理想气体的稳态定常流动;不考虑伞绳与流场之间的相互耦合;伞衣为各向同性材料。降落伞部件结构参数见表1。表 1 降落伞部件结构参数 Tab.1 Structural parameters of parachute components 伞衣 伞绳 加强带 吊带 厚度/m 6.0105 截面积/m2 9.6106 1.0105 1.5104 密度/(kg/m3)616 623 660 623 杨氏模量/Pa 5.4108 1.091010 2
22、.181010 1.091010 2.2 网格划分 2.2.1 降落伞网格 本文通过动力学方法建立降落伞初始模型,如图2所示。首先建立十字形伞平面模型,伞衣为四边形单元,伞绳为线单元(见图2(a);然后将伞衣中点固定,对连接绳一端施加远离伞衣的位移载荷,在伞绳作用下,降落伞伞衣将以十字形折叠(见图2(b)。为了保证展开后伞绳长度一致,初始时刻部分伞绳呈折线状态。(a)平面模型 (b)折叠模型(a)Planform model (b)Folding model 图 2 十字形伞网格模型 Fig.2 Simulation model of cruciform parachute 第 3 期 谢淮
23、等:十字形伞开伞充气过程数值仿真研究 37 2.2.2 流场网格 为了与试验风洞尺寸一致,从而方便与试验结果进行对比,将流场区域设置为六面体,尺寸为18 m8 m6 m。流场网格为六面体单元,在降落伞处对网格进行了十字形加密,总单元数为1.7106个,如图3所示。3 开伞过程仿真 3.1 工况设置 在建立数值计算模型的基础上,为了使仿真计算结果更加接近十字形伞实际工作过程中的状态,同时方便与已有的风洞试验结果进行对比,本文设置以下计算工况:来流速度40 m/s,气体密度1.35 kg/m3,温度258 K;流场除入口外,其余壁面均为无反射边界。此工况下获得的仿真结果具有一定的代表性,可以反映十
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 十字形 伞开伞 充气 过程 数值 仿真 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。