Describing relative motion with differences of orbital elements instead of Cartesian coordinates. Subscript 𝑝: primary (reference), 𝑠: secondary.

Relative Eccentricity and Inclination Vectors

  • Relative orbital elements include:
    • the relative eccentricity vector
    • the relative inclination vector
  • 𝛿𝑖 = angle between the orbital planes; the relative inclination vector is
    • sin(𝛿𝑖)[cos𝜃sin𝜃]
    • 𝜃: argument of latitude of the relative ascending node (the point where the second spacecraft crosses the orbital plane of the first going in the +̂𝐤 direction)

Orbital Element Representations

  • Three families:
    • Relative orbital elements
    • Orbital element differences
    • Geometric ROE (constants of motion of HCW)

Orbital Element Differences

  • Classical elements: take the difference between primary and secondary:
    • [Δ𝑎Δ𝑒Δ𝑖ΔΩΔ𝜔Δ𝑀]
  • Easy to understand, but have singularities (the same ones as the classical elements)
  • Alternatively: equinoctial element differences

Relative Orbital Elements

  • Relative semimajor axis:
    • 𝛿𝑎=𝑎𝑠−𝑎𝑝𝑎𝑝
  • Relative mean longitude, using the mean argument of latitude 𝑢𝑀=𝑀+𝜔:
    • 𝛿𝜆=𝑢𝑀𝑠−𝑢𝑀𝑝+(Ω𝑠−Ω𝑝)cos𝑖𝑝
  • Relative eccentricity vector:
    • 𝛿𝐞=[𝑒𝑠cos𝜔𝑠−𝑒𝑝cos𝜔𝑝𝑒𝑠sin𝜔𝑠−𝑒𝑝sin𝜔𝑝]
  • Relative inclination vector:
    • 𝛿𝐢=[𝑖𝑠−𝑖𝑝(Ω𝑠−Ω𝑝)sin𝑖𝑝]
  • These are commonly used for control of formation-flying or cluster-flight spacecraft

Mapping to Cartesian

  • 𝑥=𝑎𝑝(𝛿𝑎−𝛿𝑒𝑥cos𝑢𝑝−𝛿𝑒𝑦sin𝑢𝑝)𝑦=𝑎𝑝(−32𝛿𝑎𝑢𝑝+𝛿𝜆+2𝛿𝑒𝑥sin𝑢𝑝−2𝛿𝑒𝑦cos𝑢𝑝)𝑧=𝑎𝑝(𝛿𝑖𝑥sin𝑢𝑝−𝛿𝑖𝑦cos𝑢𝑝)
  • 𝑥,𝑦,𝑧: radial, along-track, cross-track; 𝑢𝑝: argument of latitude of the primary
  • −32𝛿𝑎𝑢𝑝: a semimajor-axis difference makes the secondary drift along-track

Geometric Relative Orbital Elements

  • Geometric elements/parameters that define constants of the motion (from the HCW state 𝑥,𝑦,𝑧,̇𝑥,̇𝑦,̇𝑧):
  • Relative ellipse size:
    • 𝑎ell=2√(̇𝑥𝑛)2+(3𝑥+2̇𝑦𝑛)2
  • Center of the projected ellipse:
    • 𝑐ell=[4𝑥+2̇𝑦𝑛𝑦−2̇𝑥𝑛]
  • Relative in-plane phase:
    • 𝛽=arctaṅ𝑥3𝑛𝑥+2̇𝑦
  • Cross-track amplitude and phase:
    • 𝑧max=√𝑧2+(̇𝑧𝑛)2,𝜓=arctan(𝑛𝑧̇𝑧)