A method is proposed to select the target sequence for a J 2-perturbed multiple debris rendezvous mission aimed at removing dozens of debris from several thousand debris candidates running on sun-synchronous orbits(SS...A method is proposed to select the target sequence for a J 2-perturbed multiple debris rendezvous mission aimed at removing dozens of debris from several thousand debris candidates running on sun-synchronous orbits(SSO).The solving methodology proceeds in two steps:Firstly,the variance of the right ascension of ascending node(RAAN)of the debris group is used for narrowing down the potential debris candidate;secondly,the debris of the candidate group that has closest RAAN to the current debris is chosen as the next debris.The low thrust near-minimum-fuel trajectories of each rendezvous leg are obtained by the indirect optimization method.The proposed approach is demonstrated for the problem of the 8th China Trajectory Optimization Competition(CTOC).The radar cross section(RCS)of the debris is also considered in the first step since the primary performance index of the competition is to maximize the total RCS of the debris visited.The results show that the proposed approach achieves better performance within a competition period.Of the many rendezvous sequences found,the best one submitted for the competition obtained a total RCS of 184 by accomplishing rendezvous with 70 debris within a transfer duration of one year.展开更多
High-specific-impulse electric propulsion technology is promising for future space robotic debris removal in sun-synchronous orbits.Such a prospect involves solving a class of challenging problems of low-thrust orbita...High-specific-impulse electric propulsion technology is promising for future space robotic debris removal in sun-synchronous orbits.Such a prospect involves solving a class of challenging problems of low-thrust orbital rendezvous between an active spacecraft and a free-flying debris.This study focuses on computing optimal low-thrust minimum-time many-revolution trajectories,considering the effects of the Earth oblateness perturbations and null thrust in Earth shadow.Firstly,a set of mean-element orbital dynamic equations of a chaser(spacecraft)and a target(debris)are derived by using the orbital averaging technique,and specifically a slow-changing state of the mean longitude difference is proposed to accommodate to the rendezvous problem.Subsequently,the corresponding optimal control problem is formulated based on the mean elements and their associated costate variables in terms of Pontryagin’s maximum principle,and a practical optimization procedure is adopted to find the specific initial costate variables,wherein the necessary conditions of the optimal solutions are all satisfied.Afterwards,the optimal control profile obtained in mean elements is then mapped into the counterpart that is employed by the osculating orbital dynamics.A simple correction strategy about the initialization of the mean elements,specifically the differential mean true longitude,is suggested,which is capable of minimizing the terminal orbital rendezvous errors for propagating orbital dynamics expressed by both mean and osculating elements.Finally,numerical examples are presented,and specifically,the terminal orbital rendezvous accuracy is verified by solving hundreds of rendezvous problems,demonstrating the effectiveness of the optimization method proposed in this article.展开更多
Low thrust propulsion and gravity assist (GA) are among the most promising techniques for deep space explorations.In this paper the two techniques are combined and treated comprehensively,both on modeling and numerica...Low thrust propulsion and gravity assist (GA) are among the most promising techniques for deep space explorations.In this paper the two techniques are combined and treated comprehensively,both on modeling and numerical techniques.Fuel optimal orbit rendezvous via multiple GA is first formulated as optimal guidance with multiple interior constraints and then the optimal necessary conditions,various transversality conditions and stationary conditions are derived by Pontryagin's Maximum Principle (PMP).Finally the initial orbit rendezvous problem is transformed into a multiple point boundary value problem (MPBVP).Homotopic technique combined with random searching globally and Particle Swarm Optimization (PSO),is adopted to handle the numerical difficulty in solving the above MPBVP by single shooting method.Two scenarios in the end show the merits of the present approach.展开更多
The low-thrust trajectory optimization with complicated constraints must be considered in practical engineering. In most literature, this problem is simplified into a two-body model in which the spacecraft is subject ...The low-thrust trajectory optimization with complicated constraints must be considered in practical engineering. In most literature, this problem is simplified into a two-body model in which the spacecraft is subject to the gravitational force at the center of mass and the spacecraft's own electric propulsion only, and the gravity assist (GA) is modeled as an instantaneous velocity increment. This paper presents a method to solve the fuel-optimal problem of low-thrust trajectory with complicated constraints in a full ephemeris model, which is closer to practical engineering conditions. First, it introduces various perturbations, including a third body's gravity, the nonspherical perturbation and the solar radiation pressure in a dynamic equation. Second, it builds two types of equivalent inner constraints to describe the GA. At the same time, the present paper applies a series of techniques, such as a homotopic approach, to enhance the possibility of convergence of the global optimal solution.展开更多
A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits.First,through parameterizing the control steering angles by costate variables,the s...A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits.First,through parameterizing the control steering angles by costate variables,the search space of free parameters has been decreased.Then,in order to obtain the global optimal solution effectively and robustly,the simulated annealing and penalty function strategies were used to handle the constraints,and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem,in which,a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement.Comparing to the classical direct method,this novel method has fewer free parameters,needs not initial guesses,and has higher computation precision.An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach.The results of simulation indicate that our approach is available to solve the problem of optimal multi-revolution transfer between Earth-orbits.展开更多
Spacecraft science missions to planets or asteroids have historically visited only one or several celestial bodies per mission.The research goal of this paper is to create a trajectory design algorithm that generates ...Spacecraft science missions to planets or asteroids have historically visited only one or several celestial bodies per mission.The research goal of this paper is to create a trajectory design algorithm that generates trajectory allowing a spacecraft to visit a significant number of asteroids during a single mission.For the problem of global trajectory optimization,even with recent advances in low-thrust trajectory optimization,a full enumeration of this problem is not possible.This work presents an algorithm to traverse the searching space in a practical fashion and generate solutions.The flight sequence is determined in ballistic scenario,and a differential evolution method is used with constructing a three-impulse transfer problem,then the local optimization is implemented with low-thrust propulsion on the basis of the solutions of impulsive trajectories.The proposed method enables trajectory design for multiple asteroids tour by using available low thrust propulsion technology within fuel and time duration constraints.展开更多
This paper studies the existence and stability of the artificial equilibrium points (AEPs) in the low-thrust restricted three-body problem when both the primaries are oblate spheroids. The artificial equilibrium point...This paper studies the existence and stability of the artificial equilibrium points (AEPs) in the low-thrust restricted three-body problem when both the primaries are oblate spheroids. The artificial equilibrium points (AEPs) are generated by canceling the gravitational and centrifugal forces with continuous low-thrust at a non-equilibrium point. Some graphical investigations are shown for the effects of the relative parameters which characterized the locations of the AEPs. Also, the numerical values of AEPs have been calculated. The positions of these AEPs will depend not only also on magnitude and directions of low-thrust acceleration. The linear stability of the AEPs has been investigated. We have determined the stability regions in the xy, xz and yz-planes and studied the effect of oblateness parameters A1(0A1?and ?A2(0A2<1) on the motion of the spacecraft. We have found that the stability regions reduce around both the primaries for the increasing values of oblateness of the primaries. Finally, we have plotted the zero velocity curves to determine the possible regions of motion of the spacecraft.展开更多
This paper proposes an optimal,robust,and efficient guidance scheme for the perturbed minimum-time low-thrust transfer toward the geostationary orbit.The Earth’s oblateness perturbation and shadow are taken into acco...This paper proposes an optimal,robust,and efficient guidance scheme for the perturbed minimum-time low-thrust transfer toward the geostationary orbit.The Earth’s oblateness perturbation and shadow are taken into account.It is difficult for a Lyapunov-based or trajectory-tracking guidance method to possess multiple characteristics at the same time,including high guidance optimality,robustness,and onboard computational efficiency.In this work,a concise relationship between the minimum-time transfer problem with orbital averaging and its optimal solution is identified,which reveals that the five averaged initial costates that dominate the optimal thrust direction can be approximately determined by only four initial modified equinoctial orbit elements after a coordinate transformation.Based on this relationship,the optimal averaged trajectories constituting the training dataset are randomly generated around a nominal averaged trajectory.Five polynomial regression models are trained on the training dataset and are regarded as the costate estimators.In the transfer,the spacecraft can obtain the real-time approximate optimal thrust direction by combining the costate estimations provided by the estimators with the current state at any time.Moreover,all these computations onboard are analytical.The simulation results show that the proposed guidance scheme possesses extremely high guidance optimality,robustness,and onboard computational efficiency.展开更多
Spanwise flexibility is a key factor influencing propulsion performance of pectoral foils. Performances of bionic fish with oscillating pectoral foils can be enhanced by properly selecting the spanwise flexibility. Th...Spanwise flexibility is a key factor influencing propulsion performance of pectoral foils. Performances of bionic fish with oscillating pectoral foils can be enhanced by properly selecting the spanwise flexibility. The influence law of spanwise flexibility on thrust generation and propulsion efficiency of a rectangular hydro-foil is discussed. Series foils constructed by the two-component silicon rubber are developed. NACA0015 shape of chordwise cross-section is employed. The foils are strengthened by fin rays of different rigidity to realize variant spanwise rigidity and almost the same chordwise flexibility. Experiments on a towing platform developed are carried out at low Reynolds numbers of 10 000, 15 000, and 20 000 and Strouhal numbers from 0.1 to 1. The following experimental results are achieved: (1) The average forward thrust increases with the St number increased; (2) Certain degree of spanwise flexibility is beneficial to the forward thrust generation, but the thrust gap is not large for the fins of different spanwise rigidity; (3) The fin of the maximal spanwise flexibility owns the highest propulsion efficiency; (4) Effect of the Reynolds number on the propulsion efficiency is significant. The experimental results can be utilized as a reference in deciding the spanwise flexibility of bionic pectoral fins in designing of robotic fish prototype propelled by flapping-wing.展开更多
基金We are very grateful to the organizers of the 8th China Trajectory Optimization Competition for the interesting and complex problemMost methods presented in this paper were developed under the National Natural Science Foundation of China(Nos.11172036,11572037,and 11402021)the Excellent Young Scholars Research Fund of Beijing Institute of Technology(No.2015YG0101).
文摘A method is proposed to select the target sequence for a J 2-perturbed multiple debris rendezvous mission aimed at removing dozens of debris from several thousand debris candidates running on sun-synchronous orbits(SSO).The solving methodology proceeds in two steps:Firstly,the variance of the right ascension of ascending node(RAAN)of the debris group is used for narrowing down the potential debris candidate;secondly,the debris of the candidate group that has closest RAAN to the current debris is chosen as the next debris.The low thrust near-minimum-fuel trajectories of each rendezvous leg are obtained by the indirect optimization method.The proposed approach is demonstrated for the problem of the 8th China Trajectory Optimization Competition(CTOC).The radar cross section(RCS)of the debris is also considered in the first step since the primary performance index of the competition is to maximize the total RCS of the debris visited.The results show that the proposed approach achieves better performance within a competition period.Of the many rendezvous sequences found,the best one submitted for the competition obtained a total RCS of 184 by accomplishing rendezvous with 70 debris within a transfer duration of one year.
基金supported by the National Key Research and Development Project(Grant No.2018YFB1900605)the Key Research Program of Chinese Academy of Sciences(Grant No.ZDRW-KT-2019-1).
文摘High-specific-impulse electric propulsion technology is promising for future space robotic debris removal in sun-synchronous orbits.Such a prospect involves solving a class of challenging problems of low-thrust orbital rendezvous between an active spacecraft and a free-flying debris.This study focuses on computing optimal low-thrust minimum-time many-revolution trajectories,considering the effects of the Earth oblateness perturbations and null thrust in Earth shadow.Firstly,a set of mean-element orbital dynamic equations of a chaser(spacecraft)and a target(debris)are derived by using the orbital averaging technique,and specifically a slow-changing state of the mean longitude difference is proposed to accommodate to the rendezvous problem.Subsequently,the corresponding optimal control problem is formulated based on the mean elements and their associated costate variables in terms of Pontryagin’s maximum principle,and a practical optimization procedure is adopted to find the specific initial costate variables,wherein the necessary conditions of the optimal solutions are all satisfied.Afterwards,the optimal control profile obtained in mean elements is then mapped into the counterpart that is employed by the osculating orbital dynamics.A simple correction strategy about the initialization of the mean elements,specifically the differential mean true longitude,is suggested,which is capable of minimizing the terminal orbital rendezvous errors for propagating orbital dynamics expressed by both mean and osculating elements.Finally,numerical examples are presented,and specifically,the terminal orbital rendezvous accuracy is verified by solving hundreds of rendezvous problems,demonstrating the effectiveness of the optimization method proposed in this article.
基金supported by the National Natural Science Foundation of China(Grant Nos. 10832004 and 11072122)
文摘Low thrust propulsion and gravity assist (GA) are among the most promising techniques for deep space explorations.In this paper the two techniques are combined and treated comprehensively,both on modeling and numerical techniques.Fuel optimal orbit rendezvous via multiple GA is first formulated as optimal guidance with multiple interior constraints and then the optimal necessary conditions,various transversality conditions and stationary conditions are derived by Pontryagin's Maximum Principle (PMP).Finally the initial orbit rendezvous problem is transformed into a multiple point boundary value problem (MPBVP).Homotopic technique combined with random searching globally and Particle Swarm Optimization (PSO),is adopted to handle the numerical difficulty in solving the above MPBVP by single shooting method.Two scenarios in the end show the merits of the present approach.
文摘The low-thrust trajectory optimization with complicated constraints must be considered in practical engineering. In most literature, this problem is simplified into a two-body model in which the spacecraft is subject to the gravitational force at the center of mass and the spacecraft's own electric propulsion only, and the gravity assist (GA) is modeled as an instantaneous velocity increment. This paper presents a method to solve the fuel-optimal problem of low-thrust trajectory with complicated constraints in a full ephemeris model, which is closer to practical engineering conditions. First, it introduces various perturbations, including a third body's gravity, the nonspherical perturbation and the solar radiation pressure in a dynamic equation. Second, it builds two types of equivalent inner constraints to describe the GA. At the same time, the present paper applies a series of techniques, such as a homotopic approach, to enhance the possibility of convergence of the global optimal solution.
基金Sponsored by the National Natural Science Foundation of China(Grant No.10672044)
文摘A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits.First,through parameterizing the control steering angles by costate variables,the search space of free parameters has been decreased.Then,in order to obtain the global optimal solution effectively and robustly,the simulated annealing and penalty function strategies were used to handle the constraints,and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem,in which,a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement.Comparing to the classical direct method,this novel method has fewer free parameters,needs not initial guesses,and has higher computation precision.An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach.The results of simulation indicate that our approach is available to solve the problem of optimal multi-revolution transfer between Earth-orbits.
文摘Spacecraft science missions to planets or asteroids have historically visited only one or several celestial bodies per mission.The research goal of this paper is to create a trajectory design algorithm that generates trajectory allowing a spacecraft to visit a significant number of asteroids during a single mission.For the problem of global trajectory optimization,even with recent advances in low-thrust trajectory optimization,a full enumeration of this problem is not possible.This work presents an algorithm to traverse the searching space in a practical fashion and generate solutions.The flight sequence is determined in ballistic scenario,and a differential evolution method is used with constructing a three-impulse transfer problem,then the local optimization is implemented with low-thrust propulsion on the basis of the solutions of impulsive trajectories.The proposed method enables trajectory design for multiple asteroids tour by using available low thrust propulsion technology within fuel and time duration constraints.
文摘This paper studies the existence and stability of the artificial equilibrium points (AEPs) in the low-thrust restricted three-body problem when both the primaries are oblate spheroids. The artificial equilibrium points (AEPs) are generated by canceling the gravitational and centrifugal forces with continuous low-thrust at a non-equilibrium point. Some graphical investigations are shown for the effects of the relative parameters which characterized the locations of the AEPs. Also, the numerical values of AEPs have been calculated. The positions of these AEPs will depend not only also on magnitude and directions of low-thrust acceleration. The linear stability of the AEPs has been investigated. We have determined the stability regions in the xy, xz and yz-planes and studied the effect of oblateness parameters A1(0A1?and ?A2(0A2<1) on the motion of the spacecraft. We have found that the stability regions reduce around both the primaries for the increasing values of oblateness of the primaries. Finally, we have plotted the zero velocity curves to determine the possible regions of motion of the spacecraft.
基金supported by the National Natural Science Foundation of China(No.12022214)the National Key R&D Program of China(No.2020YFC2201200)。
文摘This paper proposes an optimal,robust,and efficient guidance scheme for the perturbed minimum-time low-thrust transfer toward the geostationary orbit.The Earth’s oblateness perturbation and shadow are taken into account.It is difficult for a Lyapunov-based or trajectory-tracking guidance method to possess multiple characteristics at the same time,including high guidance optimality,robustness,and onboard computational efficiency.In this work,a concise relationship between the minimum-time transfer problem with orbital averaging and its optimal solution is identified,which reveals that the five averaged initial costates that dominate the optimal thrust direction can be approximately determined by only four initial modified equinoctial orbit elements after a coordinate transformation.Based on this relationship,the optimal averaged trajectories constituting the training dataset are randomly generated around a nominal averaged trajectory.Five polynomial regression models are trained on the training dataset and are regarded as the costate estimators.In the transfer,the spacecraft can obtain the real-time approximate optimal thrust direction by combining the costate estimations provided by the estimators with the current state at any time.Moreover,all these computations onboard are analytical.The simulation results show that the proposed guidance scheme possesses extremely high guidance optimality,robustness,and onboard computational efficiency.
基金supported by National Hi-tech Research and Development Program of China(863 Program, Grant No. 2006AA04Z252)National Natural Science Foundation of China(Grant No. 51005006)+1 种基金Research Fund for the Doctoral Program of Higher Education of China(Grand No. 20101102110022)Innovation Foundation of Beihang University for PhD Graduates, China
文摘Spanwise flexibility is a key factor influencing propulsion performance of pectoral foils. Performances of bionic fish with oscillating pectoral foils can be enhanced by properly selecting the spanwise flexibility. The influence law of spanwise flexibility on thrust generation and propulsion efficiency of a rectangular hydro-foil is discussed. Series foils constructed by the two-component silicon rubber are developed. NACA0015 shape of chordwise cross-section is employed. The foils are strengthened by fin rays of different rigidity to realize variant spanwise rigidity and almost the same chordwise flexibility. Experiments on a towing platform developed are carried out at low Reynolds numbers of 10 000, 15 000, and 20 000 and Strouhal numbers from 0.1 to 1. The following experimental results are achieved: (1) The average forward thrust increases with the St number increased; (2) Certain degree of spanwise flexibility is beneficial to the forward thrust generation, but the thrust gap is not large for the fins of different spanwise rigidity; (3) The fin of the maximal spanwise flexibility owns the highest propulsion efficiency; (4) Effect of the Reynolds number on the propulsion efficiency is significant. The experimental results can be utilized as a reference in deciding the spanwise flexibility of bionic pectoral fins in designing of robotic fish prototype propelled by flapping-wing.