CN110758775A - Multi-pulse area hovering method based on asteroid surface observation - Google Patents
Multi-pulse area hovering method based on asteroid surface observation Download PDFInfo
- Publication number
- CN110758775A CN110758775A CN201911149845.6A CN201911149845A CN110758775A CN 110758775 A CN110758775 A CN 110758775A CN 201911149845 A CN201911149845 A CN 201911149845A CN 110758775 A CN110758775 A CN 110758775A
- Authority
- CN
- China
- Prior art keywords
- hovering
- detector
- asteroid
- area
- pulse
- 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
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000012937 correction Methods 0.000 claims abstract description 30
- 230000008685 targeting Effects 0.000 claims abstract description 14
- 239000000446 fuel Substances 0.000 claims abstract description 8
- 230000001788 irregular Effects 0.000 claims abstract description 7
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 230000004048 modification Effects 0.000 claims description 5
- 238000012986 modification Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 4
- 206010048669 Terminal state Diseases 0.000 claims description 3
- 239000003380 propellant Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 claims description 2
- 238000005381 potential energy Methods 0.000 claims description 2
- 239000000725 suspension Substances 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 13
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011217 control strategy Methods 0.000 description 2
- 241001233914 Chelidonium majus Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/247—Advanced control concepts for autonomous, robotic spacecraft, e.g. by using artificial intelligence, neural networks or autonomous agents
Abstract
The invention relates to a multi-pulse area hovering method based on asteroid surface observation, and belongs to the technical field of aerospace. The invention discloses a multi-pulse area hovering method based on asteroid surface observation, which considers the asteroid irregular attraction and the sunlight pressure perturbation and designs a hovering orbit meeting the constraint of the asteroid surface hovering area based on a targeting method; the fuel consumption optimization is used as a performance index, and the obtained hovering track has the minimum total speed increment under the constraint of specific correction duration; the method has the advantages of small required speed increment, less pulse applying times, wide application range and capability of meeting hovering area constraints of different sizes, and is suitable for hovering orbit design under the constraint condition of the hovering area on the surface of the asteroid.
Description
Technical Field
The invention relates to a multi-pulse region hovering method based on asteroid surface observation, in particular to a multi-pulse region hovering method considering initial state errors and asteroid irregular gravity influence, belongs to the technical field of aerospace, and is suitable for multi-pulse region hovering orbit design of a detector under the constraint of a certain range on the asteroid surface.
Background
Asteroid surface hovering detection is a very important stage for observing the asteroid surface, and hovering detection means that a detector detects on an orbit, the relative position of which with a target asteroid is kept unchanged or the detector moves only in a very small range. Considering that a strict fixed-point observation position is not required for an actual observation task, the regional hovering application prospect of the asteroid surface is wide and easy to realize. In the hovering detection stage, the detector utilizes the infrared camera to determine the surface temperature of the asteroid, the detector is ensured to safely operate in a close range of the asteroid, the characteristic points on the surface of the asteroid are observed in detail, and the hovering orbit design is crucial in order to ensure the safety and the observation effect of the detector in the hovering detection task execution process.
The design of the regional hovering orbit in the prior art [1 ] (see Zhou Hai Pan. regional relative hovering orbit design and control [ J ]. Nanyang institute of technology academy of technology, 2015) proposes an orbit optimization design method for realizing relative hovering by using single-pulse control, which can realize the design of a single-pulse hovering orbit from a boundary to a boundary in a specified hovering area, but has a maximum flight time limit during relative motion between any two points in the area, and a detector moves without a reference point in the hovering area.
In the prior art [ 2 ] (see the president of paniculate swallowwort root, the same family of information, the korean tide, and the like. J ] relative motion control simulation research on continuous low-thrust hovering configuration [ J ]. computer simulation, 2019), a control strategy for realizing hovering by using multiple continuous constant low-thrust is provided, a constant-thrust multi-arc-segment control strategy meeting regional hovering constraint is designed according to a Gaussian perturbation equation of an orbit element based on a relative motion model of the relative orbit element, the method can effectively realize regional hovering of a mission spacecraft relative to a target spacecraft, but the control time of the method in the hovering time accounts for more than half of the total time, the control thrust acceleration is too small, and the method is difficult to realize.
Disclosure of Invention
The invention aims to provide a multi-pulse region hovering method based on asteroid surface observation, which is based on a pulse control method, realizes the region hovering orbit design of a detector under the constraint of a certain range of the asteroid surface, and has the advantages of small required speed increment, less control times and easy realization of the method.
The purpose of the invention is realized by the following technical scheme:
the invention discloses a multi-pulse area hovering method based on asteroid surface observation, which comprises the following steps of:
the method comprises the following steps: and establishing a minor planet body coordinate system, and establishing a kinetic equation by considering the influence of the gravity of the minor planet and the sunlight pressure.
The centroid of the target asteroid is taken as the origin of coordinates, the equatorial plane of the asteroid is taken as an XY plane, the X axis points to the specific direction of the surface of the asteroid, the Z axis is in right-hand relation with the X axis and the Y axis, and the Z axis points upwards along the rotation axis of the asteroid.
The position vector of the detector under the asteroid system is rB=[x,y,z]TVelocity vector is vB=[vx,vy,vz]TThen, the motion equation of the detector under the asteroid system is as follows:
wherein x, y and z are three-axis position components of the detector under the minor planet body coordinate system respectively, and vx,vy,vzThe three-axis velocity components of the detector under the minor planet body coordinate system are respectively, omega is the minor planet self-rotation angular velocity vector, aastIs the minor planet gravitational acceleration vector of that position, aSRPThe acceleration vector is perturbed by the sunlight pressure on the detector.
And establishing a gravitational field of the small planet by adopting a polyhedral model. Then there isExpressing the equation of motion as a scalar can be derived:
wherein U is asteroid potential energy, aSRPx、aSRPyAnd aSRPzThe three-axis components of the sunlight pressure perturbation force in the minor planet body coordinate system are respectively.
Step two: selecting a center point of the hovering area according to the task requirement, wherein the center point is set as rtarDetermining the radius of the hovering area and the minimum hovering duration T with the radius of the hovering area being rhof。
In the hovering stage of the area based on the asteroid surface observation, the detector needs to investigate and screen target points of subsequent tasks such as landing sampling and the like, and the thermal infrared camera needs to be used for determining the asteroid surface temperature to ensure the operation safety of the detector in the asteroid close range, so the center point of the hovering area is generally selected to be above the most concerned position point of the asteroid surface in the task and is set as rtarThe size of the hovering area is determined by comprehensively considering factors such as the geographic environment near a target point, the height of the detector from the asteroid surface, the size of the angle of view of a camera carried by the detector and the like according to task requirements, the hovering area is generally selected to be ten meters to hundred meters, and if the radius of the hovering area of the asteroid surface is determined to be rho, the movable range of the detector is determined to be centered at rtarA sphere of radius p, which can be expressed as
|rB-rtar|≤ρ (4)
Hovering duration T required by detector at taskfThe internal demand time meets the space position requirement determined by the above formula.
Step three: giving the initial speed of the detector, and integrating the initial speed to the boundary of the hovering area in the minor planet body coordinate system kinetic equation from the center position of the hovering area to obtain the initial correction time t11And a corresponding pre-correction state X (t)11)。
Given the initial time t of the detector0Velocity v at target location point 00=[vx0,vy0,vz0]Then the initial state of the detector is X0=[rtar,v0]', the dynamic equation in the minor planet body coordinate system based on the initial stateThe distance rho between the position of the detector and the central point of the target hovering area is calculated in real time by integrationt. Influenced by the relative speed of the initial hovering point and the irregular gravitational field of the asteroid, the detector gradually deviates from the central point of the hovering area, possibly reaches the boundary of the hovering area along with the time, and is set at the moment t11The detector drifts to the hover region boundary for the first time, at which time ptRho, and the state before detector correction is denoted as X (t)11)。
Step four: after the boundary of the hovering area is reached for the first time, the correction transfer time is given, pulse control is applied by adopting a targeting method, the minimum speed increment is taken as a performance index, the boundary correction pulse with the optimal fuel consumption is determined, and the detector reaches the central position of the hovering area.
The detector state at the boundary is X (t)11) Giving a corrected transfer time tmAnd correcting the pulse maneuver Deltav, the corrected terminal state r of the detector can be obtainedf1。
The pulse correction orbit is designed by adopting a targeting method, and a targeting equation of the regional hovering problem based on asteroid surface observation is established as
F=rf1-rtar(5)
At this time, the correction problem is converted into a solution problem of a nonlinear equation system, and the detector is expected to be at the moment t1=t11+tmReacquires the hover region center position rtar. Meanwhile, considering the limitation of propellant carried by the detector, the optimal fuel consumption is selected as a performance index, namely, the pulse motor with the minimum speed increment is selected as a control correction amount under the condition of sufficient target hitting times, and the pulse motor is recorded as delta v1。
Step five: and repeating the third step and the fourth step, and determining the time from the center point to the boundary of the hovering area and the pulse maneuver.
The first pulse modification makes the detector at t1=t11+tmThe moment of time reaches the hovering area center point r from the hovering area boundarytarWhen the detector state is marked as X (t)1)=[rtar,v1]. At this time | v1|≠0,There is still a relative velocity of the detector and the asteroid and therefore will continue to deviate from the desired hover point, and when the hover region boundary is reached again, the next correction pulse is applied using the targeting method.
Recording the ith slave hovering area center position rtarThe time length for reaching the boundary of the hovering area is ti1I-th arrival at rtarAt a time tiI-th time from hover region boundary back to rtarThe desired pulse maneuver is Δ vi。
Step six: and determining the required number of pulse maneuvers to be applied and the corresponding optimal total pulse speed increment size according to the minimum required hovering duration of the task.
First i times when the detector reaches the center point r of the hovering areatarIs ti
The task requires a hover time period of TfWhen t isi-1<TfAnd t isi>TfWhen the region hovers to meet the task requirement, the total speed increment required by the region hovering is
And obtaining the region hovering orbit based on the asteroid surface observation.
The invention discloses a multi-pulse area hovering method based on asteroid surface observation, which considers the asteroid irregular attraction and the sunlight pressure perturbation and designs a hovering orbit meeting the constraint of the asteroid surface hovering area based on a targeting method. And the fuel consumption optimization is taken as a performance index, and the obtained hovering track has the minimum total speed increment under the constraint of specific correction duration. The method has the advantages of small required speed increment, less pulse applying times, wide application range and capability of meeting hovering area constraints of different sizes, and is suitable for hovering orbit design under the constraint condition of the hovering area on the surface of the asteroid.
Has the advantages that:
1. according to the multi-pulse area hovering method based on asteroid surface observation disclosed by the invention, because only whether the position of the detector reaches the boundary of the hovering area is judged, and only the pulse speed increment is applied to the boundary to correct the orbit, the required engine switching times are less, and the hovering area constraints with different radiuses can be met.
2. The invention discloses a multi-pulse area hovering method based on asteroid surface observation, which is characterized in that a pulse correction orbit with the minimum speed increment is obtained by using a targeting method when the hovering area boundary is reached every time, so that the total speed increment required by the method is small.
3. The invention discloses a multi-pulse area hovering method based on asteroid surface observation, which considers asteroid irregular shape gravitational perturbation and sunlight pressure perturbation and uses a targeting method to enable a detector to reach the center point of a hovering area for many times, so that the method has wide application range.
Drawings
FIG. 1 is a flow chart of a multi-pulse region hovering method based on asteroid surface observation according to the present invention;
FIG. 2 is a schematic design diagram of a multi-pulse region hovering method based on asteroid surface observation according to the present invention;
FIG. 3 is a three-dimensional diagram of a multi-pulse region hovering orbit of the present invention based on asteroid surface observation;
FIG. 4 is an XY plane diagram of a five-region hovering orbit of the multi-pulse region hovering method based on asteroid surface observation according to the present invention;
FIG. 5 is a graph showing the distance variation from the center position of a detector in the hovering process in the step five of the multi-pulse area hovering method based on the asteroid surface observation in the invention.
Detailed Description
For a better understanding of the objects and advantages of the present invention, reference should be made to the following detailed description taken in conjunction with the accompanying drawings and examples.
Example 1:
as shown in fig. 1, the multi-pulse area hovering method based on asteroid surface observation disclosed by the present invention, taking asteroid Nereus as an example, gives the following examples including the following steps:
the method comprises the following steps: and establishing a minor planet body coordinate system, and establishing a kinetic equation by considering the influence of the gravity of the minor planet and the sunlight pressure.
The mass center of the asteroid Nereus is taken as the origin of coordinates, the equatorial plane of the Nereus is taken as an XY plane, the X axis points to the position point on the equator which is farthest away from the mass center of the asteroid, and the Z axis, the X axis and the Y axis are in a right-hand relationship and point upwards along the rotation axis of the Nereus.
Let the position vector of the detector under the minor planet Nereus system be rB=[x,y,z]TVelocity vector is vB=[vx,vy,vz]TAnd the motion equation of the detector under the asteroid Nereus system is as follows:
where ω is the asteroid Nereus rotation angular velocity vector, aastIs the minor planet gravitational acceleration vector of that position, aSRPThe acceleration is perturbed by the sunlight pressure on the detector.
And establishing a gravitational field of the minor planet Nereus by adopting a polyhedral model. Then there isExpressing the equation of motion as a scalar can be derived:
wherein, aSRPx、aSRPyAnd aSRPzThe three-axis components of the sunlight pressure perturbation force in the minor planet Nereus body coordinate system are respectively.
Selecting a pendant according to task requirementsA center point of the parking area, the center point is rtarDetermining the radius of the hovering area and the minimum hovering duration T with the radius of the hovering area being rhof。
In the hovering stage of the asteroid Nereus surface area, the detector needs to investigate and screen target points of subsequent tasks such as landing sampling and the like, and the surface temperature of the asteroid needs to be determined by utilizing a thermal infrared camera to ensure the operation safety of the detector in the short-distance range of the asteroid, wherein the central point of the hovering area is selected to be above the landing point of the asteroid surface and is selected to be rtar=[1.5,1.8,1.2]km, radius of hovering area ρ 200m, the movable range of the detector is determined as being centered at rtarA sphere of radius p, which can be expressed as
|rB-rtar|≤ρ (4)
Hovering duration T required by detector at taskfThe space position requirement determined by the above formula is satisfied at the moment of need within 2 h.
Step three: giving the initial speed of the detector, and integrating the initial speed to the boundary of the hovering area in the minor planet body coordinate system kinetic equation from the center position of the hovering area to obtain the initial correction time t11And a corresponding pre-correction state X (t)11)。
Given the initial time t of the detector0Speed deviation v at target position point of 00=[0.05,0.02,-0.03]m/s, then the initial state of the detector is X0=[rtar,v0]' integrating the initial state in a minor planet body coordinate system kinetic equation, and calculating the distance rho between the position of the detector and the central point of the target hovering area in real timet. Influenced by the state error of the initial hovering point and the irregular gravitational field of the asteroid, the detector gradually deviates from the central point of the hovering area and possibly reaches the boundary of the hovering area along with the time, and at the moment t111201.5s the first time the detector drifts to the hover region boundary, at which time ρtρ 200m, and the detector state X (t)11)=[1.67669,1.88623,1.16332,2.57e-04,1.01e-04,-3.1e-05]。
Step four: after the boundary of the hovering area is reached for the first time, the correction transfer time is given, pulse control is applied by adopting a targeting method, the minimum speed increment is taken as a performance index, the boundary correction pulse with the optimal fuel consumption is determined, and the detector reaches the central position of the hovering area.
The detector state at the boundary is X (t)11) Giving a corrected transfer time tm1000s and a corrected pulse maneuver initial value Δ v ═ 1,1]m/s, the corrected terminal state r of the detector can be obtainedf1。
The pulse correction orbit is designed by adopting a target practice method, and a target practice equation for the problem of hovering of the surface area of the asteroid Nereus is established as
F=rf1-rtar(5)
At this time, the correction problem is converted into a solution problem of a nonlinear equation system, and the detector is expected to be at the moment t1=t11+tm2201.5s again reaches the hover region center position rtar=[1.5,1.8,1.2]And km. Meanwhile, considering the limitation of propellant carried by the detector, the optimal fuel consumption is selected as a performance index, namely, the pulse motor with the minimum speed increment is selected as a control correction quantity under the condition of sufficient target hitting times, and the first minimum speed pulse is obtained as delta v1=[-0.510,-0.024,0.019]m/s。
Step five: and repeating the third step and the fourth step, and determining the time from the center point to the boundary of the hovering area and the pulse maneuver.
The first pulse modification makes the detector at t1=t11+tmTime instant 2201.5s from the hovering area boundary to the hovering area center point rtar=[1.5,1.8,1.2]km, in which the detector state is denoted X (t)1)=[rtar,v1]. At this time | v1If | not equal to 0, the detector and the minor planet Nereus still have relative speeds and therefore will continue to deviate from the desired hover point, and when the hover region boundary is reached again, the next correction pulse is applied using the targeting method.
Recording the ith slave hovering area center position rtarThe time length for reaching the boundary of the hovering area is ti1I-th arrival at rtarAt a time tiI th time from the suspensionStopping zone boundary back to rtarThe optimum velocity pulse required is Δ vi。
Step six: and determining the required number of pulse maneuvers to be applied and the corresponding optimal total pulse speed increment size according to the minimum required hovering duration of the task.
The first 4 times that the detector reaches the center point r of the hovering areatar=[1.5,1.8,1.2]The total duration of km is t3
The task requires a hover time period of Tf=2h,t3<TfAnd t is4>TfWhen the region hovers to meet the task requirement, the total speed increment required by the region hovering is
The parameters for each pulse modification are shown in table 1.
TABLE 1 parameters for pulse correction per time
Number of maneuvers | Time (h) | x(m/s) | y(m/s) | z(m/s) | Maneuver size (m/s) |
1 | 0.334 | -0.463 | -0.319 | 0.068 | 0.566 |
2 | 0.967 | -0.139 | -0.512 | -0.080 | 0.537 |
3 | 1.729 | -0.461 | -0.118 | 0.130 | 0.493 |
Fig. 3 and 4 show three-dimensional and XY plan views of the hovering orbit of the region satisfying the task constraint and having the minimum total velocity increment, and fig. 5 shows a distance change curve between a detector and the center point of the hovering region in the region hovering process, so that it can be seen that the multi-pulse region hovering method based on asteroid surface observation disclosed by the invention has strong reliability, the designed hovering orbit is smooth, and the convergence of the design method is good.
The above detailed description is intended to illustrate the objects, aspects and advantages of the present invention, and it should be understood that the above detailed description is only exemplary of the present invention and is not intended to limit the scope of the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (1)
1. A multi-pulse area hovering method based on asteroid surface observation comprises the following steps:
the method comprises the following steps: establishing a minor planet body coordinate system, and establishing a kinetic equation by considering the influence of the gravity of the minor planet and the sunlight pressure;
the centroid of the target asteroid is taken as the origin of coordinates, the equatorial plane of the asteroid is taken as an XY plane, the X axis points to the specific direction of the surface of the asteroid, the Z axis is in right-hand relation with the X axis and the Y axis, and the Z axis points upwards along the rotation axis of the asteroid; the position vector of the detector under the asteroid system is rB=[x,y,z]TVelocity vector is vB=[vx,vy,vz]TThen, the motion equation of the detector under the asteroid system is as follows:
wherein x, y and z are three-axis position components of the detector under the minor planet body coordinate system respectively, and vx,vy,vzThe three-axis velocity components of the detector under the minor planet body coordinate system are respectively, omega is the minor planet self-rotation angular velocity vector, aastIs the minor planet gravitational acceleration vector of that position, aSRPPerturbing an acceleration vector for sunlight pressure received by a detector;
establishing a gravitational field of a small planet by adopting a polyhedral model; then there isExpressing the equation of motion as a scalar:
wherein U is asteroid potential energy, aSRPx、aSRPyAnd aSRPzAre respectively provided withThe three-axis component of the sunlight pressure perturbation force in the minor planet body coordinate system is shown;
step two: selecting a center point of the hovering area according to the task requirement, wherein the center point is set as rtarDetermining the radius of the hovering area and the minimum hovering duration T with the radius of the hovering area being rhof;
In the hovering stage of the area based on the asteroid surface observation, the detector needs to investigate and screen target points of subsequent tasks such as landing sampling and the like, and the thermal infrared camera needs to be used for determining the asteroid surface temperature to ensure the operation safety of the detector in the asteroid close range, so the center point of the hovering area is generally selected to be above the most concerned position point of the asteroid surface in the task and is set as rtarThe size of the hovering area is determined by comprehensively considering factors such as the geographic environment near a target point, the height of the detector from the asteroid surface, the size of the angle of view of a camera carried by the detector and the like according to task requirements, the hovering area is generally selected to be ten meters to hundred meters, and if the radius of the hovering area of the asteroid surface is determined to be rho, the movable range of the detector is determined to be centered at rtarA sphere of radius ρ, denoted rB-rtar|≤ρ (4)
Hovering duration T required by detector at taskfThe internal demand time meets the space position requirement determined by the formula;
step three: giving the initial speed of the detector, and integrating the initial speed to the boundary of the hovering area in the minor planet body coordinate system kinetic equation from the center position of the hovering area to obtain the initial correction time t11And a corresponding pre-correction state X (t)11);
Given the initial time t of the detector0Velocity v at target location point 00=[vx0,vy0,vz0]Then the initial state of the detector is X0=[rtar,v0]' integrating the initial state in a minor planet body coordinate system kinetic equation, and calculating the distance rho between the position of the detector and the central point of the target hovering area in real timet(ii) a Influenced by the existence of relative speed at initial hovering point and the irregular gravitational field of asteroidWill gradually deviate from the center point of the hovering area and reach the boundary of the hovering area along with the time, which is set at the moment t11The detector drifts to the hover region boundary for the first time, at which time ptRho, and the state before detector correction is denoted as X (t)11);
Step four: after the boundary of the hovering area is reached for the first time, giving correction transfer time, applying pulse control by adopting a targeting method, determining a boundary correction pulse with optimal fuel consumption by taking the minimum speed increment as a performance index, and enabling the detector to reach the central position of the hovering area;
the detector state at the boundary is X (t)11) Giving a corrected transfer time tmAnd correcting the pulse maneuver Deltav to obtain the corrected terminal state r of the detectorf1;
The pulse correction orbit is designed by adopting a targeting method, and a targeting equation of the regional hovering problem based on asteroid surface observation is established as
F=rf1-rtar(5)
At this time, the correction problem is converted into a solution problem of a nonlinear equation system, and the detector is expected to be at the moment t1=t11+tmReacquires the hover region center position rtar(ii) a Meanwhile, considering the limitation of propellant carried by the detector, the optimal fuel consumption is selected as a performance index, namely, the pulse motor with the minimum speed increment is selected as a control correction amount under the condition of sufficient target hitting times, and the pulse motor is recorded as delta v1;
Step five: repeating the third step and the fourth step, and determining the time from the center point to the boundary of the hovering area and the pulse maneuver each time;
the first pulse modification makes the detector at t1=t11+tmThe moment of time reaches the hovering area center point r from the hovering area boundarytarWhen the detector state is marked as X (t)1)=[rtar,v1](ii) a At this time | v1If | ≠ 0, the detector and the asteroid still have relative speed, so that the detector can continuously deviate from the expected suspension point, and when the boundary of the hovering area is reached again, the next correction pulse is applied by using a targeting method;
remember the ith fromHovering area center position rtarThe time length for reaching the boundary of the hovering area is ti1I-th arrival at rtarAt a time tiI-th time from hover region boundary back to rtarThe desired pulse maneuver is Δ vi;
Step six: determining the number of times of pulse maneuvers required to be applied and the corresponding optimal total pulse speed increment according to the minimum hovering duration required by the task;
first i times when the detector reaches the center point r of the hovering areatarIs ti
The task requires a hover time period of TfWhen t isi-1<TfAnd t isi>TfWhen the region hovers to meet the task requirement, the total speed increment required by the region hovering is
And obtaining the region hovering orbit based on the asteroid surface observation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911149845.6A CN110758775B (en) | 2019-11-21 | 2019-11-21 | Multi-pulse area hovering method based on asteroid surface observation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911149845.6A CN110758775B (en) | 2019-11-21 | 2019-11-21 | Multi-pulse area hovering method based on asteroid surface observation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110758775A true CN110758775A (en) | 2020-02-07 |
CN110758775B CN110758775B (en) | 2021-01-08 |
Family
ID=69339124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911149845.6A Active CN110758775B (en) | 2019-11-21 | 2019-11-21 | Multi-pulse area hovering method based on asteroid surface observation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110758775B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111332498A (en) * | 2020-01-10 | 2020-06-26 | 北京理工大学 | Asteroid balance point hovering detection constant thrust threshold control method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0533489A1 (en) * | 1991-09-19 | 1993-03-24 | Nec Corporation | Satellite stabilized by gravity gradient |
RU2289533C1 (en) * | 2005-04-29 | 2006-12-20 | Государственное образовательное учреждение высшего профессионального образования Московский авиационный институт (технический университет) | Method of injection of spacecraft into inter-planetary trajectory |
US20070219813A1 (en) * | 2006-03-17 | 2007-09-20 | Moore Barrett H | Purchase option-based emergency supplies provisioning method |
US7728756B2 (en) * | 2007-08-20 | 2010-06-01 | Raytheon Company | Wide area high resolution SAR from a moving and hovering helicopter |
CN106446313A (en) * | 2015-08-10 | 2017-02-22 | 北京空间飞行器总体设计部 | Polar hovering satellite orbit-based system design method |
CN107688351A (en) * | 2017-08-31 | 2018-02-13 | 北京理工大学 | A kind of pulsion phase of spacecraft two is to method of hovering |
CN107992682A (en) * | 2017-12-05 | 2018-05-04 | 北京理工大学 | A kind of optimal multiple-pulse transfer method of interplanetary multi-body system asteroid detection |
-
2019
- 2019-11-21 CN CN201911149845.6A patent/CN110758775B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0533489A1 (en) * | 1991-09-19 | 1993-03-24 | Nec Corporation | Satellite stabilized by gravity gradient |
RU2289533C1 (en) * | 2005-04-29 | 2006-12-20 | Государственное образовательное учреждение высшего профессионального образования Московский авиационный институт (технический университет) | Method of injection of spacecraft into inter-planetary trajectory |
US20070219813A1 (en) * | 2006-03-17 | 2007-09-20 | Moore Barrett H | Purchase option-based emergency supplies provisioning method |
US7728756B2 (en) * | 2007-08-20 | 2010-06-01 | Raytheon Company | Wide area high resolution SAR from a moving and hovering helicopter |
CN106446313A (en) * | 2015-08-10 | 2017-02-22 | 北京空间飞行器总体设计部 | Polar hovering satellite orbit-based system design method |
CN107688351A (en) * | 2017-08-31 | 2018-02-13 | 北京理工大学 | A kind of pulsion phase of spacecraft two is to method of hovering |
CN107992682A (en) * | 2017-12-05 | 2018-05-04 | 北京理工大学 | A kind of optimal multiple-pulse transfer method of interplanetary multi-body system asteroid detection |
Non-Patent Citations (1)
Title |
---|
李翔宇: "小行星探测近轨操作的轨道动力学与控制", 《中国科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111332498A (en) * | 2020-01-10 | 2020-06-26 | 北京理工大学 | Asteroid balance point hovering detection constant thrust threshold control method |
CN111332498B (en) * | 2020-01-10 | 2021-10-19 | 北京理工大学 | Asteroid balance point hovering detection constant thrust threshold control method |
Also Published As
Publication number | Publication date |
---|---|
CN110758775B (en) | 2021-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105911866B (en) | The finite time full-order sliding mode control method of quadrotor unmanned vehicle | |
CN110794863B (en) | Heavy carrier rocket attitude control method capable of customizing control performance indexes | |
CN105159305B (en) | A kind of quadrotor flight control method based on sliding moding structure | |
CN109911249B (en) | Interstellar transfer limited thrust orbit-entering iterative guidance method for low thrust-weight ratio aircraft | |
CN110329546B (en) | Small celestial body landing track optimization method considering gravitational attitude and orbit coupling effect | |
CN105574261B (en) | A kind of moon borrows the ground moon libration point transfer orbit design method of force constraint | |
CN107491088B (en) | Airship track control method with saturated input | |
CN106970530B (en) | Model-free preset performance control method for autonomous sight intersection of space non-cooperative targets | |
CN111399531A (en) | Hypersonic aircraft glide phase guidance and attitude control integrated design method | |
CN110466805B (en) | Planet landing guidance method based on optimized guidance parameters | |
CN110334439B (en) | Aerodynamic force assisted superorbital speed flight technical method | |
CN113467241B (en) | Method for optimizing burn-up of convex curvature landing track | |
Peng et al. | Free return orbit design and characteristics analysis for manned lunar mission | |
CN109582039A (en) | A kind of optimal reconfiguration method under the J2 perturbation using Relative Navigation information | |
CN110758775B (en) | Multi-pulse area hovering method based on asteroid surface observation | |
CN114485672B (en) | Planar detector coupling constraint track planning method for attaching small celestial body | |
CN108459611B (en) | Attitude tracking control method of near space vehicle | |
CN108082538B (en) | Multi-body system low-energy track capturing method considering initial and final constraints | |
CN117471952A (en) | Integrated control method for backstepping supercoiled sliding mode guidance of aircraft | |
CN112560343A (en) | J2 perturbation Lambert problem solving method based on deep neural network and targeting algorithm | |
CN107688351B (en) | Two-pulse relative hovering method for spacecraft | |
CN115892519A (en) | Spacecraft control method for short-distance spacecraft orbit pulse game | |
CN108891625B (en) | Combined control method for solid micro-propeller array and magnetic torquer | |
CN107703967B (en) | Control method for controlling track of limited airship | |
CN106125757B (en) | A kind of drive lacking dirigible Track In Track control method |
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 |