A Sliding Mode Control-Based Guidance Law for a Two-Dimensional Orbit Transfer with Bounded Disturbances

: The aim of this paper is to analyze the performance of a state-feedback guidance law, which is obtained through a classical sliding mode control approach, in a two-dimensional circle-to-circle orbit transfer of a spacecraft equipped with a continuous-thrust propulsion system. The paper shows that such an inherently robust control technique can be effectively used to obtain possible transfer trajectories even when the spacecraft equations of motion are affected by perturbations. The problem of the guidance law design is ﬁrst addressed in the simpliﬁed case of an unperturbed system, where it is shown how the state-feedback control may be effectively used to obtain simple mathematical relationships and graphs that allow the designer to determine possible transfer trajectories that depend on a few control parameters. It is also shown that a suitable combination of the controller parameters may be exploited to obtain trade-off solutions between the ﬂight time and the transfer velocity change. The simpliﬁed control strategy is then used to investigate a typical heliocentric orbit raising/lowering in the presence of bounded disturbances and measurement errors.


Introduction
In a preliminary phase of mission design, the use of a state-feedback control law represents a viable option to obtain possible transfer trajectories that may be used as an initial starting point for succeeding (and more refined) analyses.In this context, an interesting approach is based on the use of a rather classical sliding mode control, which is a variable structure (control) method that alters the dynamical behaviour of a nonlinear system through the application of a suitable control signal [1].In particular, sliding mode control is a basic robust technique, which allows the system trajectory to converge towards the desired target even in the presence of significant perturbations and measurement errors [2,3].An interesting discussion about the potentialities of a sliding mode control law can be found in the review paper by Hung et al. [4] and in the work by Utkin [5].
The literature about the control of satellites by means of sliding mode techniques is rich, although essentially concentrated on attitude control and terminal guidance maneuvers [6].In this scenario, Wu et al. [7] investigated the attitude synchronization and tracking problem including model uncertainties, external disturbances, actuator failures, and control torque saturation.By proposing two decentralized sliding mode control laws, that work [7] proved that the control laws guarantee each spacecraft to approach the desired timevarying attitude and angular velocity while maintaining attitude synchronization among the other elements in a typical formation structure.Another example is offered by the work by Massey and Shtessel [8], who adopted a traditional, continuous, high-order sliding mode strategy to control a satellite formation in a robust manner (i.e., compensating for model uncertainties and external disturbances).An adaptive sliding mode tension control method was successfully proposed by Ma et al. [9] for the deployment of tethered satellites, when input tension limitations are taken into account.On the other hand, a terminal sliding mode control law was designed by Liu and Huo [10] for spacecraft rendezvous and docking while considering both model uncertainties and external disturbances, proving that the closed-loop tracking error converges to zero in a finite time.The same problem was also successfully addressed by Dong et al. [11], who constructed a nonsingular terminal sliding surface by introducing a continuous sinusoidal function to solve the inherent singularity problem.More recently, Capello et al. [12] designed two controllers, that is, a first-order sliding mode control for position tracking and a supertwisting second-order sliding mode control for attitude stability, in which the mutual influence was taken into account by the introduction of additional disturbances.Kasaeian et al. [13] presented a robust guidance algorithm to perform a rendezvous between a chaser and a target spacecraft orbiting around the Earth, revealing that sliding mode control guarantees the tracking of the required states and minimum final errors even in the presence of uncertainties and disturbances.Li et al. [14] developed a novel sliding mode control strategy to address the relative position tracking and attitude synchronization problem of spacecraft rendezvous with the requirement of collision avoidance, proving the convergence of relative position and attitude errors even in the presence of external disturbances.Finally, Bassetto et al. [15] discussed how solar sail attitude maneuvers may be designed in a collinear, artificial, equilibrium point by implementing a sliding mode control strategy that uses electrochromic devices as actuators [16][17][18].Anyway, there are many other potential feedback control techniques [19][20][21], to which the interest reader is invited to refer.
In the context of spacecraft trajectory design, a robust state-feedback control law can be used to obtain a possible transfer trajectory that is useful as an initial guess during the subsequent refinement phase [22].In that case, potential transfer trajectories can be obtained by taking into account the orbit perturbations and the model uncertainties with a reduced computational cost [23].The aim of this paper is to investigate the potentialities of a sliding mode control strategy in detecting possible trajectories in a typical circle-to-circle orbit transfer scenario, in which the spacecraft propulsion system provides a continuously adjustable and freely steerable propulsive acceleration vector.Among actuators capable of generating variable propulsive acceleration, there are variable thrust ion engines (such as NASA's Evolutionary Xenon Thruster (NEXT)), in which continuous thrust variation can be replaced by a succession of discrete thrust levels that, on average, provide the required propulsive acceleration.For example, NEXT has a total of 40 operating points, with available thrust ranging from 25.5 mN and 236 mN [24].The proposed approach uses a standard implementation of the sliding mode procedure [1] to obtain a set of preliminary results.In this way, the discussed procedure allows the designer to make a trade-off between the flight time and the required velocity change by selecting the design parameters of the controller.In particular, the simplified control strategy involves three independent parameters of the spacecraft dynamics (on which the resulting propulsive acceleration profile and the characteristics of the transfer trajectory depend), which represent tuning quantities to be selected by the designer.The main limitation of the proposed approach lies in the use of an ideal propulsion system to control the nonlinear dynamics of the spacecraft center of mass.In fact, the time-variation of the thrust vector magnitude, which is an output of the design procedure, can be used a posteriori to check whether the obtained transfer trajectory is compatible with the physical constraints of the thruster, such as the maximum thrust level.
Starting from the simplified scenario, in which the spacecraft orbital dynamics is unaffected by external disturbances or model uncertainties, we firstly discuss how the controller may be tuned by considering the flight time and total velocity change.The procedure is then used to investigate an orbit raising/lowering in the presence of bounded disturbances and measurement errors.The paper is organized as follows.Section 2 presents the mathematical model, i.e., the nonlinear differential equations describing the coplanar orbital motion of a spacecraft around an assigned primary body.Section 3 introduces the sliding mode control technique in its general form, where bounded disturbances are included in the model.Section 4 addresses the design of the sliding mode control law in the simplified case of an unperturbed system.In particular, Section 4 illustrates the time-variation of tracking errors (which can be analytically determined when no disturbance is considered in the mathematical model) and the definition of the control law parameters.The numerical simulations are described in Section 5, while the concluding remarks are drawn in Section 6.

Problem Description and Mathematical Model
Consider a spacecraft S that covers a circular parking orbit of radius r 0 around a primary body with center of mass P and gravitational parameter µ.The mission purpose is to transfer the spacecraft to a circular and coplanar target orbit of assigned radius r f = r 0 by means of a continuously adjustable (and freely steerable) propulsion system, which gives both a radial (a r ) and a transverse (a t ) component of propulsive acceleration.In this context, the spacecraft two-dimensional dynamics may be described by the classical polar equations of motion [25]: where r is the P-S distance and θ is the spacecraft polar angle measured counterclockwise from the P-S line at the initial time t 0 0, while v r (or v t ) is the radial (or transverse) component of the spacecraft velocity vector; see Figure 1.In Equations ( 3) and (4), the terms {d r , d t } represent possible unknown bounded disturbance accelerations acting along the radial and transverse directions, with where D r ≥ 0 and D t ≥ 0 are two constant parameters.Bearing in mind that the parking orbit is circular, Equations ( 1)-( 4) are completed by the initial conditions while the time-variation of the propulsive acceleration components {a r , a t } are to be found so as to bring the spacecraft states to the desired final values within a time interval t f .Note that the final polar angle θ(t f ), that is, the angle swept out by the spacecraft during the transfer, is left free.The spacecraft dynamics is more conveniently rewritten by introducing the dimensionless tracking errors in radial distance (x 1 ), radial velocity component (x 2 ), and transverse velocity component (x 3 ), defined as where ρ r f /r 0 = 1 is the dimensionless radius of the target circular orbit.In particular, ρ ∈ (0, 1) in case of an orbit lowering, whereas ρ > 1 in case of an orbit raising.Substituting Equations ( 8)- (10) into Equations ( 1), (3), and (4) yields where the prime symbol denotes a derivative taken with respect to the dimensionless time τ, defined as are the two dimensionless control variables, defined as the ratio of the propulsive acceleration components {a r , a t } to the primary body gravitational acceleration at r = r 0 .Finally, the two terms {z r , z t } in Equations ( 12) and ( 13) are the dimensionless forms of the disturbance acceleration components, defined as which, by assumption, satisfy the inequalities Equations ( 11)-( 13) are integrated with the three initial conditions from which it follows that x 1 0 < 0 and x 3 0 > 0 when ρ > 1, while x 1 0 > 0 and x 3 0 < 0 when ρ ∈ (0, 1).Finally, the target states (i.e., the conditions on the target circular orbit) are expressed in a dimensionless form as where is the dimensionless flight time.Note that the τ-variation of θ can be obtained by numerically integrating the differential equation which is not included in the dynamical system because the final polar angle is left free.However, solving Equation ( 21) is necessary to obtain the polar trajectory of the spacecraft.

State-Feedback Control Design
In this section, a classical sliding mode control law is used to determine the circle-tocircle orbit transfer trajectory.The spacecraft states are brought and maintained on two sliding surfaces, where the system exhibits the desired dynamics of reduced order or one of the states is at its final equilibrium point.More specifically, the first sliding surface is described by the equation where λ > 0 is a dimensionless design parameter, so that the states {x 1 , x 2 } exhibit a first-order dynamics when the system is on that sliding surface.In fact, bearing in mind Equation ( 11), the condition s = 0 implies from which the τ-variations of the tracking errors {x 1 , x 2 } turn out to be proportional to e −λ τ , viz.
In other terms, when the system is brought to the sliding surface s = 0, both x 1 and x 2 converge exponentially to zero with a convergence rate equal to λ.Now, in order to bring the system on the sliding surface s = 0, it is required that s < 0 when s > 0, and s > 0 when s < 0. To this end, differentiating s with respect to τ yields from which selecting u r according to the law where sign( ) is the signum function and δ is given by in which K > 0 is a dimensionless design parameter, one obtains In this case, s < 0 implies s = z r + Z r + K > 0 and s > 0 implies s = z r − Z r − K < 0, while s = 0 implies s = z r , that is, the perturbative term z r forces the system to leave the sliding surface s = 0 once it has been reached.
The second sliding surface is the plane x 3 = 0; see Equation (10).Note that the system is driven to the sliding surface x 3 = 0 if x 3 < 0 when x 3 > 0, and if x 3 > 0 when x 3 < 0. In this context, if the control parameter u t is selected as where c > 0 is a dimensionless design parameter, then In this case, x 3 < 0 implies that is, the perturbative term z t moves the system away from the sliding surface x 3 = 0 once it has been reached.

Disturbance Modeling
This section gives a brief description of the source of disturbances (or uncertainties) that will be included in the numerical simulations.The first one is related to the state measurement.In fact, when applying a state-feedback control law, it is necessary to verify whether and how measurement errors or low-frequency sampling affect the control effectiveness.The measured states, denoted as { x1 , x2 , x3 }, are the sum of a true value x i and a measurement error X i , that is, where X i is a zero-mean random variable with normal distribution and standard deviation σ i .
The second source of disturbance is related to the approximation of the signum function in Equations ( 26) and ( 29) with a sigmoid-like function.Note, in fact, that the change in sign of u r (or u t ) each time the system crosses the sliding surface s = 0 (or x 3 = 0) gives rise to a chattering behaviour, which is typical of the sliding mode control.Such a phenomenon must be mitigated to prevent the switching frequency of the control signals from being too high and, therefore, not applicable.A viable option is to implement a pseudo-sliding mode control [26], which consists of smoothing the discontinuity in the signum function to obtain an arbitrarily close but continuous approximation.One possibility is to approximate the signum function with the sigmoid-like function [26] S = S(x) where κ is an arbitrarily small positive scalar.Note that S(x) → sign(x) as κ → 0. Using such a pseudo-sliding mode control, however, causes sliding to no longer take place because the (continuous) control only drives the states to a neighbourhood of the switching surfaces [26].Accordingly, introducing the measured states in Equations ( 26) and ( 29) and substituting sign(x) with S(x), the control variables become where s x2 + λ x1 .Measurement errors, low-frequency sampling, and the approximation of the signum function with the sigmoid-like function S are all treated as disturbance sources.In practice, this situation is equivalent to using ideal sensors and actuators (i.e., sensors capable of measuring the actual states with continuity and actuators capable of adjusting their control signals with continuity) and to perturbing the system with the following (bounded) disturbance accelerations:

Case of an Unperturbed System
The control law described by Equations ( 26) and (29) takes a simpler form in the case of an unperturbed system, which allows some useful analytical relationships to be found in such a simplified mission scenario.Accordingly, in this section, we analyze the evolution of the tracking errors and address the control law design problem with the significant assumption that Z r = Z t = 0.In this simplified case, Equations ( 28) and (31) become This means that in the absence of perturbative terms, the value of s (or x 3 ) approaches zero linearly with respect to τ, and once the sliding surface s = 0 (or x 3 = 0) is reached for the first time, the term s (or x 3 ) remains stationary at zero.When Equations (40) and (41) are integrated with respect to the dimensionless time τ, one obtains the τ-variations of s and x 3 before reaching the sliding surfaces s = 0 and x 3 = 0, respectively.The result is where s 0 s(τ 0 ) can be written, according to Equations ( 18) and ( 22), as while x 3 0 is given by the last of Equation ( 18) as a function of ρ.Note that {K, c} represent a sort of approach speed to the two sliding surfaces.The value of τ at which the system reaches the sliding surface s = 0 (i.e., τ = τ s ) or the sliding surface x 3 = 0 (i.e., τ = τ x 3 ) can be expressed in a compact form using Equations ( 42) and (43).In fact, enforcing the condition s = 0 in Equation (42) gives while the condition x 3 = 0 in Equation (43) gives from which The value of τ x 3 may be written as a function of τ s in a more convenient way by introducing the dimensionless parameter β > 0 such that Observing that β is a redundant parameter, it may be used in place of c, which can be expressed as a function of {K, λ, β, ρ} as (50)

The τ-Variation of Tracking Errors and Controls
The τ-variation of the tracking errors is now calculated, thus allowing the expressions of u r and u t to be determined through Equations ( 26) and (29) by simply setting δ = K and γ = c.To that end, the differential equation governing the τ-evolution of x 1 is found by substituting Equations ( 11) and (22) into Equation (42) and bearing in mind Equations ( 44)-(46), viz.
Integrating Equation (51) with respect to τ with the initial condition x 1 (τ 0 ) = 1 − ρ (see Equation ( 18)) gives the τ-variation of x 1 , that is, The τ-variation of x 2 is instead obtained by deriving Equation (52) with respect to τ (see Equation ( 11)), that is, Finally, the τ-variation of x 3 is governed by the differential equation which must be solved recalling the initial condition x 3 (τ 0 ) = 1 − 1/ √ ρ (see Equation ( 18)), and the result is Note that in the absence of perturbative terms, the maximum values of |x 1 | and |x 3 | occur when τ = τ 0 , that is, while the maximum value of |x 2 | (which corresponds to the maximum of |v r |) is reached when τ = τ s , viz.
The dimensionless propulsive acceleration components {u r , u t } in absence of perturbative terms are simply obtained by substituting Equations ( 52), (53), and (55) into Equations ( 26) and (29) and setting δ = K and γ = c.Those expressions, which are here omitted for the sake of conciseness, change according to whether β < 1, β = 1, or β > 1.In particular, u r exhibits a discontinuity equal to K sign(1 − ρ) when τ = τ s , whereas u t exhibits a discontinuity equal to c sign 1 − √ ρ when τ = β τ s .Accordingly, if β = 1, the profile of the magnitude u u 2 r + u 2 t has two discontinuities (one at τ = τ s , the other at τ = β τ s ).Otherwise (i.e., when β = 1), the profile of u presents a single discontinuity at τ = τ s .

Control Parameter Selection
For a given value of ρ, the design of the sliding mode control law amounts to selecting the values of the triplet {K, λ, β}.The previous expressions allow the flight time and the total velocity change to be determined and the dimensionless parameters in the control law to be established, according to arbitrary criteria.More precisely, when the spacecraft orbital dynamics is unaffected by external disturbances or model uncertainties, the dimensionless flight time τ f and the total velocity change ∆v can be calculated with analytical expressions or graphic plots that only depend on the design parameters {K, λ, β}.
For example, the value of τ f can be obtained by assuming that the orbit transfer terminates when the tracking errors x 1 and x 2 are sufficiently close to zero.To that end, the value of τ f is defined as the instant at which the exponent λ (τ − τ s ) in Equations ( 52) and ( 53) satisfies the equality λ (τ f − τ s ) = n, for an assigned value of n ∈ R + .In this context, using Equation ( 46), one obtains Note that τ f can be minimized with respect to λ by enforcing the necessary condition in Equation ( 59), from which so that, by assuming λ = λ , the expression of the dimensionless flight time becomes A suitable value of n may be chosen by evaluating the tracking error x 1 at the final time τ = τ f , that is, Figure 2, which describes the variation of x 1 (τ f )/x 1 0 with n when λ = λ , shows that a value of n = 4 (when x 1 (τ f )/x 1 0 0.0045) is reasonable from a practical point of view.In fact, the percentage error in orbital radius, that is, the function is less than 1% when n = 4 and ρ > 0.310; see Figure 3. Therefore, it is assumed that n = 4 in the rest of the paper.According to Equations ( 61) and ( 62), the expressions of λ and τ f when n = 4 become Note that Equation (66) relates the flight time τ f to the value of λ necessary to minimize the flight time for fixed values of K and ρ.Such a value of λ is a function of K and ρ, as described by Equation (65).This means that for given values of ρ and τ f , λ can be chosen by reversing Equation (66), that is, by setting λ = λ ≡ 8/τ f .In this case, the value of K is related to ρ and λ (or to ρ and τ f ) through Equation (65), and Figure 4 shows the variation in τ f with K and ρ.The value of β can be chosen with the aid of another parameter that usually determines the transfer performance.More precisely, β may be related to the total velocity change ∆v of the transfer, defined as Note that β ∈ (0, 2], since β → 0 (or β = 2) means that the sliding surface x 3 = 0 is reached at the beginning (or at the end) of the transfer; see Equations ( 49) and (62).Figure 5 shows the values of β (referred to as β ) that minimize the total velocity change when λ = λ , n = 4, and ρ = {0.723,1.524} (the same values of ρ that will be used for some numerical applications of the proposed control law) as a function of K ∈ (0, 1].In fact, K corresponds to the magnitude of the discontinuity of u r when τ = τ s (see Section 4.1), and a value of K greater than 1 would imply a discontinuity of |a r | greater than the gravitational acceleration on the parking orbit.Figure 6, instead, shows the variation in ∆v with K when λ = λ , n = 4, and β = β .Note that the function ∆v| β=β (K) exhibits a global minimum.When ρ = 0.723, such a minimum is reached when K 0.097 and the corresponding values of ∆v and β are ∆v 0.357 and β 1.368, respectively.If, instead, ρ = 1.524, such a minimum is reached when K 0.032 and the corresponding values of ∆v and β are ∆v 0.324 and β 1.242, respectively.In essence, the design of the control law only requires the choice of the single parameter K.For example, we have already seen that a value of K may be determined using Equations ( 65) and (66) by fixing the total flight time τ f , or using Figure 6 by looking for the value of K that minimizes the function ∆v| β=β (K), so that K may be thought of as sort of trade-off parameter, as discussed in the next section.

Numerical Simulations and Mission Application
The proposed control strategy is now used to analyze two classical circle-to-circle interplanetary transfers.In particular, the radius of the circular parking orbit is r 0 = r ⊕ 1 au, which is consistent with a spacecraft that leaves the Earth's sphere of influence using a parabolic escape trajectory, with the simplifying assumption that the Earth's heliocentric orbit is circular.The radii of the target orbits are r f = {0.723,1.524} au, so the analyzed mission scenarios describe simplified ephemeris-free Earth-Venus and Earth-Mars orbit transfers.
Bearing in mind Equations ( 32)-(34), it is assumed that (i) the sensors measure the states once per day; (ii) σ 1 = σ 2 = σ 3 = 10 −4 , which means that the measurement error in {x 1 , x 2 , x 3 } is less than 0.01% with a probability of 68.3%; and (iii) the sigmoid-like function S = S(x) of Equation ( 35) is obtained with κ = 10 −2 .Although the numerical simulations consider measurement errors, low-frequency sampling, and the approximation of the signum function with the sigmoid-like function, the parameters used in the control law can be those found in Section 4.2 thanks to the robustness of the proposed approach.
For example, assume that λ = λ , n = 4, β = β , and select K such that ∆v| β=β is minimized (the corresponding value of K will be referred to as K v ), so that according to Figure 6, one has K v 0.0969 (or K v 0.0320) when ρ = 0.723 (or ρ = 1.524).The numerical simulations give a flight time of about 394 days (or 949 days) in the Earth-Venus (or Earth-Mars) mission scenario.Moreover, Figure 7 shows the corresponding (twodimensional) heliocentric trajectories, while Figure 8 shows the time-variations of the propulsive acceleration components {a r , a t }.In particular, each black dot in Figure 8 corresponds to one day (i.e., to the sampling period of the states), while the red lines show the propulsive acceleration components in case of ideal sensors and actuators (that is, when only the approximation of the signum function with the sigmoid-like function is taken into account).A second case considered in the simulations is when the flight time is assigned, λ = λ ≡ 8/τ f , n = 4, and β = β .For example, by assuming that the transfer time coincides with the Hohmann transfer one, the values of τ f become while the values of β are chosen by using Figure 5 to minimize the total velocity change, viz.
In this context, Figure 9 shows the interplanetary transfer trajectories, while Figure 10 collects the time-variations of the propulsive acceleration components {a r , a t } for the two mission scenarios.Once the control parameters are selected and the transfer trajectory is obtained, it is possible to evaluate the time-variation of the magnitude of the propulsive acceleration vector a a 2 r + a 2 t during the transfer.Figure 11 shows the values of a as a function of time in the four cases previously described.The curves depicted in that figure can be used to evaluate, a posteriori, the feasibility of the obtained transfer trajectory according to the actual thruster installed on board.In particular, Figure 11a (or Figure 11b) indicates that the maximum value of a during an Earth-Venus (or Earth-Mars) transfer with K = K v is about 0.6 mm/s 2 (or 0.25 mm/s 2 ), while Figure 11c (or Figure 11d) shows that the maximum value of a is roughly 4.2 mm/s 2 (or 2.5 mm/s 2 ) for an Earth-Venus (or Earth-Mars) case when τ f = τ H . Therefore, if, for example, the installed thruster gives a maximum propulsive acceleration of 0.3 mm/s 2 , when K = K v , one concludes that the transfer trajectory obtained in the Earth-Mars case can be theoretically flown, while the result in the Earth-Venus scenario gives a trajectory that violates the propulsive constraint.In the latter case (that is, in the Earth-Venus scenario with K = K v ), the designer could suitably change the control law parameters in order to reduce the maximum value of a reached during the transfer.For example, when K = 0.03 and β = 1.48, the maximum value of a reduces to about 0.29 mm/s 2 , while the flight time rises to roughly 1413 days.The previous results may be easily extended to trade-off solutions between the flight time and the total velocity change necessary to complete the transfer.Recall in fact that, for a given value of ρ, the flight time is a function of K according to Equation (66), while the total velocity change depends on K as Figure 6 shows.Therefore, it is possible to plot the value of ∆v| β=β as a function of τ f .The results are shown in Figure 12, where K ranges within the interval [K v , 1] and the black squares correspond to the cases in which the flight time equals the Hohmann transfer time.Figure 12 represents a simple and effective means to identify reasonable compromise solutions, which are useful in a preliminary analysis of the trajectory design.

Conclusions
In this paper, we have investigated the capabilities of a sliding mode control technique of rapidly generating a possible trajectory for a spacecraft in a typical circle-to-circle orbit transfer mission scenario.Such a control strategy, which is a robust technique usually used for controlling nonlinear systems affected by disturbances, has been employed here to obtain simple mathematical relations and graphs that allow the designer to estimate the (possible) transfer trajectory characteristics as a function of few tuning parameters.
The proposed approach has shown to be effective even in the presence of bounded disturbances due to measurement errors and low-frequency sampling.It may be effectively employed in an early phase of trajectory planning, that is, just before the usual refinement phase that provides the nominal spacecraft trajectory to be tracked during the transfer.In particular, the discussed approach has the scope of reducing the complexity of the mathematical model and the computational cost required to obtain a possible solution to the transfer problem.However, rather strong simplifying assumptions have been adopted, such as the use of a sort of ideal thruster with a freely steerable thrust vector, or the definition of the error in terms of the desired states.A more accurate estimate of the actual propulsive acceleration profile can be obtained by relaxing some of those assumption.This aspect represents the natural extension of this work.

Figure 1 .
Figure 1.Reference frame and conceptual scheme of the two-dimensional mission scenario.

Figure 7 .
Figure 7. Transfer trajectories when K = K v in the two interplanetary mission scenarios.(a) Earth-Venus case; (b) Earth-Mars case.

Figure 8 .
Figure 8. Time-variations of a r and a t when K = K v in two typical interplanetary mission scenarios.(a) Radial component, Earth-Venus case; (b) transverse component, Earth-Venus case; (c) radial component, Earth-Mars case; (d) transverse component, Earth-Mars case.

Figure 12 .
Figure 12.Trade-off solution between flight time and total velocity change in two typical interplanetary mission scenarios.(a) Earth-Venus case; (b) Earth-Mars case.
Mars case.In these cases, by using Equation (66), the values of λ and K, respectively referred to as λ H and K H , are given by convergence rate of x 1 and x 2 λ H value of λ corresponding to τ f = τ H µ primary body gravitational parameter (km 3 /s 2 ) ρ ratio of r f to r 0 λ