火星进入舱配平翼机构展开冲击动力学分析_杨智杰.pdf

http:/DOI:10.13700/j.bh.1001-5965.2021.0234火星进入舱配平翼机构展开冲击动力学分析杨智杰1，王刚2，赵瑞杰1，王春洁3，赵军鹏1，*(1.北京航空航天大学机械工程及自动化学院，北京100191；2.北京空间飞行器总体设计部，北京100094；3.北京航空航天大学虚拟现实技术与系统国家重点实验室，北京100191)摘要：在着陆巡视器进入火星大气的过程中，配平翼机构会根据指令由收拢状态展开，并在到达指定位置后锁定，进而将进入舱配平攻角降至合理范围内，因此，其展开动力学性能对后续任务的成败至关重要。以配平翼机构功能的顺利实现为背景，研究复合材料构件的冲击动力学分析方法，建立了其有限元分析模型并基于隐式动力学算法对其展开过程进行了仿真。通过与地面试验结果对比，验证了展开动力学分析模型的正确性。在此基础上，对配平翼机构进入火星大气过程中的 2 种气动载荷工况下的展开过程进行了分析；基于 Hashin 理论对碳纤维蒙皮翼板的强度进行了校核，在 2 种气动载荷工况下各铺层的纤维拉伸、纤维压缩、基体拉伸及基体压缩 4 种失效模式所对应的失效因子均处于安全范围。可对类似机构的展开冲击问题研究提供参考。关键词：配平翼机构；有限元分析；展开锁定；冲击；Hashin 理论中图分类号：V415.4文献标志码：A文章编号：1005-5965（2023）02-0422-08配平翼机构是中国火星着陆巡视器的重要组成部分，其主要作用是在进入舱进入火星大气过程中，通过产生反向气动力矩来抵消质心偏移所产生的气动力矩，从而将进入舱攻角调整至合适范围内，进而为降落伞开伞提供先决条件，因此，火星进入舱配平翼机构的展开动力学性能非常重要。配平翼在由收拢状态展开并锁定的过程中可能会产生较大的冲击载荷，从而引起翼板结构上复合材料的损伤，因此，有必要对配平翼机构展开冲击过程进行分析，验证机构设计的合理性。由于地面展开试验的时间成本和经济成本较高，且难以对配平翼机构进入火星大气过程中的实际工况进行准确模拟，因此，通过计算机仿真分析便成了对配平翼机构展开锁定冲击问题进行分析的重要手段之一。现关于航天器机构展开冲击过程仿真的研究主要包括 2 种途径：采用多体动力学理论建立机构的展开动力学分析模型，并基于多体动力学方法进行仿真，其中根据是否考虑机构部件的变形可分为多刚体动力学、多柔体动力学及刚柔耦合动力学分析；采用有限元理论建立机构的展开动力学分析模型，并基于隐式或显式动力学算法进行展开过程仿真。针对太阳翼展开锁定问题，任守志和刘立平1建立了太阳翼多柔体动力学仿真模型，对太阳翼展开过程进行了仿真分析，并给出影响太阳翼展开时间和展开结束角速度的主要因素。王晛等2采用多柔体动力学分析方法，对太阳翼地面展开锁定过程进行了仿真分析，并验证了仿真模型的正确性。濮海玲等3采用类似方法，对太阳翼在轨和地面试验展开过程进行了仿真分析，通过对 2 种仿真结果的对比，评估了阻尼器对太阳翼展开的影响。Zhang 等4针对太阳翼的展开冲击问题，建立了多柔体动力学仿真模型，得到了较为准确的冲击载荷收稿日期：2021-05-07；录用日期：2021-07-09；网络出版时间：2021-07-3009：03网络出版地址： J.北京航空航天大学学报，2023，49（2）：422-429.YANG Z J，WANG G，ZHAO R J，et al.Dynamic analysis of deployment impact of trim-wing mechanism of Mars entry capsulesJ.Journal of Beijing University of Aeronautics and Astronautics，2023，49（2）：422-429（in Chinese）.2023年2月北京航空航天大学学报February2023第49卷第2期JournalofBeijingUniversityofAeronauticsandAstronauticsVol.49No.2结果。荣吉利等5采用有限元方法，以圆形薄膜太阳翼为研究对象，对其展开过程进行了动力学仿真分析，研究了转角驱动函数对太阳翼展开稳定性的影响。针对着陆器着陆冲击问题，吴宏宇等6-7基于多刚体动力学分析方法，建立了着陆器软着陆过程的动力学仿真分析模型，并进行了仿真优化；采用相同方法分别建立了着陆器 2 种模式下的动力学仿真模型，并对其着陆冲击性能进行了分析。逯运通等8提出了一种刚柔耦合动力学分析方法，基于该方法建立了月球着陆器着陆动力学仿真模型并进行了分 析 计 算。吴 宏 宇 等9建 立 了 新 型 腿 式 着 陆器的刚柔耦合动力学分析模型，并对着陆器软着陆过程进行了仿真和优化。梁东平等10-11基于有限元理论，对着陆器地面着陆冲击试验进行了有限元建模和仿真分析；采用非线性有限元理论建立了单腿着陆冲击有限元模型，并基于显式动力学算法进行了着陆冲击动力学仿真。针对天线展开锁定和冲击问题，Li 和 Wang12采用多刚体动力学方法，建立了环形桁架式可展开天线的展开动力学分析模型，并对天线展开过程进行了仿真模拟，但无法得到展开过程中天线变形和应力分布情况。李培13基于绝对坐标框架，建立了大型星载环形桁架天线有限元模型,并对其展开过程进行了动力学分析。李团结等14结合刚柔耦合动力学和有限元分析，实现了可展开天线展开过程的仿真分析。对于太阳翼展开锁定、着陆器着陆冲击和天线展开锁定冲击等问题，研究者已做了大量仿真和研究工作，但针对配平翼机构展开冲击问题的研究较少，且复合材料构件冲击损伤的仿真及校核难度较大，尚没有形成系统性方法。而相比于多体动力学分析方法，有限元仿真分析可以在获得机构展开性能的同时十分方便地对其结构强度进行校核，并获得各部件在给定损伤容限下的潜在损伤模式，故本文基于有限元仿真分析，研究含复合材料构件的火星着陆巡视器配平翼机构展开过程分析方法。首先，根据火星着陆巡视器配平翼机构的组成和结构特点研究其有限元建模方法，然后基于隐式动力学算法研究其展开过程仿真分析方法，通过与地面试验结果对比验证了所建立展开动力学分析模型的正确性。在此基础上，对配平翼机构在着陆巡视器进入火星大气过程中 2 种气动载荷工况下的展开过程进行了分析，并基于 Hashin 理论对碳纤维蒙皮翼板的强度进行了校核，有效验证了配平翼机构设计的合理性。本文所提基于复合材料的结构有限元建模、仿真和失效分析方法可为类似航天器机构展开冲击动力学分析问题的研究提供参考。1配平翼机构的原理与有限元建模如图 1 所示，火星着陆巡视器配平翼机构主要由复合材料翼板、展开臂组件、翼板连接架、连杆、曲柄和阻尼器组成，其中翼板与展开臂组件固连，展开臂组件、连杆和曲柄共同构成曲柄摇杆机构，展开臂组件与翼板固连并通过绕根部转轴转动将翼板展开至指定位置。阻尼器位于曲柄末端，用于限制翼板展开过程中的最大速度。在火星进入舱进入火星大气前，配平翼机构处于收拢状态；配平翼机构在接收到展开指令后，在驱动力矩的作用下展开至指定位置。左右展开臂组件连杆曲柄阻尼器翼板前端测点展开臂组件旋转中心翼板连接架复合材料翼板传感器测点xyz图1配平翼机构有限元模型Fig.1Finiteelementmodeloftrim-wingmechanism该配平翼机构需要满足以下动力学指标：机构展开时间小于等于 800ms；复合材料翼板各铺层不发生失效。为了进行展开冲击动力学分析，首先根据火星进入舱配平翼机构的组成和结构特点建立其有限元模型，将复合材料翼板、展开臂组件和翼板连接架等效为壳单元，将曲柄和连杆等效为梁单元。复合材料翼板结构采用常规壳单元(conventionalshell，CS)进行模拟，并对各铺层的材料、区域、层数、厚度和角度进行定义。展开臂组件、翼板连接架与复合材料翼板之间的螺栓连接均采用多点约束(multiplepointconstraint,MPC)模拟，连杆与翼板连接架、连杆与曲柄之间的转动副均采用连接器单元(connectorelement，CE)模拟。在分析中，对曲柄末端和展开臂组件旋转中心处除绕 y 轴旋转以外的自由度进行约束，驱动力矩加载在曲柄末端，作为翼板展开的动力。在配平翼机构中，展开臂组件、翼板连接架、连杆和曲柄等结构的材料为铝合金，翼板使用铝蜂窝夹芯复合材料，该材料由上下对称铺设的 4 层碳纤维复合材料蒙皮和中间的铝蜂窝芯层组成，上层蒙皮角度分别为 0、45、45和 90。阻尼器根据试验中不同速度下阻尼力的测试第2期杨智杰，等：火星进入舱配平翼机构展开冲击动力学分析423结果进行建模，阻尼力矩与扭转角速度之间的关系曲线，如图 2 所示。123456050100150200250300350400阻尼力矩/(Nm)扭转角速度/(rads1)图2阻尼力矩与扭转角速度关系曲线Fig.2Relationshipcurvebetweendampingtorqueandtorsionalangularvelocity2展开过程冲击动力学分析2.1隐式动力学分析方法基于所建立的有限元模型，利用隐式动力学分析方法，对配平翼机构进行展开动力学分析，以获取配平翼机构在各个工况下由收拢状态到展开锁定过程的展开时间、展开性能、翼板应力状况、强度裕度与在给定损伤容限下的潜在损伤模式。根据动力学分析的基本理论，对于任意结构均可列出其动力学方程15：MU+CU+KU=R（1）UUU式中：M 为质量矩阵；C 为阻尼矩阵；K 为刚度矩阵；R 为外部作用载荷矢量，主要包括与时间相关的外力矢量及与位移和速度相关的非线性外力矢量和边界约束反力矢量；、分别为对应有限元组合体的位移、速度和加速度矢量。t+t对于动力学分析，利用隐式时间积分，分析时刻的平衡方程：Mt+tU+Ct+tU+Kt+tU=t+tR（2）根据 Newmark 积分法，系统的位移和速度向量可利用式(3)和式(4)进行求解：t+tU=tU+tUt+(12)tU+t+tUt2（3）t+tU=tU+(1)tU+t+tUt（4）式中：依据 Newmark 积分法的特点，和 分别取1/2 和 1/4。将求得的位移和速度向量代入式(2)，即可求解结构在不同时刻的加速度。为获得翼板在给定损伤容限下潜在的损伤模式，使用 Hashin 理论作为纤维复合材料的失效准则，该准则考虑纤维拉伸、纤维压缩、基体拉伸及基体压缩 4 种破坏模式，4 种破坏模式下的损伤因子计算表达式如下16-17：11 0纤维拉伸失效（）：Ftf=(11XT)2+(12SL)2（5）11 0纤维压缩失效（）：Fcf=(11XC)2（6）22 0基体拉伸失效（）：Ftm=(22YT)2+(12SL)2（7）22 0基体压缩失效（）：Fcm=(222ST)2+(YC2ST)21 22YC+(12SL)2（8）XTXCYTYCSLST 11 22 12式中：为纵向拉伸强度；为纵向压缩强度；为横向拉伸强度；为横向压缩强度；为纵向剪切强度；为横向剪切强度；为纤维压缩剪切应力分布系数；、为有效应力张量。对翼板每个铺层每个时刻的失效因子进行计算，若存在某种破坏模式对应失效因子的值大于 1，则该铺层会发生破坏；否则，在 4 种失效模式下翼板铺层均不会发生失效。2.2仿真模型试验验证为验证有限元模型的准确性，以地面试验中的展开时间和翼板传感器测点处 x 方向最大加速度数值作为目标，对所建立的有限元分析模型进行修正。配平翼机构在实际展开时会受到飞行气动载荷的作用，阻碍配平翼的展开。在地面试验中，在曲柄末端施加大小为 240Nm 的驱动力矩，同时在展开臂组件旋转中心处加载如表 1 所示的等效阻力矩来模拟试验气动载荷的作用。表1试验气动载荷数据Table1Dataoftheaerodynamicloadtest转角/()试验气动载荷/(Nm)00243.79548.538413.2618114同时，基于隐式动力学分析方法对配平翼机构展开过程进行分析，驱动力矩及气动载荷模拟阻力矩的设