CN103606332A - Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof - Google Patents

Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof Download PDF

Info

Publication number
CN103606332A
CN103606332A CN201310547320.4A CN201310547320A CN103606332A CN 103606332 A CN103606332 A CN 103606332A CN 201310547320 A CN201310547320 A CN 201310547320A CN 103606332 A CN103606332 A CN 103606332A
Authority
CN
China
Prior art keywords
max
acceleration
guide rail
docking
vertical beam
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.)
Pending
Application number
CN201310547320.4A
Other languages
Chinese (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201310547320.4A priority Critical patent/CN103606332A/en
Publication of CN103606332A publication Critical patent/CN103606332A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Navigation (AREA)

Abstract

The invention provides a spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and a device thereof. The device comprises a machine body, a cross beam, a vertical beam, a three-shaft turntable, a motion simulation device controller, a GNC system and a measurement system. The machine body is mounted on a ground base. The cross beam is placed on the guide rail of the machine body. The vertical beam is installed on the guide rail of the cross beam. The three-shaft turntable is installed on a vertical beam guide rail. The motion simulation device controller, the GNC system and the measurement system are placed next to a system. Based on a similarity theorem and a scaling ratio criterion, the effective travel range of a simulation device is determined according to an aircraft actual flight range, and assuming that a scaling parameter is k, according to the actual flight condition of the aircraft, assuming that an aircraft true maximum speed, an acceleration, an angular speed and an angular acceleration are v[max], a[max], omega[max], and alpha[max] respectively, the maximum speed, the acceleration, the angular speed and the angular acceleration of the simulation device are determined as k*v[max], k*a[max], omega[max], and alpha[max]. The method and the device have the advantages of a wide application range, a simple structure, and multi degrees of freedom.

Description

Space Rendezvous And Docking docking multiple degrees of freedom semi-physical simulation method and device thereof
Technical field
The present invention relates to emulation technology, is exactly semi-physical simulation method and device thereof to section tracking aircraft far away and target aircraft relative motion and control method in Technique in Rendezvous and Docking specifically.
Background technology
Space Rendezvous And Docking docking technique is complicated, because environment of living in is special, need to guarantee its reliability and precision, therefore must carry out on ground emulation.Although since the sixties in last century, the numerous and confused Technique in Rendezvous and Docking Simulation system of setting up separately in various countries, at present in the research of emulation mode of section far away rare.But in the final intersection docking of the result direct relation result of section far away, therefore in the research of emulation mode of section far away there is important practical significance and using value.
Through searching document, find, Chinese invention patent application number: 200910243276.1, patent name Manual Control Rendezvous and Docking semi-physical simulation system, this invention adopts the motion of Three-degree of Freedom Rotational Platform simulated target aircraft, and the present invention adopts the motion of a Three-degree of Freedom Rotational Platform and a common simulated target aircraft of Three Degree Of Freedom translational motion simulator.Two aircraft have ten two degrees of freedom in real space, and former invention is also adopting the design of nine-degree of freedom to carry out analog simulation, so former invention is not high to the simulation degree of Technique in Rendezvous and Docking, precision is limited.And this invention adopts the applicable measuring distance minimum of laser rendezvous radar at 50 meters, so application scenario is very limited.The present invention adopts the sensor of vision guided navigation camera and laser range finder to substitute laser rendezvous radar, has improved precision and the range of application of measuring emulation.
Chinese invention patent application number: 200910243277.6, patent name is a kind of human control intersection docking operation method, this patent has mainly designed in a kind of intersection docking test, experimenter's method of operating, mainly to use for personnel training, be not emulation test method involved in the present invention, the two is not same class invention.
Summary of the invention
The object of the present invention is to provide a kind of multiple degrees of freedom Space Rendezvous And Docking docking semi-physical simulation method and device thereof.
Device of the present invention is achieved in that a kind of Space Rendezvous And Docking docking multiple degrees of freedom semi-physical simulation device, comprise lathe bed, crossbeam, vertical beam, three-axle table, motion simulator controller, GNC system and measuring system, lathe bed is installed on ground, crossbeam is placed on the guide rail of lathe bed, vertical beam is arranged on the guide rail of crossbeam, three-axle table is installed on vertical beam guide rail, and motion simulator controller, GNC system and measuring system are placed in system side.
Method of the present invention is: based on correspondence theorem and contracting, compare criterion, according to aerocraft real flight range, determine the effective travel of simulator, suppose that contracting is k than parameter, according to aerocraft real flight situation, suppose that the true maximal rate of aircraft, acceleration, angular velocity and angular acceleration are respectively v max, a max, ω max, α max, determine that maximal rate, acceleration, angular velocity and the angular acceleration of simulator is kv max, ka max, ω max, α max.
The invention has the beneficial effects as follows: applied range, simple in structure, multiple degrees of freedom.
Accompanying drawing explanation
Fig. 1 is the key diagram to respective angles in coordinate system;
Fig. 2 is ground experiment system principle diagram;
Fig. 3 is for adopting cross beam frame in the schematic diagram of earth construction.
Embodiment
Below in conjunction with Figure of description, the invention will be further described:
Embodiment 1
In conjunction with Fig. 1, this example is set up the process of kinetic model mainly for illustrative system, after system modelling, and stressed in the time of can clear and definite flight simulator motion, and people is for applying this power, and the truth to aircraft in space is simulated.In addition, in the computation process of the topworks's index to analogue system, also need the kinetic model of the system that is applied to.
According to spatial dynamics principle, rendezvous and docking system is carried out to modeling, obtains Technique in Rendezvous and Docking kinetic model:
x . . - 2 ω y . = a x y . . + 2 ω x . - 3 ω 2 y = a y z . . + ω 2 z = a z
Intersection Docking simulation model:
x . . - 2 ω y . = a x k y . . + 2 ω x . - 3 ω 2 y = a y k z . . + ω 2 z = a z k
Wherein k is emulation contracting ratio.
And intersection docking measurement model:
ρ = x 4 2 + y 4 2 + z 4 2 α = arctan z 4 x 4 β = arctan y 4 x 4 2 + y 4 2
ρ: passive space vehicle for the distance of pursuit spacecraft (can be understood as in Fig. 1 | O 1o 2|);
α: passive space vehicle is the angle of pitch in the coordinate system of spaceborne measuring system for pursuit spacecraft.Its angle of pitch be defined as passive space vehicle for pursuit spacecraft the projection in pursuit spacecraft body coordinate system X4Z4 plane and the angle between X4 axle positive axis, X4 axle positive dirction is just (can with reference to figure 1);
β: passive space vehicle is the position angle in the coordinate system of spaceborne measuring system for pursuit spacecraft.Its position angle be defined as passive space vehicle for pursuit spacecraft line the angle between pursuit spacecraft body coordinate system X4Z4 plane projection, face upward as just.
Embodiment 2
Definite method of the main explanation of this example emulation contracting ratio.Emulation contracting represents with k than in this article, refers to that laboratory simulations environmental field is to Technique in Rendezvous and Docking Simulation scope a ratio, and generally it is one and is less than 1 positive number.
First according to simulation requirements, determine the selection scheme of sensor; Meanwhile, according to sensor accuracy, want in the derivation of equation and Technique in Rendezvous and Docking section real data far away and in conjunction with self laboratory condition, determine that suitable emulation contracting compares.
For example: in space, section tracking aircraft far away and the true relative distance of target aircraft are probably 10km.The precision satisfied with analogue system compared in emulation contracting, should meet following relation:
ρ 2 ≥ kρ 1 Δ ρ 2 ≤ kΔ ρ 1
Wherein, Δ ρ 1, Δ ρ 2be respectively the measuring error of sensor replacement scheme and laser rendezvous radar, ρ 1, ρ 2be respectively the measurement maximum magnitude of sensor replacement scheme and laser rendezvous radar.Table 4-1 has provided the measurement of correlation parameter of the laser rendezvous radar of aerocraft real employing, according to the precision of emulation demand (ground short distance cannot be used laser rendezvous radar) and actual available sensors, determine that the sensor plan by vision camera and laser range finder combination substitutes the laser rendezvous radar sensor plan on star.Table 4-2 has provided the laser range finder of analogue system plan employing and the measurement of correlation parameter of vision sensor.
Table 1 laser rendezvous radar measurement of correlation parameter list
Figure BSA0000097364610000041
Table 2 laser range finder and vision sensor measurement of correlation parameter list
Figure BSA0000097364610000042
With reference to the vision camera providing above, laser range finder, the precision parameter of laser rendezvous radar, and the formula proposing, can determine that contracting is 0.003≤k≤0.02 than k scope.Now, can also be according to actual conditions, and formula
Figure BSA0000097364610000043
determine the three-dimensional scope of simulator.For example according to our actual conditions in laboratory, get k=0.0043, obtain crossbeam length 39m, vertical beam length 5m, lathe bed length 17m, can simulate in the middle of real space x = 9304 m y = 3489 m z = 1163 m Motion in scope.
Embodiment 3
This routine article is for the computing method of analogue system topworks index.
Introduced the kinetic model of system above, it meets Hill's equation.Separate Hill's equation, can obtain the movement position of system and the ideal solution of speed.
Separate position:
x = x 0 + 6 ( θ - sin θ ) y 0 + ( 4 sin θ - 3 θ ) x 0 . n + 2 ( 1 - cos θ ) y 0 . n + [ 4 ( 1 - cos θ ) - 1.5 θ 2 ] a x n 2 + 2 ( θ - sin θ ) a y n 2 y = ( 4 - 3 cos θ ) y 0 - 2 ( 1 - cos θ ) x 0 . n + sin θ n y 0 + 2 ( sin θ - θ ) a x n 2 + ( 1 - cos θ ) a y n 2 z = z 0 cos θ + sin θ n z 0 . + ( 1 - cos θ ) a z n 2
Velocity solution:
x . = 6 n ( 1 - cos θ ) y 0 + ( 4 cos θ - 3 ) x 0 + 2 y 0 sin θ + ( 4 sin θ - 3 θ ) a x n + 2 ( 1 - cos θ ) a y n y . = 3 n y 0 sin θ - 2 x 0 . sin θ + y 0 . cos θ + 2 ( cos θ - 1 ) a x n + a y n sin θ z . = - nz 0 sin θ + z . 0 cos θ + a z n sin θ
The angle that θ in formula---pursuit spacecraft turns over around passive space vehicle (rad) is π;
The rotational angular velocity of n---pursuit spacecraft (rad/s).
The moment cartesian coordinate system upper/lower positions for t is separated in position, and velocity solution is corresponding speed under t moment cartesian coordinate system.Can be according to kepler's third law for rotational angular velocity n:
2 π T = μ a 3
A in formula---shift elliptical orbit center to apogean distance, i.e. a=(r 1+ r 0)/2;
R 1---the earth's core is to the distance of passive space vehicle track, 6768km;
R 0---the earth's core is to the distance of the motor-driven front track of pursuit spacecraft, 6778km;
T---the orbital period;
μ---Gravitational coefficient of the Earth, numerical value is 3.986 * 105km3/s 2.
Initial error x while simultaneously considering Technique in Rendezvous and Docking, carries out error y and measuring error z, adopts the geometric mean counting method in statistical mathematics can obtain total error,
Figure BSA0000097364610000061
can then to position solution and velocity solution, carry out mathematical computations respectively by maximum positive error and the movement position of negative error substitution system and the ideal solution of speed, can calculate the extreme value of system motion parameter.By the contrast of the system motion parameter extreme value in two kinds of situations, get the motion index that maximal value can obtain system actuator.
Such as, for certain space flight intersection docking mission, known initial error x and execution error y are as follows: X axis position initial error is x x=4.3mm, Y direction position initial error is x y=8.6mm, X-direction speed initial error
Figure BSA0000097364610000062
y direction speed initial error
Figure BSA0000097364610000063
x axis executing location error is y x=8.6mm, velocity deviation is
Figure BSA0000097364610000064
y-axis position deviation is y y=4.3mm, velocity deviation is
Figure BSA0000097364610000065
and the survey sensor that measuring error z can select according to oneself is determined, the survey sensor positional precision that this example is selected is z=21.5mm, and velocity accuracy is
Figure BSA0000097364610000066
adopt the geometric mean counting method in statistical mathematics can obtain total error,
Figure BSA0000097364610000067
respectively by maximum positive error (+σ=13.59mm,
Figure BSA00000973646100000615
) and negative error (σ=-13.59mm,
Figure BSA0000097364610000069
) in the movement position of substitution system and the ideal solution of speed.
Known starting condition and final state condition:
y = y 1 = 0 y . = y 1 . = 0
Can obtain the starting condition with error:
y = y ± σ = ± 13.59 mm y . = y . ± σ . = ± 2.13 mm / s
By MATLAB, programme respectively in the ideal case again ( y = y 1 = 0 y · = y · 1 = 0 ) With with under error condition ( y = y ± σ = ± 13.59 mm y · = y · ± σ · = ± 2.13 mm / s ) Position solution and velocity solution are carried out to mathematical computations, can calculate the extreme value of system motion parameter.By the contrast of the system motion parameter extreme value in two kinds of situations, get the motion index that maximal value can obtain system actuator.Specific targets are as following table:
Table 3
Figure DEST_PATH_GSB0000120340110000071
Note: speed and unit of error: mm/s; Acceleration and unit of error: mm/s 2.
Embodiment 4
In conjunction with Fig. 2-Fig. 3, the frame mode that analogue system adopts is set up in the main explanation of this example, it adopts console mode driving scheme, by lathe bed (1), crossbeam (2), vertical beam (3), three-axle table (4), motion simulator controller (5), GNC system (6), measuring system (7) forms, it is characterized in that: lathe bed (1) is installed on ground, crossbeam (2) is placed on the guide rail of lathe bed (1), vertical beam (3) is arranged on the guide rail of crossbeam (2), three-axle table (4) is installed on vertical beam (3) guide rail, motion simulator controller (5), GNC system (6) is placed in system side with measuring system (7).
The sensor such as radar, photoelectric follow-up is measured the motion state of spacecraft simulation mechanism in real time, and is sent to GNC control analogue system.GNC control analogue system according to the current motion state of the dynamics of aircraft, kinematics and spacecraft simulation mechanism according to the navigation algorithm of implanting, send movement instruction to motion simulator controller, motion simulator controller is controlled output by execution units such as its drive motor according to the current pose of analogue means and expection pose through calculating, and drives the motion of multiple degrees of freedom flight simulator.Ground observing and controlling system by system significant data show, storage etc., thereby facilitate designer can be in real time, the motion state data of recovering and analysis process of the test.

Claims (2)

1. a Space Rendezvous And Docking docks multiple degrees of freedom semi-physical simulation device, comprise lathe bed (1), crossbeam (2), vertical beam (3), three-axle table (4), motion simulator controller (5), GNC system (6) and measuring system (7), it is characterized in that: lathe bed (1) is installed on ground, crossbeam (2) is placed on the guide rail of lathe bed (1), vertical beam (3) is arranged on the guide rail of crossbeam (2), three-axle table (4) is installed on vertical beam (3) guide rail, motion simulator controller (5), GNC system (6) is placed in system side with measuring system (7).
2. a Space Rendezvous And Docking docks multiple degrees of freedom semi-physical simulation method, it is characterized in that: based on correspondence theorem and contracting, compare criterion, according to aerocraft real flight range, determine the effective travel of simulator, suppose that contracting is k than parameter, according to aerocraft real flight situation, suppose that the true maximal rate of aircraft, acceleration, angular velocity and angular acceleration are respectively v max, a msx, ω max, α max, determine that maximal rate, acceleration, angular velocity and the angular acceleration of simulator is kv max, ka max, ω max, α max.
CN201310547320.4A 2013-10-30 2013-10-30 Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof Pending CN103606332A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310547320.4A CN103606332A (en) 2013-10-30 2013-10-30 Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310547320.4A CN103606332A (en) 2013-10-30 2013-10-30 Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof

Publications (1)

Publication Number Publication Date
CN103606332A true CN103606332A (en) 2014-02-26

Family

ID=50124551

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310547320.4A Pending CN103606332A (en) 2013-10-30 2013-10-30 Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof

Country Status (1)

Country Link
CN (1) CN103606332A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104077490A (en) * 2014-07-03 2014-10-01 哈尔滨工业大学 Aircraft navigation guidance and control ground simulation system performance evaluating method
CN110826189A (en) * 2019-10-14 2020-02-21 中国科学院力学研究所 Method for determining aircraft scale model experiment system
CN112857844A (en) * 2021-01-06 2021-05-28 安徽省湘凡科技有限公司 Triaxial simulation test rotary table

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003306200A (en) * 2002-04-15 2003-10-28 Natl Space Development Agency Of Japan Image navigation for rendezvous docking and navigation device
CN101794527A (en) * 2009-12-30 2010-08-04 北京控制工程研究所 Manual control rendezvous and docking semi-physical simulation testing system
CN203134243U (en) * 2013-01-15 2013-08-14 北京化工大学 Three degree-of-freedom motion control system teaching practical training robot platform
CN103268070A (en) * 2013-04-24 2013-08-28 哈尔滨工业大学 Space multi-motion-body relative motion scaling semi-physical simulation system
CN103323823A (en) * 2013-05-30 2013-09-25 北京控制工程研究所 Method for analyzing navigation error of rendezvous radar in rendezvous and docking

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003306200A (en) * 2002-04-15 2003-10-28 Natl Space Development Agency Of Japan Image navigation for rendezvous docking and navigation device
CN101794527A (en) * 2009-12-30 2010-08-04 北京控制工程研究所 Manual control rendezvous and docking semi-physical simulation testing system
CN203134243U (en) * 2013-01-15 2013-08-14 北京化工大学 Three degree-of-freedom motion control system teaching practical training robot platform
CN103268070A (en) * 2013-04-24 2013-08-28 哈尔滨工业大学 Space multi-motion-body relative motion scaling semi-physical simulation system
CN103323823A (en) * 2013-05-30 2013-09-25 北京控制工程研究所 Method for analyzing navigation error of rendezvous radar in rendezvous and docking

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
梁刚刚: "《空间交会对接中远段相对导航仿真验证研究》", 《中国优秀硕士学位论文全文数据库》, 15 May 2011 (2011-05-15) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104077490A (en) * 2014-07-03 2014-10-01 哈尔滨工业大学 Aircraft navigation guidance and control ground simulation system performance evaluating method
CN110826189A (en) * 2019-10-14 2020-02-21 中国科学院力学研究所 Method for determining aircraft scale model experiment system
CN112857844A (en) * 2021-01-06 2021-05-28 安徽省湘凡科技有限公司 Triaxial simulation test rotary table
CN112857844B (en) * 2021-01-06 2023-08-01 安徽省湘凡科技有限公司 Triaxial simulation test turntable

Similar Documents

Publication Publication Date Title
CN100495261C (en) Half-physical emulation test system for controlling and guiding, navigating and controlling soft landing for moon
CN103901907B (en) Soft landing obstacle avoidance simulation test system
CN102519441B (en) Method for measuring positioning points based on laser tracker in docking process of airplane parts
CN105974822B (en) A kind of spacecraft, which is independently diversion, intersects the verification method of control system ground validation device
CN102393200B (en) General inertial navigation test method based on flight simulation
CN109592079A (en) A kind of spacecraft coplanar encounter of limiting time becomes rail strategy and determines method
CN104077490A (en) Aircraft navigation guidance and control ground simulation system performance evaluating method
CN101794527B (en) Manual control rendezvous and docking semi-physical simulation testing system
CN101499220B (en) Method and apparatus for simulating large thruster on spacecraft
CN103674034B (en) Multi-beam test the speed range finding revise robust navigation method
CN104718508A (en) Three-dimensional manipulation of teams of quadrotors
CN104236546A (en) Satellite starlight refraction navigation error determination and compensation method
CN104298128A (en) Ground simulation method for spacecraft navigation guidance technology
CN106094565A (en) A kind of spacecraft autonomous rendezvous control system ground simulation test method
CN103496449A (en) Pose adjustment track planning method for plane side wall component assembling
CN103868648A (en) Barycenter measuring method for three-axis air floatation simulation experiment platform
CN112650076B (en) Constellation cooperative control ground simulation system
CN114936471B (en) Spacecraft collision early warning layered rapid screening method based on parallel computing
CN114625027A (en) Multi-spacecraft attitude and orbit control ground full-physical simulation system based on multi-degree-of-freedom motion simulator
CN103606332A (en) Spacecraft rendezvous and docking multi-degree-of-freedom semi-physical simulation method and device thereof
CN109189079A (en) Mobile Robotics Navigation control method based on GPS positioning
CN203966431U (en) The target detection of a kind of dexterous ammunition device and control analogue system
CN105737848B (en) System-level star sensor star viewing system and star viewing method
CN105973237B (en) Emulation dynamic trajectory based on practical flight data interpolating parses generation method
CN111637902A (en) Ground demonstration verification system and method for remote approach of small deep space celestial body

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20140226

RJ01 Rejection of invention patent application after publication