CN103778290A - 一种腿式着陆器着陆缓冲机构的建模方法 - Google Patents
一种腿式着陆器着陆缓冲机构的建模方法 Download PDFInfo
- Publication number
- CN103778290A CN103778290A CN201410018124.2A CN201410018124A CN103778290A CN 103778290 A CN103778290 A CN 103778290A CN 201410018124 A CN201410018124 A CN 201410018124A CN 103778290 A CN103778290 A CN 103778290A
- Authority
- CN
- China
- Prior art keywords
- urceolus
- outer tube
- inner core
- motion
- point
- 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
Links
Images
Landscapes
- Vibration Dampers (AREA)
Abstract
本发明涉及一种腿式着陆器着陆缓冲机构的建模方法，属于腿式星球着陆器技术领域。本发明方法首先建立缓冲着陆腿构形；然后建立铝蜂窝材料等效模型，应力应变关系等效转化成内外筒之间的载荷行程关系，采用连接单元模拟载荷行程关系；最后建立缓冲着陆腿的连接关系，其优点在于采用有限元建模，充分考虑了结构柔性对响应的影响；考虑了缓冲着陆腿柔性对着陆腿内外筒之间的相互作用；对铝蜂窝力学特性的简化建模，在保证准确度的前提下有效地提高了计算效率。
Description
技术领域
本发明涉及一种腿式着陆器着陆缓冲机构的建模方法，属于腿式星球着陆器技术领域。
背景技术
腿式着陆器通常使用一次性缓冲器，该缓冲器具有工作可靠性高、结构简单、质量轻等诸多优势，广泛应用于星球探测领域。我国开展的探月工程项目即使用腿式着陆器，腿式着陆器动力学仿真对着陆器设计具有极大的指导意义；而腿式着陆器动力学仿真建模的一项关键技术便是对缓冲着陆腿的合理建模。
传统的缓冲着陆腿建模基于多刚体动力学模型，模型比较简单，计算效率高，但没能考虑结构的柔性，在进行着陆动力学分析时不能准确反映柔性对结构响应的影响；另一类模型采用部件模态综合方法，通过有限元分析软件生成着陆腿的模态中性文件，并导入多体动力学软件中，建立柔性模型。但该建模方法由于引入浮动坐标系，在数值计算中，大范围运动常常会吞噬掉柔性体的小变形运动，因而数值精度比较差。还有一类模型，采用非线性有限元方法，使用可压溃单元模拟缓冲器的动力特性，该方法全面考虑了结构的柔性，但计算耗时巨大。
近年来，国外开发了诸多有限元分析软件，可以通过调用该类软件实施本发明的方法。例如：ABAQUS是国际先进的大型通用有限元计算分析软件，可以进行复杂的动力学分析计算。
发明内容
本发明的目的是为解决数值精度差等问题，提供一种腿式着陆器着陆缓冲机构的建模方法，能够很好反应缓冲腿的力学特性，同时考虑结构的柔性影响，计算耗时小、效率高。
一种腿式着陆器着陆缓冲机构的建模方法，具体步骤如下：
步骤一，建立缓冲着陆腿构形。
缓冲着陆腿由内筒、外筒以及铝蜂窝芯子构成，内外筒之间能相对运动。内筒同轴安装于外筒内，内外筒沿轴向平行安装，安装后外筒远离内筒的一端称为外筒上端，内筒位于外筒内的一端称为内筒上端，铝蜂窝芯子填满于外筒上端面、外筒内壁以及内筒上端面形成的中空部分。内外筒的几何特性参数包括：内筒的长度lin、内筒的外半径Rin、内筒壁厚度hin、外筒的长度lout、外筒的外半径Rout、外筒壁厚度hout以及内筒安装到外筒内后缓冲着陆腿的总长度lw，为满足安装要求其中Rout-hout≥Rin，lin+lout≥lw。
内外筒采用壳单元模拟，赋予内外筒各自的截面属性，并对内外筒分别进行网格划分。截面属性包含壳单元的材料属性与壳单元的厚度，材料属性为杨氏模量E与材料密度ρ，壳单元的厚度值即为内外筒壁的厚度hin和hout。
步骤二，建立铝蜂窝材料等效模型。
铝蜂窝材料具有特殊的应力应变σ～ε关系，其中σ表示应力，ε表示应变，应力应变关系由对铝蜂窝材料进行应力应变实验获得。应力应变关系等效转化成内外筒之间的载荷行程关系为：当内外筒之间相对运动距离s，内外筒之间会产生沿轴向方向阻碍其相对运动的载荷力F。相对运动距离s＝lwε，载荷力
本方法采用连接单元模拟上述载荷行程关系，连接单元建立过程为：先在缓冲着陆腿两端部分别建立参考点RA、RB，RB点处于内筒下端面圆心处（内筒位于外筒外的一端称为内筒下端），RA点处于外筒上端面圆心处。在RB点建立随体坐标系，随体坐标系中只保留RA点沿轴线方向的自由度。接着在RA、RB之间建立载荷行程关系。具体设置为：连接单元的塑性屈服力设置为载荷力F，连接单元的塑性变形值设置为相对运动距离s，连接单元的弹性刚度设置为铝蜂窝材料的杨氏模量Eal。
步骤三，建立缓冲着陆腿的连接关系。
通过设置摩擦系数μ表示摩擦力，在内筒外壁与外筒内壁之间建立接触属性。将步骤二中建立的RA点和RB点通过多点约束方式分别约束固定到外筒上端面和内筒下端面，使得RA点的运动与外筒上端面运动一致，RB点的运动与内筒下端面的运动一致。
至此，缓冲着陆腿模型建立完成。在下一步着陆器整体分析中，只需将RA点和RB点分别与着陆器主体以及足垫再次按照多点约束的方式进行连接。
有益效果
本发明方法与现有技术相比，具有以下优点：
（1）采用有限元建模，充分考虑了结构柔性对响应的影响；
（2）考虑了缓冲着陆腿柔性对着陆腿内外筒之间的相互作用；
（3）对铝蜂窝力学特性的简化建模，在保证准确度的前提下有效地提高了计算效率。
附图说明
图1为本发明的缓冲着陆腿结构示意图；
图2为具体实施方式中缓冲着陆腿连接单元示意图；
图3为具体实施方式中缓冲着陆腿载荷行程示意图；
图4为具体实施方式中着陆器整体示意图；
标号说明：1-外筒，2-铝蜂窝芯子，3-内筒，4-缓冲腿，5-着陆器主体，6-足垫。
具体实施方式
下面将结合附图对本发明作进一步的详细说明。
步骤一、缓冲着陆腿构形建立
按照内外筒的几何特性参数建立几何模型，缓冲着陆腿结构示意图如图1所示。几何参数包括：内筒的长度lin＝875mm、内筒的外半径Rin＝49mm、内筒壁厚度hin＝2.5mm、外筒的长度lout＝1424mm、外筒的外半径Rout＝52.5mm、外筒壁厚度hout＝2.5mm以及内筒安装到外筒内后的总长度lw＝1880mm。
内外筒采用壳单元模拟，赋予内外筒各自的截面属性，并对内外筒分别进行网格划分。截面属性包含壳单元的材料属性与壳单元的厚度，材料属性为杨氏模量E与材料密度ρ，壳单元的厚度值即为内外筒壁的厚度hin＝2.5mm和hout＝2.5mm。
步骤二、铝蜂窝材料缓冲特性模型建立
采用连接单元模拟上述载荷行程关系，连接单元建立过程为：先在缓冲着陆腿两端部建立参考点RA、RB，如图2所示，RB点处于内筒下端面圆心处，RA点处于外筒上端面圆心处，在RB点建立随体坐标系，在随体坐标系中只保留RA点沿轴线方向的自由度，即图中x方向的自由度。
接着在RA、RB之间建立载荷行程关系。载荷行程关系如图3所示。在ABAQUS平台下，可以通过非线性移动连接单元Translator实现，连接单元的各级塑性屈服力YieldForce设置为载荷力F1与F2，连接单元的塑性变形值PlasticMotion设置为相对运动量S1与S2，连接单元的弹性刚度Elasticity设置为铝蜂窝材料的杨氏模量Eal。
步骤三、缓冲着陆腿各部分之间连接关系的建立
内筒外壁与外筒内壁之间建立接触，通过设置摩擦系数μ＝0.4表达筒壁之间的摩擦特性。将步骤二中建立的RA点和RB点通过多点约束方式分别约束固定到步骤一所建立的外筒上端面和内筒下端面，在ABAQUS平台下可通过Coupling建立多点约束，使得RA点的运动与外筒上端面的运动一致，RB点的运动与内筒下端面的运动一致。
将各个缓冲腿结构与着陆器主体以及足垫连接，如图4所示，并建立相应的边界及初始条件，便可以进行着陆动力学仿真分析。该建模方法在保证一定准确度的前提下节省了计算耗时，特别适用于多工况着陆仿真分析。
Claims (4)
1.一种腿式着陆器着陆缓冲机构的建模方法，其特征在于：具体步骤如下：
步骤一，建立缓冲着陆腿构形；
缓冲着陆腿由内筒、外筒以及铝蜂窝芯子构成，内外筒之间能相对运动；内筒同轴安装于外筒内，内外筒沿轴向平行安装，安装后外筒远离内筒的一端称为外筒上端，内筒位于外筒内的一端称为内筒上端，内筒位于外筒外的一端称为内筒下端；铝蜂窝芯子填满于外筒上端面、外筒内壁以及内筒上端面形成的中空部分；内外筒的几何特性参数包括：内筒的长度lin、内筒的外半径Rin、内筒壁厚度hin、外筒的长度lout、外筒的外半径Rout、外筒壁厚度hout以及内筒安装到外筒内后缓冲着陆腿的总长度lw，其中Rout-hout≥Rin，lin+lout≥lw；
内外筒分别采用壳单元模拟，赋予内外筒各自的截面属性，并对内外筒分别进行网格划分；截面属性包含壳单元的材料属性与壳单元的厚度，材料属性为杨氏模量E与材料密度ρ，壳单元的厚度值为内外筒壁的厚度hin和hout；
步骤二，建立铝蜂窝材料等效模型；
铝蜂窝材料的应力应变关系等效转化成内外筒之间的载荷行程关系为：当内外筒之间相对运动距离s，内外筒之间产生沿轴向方向阻碍其相对运动的载荷力F；相对运动距离s＝lwε，载荷力
其中σ表示应力，ε表示应变；
采用连接单元模拟载荷行程关系，连接单元建立过程为：在缓冲着陆腿两端分别建立参考点RA、RB，RB点处于内筒下端面圆心处，RA点处于外筒上端面圆心处；在RB点建立随体坐标系，随体坐标系中只保留RA点沿轴线方向的自由度；接着在RA、RB之间建立载荷行程关系；具体设置为：连接单元的塑性屈服力设置为载荷力F，连接单元的塑性变形值设置为相对运动距离s，连接单元的弹性刚度设置为铝蜂窝材料的杨氏模量Eal；
步骤三，建立缓冲着陆腿的连接关系；
通过设置摩擦系数μ表示摩擦力，在内筒外壁与外筒内壁之间建立接触属性；将步骤二中建立的RA点和RB点通过多点约束方式分别约束固定到外筒上端面和内筒下端面，使得RA点的运动与外筒上端面运动一致，RB点的运动与内筒下端面的运动一致；
至此，缓冲着陆腿模型建立完成。
2.根据权利要求1所述的一种腿式着陆器着陆缓冲机构的建模方法，其特征在于：应力应变σ～ε关系由对铝蜂窝材料进行应力应变实验获得。
3.根据权利要求1所述的一种腿式着陆器着陆缓冲机构的建模方法，其特征在于：在ABAQUS平台下通过Coupling建立多点约束，使得RA点的运动与外筒上端面的运动一致，RB点的运动与内筒下端面的运动一致。
4.根据权利要求1所述的一种腿式着陆器着陆缓冲机构的建模方法，其特征在于：在ABAQUS平台下，通过非线性移动连接单元Translator实现连接单元模拟载荷行程关系。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410018124.2A CN103778290B (zh) | 2014-01-15 | 2014-01-15 | 一种腿式着陆器着陆缓冲机构的建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410018124.2A CN103778290B (zh) | 2014-01-15 | 2014-01-15 | 一种腿式着陆器着陆缓冲机构的建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103778290A true CN103778290A (zh) | 2014-05-07 |
| CN103778290B CN103778290B (zh) | 2015-02-25 |
Family
ID=50570520
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410018124.2A Expired - Fee Related CN103778290B (zh) | 2014-01-15 | 2014-01-15 | 一种腿式着陆器着陆缓冲机构的建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103778290B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112597641A (zh) * | 2020-12-10 | 2021-04-02 | 上海宇航系统工程研究所 | 一种运载器着陆稳定性优化方法 |
| CN113609725A (zh) * | 2021-07-24 | 2021-11-05 | 上海大学 | 一种着陆缓冲机构的刚柔耦合动力学仿真方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100010782A1 (en) * | 2008-07-09 | 2010-01-14 | Livermore Software Technology Corporation | Systems and Methods of Selecting a CAE Analysis Solver with Appropriate Numerical Precision in Each of a Series of Hierarchically Related Engineering Simulations |
| CN102981507A (zh) * | 2012-11-29 | 2013-03-20 | 北京理工大学 | 一种软着陆自主障碍规避常推力器控制方法 |
-
2014
- 2014-01-15 CN CN201410018124.2A patent/CN103778290B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100010782A1 (en) * | 2008-07-09 | 2010-01-14 | Livermore Software Technology Corporation | Systems and Methods of Selecting a CAE Analysis Solver with Appropriate Numerical Precision in Each of a Series of Hierarchically Related Engineering Simulations |
| CN102981507A (zh) * | 2012-11-29 | 2013-03-20 | 北京理工大学 | 一种软着陆自主障碍规避常推力器控制方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112597641A (zh) * | 2020-12-10 | 2021-04-02 | 上海宇航系统工程研究所 | 一种运载器着陆稳定性优化方法 |
| CN113609725A (zh) * | 2021-07-24 | 2021-11-05 | 上海大学 | 一种着陆缓冲机构的刚柔耦合动力学仿真方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103778290B (zh) | 2015-02-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN203982627U (zh) | 多功能桁架结构实验装置 | |
| CN103425833B (zh) | 一种基于熵格子波尔兹曼模型的并行流体计算实现方法 | |
| Zhang et al. | Hierarchical anisogrid stiffened composite panel subjected to blast loading: Equivalent theory | |
| Dang et al. | Accelerated loosely-coupled CFD/CSD method for nonlinear static aeroelasticity analysis | |
| CN103488812A (zh) | 卫星在轨微振动仿真建模修正方法 | |
| CN103778290B (zh) | 一种腿式着陆器着陆缓冲机构的建模方法 | |
| CN103407516B (zh) | 一种用于月面上升的具有姿态调整功能的腿式缓冲装置 | |
| CN105279325A (zh) | 考虑整体空间作用的钢管混凝土锥形柱等效计算长度及稳定性的计算方法 | |
| Zhao et al. | Investigation on sloshing and vibration mitigation of water storage tank of AP1000 | |
| Cheng et al. | A second-order cell-centered Lagrangian method for two-dimensional elastic-plastic flows | |
| Zhou et al. | Improving the buffer energy absorption characteristics of movable lander-numerical and experimental studies | |
| CN104537205B (zh) | 一种被动约束阻尼旋转体结构的振动分析方法 | |
| CN106909755B (zh) | 一种车辆荷载下沉管隧道管节竖向位移的计算方法 | |
| Xiao et al. | Julia language in computational mechanics: A new competitor | |
| CN103955578A (zh) | 一种直升机全机抗坠毁分级仿真方法 | |
| Duan et al. | Partition method for impact dynamics of flexible multibody systems based on contact constraint | |
| CN111859738B (zh) | 一种暴雨环境下管道多物理场作用致灾机理演示方法 | |
| Lihua et al. | Dynamic analysis of aircraft arresting gear based on finite element method | |
| CN106951585B (zh) | 一种弹塑性弯曲变形限力结构的建模方法 | |
| Kang et al. | Dynamic analysis of communication transmission tower based on finite element method | |
| Yi et al. | Analysis on Seismic Performance of Steel Reinforced High Strength Concrete Columns | |
| Wang et al. | The Rigid-Flexible Coupling Dynamic Condition Analysis of New Type of Pipe Racking Manipulator | |
| Chen et al. | Efficient multidisciplinary aerodynamic optimization design based on discrete adjoint method | |
| Portioli et al. | Three-dimensional limit analysis of dry-joint masonry walls using rigid block modelling | |
| Qu et al. | Energy analysis of water impact of an elastic cylindrical shell |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150225 Termination date: 20220115 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |