中文核心期刊

中国科学引文数据库(CSCD)来源期刊

中国高校优秀科技期刊

中国宇航学会深空探测技术专业委员会会刊

高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

着陆器足垫冲击月壤动态行为离散元仿真分析

林云成 李立犇 赵振家 张荣荣 金聪 邹猛

林云成, 李立犇, 赵振家, 张荣荣, 金聪, 邹猛. 着陆器足垫冲击月壤动态行为离散元仿真分析[J]. 深空探测学报(中英文), 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
引用本文: 林云成, 李立犇, 赵振家, 张荣荣, 金聪, 邹猛. 着陆器足垫冲击月壤动态行为离散元仿真分析[J]. 深空探测学报(中英文), 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
LIN Yuncheng, LI Liben, ZHAO Zhenjia, ZHANG Rongrong, JIN Cong, ZOU Meng. Simulation Analysis of Dynamic Behavior of Lander Footpad Impact on Lunar Regolith[J]. Journal of Deep Space Exploration, 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
Citation: LIN Yuncheng, LI Liben, ZHAO Zhenjia, ZHANG Rongrong, JIN Cong, ZOU Meng. Simulation Analysis of Dynamic Behavior of Lander Footpad Impact on Lunar Regolith[J]. Journal of Deep Space Exploration, 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002

着陆器足垫冲击月壤动态行为离散元仿真分析

doi: 10.15982/j.issn.2095-7777.2020.20190313002
基金项目: 国家自然科学基金资助项目(51775233);中国空间技术研究院预研项目(5010120160759)
详细信息
    作者简介:

    林云成(1987–),男,工程师,主要研究方向:航天器设计与试验研究。通讯地址:北京市海淀区友谊路航天城北京空间飞行器总体设计部(100094)E-mail:gentlycloud@126.com

    李立犇(1992–),男,硕士研究生,主要研究方向:深空探测车辆地面力学。通讯地址:吉林省长春市吉林大学南岭校区工程仿生教育部重点实验室(130022)E-mail:liliben86@163.com

    通讯作者:

    邹猛(1978–),男,教授,博士生导师,主要研究方向:深空探测车辆地面力学。本文通讯作者。通讯地址:吉林省长春市吉林大学南岭校区工程仿生教育部重点实验室(130022)E-mail:zoumeng@jlu.edu.cn

  • ● A dynamic mechanical model of angled footpads impacting on lunar regolith is established. ● Considering the characteristics of low cohesion between lunar soil particles and zero water content, a linear contact stiffness model is used to analyze and establish lunar soil particles. ● Through the simulation software, the changes of the internal force field of the lunar soil after the impact of the footpad can be seen directly, and the changes of the lunar regolith after the impact can be qualitatively analyzed.
  • 中图分类号: V42

Simulation Analysis of Dynamic Behavior of Lander Footpad Impact on Lunar Regolith

  • 摘要: 研究着陆器与月壤的相互作用对着陆稳定性、安全性和着陆机构设计具有重要的意义。建立了具有角度的着陆器足垫冲击月壤的动态作用力学模型,并通过离散单元法分析了足垫在0°、8°、15°和22°冲击角度和1 m/s、2 m/s和4 m/s冲击速度下月壤的动态变化情况,得出了在冲击速度4 m/s、冲击角度为0°时,冲击力最大为5 080.6 N;在竖向位移相同的情况下,冲击角度越大冲击力越小。另外,足垫冲击月壤角度越大对月壤的压实作用越小,其中22°时,前后孔隙率变化为26.5%。冲击角度一定时,速度越大对月壤的压实较滞后,而压实程度较大,冲击速度为4 m/s时,前后孔隙率变化为78.1 %。研究结果对着陆机构设计和着陆控制具有参考意义。
    Highlights
    ● A dynamic mechanical model of angled footpads impacting on lunar regolith is established. ● Considering the characteristics of low cohesion between lunar soil particles and zero water content, a linear contact stiffness model is used to analyze and establish lunar soil particles. ● Through the simulation software, the changes of the internal force field of the lunar soil after the impact of the footpad can be seen directly, and the changes of the lunar regolith after the impact can be qualitatively analyzed.
  • 图  1  足垫冲击月壤力学模型

    Fig.  1  Mechanical model of footpad impacting on lunar soil

    图  2  线性接触刚度模型

    Fig.  2  Linear contact stiffness model

    图  3  离散元模拟月壤双轴试验结果

    Fig.  3  Biaxial test results from DEM simulation

    图  4  完成4 m/s不同冲击角度足垫受力曲线

    Fig.  4  Force curve of 4 m/s footpad cushion at different impact angles

    图  5  冲击角度与冲击力的关系

    Fig.  5  Relation curve between impact Angle and impact force

    图  6  不同刺入深度月壤破坏图

    Fig.  6  Failure maps of lunar soil at different penetration depths

    图  7  不同冲击角度对月壤的破坏情况

    Fig.  7  Destruction of lunar soil at different impact angles

    图  8  不同冲击角度下月壤内部接触力场图

    Fig.  8  Contact field diagram of lunar soil at different impact angles

    图  9  足垫冲击月壤测量圆示意图

    Fig.  9  Schematic diagram of measurement circle of footpad impacting lunar soil

    图  10  4 m/s下不同角度对孔隙率的影响

    Fig.  10  Effect of different angles on porosity at 4 m/s

    图  11  0°下不同冲击速度对月壤孔隙率的影响

    Fig.  11  Effects of different impact velocities on the porosity of lunar soil at 0°

    表  1  月壤物理力学性质

    Table  1  Mechanical characters of lunar soil

    颗粒密度/
    (g.cm–3
    孔隙率/
    %
    内摩
    擦角 /(°)
    凝聚力 /
    kPa
    容重/
    (g.cm–3
    粒径 /
    μm
    2.9~3.540~5025~450.26~1.81.36~1.830~1 000
    下载: 导出CSV

    表  2  月壤离散元接触力学模型参数

    Table  2  Contact mechanics model parameters of discrete elements in soil

    参数参数值
    颗粒数/个3万
    颗粒半径/mm5~12.5
    法向刚度/(N·m–11×106
    切向刚度/(N·m–12.5×106
    颗粒间摩擦系数0.75
    密度/(g·cm–31.63
    初始孔隙率0.21
    重力g/(m·cm–21.63
    法向临界阻尼比0.84
    切向临界阻尼比0.62
    下载: 导出CSV
  • [1] 曾福明, 杨建中. 月球及行星探测器软着陆缓冲器技术综述[C]//中国宇航学会深空探测技术专业委员会会议论文集. 哈尔滨: CDSET, 2005.

    ZENG F M, YANG J Z. Technical overview of soft landing buffers for lunar and panetary probes[C]//Proceedings of the Conference of the Professional Committee of Deep Space Exploration Technology of the Chinese Astronautical Society. Harbin: CDSET, 2005.
    [2] 丁宗茂,王春洁. 探测器触地机关软着陆稳定性分析[J]. 北京航空航天大学学报,2018(3):614-620.

    DING Z M, WANG C J. Stability analysis of explorer in soft landing mode of engine shutdown at touchdown[J]. Journal of Beijing University of Aeronautics and Astronautics,2018(3):614-620.
    [3] WU H Y, WANG C J, DING J Z, et al. Dynamics simulation analysis for novel lander based on two kinds of landing mode[C]//9th International Conference on Measuring Technology and Mechatronics Automation(ICMTMA). Piscataway: ICMTMA, 2017.
    [4] 王少纯,邓宗全. 微型月球着陆器冲击隔离着陆腿研究[J]. 哈尔滨工业大学学报,2004,26(2):180-182.

    WANG S C,DENG Z Q. Research on impact isolation landing leg of miniature lunar lander[J]. Journal of Harbin Institute of Technology,2004,26(2):180-182.
    [5] 丁建中,王春洁,王家俊,等. 着陆工况对月球探测器着陆缓冲性能影响分析[J]. 载人航天,2016,22(1):132-137.

    DING J Z, WANG C J,WANG J J,et al. Analysis of the impact of landing conditions on the lunar probe landing buffer performance[J]. Manned Spaceflight,2016,22(1):132-137.
    [6] 吴晓君,钟世英. 着陆器足垫垂直冲击模型试验研究[J]. 岩土力学,2012(4):1045-1050.

    WU X J,ZHONG S Y. Model test research on vertical impact of lander foot pad[J]. Rock and Soil Mechanics,2012(4):1045-1050.
    [7] 黄晗,李建桥. 着陆器足垫冲击特性模拟试验[J]. 吉林大学学报,2017(4):1194-1200.

    HUANG H,LI J Q. Simulation test of impact characteristics of lander footpads[J]. Journal of Jilin University,2017(4):1194-1200.
    [8] 郝飞,蒋明镜. TJ-1模拟月壤承载特性室内试验研究[J]. 岩土力学,2015(7):1953-1958.

    HAO F,JIANG M J. Experimental mirror. TJ-1 lunar soil bearing characteristics simulation chamber[J]. Rock and Soil Mechanics,2015(7):1953-1958.
    [9] 朱汪,杨建中. 月球探测器软着陆机构着陆稳定性仿真分析[J]. 宇航学报,2009,30(5):1792-1796.

    ZHU W,YANG J Z. Simulation analysis of landing stability of lunar probe soft landing mechanism[J]. Journal of Astronautics,2009,30(5):1792-1796.
    [10] NOHMI M. Modeling for lunar lander by mechanical dynamicssoftware[C]//AIAA Modeling and Simulation Technologies Conference and Exhibit. Reston: AIAA, 2005.
    [11] 龙铝波,卿启湘,文桂林,等. 基于ADAMS的着陆器软着陆稳定性仿真分析[J]. 工程设计学报,2010,17(5):334-338,387.

    LONG L B,QING Q X,WEN G L,et al. Simulation analysis of soft landing stability of lander based on ADAMS[J]. Chinese Journal of Engineering Design,2010,17(5):334-338,387.
    [12] 刘焕焕. 多腿式月球探测软着陆器着陆动力学建模与仿真研究[D]. 南昌: 南昌大学, 2008.

    LIU H H. Modeling and simulation of landing dynamics of multi-legged lunar exploration soft lander[D]. Nanchang: Nanchang University, 2008.
    [13] HUANG H,LI J Q,CHEN B C,et al. Experiment of impact performance of space lander footpad[J]. Journal of Jinlin University,2017,47(4):1194-1200.
    [14] SHUAI Y,HONG N,MING Z,et al. Analysis on landing dynamics and influence factors of near space manned capsule[J]. Journal of Astronautics,2018,39(3):264-274.
    [15] ZHENG G,NIE H,CHEN J,et al. Dynamic analysis of lunar lander during soft landing using explicit finite element method[J]. Acta Astronautica,2018,148:69-81. doi:  10.1016/j.actaastro.2018.04.014
    [16] 梁绍敏,王永滨,季顺迎. 基于离散元方法的月球着陆器冲击月壤过程分析[J]. 航天返回与遥感,2017(4):55-63.

    LIANG S M, WANG Y B, JI S Y. Analysis of landing impact process of lunar landing lander based on discrete element method[J]. Spacecraft Recovery & Remote Sensing,2017(4):55-63.
    [17] 梁绍敏,王永滨,王立武,等. 月球着陆器着陆过程的DEM-FEM耦合分析[J]. 固体力学学报,2019,40(1):39-50.

    LIANG S M,WANG Y B,WANG L W,et al. DEM-FEM coupling analysis of the landing process of the lunar lander[J]. Chinese Journal of Solid Mechanics,2019,40(1):39-50.
  • [1] 杨霞, 李少腾, 贾贺, 刘乃彬, 王垒.  深空大载荷航天器软着陆气囊制造技术 . 深空探测学报(中英文), 2021, 8(3): 324-332. doi: 10.15982/j.issn.2096-9287.2021.20210010
    [2] 王迎春, 王国欣, 赵帆, 张鼐, 王书超.  深层月壤钻取冗余绕组电机控制方法 . 深空探测学报(中英文), 2021, 8(3): 259-268. doi: 10.15982/j.issn.2096-9287.2021.20210014
    [3] 庞勇, 王国欣, 汤滨, 王书超, 赵忠贤.  月球浅层月壤螺旋钻具进芯机制设计 . 深空探测学报(中英文), 2021, 8(3): 252-258. doi: 10.15982/j.issn.2096-9287.2021.20210019
    [4] 陈化智, 刘荣凯, 马超, 姜生元.  基于颗粒单向流动效应的手持式月壤取心装置设计 . 深空探测学报(中英文), 2020, 7(2): 164-170. doi: 10.15982/j.issn.2095-7777.2020.20171124001
    [5] 马超, 刘飞, 曾婷, 赵曾, 庞勇, 樊明旭, 姜生元.  无轴螺旋式模拟月壤主动填充装置研制 . 深空探测学报(中英文), 2019, 6(1): 57-62. doi: 10.15982/j.issn.2095-7777.2019.01.009
    [6] 任德鹏, 李青, 刘振春, 张旭辉.  月面着陆起飞试验技术研究 . 深空探测学报(中英文), 2018, 5(3): 281-285,298. doi: 10.15982/j.issn.2095-7777.2018.6.011
    [7] 孟治国, 李翠, 平劲松, 黄倩, 蔡占川, AlexanderGusev.  月面冯·卡门撞击坑的着陆选址和科学探测目标浅析 . 深空探测学报(中英文), 2018, 5(1): 3-11. doi: 10.15982/j.issn.2095-7777.2018.01.001
    [8] 黄倩, 王太茂, 赵健楠, 孟治国.  “嫦娥4号”月球背面着陆区月壳及深部结构特征 . 深空探测学报(中英文), 2018, 5(1): 41-49. doi: 10.15982/j.issn.2095-7777.2018.01.006
    [9] 申智春, 林小艳, 程坤, 王海鹏.  火星探测器器箭分离冲击响应影响分析与评价 . 深空探测学报(中英文), 2018, 5(5): 483-487. doi: 10.15982/j.issn.2095-7777.2018.05.012
    [10] 鲁亮, 胡宇鹏, 欧阳智江, 向延华, 黄含军, 李思忠.  同位素热源高温-冲击复合环境试验 . 深空探测学报(中英文), 2017, 4(5): 466-470. doi: 10.15982/j.issn.2095-7777.2017.05.007
    [11] 吴伟仁, 王琼, 唐玉华, 于国斌, 刘继忠, 张玮, 宁远明, 卢亮亮.  “嫦娥4号”月球背面软着陆任务设计 . 深空探测学报(中英文), 2017, 4(2): 111-117. doi: 10.15982/j.issn.2095-7777.2017.02.002
    [12] 苏奕粼, 侯绪研, 曹盼, 薛萍萍, 唐天峰, 李龙, 陈涛, 姜生元.  针对在轨维护的爬行机器人足粘附机理研究 . 深空探测学报(中英文), 2017, 4(3): 293-298. doi: 10.15982/j.issn.2095-7777.2017.03.015
    [13] 赖小明, 杜志豪, 王国峰, 王国欣, 莫桂冬.  月壤取芯钻具热特性有限元分析 . 深空探测学报(中英文), 2017, 4(6): 544-551. doi: 10.15982/j.issn.2095-7777.2017.06.007
    [14] 闫晓鹏, 孙海滨, 郭雷.  火星着陆器的大气进入段有限时间抗干扰制导律设计 . 深空探测学报(中英文), 2016, 3(1): 61-67. doi: 10.15982/j.issn.2095-7777.2016.01.010
    [15] 王永滨, 蒋万松, 王磊, 黄伟.  载人登月舱月面着陆缓冲装置设计与研制 . 深空探测学报(中英文), 2016, 3(3): 262-267. doi: 10.15982/j.issn.2095-7777.2016.03.010
    [16] 陈化智, 姜生元, 张明, 张玉良, 沈毅, 全齐全, 梁鲁, 张志贤.  宇航员手持式高频冲击采样装置的设计研究 . 深空探测学报(中英文), 2015, 2(2): 131-136. doi: 10.15982/j.issn.2095-7777.2015.02.005
    [17] 沈毅, 王冬, 姜生元, 刘杰, 张伟伟, 陈化智, 邓宗全.  月壤剖面冲击贯入式探测方案研究 . 深空探测学报(中英文), 2015, 2(3): 213-217. doi: 10.15982/j.issn.2095-7777.2015.03.004
    [18] 黄翔宇, 张洪华, 王大铁, 李骥, 关软峰, 王鹏基.  “嫦娥三号”探测器软着陆自主导航与制导技术 . 深空探测学报(中英文), 2014, 1(1): 52-59.
    [19] 吴伟仁, 王大轶, 毛晓艳, 黄翔宇.  基于月面单幅图像的软着陆障碍识别与安全区选取方法 . 深空探测学报(中英文), 2014, 1(4): 262-268. doi: 10.15982/j.issn.2095-7777.2014.04.003
    [20] 张大鹏, 雷勇军.  深空探测返回舱着陆冲击动力学分析 . 深空探测学报(中英文), 2014, 1(2): 150-155.
  • 加载中
图(11) / 表 (2)
计量
  • 文章访问数:  313
  • HTML全文浏览量:  102
  • PDF下载量:  41
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-13
  • 修回日期:  2019-05-21
  • 网络出版日期:  2020-08-10
  • 刊出日期:  2020-08-10

着陆器足垫冲击月壤动态行为离散元仿真分析

doi: 10.15982/j.issn.2095-7777.2020.20190313002
    基金项目:  国家自然科学基金资助项目(51775233);中国空间技术研究院预研项目(5010120160759)
    作者简介:

    林云成(1987–),男,工程师,主要研究方向:航天器设计与试验研究。通讯地址:北京市海淀区友谊路航天城北京空间飞行器总体设计部(100094)E-mail:gentlycloud@126.com

    李立犇(1992–),男,硕士研究生,主要研究方向:深空探测车辆地面力学。通讯地址:吉林省长春市吉林大学南岭校区工程仿生教育部重点实验室(130022)E-mail:liliben86@163.com

    通讯作者: 邹猛(1978–),男,教授,博士生导师,主要研究方向:深空探测车辆地面力学。本文通讯作者。通讯地址:吉林省长春市吉林大学南岭校区工程仿生教育部重点实验室(130022)E-mail:zoumeng@jlu.edu.cn
  • ● A dynamic mechanical model of angled footpads impacting on lunar regolith is established. ● Considering the characteristics of low cohesion between lunar soil particles and zero water content, a linear contact stiffness model is used to analyze and establish lunar soil particles. ● Through the simulation software, the changes of the internal force field of the lunar soil after the impact of the footpad can be seen directly, and the changes of the lunar regolith after the impact can be qualitatively analyzed.
  • 中图分类号: V42

摘要: 研究着陆器与月壤的相互作用对着陆稳定性、安全性和着陆机构设计具有重要的意义。建立了具有角度的着陆器足垫冲击月壤的动态作用力学模型,并通过离散单元法分析了足垫在0°、8°、15°和22°冲击角度和1 m/s、2 m/s和4 m/s冲击速度下月壤的动态变化情况,得出了在冲击速度4 m/s、冲击角度为0°时,冲击力最大为5 080.6 N;在竖向位移相同的情况下,冲击角度越大冲击力越小。另外,足垫冲击月壤角度越大对月壤的压实作用越小,其中22°时,前后孔隙率变化为26.5%。冲击角度一定时,速度越大对月壤的压实较滞后,而压实程度较大,冲击速度为4 m/s时,前后孔隙率变化为78.1 %。研究结果对着陆机构设计和着陆控制具有参考意义。

注释:
1)  ● A dynamic mechanical model of angled footpads impacting on lunar regolith is established. ● Considering the characteristics of low cohesion between lunar soil particles and zero water content, a linear contact stiffness model is used to analyze and establish lunar soil particles. ● Through the simulation software, the changes of the internal force field of the lunar soil after the impact of the footpad can be seen directly, and the changes of the lunar regolith after the impact can be qualitatively analyzed.

English Abstract

林云成, 李立犇, 赵振家, 张荣荣, 金聪, 邹猛. 着陆器足垫冲击月壤动态行为离散元仿真分析[J]. 深空探测学报(中英文), 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
引用本文: 林云成, 李立犇, 赵振家, 张荣荣, 金聪, 邹猛. 着陆器足垫冲击月壤动态行为离散元仿真分析[J]. 深空探测学报(中英文), 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
LIN Yuncheng, LI Liben, ZHAO Zhenjia, ZHANG Rongrong, JIN Cong, ZOU Meng. Simulation Analysis of Dynamic Behavior of Lander Footpad Impact on Lunar Regolith[J]. Journal of Deep Space Exploration, 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
Citation: LIN Yuncheng, LI Liben, ZHAO Zhenjia, ZHANG Rongrong, JIN Cong, ZOU Meng. Simulation Analysis of Dynamic Behavior of Lander Footpad Impact on Lunar Regolith[J]. Journal of Deep Space Exploration, 2020, 7(2): 171-177. doi: 10.15982/j.issn.2095-7777.2020.20190313002
    • 月球探测及其资源的开发利用离不开着陆器,而足式着陆机构不仅能起到减少冲击能量的作用,更能够起到保护其内部精密仪器的作用,从而保证月球探测器的安全性和稳定性[1]。足式着陆器的软着陆可简化为3个过程:首先,着陆腿竖向冲击月壤过程;其次,着陆腿水平滑移过程;最终,达到力平衡。着陆过程中,足垫是直接与月壤接触的部件。因此,针对着陆器足垫冲击月壤的研究,国内外学者从缓冲器材料、缓冲腿结构、等多方面分析了着陆器着陆过程的冲击力学问题。针对缓冲腿、足垫与月壤相互作用方面,试验表明了冲击能量越大,刺入深度和轴力也越大,月壤密实度对冲击加速度影响峰值较大。

      着陆器在着陆过程中,由于月面凸凹不平且具有一定坡度,加之着陆器在着陆时有一定的俯仰角,因此足垫在接触月壤的实际情况中,往往与月面呈一定的倾角。目前,研究足垫与月壤二者相互之间的力学关系较多是足垫以垂直方向冲击月壤。

      文献[2]提出针对触地机关软着陆稳定性分析,建立了新型探测器的ADAMS(Automatic Dynamic Analysis of Mechanical Systems)动力学分析模型,并引入变推力主发动机及姿控推进器的推力控制模型。文献[3]讨论了新型着陆器在自由落体着陆和低速着陆两种着陆方式下的着陆稳定性,建立了新型着陆器的动力学分析模型,对着陆过程进行了仿真。文献[3-4]讨论了着陆腿和着陆工况对着陆性能的影响。吴晓君等[6]基于对着陆过程的研究,将问题分解为垂直冲击和水平滑动两个方向的运动,并通过力的平衡和动量定理建立了垂直冲击的动力模型。利用独立研发的垂直冲击模型试验装置,探讨了冲击过程中土体密实度、冲击速率和冲击质量等对足垫刺入深度、最大轴力和最大加速度等的影响。试验结果与动力模型进行了对比,结果表明:动力模型能合理地反映出着陆器在垂直冲击阶段的动力特性。文献[6-7]讨论了着陆器足垫冲击月壤的模拟试验。文献[8]给出了着陆动力学分析的数学模型,并在分析铝蜂窝材料缓冲特性的基础上,建立了铝蜂窝缓冲器模型,将其应用于软着陆机构单条着陆腿冲击仿真。文献[9]讲述了机械动力学软件在在月球着陆器建模中的应用。文献考虑月球1/6 g的低重力环境和月壤物理力学性质的影响,以及真实地面环境下无法进行大量的冲击试验来研究足垫冲击月壤后对月壤的影响情况。龙铝波等[10]建立了着陆器足垫–月壤相互作用模型及铝蜂窝的力学模型,在此基础上应用ADAMS软件构建准确的着陆器整体的多体动力学模型,并结合力–角稳定性度量法,为着陆器安全软着陆提出了合理的稳定性判据。该文章着重分析了月面着陆环境和初始条件对软着陆稳定性的影响,并获得了着陆器发生倾覆时不同着陆偏航角与极限水平速度关系的稳定性边界曲线。刘焕焕[11]建立了包括三腿、四腿月球探测器下坡软着陆的二维动力学方程,并对着陆过程中缓冲器的力学特性、足垫与月表的相互作用和着陆过程进行了数学描述。借助所编制的仿真软件,对月球探测器的二维软着陆动力学进行了仿真,对各种着陆条件参数,以及探测器的总体参数,如探测器的回转半径对着陆过载系数、稳定性和着陆腿压溃长度的影响进行了参数敏度分析。黄晗等[12]研究了土壤类型和密实度,脚垫直径和冲击高度对脚垫冲击性能的影响,结果表明,土壤类型和压实度对冲击深度的影响最大,其次是脚垫直径,然后是冲击高度,验证了所提出的模型能够准确预测脚垫峰值冲击力。Shuai等[13]建立了考虑地面弹性变形的联合仿真动力学模型,用于着陆器着陆冲击过程,并通过单腿着陆冲击试验验证了该模型的有效性,结果表明,降低水平着陆速度可以有效地降低着陆器水平着陆过载,提高着陆稳定性。Zheng等[14]提出了模型的能量和质量平衡准则,并结合理论和计算方法,该模型中的缓冲能量吸收特性具有良好的能量吸收能力,其中可以消散高达84%的初始能量。梁绍敏等[15-16]对着陆过程中足垫与月壤的相互作用进行了研究。

      本文建立了具有角度的足垫冲击月壤的动态作用力学模型,通过离散元仿真分析对足垫冲击月壤的动态过程做了分析,研究冲击速度、冲击角度对足垫受力和月壤孔隙率变化情况。

    • 足垫与月壤作用力学模型是基于着陆器足垫以一定速度垂直或倾斜冲击月壤。在足垫冲击月壤时,着陆器的动量和势能被着陆器的吸能结构和月壤的变形吸收,足垫冲击月壤时受力模型如图1所示。其中:${M_1}$为冲击体总质量的1/4;$\theta $为冲击方向与月壤表面法向的夹角;${r_{\rm m}}$为足垫的最大半径;$r$为足垫沿冲击方向接触月壤面积的半径($0 \text{≤} r \text{≤} {r_{\rm m}}$);${K_{\rm s}}$为月壤的刚度系数;$\delta $为足垫沿冲击方向的轴向刺入位移;$\delta_1 $为足垫沿冲击方向的轴向刺入速度;$\delta_2 $为锥形月壤沿着冲击方向的轴向位移;$\dot{\delta}$为锥形月壤沿着轴向运动速度;$\varphi $为月壤的内摩擦角;$\text{π} /2 - \varphi $为锥形体月壤顶角;${M_{\rm s}}$为锥形体有效域月壤质量,其中${M_{\rm s}}{\rm{ = }}\eta \rho {{\rm{r}}^3}\left( {{{ \delta} _1}} \right)$

      图  1  足垫冲击月壤力学模型

      Figure 1.  Mechanical model of footpad impacting on lunar soil

      根据冲击机械部分得出受力平衡方程如下

      $$ {M_1}{{\delta} _1} + \frac{r}{{{r_{\rm m}}}}{K_{\rm s}}\left( {{{\delta} _1} - {{\delta} _2}} \right) - {M_1}g\cos \theta = 0 $$ (1)

      其中:$g$为当地实际重力加速度。

      根据锥形月壤的动量定理得出微分方程如下

      $$ \frac{{\rm d}}{{{{\rm d}_t}}}\left( {\eta \rho {r^3}{{\delta^. }_2}} \right) + {C_{{\rm s}2}}r\left| {{{\mathop {\delta} \limits^. }_2}} \right|\mathop {\delta} \limits^. + {a_1}{{\delta} _2} - \frac{r}{{{r_{\rm m}}}}{K_{\rm s}}\left( {{{{\delta}} _1} - {{\delta} _2}} \right) = 0 $$ (2)

      冲击动态过程为

      $$ \frac{\rm d}{{{{\rm d}_t}}}\left( {\eta \rho {r^3}{{\mathop{\delta}^. }_2}} \right){\rm{ = }}\eta \rho \left( {{r^3}{{\mathop{\delta}^. }_2} + 3{r^2}{{\mathop{\delta}^. }_2} \mathop r\limits^. } \right) $$ (3)

      足垫接触土壤半径与入土深度变化为

      $$ \mathop r\limits^. = \frac{{{{\rm d}_r}}}{{{\rm d}{\delta _1}}}{{\mathop {\delta} \limits^.} _1} $$ (4)

      将式(3)~(4)带入式(2)得到

      $$ \begin{array}{l} \eta \rho \left( {{r^3}{{\mathop {\delta} \limits^. }_2} + 3{r^2}\dfrac{{{{\rm d}_r}}}{{{\rm d}{{\delta} _1}}}{{\mathop {\delta} \limits^. }_1}{{\mathop {\delta} \limits^. }_2}} \right) + {C_{{\rm S}2}}r\left| {{{\mathop {\delta} \limits^. }_2}} \right|{{\mathop {\delta} \limits^. }_2} +\\ {a_1}{{\delta} _2} - \dfrac{r}{{{r_{\rm m}}}}{K_{\rm s}}\left( {{\delta _1} - {{\delta} _2}} \right) = 0 \\ \end{array} $$ (5)

      其中:$\rho $为月壤密度;${C_{{\rm S}2}}$为表征足垫刺入月壤作用的锥形土体与周围月壤发生动量转移现象的系数;${a_1}$为月壤强度系数,与锥形土体有效域质量成正比。这里${C_{{\rm S}2}}$${a_1}$都为经验值。联立方程(1)、(5)可求得足垫在冲击中的运动与受力情况。从而得出,影响足垫与月壤相互作用力学关系的因素有4个:冲击物重量${M_1}$、冲击速度${{V}}$、月壤密度$\rho $以及着陆倾角$\theta $

    • 公式(5)表明,冲击角度和月壤密实度对冲击的能量的吸收具有重要影响,但通过公式计算难以准确获取月壤的动态变形和其性能变化。因此,文章采用离散元软件——PFC(Particle Flow Code)进行仿真分析。

      月壤是指月球表面覆盖着的土壤,一般呈现淡褐、暗灰色。绝大多数月壤元素与地球土壤元素一致,包括氧、硅和铝等。月壤大部分由小于1 mm的具有粘聚性的细颗粒组成,颗粒大小的中间值为40~130 μm,平均大小为70 μm。月壤物理力学特性对月面探测车辆运动影响很大,包括颗粒大小、孔隙率、凝聚力、内摩擦角、容积密度、颗粒比重、颗粒形态等。其基本物理力学性质如表1所示。

      表 1  月壤物理力学性质

      Table 1.  Mechanical characters of lunar soil

      颗粒密度/
      (g.cm–3
      孔隙率/
      %
      内摩
      擦角 /(°)
      凝聚力 /
      kPa
      容重/
      (g.cm–3
      粒径 /
      μm
      2.9~3.540~5025~450.26~1.81.36~1.830~1 000

      文章采用PFC软件分析足垫与月壤作用。根据表1月壤的物理力学性质,选用线性刚度模型建立月壤离散元模型。线性模型如图2所示,图中${{F}}_i^s$为切向接触力;${{F}}_i^n$为法向接触力;${{F}}$为接触合力,并认为变形只发生在接触点处。根据力–位移接触定律可以得出两颗粒间的接触力

      图  2  线性接触刚度模型

      Figure 2.  Linear contact stiffness model

      $$ {{F}}_i^n = {k^n}{U^n}{n_i} $$ (6)

      其中:${k^n}$为法向接触刚度;${{{U}}^n}$为法向接触位移;${{{n}}_i}$为单位法向量。

      研究通过双轴模拟试验,确定颗粒接触力学参数,模拟时设定颗粒在1/6重力加速度条件下达到平衡。通过反复调整颗粒输入参数,直到模拟试验结果与月壤力学性质接近。此时,参数如表2所示。

      表 2  月壤离散元接触力学模型参数

      Table 2.  Contact mechanics model parameters of discrete elements in soil

      参数参数值
      颗粒数/个3万
      颗粒半径/mm5~12.5
      法向刚度/(N·m–11×106
      切向刚度/(N·m–12.5×106
      颗粒间摩擦系数0.75
      密度/(g·cm–31.63
      初始孔隙率0.21
      重力g/(m·cm–21.63
      法向临界阻尼比0.84
      切向临界阻尼比0.62

      此时,主应力差–轴向应变关系曲线如图3(a)所示,体积应变–轴向应变关系曲线如图3(b)所示。

      图  3  离散元模拟月壤双轴试验结果

      Figure 3.  Biaxial test results from DEM simulation

    • 着陆器着陆区域往往选择较平坦的月面,月面的坡度一般不超过10°,考虑着陆区域特殊情况,仿真过程设定4个冲击角度进行研究,角度分别为0°、8°、15°、22°。

      着陆器足垫在冲击月壤时,受到月壤的反作用力,足垫会产生较大的反向加速度并在较短的时间内停止。通过PFC软件模拟时,假定受到的加速度为定值,来分析不同冲击角度足垫的受力分析。如图4为冲击速度达到4 m/s条件下不同冲击角度足垫的受力情况。

      图  4  完成4 m/s不同冲击角度足垫受力曲线

      Figure 4.  Force curve of 4 m/s footpad cushion at different impact angles

      图4分析可知,足垫在相同的沉陷深度下,冲击角度越大,足垫与月壤接触面积越小,故而足垫对颗粒的作用力较小,最终足垫受到的力也就越小。其中0°、8°、15°、22°峰值冲击力分别为5 080.6 N、4 478.8 N、4 224.9 N、3 605.9 N。冲击力峰值在0°时比8°时增加了13.4%、比15°时增加了20.2%、比22°时增加了40.8%。

      图5为4个冲击角度在3个冲击速度下的冲击力峰值图,从图中可以看出,足垫受到的力随着冲击角度增大而减少。速度不同,减小的幅度不同,其中最大减幅出现在4 m/s时,0°比22°的情况下,减幅为40.8%。

      图  5  冲击角度与冲击力的关系

      Figure 5.  Relation curve between impact Angle and impact force

    • 足垫在冲击月壤的过程中,会对月壤造成破坏,文章通过对月壤进行区域标记来研究足垫对月壤的破坏情况。以冲击角度22°为例,图6为不同深度下足垫刺入月壤对月壤的破坏情况。

      图  6  不同刺入深度月壤破坏图

      Figure 6.  Failure maps of lunar soil at different penetration depths

      图6可以看出,随着足垫向下的刺入,对月壤的剪切破坏作用越来越大,主要表现沿着足垫垂直作用方向,对月壤进行不同程度的压实,刺入越深,压实越严重,而沿水平方向,月壤并没有切向剪切作用。

      图7为0°、8°、15°、22° 4个冲击角度在相同刺入位移下的月壤破坏情况。

      图  7  不同冲击角度对月壤的破坏情况

      Figure 7.  Destruction of lunar soil at different impact angles

      图7可知,在沿垂直方向足垫刺入月壤位移相同情况下,随着足垫冲击月壤角度的增大,足垫与月壤的接触面积减小,而足垫刺入月壤位移大的一侧对月壤的剪切破坏程度增大。无论角度大小,对月壤的破坏主要是沿冲击方向做了不同程度的压实,而没有出现较大的沿切向方向的剪切破坏情况。

    • 根据本文模型理论推导,足垫冲击月壤后对月壤作用形成一个圆锥形区域,从足垫作用圆型区域纵剖面来看应是一个楔形区域。通过DEM(Discrete Element Method)软件可以直观地看出足垫冲击后月壤作用的内部力场的变化情况。从而定性地分析月壤受冲击后的变化情况。图8为足垫以不同角度冲击月壤其颗粒内部接触力场图,可以看出形状接近楔形,在垂直方向相同沉陷位移情况下,冲击角度越大,足垫接触面积越小,颗粒内部接触力也越小。

      图  8  不同冲击角度下月壤内部接触力场图

      Figure 8.  Contact field diagram of lunar soil at different impact angles

    • 着陆器在着陆过程中,足垫接触月壤后会对月壤产生破坏并对其有一定的压实作用,月壤的固有孔隙率也随着发生变化,反映了月壤对冲击能量吸收的情况。在PFC2D中对孔隙率的定义按以下公式定义

      $$ {A_p} = \sum {\text{π} {R^2}} $$ (7)
      $$ R = \frac{{{R_{\min }} + {R_{\max }}}}{2} $$ (8)
      $$ \sum {{R^2}} {\rm{ = }}\sum {{{\bar R}^2}} = N{\bar R^2} $$ (9)
      $$ n = 1 - \frac{{{A_p}}}{A} $$ (10)

      其中:A为颗粒所占据区域面积;Ap为所有颗粒面积之和;R为颗粒半径;N为颗粒数。

      通过PFC2D软件对足垫作用下方定义一个半径为60 mm的测量圆进行检测分析,该测量圆圆心距离地表面120 mm。图9为足垫冲击模拟月壤测量圆示意图。

      图  9  足垫冲击月壤测量圆示意图

      Figure 9.  Schematic diagram of measurement circle of footpad impacting lunar soil

      足垫冲击月壤的角度不同时,以垂直方向沉陷位移为标准,足垫对月壤的压实面积也就不同,故而足垫下月壤的孔隙率随足垫冲击角度的变化应有一定的规律。图10为冲击速度4 m/s的情况下,不同冲击角度对月壤孔隙率的影响情况。

      图  10  4 m/s下不同角度对孔隙率的影响

      Figure 10.  Effect of different angles on porosity at 4 m/s

      图10可知,角度越大孔隙率变化情况越缓慢,孔隙率变化率较小。分析得出0°、8°、15°、22°,4个不同冲击角度冲击前后孔隙率变化率分别为78.1%、64.3%、42%、26.5%,这是因为随着冲击角度的增加,单位沉陷位移下足垫接触月壤的面积减少,故而对月壤压实程度较小,孔隙率变化情况较小。

      图11为冲击角度为0°时,不同冲击速度对月壤孔隙率的影响情况。

      图11可知,冲击速度1 m/s时,足垫沉陷约10 mm后,孔隙率开始变化,在速度为2 m/s情况下,沉陷约为13 mm后,孔隙率改变,而在4 m/s的情况下,则是在沉陷21 mm后,孔隙率较明显改变。分析认为:冲击速度越大,对孔隙率的改变程度越滞后。这是由于冲击速度越大,短时间内对表层月壤迅速造成破坏,月壤颗粒相互作用时间极短,而在这段时间内,较深层的月壤并没有受到破坏,故而孔隙率变化较滞后。同时,从图11中可知,冲击速度分别为1 m/s、2 m/s、4 m/s时,月壤孔隙率的变化率分别为50%、61.6%、78.1%,可见速度越大对孔隙率的影响也越大,这是因为冲击速度越大,冲击能量也越大,对月壤的压实破坏程度也越大,故而造成孔隙率变化较大。

      图  11  0°下不同冲击速度对月壤孔隙率的影响

      Figure 11.  Effects of different impact velocities on the porosity of lunar soil at 0°

    • 1)通过受力分析,给出了足垫冲击月壤的力学模型,该模型能够反映冲击角度和月壤密度对冲击受力的影响。

      2)冲击速度越大、足垫受到的冲击力越大,在冲击速度4 m/s、冲击角度为0°时,冲击力最大为5 080.6 N。另外,在竖向位移相同的情况下,冲击角度越大冲击力越小。

      3)足垫冲击月壤角度越大对月壤的压实作用越小,得出0°、8°、15°、22° 4个不同冲击角度冲击前后孔隙率变化率分别为78.1%、64.3%、42%、26.5%,其中22°时,前后孔隙率变化仅为26.5%。

      4)冲击角度一定时,速度越大对月壤的压实较滞后,而压实程度较大,冲击速度分别为1 m/s、2 m/s、4 m/s时,月壤孔隙率的变化率分别为:50%、61.6%、78.1%。

      5)对着陆器在着陆时选择冲击角度和冲击速度有一定的参考意义,选择合适的角度和速度来减小冲击力和孔隙率的变化,以增加着陆时的安全性和稳定性。

参考文献 (17)

目录

    /

    返回文章
    返回