2. 吉林大学工程仿生教育部重点实验室, 吉林 长春 130022
2. Key Lab for Bionics Engineering of Education Ministry, Jilin University, Changchun 130022, China
21世纪是人类全面开展深空探测和开发的新时代,火星与地球的构造、环境具有相似性,更是成为世界各国探测的重点[1-3]。2020年是火星探测的超级年,阿联酋发射“希望号”火星探测器并成功实现环火探测;美国国家航空航天局(NASA)发射“毅力号”火星探测器实现火星表面巡视探测[4-5];中国发射的“天问一号”火星探测器也已成功实现火星环绕和着陆,目前“祝融号”火星车已成功驶离着陆平台至火星表面,并开展了火星表面巡视探测。
迄今为止,成功实现火星表面巡视探测的火星车有美国的“索杰纳”“勇气号”“机遇号”“好奇号”“毅力号”以及刚刚登陆火星的中国“祝融号”。由于火星表面部分区域土壤松软,且部分区域环境较为恶劣,因此火星车在执行火星探测任务的过程中,其不确定性较高[6]。其中,“勇气号”多次在松软地面上出现沉陷过大的情况,且脱困难度大,只能采用车轮前后移动以及转向辅助的方式进行脱困,耗时过长。“机遇号”在2006年5月探测过程中,后轮陷入了松散的火星表面土壤中,其移动能力下降了95% 左右,导致其不能正常行驶[7]。因此,火星车沉陷机理以及脱困策略已经成为现如今火星车移动系统研制中的关键技术之一[8]。
部分学者通过仿真分析的方法来模拟火星表面可能出现的异常情况。其中,Iagnemma[9]通过多刚体动力学软件模拟了火星车在火星表面松软环境下的运动情况,开发了一种将遥测技术与轨道卫星观测数据相结合的技术,该技术可用于预测火星车在火星表面开阔区域的行驶路径。此外,该方法可协助工程师制定火星车相关在轨运行策略,同时该结果也被应用到“勇气号”火星车的在轨运行当中。Johnson等[10]通过离散元方法研究了高滑转情况下火星车沉陷以及挂钩牵引力的情况,可对火星车路径规划起到指导作用。NASA通过火星车试验对滑转率进行了标定,并对滑转率阈值进行了定义,但在火星车实际工作当中,滑转率情况难以进行实时测量[11]。
本文以“祝融号”火星车为研究对象,通过理论分析与试验验证相结合的方式对火星车沉陷机理、沉陷判别准则以及沉陷后脱困策略展开研究,可为“祝融号”火星车在轨使用提供参考。
2 沉陷机理分析与判别准则(Mechanism analysis and determination criterion of sinkage) 2.1 火星车移动系统简介“祝融号”火星车采用六轮驱动主副摇臂悬架形式,且为主动悬架形式,相较被动悬架增加夹角调整机构和离合器机构,如图 1所示。移动系统机构主要包括6套车轮驱动机构、6套车轮转向机构、2套夹角调整机构、2套离合器机构和1套差动机构,以及由2套主摇臂、2套副摇臂、6套转向臂组成的悬架结构。移动系统的17套机构中除差动机构为被动机构外,其余16套机构均为主动机构。
“祝融号”火星车由于使用主动悬架的形式,因此与传统火星车相比具有一些优势。首先,其车厢高度可调,可在运输过程中收紧车体,节省一定空间;在着陆后可将车体升回指定位置,继而进行正常探测。其次,该种设计增加了一种特殊的行走方式(蠕动行走),为“祝融号”火星车移动策略提供了新选择;第三,祝融号可通过改变其质心位置,减小轮压,从而使其可以更好地面对异常路况;同时,主动悬架可在“祝融号”6个轮子的某一轮发生故障时,单独将该轮抬起,其余车轮正常行驶,继续完成火星探测任务。
2.2 沉陷机理分析火星车车轮是火星车直接与火壤接触的部分。研究火星车沉陷机理,首先要对车轮与火壤间的相互作用进行分析。正常行驶时,其轮/壤作用关系如图 2所示,其中:
当火星车正常行驶时,其沉陷深度应小于火星车设计时所规定的临界沉陷量,图 2所示即为火星车正常行驶时的车轮沉陷情况。当其沉陷量大于火星车设计时所规定的临界沉陷量时,可认为火星车出现沉陷异常的状况,如图 3所示。
火星车的异常沉陷主要由两部分组成,一部分是车轮所受纵向载荷
$ \begin{align} z_{0} = \Bigg[{\frac{3W}{ b \left({{k_{\rm c}} /{b}+k_{\psi}}\right)\sqrt{D}({3-n})}}\Bigg]^{ \frac{2}{2n+1}} \end{align} $ | (1) |
其中
由表 1可知,密实火壤的摩擦模量
根据表 1中数据、式(1) 以及“祝融号”火星车车轮及载荷参数可知,在不同土壤状态下其静态沉陷情况如表 2所示。
除静态沉陷外,另一部分沉陷是由于火星车在火星表面行驶时,火星车车轮对模拟火壤具有一定的切削作用从而产生的车轮滑转沉陷
$ \begin{align} z_{\rm s} =\frac{h_{\rm{gr}} \times s}{1-s} \end{align} $ | (2) |
其中
综上可得火星车在松软环境下总沉陷如下:
$ \begin{align} z=z_{0} +z_{\rm s} \end{align} $ | (3) |
另外,可从轮/壤作用关系的角度来分析产生沉陷的原因。由图 4(a)可知,对于火星车来说,当其匀速正常行驶时,轮部纵向载荷
根据Bekker和Janosi等的地面力学理论[16-17]可知,当土壤种类一定时,正应力和切应力随着沉陷深度的增加而逐渐增大。因此,当火星车行驶在松软火壤表面时,为了平衡车轮所受纵向载荷,火星车沉陷深度会大大增加,如图 4(b)所示。同时,当沉陷过大时火星车所受切向阻力也随之增大,因此火星车车轮上所施加的扭矩
由上文分析可知,当火星车行驶在火星表面时,由于各种原因,其车轮可能出现不同的沉陷情况。“祝融号”火星车的轮部接地压强临界值设计为5 kPa。其压强计算公式为
$ \begin{align} p=\frac{F}{S} \end{align} $ | (4) |
其中
根据式(4) 以及车轮尺寸计算可知,当接地压强为5 kPa时,其车轮沉陷深度为20.15 mm。综上可知“祝融号”火星车的临界沉陷量设计为20.15 mm。在此,对不同沉陷情况加以定义,如表 3所示。
普通沉陷对“祝融号”火星车正常巡视不会造成过大影响;当沉陷量达到临界沉陷值时,则认定“祝融号”火星车达到沉陷临界点;超出该点则进入极端沉陷状态,“祝融号”火星车需通过更改运动方式进行脱困。
由表 2中数据、临界沉陷量以及式(3) 可知,“祝融号”火星车在松软火壤环境下的最大滑转沉陷
由于静态沉陷为定值,因此判别“祝融号”火星车是否超过临界沉陷量的主要判据为其滑转沉陷量。由式(3) 可知,对于“祝融号”火星车来说,当其滑转沉陷量大于9.76 mm时,总的沉陷量大于所设计的临界沉陷量20.15 mm,可判断此时“祝融号”火星车沉陷量大于临界沉陷量。此外,由式(2) 可知火星车滑转沉陷量主要与火星车在行驶时的车轮滑转率有关,将滑转沉陷临界值以及“祝融号”火星车轮刺高度代入式(2) 可知,当滑转率大于66.1% 时,其滑转沉陷值大于9.76 mm,总沉陷大于“祝融号”火星车临界沉陷值。此时,可判断“祝融号”火星车陷入极端沉陷情况。
2.4 沉陷判别依据滑转率的大小可以作为判断“祝融号”火星车是否出现极端沉陷的重要判定依据。通过上文的理论分析可知,当行驶过程中滑转率大于66.1% 时,则可认为滑转沉陷大于9.76 mm,总沉陷量大于20.15 mm。此时,沉陷状态为极端沉陷,超出火星车设计范围,需更改运动方式。
但对于火星车来说,由于其信号传输时间较长,若通过这种方式进行测量,则易出现火星车极端沉陷值过大、无法脱困的情况。因而需要提出与滑转率相关的判断方式,使其可以更快作出判断。
在此,提出另一种运动异常状态判断方式,即通过驱动电流进行判别。当火星车陷入不同深度的火壤时,其前进所需驱动电流会随着下陷深度的变化而产生一定的变化。故可通过监测不同下陷深度下驱动电流的峰值,对其下陷深度进行判断。当驱动电流超过相应沉陷条件的峰值时,则可认为“祝融号”火星车出现沉陷量大于临界沉陷的情况,需更改运动方式。
2.5 脱困策略分析对于“祝融号”火星车来说,其脱困方法主要有2种。第1种方法通过增大驱动电流来增大轮部驱动力,从而克服由于沉陷等情况所产生的阻力。但当滑转率过大时,单纯增加驱动能力也无法实现脱困。
第2种方法通过改变运动方式进行脱困,主要指将正常轮式行驶方式改为蠕动式行驶进行脱困。其运动方式示意图如图 5所示。
在进行蠕动脱困时,其运动过程分为车体抬升以及车体下降2个阶段。在车体抬升阶段,火星车通过主副摇臂主动向上运动,使车前轮与中后轮轮间距变小,其中中轮和后轮通过电机驱动同步向前转动,前轮锁死不动。在车体下降阶段,火星车主副摇臂主动运动,导致火星车车体向下运动,此时中后轮锁死不动,前轮通过电机驱动向前转动,导致前轮与中后轮的轮间距变大。当中轮以及后轮锁死时,车轮受力
为了验证“祝融号”火星车沉陷机理与其脱困策略,进行了火星车沉陷及脱困试验。为了模拟火星表面的环境,在土槽中铺设模拟火壤。其中松软火壤为火星表面的飘浮物和沉积沙丘,该部分为火壤最松软状态,强度最小,粒径最细,是火星车最难通过的松软路面。因此,依据飘浮物火壤的参数,确定细粒径的中值粒径为40 μm,内聚强度为0.5 kPa,内摩擦角为18°,承压参数如表 1所示。
密实火壤中值粒径为200 μm,该中值粒径为火壤的平均值,能代表火星上大部分的松软沙地,其力学参数强于浮沙;内摩擦角(35°)和内聚强度(0.24 kPa)采用火壤的平均值,承压参数如表 1所示。
被试车为“祝融号”火星车初样鉴定件,其系统质量为90 kg。其结构增加主动悬架,使其具有蠕动脱困模式,蠕动模式是轮步复合式运动模式,可有效解决大滑转率情况下轮式移动无法通过的问题。试验环境及试验车如图 6所示。
在试验过程中,通过在“祝融号”火星车后施加一定的挂钩牵引力来使火星车在行驶过程中出现沉陷,通过改变挂钩牵引力数值来改变火星车沉陷量。首先,通过直接行驶的方式前进,当沉陷小于15 mm时则认定其脱困成功。若沉陷不变或沉陷继续增加则认定脱困失败。当直接行驶脱困失败时,则采用蠕动方式进行脱困,观察能否脱困成功。在试验过程中,记录火星车滑转率、沉陷深度、电流等数据。
4 结果与分析(Results and analysis) 4.1 不同沉陷深度下脱困方式分析为研究沉陷深度对脱困方式的影响,对车轮在不同沉陷深度情况下进行了试验,试验过程如图 7所示。图 8所示为沉陷深度为20 mm时,脱困前后的对比图。试验时,对不同沉陷深度下火星车脱困时间进行记录,得到不同沉陷情况下的脱困所需时间,如表 4所示。
由表 4可知,普通沉陷情况下,直接方式与蠕动方式均可使火星车成功脱困。其中直接方式脱困所需的时间远小于蠕动方式。同时发现,随着沉陷深度的增加,直接方式的脱困速度逐渐降低,这是由于随着沉陷深度的增加,火星车所受火壤阻力增大,导致其脱困速度降低。当沉陷量大于临界沉陷量时,使用直接方式无法使火星车成功脱困。通过分析发现,此时车轮所提供的驱动力无法克服火壤对车轮所施加的阻力,若继续增大驱动力则可能导致驱动电机烧毁,因此需改变运动方式进行脱困。由表 4可知,选用蠕动方式可使火星车在沉陷过大的情况下顺利脱困。此外,与直接脱困大体相同,随着沉陷深度的增加,其脱困效率也逐渐降低。
4.2 沉陷异常判据分析 4.2.1 滑转率判据在试验过程中分别记录不同沉陷深度条件下火星车在运动过程中的滑转率,绘制沉陷深度与滑转率之间的关系,如图 9所示。
由图 9,随着沉陷深度的增加,滑转率逐渐增大,且沉陷深度与滑转率呈现一定的函数关系。通过曲线拟合,可得沉陷深度与滑转率的关系:
$ \begin{align} s=88.86-\frac{722.99}{\sqrt{2\pi} z}{\rm e}^{ {\frac{\ln^{2} {\frac{z}{6.4}}}{2}}} \end{align} $ | (5) |
其中
在试验过程中分别记录不同下陷深度条件下火星车运动过程中的峰值电流情况,绘制下陷深度与电流之间的关系图,如图 10所示。
由图 10可知,随着下陷深度的增加,“祝融号”火星车前进所需驱动电流也逐渐增大,且二者呈现一定的函数关系。通过曲线拟合可得驱动电流峰值与火星车沉陷深度之间的关系式:
$ \begin{align} I=2.1-\frac{28.57}{z}{\rm e}^{\ln^{2} {\frac{z}{10.5}}} \end{align} $ | (6) |
其中
通过上文试验以及分析可制定“祝融号”火星车沉陷的判别标准。当火星车达到临界沉陷量20.15 mm时,则认为其无法正常前进。
基于试验结果,由拟合公式(5)(6) 可知,当“祝融号”火星车陷入极端沉陷情况时,对应的临界滑转率以及峰值电流分别大于61.25% 和1.9 A。
综上所述,可给出“祝融号”火星车脱困策略:在其行驶过程中,下陷深度、滑转率和峰值电流任一指标超过表 5中所列阈值时,需采用蠕动方式使其继续前进;当全部指标小于阈值且下陷深度小于15 mm时,可认定火星车脱离极端沉陷状态,可采用正常行驶方式前进。
另外,“祝融号”火星车的一个蠕动行驶周期为140 s,前进距离为318 mm。以此作为参考标准,当行驶时间为140 s且前进距离小于318 mm时,可认为火星车出现极端沉陷情况,需要采取蠕动方式脱困。此外,为了火星车安全性考虑,当上述条件同时发生时,即在表 6中序号5情况下,火星车不宜继续前进,应选择后退或转向。
1) 火星车在松软模拟火壤上产生的沉陷主要包括其自身载荷所产生的静态沉陷以及运动过程中滑转引起的滑转沉陷两部分。
2) 根据沉陷量的大小,可将火星车运行过程中所产生的沉陷分为普通沉陷、临界沉陷、极端沉陷3类。对于“祝融号”火星车,普通沉陷的沉陷量小于20.15 mm,临界沉陷的沉陷量为20.15 mm,极端沉陷的沉陷量大于20.15 mm。
3) 通过试验验证理论分析的准确性,得出了具体判别准则。当沉陷深度超过临界沉陷、滑转率大于61.25%、临界电流值大于1.90 A以及行驶318 mm时间大于140 s等4种情况发生任意一种时,应采取蠕动方式进行脱困。此外,为了保证“祝融号”火星车行驶安全,当上述条件同时发生时,应使“祝融号”火星车后退或转向。
[1] |
Poo M M. Mars exploration on the move[J]. National Science Review, 2020, 7(9): 1413. DOI:10.1093/nsr/nwaa181 |
[2] |
Biswal M M K, Das N B, Annavarapu R N. Orbital and plan etary challenges for human Mars exploration[DB/OL]. (2021-01-12)[2021-06-01]. https://arxiv.org/abs/2101.04725.
|
[3] |
Chen Y S, Han J J, Fan Y, et al. Overview of the latest scientific results of China's lunar exploration program[J]. Chinese Journal of Space Science, 2020, 40(5): 626-642. |
[4] |
司马光. 灵活高效的玉兔号月球车[J]. 国际太空, 2013(12): 15-19. Sima G. The flexible and efficient YuTu lunar rover[J]. International Space, 2013(12): 15-19. |
[5] |
Sebastián E, Martínez G, Ramos M, et al. Thermal calibration of the MEDA-TIRS radiometer onboard NASA's Perseverance rover[J]. Acta Astronautica, 2021, 182: 144-159. DOI:10.1016/j.actaastro.2021.02.006 |
[6] |
王琼, 贾阳, 陶灼, 等. 火星移动智能体技术探讨[J]. 航天器工程, 2015, 24(4): 27-32. Wang Q, Jia Y, Tao Z, et al. Discuss on Mars mobile agent technologies[J]. Spacecraft Engineering, 2015, 24(4): 27-32. DOI:10.3969/j.issn.1673-8748.2015.04.005 |
[7] |
Ding L, Deng Z Q, Gao H B, et al. Planetary rovers' wheelsoil interaction mechanics: New challenges and applications for wheeled mobile robots[J]. Intelligent Service Robotics, 2011, 4(1): 17-38. |
[8] |
黄晗, 李建桥, 陈百超, 等. 滑转条件下星球车坡面通过性评估试验[J]. 农业工程学报, 2016, 32(16): 40-44. Huang H, Li J Q, Chen B C, et al. Experiment of slope trafficability assessment of planetary rover under slip condition[J]. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(16): 40-44. DOI:10.11975/j.issn.1002-6819.2016.16.006 |
[9] |
Iagnemma K, Senatore C, Trease B P. Terramechanics modeling of Mars surface exploration rovers for simulation and parameter estimation[C]//ASME International Design Engineering Technical Conferences. New York, USA: ASME, 2011: 805-812.
|
[10] |
Johnson J B, Kulchitsky A V, Duvoy P, et al. Discrete element method simulations of Mars Exploration Rover wheel performance[J]. Journal of Terramechanics, 2015, 62: 31-40. |
[11] |
Heverly M, Matthews J, Lin J, et al. Traverse performance characterization for the Mars Science Laboratory rover[J]. Journal of Field Robotics, 2013, 30(6): 835-846. |
[12] |
Meirion-Griffith G, Spenko M. A modified pressure-sinkage model for small, rigid wheels on deformable terrains[J]. Journal of Terramechanics, 2011, 48(2): 149-155. |
[13] |
Lyasko M. Slip sinkage effect in soil-vehicle mechanics[J]. Journal of Terramechanics, 2010, 47(1): 21-31. |
[14] |
Huang H, Xu S C, Meng Z, et al. The sinkage characteristics and prediction of a planetary rover based on a similarity model experiment[J]. Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, 2019, 233(10): 3762-3774. |
[15] |
Richter L, Ellery A, Gao Y, et al. A predictive wheel-soil interaction model for planetary rovers validated in testbeds and against MER Mars rover performance data[C]//European Planetary Science Congress. 2006.
|
[16] |
Bekker M G. Theory of land locomotion[M]. Ann Arbor, USA: University of Michigan Press, 1956.
|
[17] |
Wong J Y. Terramechanics and off-road vehicle engineering: Terrain behaviour, off-road vehicle performance and design[M]. 2nd ed. Oxford, UK: Butterworth-Heinemann, 2009.
|