基于多工况的新型着陆器软着陆性能优化-吴宏宇.pdf

《基于多工况的新型着陆器软着陆性能优化-吴宏宇.pdf》由会员分享，可在线阅读，更多相关《基于多工况的新型着陆器软着陆性能优化-吴宏宇.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2017年 4月第43卷第4期北京航空航天大学学报Journal of Beijing University of Aeronautics and AstronauticsApril 2017V0143 No4http：bhxbbuaaeduCB jbuaabuaaedu，cnDOI：1013700jbh1001596520160296基于多工况的新型着陆器软着陆性能优化吴宏宇1，王春洁12一，丁建中1，满剑锋13，罗敏3(1北京航空航天大学机械工程及自动化学院，北京100083；2北京航空航天大学虚拟现实技术与系统国家重点实验室，北京100083；3北京空间飞行器总体设计部，北京100094

2、)摘 要：以新型腿式着陆器为研究对象，建立其刚柔耦合动力学分析模型，实现着陆器软着陆过程的仿真。通过仿真计算，确定着陆器最易翻倒、底面最易与星球表面岩石碰撞、主体承受最大碰撞力的3组恶劣着陆工况。分析着陆器缓冲机构构型选取设计变量，基于仿真得到的3组恶劣工况，应用第二代非劣排序遗传算法(NSGAlI)实现着陆器软着陆性能的优化，优化目标为增强着陆器抗翻倒能力、降低着陆器底面与星球表面岩石碰撞的可能性、降低着陆器主体最大受力值。将优化所得参数代入模型重新进行仿真，着陆器不再发生翻倒，着陆平台底面与星球表面最小距离提高42，主体最大受力值降低121。关 键 词：腿式着陆器；刚柔耦合模型；动力学分析

3、；软着陆性能；多目标优化中图分类号：V4236文献标识码：A 文章编号：10015965(2017)04077606我国已完成着陆器软着陆月球表面并利用自动探测机器人对月球进行探测。作为航天大国，我国未来将陆续开展其他星球的探测工程嵋。腿式着陆器作为常用的着陆缓冲装置，对月球之外其他星球进行探索还可能用到该装置。31，因此必须对腿式着陆器持续进行研究。为了检验腿式着陆器的软着陆性能并减少物理样机试验的成本，开展着陆模拟试验的动力学仿真显得尤为重要。由于着陆器的软着陆性能直接决定探测器能否安全抵达星球表面”o，对着陆器的软着陆性能进行优化成为设计阶段必不可少的环节。文献5-9提出了着陆器的动力学

4、分析方法，给出对软着陆性能影响较大的因素。文献10提出月球着陆缓冲器刚柔耦合模型的建模与分析方法。文献11运用响应面法，对单一工况下月球着陆器的软着陆性能进行了多目标优化。以上研究内容局限于着陆器单个性能指标或单一工况下多个性能指标的改进。相比于月球，其他星球的环境更加难以预测，着陆器各性能指标可能不在同一工况下达到最差，有必要在多组工况下，对着陆器多个性能指标进行优化。本文建立了新型腿式着陆器的刚柔耦合动力学分析模型。探索着陆器最易翻倒、底面最易与星球表面岩石碰撞、主体承受最大碰撞力时对应的着陆工况。基于以上3组工况，以增强着陆器抗翻倒能力、降低着陆器主体底面与星球表面岩石碰撞的可能性、降低

5、着陆器主体受力为目标进行优化设计。应用第二代非劣排序遗传算法(NSGA-II)寻优，得到的优化结果明显提高了着陆器的软着陆性能。1 着陆器刚柔耦合动力学模型11新型腿式着陆器构型本文研究的新型着陆器包括主体结构和4套收稿日期：2016-04-14；录用日期：2016-0513；网络出版时间：2016-05，18 17：28网络出版地址：WWWcnkinetkcmsdetail112625V201605181728001html通讯作者：E-mail：wangcjbuaaeducn引用格式：吴宏字，王春洁，丁建中，等基于多工况的新型着陆器欺着陆性能优化EJ北京航空航天大学学报，2017，43似)

6、：776781，WU H Y，WANG C J，DING J Z，et a1Soft landing performance optimization for力ovej lander based O1 muhiple working conditions I】Journal of Beijing University of Aeronautics and Astronautics，2017，43(4)：776781(in Chinese)万方数据第4期 吴宏宇，等：基于多工况的新型着陆器软着陆性能优化 777对称分布的着陆缓冲机构。主体在动力学分析中作为负载，着陆缓冲机构用于吸收着陆器着陆过程

7、产生的冲击力，每套缓冲机构由主支柱、辅助支柱、足垫和外部缓冲件组成，如图1所示。主支柱起支撑作用；辅助支柱由内、外简和缓冲元件组成，内、外筒之间存在缓冲元件，用于吸收压缩载荷；足垫用来增大着陆器与星球土壤的接触面积，防止着陆时着陆腿过度下陷；外部缓冲件与主体固连，当主支柱受到较大冲击力时，外部缓冲件间接受力，发生弯曲变形，进而产生缓冲弯矩，吸收冲击能量，防止主体受力过大。辅助支柱外筒辅助支柱内筒足垫图1 着陆缓冲机构Fig1 Landing buffer mechanism12着陆缓冲机构力学模型辅助支柱内部的缓冲元件是着陆器着陆过程吸收冲击能量的主要结构，元件只能承受压缩载荷。在辅助支柱内外

8、筒之间施加与内、外筒相对位移(D)有关的力(F。)代替缓冲元件的作用效果，辅助腿缓冲机构示意图及D与F。关系曲线如图2所示。辅助支柱缓冲能力较强，外部缓冲件发生侧向弯曲的可能性很小。由于使用较多柔体构件导致动力学模型运算速度缓慢，使用如图3方法对萋誊 薯一-一-毫鞫外筒 缓冲力 内筒(a)辅助支柱缓冲机构Dm(b)D-关系曲线图2辅助支柱缓冲机构及力学特性Fig2 Buffer mechanism and mechanicalbehavior of assistant pillar外部缓冲件进行刚体等效。建立两段刚体模型，质量特性参照外部缓冲件柔体模型设置，冈0体1绕转轴O转动，刚体2与刚体1

9、只能沿轴线相对移动。刚体1转轴处施加与转角(0)相关的转矩(M)，模拟外部缓冲件受弯时产生的缓冲弯矩；刚体2上施加与0相关，沿刚体1轴线方向的运动(S)，模拟外部缓冲件自由端位置变化。在有限元软件中建立外部缓冲件模型，施加转矩，测得0与M、0与s关系曲线如图4所示，输出数据文件，导入刚体模型。考虑等效模型质量、运动与力学特性均由有限元软件计算得到，即等效模型力学性能与柔体模型一致。足垫触地时与星球表面的相互作用力分解为法向接触力F。和切向摩擦力F，接触力F。采用非线性阻尼弹簧碰撞方程描述：F。=硒。+C6 (1)式中：K为接触刚度；6为接触变形；e为非线性指图3 外部缓冲件等效方法Fig3 E

10、quivalent method of external buffer0(o)(a)良M关系曲线o(。)(b)0-S关系曲线图4外部缓冲件力学特性Fig4 Mechanical behavior of external buffer万方数据778 北京航空航天大学学报 2017年数；C为阻尼系数。以上系数均与星球土壤特性 表1着陆工况各参数取值密切相关。Table 1 Value of parameters of landing condition摩擦力F，采用库伦摩擦方程描述：Ff=pF。 (2)式中：肚为足垫与星球表面间的摩擦系数。13刚柔耦合动力学模型建模刚柔耦合多体系统中，构件既存在大

11、范围的刚体运动，也存在局部的弹性形变引。着陆器动力学模型中，足垫作为柔性构件可发生弹性形变，其余构件均为刚体构件。采用有限元软件Patran建立足垫有限元模型，导人多体动力学分析软件ADAMS中作为柔体构件；采用ADAMS软件建立其余构件的刚体模型；应用ADAMS软件对着陆器着陆过程的动力学问题进行求解。2 着陆器软着陆极限工况的确定21软着陆性能评判原则根据着陆器安全着陆的性能指标，结合着陆器的构型，提出以下3种软着陆性能评判原则：1)着陆器在着陆过程中不发生翻倒，即着陆器质心距包含任意两足垫中心点的竖直翻倒墙平面的最小距离大于设计许用值“。2)着陆过程中辅助支柱外筒、外部缓冲件与主体连接点

12、处的最大碰撞力不超过允许值，保证着陆器主体不会因为承受过大的碰撞力而受损。3)当着陆地点摩擦系数较小时，着陆器在着陆过程中，其着陆平台底面中心点距离星球表面的最小距离应大于预留间距值，避免着陆器底面与星球表面岩石发生碰撞而导致着陆器主体受损。22着陆工况参数着陆工况的参数包括：式(2)中着陆器足垫与星球表面间的摩擦系数(t)、星球表面坡度(d)与凹坑分布情况(n)、着陆器触地瞬间的竖直速度(移。)、水平速度(叱)、着陆偏航角(0。)、着陆俯仰角(0，)以及着陆器整机质量(m)，如图5所示。着陆工况各参数的取值如表1所示。图5着陆工况的参数Fig5 Parameters of landing c

13、ondition参数 取值弘a(o)叱(ms。)vz(ms。)Op(。)Of(o)mkg03，080，1，2，3351+10453+3700800表1中，n取0，1，2，3分别表示：着陆器没有着陆腿落入凹坑，2号与4号着陆腿落入坑中，1号与4号着陆腿落入坑中，1号、2号与4号着陆腿落入坑中；口。取正值表示水平速度方向背离坡面，负值表示水平速度朝向坡面；0，取正值表示着陆器俯仰角为仰角，负值表示俯仰角为俯角；遍历过程偏航角以1。为增量从0。增加至45。将表1中各参数取值进行组合，得到的着陆工况个数为214124622，即2 944种着陆工况。23软着陆极限工况参数取值着陆器在特定工况下着陆，着陆

14、过程中着陆器质心至翻倒墙距离取得最小值(k)；或着陆过程中主体受碰撞力取得最大值(FM)；或着陆过程中着陆平台底面中心点至星球表面距离取得最小值(。)，则称该工况为一种极限工况，对应的最值称为一种极限指标值。基于2944种着陆工况，通过仿真计算，得到3组极限工况及对应极限指标如表2所示。根据表2结果可知：着陆器质心距翻倒墙最小距离低于许用值(30 mm)，该情况认为着陆器已经翻倒；星球表面摩擦系数较小时，着陆平台底面中心点距星球表面最小距离低于预留间距值(760 mm)，着陆器极易与星球表面岩石碰撞；由于辅助支柱缓冲能力较强，主体所受的最大碰撞力位于主体与外部缓冲件连接处，接近允许值(30 k

15、N)。为了提升着陆器的软着陆性能，下文利用多目标优化方法对缓冲机构进行优化设计。表2极限工况计算结果Table 2 Calculation results of the worst landing condition万方数据第4期 吴宏宇，等：基于多工况的新型着陆器软着陆性能优化 7793着陆器软着陆性能多目标优化31设计变量的选取根据缓冲机构构型可知，增强外部缓冲件的缓冲力矩，可以提升着陆平台底面至星球表面距离的最小值，但主体受力可能超过许用值；增强辅助腿的缓冲力，可以提升着陆平台底面至星球表面距离的最小值，降低外部缓冲件与主体问的作用力，但缓冲性能降低导致着陆器更易翻倒；如图6所示，增加足

16、垫至着陆器中心轴线的距离(d：)，可使质心至翻倒墙初始距离增加，提升抗翻倒能力，但足垫受力点与外部缓冲件转轴水平距离(d)随之增加。根据力矩计算方程：M，=Fdcos dF+Fhsin d， (3)若着陆器足垫承受同样的冲击力(F)，足垫受力点与外部缓冲件转轴竖直距离h不变，冲击力与竖直方向夹角O，不变，d值增大，冲击转矩M，增大，着陆腿向外扩张能力增强，进而导致着陆平台底面中心点至星球表面最小距离减小，着陆器底面更易与星球表面岩石相撞。在着陆器主体结构与缓冲机构连接位置确定的情况下，同时避免着陆器竖直方向整体尺寸的增加，选取外部缓冲件缓冲转矩放大倍数(|j，)、辅助支柱缓冲力放大倍数(k，)

17、、足垫至着陆器中心轴线的距离(d：)作为设计变量进行多目标优化。根据图6可知，d：值增加导致着陆腿变长，着陆器整体尺寸随之增加。本文提到的缓冲元件提供的缓冲力与其尺寸成正比。为了保证着陆器整体尺寸不超出运载火箭可提供的最大安置空间，同时保证缓冲元件与着陆器其他结构不存在干涉，确定设计变量取值范围如表3所示。图6着陆器性能分析Fig6 Performance analysis of lander表3设计变量取值范围Table 3 Value range of design variables设计变量 取值范围k，k：d：m32软着陆性能多目标优化的实现选取。三。、FM作为优化目标，约束”L。、F

18、。满足23节给定的许用范围，结合设计变量的取值，建立数学模型如下：工(L)X工(u)minLT，一三H，FMSt30一LT0760一LH0FM一350 (4)式中：X=(k。，k：，d；)7代表设计变量；工。与xo分别为设计变量的下界与上界。本文采用NSGAII进行优化计算，算法基于帕累托最优理论，具有迭代次数少、寻优能力强等特点。“。算法参数设置如表4所示1。综上所述，着陆器软着陆性能多目标优化具体流程如图7所示。经过优化计算，得到帕累托最优解集如表5所示。根据图6可知，d：值增加导致着陆腿变长，着陆器整体尺寸随之增加，为了减小着陆器整体尺寸，降低着陆器在运载火箭中的安置空间，选取帕累托解集

19、中d，值最小的结果作为最优解，重新进行2 944次仿真计算，获得极限指标值与优化前结果对比如表6所示。表4优化参数取值1Table 4 Value of optimal parameters刚柔耦合动力学模型的建立3组极限工况的选取优化数学模型的建立NSGA一1I寻优最易翻倒工况l l最易碰撞岩石计算三， l工况计算。主体易被刺穿工况计算民迭代是否终止?Y帕累托最优解的选取验证优化结果N图7多目标优化流程图Fig7 Flowchart of multiobjcctive optimization万方数据780 北京航空航天大学学报 2017年表5帕累托最优解集 nese)Table 5 Par

20、eto optimal solution aggregate 2李萌腿式着陆缓冲装置吸能特性及软着陆过程动力学仿表6优化前后性能对比Table 6 Comparison of performance before andafter optimization根据表6可知，相比于初始设计，优化后的参数使三。增加42，F。削弱121，L，显著增加，着陆器在着陆过程中不再翻倒，软着陆性能获得明显提升。4 结 论1)建立了新型腿式着陆器刚柔耦合动力学分析模型，实现着陆器软着陆过程的仿真，提出新型着陆器软着陆性能评判原则。2)根据着陆工况各参数的取值，确定2944种着陆工况，针对所有工况进行仿真计算，得到

21、着陆器着陆过程中最易受损的3组恶劣工况。3)应用NSGA一算法，基于仿真得到的3组恶劣工况，对着陆器的软着陆性能进行多目标优化。优化后，着陆器在2 944种工况下着陆不存在翻倒现象，着陆平台底面中心点与星球表面的最小距离大于预留间距值，且主体受力得到明显地减弱。参考文献(References)1杨建中，满剑锋，曾福明，等“嫦娥三号”着陆缓冲机构的研究成果及其应用J航天返回与遥感，2014，35(6)：20_27YANG J Z，MAN J F，ZENG F M，et a1Achievements and ap-plieations of landing gear for Change-3 la

22、nderJSpacecraftRecoveryRemote Sensing，2014，35(6)：20-27(in Chi一真研究D哈尔滨：哈尔滨工业大学，2013：110LI MResearch on energy absorbers of leggedtype lander anddynamic simulation on its soft landing processDHarbin：Harbin Institute of Technology，2013：1-10(in Chinese)3蒋万松，黄伟，沈祖炜，等月球探测器软着陆动力学仿真J宇航学报，2011，32(3)：462-469J

23、IANG W S，HUANG W，SHEN Z W，et a1Soft landing dynamic simulation for lunar explorerJJournal of Astronautics，2011，32(3)：462-469(in Chinese)4ADAM D S，MIGUEL S M，TOMMASO P R，et a1Mars sciencelaboratory entry，descent，and landing system development chal-lengesJJournal of Spacecraft and Rockets，2014，51(4)：9

24、9410035DING J Z，WANG C JFast modeling for lunar landing dynam-ics analysis：AIAA-2016-1191RReston：AIAA，20166刘晓宇着陆器软着陆缓冲稳定性仿真分析D长沙：湖南大学，2013：1928LIU X YSimulation analysis on soft landingS buffer and stability of lunar explorationDChangsha：Hunan University，2013：1928(in Chinese)7陈金宝，聂宏，万峻麟深空探测着陆器数字化设计及

25、着陆性能影响因素J航空学报，2014，35(2)：541554CHEN J BNIE H，WAN J LDigital design and landing per-formance influence factors of deep space landerJActa Aero-nautica et Astronautica Sinica，2014，35(2)：541554(in Chi-nese)8罗昌杰，邓宗全，刘荣强，等基于零力矩点理论的腿式着陆器着陆稳定性研究J机械工程学报，2010，46(9)：38-45LUO C J，DENG Z Q，LIU R Q，et a1Landing s

26、tability investigation of legged-type spacecraft lander based on zero momentpoint theoryJJournal of Mechanical Engineering，2010，46(9)：38-45(in Chinese)9罗松柏，赵永嘉月球软着陆动力学分析与仿真J北京航空航天大学学报，2012，38(2)：185190LUO S B，ZHAO Y JDynamic analysis and simulation of softlanding for lunar landerJJournal of Beijing

27、University ofAeronautics and Astronautics，2012，38(2)：185-190(in Chinese)10逯运通，宋顺广，王春洁基于刚柔耦合模型的月球着陆器动力学分析J北京航空航天大学学报，2010，36(11)：13481352LU Y T，SONG S G，WANG C JDynamic analysis for lunarlander based on rigid-flexible coupled modelJJournal ofBeijing University of Aeronautics and Astronautics，2010，36(

28、11)：13481352(in Chinese)11王家俊，王春洁，宋顺广基于响应面法的月球着陆器软着陆性能优化J北京航空航天大学学报，2014，40(5)：707711WANG J J，WANG C J，SONG S GPerformance optimization oflunar lander based on response surface methodologyJJournalof Beijing University of Aeronautics and Astronautics，2014，40(5)：707-711(in Chinese)12BRIAN C S，MELISSA

29、S，PETER ITelecommunications per-万方数据第4期 吴宏宇，等：基于多工况的新型着陆器软着陆性能优化 781formance during entry，descent，and landing of the Mars science HUANG J JZHENG L X，LIU G Q，et a1Optimization algolaboratoryJJournal of Spacecraft and Rockets，2014，5 1 rithm and experiment of two-disk rotor systemJJoumal of(4)：12371250

30、 Aerospace Power，2016，31(1)：65-71(in Chinese)1 13ZUPP G A，DOIRON H HA mathematical proeedure for Dre1415dicting the touchdown dynamics of a soft1anding vehicl。：NASATND-7045RWashington，DC：NASA1971MARLER R T，ARORA J SSurvey of multiobjective optimization methods for engineeringJStructural and Muhidisc

31、iplinary Optimization，2004，26(6)：369-395黄晶晶，郑龙席，刘钢旗，等双盘转予系统优化算法与试验J航空动力学报，2016，31(1)：6571作者简介：吴宏宇男，硕士研究生。主要研究方向：航天机构的设计与动力学分析。王春洁女，博士，教授，博士生导师。主要研究方向：机械设计及理论。Soft landing performance optimization for novel lander based onmultiple working conditionsWU Hongyul，WANG Chunjiel2一，DING Jianzhon91，MAN Jianf

32、en91，LUO Min3(1School of Mechanical Engineering and Automation，Beijing University ofAeronautics and Astronautics，Beijing 100083，China；2State Key Laboratory of Virtual Reality and Systems，Beijing University of Aeronautics and Astronautics，Beijing 100083，China；3Beijing Institute of Spacecraft System E

33、ngineering，Beijing 100094，China)Abstract：A rigidflexible coupled dynamics model was established for simulating and analyzing the softlanding process of the novel leg type landerThree groups of bad landing conditions，in which the lander mosteasily overturns，bottom surface of the lander most easily co

34、llides with rocks on the surface of the planet，andbody of the lander bears the greatest impact force，were found by iterating over landing parametersAccordingto the configuration of buffering mechanism，design variables of optimization were selectedBased on the threegroups of bad landing conditions an

35、d nondominated sorting genetic algorithm(NSGAH)，a multiobjective optimization method was applied to enhance the ability of the lander to resist overturning，reduce the possibility that the bottom surface of lander collides with rocks on the surface of the planet，and reduce the impacton the body of la

36、nderIn the simulation using optimized parameters，the model does not overturn any longerThe minimum distance between bottom surface of the landing platform of the lander and surface of the planetincreases by 42。and the impact on the body of lander reduces by 121Keywords：leg type lander；rigidflexible coupled model；dynamical analysis；soft landing performance；multiobjective optimizationReceived：2016-0414；Accepted：2016_0513；Published online：2016_05一18 17：28URL!WWWcnkinetkemsdetail112625V201605181728001htmlCorresponding authorEmail：wangcjbnaaeduca万方数据