CN110287575A - 一种空间采样回转机构的扭矩有限元仿真方法 - Google Patents

一种空间采样回转机构的扭矩有限元仿真方法 Download PDF

Info

Publication number
CN110287575A
CN110287575A CN201910535855.7A CN201910535855A CN110287575A CN 110287575 A CN110287575 A CN 110287575A CN 201910535855 A CN201910535855 A CN 201910535855A CN 110287575 A CN110287575 A CN 110287575A
Authority
CN
China
Prior art keywords
swing mechanism
standard shaft
spatial sampling
finite element
rigid
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
CN201910535855.7A
Other languages
English (en)
Other versions
CN110287575B (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.)
Shenyang Institute of Automation of CAS
Original Assignee
Shenyang Institute of Automation of CAS
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 Shenyang Institute of Automation of CAS filed Critical Shenyang Institute of Automation of CAS
Priority to CN201910535855.7A priority Critical patent/CN110287575B/zh
Publication of CN110287575A publication Critical patent/CN110287575A/zh
Application granted granted Critical
Publication of CN110287575B publication Critical patent/CN110287575B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

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

Abstract

本发明涉及一种空间采样回转机构的扭矩有限元仿真方法,建立采样器的初始模型,并将初始模型导入HYPERMESH进行网格划分,并设置刚性构件和柔性构件;在LSDYNA中,根据DEM离散元方法建立土壤模型;在采样器初始模型的基轴转轴质心位置建立刚性转轴,使刚性转轴与基轴同轴线;在LSDYNA中,采用同轴线的弹簧阻尼单元连接刚性转轴与基轴,设置弹簧阻尼单元的弹簧属性和扭转属性,设置工况参数,计算得出仿真结果。本发明可实现对对称回转机构采样的优化设计,针对微重力环境中空间采样机构进行过程仿真,提高了机械的可靠性,确保了任务的完成,给空间采样机构的设计提供了重要参考。

Description

一种空间采样回转机构的扭矩有限元仿真方法
技术领域
本发明涉及深空探测仿真分析技术领域,具体地说是一种空间采样回转机构的扭矩有限元仿真方法。
背景技术
深空探测是研究行星演化和形成规律的重要手段。为了了解太阳系的起源与演变过程,小行星的样本是蕴藏这些信息的重要载体。小行星是46亿年前太阳系初期形成的行星体,内部蕴含了大量太阳系早期形成和演化的信息,同时小行星含有丰富的贵金属及稀有元素,具有巨大的利用价值。
行星表面基本都是微重力环境,在小行星表面实行长期采样无法实现,采样过程所产生的反作用力可能会将采样机器推离小行星,这一特征使得合理的设计,优化机械结构是任务成功的必要保证。另外,小行星距地球遥远,且目前的技术水平限制了探测器所带燃料的质量,所以取样系统还应轻型化和小型化。
空间中的机械器件能否正常运行是空间任务成功的关键所在。由于将器件带上空间的成本巨大,而且对其重量的控制也相对严格,所以在设计时,合理的型号选择显得十分重要。电机是机器的重要动力来源,电机型号的选择对于空间探测任务至关重要。由于无重力的太空实验环境难以在地球模拟,空间中的运动预测主要通过运动仿真的方法进行验证。现阶段的扭矩仿真多为线性动力学方法,求解过程简化,对于结果的精度有所影响。扭矩的非线性动力学仿真方法能够提供更可靠设计建议。
发明内容
针对现有技术的不足,本发明提供一种空间采样回转机构的扭矩有限元仿真方法,用于预测在微重力的太空环境中回转机构的扭矩,作为空间装备设计的重要设计指标与优化目标。
本发明为实现上述目的所采用的技术方案是:
一种空间采样回转机构的扭矩有限元仿真方法,
步骤1:建立采样器的初始模型,并将初始模型导入HYPERMESH进行网格划分,并设置刚性构件和柔性构件;
步骤2:在LSDYNA中,根据DEM离散元方法建立土壤模型;
步骤3:在采样器初始模型的基轴转轴质心位置建立刚性转轴,使刚性转轴与基轴同轴线;
步骤4:在LSDYNA中,采用同轴线的弹簧阻尼单元连接刚性转轴与基轴,设置弹簧阻尼单元的弹簧属性和扭转属性,设置工况参数,计算得出仿真结果。
所述采样器的初始模型为对称回转机构,包括基轴和毛刷,基轴为两根,在基轴外表面上均布圆柱体毛刷,且两根基轴之间的毛刷交错排布。
所述土壤模型为根据月壤的材料参数建立,粒子大小范围为4~25mm。
所述刚性转轴采用壳单元建立,刚性转轴的直径小于基轴直径,位移与基轴同步,转速为预设采样转速。
所述刚性构件为基轴。
所述柔性构件为毛刷。
所述工况参数包括:下降速度,刚性转轴的转速,重力加速度,材料属性,厚度属性。
本发明具有以下有益效果及优点:
1.本发明提出了一种针对空间采样机构的非线性动力学扭矩仿真方法,给空间采样机构的设计提供了重要参考。
2.本发明所提的方法可实现对对称回转机构采样的优化设计。
3.本发明针对微重力环境中空间采样机构进行过程仿真,提高了机械的可靠性,确保了任务的完成。
附图说明
图1是本发明的方法流程图;
图2是空间采样回转机构的简化三维几何模型的等轴侧结构示意图;
图3是空间采样回转机构的有限元模型;
图4是空间采样回转机构的弹簧单元连接形式示意图;
其中,1为采样区三维模型、2为回转基轴、3为采样毛刷、4为土壤离散元节点、5为刚性转轴、6为弹簧单元。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但本发明能够以很多不同于在此描述的其他方式来实施,本领域技术人员可以在不违背发明内涵的情况下做类似改进,因此本发明不受下面公开的具体实施的限制。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。
如图1所示为本发明的方法流程图。
方法包括以下步骤:
1)将初始模型导入HYPERMESH进行网格划分与初始设置,建立有效的土壤模型,合理设置刚性构件。
2)在采样回转机构的回转轴线质心处建立一刚性转轴,约束其位移与转速。
3)应用弹簧阻尼单元将刚性转轴与回转零件连接,给弹簧单元赋予相应的section与材料,最后在LSDYNA中进行仿真求解。
所述土壤模型采用DEM离散元方法,建立土壤模型,采用无反射边界条件与微重力设置模拟太空环境。将回转机构的转轴设置为刚性体,与其相连的采样刷毛设置为柔性体。
定义的约束,对回转零件与刚性转轴施加相同的位移约束,只对刚性转轴施加转动约束。
弹簧阻尼单元特性为弹簧特性,弹簧只有扭转特性,其材料特性设置为线性,其大小可以带动回转基轴,并且基轴转速波动较小。
如图2所示为本发明的空间采样回转机构的简化三维几何模型的等轴侧结构示意图。
图中1为采样区三维模型、2为回转基轴、3为采样毛刷。
图3是空间采样回转机构的有限元模型。
在本实例中模拟土壤材料为月壤,重力环境为1/6g,利用LSDYNA软件模拟月球采样机构的运动过程与扭矩求解。
对采样回转机构进行模型简化,将采样毛刷3与回转基轴2均简化为壳网格,划分网格。其中回转基轴2设置为刚性体,采样毛刷3设置为柔性体,采用关键字*CONSTRAINED_EXTRA_NODES_SET进行连接。4为土壤离散元节点。
对于采样区土壤模型的建立采用DEM离散元模型,DEM模型的粒子大小可设置为4到25mm,本实例采用统一直径大小为4mm的粒子。材料应用LSDYNA内置的土壤材料模型*MAT_SOLD_AND_FOAM,设置材料参数包括密度、剪切模量等,应用关键字*CONTROL_DISCRETE_ELEMENT定义粒子间的刚度系数、摩擦系数、阻尼系数等,以上参数可通过指导手册查询和仪器测量等途径获得。对于土壤模型的边界条件采用无反射边界条件,此方法可以模拟土壤的无限域的分析情况。
如图4所示为空间采样回转机构的弹簧单元连接形式示意图。
其中5为刚性转轴、6为弹簧单元。在回转基轴的质心位置建立一刚性转轴5,使回转基轴2轴线与刚性转轴5的回转轴线共线。对于约束的施加,应用*BOUNDARY_PRESCRIBED_MOTION_RIGID关键字对回转基轴2与建立的刚性转轴5施加相同的位移约束,只对刚建立的刚性转轴施加回转约束。
在刚性转轴5与回转基轴2间采用弹簧单元6进形连接,刚性转轴5的转动通过弹簧单元5传递给回转基轴2,弹簧单元6的力矩即为采样过程的力矩需求。最后设置弹簧单元的相关输出选项,在LSDYNA中进行仿真求解。

Claims (7)

1.一种空间采样回转机构的扭矩有限元仿真方法,其特征在于:
步骤1:建立采样器的初始模型,并将初始模型导入HYPERMESH进行网格划分,并设置刚性构件和柔性构件;
步骤2:在LSDYNA中,根据DEM离散元方法建立土壤模型;
步骤3:在采样器初始模型的基轴转轴质心位置建立刚性转轴,使刚性转轴与基轴同轴线;
步骤4:在LSDYNA中,采用同轴线的弹簧阻尼单元连接刚性转轴与基轴,设置弹簧阻尼单元的弹簧属性和扭转属性,设置工况参数,计算得出仿真结果。
2.根据权利要求1所述的空间采样回转机构的扭矩有限元仿真方法,其特征在于:所述采样器的初始模型为对称回转机构,包括基轴和毛刷,基轴为两根,在基轴外表面上均布圆柱体毛刷,且两根基轴之间的毛刷交错排布。
3.根据权利要求1所述的空间采样回转机构的扭矩有限元仿真方法,其特征在于:所述土壤模型为根据月壤的材料参数建立,粒子大小范围为4~25mm。
4.根据权利要求1所述的空间采样回转机构的扭矩有限元仿真方法,其特征在于:所述刚性转轴采用壳单元建立,刚性转轴的直径小于基轴直径,位移与基轴同步,转速为预设采样转速。
5.根据权利要求1所述的空间采样回转机构的扭矩有限元仿真方法,其特征在于:所述刚性构件为基轴。
6.根据权利要求1所述的空间采样回转机构的扭矩有限元仿真方法,其特征在于:所述柔性构件为毛刷。
7.根据权利要求1所述的空间采样回转机构的扭矩有限元仿真方法,其特征在于:所述工况参数包括:下降速度,刚性转轴的转速,重力加速度,材料属性,厚度属性。
CN201910535855.7A 2019-06-20 2019-06-20 一种空间采样回转机构的扭矩有限元仿真方法 Expired - Fee Related CN110287575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910535855.7A CN110287575B (zh) 2019-06-20 2019-06-20 一种空间采样回转机构的扭矩有限元仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910535855.7A CN110287575B (zh) 2019-06-20 2019-06-20 一种空间采样回转机构的扭矩有限元仿真方法

Publications (2)

Publication Number Publication Date
CN110287575A true CN110287575A (zh) 2019-09-27
CN110287575B CN110287575B (zh) 2021-08-31

Family

ID=68004882

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910535855.7A Expired - Fee Related CN110287575B (zh) 2019-06-20 2019-06-20 一种空间采样回转机构的扭矩有限元仿真方法

Country Status (1)

Country Link
CN (1) CN110287575B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111473901A (zh) * 2020-05-12 2020-07-31 中国科学院沈阳自动化研究所 一种回转式对称采样刷机构的扭矩测试装置
CN113946914A (zh) * 2020-07-16 2022-01-18 中国科学院沈阳自动化研究所 一种回转式小行星采样装置的扭矩动力学仿真方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201123795Y (zh) * 2007-12-14 2008-10-01 朱秦颎 一种拖把及与其配套的甩干脱水装置
CN104155141A (zh) * 2014-07-29 2014-11-19 北京空间飞行器总体设计部 一体式行星表面土壤采集装置
CN105608297A (zh) * 2016-02-17 2016-05-25 中国科学院沈阳自动化研究所 一种角接触球轴承动刚度的数值模拟方法
CN205946880U (zh) * 2016-08-29 2017-02-15 凯里市鸿睿畜禽养殖合作社有限公司 便于清粪的鸡舍
CN109015255A (zh) * 2018-08-13 2018-12-18 张文英 一种抛光机砂辊

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201123795Y (zh) * 2007-12-14 2008-10-01 朱秦颎 一种拖把及与其配套的甩干脱水装置
CN104155141A (zh) * 2014-07-29 2014-11-19 北京空间飞行器总体设计部 一体式行星表面土壤采集装置
CN105608297A (zh) * 2016-02-17 2016-05-25 中国科学院沈阳自动化研究所 一种角接触球轴承动刚度的数值模拟方法
CN205946880U (zh) * 2016-08-29 2017-02-15 凯里市鸿睿畜禽养殖合作社有限公司 便于清粪的鸡舍
CN109015255A (zh) * 2018-08-13 2018-12-18 张文英 一种抛光机砂辊

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHRISTOPH M. SCHÄFER 等: "Numerical simulations of regolith sampling processes", 《PLANETARY AND SPACE SCIENCE》 *
ROBERT BONITZ 等: "The Brush Wheel Sampler-a Sampling Device for Small-body Touch-and-Go Missions", 《2012 IEEE AEROSPACE CONFERENCE》 *
庞彧 等: "月面钻取式自动采样机构的设计与分析", 《中国空间科学技术》 *
桑楠 等: "矩形截面弹簧应力及变形分析", 《国外建材科技》 *
谢宇明 等: "月壤的钻取采样离散元动态行为研究", 《工程设计学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111473901A (zh) * 2020-05-12 2020-07-31 中国科学院沈阳自动化研究所 一种回转式对称采样刷机构的扭矩测试装置
CN113946914A (zh) * 2020-07-16 2022-01-18 中国科学院沈阳自动化研究所 一种回转式小行星采样装置的扭矩动力学仿真方法
CN113946914B (zh) * 2020-07-16 2024-04-12 中国科学院沈阳自动化研究所 一种回转式小行星采样装置的扭矩动力学仿真方法

Also Published As

Publication number Publication date
CN110287575B (zh) 2021-08-31

Similar Documents

Publication Publication Date Title
Vršnak et al. Influence of the aerodynamic drag on the motion of interplanetary ejecta
Dobbs-Dixon et al. Spin-orbit evolution of short-period planets
Taktakishvili et al. Validation of the coronal mass ejection predictions at the Earth orbit estimated by ENLIL heliosphere cone model
Jia et al. MESSENGER observations and global simulations of highly compressed magnetosphere events at Mercury
CN103473417B (zh) 一种针对轮体结构伞形振动的颗粒阻尼减振设计方法
CN110287575A (zh) 一种空间采样回转机构的扭矩有限元仿真方法
Usmanov et al. A tilted‐dipole MHD model of the solar corona and solar wind
Lavraud et al. Asymmetry of magnetosheath flows and magnetopause shape during low Alfvén Mach number solar wind
Yao et al. Experimental study on broadband bistable energy harvester with L-shaped piezoelectric cantilever beam
Harne et al. Prospects for nonlinear energy harvesting systems designed near the elastic stability limit when driven by colored noise
Baker et al. Solar wind forcing at Mercury: WSA‐ENLIL model results
Furuya et al. Relativistic turning acceleration of radiation belt electrons by whistler mode chorus
McKee The dynamical structure and evolution of giant molecular clouds
Samsonov et al. Do we know the actual magnetopause position for typical solar wind conditions?
Ferrari et al. N-body gravitational and contact dynamics for asteroid aggregation
Chen et al. The cosmic neutrino background as a collection of fluids in large-scale structure simulations
Kerin et al. Mountain formation by repeated, inhomogeneous crustal failure in a neutron star
Cook et al. Magnetic braking and viscous damping of differential rotation in cylindrical stars
Fowler et al. MAVEN observations of low frequency steepened magnetosonic waves and associated heating of the Martian nightside ionosphere
Caballero‐Lopez et al. Galactic cosmic ray modulation: effects of the solar wind termination shock and the heliosheath
Taktakishvili et al. Model uncertainties in predictions of arrival of coronal mass ejections at Earth orbit
Wang et al. A new dynamic fluid‐kinetic model for plasma transport within the plasmasphere
Li et al. Numerical considerations in the simulation of equatorial spread F
Wang et al. Aeolian creep transport: Theory and experiment
Escultura The mathematics of the grand unified theory

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210831

CF01 Termination of patent right due to non-payment of annual fee