Clohessy–Wiltshire equations explained

The Clohessy–Wiltshire equations describe a simplified model of orbital relative motion, in which the target is in a circular orbit, and the chaser spacecraft is in an elliptical or circular orbit. This model gives a first-order approximation of the chaser's motion in a target-centered coordinate system. It is used to plan the rendezvous of the chaser with the target.[1] [2]

History

Early results about relative orbital motion were published by George William Hill in 1878.[3] Hill's paper discussed the orbital motion of the moon relative to the Earth.

In 1960, W. H. Clohessy and R. S. Wiltshire published the Clohessy–Wiltshire equations to describe relative orbital motion of a general satellite for the purpose of designing control systems to achieve orbital rendezvous.[1]

System Definition

Suppose a target body is moving in a circular orbit and a chaser body is moving in an elliptical orbit. Let

x,y,z

be the relative position of the chaser relative to the target with

x

radially outward from the target body,

y

is along the orbit track of the target body, and

z

is along the orbital angular momentum vector of the target body (i.e.,

x,y,z

form a right-handed triad).Then, the Clohessy–Wiltshire equations are\begin\ddot &= 3n^x + 2n\dot \\\ddot &= -2n\dot \\\ddot &= -n^z \endwhere n = \sqrt is the orbital rate (in units of radians/second) of the target body,

a

is the radius of the target body's circular orbit,

\mu

is the standard gravitational parameter,

If we define the state vector as

x=(x,y,z,

x,
y,
z)
, the Clohessy–Wiltshire equations can be written as a linear time-invariant (LTI) system,[4] \dot \mathbf = A\mathbfwhere the state matrix

A

isA = \begin0 & 0 & 0 & 1 & 0 & 0 \\0 & 0 & 0 & 0 & 1 & 0 \\0 & 0 & 0 & 0 & 0 & 1 \\3n^2 & 0 & 0 & 0 & 2 n & 0 \\0 & 0 & 0 & -2 n & 0 & 0 \\0 & 0 & -n^2 & 0 & 0 & 0\end.

For a satellite in low Earth orbit,

\mu=3.986 x 1014

m3/s2
and

a=6,793,137m

, implying

n=0.00113

s-1
, corresponding to an orbital period of about 93 minutes.

If the chaser satellite has mass

m

and thrusters that apply a force

F=(Fx,Fy,Fz),

then the relative dynamics are given by the LTI control system[4] \dot \mathbf = A \mathbf + B \mathbfwhere

u=F/m

is the applied force per unit mass and \boldsymbol = \begin0 & 0 & 0 \\0 & 0 & 0 \\0 & 0 & 0 \\1 & 0 & 0 \\0 & 1 & 0 \\0 & 0 & 1\end.

Solution

We can obtain closed form solutions of these coupled differential equations in matrix form, allowing us to find the position and velocity of the chaser at any time given the initial position and velocity.[5] \begin\delta\vec(t) &= [\Phi_{rr}(t)]\delta\vec + [\Phi_{rv}(t)]\delta\vec \\\delta\vec(t) &= [\Phi_{vr}(t)]\delta\vec + [\Phi_{vv}(t)]\delta\vec\endwhere:\begin\Phi_(t) &= \begin 4-3\cos & 0 & 0 \\ 6(\sin-nt) & 1 & 0 \\ 0 & 0 & \cos \end \\\Phi_(t) &= \begin \frac\sin & \frac(1-\cos) & 0 \\ \frac(\cos-1) & \frac(4\sin-3nt) & 0 \\ 0 & 0 & \frac\sin \end \\\Phi_(t) &= \begin 3n\sin & 0 & 0 \\ 6n(\cos-1) & 0 & 0 \\ 0 & 0 & -n\sin \end \\\Phi_(t) &= \begin \cos & 2\sin & 0 \\ -2\sin & 4\cos-3 & 0 \\ 0 & 0 & \cos \end\endNote that

\Phivr(t)=

\Phi

rr(t)

and

\Phivv(t)=

\Phi

rv(t)

. Since these matrices are easily invertible, we can also solve for the initial conditions given only the final conditions and the properties of the target vehicle's orbit.

See also

References

  1. Clohessy . W. H. . Wiltshire . R. S. . 1960 . Terminal Guidance System for Satellite Rendezvous . Journal of the Aerospace Sciences . 27 . 9 . 653–658 . 10.2514/8.8704.
  2. Web site: Clohessy-Wiltshire equations. University of Texas at Austin. 12 September 2013.
  3. ((Hill, G. W.)) . American Journal of Mathematics . Researches in the Lunar Theory . 1 . 1 . 5–26 . Johns Hopkins University Press . 1878 . 0002-9327 . 10.2307/2369430. 2369430 .
  4. ((Starek, J. A.)), ((Schmerling, E.)), ((Maher, G. D.)), ((Barbee, B. W.)), ((Pavone, M.)) . Journal of Guidance, Control, and Dynamics . Fast, Safe, Propellant-Efficient Spacecraft Motion Planning Under Clohessy–Wiltshire–Hill Dynamics . 40 . 2 . 418–438 . American Institute of Aeronautics and Astronautics . February 2017 . 0731-5090 . 10.2514/1.G001913. 1601.00042 . 2017JGCD...40..418S . 4956601 .
  5. Book: Curtis, Howard D.. Orbital Mechanics for Engineering Students. Elsevier. 2014. 9780080977478. Oxford, UK. 383–387. 3rd.

Further reading

External links