Longitudinal maneuver simulation of an X-51A-like aircraft based on numerical virtual flight
-
摘要:
吸气式高超声速飞行器在机动过程中,由于构型复杂,气动特性呈现强烈的非定常特性。传统基于数据库或气动力模型的飞行仿真不能准确描述机动过程中复杂的气动特性和运动规律。针对这一问题,基于现代软件分布式、模块化的发展趋势,建立了一个高效的数值虚拟飞行仿真平台。利用该平台,对一种类X-51A外形的吸气式高超声速飞行器开展了纵向机动闭环数值仿真,并与工程方法的结果进行了对比。研究发现:对于类X-51A外形的吸气式高超声速飞行器,在纵向拉起时,工程方法给出的结果可能不能完全反映非定常效应的影响。此时,应该采用更为精确的虚拟飞行方法来研究飞行器的闭环响应特性。此外,借助该仿真平台还研究了舵回路时间常数对控制系统的影响,为控制律设计提供了一定的参考。
Abstract:In the maneuvering process of air-breathing hypersonic vehicles, the aerodynamic characteristics show strong unsteady characteristics due to the complicated configuration. Traditional flight simulations based on database or aerodynamic models cannot exactly describe the complex aerodynamic characteristics and motion characteristics in the maneuvering process. To solve this problem, a numerical virtual flight simulation platform is established based on the distributed and modular development trend of modern software. With this coupled platform, longitudinal maneuvering simulations of an X-51A-like aircraft is carried out and the results are compared with those of the engineering method. It is found that For the X-51A-like air breathing hypersonic vehicle, the results given by the engineering method may not fully reflect the influence of unsteady effects when it pitches up longitudinally. In this case, the more accurate virtual flight method should be used to study the closed-loop response characteristics of the vehicle.In addition, the influence of the rudder loop time constant on the control system is also studied with the aid of the simulation platform, which provides a certain reference for the design of the control law.
-
我国已完成着陆器软着陆月球表面并利用自动探测机器人对月球进行探测[1]。作为航天大国,我国未来将陆续开展其他星球的探测工程[2]。腿式着陆器作为常用的着陆缓冲装置,对月球之外其他星球进行探索还可能用到该装置[3], 因此必须对腿式着陆器持续进行研究。为了检验腿式着陆器的软着陆性能并减少物理样机试验的成本,开展着陆模拟试验的动力学仿真显得尤为重要。由于着陆器的软着陆性能直接决定探测器能否安全抵达星球表面[4],对着陆器的软着陆性能进行优化成为设计阶段必不可少的环节。
文献[5-9]提出了着陆器的动力学分析方法,给出对软着陆性能影响较大的因素。文献[10]提出月球着陆缓冲器刚柔耦合模型的建模与分析方法。文献[11]运用响应面法,对单一工况下月球着陆器的软着陆性能进行了多目标优化。以上研究内容局限于着陆器单个性能指标或单一工况下多个性能指标的改进。相比于月球,其他星球的环境更加难以预测[12],着陆器各性能指标可能不在同一工况下达到最差,有必要在多组工况下,对着陆器多个性能指标进行优化。
本文建立了新型腿式着陆器的刚柔耦合动力学分析模型。探索着陆器最易翻倒、底面最易与星球表面岩石碰撞、主体承受最大碰撞力时对应的着陆工况。基于以上3组工况,以增强着陆器抗翻倒能力、降低着陆器主体底面与星球表面岩石碰撞的可能性、降低着陆器主体受力为目标进行优化设计。应用第二代非劣排序遗传算法 (NSGA-Ⅱ) 寻优,得到的优化结果明显提高了着陆器的软着陆性能。
1. 着陆器刚柔耦合动力学模型
1.1 新型腿式着陆器构型
本文研究的新型着陆器包括主体结构和4套对称分布的着陆缓冲机构。主体在动力学分析中作为负载,着陆缓冲机构用于吸收着陆器着陆过程产生的冲击力,每套缓冲机构由主支柱、辅助支柱、足垫和外部缓冲件组成,如图 1所示。
主支柱起支撑作用;辅助支柱由内、外筒和缓冲元件组成,内、外筒之间存在缓冲元件,用于吸收压缩载荷;足垫用来增大着陆器与星球土壤的接触面积,防止着陆时着陆腿过度下陷;外部缓冲件与主体固连,当主支柱受到较大冲击力时,外部缓冲件间接受力,发生弯曲变形,进而产生缓冲弯矩,吸收冲击能量,防止主体受力过大。
1.2 着陆缓冲机构力学模型
辅助支柱内部的缓冲元件是着陆器着陆过程吸收冲击能量的主要结构,元件只能承受压缩载荷。在辅助支柱内外筒之间施加与内、外筒相对位移 (D) 有关的力 (FD) 代替缓冲元件的作用效果,辅助腿缓冲机构示意图及D与FD关系曲线如图 2所示。
辅助支柱缓冲能力较强,外部缓冲件发生侧向弯曲的可能性很小。由于使用较多柔体构件导致动力学模型运算速度缓慢,使用如图 3方法对外部缓冲件进行刚体等效。
建立两段刚体模型,质量特性参照外部缓冲件柔体模型设置,刚体1绕转轴O转动,刚体2与刚体1只能沿轴线相对移动。刚体1转轴处施加与转角 (θ) 相关的转矩 (M),模拟外部缓冲件受弯时产生的缓冲弯矩;刚体2上施加与θ相关,沿刚体1轴线方向的运动 (S),模拟外部缓冲件自由端位置变化。在有限元软件中建立外部缓冲件模型,施加转矩,测得θ与M、θ与S关系曲线如图 4所示,输出数据文件,导入刚体模型。考虑等效模型质量、运动与力学特性均由有限元软件计算得到,即等效模型力学性能与柔体模型一致。
足垫触地时与星球表面的相互作用力分解为法向接触力Fn和切向摩擦力Ff,接触力Fn采用非线性阻尼弹簧碰撞方程描述:
(1) 式中:K为接触刚度;δ为接触变形;e为非线性指数;C为阻尼系数。以上系数均与星球土壤特性密切相关。
摩擦力Ff采用库伦摩擦方程描述:
(2) 式中:μ为足垫与星球表面间的摩擦系数。
1.3 刚柔耦合动力学模型建模
刚柔耦合多体系统中,构件既存在大范围的刚体运动,也存在局部的弹性形变[10]。着陆器动力学模型中,足垫作为柔性构件可发生弹性形变,其余构件均为刚体构件。
采用有限元软件Patran建立足垫有限元模型,导入多体动力学分析软件ADAMS中作为柔体构件;采用ADAMS软件建立其余构件的刚体模型;应用ADAMS软件对着陆器着陆过程的动力学问题进行求解。
2. 着陆器软着陆极限工况的确定
2.1 软着陆性能评判原则
根据着陆器安全着陆的性能指标,结合着陆器的构型,提出以下3种软着陆性能评判原则:
1) 着陆器在着陆过程中不发生翻倒,即着陆器质心距包含任意两足垫中心点的竖直翻倒墙平面的最小距离大于设计许用值[13]。
2) 着陆过程中辅助支柱外筒、外部缓冲件与主体连接点处的最大碰撞力不超过允许值,保证着陆器主体不会因为承受过大的碰撞力而受损。
3) 当着陆地点摩擦系数较小时,着陆器在着陆过程中,其着陆平台底面中心点距离星球表面的最小距离应大于预留间距值,避免着陆器底面与星球表面岩石发生碰撞而导致着陆器主体受损。
2.2 着陆工况参数
着陆工况的参数包括:式 (2) 中着陆器足垫与星球表面间的摩擦系数 (μ)、星球表面坡度 (α) 与凹坑分布情况 (n)、着陆器触地瞬间的竖直速度 (vx)、水平速度 (vz)、着陆偏航角 (θp)、着陆俯仰角 (θf) 以及着陆器整机质量 (m),如图 5所示。
着陆工况各参数的取值如表 1所示。
表 1 着陆工况各参数取值Table 1. Value of parameters of landing condition参数 取值 μ 0.3,0.8 α/(°) 7 n 0,1,2,3 vx/(m·s-1) 3.5 vz/(m·s-1) -1,+1 θp/(°) 0~45 θf/(°) -3,+3 m/kg 700,800 表 1中,n取0,1,2,3分别表示:着陆器没有着陆腿落入凹坑,2号与4号着陆腿落入坑中,1号与4号着陆腿落入坑中,1号、2号与4号着陆腿落入坑中;vz取正值表示水平速度方向背离坡面,负值表示水平速度朝向坡面;θf取正值表示着陆器俯仰角为仰角,负值表示俯仰角为俯角;遍历过程偏航角以1°为增量从0°增加至45°。将表 1中各参数取值进行组合,得到的着陆工况个数为2×1×4×1×2×46×2×2,即2 944种着陆工况。
2.3 软着陆极限工况参数取值
着陆器在特定工况下着陆,着陆过程中着陆器质心至翻倒墙距离取得最小值 (LT);或着陆过程中主体受碰撞力取得最大值 (FM);或着陆过程中着陆平台底面中心点至星球表面距离取得最小值 (LH),则称该工况为一种极限工况,对应的最值称为一种极限指标值。
基于2 944种着陆工况,通过仿真计算,得到3组极限工况及对应极限指标如表 2所示。
表 2 极限工况计算结果Table 2. Calculation results of the worst landing condition极限指标值 μ n vz/(m·s-1) θp/(°) θf/(°) m/kg LT=1.11 mm 0.8 2 1 9 3 700 LH=734.23 mm 0.3 3 -1 0 -3 800 FM=27.22 kN 0.8 3 -1 0 3 800 根据表 2结果可知:着陆器质心距翻倒墙最小距离低于许用值 (30 mm),该情况认为着陆器已经翻倒;星球表面摩擦系数较小时,着陆平台底面中心点距星球表面最小距离低于预留间距值 (760 mm),着陆器极易与星球表面岩石碰撞;由于辅助支柱缓冲能力较强,主体所受的最大碰撞力位于主体与外部缓冲件连接处,接近允许值 (30 kN)。为了提升着陆器的软着陆性能,下文利用多目标优化方法对缓冲机构进行优化设计。
3. 着陆器软着陆性能多目标优化
3.1 设计变量的选取
根据缓冲机构构型可知,增强外部缓冲件的缓冲力矩,可以提升着陆平台底面至星球表面距离的最小值,但主体受力可能超过许用值;增强辅助腿的缓冲力,可以提升着陆平台底面至星球表面距离的最小值,降低外部缓冲件与主体间的作用力,但缓冲性能降低导致着陆器更易翻倒;如图 6所示,增加足垫至着陆器中心轴线的距离 (dz),可使质心至翻倒墙初始距离增加,提升抗翻倒能力,但足垫受力点与外部缓冲件转轴水平距离 (d) 随之增加。根据力矩计算方程:
(3) 若着陆器足垫承受同样的冲击力 (F),足垫受力点与外部缓冲件转轴竖直距离h不变,冲击力与竖直方向夹角αF不变,d值增大,冲击转矩Mv增大,着陆腿向外扩张能力增强,进而导致着陆平台底面中心点至星球表面最小距离减小,着陆器底面更易与星球表面岩石相撞。
在着陆器主体结构与缓冲机构连接位置确定的情况下,同时避免着陆器竖直方向整体尺寸的增加,选取外部缓冲件缓冲转矩放大倍数 (kx)、辅助支柱缓冲力放大倍数 (kz)、足垫至着陆器中心轴线的距离 (dz) 作为设计变量进行多目标优化。
根据图 6可知,dz值增加导致着陆腿变长,着陆器整体尺寸随之增加。本文提到的缓冲元件提供的缓冲力与其尺寸成正比。为了保证着陆器整体尺寸不超出运载火箭可提供的最大安置空间,同时保证缓冲元件与着陆器其他结构不存在干涉,确定设计变量取值范围如表 3所示。
表 3 设计变量取值范围Table 3. Value range of design variables设计变量 取值范围 kx 1~2 kz 1~1.5 dz/m 1.37~1.45 3.2 软着陆性能多目标优化的实现
选取LT、LH、FM作为优化目标,约束LT、LH、FM满足2.3节给定的许用范围,结合设计变量的取值,建立数学模型如下:
(4) 式中:x=(kx,kz,dz)T代表设计变量;x(L)与x(U)分别为设计变量的下界与上界。
本文采用NSGA-Ⅱ进行优化计算,算法基于帕累托最优理论,具有迭代次数少、寻优能力强等特点[14-15]。算法参数设置如表 4所示[11]。
综上所述,着陆器软着陆性能多目标优化具体流程如图 7所示。
经过优化计算,得到帕累托最优解集如表 5所示。
表 5 帕累托最优解集Table 5. Pareto optimal solution aggregate编号 kz kx dz/m LH/mm FM/kN LT/mm 1 1.46 1.07 1.45 757.91 24.07 413.92 2 1.46 1.23 1.45 757.97 26.31 416.59 3 1.49 1.13 1.44 761.38 25.57 378.00 4 1.49 1.03 1.43 763.07 23.36 349.08 5 1.49 1.05 1.41 765.81 23.95 258.71 6 1.49 1.22 1.42 764.49 27.31 306.33 7 1.49 1.08 1.43 762.92 24.07 352.49 8 1.49 1.04 1.43 762.57 23.50 360.83 9 1.49 1.03 1.42 764.15 23.46 311.13 10 1.48 1.19 1.43 761.46 25.72 365.11 根据图 6可知,dz值增加导致着陆腿变长,着陆器整体尺寸随之增加,为了减小着陆器整体尺寸,降低着陆器在运载火箭中的安置空间,选取帕累托解集中dz值最小的结果作为最优解,重新进行2 944次仿真计算,获得极限指标值与优化前结果对比如表 6所示。
表 6 优化前后性能对比Table 6. Comparison of performance before and after optimization对比项 kx kz dz/m LT/mm LH/mm FM/kN 优化前 1 1 1.36 1.11 734.23 27.22 优化后 1.05 1.49 1.41 234.36 765.81 23.95 根据表 6可知,相比于初始设计,优化后的参数使LH增加4.2%,FM削弱12.1%,LT显著增加,着陆器在着陆过程中不再翻倒,软着陆性能获得明显提升。
4. 结论
1) 建立了新型腿式着陆器刚柔耦合动力学分析模型,实现着陆器软着陆过程的仿真,提出新型着陆器软着陆性能评判原则。
2) 根据着陆工况各参数的取值,确定2 944种着陆工况,针对所有工况进行仿真计算,得到着陆器着陆过程中最易受损的3组恶劣工况。
3) 应用NSGA-Ⅱ算法,基于仿真得到的3组恶劣工况,对着陆器的软着陆性能进行多目标优化。优化后,着陆器在2 944种工况下着陆不存在翻倒现象,着陆平台底面中心点与星球表面的最小距离大于预留间距值,且主体受力得到明显地减弱。
-
表 1 不同Td下控制系统的性能指标
Table 1. Performance indexes of control system under different Td
Td td tr tp ts(±5%) σ% 0 0.355 0.430 0.808 0.59 4.20 0.05 0.366 0.361 0.749 1.44 5.96 0.1 0.388 0.326 0.856 2.01 16.0 0.2 0.382 0.320 0.787 6.02 31.6 -
[1] HANK J M, MURPHY J S, MUTZMAN R C.The X-51A scramjet engine flight demonstration program: AIAA-2008-2540[R].Reston: AIAA, 2008. [2] 余华峰, 刘宏康, 陈树生, 等.类X-51A飞行器非定常湍流精细模拟[J].北京航空航天大学学报, 2019, 45(3):624-632. doi: 10.13700/j.bh.1001-5965.2018.0409YU H F, LIU H K, CHEN S S, et al.High-resolution unsteady turbulence simulation of an X-51A-like aircraft[J].Journal of Beijing University of Aeronautics and Astronautics, 2019, 45(3):624-632(in Chinese). doi: 10.13700/j.bh.1001-5965.2018.0409 [3] SALAS M D.Digital flight:The last CFD aeronautical grand challenge[J].Journal of Scientific Computing, 2006, 28(2-3):479-505. doi: 10.1007/s10915-006-9087-7 [4] CAI H M.Virtual flight simulation of a dual rotor micro air vehicle[J].International Journal of Computational Fluid Dynamics, 2015, 29(2):192-198. doi: 10.1080/10618562.2015.1016426 [5] CHEN Q, CHEN J Q, XIE Y F, et al.Study and application of virtual flight simulation for rolling control of vehicles[J].Journal of Computational Science, 2017, 21:77-85. doi: 10.1016/j.jocs.2017.05.009 [6] 张来平, 马戎, 常兴华, 等.虚拟飞行中气动、运动和控制耦合的数值模拟技术[J].力学进展, 2014, 44:376-417.ZHANG L P, MA R, CHANG X H, et al.Review of aerodynamics/kinematics/flight control coupling method in virtual flight simulations[J].Advances in Mechanics, 2014, 44:376-417(in Chinese). [7] 达兴亚, 陶洋, 赵忠良.基于预估校正和嵌套网格的虚拟飞行数值模拟[J].航空学报, 2012, 33(6):977-983.DA X Y, TAO Y, ZHAO Z L.Numerical simulation of virtual flight based on prediction-correction coupling method and chimera grid[J].Acta Aeronautica et Astronautica Sinica, 2012, 33(6):977-983(in Chinese). [8] ZHANG L P, CHANG X H, MA R, et al.A CFD-based numerical virtual flight simulator and its application in control law design of a maneuverable missile model[J].Chinese Journal of Aeronautics, 2019, 32(12):2577-2591. doi: 10.1016/j.cja.2019.07.001 [9] 陈琦, 郭勇颜, 谢昱飞, 等.PID控制器与CFD的耦合模拟技术研究及应用[J].航空学报, 2016, 37(8):2507-2516.CHEN Q, GUO Y Y, XIE Y F, et al.Research and application of coupled simulation techniques of PID controller and CFD[J].Acta Aeronautica et Astronautica Sinica, 2016, 37(8):2507-2516(in Chinese). [10] 李锋, 杨云军, 刘周, 等.飞行器气动/飞行/控制一体化耦合模拟技术[J].空气动力学学报, 2015, 33(2):156-161. doi: 10.7638/kqdlxxb-2014.0097LI F, YANG Y J, LIU Z, et al.Integrative simulation technique of coupled aerodynamics and flight dynamics with control law on a vehicle[J].Acta Aerodynamica Sinica, 2015, 33(2):156-161(in Chinese). doi: 10.7638/kqdlxxb-2014.0097 [11] MORTON S A, EYMANN T A, LAMBERSON S, et al.Relative motion simulations using an overset multi-mesh paradigm with kestrel v3[C]//Proceedings of the 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition.Reston: AIAA, 2012. [12] 张来平, 邓小刚, 张涵信.动网格生成技术及非定常计算方法进展综述[J].力学进展, 2010, 40:424-447. doi: 10.6052/1000-0992-2010-4-J2009-123ZHANG L P, DENG X G, ZHANG H X.Reviews of moving grid generation techniques and numerical methods for unsteady flow[J].Advances in Mechanics, 2010, 40:424-447(in Chinese). doi: 10.6052/1000-0992-2010-4-J2009-123 [13] CHEN Q, YUAN X X, WANG X, et al.Study of saddle-node bifurcation for longitudinal flight with CFD/RBD technique[J].Journal of Computational Science, 2018, 29:153-162. doi: 10.1016/j.jocs.2018.10.006 [14] CHATURVEDI N A, SANYAL A K, MCCLAMROCH N H.Rigid-body attitude control[J].Control Systems IEEE, 2011, 31(3):30-51. doi: 10.1109/MCS.2011.940459 [15] WALLNER E M, WELL K H.Attitude control of a reentry vehicle with internal dynamics[J].Journal of Guidance, Control, and Dynamics, 2003, 26(6):846-854. doi: 10.2514/2.6928 [16] 席柯, 袁武, 阎超, 等.基于闭环控制的带翼导弹虚拟飞行数值模拟[J].航空学报, 2014, 35(3):634-642.XI K, YUAN W, YAN C, et al.Virtual flight numerical simulation of the basic finner projectile with closed loop[J].Acta Aeronautica et Astronautica Sinica, 2014, 35(3):634-642(in Chinese). [17] WANG W, YAN C, WANG S, et al.An efficient, robust and automatic overlapping grid assembly approach for partitioned multi-block structured grids[J].Proceedings of the Institution of Mechanical Engineers, Part G:Journal of Aerospace Engineering, 2018, 233(4):1217-1236. [18] ZHONG K, YAN C, CHEN S S, et al.Numerical study on the aerothermodynamics of different heatshield configurations for Mars entry capsules[J].Acta Astronautica, 2019, 157:189-198. doi: 10.1016/j.actaastro.2018.12.025 [19] 王吉阳, 杨学良.远程过程调用的设计与实现[J].计算机工程与设计, 1991(2):71-79.WANG J Y, YANG X L.The design and implementation of remote procedure call[J].Computer Engineering and Design, 1991(2):71-79(in Chinese). [20] BISWAS R, LU X Y, PANDA D K.Accelerating tensor flow with adaptive RDMA-based gRPC[C]//Proceedings of the 25th International Conference on High Performance Computing.Piscataway: IEEE Press, 2018: 2-11. [21] 刘君, 刘瑜, 陈泽栋.非结构变形网格和离散几何守恒律[J].航空学报, 2016, 37(8):2395-2407.LIU J, LIU Y, CHEN Z D.Unstructured deforming mesh and discrete geometric conservation law[J].Acta Aeronautica et Astronautica Sinica, 2016, 37(8):2395-2407(in Chinese). [22] 李跃军.飞行器强迫运动和自由飞的非定常流动数值模拟[D].北京: 北京航空航天大学, 2007: 56-63.LI Y J.Numerical simulation of unsteady flows around aircrafts experiencing forced motion and free flight[D].Beijing: Beihang University, 2007: 56-63(in Chinese). [23] 岳瑞华, 徐中英, 周涛.导弹控制原理[M].北京:北京航空航天大学出版社, 2016:17-41.YUE R H, XU Z Y, ZHOU T.Control principle of missiles[M].Beijing:Beihang University Press, 2016:17-41(in Chinese). [24] 吴森堂.飞行控制系统[M].2版.北京:北京航空航天大学出版社, 2013:180-191.WU S T.Flight control system[M].2nd ed.Beijing:Beihang University Press, 2013:180-191(in Chinese). [25] HALL L H, PARTHASARATHY V.Validation of an automated Chimera/6-DOF methodology for multiple moving body problems: AIAA-1998-0767[R].Reston: AIAA, 1998. [26] 刘雪松.X-51A高超声速飞行器三维重建及气动/隐身特性分析[D].南京: 南京航空航天大学, 2015: 9-18.LIU X S.Three-dimensional reconstruction and analysis aerodynamic and stealth characteristic of X-51A hypersonic vehicle[D].Nanjing: Nanjing University of Aeronautics and Astronautics, 2015: 9-18(in Chinese). 期刊类型引用(11)
1. 钟世英,刘文昊,凌道盛,岳茜茜. 足垫构型对着陆器软着陆稳定性影响研究. 计算机仿真. 2024(12): 81-87 . 百度学术
2. 陈含智,孙蕊,邱明,毛继志,胡浩亮,张立东. 基于自适应噪声方差的卫星定位故障检测法. 北京航空航天大学学报. 2023(02): 406-421 . 本站查看
3. 杨智杰,王刚,赵瑞杰,王春洁,赵军鹏. 火星进入舱配平翼机构展开冲击动力学分析. 北京航空航天大学学报. 2023(02): 422-429 . 本站查看
4. Bo LEI,Ming ZHANG,Hanyu LIN,Hong NIE. Optimization design containing dimension and buffer parameters of landing legs for reusable landing vehicle. Chinese Journal of Aeronautics. 2022(03): 234-249 . 必应学术
5. 贾山,赵建华,胡汝洁,陈金宝,周向华,张胜. 可复用小型月表着陆器设计优化及仿真分析. 宇航学报. 2022(03): 356-365 . 百度学术
6. Zongmao Ding,Hongyu Wu,Chunjie Wang,Jianzhong Ding. Hierarchical Optimization of Landing Performance for Lander with Adaptive Landing Gear. Chinese Journal of Mechanical Engineering. 2019(02): 24-35 . 必应学术
7. 关俊山,杨大鹏,张艳兵,王进军. 基于折衷规划的客车电池箱结构多工况拓扑优化. 武汉理工大学学报(交通科学与工程版). 2019(02): 363-367 . 百度学术
8. 丁宗茂,王春洁,吴宏宇,丁建中,郭璠. 探测器触地关机软着陆稳定性分析. 北京航空航天大学学报. 2018(03): 614-620 . 本站查看
9. 刘学翱,吴宏宇,王春洁,丁宗茂,丁建中. 着陆器变阻尼缓冲器性能分析及参数优化. 北京航空航天大学学报. 2018(10): 2149-2155 . 本站查看
10. 吴宏宇,王春洁,丁宗茂,丁建中,董洋,满剑锋. 着陆姿态不确定下的着陆器缓冲机构优化设计. 宇航学报. 2018(12): 1323-1331 . 百度学术
11. 吴宏宇,王春洁,丁宗茂,丁建中,刘学翱. 两种着陆模式下的着陆器缓冲机构构型优化. 宇航学报. 2017(10): 1032-1040 . 百度学术
其他类型引用(3)
-