CN112597641A - 一种运载器着陆稳定性优化方法 - Google Patents

一种运载器着陆稳定性优化方法 Download PDF

Info

Publication number
CN112597641A
CN112597641A CN202011455291.5A CN202011455291A CN112597641A CN 112597641 A CN112597641 A CN 112597641A CN 202011455291 A CN202011455291 A CN 202011455291A CN 112597641 A CN112597641 A CN 112597641A
Authority
CN
China
Prior art keywords
landing
carrier
design
buffer system
parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011455291.5A
Other languages
English (en)
Other versions
CN112597641B (zh
Inventor
张华�
刘汉武
王栋梁
丁秀峰
林剑锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Aerospace System Engineering Institute
Original Assignee
Shanghai Aerospace System Engineering Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Aerospace System Engineering Institute filed Critical Shanghai Aerospace System Engineering Institute
Priority to CN202011455291.5A priority Critical patent/CN112597641B/zh
Publication of CN112597641A publication Critical patent/CN112597641A/zh
Application granted granted Critical
Publication of CN112597641B publication Critical patent/CN112597641B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明一种运载器着陆稳定性优化方法,涉及运载器支腿式着陆缓冲技术领域。本发明以运载器可垂直起降的支腿式着陆缓冲系统为研究对象,通过将运动系统之间复杂的交互作用采用数学模型等效处理,并将数学模型中的设计参数实现参数化建模,然后针对经过参数化处理的结构设计参数进行优化设计,获取影响着陆稳定性的关键因素。基于获取的关键因素物理量特点,确立其概率分布准则,进行随机打靶分析,评估关键因素对着陆稳定性影响的权重,解决运载器支腿式着陆缓冲系统的结构优化设计问题和无法开展地面有效试验验证的设计难题。

Description

一种运载器着陆稳定性优化方法
技术领域
本发明涉及运载器的着陆缓冲技术领域,尤其涉及一种运载器着陆稳定性优化方法,解决运载器支腿式着陆缓冲系统的结构优化设计问题和无法开展地面有效试验验证的设计难题,通过快速、有效地动力学分析、迭代优化、参数评估,能极大地缩短运载器支腿式着陆缓冲系统的设计研发周期。
背景技术
运载器着陆缓冲系统在返回地面前距地面一定高度时展开,通过与地面的接触碰撞,利用着陆缓冲系统内部的耗能部件结构实现运载器稳定着陆。运载器在着陆前,会受地面着陆场风载荷,空气阻力,地面摩擦阻力,着陆水平速度、垂直速度、滚转角速度、偏航角速度、俯仰角速度及横向初始姿态等因素影响,这使得运载器支腿式着陆缓冲系统在着陆过程中会出现多种情况:单腿着陆、两腿着陆及四腿同时着陆。上述每种着陆情况,运载器的着陆姿态、着陆缓冲系统的承载有很大差异,有些初始条件会使得运载器出现翻转,导致着陆失败。
为解决上述运载器着陆过程中的稳定性问题,有必要进行运载器着陆稳定性分析及优化设计。
发明内容
本发明的目的在于:提供一种运载器着陆稳定性优化方法,为运载器着陆缓冲系统快速、高效的迭代设计提供方法支撑和理论指导,解决运载器着陆缓冲系统的着陆稳定性问题,保证运载器可重复使用领域的工程实现。
本发明的技术方案是:一种运载器着陆稳定性分析与优化设计方法,步骤如下:
步骤一:建立运载器着陆缓冲系统结构的几何模型,并对结构设计参数中的关注变量实现参数化;根据运载器着陆缓冲系统结构的工作原理,建立连接参数、约束参数、运动参数及力元素等效模型,并对所关注的设计参数变量实现参数化;
步骤二:基于步骤一建立的参数化几何模型、参数化数学等效模型,建立参数化的着陆缓冲系统动力学模型;
步骤三:基于步骤二建立的动力学模型,针对所有可能会影响着陆性能的设计变量,进行动力学分析和优化设计,确定主要关键影响参数,并结合其物理量意义或特点建立其参数概率分布准则;
步骤四:基于步骤二建立的动力学模型,结合步骤三关键影响参数的概率分布准则,进行随机打靶分析,并基于分析结果,综合评估主要关键参数对着陆缓冲性能的影响权重。
进一步的,步骤一中将运载器本体、着陆缓冲系统中的部件均视作刚体,建立参数化几何模型。运载器着陆缓冲系统由运载器本体和着陆缓冲系统通过转动副、转动副连接。着陆缓冲系统由脚撑、主支柱、辅助支撑连接,且沿运载器本体周向正交布局,主支柱与辅助支撑由球铰副连接,脚撑与主支柱由球铰副连接。
进一步的,步骤一中运动参数包含运载器本体横向姿态偏差Θ,着陆滚转角速度ωx、偏航角速度ωy、俯仰角速度ωz,着陆垂直速度Vx、水平速度Vy,水平速度Vz。力元素参数包含地面着陆场风载荷W,主支柱内部缓冲器缓冲力f、摩擦力ff、耗能部件fw
进一步的,步骤一中主支柱包含油气结构和耗能部件,用于吸收运载器着陆缓冲过程中的能量。依据运载器着陆缓冲系统连接关系及着陆工作原理,对所有部件间的连接、约束、运动参数及力元素采用数学等效模型建模并对所有参数变量实现参数化,结合几何模型一起搭建整个运载着陆缓冲系统的动力学模型。
进一步的,步骤二中以运载器本体质心姿态θ、w、v及脚撑与地面的接触碰撞力Fc作为设计目标,对建立的动力学模型进行单因素的动力学分析和优化设计,确定所有参数变量对设计目标的敏感度,获取影响设计目标的主要关键因素。上述关键因素中与长度相关量、摩擦系数一般假设为均匀分布,与力元素相关量一般假设为正态分布且置信度一般设为3σ。
进一步的,步骤四中基于步骤二建立的动力学模型,结合步骤三参数概率分布情况,进行多工况(单腿着陆、两腿着陆、四腿着陆)、多参数、多目标(姿态角、着陆缓冲碰撞力)的动力学随机打靶分析。根据分析结果计算各个参数化变量与优化目标之间的相关性,并对相关性进行降序或升序排列,综合评估主要关键参数对着陆缓冲性能的影响权重。
本发明的有益效果:本发明为运载器着陆缓冲系统快速、高效的迭代设计提供方法支撑和理论指导,解决运载器着陆缓冲系统的着陆稳定性问题,保证运载器可重复使用领域的工程实现。本发明以运载器可垂直起降的支腿式着陆缓冲系统为研究对象,通过将运动系统之间复杂的交互作用采用数学模型等效处理,并将数学模型中的设计参数实现参数化建模,然后针对经过参数化处理的结构设计参数进行优化设计,获取影响着陆稳定性的关键因素。基于获取的关键因素物理量特点,确立其概率分布准则,进行随机打靶分析,评估关键因素对着陆稳定性影响的权重,解决运载器支腿式着陆缓冲系统的结构优化设计问题和无法开展地面有效试验验证的设计难题。
附图说明
图1为本发明的稳定性分析与优化设计方法流程图;
图2为本发明的运载器着陆缓冲系统示意图,其中2(a)为侧视图,2(b)为俯视图;
图3为本发明的运载器着陆缓冲系统与地面相互作用工况示意图;
图4为本发明的主支柱耗能部件位移-载荷作用曲线示意图;
图5为本发明的脚撑与地面接触碰撞力变化曲线示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
一种运载器着陆稳定性分析与优化设计方法的流程见图1,包括以下步骤:
步骤一:建立运载器着陆缓冲系统结构的几何模型,并对结构设计参数中的预先关注变量实现参数化;运载器着陆缓冲系统由运载器本体1和着陆缓冲系统通过转动副5、转动副6连接。着陆缓冲系统由脚撑4、主支柱2、辅助支撑3连接,且沿运载器本体周向正交布局,主支柱2与辅助支撑3由球铰副连接,脚撑4与主支柱2由球铰副连接。主支柱2包含油气结构和耗能部件,用于吸收运载器着陆缓冲过程中的能量,见图2。根据运载器着陆缓冲系统结构的工作原理,建立连接参数、约束参数、运动参数及力元素等效模型,并对所预先所关注的设计参数变量实现参数化;依据运载器着陆缓冲系统连接关系及着陆工作原理示意图,将运载器本体1、着陆缓冲系统中的部件均视作刚体,所有设计变量实现参数化;运载器本体1与着陆缓冲系统的连接约束5、6以及主支柱2与辅助支撑3的连接约束7,均采用数学等效模型,约束三个方向平动自由度,释放三个转动自由度;脚撑4与主支柱2的连接约束采用固定副,约束所有的平动、转动自由度。主支柱与运载器本体的连接点(PUp_X,PUp_Y,PUp_Z)、辅助支撑与运载器本体连接点(PMid_X,PMid_Y,PMid_Z)、主支柱与脚撑连接点(PDown_X,PDown_Y,PDown_Z)、主支柱外筒半径及长度(ROuter,LOuter)、主支柱活塞杆半径及长度(RInner,LInner)、辅助支撑半径及长度(RFuzhu,LFuzhu)、脚撑半径RFoot、油气最大行程Smax等几何模型参数实现变量参数化。
运载器本体横向姿态偏差Θ、着陆滚转角速度ωx、偏航角速度ωy、俯仰角速度ωz、着陆垂直速度Vx、水平速度Vy,水平速度Vz等运动参数实现变量参数化,运载器与地面的相互作用方式分为三种:四腿着陆、两腿着陆、单腿着陆,见图3。
步骤二:基于步骤一建立的参数化几何模型、参数化数学等效模型,建立参数化的着陆缓冲系统动力学模型;运载器着陆缓冲过程中,施加力元素如当地重力加速度a,方向为竖直向下,着陆场风载荷Fw,作用于运载器本体质心,方向垂直于运载器轴线。主支柱内部缓冲器力学模型可分解为四个部分:
f=fh+fa+ff+fs+fw (1)
其中,fs为缓冲器结构限制力,fa为缓冲器弹性力,fh为缓冲器油液阻尼力,ff为缓冲器内部摩擦力,fw为耗能部件缓冲力。
结构限制力fs可以由下式表达:
Figure BDA0002828556100000051
其中,KS缓冲器轴向拉压刚度;S0为缓冲器的初始行程,Smax为缓冲器的最大压缩行程。
缓冲器弹性力fa可以由下式表达:
Figure BDA0002828556100000052
其中,V0为缓冲器空腔初始体积;Aa为有效压气面积;N为气体多变系数;Patm为大气常压。
缓冲器油液阻尼力fh可以由下式表达:
Figure BDA0002828556100000053
其中,ρ为油液密度;Ah为压油面积;A+为正行程油孔面积;A-为反行程油孔面积;Cd为油孔的缩流系数。
缓冲器内部摩擦力ff可以由下式表达:
Figure BDA0002828556100000061
其中,μm为摩擦系数。
耗能部件缓冲力fw为单向过程,通过耗能部件压缩吸收能耗,被压缩后不可恢复逆转,其位移与载荷曲线见图4。
步骤三:基于步骤二建立的动力学模型,针对所有预估可能会影响着陆性能的设计变量,进行动力学分析和优化设计,确定主要关键影响参数,并结合其物理量意义或特点获得的影响运载器着陆稳定性的关键参数概率分布准则,与长度相关量、摩擦系数一般假设为均匀分布,与力元素相关量一般假设为正态分布且置信度一般设为3σ。
步骤四:基于步骤二建立的动力学模型,结合步骤三关键影响参数的概率分布准则,进行随机打靶分析,并基于分析结果,综合评估主要关键参数对着陆缓冲性能的影响权重。采用式(6)分析参化变量与优化目标(如脚撑与地面接触碰撞力见图5)之间的相关性,并对相关性进行降序或升序排列,综合评估主要关键参数对着陆缓冲性能的影响权重。
Figure BDA0002828556100000062
式中xi为参数化变量,yi为优化设计目标变量,
Figure BDA0002828556100000063
均为变量的加权平均值。
通过本发明,可实现运载器着陆缓冲系统的稳定性分析与优化设计,本发明的方法已运用于我国可垂直起降运载器的缓冲着陆系统的结构设计中,为运载器着陆缓冲系统的实施应用奠定了基础。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (7)

1.一种运载器着陆稳定性优化方法,其特征在于包括以下步骤:
步骤一:建立运载器着陆缓冲系统结构的几何模型,并对结构设计参数中的预先关注变量实现参数化;根据运载器着陆缓冲系统结构的工作原理,建立连接参数、约束参数、运动参数及力元素等效模型,并对所预先所关注的设计参数变量实现参数化;
步骤二:基于步骤一建立的参数化几何模型、参数化数学等效模型,建立参数化的着陆缓冲系统动力学模型;
步骤三:基于步骤二建立的动力学模型,针对所有预估可能会影响着陆性能的设计变量,进行动力学分析和优化设计,确定主要关键影响参数,并结合其物理量意义或特点建立其参数概率分布准则;
步骤四:基于步骤二建立的动力学模型,结合步骤三关键影响参数的概率分布准则,进行随机打靶分析,并基于分析结果,综合评估主要关键参数对着陆缓冲性能的影响权重。
2.如权利要求1所述的一种运载器着陆稳定性优化方法,其特征在于,步骤一中所述的几何模型包括运载器几何模型和着陆缓冲系统几何模型,步骤一中所述的参数化指对所述模型的设计赋予参数变量;
所述的运载器几何模型包含质量、惯量、外包络长度及直径;
所述的着陆缓冲系统几何模型包含支腿张开角度、主支柱套筒及级数、辅着陆腿关节、油液缓冲器、压溃式缓冲器。
3.如权利要求1所述的一种运载器着陆稳定性优化方法,其特征在于,在步骤一几何模型基础上,根据结构部件间的实际运动状态和工作原理,建立连接、运动、约束及外力。
4.如权利要求1所述的一种运载器着陆稳定性优化方法,其特征在于,步骤二中建立的动力学模型的所有参数变量完全实现参数化,根据初始条件的变化实现模型自动更新。
5.如权利要求1所述的一种运载器着陆稳定性优化方法,其特征在于,步骤三中对所有设计参数变量进行单因素动力学分析和优化设计,评估每个因素对运载器着陆性能的影响程度,以此确定影响着陆稳定性的关键因素。
6.如权利要求1所述的一种运载器着陆稳定性优化方法,其特征在于,步骤四中的随机打靶分析,分析内容包括多参数、多目标优化设计。
7.如权利要求1所述的一种运载器着陆稳定性优化方法,其特征在于,步骤四中所述的评估影响权重指关键影响因素对运载器着陆姿态及在该条件下与地面接触碰撞力大小的影响程度的相对比例。
CN202011455291.5A 2020-12-10 2020-12-10 一种运载器着陆稳定性优化方法 Active CN112597641B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011455291.5A CN112597641B (zh) 2020-12-10 2020-12-10 一种运载器着陆稳定性优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011455291.5A CN112597641B (zh) 2020-12-10 2020-12-10 一种运载器着陆稳定性优化方法

Publications (2)

Publication Number Publication Date
CN112597641A true CN112597641A (zh) 2021-04-02
CN112597641B CN112597641B (zh) 2022-07-22

Family

ID=75192568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011455291.5A Active CN112597641B (zh) 2020-12-10 2020-12-10 一种运载器着陆稳定性优化方法

Country Status (1)

Country Link
CN (1) CN112597641B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110098835A1 (en) * 2009-10-23 2011-04-28 Yuecel Sermet Parametric configurator for product design: system and method
US20140107990A1 (en) * 2011-03-17 2014-04-17 Messler-Dowty Limited METHOD AND SYSTEM FOR DETERMINING FRICTION COEFFICIENT μ FOR AN AIRCRAFT LANDING EVENT
CN103778290A (zh) * 2014-01-15 2014-05-07 北京理工大学 一种腿式着陆器着陆缓冲机构的建模方法
US8725470B1 (en) * 2010-05-17 2014-05-13 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) Co-optimization of blunt body shapes for moving vehicles
US20140214370A1 (en) * 2013-01-30 2014-07-31 Honda Research Institute Europe Gmbh Optimizing the design of physical structures/objects
CN109783953A (zh) * 2019-01-22 2019-05-21 中国人民解放军空军工程大学 一种新型飞机着陆动荷载的计算方法
CN110442117A (zh) * 2019-08-26 2019-11-12 北京理工大学 一种火星探测器大底气动耦合分离过程安全性分析方法
CN111783232A (zh) * 2020-07-23 2020-10-16 中国人民解放军国防科技大学 一种基于聚类分析的可回收火箭返回段自适应优化方法
CN112001029A (zh) * 2020-07-28 2020-11-27 清华大学 一种基于凸优化的火箭在线轨迹优化定制化求解器

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110098835A1 (en) * 2009-10-23 2011-04-28 Yuecel Sermet Parametric configurator for product design: system and method
US8725470B1 (en) * 2010-05-17 2014-05-13 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) Co-optimization of blunt body shapes for moving vehicles
US20140107990A1 (en) * 2011-03-17 2014-04-17 Messler-Dowty Limited METHOD AND SYSTEM FOR DETERMINING FRICTION COEFFICIENT μ FOR AN AIRCRAFT LANDING EVENT
US20140214370A1 (en) * 2013-01-30 2014-07-31 Honda Research Institute Europe Gmbh Optimizing the design of physical structures/objects
CN103778290A (zh) * 2014-01-15 2014-05-07 北京理工大学 一种腿式着陆器着陆缓冲机构的建模方法
CN109783953A (zh) * 2019-01-22 2019-05-21 中国人民解放军空军工程大学 一种新型飞机着陆动荷载的计算方法
CN110442117A (zh) * 2019-08-26 2019-11-12 北京理工大学 一种火星探测器大底气动耦合分离过程安全性分析方法
CN111783232A (zh) * 2020-07-23 2020-10-16 中国人民解放军国防科技大学 一种基于聚类分析的可回收火箭返回段自适应优化方法
CN112001029A (zh) * 2020-07-28 2020-11-27 清华大学 一种基于凸优化的火箭在线轨迹优化定制化求解器

Also Published As

Publication number Publication date
CN112597641B (zh) 2022-07-22

Similar Documents

Publication Publication Date Title
CN105912772A (zh) 一种集频率、振型和频响函数为一体的模型修正方法
Yue et al. Design and landing dynamic analysis of reusable landing leg for a near-space manned capsule
Wu et al. Robust optimization design method for powertrain mounting systems based on six sigma quality control criteria
CN112009577A (zh) 一种用于重卡驾驶室半主动悬置的控制方法
Li et al. A new dynamic modelling methodology of a force measuring system for hypersonic impulse wind tunnel
CN112597641B (zh) 一种运载器着陆稳定性优化方法
Xu et al. Semi-active control of a new quasi-zero stiffness air suspension for commercial vehicles based on event-triggered H∞ dynamic output feedback
CN117113725A (zh) 参数摆式波浪能装置的获能分析方法
Legeza Dynamics of vibration isolation system with a quasi-isochronous roller shock absorber
Rong1a et al. A novel nonlinear gas-spring TMD for the seismic vibration control of a MDOF structure
Lin et al. Dynamic modeling and parameter identification of a gun saddle ring
Zhang et al. Modal and natural frequency sensitivity analysis of electrohydraulic Stewart platform
CN112793805B (zh) 全机落震缩比模型试验方法
Gao et al. Programmed load spectrum for fatigue bench test of a vehicle body
Song et al. Aiming error analysis of guns in ground combat vehicles operating on bumpy roads
CN108647383A (zh) 一种结构增强调谐质量阻尼器优化设计方法
Shangguan et al. Dynamic analysis and design calculation methods for powertrain mounting systems
Andersen Multibody dynamics modeling and system identification for a quarter-car test rig with McPherson strut suspension
Martynowicz et al. Wind turbine tower-nacelle system with MR tuned vibration absorber: modelling, test rig, and identification
Wang et al. Dynamics analysis of a missile vehicle considering the pavement roughness
Preissner et al. Nostradamus and the synchrotron engineer: Key aspects of predicting accelerator structural response
Zhou et al. Identification of Planar Double-wishbone Suspension Mechanism Using Jacobian Approach
Wang et al. Analysis of Vibration Characteristics of Electro-hydraulic 3-UPS/S Parallel Stabilized Platform
Motahari et al. Aerodynamic Load Maps of Bluff-Body Combinations in Incompressible Flow
Zhang et al. The modeling and analysis of vibration response with airdrop vehicle in landing process

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant