CN103678822A - Mechanical environment prediction method of lunar probe soft landing impact - Google Patents

Mechanical environment prediction method of lunar probe soft landing impact Download PDF

Info

Publication number
CN103678822A
CN103678822A CN201310727234.1A CN201310727234A CN103678822A CN 103678822 A CN103678822 A CN 103678822A CN 201310727234 A CN201310727234 A CN 201310727234A CN 103678822 A CN103678822 A CN 103678822A
Authority
CN
China
Prior art keywords
detector
environment
landing
lunar
spectrum
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
CN201310727234.1A
Other languages
Chinese (zh)
Other versions
CN103678822B (en
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201310727234.1A priority Critical patent/CN103678822B/en
Publication of CN103678822A publication Critical patent/CN103678822A/en
Application granted granted Critical
Publication of CN103678822B publication Critical patent/CN103678822B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention relates to a mechanical environment prediction method of lunar probe soft landing impact, and belongs to the technical field of deep space exploration. According to the method, a non-linear finite element model of probe lunar surface soft landing is established, model reduction is conducted on a probe according to the generalized dynamic condensation method, and therefore the purpose of improving the solving efficiency is achieved; then, non-linear finite element solving is conducted on the reduced model to obtain a dynamic response of the structure of the probe; after an acceleration speed response at a concerted position is obtained, the shock spectrum of the acceleration speed response of the structure is calculated according to a recursion digital filtering method, wherein the shock spectrum is used for describing the dynamic environment of probe landing impact; at last, a maximum expected environment is obtained by the adoption of a spectrum enveloping method, wherein the maximum expected environment is used for setting the designing and testing conditions of equipment. By means of the mechanical environment prediction method, the solving efficiency is improved and the product designing period is shortened.

Description

Moon detector in flexible landing ballistics environment predicting method
Technical field
The present invention relates to a kind of moon detector in flexible landing ballistics environment predicting method, belong to survey of deep space technical field.
Background technology
For realizing lunar surface soft landing detection mission, lunar orbiter system will experience two large processes: flight landing process and the lunar surface course of work.According to current achievement in research, the moon does not have geology activity substantially, so after safe landing, except a small amount of disturbance that each mechanism kinematic of detector system produces, without typical mechanical environment.So the main mechanical environment of detector system is all embodied in flight landing process.Flight landing process can be divided into again carrier rocket transmitter section, month transfer leg, ring month section and four parts of power descending branch.Wherein first three section and conventional satellite, particularly the mechanical environment of ring month explorer satellite experience is basically identical, mainly comprise the mechanical environment that transmitter section is transmitted by carrier rocket, the noise circumstance of endoatmosphere flight, and the detector shock environment that engine ignition produces while becoming rail.
Power descending branch is the distinctive mission phase of soft lunar landing task, and at its latter end, detector can be from the height free-falling of about 4m to moonscape.Can, in landing shock process, detector need to experience a transient impact load, hold out against the key that this mechanical environment is soft landing success.Therefore be necessary just this environment to be made to indication in the schematic design phase, to formulate design and the test condition of equipment.Three kinds of methods of the general employing of indication of mechanical environment: testing method, analog structure extrapolation method and mathematical model prediction method.The method of experimental test just can be carried out after need to waiting until the detector first sample stage, in the schematic design phase, is infeasible.Detector planet landing task still belongs to the first time in China's survey of deep space activity, and the mechanical environment that its landing shock section experiences does not have relevant experience and measured data accumulation, and the method for analog structure extrapolation is also infeasible.Therefore, indicate that its landing shock mechanics environment just need to adopt the method for mathematical model prediction in the detector scheme design phase.
So far, conventional lunar orbiter mathematical model is broadly divided into three classes in the world: Rigid Body Dynamics Model, Rigid-flexible Coupling Dynamics model and nonlinear finite element model.Rigid Body Dynamics Model counting yield is high, but does not consider structural flexibility, shock response that can not actual response detector useful load.Rigid-flexible Coupling Dynamics model cannot be described the elastic-plastic deformation of lunar soil, and the interaction of detector buffer gear and lunar soil depends on experimental formula.Nonlinear finite element model can be good at considering the impact of detector plastic deformation and various non-linear factors, and testing surface adopts the analysis result of nonlinear finite element model to be better than the above two.But among three kinds of models, nonlinear finite element model solution efficiency is minimum, affects designer's work efficiency, therefore how obtaining an accurate and efficient mathematical model becomes one of gordian technique of detector soft landing ballistics environment indication.
Summary of the invention
The object of the invention is must not the shock response of actual response detector useful load for solution problem, a kind of moon detector in flexible landing ballistics environment predicting method based on non linear finite element analysis is proposed, be applicable to the detector scheme design phase.
Technical scheme of the present invention specifically comprises lunar orbiter modeling method and mechanical environment analytical approach: first, set up the nonlinear finite element model of detector lunar surface soft landing, utilize broad sense dynamic condensation method to carry out model reduction to detector, reach the object that improves solution efficiency; Then reduced-order model is carried out to nonlinear finite element and solve, obtain the dynamic response of sounder structure; Obtaining being concerned about after the acceleration responsive of position, adopting the shock response spectrum of recurrence digital filtering method computation structure acceleration responsive, for describing the mechanical environment of detector landing shock; Finally adopt the method for envelope spectrum to obtain greatest hope environment, for formulating design and the test condition of equipment.
A moon detector in flexible landing ballistics environment predicting method, specifically comprises the steps:
Step 1: lunar orbiter is encouraged to prediction, the dynamic response while selecting finite element prediction method prediction detector to land according to the feature of excitation.
Step 2: the method for employing mathematical model prediction is obtained the acceleration responsive of measuring point on detector, and described measuring point is the whole finite element nodes on the plate of cabin.
Step 2.1: set up the dynamic (dynamical) nonlinear finite element model of moon detector in flexible landing.Described nonlinear finite element model comprises centrosome, fuel tank, lunar rover, buffer gear and lunar soil finite element model.
In the present lunar soil material of model non-linear body and buffer gear in the moulding constitutive relation of bullet of padded coaming and the contact relation of buffer gear and lunar soil.
Step 2.2: utilize broad sense dynamic condensation method to carry out model reduction to detector centrosome finite element model.
Step 2.3: to the detector centrosome finite element model loaded load after depression of order, carry out Nonlinear FEM Simulation, obtain the acceleration responsive of all measuring points.
Step 3: ballistics environment when detector is landed is described.The present invention uses the ballistics environment of the shock response profiling detector of acceleration.
Step 3.1: the acceleration responsive that extracts all measuring points on the plate of the same cabin of detector centrosome.
Step 3.2: the acceleration responsive of all measuring points on the plate of same cabin is converted to shock response spectrum with recurrence digital filtering algorithm.
Step 4: determine greatest hope environment (Maximum expected environment, MEE).Greatest hope environment is the upper limit of an expection of detector landing shock mechanics environment, and it can retrain the prediction frequency spectrum of all putting in region, guarantees that any point can seriously not surpass spectral range.In order to guarantee that test can produce a conservative result, the present invention adopts maximum spectrum as greatest hope environment.
Concrete definite method is: on the same cabin of the detector centrosome plate that step 3.2 is obtained, the shock response spectrum of all measuring point acceleration carries out envelope.Then the multiple typical condition running in landing according to detector changes the landing attitude of detector, and operating mode repeating step 2.3 of every change is to step 3.2, obtains under every kind of operating mode the envelope spectrum of the shock response spectrum of all measuring point acceleration on the plate of same cabin.Finally the envelope spectrum of the shock response spectrum of all measuring point acceleration on same cabin plate under above-mentioned all typical conditions is got to envelope again, the greatest hope environment of this cabin plate during as detector landing shock.
Described multiple typical condition comprises that detector lands and lands in horizontal lunar surface and have the landing leg of initial level speed, detector to land and land successively in slope lunar surface in the landing leg of slope lunar surface, detector with symmetric mode in horizontal lunar surface and without initial level speed, detector.
Step 5: formulate the test condition of detector landing shock mechanics environment, for Change In Design and the optimization of useful load on detector provides decision information.The testing standard of the mechanical environment of timing considerations same frequency segment processed, other stage of existing detector (carrier rocket transmitter section, month transfer leg, ring month section).Method is that the mechanical environment condition in other stage is converted to corresponding shock response spectrum, sees the greatest hope environment that whether can cover landing impact section.If can cover, do not need to formulate extra test condition.If can not cover, consider improve original test condition or formulate extra test condition.
Beneficial effect
The present invention utilizes the mathematical model of lunar orbiter to indicate the mechanical environment of its landing shock, can provide the restrictive condition in design for numerous useful load such as the radar on detector, camera, antenna, vector engine, masts in starting stage of lunar orbiter design, do not rely on any test figure, saved product design costs; Before moon detector in flexible landing dynamics is analyzed, adopt the method for broad sense dynamic condensation to carry out model reduction to detector central body model, improved solution efficiency, shortened the product design cycle.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of moon detector in flexible landing ballistics environment predicting method of the present invention;
Fig. 2 is the assembly finite element model schematic diagram of lunar orbiter in embodiment;
Fig. 3 is the amplitude-versus-frequency curve of detector damper leg to the exciting force of centrosome in embodiment;
Fig. 4 is shock response spectrum concept map in embodiment;
Fig. 5 is the envelope again of all measuring point acceleration of centrosome top board shock response spectrum envelope spectrum under 5 kinds of typical conditions in embodiment;
Fig. 6 is the comparison of the acceleration shock response spectrum of sinusoidal vibration environment equivalence in embodiment and the centrosome top board shock response spectrum that finite element solving obtains;
Fig. 7 is the comparison of the acceleration shock response spectrum of random vibration environment equivalence in embodiment and the centrosome top board shock response spectrum that finite element solving obtains.
Label declaration: 1-lunar rover, 2-centrosome, 3-damper leg, 4-foot pad, 5-lunar soil.
Embodiment
The centrosome top board (panel detector structure finite element model is shown in Fig. 2) of take is example, formulates its ballistics Environment Design and test condition.
Step 1: excitation is predicted, selects the dynamic response of suitable method predict according to the feature of excitation.
For soft lunar landing detector, amplitude limit effect due to Landing Buffer device, the excitation energy acting in detector matrix main structure mainly concentrates on low-frequency range, therefore the dynamic response of detector main structure also concentrates on low-frequency range, according to NASA handbook < < NASA-HDBK-7005 > >, the prediction that can respond by Finite Element Method, i.e. detector soft landing dynamic response is the Transient Dynamics problem of a low frequency.
Step 1.1: the Rigid Body Dynamics Model of setting up the soft landing of lunar orbiter lunar surface.(most conventional rigid dynamics software and the rigid dynamics solver of Nonlinear Finite meta software can be realized, as MSC.ADAMS, Abaqus, HyperWorks etc.Setting up Rigid Body Dynamics Model is in order to improve the speed of solving.)
Step 1.2: solve and obtain the exciting force-time curve of damper leg to detector centrosome.
Step 1.3: the load-time curve in step 1.2 is carried out to Fourier transform, obtain the amplitude-versus-frequency curve (as shown in Figure 3) of exciting force.
Step 1.4: the amplitude-versus-frequency curve that analytical procedure 1.3 obtains, determine the main frequency component of exciting force, scope is [0,500] Hz.This 500Hz is using the mode cutoff frequency during as model reduction for step 2.2.
Step 2: response prediction, adopts the method for mathematical model prediction to obtain the acceleration responsive of measuring point on detector.
Step 2.1: utilize Abaqus software to set up the dynamic (dynamical) nonlinear finite element model of moon detector in flexible landing.
Step 2.2: utilize " broad sense dynamic condensation " method that in Abaqus software, Substructure module provides to carry out model reduction to detector centrosome.
Step 2.2.1: detector centrosome is carried out to free vibration model analysis, obtain the 900th rank model frequency and be about 500Hz.(most business finite element softwares can carry out model analysis, and the present invention utilizes the Linear perturbation analysis module of Abaqus for simplicity)
Step 2.2.2: the cutoff frequency of broad sense dynamic condensation method is got 500Hz, main degree of freedom is got 72 degree of freedom on 12 nodes that damper leg is connected with centrosome.In the Substructure of Abaqus software module, centrosome is carried out to dynamic condensation, then assemble with Nonlinear Free degree (snubber assembly and lunar soil).
Step 2.3: utilize Abaqus/Standard solver to carry out Nonlinear FEM Simulation to detector with 5 kinds of typical conditions of different attitude soft landings, wherein moon acceleration is got 1.63m/s 2.
Step 3: mechanical environment is described, ballistics environment is described by the frequency spectrum of certain type conventionally, China GB GB2423, IEC68-2-27 and IEC721-3, American army mark MIL-STD-810D and the mark DEF STAN07-55 of the British army require to adopt acceleration as test data.So use the ballistics environment of the shock response profiling detector of acceleration in the present invention, shock response spectrum principle as shown in Figure 4.
Step 3.1: the acceleration responsive that extracts all finite element nodes on centrosome top board.
Step 3.2: the acceleration responsive of all nodes on top board is converted to shock response spectrum with recurrence digital filtering algorithm.(amplification factor Q gets 10)
Step 4: the determining of greatest hope environment (Maximum expected environment, MEE).Greatest hope environment is the upper limit of an expection of detector landing shock mechanics environment, and it can retrain the prediction frequency spectrum of all putting in region, guarantees that any point can seriously not surpass spectral range.In order to guarantee that test can produce a conservative result, the present invention has adopted maximum spectrum as greatest hope environment.
Step 4.1: respectively under 5 kinds of typical conditions to centrosome top board on the acceleration shock response spectrum of all nodes carry out envelope, corresponding 1 envelope spectrum of every 1 operating mode.
Described 5 kinds of typical conditions are respectively:
Operating mode 1 detector lands in horizontal lunar surface, and detector is without initial level speed.
Operating mode 2 detectors land in horizontal lunar surface, and detector has the initial level speed of 1m/s.
Operating mode 3 detectors are so that " 2-2 " land in 15 ° of slope lunar surfaces, has 2 to land in lunar surface in 4 of detector foot pads, and then all the other 2 are landed in lunar surface again by symmetrical mode simultaneously simultaneously.
Operating mode 4 detectors land in 15 ° of slope lunar surfaces in " 1-2-1 " symmetrical mode, in 4 of detector foot pads, have at first 1 to land in lunar surface, and then other 2 foot pads land in lunar surface simultaneously, and remaining 1 is finally landed in lunar surface.
Operating mode 5 detectors land in 15 ° of slope lunar surfaces in the mode of " 1-2-3-4 ", and 4 of detector sufficient pads successively (during difference) are landed in lunar surface.
Step 4.2: to the envelope again of 5 envelope spectrums in step 4.1, as the greatest hope environment of centrosome top board.The envelope again of 5 kinds of envelope spectrums as shown in Figure 5.
Step 5: formulate design and test condition.During formulation, can take into full account the testing standard of other stage mechanical environment of the existing detector of same frequency segment, for example, deliver sinusoidal vibration (as shown in table 1) and the random vibration mechanical environment condition (as shown in table 2) of launching phase.Method is that the mechanical environment condition in other stage is converted to corresponding shock response spectrum, sees the greatest hope environment that whether can cover landing impact section.If can cover, do not need to formulate extra test condition.If can not cover, consider improve original test condition or formulate extra test condition.
Step 5.1: the centrosome top board shock response spectrum comparison that sinusoidal vibration environment is equivalent to acceleration shock response spectrum and obtains with step 4.2.The sinusoidal vibration environment of centrosome top board is provided by table 1.Equivalent method is as follows:
In little damping situation, when excitation frequency equals the natural frequency of the system that is energized, acceleration responsive can be determined by following formula
u &CenterDot; &CenterDot; max x &CenterDot; &CenterDot; max = Q ( 1 - e - &pi;N Q ) - - - ( 1 )
In formula,
Figure BDA0000446560100000062
the maximum absolute acceleration response of system,
Figure BDA0000446560100000063
excitation acceleration peak value, (ξ is damping ratio), N is the sinusoidal excitation cycle, when N is enough large, has
u &CenterDot; &CenterDot; max x &CenterDot; &CenterDot; max = Q - - - ( 2 )
The input amplitude (shown in table 1) of supposing the sinusoidal excitation in regulation frequency range is
Figure BDA0000446560100000075
, according to formula, the shock environment spectrum value of its equivalence is if get Q=10, the centrosome top board shock response spectrum being obtained by acceleration shock response spectrum and the step 4.2 of table 1 equivalence more as shown in Figure 6.As seen from the figure, in 0~100Hz frequency range, sinusoidal vibration environment has covered the landing shock environment of centrosome top board completely.
Step 5.2: the centrosome top board shock response spectrum comparison that random vibration environment is equivalent to acceleration shock response spectrum and obtains with step 4.2.The sinusoidal vibration environment of centrosome top board is provided by table 2.Equivalent method is as follows:
The input amplitude (as shown in table 2) of supposing the arbitrary excitation in regulation frequency range is G (f), and the shock environment spectrum value of its equivalence is
a max = 3 &pi; 2 fQG ( f ) - - - ( 3 )
In formula, f is frequency (Hz of unit), if get Q=10, the centrosome top board shock response spectrum being obtained by acceleration shock response spectrum and the step 4.2 of table 2 equivalence more as shown in Figure 7.As seen from the figure, in the frequency range more than 100Hz, random vibration environment has covered the landing shock environment of centrosome top board.
Step 5.3: because sinusoidal vibration environment in 5Hz~100Hz has covered the landing shock environment of centrosome top board completely, 100Hz has covered the landing shock environment of centrosome top board with random vibration environment in super band, so do not need to formulate extra test condition for centrosome top board.
Table 1 sine vibration test condition
Figure BDA0000446560100000073
Table 2 random vibration test condition
Figure BDA0000446560100000074

Claims (4)

1. moon detector in flexible landing ballistics environment predicting method, is characterized in that: specifically comprise the steps:
Step 1: lunar orbiter is encouraged to prediction, the dynamic response while selecting finite element prediction method prediction detector to land according to the feature of excitation;
Step 2: the method for employing mathematical model prediction is obtained the acceleration responsive of measuring point on detector, and described measuring point is the whole finite element nodes on the plate of cabin;
Step 2.1: set up the dynamic (dynamical) nonlinear finite element model of moon detector in flexible landing; Described nonlinear finite element model comprises centrosome, fuel tank, lunar rover, buffer gear and lunar soil finite element model;
Step 2.2: utilize broad sense dynamic condensation method to carry out model reduction to detector centrosome finite element model;
Step 2.3: to the detector centrosome finite element model loaded load after depression of order, carry out Nonlinear FEM Simulation, obtain the acceleration responsive of all measuring points;
Step 3: ballistics environment when detector is landed is described;
Step 3.1: the acceleration responsive that extracts all measuring points on the plate of the same cabin of detector centrosome;
Step 3.2: the acceleration responsive of all measuring points on the plate of same cabin is converted to shock response spectrum with recurrence digital filtering;
Step 4: determine greatest hope environment; Greatest hope environment is the expection upper limit of detector landing shock mechanics environment, can retrain the prediction frequency spectrum of all putting in region, guarantees that any point can seriously not surpass spectral range; Adopt maximum spectrum as greatest hope environment;
Concrete definite method is: on the same cabin of the detector centrosome plate that step 3.2 is obtained, the shock response spectrum of all measuring point acceleration carries out envelope; Then the multiple typical condition running in landing according to detector changes the landing attitude of detector, and operating mode repeating step 2.3 of every change is to step 3.2, obtains under every kind of operating mode the envelope spectrum of the shock response spectrum of all measuring point acceleration on the plate of same cabin; Finally the envelope spectrum of the shock response spectrum of all measuring point acceleration on same cabin plate under above-mentioned all typical conditions is got to envelope again, the greatest hope environment of this cabin plate during as detector landing shock;
Step 5: formulate the test condition of detector landing shock mechanics environment, for Change In Design and the optimization of useful load on detector provides decision information; Timing considerations same frequency segment processed, existing detector carrier rocket transmitter section, the testing standard of mechanical environment of month transfer leg, ring month section; The mechanical environment condition in other stage is converted to corresponding shock response spectrum, if can cover the greatest hope environment of landing impact section, does not need to formulate extra test condition; If can not cover the greatest hope environment of landing impact section, improve original test condition or formulate extra test condition.
2. moon detector in flexible landing ballistics environment predicting method according to claim 1, is characterized in that: in the present lunar soil material of model non-linear body and buffer gear in the moulding constitutive relation of bullet of padded coaming and the contact relation of buffer gear and lunar soil.
3. moon detector in flexible landing ballistics environment predicting method according to claim 1, is characterized in that: the present invention uses the ballistics environment of the shock response profiling detector of acceleration.
4. moon detector in flexible landing ballistics environment predicting method according to claim 1, is characterized in that: described multiple typical condition comprises that detector lands and lands in horizontal lunar surface and have the landing leg of initial level speed, detector to land and land successively in slope lunar surface in the landing leg of slope lunar surface, detector with symmetric mode in horizontal lunar surface and without initial level speed, detector.
CN201310727234.1A 2013-12-25 2013-12-25 Mechanical environment prediction method of lunar probe soft landing impact Expired - Fee Related CN103678822B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310727234.1A CN103678822B (en) 2013-12-25 2013-12-25 Mechanical environment prediction method of lunar probe soft landing impact

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310727234.1A CN103678822B (en) 2013-12-25 2013-12-25 Mechanical environment prediction method of lunar probe soft landing impact

Publications (2)

Publication Number Publication Date
CN103678822A true CN103678822A (en) 2014-03-26
CN103678822B CN103678822B (en) 2015-02-25

Family

ID=50316358

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310727234.1A Expired - Fee Related CN103678822B (en) 2013-12-25 2013-12-25 Mechanical environment prediction method of lunar probe soft landing impact

Country Status (1)

Country Link
CN (1) CN103678822B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104809300A (en) * 2015-05-07 2015-07-29 北京理工大学 Dynamic simulation method for pulse substructure and finite element rigid-elastic mixed connection
CN105243246A (en) * 2015-11-18 2016-01-13 北京理工大学 Landing mechanical environment analysis method based on returner component
CN106773671A (en) * 2016-11-21 2017-05-31 上海卫星工程研究所 Deep space probe MDO methods stage by stage
CN113158418A (en) * 2021-02-25 2021-07-23 北京空间飞行器总体设计部 Multi-band test simulation method for lunar surface soft landing mechanical environment
CN113656935A (en) * 2021-07-05 2021-11-16 清华大学 Asteroid rigid-flexible coupling detector and calculation method of simulation motion parameters thereof
CN114239131A (en) * 2021-11-24 2022-03-25 沈阳航天新光集团有限公司 Attitude control power system space flight mechanics environment prediction method
CN118504360A (en) * 2024-07-17 2024-08-16 烟台哈尔滨工程大学研究院 Impact load path identification method based on impact energy fitting

Citations (2)

* Cited by examiner, † Cited by third party
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 (en) * 2012-11-29 2013-03-20 北京理工大学 Control method of soft landing autonomic obstacle avoidance fixed-thruster

Patent Citations (2)

* Cited by examiner, † Cited by third party
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 (en) * 2012-11-29 2013-03-20 北京理工大学 Control method of soft landing autonomic obstacle avoidance fixed-thruster

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104809300A (en) * 2015-05-07 2015-07-29 北京理工大学 Dynamic simulation method for pulse substructure and finite element rigid-elastic mixed connection
CN104809300B (en) * 2015-05-07 2016-11-23 北京理工大学 Pulse minor structure and finite element have just played the dynamics simulation method of Hybrid connections
CN105243246A (en) * 2015-11-18 2016-01-13 北京理工大学 Landing mechanical environment analysis method based on returner component
CN105243246B (en) * 2015-11-18 2017-03-29 北京理工大学 Based on the landing mechanical environment analysis method for returning device assembly
CN106773671A (en) * 2016-11-21 2017-05-31 上海卫星工程研究所 Deep space probe MDO methods stage by stage
CN106773671B (en) * 2016-11-21 2019-08-23 上海卫星工程研究所 Deep space probe MDO method stage by stage
CN113158418A (en) * 2021-02-25 2021-07-23 北京空间飞行器总体设计部 Multi-band test simulation method for lunar surface soft landing mechanical environment
CN113656935A (en) * 2021-07-05 2021-11-16 清华大学 Asteroid rigid-flexible coupling detector and calculation method of simulation motion parameters thereof
CN113656935B (en) * 2021-07-05 2024-03-15 清华大学 Asteroid rigid-flexible coupling detector and calculation method of simulation motion parameters thereof
CN114239131A (en) * 2021-11-24 2022-03-25 沈阳航天新光集团有限公司 Attitude control power system space flight mechanics environment prediction method
CN114239131B (en) * 2021-11-24 2024-10-01 沈阳航天新光集团有限公司 Space flight mechanical environment prediction method of attitude control power system
CN118504360A (en) * 2024-07-17 2024-08-16 烟台哈尔滨工程大学研究院 Impact load path identification method based on impact energy fitting
CN118504360B (en) * 2024-07-17 2024-09-24 烟台哈尔滨工程大学研究院 Impact load path identification method based on impact energy fitting

Also Published As

Publication number Publication date
CN103678822B (en) 2015-02-25

Similar Documents

Publication Publication Date Title
CN103678822B (en) Mechanical environment prediction method of lunar probe soft landing impact
Ding et al. Influence of the North Pacific Victoria mode on the Pacific ITCZ summer precipitation
Liu et al. Hydrodynamic performance in a sloshing liquid oxygen tank under different initial liquid filling levels
Cooke et al. Assessing the work budget and efficiency of fault systems using mechanical models
CN110489904B (en) Dynamic model and parameter identification method of metal rubber vibration isolation system
CN109828033A (en) Damnification recognition method and system based on vibratory response similarity analysis
Hurlburt The potential for ocean prediction and the role of altimeter data
CN113051787A (en) Air-to-air missile hanging fatigue life estimation method and system based on short-time dynamic stress measurement
CN102730205B (en) Low-rigidity suspension system and implementation method thereof
Luo et al. Design and experimental research on buffer protection of high-g penetrator for deep space exploration
Wang et al. A nonlinear dynamic model for 2D deepwater steel lazy-wave riser subjected to top–end imposed excitations
CN102787677A (en) Engineering structure seismic intensity prediction method
Li et al. Experimental and numerical seismic investigations of the Three Gorges dam
Gücüyen Analysis of offshore wind turbine tower under environmental loads
CN105160084A (en) Ejection escape simulation method
Pham et al. Landing stability simulation of a 1/6 lunar module with aluminum honeycomb dampers
Ding et al. Numerical and experimental investigation on the shock mitigation of satellite-rocket separation
CN109446702B (en) Passive vibration reduction method of space science experiment cabinet
Shi et al. Megathrust locking and viscous mantle flow induce continental shortening in Central Andes
Wang et al. 3D geodynamic models for HP‐UHP rock exhumation in opposite‐dip double subduction‐collision systems
CN106202747B (en) A kind of scale-span analysis method of silicon-carbon base ceramic coating interface mechanics characteristic
Damgaard et al. Assessment of dynamic substructuring of a wind turbine foundation applicable for aeroelastic simulations
CN105570377B (en) Control-moment gyro micro-vibration passive vibration isolation device
CN210571263U (en) Dam body anti-seismic test simulation device with controllable vibration frequency partitions
Wan et al. Modeling and simulation of lunar lander soft-landing using transient dynamics approach

Legal Events

Date Code Title Description
PB01 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: 20191225

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