4.2. Spherical shaping¶
The spherical shapebased method shapes the position of the spacecraft using spherical coordinates, and uses the azimuth angle as independent variable. This shaping method assumes an inverse function for the radial distance, and a polynomial function for the elevation angle. It only requires two shaping functions, as one of the spherical coordinates is used as independent variable. As opposed to hodographic shaping, the current implementation of the spherical shaping method does not let the choice of base functions free. The ones proposed in Novak, 2012 (Methods and tools for preliminary low thrust mission analysis), where this shaping method has been initially developed, are used directly. The weighting coefficients of those base functions are calculated so that the boundary conditions of the lowthrust trajectory are fulfilled. The spherical shaping method has no free parameter (on the contrary to the hodographic shaping method), so that the shaped trajectory is entirely defined by its boundary conditions and timeofflight and cannot be optimised in terms of deltaV budget.
One of the main drawbacks of this shaping method is that it assumes the thrust is always tangential to the orbit. Consequently, there is no guarantee this shapebased method always leads to a feasible trajectory, depending on the boundary conditions and required timeofflight. This is especially the case for trajectories implying large inclination changes.
Moreover, because this method shapes the position of the spacecraft and not its velocity, matching the required time of flight is not inherently ensured by the implementation of the shapebased method itself (as it is the case for hodographic shaping for instance). An iterative process is required to match the expected timeofflight. The coefficient associated to the second base function of the inverse polynomial radial distance function (referred to as the “free coefficient” is the following) is thus not determined by the departure and arrival boundary conditions, but its value is tuned until the expected timeofflight is matched.
4.2.1. Setting up a spherically shaped trajectory¶
A spherically shaped trajectory is defined with the SphericalShaping
class. It is a derived class from ShapeBasedMethodLeg
. A corresponding settings class (SphericalShapingLegSettings
) has been implemented to facilitate the creation of SphericalShaping
objects (see Setting up a lowthrust trajectory for more details on the SphericalShapingLegSettings
).

class
SphericalShaping
¶
The SphericalShaping class is defined as follows:
SphericalShaping(Eigen::Vector6d initialState,
Eigen::Vector6d finalState,
double requiredTimeOfFlight,
int numberOfRevolutions,
simulation_setup::NamedBodyMap& bodyMap,
const std::string bodyToPropagate,
const std::string centralBody,
double initialValueFreeCoefficient,
std::shared_ptr< root_finders::RootFinderSettings >& rootFinderSettings,
const double lowerBoundFreeCoefficient = TUDAT_NAN,
const double upperBoundFreeCoefficient = TUDAT_NAN,
std::shared_ptr< numerical_integrators::IntegratorSettings< double > > integratorSettings
= std::shared_ptr< numerical_integrators::IntegratorSettings< > >( ) )
where the inputs are:
initialState
State of the spacecraft at departure.
finalState
State of the spacecraft at arrival.
requiredTimeOfFlight
Time of flight required for the shapebased trajectory.
numberOfRevolutions
Required number of revolutions before the spacecraft reaches its final state.
bodyMap
Map of pointers to
Body
objects involved in the lowthrust trajectory.
bodyToPropagate
Name of the spacecraft to be propagated.
centralBody
Name of the central body of the lowthrust trajectory.
initialValueFreeCoefficient
Initial value for the free coefficient, to be tuned until the required timeofflight is matched.
rootFinderSettings
Settings that define the root finder algorithm, used to find the proper free coefficient value for which the required timeofflight would be achieved.
lowerBoundFreeCoefficient
Lower bound for the free coefficient search.
upperBoundFreeCoefficient
Upper bound for the free coefficient search.
integratorSettings
Integrator settings (empty by default), used to propagate the spacecraft mass or thrust profiles, or to numerically propagate the fully perturbed trajectory (as a means to assess the quality of the analytical shapedbased preliminary design).
Note
The creation of a spherically shaped trajectory requires much less inputs from the user than hodographic shaping does. The shaping functions are indeed fixed in spherical shaping, while they are userdefined in hodogaphic shaping. Morever, spherical shaping does not allow for free parameter, as opposed to hodographic shaping. This makes the definition of a spherically shaped trajectory more straightforward (no need to set up the base functions) from the user’s point of view. The downside of it is that spherical is less flexible than hodographic shaping, and cannot be optimised to get a lower deltaV.