CL_man_dvBiElliptic

Description

• Computes the maneuvers of a bi-elliptical transfer from a circular orbit with semi-major axis ai to a circular orbit with semi-major axis af.

  The apogee radius of the elliptical transfer orbit is rt.

  deltav is the sum of the norms of the velocity increments.

  Velocity increments are expressed in cartesian coordinates in the "qsw" local frame.

  If the argument res is present and is equal to "s", all the output data are returned in a structure.

  

Parameters

ai :

Semi-major axis of initial circular orbit [m] (1xN or 1x1)

af :

Semi-major axis of final circular orbit [m] (1xN or 1x1)

rt :

Radius at the position of the second maneuver [m] (1xN or 1x1)

mu :

(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)

res :

(string, optional) Type of output: "d" or "s" for . Default is "d".

deltav :

Som of norms of velocity increments (=|dv1|+|dv2|+|dv3]) [m/s] (1xN)

dv1 :

First velocity increment in cartesian coordinates in the "qsw" frame [m/s] (3xN)

dv2 :

Second velocity increment in cartesian coordinates in the "qsw" frame [m/s] (3xN)

dv3 :

Third velocity increment in cartesian coordinates in the "qsw" frame [m/s] (3xN)

anv1 :

True anomaly at the position of the 1st maneuver. The initial orbit is circular so the value is arbitrary and set to 0 (1xN)

anv2 :

True anomaly at the position of the 2nd maneuver (either 0 or %pi) [rad] (1xN)

anv3 :

True anomaly at the position of the 3rd maneuver (either 0 or %pi) [rad] (1xN)

man :

Structure containing all the output data.

Authors

• CNES - DCT/SB

Bibliography

• 1) Orbital Mechanics for engineering students, H D Curtis, Chapter 6

See also

• CL_man_dvHohmann

• CL_man_dvHohmannG

• CL_man_dvSma

Examples

// 7000 km to 98 000km through a 280 000 transfer orbit: ai = 7000.e3; af = 98000.e3; rt = 280000.e3; [deltav,dv1,dv2,dv3,anv1,anv2,anv3] = CL_man_dvBiElliptic(ai,af,rt) // Check results: kep = [ai ; 0 ; %pi/2 ; 0 ; 0 ; anv1]; kep1 = CL_man_applyDvKep(kep,dv1); kep1(6) = anv2; kep2 = CL_man_applyDvKep(kep1,dv2); kep2(6) = anv3; kep3 = CL_man_applyDvKep(kep2,dv3) // Same example with a Hohmann transfer: // more expensive ! [deltav,dv1,dv2,anv1,anv2] = CL_man_dvHohmann(ai,af)