Earth-moon Trajectory Optimization Using Solar Electric Propulsion
|
|
- Simon Ellis
- 5 years ago
- Views:
Transcription
1 Chinese Journal o Aeronautics (7) Chinese Journal o Aeronautics arth-moon rajectory Optimization Using Solar lectric Propulsion Gao Yang* Academy o Opto-lectronics Chinese Academy o Sciences Beijing 18 China Received 5 April 7; accepted 17 August 7 Abstract he optimization o the arth-moon trajectory using solar electric propulsion is presented. A easible method is proposed to optimize the transer trajectory starting rom a low arth circular orbit (5 km altitude) to a low lunar circular orbit ( km altitude). Due to the use o low-thrust solar electric propulsion the entire transer trajectory consists o hundreds or even thousands o orbital revolutions around the arth and the moon. he arth-orbit ascending (rom low arth orbit to high arth orbit) and lunar descending (rom high lunar orbit to low lunar orbit) trajectories in the presence o J perturbations and shadowing eect are computed by an analytic orbital averaging technique. A direct/indirect method is used to optimize the control steering or the trans-lunar trajectory segment a segment rom a high arth orbit to a high lunar orbit with a ixed thrust-coast-thrust engine sequence. For the trans-lunar trajectory segment the equations o motion are expressed in the inertial coordinates about the arth and the moon using a set o nonsingular equinoctial elements inclusive o the gravitational orces o the sun the arth and the moon. By way o the analytic orbital averaging technique and the direct/indirect method the arth-moon transer problem is converted to a parameter optimization problem and the entire transer trajectory is ormulated and optimized in the orm o a single nonlinear optimization problem with a small number o variables and constraints. Finally an example o an arth-moon transer trajectory using solar electric propulsion is demonstrated. Keywords: trajectory optimization; solar electric propulsion; analytic orbital averaging technique; direct/indirect method 1 Introduction * Since the Deep Space 1 spacecrat demonstrated the irst use o solar electric propulsion (SP) or an interplanetary mission [1] the application o low-thrust propulsion or uture space missions has been a popular research subject. More recently the SA Smart-1 spacecrat [] launched in 3 successully perormed the irst lunar mission using solar electric propulsion. It is well known that the spacecrat propelled by low-thrust SP engines is capable o delivering a greater payload raction *Corresponding author. el.: mail address: gaoy@aoe.ac.cn Foundation item: National Natural Science Foundation o China (1635) compared to the spacecrat using conventional chemical propulsion. However the low-thrust continuous control proiles can not be approximated by the impulsive velocity dierences. How to optimize low-thrust transer trajectories becomes a new challenge to the mission designers who want to use SP as the primary spacecrat propulsion. In the past two decades lots o researches on low-thrust arth-moon trajectories were conducted. For instance Golan and Breakwell [3] investigated minimum-uel lunar trajectories with the ixed transer time by patching arth- and moon-centered spirals at an intermediate point. Huelman [4] introduced a power-limited optimal guidance law or arth-moon transers in a planar restricted three-
2 Gao Yang / Chinese Journal o Aeronautics (7) body dynamics. Piersion and Kluever [5] solved optimal arth-moon trajectories using a three-stage approach. Later they obtained a series o arthmoon trajectories [6-8] in the classical restricted threebody dynamics. Herman and Conway [9] employed a parallel Runge-Kutta method to solve optimal arth-moon transers which start rom a circular orbit with 6.3 arth radii and end at a two-lunar-radii orbit. he initial thrust acceleration o 1 4 g (g is the arth sea-level gravitational acceleration) yields a 3-day light journey. Recently Betts [1] obtained an optimal trajectory rom a geostationary transer orbit (GO) to a high elliptic lunar orbit with the perilune radius o 378 km and the apolune radius o km using the direct transcription method or collocation method which creates a large-scale nonlinear optimization problem including variables and constraints. A less accurate approximate solution to this problem was obtained with the computational time o s (11.67 h) and the longer time is needed to obtain a more accurate solution. However complete transers rom low arth orbits (LO) to low lunar orbits (LLO) using solar electric propulsion with or without consideration o signiicant shadowing eect still have not been solved in above-mentioned papers. In Res.[7] and [8] Kluever and Piersion used the delbaum s method to approximate the spirals rom LO to a high arth orbit (HO) and rom a high lunar orbit (HLO) to LLO. Nevertheless the delbaum s method does not consider J perturbations and shadowing condition that signiicantly aects trajectory evolution in relatively low-altitude orbits around the arth and the moon. I Betts direct transcription and the parallel Runge-Kutta method [9] are utilized to solve a long-duration transer trajectory rom LO to LLO the dimension o the nonlinear optimization problem would be much larger. In addition it is diicult to consider the shadowing eect in the existing direct methods. his paper demonstrates a easible method to compute the low-thrust arth-moon trajectory rom a low earth circular orbit (5 km altitude) to a low lunar circular orbit ( km altitude). he transer rom HO to HLO is computed by a direct/indirect method utilized by Kluever and Piersion [5-8] but the equinoctial elements in a our-body (the arth the moon the sun and the spacecrat) dynamics including J perturbations are employed. Furthermore the low-thrust spirals rom LO to HO and rom HLO to LLO are computed by an analytic orbital averaging technique (AOA) in which the thirdbody perturbations are ignored but the signiicant J perturbations and shadowing eect are taken into account. he ephemeris o the arth the moon and the sun are obtained using JPL planetary and lunar ephemerides [11]. he moon s orbit around the arth is not a planar circular orbit; the moon s inclination varies between 18 and 8 every 15 years. he entire trajectory is assumed to be a burn-coast-burn sequence which is a good tradeo between uel consumption and transer time. he optimal control problem or the arth-moon transer trajectory is converted to a parameter optimization problem that is in turn solved by nonlinear programming sequential quadratic programming (SQP). he entire transer trajectory is optimized by a single nonlinear SQP problem including only a ew variables and constraints. quations o Motion.1 Deinitions o CI and MCI coordinates he CI coordinate (arth-centered inertial coordinate O e x e y e z e see Fig.1) reers to the J arth equatorial rame that is deined by the mean orientation o the arth s equator and ecliptic orbit at the beginning o the year. he O e x e y e plane is parallel to the mean arth s equator. he line ormed by intersection o the ecliptic orbit plane and the arth s equatorial plane deines the axis x e. he line on the irst day o autumn (starting rom the sun to the center o the arth) deines the positive direction o the axis x e which is called the vernal equinox. he axis z e is perpendicular to the O e x e y e plane and points to the north. he axis y e completes the Cartesian coordinate using the right-handed
3 454 Gao Yang / Chinese Journal o Aeronautics (7) principle. he MCI coordinate i.e. moon-centered inertial coordinate O m x m y m z m (Fig.1) is parallel to O e x e y e z e with dierent origins. In this paper the equations o motion in both CI and MCI coordinates are expressed by a set o equinoctial elements [1-14] which avoid singularities when the orbit s inclination and eccentricity equal to zero. Fig.1 Illustrations o CI MCI coordinates and thrust steering angles in the local rotating coordinate O R x R y R z R where the axis x R is rom the arth to the spacecrat the axis z R is perpendicular to the orbital plane and the axis y R completes the Cartesian coordinate O R x R y R z R ollowing the right-handed principle which is in the orbital plane and points to the velocity direction.. quations o motion in the CI coordinate he equations o motion in the CI coordinate or the arth-departure trajectory are x M( Jarth moon sun ) D (1) 1 m m /( gisp) where P/( gisp) () where x [ pghkl] is the vector o equinoctial elements in the CI coordinate and J arth moon and sun 1 are arth J oblateness acceleration moon s perturbation and sun s perturbation respectively in the CI coordinate. he elements in the matrices M and D are given in Appendix A. q.() is the mass low rate where is thrust amplitude and P I sp and are power speciic impulse and eiciency o the SP system respectively. he thrusting direction unit vector can be expressed in terms o local pitch and yaw steering angles sin cos cos cos sin (3) where the pitch angle () is measured rom the local horizon (the axis y R ) to the projection o the thrust vector onto the orbit plane and the yaw angle () is measured rom the orbit plane to the thrust vector (see Fig.1). he arth gravitational parameter and arth radius ( e and R e ) deine the non-dimensional units distance velocity time and acceleration in the CI coordinate as ollows e Re d e Re ve te ae Re / te R v (4) e he arth J perturbations can be expressed in the rotating coordinate O R x R y R z R J arth 3eJeR e 1( hsin L kcos L) 4 1 r (1 h k ) 1 ejere ( hsin Lkcos L)( hcos Lksin L) 4 r (1 h k ) 6 ejere (1 h k )( hsin Lkcos L) 4 r (1 h k ) (5) where r is the distance rom the arth center to the spacecrat. he gravitational accelerations o the sun and the moon expressed in the CI coordinate can be written as rsunsc r esun sun 1 sun (6) 3 3 rsunsc re sun rmsc r em moon m (7) 3 3 rmsc rem where r sunsc is the position vector rom the sun to the spacecrat r esun rom the arth to the sun r mscrom the moon to the spacecrat and r em rom the arth to the moon. he gravitational parameters o the sun and the moon are denoted by sun and m respectively. he vector sun 1 and moon in the CI coordinate should be transormed to and moon in the local rotating coordinate. sun 1.3 quations o motion in the MCI coordinate Likewise the equations o motion in the MCI coordinate or the moon-capture trajectory are e
4 Gao Yang / Chinese Journal o Aeronautics (7) where x M( Jmoon arth sun ) D (8) m x [ p g h k L] is the vector o equinoctial elements in the MCI coordinate and J moon arth and sun are moon J oblateness acceleration arth s perturbation and sun s perturbation respectively in the MCI coordinate. he elements in the matrices M and D are in the expressions as in the matrices M and D and the mass low rate is the same as in q.(). he local thrusting pitch and yaw steering angles have the same deinitions as in q.(3) but reerenced in the MCI coordinate. he moon s gravitational parameter and moon s radius ( m and R m ) also deine the non-dimensional units in the MCI coordinate as ollows m Rm d m Rm vm tm am Rm / tm R v (9) m he moon s J perturbations have the same ormulation as in q.(5) with the substitutions o J m m and R m or J e e and R e. he gravitational accelerations o the sun and the arth expressed in the MCI coordinate can be written as rsunsc r msun sun sun (1) 3 3 rsunsc rm sun m resc r me arth m (11) 3 3 resc rme where r msun is the position vector rom the moon to the sun r esc rom the arth to the spacecrat and r me rom the moon to the arth. Also sun and in the MCI coordinate should be transormed arth to and arth in the local rotating coordinate. sun 3 Five Segments o the arth-moon rajectory he entire arth-moon trajectory is divided into ive segments: (1) Burn arc rom LO to HO; () Burn arc rom HO to the start o the trans-lunar coast arc; (3) rans-lunar coast arc; (4) Burn arc rom the end o trans-lunar coast arc to HLO; (5) Burn arc rom HLO to LLO. Fig. illustrates trajectory segments and able 1 summarizes perturbations trajectory propagation methods and inertial coordinates in dierent trajectory segments. Segment (1) () (3) (4) (5) Fig. Illustration o ive trajectory segments. able 1 Dierent trajectory segments Perturbations Propagation method Coordinate arth shadow AOA CI arth J arth J moon Numerical CI gravity sun gravity integration Moon J arth gravity sun gravity Numerical integration MCI Moon shadow moon J AOA MCI It is well known that the low-thrust spirals rom LO to HO and those rom HLO to LLO consist o hundreds or even thousands o orbital revolutions which requires a ormidable amount o time to compute i trajectories are numerically integrated. Segments ()-(4) having relatively ewer orbital revolutions are easier to be numerically integrated or moderate duration. hus the AOA a ast trajectory propagation algorithm is used to compute the segments (1) and (5). For the segments (1)-(3) the trajectories in the CI coordinate are propagated orward in time while those in segments (4) and (5) in the MCI coordinate backward in time. Since the gravitational orces o the sun the arth and the moon in both CI and MCI coordinates are considered the optimization subroutine would automatically ind the proper location to switch rom the CI coordinate to the MCI coordinate. hus it is not necessary to speciy a ixed intermediate point requisite or the conventional patched conic method.
5 456 Gao Yang / Chinese Journal o Aeronautics (7) Computing arth-orbit Ascending and Lunar Descending rajectories Using AOA he tangential steering along the velocity direction is employed or the arth-orbit ascending trajectory rom LO to HO while the anti-tangential steering or the lunar descending trajectory rom HLO to LLO. ither the tangential or anti-tangential steering is the optimal control strategy to change the instantaneous rate o semi-major axis. he trajectory segments (1) and (5) computed by the AOA are accurate only or relatively low orbits about an attracting body where the J perturbations and signiicant shadowing eect are included but the insigniicant third-body perturbations are ignored. he AOA algorithm is described in Re.[15] and also briely presented in Appendix B. he primary advantage o the AOA lies in the multirevolution trajectories which can be quickly propagated while maintaining satisactory accuracy [15]. Note that the computational procedure o the mooncapture trajectory using anti-tangential thrust is opposite to that using tangential-thrust trajectory. Let the initial condition o classical orbit elements mass and time at LO be (1) (1) (1) (1) (1) (1) (1) a e i m t where the superscripts represent the segment index number and the subscript denotes the initial time o the corresponding trajectory segment. he AOA propagates the tangential-thrust trajectory and determines the orbital elements spacecrat mass and light time at a point where the semi-major axis is pre-deined by a (1) with the subscript denoting the terminal time o the corresponding trajectory segment (1) (1) (1) (1) (1) (1) (1) a e i m t Likewise or the segment (5) that is propagated backward in time i the terminal condition at LLO is (5) (5) (5) (5) (5) (5) (5) a e i m t he initial condition o the segment (5) can be obtained by backward integrating anti-tangentialthrust trajectory to a point where the semi-major axis is pre-deined by (5) a (5) (5) (5) (5) (5) (5) (5) a e i m t Note that the AOA does not consider true anomaly which is a variable to be optimized. 5 Optimization Method or rans-lunar rajectory Segments A direct/indirect method is used to obtain near optimal control steering or trans-lunar trajectory rom HLO to HLO. Note that this method is employed to solve an arth-orbit transer (Zondervan Wood and Caughey [16] ) low-thrust arth-moon transers in the classical restricted three-body dynamics (Kluever and Pierson [5-8] ) as well as interplanetary transers (Gao and Kluever [17] ). In this paper the direct/indirect method is used in terms o the equinoctial elements. According to the calculus o variation theory the Hamiltonian o optimal control problem takes the ollowing orm H M p D m (1) m g I sp where p includes J perturbations and third-body perturbations (in both CI and MCI coordinates). p g h k L [ ] is the costate vector associated with the corresponding equinoctial elements and m the costate variable associated with the spacecrat mass. he optimal control steering direction unit vector is obtained by setting H / with the constraint 1. * [ ] M M (13) aking the partial derivative o the Hamiltonian with respect to the states the costate equations are determined. H M * D x x m x M p p M (14) x x H * m M M (15) m m m It indicates that the optimal control is governed by costate variables whose dynamics are given by
6 Gao Yang / Chinese Journal o Aeronautics (7) q.(19). For simplicity the costate dynamics associated with the two-body dynamic model is used to govern the control steering to make deriving derivatives o p with respect to the states unnecessary. In act the third-body perturbations are time-varying unction without explicit derivatives with respect to equinoctial elements. However the equations o motion do include perturbations. he initial costate variables or terminal costate variables need to be guessed to satisy the necessary boundary constraints without regard to the transversality condition and time-varying condition which should be considered in a classical two-point boundary-value problem. he optimal objective unction and boundary conditions are all treated by nonlinear optimization programming. Since the costate variable associated with mass do not aect the optimal control the irst six costate variables ( m is not used) could be used. hus or the segment () the state equations are q.(1) and q.(16) and the costate equation is q.(17). P m where gisp gisp (16) H M * D x x x (17) Likewise or the segment (4) the state equations are q.(9) and q.(16) and the costate equation has the same orm as in q.(17) but reerenced in the MCI. he advantage o the direct/indirect method over the direct transcription method lies in ewer variables and constraints leading to saving on considerable time or computation. However the costate variables needed to be integrated have no intuitive physical meaning which makes optimization problem much harder to converge although some techniques such as the adjoint-control transormation [18] and the multiple-shooting technique [17] are helpul in improving convergence robustness. he initial condition o the segment () can be denoted by state variables x () and costate variables where x () expressed in equinoctial ele- () () () () () () ments [ p g h k ] is transormed rom (1) (1) (1) (1) (1) () (1) () (1) a e i and m m t t. he initial values o true longitude L () and costate variables need to be guessed. he segment (3) does not integrate the costate equations since there is no thrust during the coast arc. he initial condition o the segment (3) is the terminal condition o the segment (). In contrast the segment (4) is integrated backward and the terminal condition o the segment (4) is state variables x (4) and costate variables (4) (4) where x in equinoctial elements (4) (4) (4) (4) (4) [ p g h k ] are transormed rom (5) (5) (5) (5) (5) (4) (5) (4) (5) a e i and m m t t.he terminal values o true longitude L (4) and costate variables need to be guessed. Obviously a complete trajectory should satisy the ollowing constraint (3) (3) (3) (3) (3) (3) (3) (3) (4) (4) (4) (4) (4) (4) (4) (4) [ p g h k L m t ] [ p g h k L m t ] 6 Formulation o Nonlinear Optimization Problem In the Section 4 and Section 5 the AOA and the direct/indirect method are described and the initial/ terminal conditions are speciied or each burn segment. he next step is to ormulate a parameter optimization problem that is in turn solved by nonlinear programming sequential quadratic programming (SQP) in which the overall mission objective is to minimize the propellant consumption during the transer. he SQP variables and constraints are summarized as ollows SQP variables: 11 variables in the orward trajectory propagation or segments (1)-(3): Segment (1): LO departure date (1) (1) t LO s ; () Segment (): six initial costate variables or the () () arth-departure burn arc burn arc duration t -t () HO s longitude angle L ; (3) (3) Segment (3): coast arc duration t -t. 11 variables in the backward trajectory propagation or segments (4) and (5): Segment (4): six terminal costate variables (4) or the moon-capture burn arc burn arc duration t (4) - (4) t HLO s longitude angle (4) L ;
7 458 Gao Yang / Chinese Journal o Aeronautics (7) Segment (5): LLO departure date t (5) spacecrat (5) (5) mass at LLO m LLO s. 8 SQP equality constraints: 6 constraints: terminal states o segment (3) = initial states o segment (4); 1 constraint: date at the terminal time o segment (3) = date at the initial time o segment (4); 1 constraint: mass at the terminal time o segment (3) = mass at the initial time o segment (4). 7 Comparison with Direct ranscription Method he direct transcription method parameterizes the control history using discrete nodes which inevitably results in a large number o variables and constraints i the transer contains much more revolutions. In this paper with the help o costate equations the ormulated SQP problem involves only variables and 8 constraints which is generally regarded as a smaller dimensional nonlinear optimization. Another ocused point is the use o the AOA a semi-analytic method or trajectory propagation. Capable o propagating hundreds or even thousands o revolutions it is much aster than precise numerical integration. Compared with the direct transcription method the proposed method reduces computational time by orders especially or long-duration multi-revolution transers. his makes it possible to optimize the entire transer trajectory as a single nonlinear optimization problem. 8 Numerical Results An arth-moon trajectory rom a LO with 5 km altitude to a LLO with km altitude is presented. First the parameters o the spacecrat to be deined are: the input power o the solar electric propulsion P = 1 kw the eiciency =.65 the speciic impulse I sp = 3 3 s and the initial spacecrat mass set to be 1 kg. hus the calculated initial thrust-to-weight ratio (/m g ) equals to he parameters o the initial LO and the terminal LLO are speciied in able. he initial arth departure date is set to be Jan Note that the 9 inclination o the LLO is reerenced in the MCI coordinate deined in this paper. In act 9 might not be an ideal inclination o a LLO parking orbit or the spacecrat which should take into account the moon s obliquity angle angle between the moon s rotation axis and the arth equatorial plane or ecliptic plane and urther detailed analysis o lunar gravity ield. Falling out o the scope o this paper this problem will not be discussed. able Initial LO and terminal LLO Orbital elements Initial LO erminal LLO Semi-major axis R e R m ccentricity.1.1 Inclination/( ) Ascension o ascending node ree ree Argument o periapsis/( ) rue anomaly ree ree As described in previous sections the AOA computes trajectory segment rom the LO to a HO (a HO = 1R e ) in the CI coordinate. It goes the same way or the AOA to compute the trajectory segment rom the LLO to a HLO (a HLO = R m ) in the MCI. As shown in Re.[15] the lower the altitudes o HO and HLO the more accurate the solutions. However the low attitudes might result in more revolutions or the trans-lunar trajectory segment which will conceptually cause diiculties in converging i the indirect/direct method is used. But i the number o revolution is not excessively large or instance under 5 the indirect/direct method still works well without too many troubles guessing initial costate variables. he indirect/direct method was once used to solve a lot o transer problems with small numbers o revolutions [6-8]. On the base o previous researches and author s experience in trajectory optimization a HO = 1R e and a HLO = R m are selected in the ollowing example. he ormulated parameter optimization problem has a small number o variables and constraints. However it appears not easy to make initial guesses or SQP variables to obtain converged solutions at the irst run. As a result the ollowing steps are
8 Gao Yang / Chinese Journal o Aeronautics (7) taken to ind out the inal optimal solutions: Step 1 o generate a trajectory to rendezvous moon in the CI coordinate with a burn-coast-burn engine sequence exclusive o the lunar perturbation. In this step the lunar capture is not considered. Step o remove the second burn and adjust the duration o the irst burn and coast arcs to make the end o the coast arc in the vicinity o moon. Step 3 o generate a trajectory in the MCI coordinate to escape moon s gravity and try dierent values or SQP variables in the segments (4) and (5) to make 8 SQP equality constraints as close as possible. Step 4 o optimize the trajectories rom the LO to the LLO in a single SQP problem based on the solutions obtained rom Step 3. Ater perorming the Step 4 the obtained LO- LLO trajectories are summarized in able 3 which shows about 8 days let or the trans-lunar trajectory rom the HO to the HLO. he arth-orbit ascending trajectory takes about 161 days and the moon descending trajectory about 1 days. hrough AOA are obtained about 1 9 orbital revolutions rom the LO to the HO and about 94 revolutions rom the HLO to the LLO. It indicates that the trajectories computed by the AOA have majority o revolutions o the entire transer trajectory and take about almost hal the transer time. he time histories o semi-major axis eccentricity and inclination in both CI and MCI coordinates are presented in Figs.3-5 respectively rom which it is clear that semi-major axis increases in the CI coordinate and decreases in the MCI coordinate. From the LO to the HO the eccentricity is raised at irst to about.16 ollowed by slightly alling and then quickly increasing. In the mooncapture phase the eccentricity is about.14 at the HLO. he eccentricity does not change monotonically mainly because o the eect o shadow. Without the shadowing eect the trajectory rom the LO to the HO or rom HLO to the LLO should be a near-circular transer. he inclination change in the MCI coordinate is not signiicant since the spacecrat enters moon-capture trajectory with an optimally chosen orientation. Additionally the transer trajectories in the CI coordinate are presented in Fig.6 and the lunar capture trajectory in Fig.7. he control steering direction angles in both CI and MCI coordinates are presented in Fig.8 and Fig.9. able 3 Solutions or the LO-LLO trajectories Mission parameters Solutions LO departure date Dec Fig.3 ime history o semi-major axis rom LO to LLO. HO arrival date June 18 Duration rom LO to HO/day st burn arc duration (trans-lunar)/day 5.3 coast arc duration (trans-lunar)/day 9.49 nd burn arc duration (trans-lunar)/day 19.4 HLO arrival date Aug. 8 8 LLO arrival date Sept. 9 8 Duration rom HLO to LLO/day 1.37 otal transer time/day Final mass at LLO/kg Fig.4 ime history o eccentricity rom LO to LLO.
9 46 Gao Yang / Chinese Journal o Aeronautics (7) Fig.5 ime history o inclination rom HO to HLO. Fig.9 ime histories o control steering or trajectory rom capture to HO (MCI). 9 Conclusions Fig.6 rajectory rom HO to HLO in the CI coordinate. Fig.7 Moon-capture trajectory to HLO in the MCI coordinate. A easible method has been proposed to ind preliminary solutions or arth-moon transer trajectories rom LO to LLO using SP. he shadow eects oblateness o the arth and the moon the gravitational orces o the sun the arth and the moon are all taken into consideration and the trajectories are solved using a practical complex threedimensional dynamic model. he JPL planetary and lunar ephemerides are utilized to compute precise positions o involved celestial bodies. he signiicant shadowing eect at relatively low altitude orbits is considered to highlight the use o solar electric propulsion. By means o AOA and the direct/indirect method the optimal orbit transer problem can be converted to a parameter optimization problem that only involves a small number o SQP variables and constraints. he main advantage o the proposed method lies in its provision o a more eicient algorithm to optimize long-duration multi-revolution transer trajectories while keeping the solution satisactorily accurate. he obtained solutions can be used as a good initial guess or urther trajectory optimization problems in highidelity dynamic models. Reerences [1] Rayman M D Williams S N. Design o the irst interplanetary Fig.8 ime histories o control steering or trajectory rom HO to escape (CI). solar electric propulsion mission. Journal o Spacecrat and Rockets ; 39(4): [] Koppel C R Marchandise F stublier D et al. he SMAR-1
10 Gao Yang / Chinese Journal o Aeronautics (7) electric propulsion subsystem in light experience. AIAA Paper [3] Golan O M Breakwell J V. Minimum uel lunar trajectories or a low-thrust power-limited spacecrat. AIAA Paper [4] Guelman M. arth-to-moon transer with a limited power engine. Journal o Guidance Control and Dynamics 1995; 18(5): [5] Pierson B L Kluever C A. hree-stage approach to optimal low-thrust earth-moon trajectories. Journal o the Astronautical Sciences 1994; 4(3): [6] Kluever C A Pierson B L. Vehicle-and-trajectory optimization o nuclear electric spacecrat or lunar missions. Journal o Spacecrat and Rockets 1995; 3(1): [7] Kluever C A Pierson B L. Optimal low-thrust three-dimensional earth-moon trajectories. Journal o Guidance Control and Dynamics 1995; 18(4): [8] Kluever C A Pierson B L. Optimal earth-moon trajectories using nuclear electric propulsion. Journal o Guidance Control and Dynamics 1997; (): [9] Herman A L Bruce B A. Optimal low-thrust earth-moon orbit transer. Journal o Guidance Control and Dynamics 1998; 1(1): [1] Betts J rb S O. Optimal low thrust trajectory to the moon. SIAM Journal o Applied Dynamical System 3; (): [11] Standish M. he JPL planetary and lunar ephemerides D4/ L4. Bulletin o the American Astronomical Society 1995; 7: 13. [1] Walker M J H Ireland B Owens J. A set o modiied equinoctial orbit elements. Celestial Mechanics 1985; 36(4): [13] Betts J. Optimal interplanetary orbit transers by direct transcription. Journal o the Astronautical Sciences 1994; 4(3): [14] Kluever C A. Optimal low-thrust interplanetary trajectories by direct method techniques. Journal o the Astronautical Sciences 1997; 45(3): [15] Gao Y Kluever C A. Analytic orbital averaging technique or computing tangential-thrust trajectories. Journal o Guidance Control and Dynamics 5; 8(6): [16] Zondervan K P Wood L J Caughey K. Optimal low-thrust three-burn orbit transers with large plane changes. Journal o the Astronautical Sciences 1984; 3(3): using hybrid trajectory optimization method with multiple shooting. AIAA Paper [18] Dixon L C Bartholomew-Biggs M C. Adjoint control transormations or solving practical optimal control problems. Optimal Control Applications and Methods 1981; (4): [19] Battin R H. An introduction to the mathematics and methods o astrodynamics. AIAA ducation Serials. Washington DC: AIAA 1987: [] Neta B Vallado D. On satellite umbra/penumbra entry and exit positions. Journal o the Astronautical Sciences 1998; 46(1): Biography: Gao Yang Born in 1974 he received B.S. degree rom Beijing University o Aeronautics and Astronautics in 1997 and M.S. degree rom Chinese Academy o Sciences in. He received his doctoral degree in 3 rom the University o Missouri-Columbia USA. From 4 to 5 he was a postdoctoral research ellow in the University o Missouri-Columbia. Since June 5 he has been an associate research ellow at Academy o Opto-lectronic Chinese Academy o Sciences. His research interests include orbital mechanics spacecrat dynamics and control guidance and navigation trajectory optimization and space mission design. -mail: gaoy@aoe.ac.cn Appendix A: Matrices M and D he elements o matrices M D expressed by equinoctial elements are as ollows M M M M M M M M M M M M M M M M M M M M11 M p p w M 13 p M1 sin L p 1 M [( w1)cos L ] w [17] Gao Y Kluever C A. Low-thrust interplanetary orbit transers
11 46 Gao Yang / Chinese Journal o Aeronautics (7) p g M3 ( hsin Lkcos L) w p p 1 M31 cos L M3 [( w1)sin Lg] w p M33 ( hsin Lkcos L) w M M M 41 M4 51 M5 61 M6 p s M43 cos L w p s M53 sin L w 1 p M63 ( hsin Lkcos L) w w D [ D6] D6 p p where w1 cos L gsin L is the arth gravitational parameter and s 1h k. he equinoctial elements can be obtained in terms o classical orbital elements (a e i ): p a(1 e ) ecos( ) g esin( ) h tan( i/ )cos k tan( i/ )sin L where is true anomaly. Appendix B: Analytic Orbital Averaging echnique In this section equations are deduced in the CI coordinate which in the MCI coordinate can be deduced likewise. In order to obtain the analytic incremental changes using tangential-thrust in classical orbital elements it is preerable to start with the Gauss planetary equations [17] : da a esin a p r (B1) d t h* h* r de 1 1 psin r [( pr)cos re] (B) d t h* h* d pcos ( p r)sin r (B3) d t h* e h* e d na 1 r [ r (cos e ) (1 )sin ] (B4) dt r nae a where p a(1 e ) h* p n / a and r p/(1 ecos ). When the thrust acceleration ( r and ) is much smaller than the gravitational acceleration the derivative o eccentric anomaly with respect to time is approximated by removing the 3 thrust term d na (B5) dt r Ater transorming eccentric anomaly to true anomaly by sin 1 e sin 1 ecos cos e cos (B6) 1 ecos the derivatives o the irst ive classical orbital elements (a e i ) with respect to eccentric anomaly can be computed by dividing the Gauss planetary equations (B1)-(B4) by q.(b5): 3 da a e r sin 1 e (B7) d de a r (1 e ) sin ( cos d cos ) 1 (B8) e e e d d a cos i r (cos e) 1e d d e ( e ecos )sin (B9) he thrust is so low that it can be assumed that over an orbital arc the classical orbital elements and acceleration are kept constant. he steering that maximizes the rate o semi-major axis (i.e. tangential steering) can be derived by setting (d a/ d ) / and (d a/ d ) / sin esin 1 e cos 1e cos 1e cos (B1) he analytic expressions or integrals with respect to eccentric anomaly are: da a d 3 d in 1 cos d e (B11) de a d (1 e ) in 1 e cos d d e 1 1e cos d [ln(sin 1e cos e)] d a d 1e d e 1 1e cos sin ( ecos ) in (B1) (B13)
12 Gao Yang / Chinese Journal o Aeronautics (7) where in is acceleration amplitude. Note that it is unable to ind analytic expressions or the terms 1 1 e cos d and d. he 1 e cos ollowing approximations can be utilized or these unctions to be integrated. 1e cos 1 e (1 1 e )sin (B14) 1 1e cos 1 e cos e cos e cos 1e cos 1e c (B15) o make the unctions and their approximations as close as possible c =.8 is chosen by way o trial-and-error tests. he corresponding integrals are then obtained in the ollowing analytic orms: 1e cos d 1 e (1 1 e )(.5.5sin ) 1 e 1 e cos 1 cos d d e 1 ec.5.5sin (B16) (B17) Furthermore the averaging changes o orbital elements due to J perturbation are as ollows: da de di (B18) dt dt dt Re d 3 J ncos i dt a (1 e ) Re d 3 J n(5cos i 1) dt 4 a (1 e ) (B19) (B) he increments in classical orbital elements and time including J perturbations and cylindrical shadow [] per revolution are approximated as: dz dz d dt n t n (B) en z d (B1) ex where z = [ a e i ]. he mass loss per revolution is computed by arth shadow entrance and exit angles. P 1 m [ ( en e sin en ex e sin ex )] ( gi ) n sp (B3) With y [ z t m] and y [ z t m] the elements in the (i+1)th revolution is computed only in terms o the elements in the (i)th revolution y y y (B4) i1 i i A terminal semi-major axis can be speciied or the stop condition o the orbital averaging. he states at the terminal semi-major are obtained by interpolating the orbital elements during the last revolution a a y y ( y y ) (B5) i i i1 i ai1 ai
Minimum-Time Trajectory Optimization of Multiple Revolution Low-Thrust Earth-Orbit Transfers
Minimum-Time Trajectory Optimization o Multiple Revolution Low-Thrust Earth-Orbit Transers Kathryn F. Graham Anil V. Rao University o Florida Gainesville, FL 32611-625 Abstract The problem o determining
More informationThe Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (014) 000 000 www.elsevier.com/locate/procedia APISAT014, 014 Asia-Paciic International Symposium on Aerospace Technology,
More informationOptimal Control based Time Optimal Low Thrust Orbit Raising
Optimal Control based Time Optimal Low Thrust Orbit Raising Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P., India 2
More informationStrathprints Institutional Repository
Strathprints Institutional Repository Docherty, Stephanie and Macdonald, Malcolm (2012) Analytical sun synchronous low-thrust manoeuvres. Journal of Guidance, Control and Dynamics, 35 (2). pp. 681-686.
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationPRECISION ZEM/ZEV FEEDBACK GUIDANCE ALGORITHM UTILIZING VINTI S ANALYTIC SOLUTION OF PERTURBED KEPLER PROBLEM
AAS 16-345 PRECISION ZEM/ZEV FEEDBACK GUIDANCE ALGORITHM UTILIZING VINTI S ANALYTIC SOLUTION OF PERTURBED KEPLER PROBLEM Jaemyung Ahn, * Yanning Guo, and Bong Wie A new implementation o a zero-eort-miss/zero-eort-velocity
More informationInterplanetary Mission Opportunities
Interplanetary Mission Opportunities Introduction The quest for unravelling the mysteries of the universe is as old as human history. With the advent of new space technologies, exploration of space became
More informationMinimum-Fuel Trajectory Optimization of Many Revolution Low-Thrust Earth-Orbit Transfers
Minimum-Fuel Trajectory Optimization o Many Revolution Low-Thrust Earth-Orbit Transers Kathryn F. Graham 1, Anil V. Rao 2 University o Florida, Gainesville, FL 32611-625 Abstract The problem o determining
More informationAPPLYING FLYWHEEL ENERGY STORAGE TO SOLAR ELECTRIC ORBITAL TRANSFERS THESIS. Presented to the Faculty of the Graduate School of Engineering
AFIT/GA/ENY/97D-02 APPLYING FLYWHEEL ENERGY STORAGE TO SOLAR ELECTRIC ORBITAL TRANSFERS THESIS Presented to the Faculty o the Graduate School o Engineering o the Air Force Institute o Technology Air University
More informationAUTOMATIC MINIMUM PRINCIPLE FORMULATION FOR LOW THRUST OPTIMAL CONTROL IN ORBIT TRANSFERS USING COMPLEX NUMBERS
AUOMAIC MINIMUM PRINCIPLE FORMULAION FOR LOW HRUS OPIMAL CONROL IN ORBI RANSFERS USING COMPLEX NUMBERS hierry Dargent hales Alenia Space, 1 Bd du Midi, BP 99, 6322 Cannes La Bocca Cedex France, E-mail:
More informationOptElec: an Optimisation Software for Low-Thrust Orbit Transfer Including Satellite and Operation Constraints
OptElec: an Optimisation Software for Low-Thrust Orbit Transfer Including Satellite and Operation Constraints 7th International Conference on Astrodynamics Tools and Techniques, DLR, Oberpfaffenhofen Nov
More informationProbabilistic Optimisation applied to Spacecraft Rendezvous on Keplerian Orbits
Probabilistic Optimisation applied to pacecrat Rendezvous on Keplerian Orbits Grégory aive a, Massimiliano Vasile b a Université de Liège, Faculté des ciences Appliquées, Belgium b Dipartimento di Ingegneria
More informationAUTONOMOUS AND ROBUST RENDEZVOUS GUIDANCE ON ELLIPTICAL ORBIT SUBJECT TO J 2 PERTURBATION.
AUTONOMOUS AND ROBUST RENDEZVOUS GUIDANCE ON ELLIPTICAL ORBIT SUBJECT TO J 2 PERTURBATION Emmanuel GOGIBUS (1), Hervé CHARBONNEL (2), Patrick DELPY (3) (1) Astrium Space Transportation, 66 route de Verneuil,
More informationResearch Article Generalized Guidance Scheme for Low-Thrust Orbit Transfer
Mathematical Problems in Engineering, Article ID 4787, 9 pages http://dx.doi.org/1.1155/214/4787 Research Article Generalized Guidance Scheme for Low-Thrust Orbit Transfer Henzeh Leeghim, 1 Dong-Hyun Cho,
More informationMinimum-Time Trajectory Optimization of Low-Thrust Earth-Orbit Transfers with Eclipsing
Minimum-Time Trajectory Optimization of Low-Thrust Earth-Orbit Transfers with Eclipsing Kathryn F. Graham Anil V. Rao University of Florida Gainesville, FL 32611-625 Abstract The problem of determining
More informationCross-Range Error in the Lambert Scheme
Proceedings o the Tenth National Aeronautical Conerence, Edited by Sheikh SR, Khan AM, College o Aeronautical Engineering, PAF Academy, Risalpur, NWFP, Pakistan, April -, 006, pp 55-6 Cross-Range Error
More informationASTOS for Low Thrust Mission Analysis
ASTOS for Low Thrust Mission Analysis 3rd Astrodynamics Workshop, Oct. 26, ESTEC Overview Low Thrust Trajectory Computation Description of the Optimal Control Problem Trajectory Optimization and Mission
More informationPowered Space Flight
Powered Space Flight KOIZUMI Hiroyuki ( 小泉宏之 ) Graduate School of Frontier Sciences, Department of Advanced Energy & Department of Aeronautics and Astronautics ( 基盤科学研究系先端エネルギー工学専攻, 工学系航空宇宙工学専攻兼担 ) Scope
More informationCross-Range Error in the Lambert Scheme
Proceedings o the Tenth National Aeronautical Conerence, dited by Sheikh SR, Khan AM, College o Aeronautical ngineering, PAF Academy, Risalpur, KP, Pakistan, April -, 006, pp 55-6 Cross-Range rror in the
More informationOptimal Guidance Strategy for Low Earth Orbit Based on the Laser Solutions
Optimal Guidance Strategy for Low Earth Orbit Based on the Laser Solutions Ting Jin, Jian Cai Institute of Microelectronics of Chinese Academy of Science, Beijing, China 28 April 2015 Laboratory Research
More informationJames E. Pollard The Aerospace Corporation, P-0. Box 92957, M5-754, Los Angeles, CA 90009
IEPC-97-160 979 Simplified Approach for Assessment of Low-Thrust Elliptical Orbit Transfers James E. Pollard The Aerospace Corporation, P-0. Box 92957, M5-754, Los Angeles, CA 90009 Abstract Low-thrust
More informationFundamentals of Astrodynamics and Applications
Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA
More informationLecture D30 - Orbit Transfers
J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we
More informationA Concept Study of the All-Electric Satellite s Attitude and Orbit Control System in Orbit Raising
Journal of Automation and Control Engineering Vol., No., December A Concept Study of the All-Electric Satellite s Attitude and Orbit Control System in Orbit Raising Yoshinobu Sasaki Japan Aerospace Exploration
More informationEscape Trajectories from Sun Earth Distant Retrograde Orbits
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists30, pp. Pd_67-Pd_75, 06 Escape Trajectories from Sun Earth Distant Retrograde Orbits By Yusue OKI ) and Junichiro KAWAGUCHI ) ) Department of Aeronautics
More informationEasyChair Preprint. Retrograde GEO Orbit Design Method Based on Lunar Gravity Assist for Spacecraft
EasyChair Preprint 577 Retrograde GEO Orbit Design Method Based on Lunar Gravity Assist for Spacecraft Renyong Zhang EasyChair preprints are intended for rapid dissemination of research results and are
More informationEarth-Mars Halo to Halo Low Thrust
Earth-Mars Halo to Halo Low Thrust Manifold Transfers P. Pergola, C. Casaregola, K. Geurts, M. Andrenucci New Trends in Astrodynamics and Applications V 3 June / -2 July, 28 Milan, Italy Outline o Introduction
More informationA SURVEY OF BALLISTIC TRANSFERS TO LOW LUNAR ORBIT
AAS 11-277 A SURVEY OF BALLISTIC TRANSFERS TO LOW LUNAR ORBIT Jeffrey S. Parker Rodney L. Anderson and Andrew Peterson A simple strategy is identified to generate ballistic transfers between the Earth
More informationCOUPLED OPTIMIZATION OF LAUNCHER AND ALL-ELECTRIC SATELLITE TRAJECTORIES
COUPLED OPTIMIZATION OF LAUNCHER AND ALL-ELECTRIC SATELLITE TRAJECTORIES M. Verlet (1), B. Slama (1), S. Reynaud (1), and M. Cerf (1) (1) Airbus Defence and Space, 66 Route de Verneuil, 78133 Les Mureaux,
More informationThe B-Plane Interplanetary Mission Design
The B-Plane Interplanetary Mission Design Collin Bezrouk 2/11/2015 2/11/2015 1 Contents 1. Motivation for B-Plane Targeting 2. Deriving the B-Plane 3. Deriving Targetable B-Plane Elements 4. How to Target
More informationOPTIMIZING PERIAPSIS-RAISE MANEUVERS USING LOW-THRUST PROPULSION
AAS 8-298 OPTIMIZING PERIAPSIS-RAISE MANEUVERS USING LOW-THRUST PROPULSION Brenton J. Duffy and David F. Chichka This study considers the optimal control problem of maximizing the raise in the periapsis
More informationSteered Spacecraft Deployment Using Interspacecraft Coulomb Forces
Steered Spacecrat Deployment Using Interspacecrat Coulomb Forces Gordon G. Parker, Lyon B. King and Hanspeter Schaub Abstract Recent work has shown that Coulomb orces can be used to maintain ixed-shape
More informationLunar Landing Trajectory and Abort Trajectory Integrated Optimization Design.
Lunar Landing Trajectory and Abort Trajectory Integrated Optimization Design Bai Chengchao (1), Guo Jifeng (2), and Xu Xibao (3) (1)(2) School of Astronautics, Harbin Institute of Technology, (451)864128766
More information5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...
1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...
More informationACCURACY ASSESSMENT OF GEOSTATIONARY-EARTH-ORBIT WITH SIMPLIFIED PERTURBATIONS MODELS
ARTIFICIAL SATELLITES, Vol. 51, No. 2 2016 DOI: 10.1515/arsa-2016-0005 ACCURACY ASSESSMENT OF GEOSTATIONARY-EARTH-ORBIT WITH SIMPLIFIED PERTURBATIONS MODELS Lihua Ma, Xiaojun Xu, Feng Pang National Astronomical
More informationAstromechanics. 6. Changing Orbits
Astromechanics 6. Changing Orbits Once an orbit is established in the two body problem, it will remain the same size (semi major axis) and shape (eccentricity) in the original orbit plane. In order to
More informationLow-Thrust Trajectories to the Moon
3rd WSEAS International Conference on APPLIED and THEORETICAL MECHANICS, Spain, December 14-16, 7 257 Low-Thrust Trajectories to the Moon ANTONIO F. B. A. PRADO Space Mechanics and Control Division INPE
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationExpanding opportunities for lunar gravity capture
Expanding opportunities for lunar gravity capture Keita Tanaka 1, Mutsuko Morimoto 2, Michihiro Matsumoto 1, Junichiro Kawaguchi 3, 1 The University of Tokyo, Japan, 2 JSPEC/JAXA, Japan, 3 ISAS/JAXA, Japan,
More informationA Comparison of Low Cost Transfer Orbits from GEO to LLO for a Lunar CubeSat Mission
A Comparison of Low Cost Transfer Orbits from GEO to LLO for a Lunar CubeSat Mission A presentation for the New Trends in Astrodynamics conference Michael Reardon 1, Jun Yu 2, and Carl Brandon 3 1 PhD
More informationRigorous Global Optimization of Impulsive Space Trajectories
Rigorous Global Optimization of Impulsive Space Trajectories P. Di Lizia, R. Armellin, M. Lavagna K. Makino, M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16 19, 2006 Motivation
More informationOptimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming
Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationLaunch strategy for Indian lunar mission and precision injection to the Moon using genetic algorithm
Launch strategy for Indian lunar mission and precision injection to the Moon using genetic algorithm VAdimurthy, R V Ramanan, S R Tandon and C Ravikumar Aeronautics Entity, Vikram Sarabhai Space Centre,
More information1 The Problem of Spacecraft Trajectory Optimization
1 The Problem of Spacecraft Trajectory Optimization Bruce A. Conway Dept. of Aerospace Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 1.1 Introduction The subject of spacecraft trajectory
More informationTHIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL
INPE-1183-PRE/67 THIRD-BODY PERTURBATION USING A SINGLE AVERAGED MODEL Carlos Renato Huaura Solórzano Antonio Fernando Bertachini de Almeida Prado ADVANCES IN SPACE DYNAMICS : CELESTIAL MECHANICS AND ASTRONAUTICS,
More informationASTRIUM. Interplanetary Path Early Design Tools at ASTRIUM Space Transportation. Nathalie DELATTRE ASTRIUM Space Transportation.
Interplanetary Path Early Design Tools at Space Transportation Nathalie DELATTRE Space Transportation Page 1 Interplanetary missions Prime approach: -ST has developed tools for all phases Launch from Earth
More informationMAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30.
MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details that
More informationMission Scenarios for a Controlled Lunar Impact of a Small Satellite
IAC-4-IAA.4.11.P.5 Mission Scenarios for a Controlled Lunar Impact of a Small Satellite Nikolas Trawny, Michael Graesslin, Rene Laufer and Hans-Peter Roeser Email: n.trawny@gmx.de, {graesslin,laufer,roeser}@irs.uni-stuttgart.de
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem
More informationA Study of the Close Approach Between a Planet and a Cloud of Particles
A Study of the Close Approach Between a Planet a Cloud of Particles IIAN MARTINS GOMES, ANTONIO F. B. A. PRADO National Institute for Space Research INPE - DMC Av. Dos Astronautas 1758 São José dos Campos
More informationPLANETARY MISSIONS FROM GTO USING EARTH AND MOON GRAVITY ASSISTS*
. AIAA-98-4393 PLANETARY MISSIONS FROM GTO USING EARTH AND MOON GRAVITY ASSISTS* Paul A. Penzo, Associate Fellow AIAA+ Jet Propulsion Laboratory California Institute of Technology 4800 Oak Grove Dr. Pasadena,
More informationResearch Article Low-Thrust Transfer Design of Low-Observable Geostationary Earth Orbit Satellite
International Journal o Aerospace Engineering Volume 215, Article ID 439815, 12 pages http://dx.doi.org/1.1155/215/439815 Research Article Low-Thrust Transer Design o Low-Observable Geostationary Earth
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) L06: Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 Problem Statement? Hint #1: design the Earth-Mars transfer using known concepts
More informationASEN 6008: Interplanetary Mission Design Lab Spring, 2015
ASEN 6008: Interplanetary Mission Design Lab Spring, 2015 Lab 4: Targeting Mars using the B-Plane Name: I d like to give credit to Scott Mitchell who developed this lab exercise. He is the lead Astrodynamicist
More informationINDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT ABOUT LIBRATION POINTS
6 th International Conference on Astrodynamics Tools and Technique (ICATT) INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT LI Xiangyu, Qiao Dong, Cui Pingyuan Beijing Institute of Technology Institute of
More information1. (a) Describe the difference between over-expanded, under-expanded and ideallyexpanded
Code No: R05322106 Set No. 1 1. (a) Describe the difference between over-expanded, under-expanded and ideallyexpanded rocket nozzles. (b) While on its way into orbit a space shuttle with an initial mass
More informationGRAIL Takes A Roundabout Route To Lunar Orbit
The Gravity Recovery and Interior Laboratory (GRAIL) mission departed Earth from Cape Canaveral Air Force Station Space Launch Complex 17B on 2011 September 10 at 13:08 UTC, carried by a Delta II Heavy
More informationA Simple Semi-Analytic Model for Optimum Specific Impulse Interplanetary Low Thrust Trajectories
A Simple Semi-Analytic Model for Optimum Specific Impulse Interplanetary Low Thrust Trajectories IEPC-2011-010 * Presented at the 32nd International Electric Propulsion Conference, Wiesbaden Germany David
More informationRELATIVISTIC EFFECTS O N SV CLOCKS DUE TO ORBIT CHANGES, AND DUE TO EARTH S OBLATENESS
33rdAnnual Precise Time and Time Interval (PTTI)Meeting RELATIVISTIC EFFECTS O N SV CLOCKS DUE TO ORBIT CHANGES, AND DUE TO EARTH S OBLATENESS Neil Ashby Department o Physics, UCB 390 University o Colorado,
More informationGlobal Optimization of Impulsive Interplanetary Transfers
Global Optimization of Impulsive Interplanetary Transfers R. Armellin, Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano Taylor Methods and Computer Assisted Proofs Barcelona, June, 3 7, 2008
More informationOptimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission
Optimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P.,
More information: low-thrust transfer software, optimal control problem, averaging techniques.
J. Fourcade S. Geffroy R.Epenoy Centre National d Etudes Spatiales 8 avenue Edouard Belin 4 Toulouse cedex 4 France e-mail : Jean.Fourcade@cnes.fr Sophie.Geffroy@cnes.fr Richard.Epenoy@cnes.fr Low thrust
More informationChapter 4. Integrated Algorithm for Impact Lunar Transfer Trajectories. using Pseudo state Technique
Chapter 4 Integrated Algorithm for Impact Lunar Transfer Trajectories using Pseudo state Technique Abstract A new integrated algorithm to generate the design of one-way direct transfer trajectory for moon
More informationAUGMENTED POLYNOMIAL GUIDANCE FOR TERMINAL VELOCITY CONSTRAINTS
AUGMENTED POLYNOMIAL GUIDANCE FOR TERMINAL VELOCITY CONSTRAINTS Gun-Hee Moon*, Sang-Woo Shim*, and Min-Jea Tah* *Korea Advanced Institute Science and Technology Keywords: Polynomial guidance, terminal
More informationFundamentals of Satellite technology
Fundamentals of Satellite technology Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Orbital Plane All of the planets,
More informationTRAJECTORY SIMULATIONS FOR THRUST-VECTORED ELECTRIC PROPULSION MISSIONS
RAJECORY SIMULAIONS FOR HRUS-VECORED ELECRIC PROPULSION MISSIONS Abstract N. Leveque, C. Welch, A. Ellery, A. Curley Kingston University, Astronautics and Space Systems Group School of Engineering Friars
More informationFeedback Optimal Control for Inverted Pendulum Problem by Using the Generating Function Technique
(IJACSA) International Journal o Advanced Computer Science Applications Vol. 5 No. 11 14 Feedback Optimal Control or Inverted Pendulum Problem b Using the Generating Function echnique Han R. Dwidar Astronom
More informationOrbital Transfer Trajectory Optimization. James K Whiting
Orbital Transfer Trajectory Optimization by James K Whiting Submitted to the Department of Aeronautical and Astronautical Engineering in partial fulfillment of the requirements for the degree of Master
More informationORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM
ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM JORGE K. S. FORMIGA 1,2 and ANTONIO F B A PRADO 2 National Institute for Space Research -INPE 1 Technology Faculty-FATEC-SJC
More informationList of Tables. Table 3.1 Determination efficiency for circular orbits - Sample problem 1 41
List of Tables Table 3.1 Determination efficiency for circular orbits - Sample problem 1 41 Table 3.2 Determination efficiency for elliptical orbits Sample problem 2 42 Table 3.3 Determination efficiency
More informationEscape Trajectories from the L 2 Point of the Earth-Moon System
Trans. Japan Soc. Aero. Space Sci. Vol. 57, No. 4, pp. 238 244, 24 Escape Trajectories from the L 2 Point of the Earth-Moon System By Keita TANAKA Þ and Jun ichiro KAWAGUCHI 2Þ Þ Department of Aeronautics
More informationLOW-COST LUNAR COMMUNICATION AND NAVIGATION
LOW-COST LUNAR COMMUNICATION AND NAVIGATION Keric Hill, Jeffrey Parker, George H. Born, and Martin W. Lo Introduction Spacecraft in halo orbits near the Moon could relay communications for lunar missions
More informationVerified High-Order Optimal Control in Space Flight Dynamics
Verified High-Order Optimal Control in Space Flight Dynamics R. Armellin, P. Di Lizia, F. Bernelli-Zazzera K. Makino and M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16
More informationTrajectory Optimization for Ascent and Glide Phases Using Gauss Pseudospectral Method
Trajectory Optimization for Ascent and Glide Phases Using Gauss Pseudospectral Method Abdel Mageed Mahmoud, Chen Wanchun, Zhou Hao, and Liang Yang Abstract The trajectory optimization method for ascent
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More informationAn Architecture for a Generalized Spacecraft Trajectory Design and Optimization System
An Architecture for a Generalized Spacecraft Trajectory Design and Optimization System Cesar Ocampo Assistant Professor, Department of Aerospace Engineering and Engineering Mechanics, The University of
More informationChapter 5 - Part 1. Orbit Perturbations. D.Mortari - AERO-423
Chapter 5 - Part 1 Orbit Perturbations D.Mortari - AERO-43 Orbital Elements Orbit normal i North Orbit plane Equatorial plane ϕ P O ω Ω i Vernal equinox Ascending node D. Mortari - AERO-43 Introduction
More informationAn Analysis of N-Body Trajectory Propagation. Senior Project. In Partial Fulfillment. of the Requirements for the Degree
An Analysis of N-Body Trajectory Propagation Senior Project In Partial Fulfillment of the Requirements for the Degree Bachelor of Science in Aerospace Engineering by Emerson Frees June, 2011 An Analysis
More informationOrbital Mechanics MARYLAND. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design
Lecture #08 September 22, 2016 Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit
More informationGravity Assisted Maneuvers for Asteroids using Solar Electric Propulsion
Gravity Assisted Maneuvers for Asteroids using Solar Electric Propulsion Denilson P. S. dos Santos, Antônio F. Bertachini de A. Prado, Division of Space Mechanics and Control INPE C.P. 515, 17-310 São
More informationSpace Travel on a Shoestring: CubeSat Beyond LEO
Space Travel on a Shoestring: CubeSat Beyond LEO Massimiliano Vasile, Willem van der Weg, Marilena Di Carlo Department of Mechanical and Aerospace Engineering University of Strathclyde, Glasgow 5th Interplanetary
More informationEffect of Coordinate Switching on Translunar Trajectory Simulation Accuracy
Effect of Coordinate Switching on Translunar Trajectory Simulation Accuracy Mana P. Vautier Auburn University, Auburn, AL, 36849, USA This paper focuses on the affect of round-off error in the accurate
More informationIMPROVED DESIGN OF ON-ORBIT SEPARATION SCHEMES FOR FORMATION INITIALIZATION BASED ON J 2 PERTURBATION
IAA-AAS-DyCoSS- IMPROVED DESIGN OF ON-ORBIT SEPARATION SCHEMES FOR FORMATION INITIALIZATION BASED ON J PERTURBATION Jiang Chao, * Wang Zhaokui, and Zhang Yulin INTRODUCTION Using one Multi-satellite Deployment
More informationL eaving Earth and arriving at another planet or asteroid requires
Designing Interplanetary Transfers L eaving Earth and arriving at another planet or asteroid requires a spacecraft to implement a sequence of manoeuvres. These include changes of velocity needed to escape
More informationDEVELOPMENT OF A MULTIPURPOSE LOW THRUST INTERPLANETARY TRAJECTORY CALCULATION CODE
AAS 03-667 DEVELOPMENT OF A MULTIPURPOSE LOW THRUST INTERPLANETARY TRAJECTORY CALCULATION CODE Tadashi Sakai 1 and John R. Olds 2 A multipurpose low thrust interplanetary trajectory calculation code has
More informationOrbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design
Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit Interplanetary trajectories
More informationSUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS. Sandro da Silva Fernandes. Cleverson Maranhão Porto Marinho
SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS Sandro da Silva Fernandes Instituto Tecnológico de Aeronáutica, São José dos Campos - 12228-900 - SP-Brazil, (+55) (12) 3947-5953 sandro@ita.br Cleverson
More informationSatellite meteorology
GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading
More informationOrbital Mechanics. John E. Prussing & Bruce A.Conway. Published by Oxford University Press 2013
Orbital Mechanics by John E. Prussing & Bruce A.Conway Published by Oxford University Press 2013 Copyright 2013, 1993 by Oxford University Press ISBN: 978-0-19-983770-0 Chapter 8 Continuous-Thrust Orbit
More informationFlight and Orbital Mechanics
Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and Orbital Mechanics AE-104, lecture hours 1-4: Interplanetary flight Ron Noomen October 5, 01 AE104 Flight and Orbital Mechanics
More informationESMO Mission Analysis
Changing the economics of space ESMO Mission Analysis SRR Workshop Alison Gibbings 22 nd 26 th March 2010 Review of the existing baseline Sensitivity analysis Contents At lunar Injection Along the WSB-Moon
More informationRAPID GEOSYNCHRONOUS TRANSFER ORBIT ASCENT PLAN GENERATION. Daniel X. Junker (1) Phone: ,
RAPID GEOSYNCHRONOUS TRANSFER ORBIT ASCENT PLAN GENERATION Daniel X. Junker (1) (1) LSE Space GmbH, Argelsrieder Feld 22, 82234 Wessling, Germany, Phone: +49 160 9111 6696, daniel.junker@lsespace.com Abstract:
More informationAn adaptive model predictive controller for turbofan engines
American Journal o Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-4, Issue-12, pp-170-176 www.ajer.org Research Paper Open Access An adaptive model predictive controller or turboan
More informationProblem Description for the 8 th Global Trajectory Optimisation Competition. Anastassios E. Petropoulos 1
1 Background Problem Description for the 8 th Global Trajectory Optimisation Competition Anastassios E. Petropoulos 1 Outer Planet Mission Analysis Group Mission Design and Navigation Section Jet Propulsion
More informationEarth-Centered, Earth-Fixed Coordinate System
Fundamentals of Global Positioning System Receivers: A Software Approach James Bao-Yen Tsui Copyright 2000 John Wiley & Sons, Inc. Print ISBN 0-471-38154-3 Electronic ISBN 0-471-20054-9 CHAPTER FOUR Earth-Centered,
More informationIdentifying Safe Zones for Planetary Satellite Orbiters
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 16-19 August 2004, Providence, Rhode Island AIAA 2004-4862 Identifying Safe Zones for Planetary Satellite Orbiters M.E. Paskowitz and D.J. Scheeres
More informationOrbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design
Lecture #05 September 15, 2015 Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit
More informationExperimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations
Experimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations Aman Saluja #1, Manish Bansal #2, M Raja #3, Mohd Maaz #4 #Aerospace Department, University of Petroleum and Energy
More informationSession 6: Analytical Approximations for Low Thrust Maneuvers
Session 6: Analytical Approximations for Low Thrust Maneuvers As mentioned in the previous lecture, solving non-keplerian problems in general requires the use of perturbation methods and many are only
More information