Appell's equation of motion explained

In classical mechanics, Appell's equation of motion (aka the Gibbs–Appell equation of motion) is an alternative general formulation of classical mechanics described by Josiah Willard Gibbs in 1879[1] and Paul Émile Appell in 1900.[2]

Statement

The Gibbs-Appell equation reads

Qr=

\partialS
\partial\alphar

,

where

\alphar=\ddot{q}r

is an arbitrary generalized acceleration, or the second time derivative of the generalized coordinates

qr

, and

Qr

is its corresponding generalized force. The generalized force gives the work done

dW=

D
\sum
r=1

Qrdqr,

where the index

r

runs over the

D

generalized coordinates

qr

, which usually correspond to the degrees of freedom of the system. The function

S

is defined as the mass-weighted sum of the particle accelerations squared,

S=

1
2
N
\sum
k=1

mk

2
a
k

,

where the index

k

runs over the

K

particles, and

ak=\ddot{r

}_k = \frac

is the acceleration of the

k

-th particle, the second time derivative of its position vector

rk

. Each

rk

is expressed in terms of generalized coordinates, and

ak

is expressed in terms of the generalized accelerations.

Relations to other formulations of classical mechanics

Appell's formulation does not introduce any new physics to classical mechanics and as such is equivalent to other reformulations of classical mechanics, such as Lagrangian mechanics, and Hamiltonian mechanics. All classical mechanics is contained within Newton's laws of motion. In some cases, Appell's equation of motion may be more convenient than the commonly used Lagrangian mechanics, particularly when nonholonomic constraints are involved. In fact, Appell's equation leads directly to Lagrange's equations of motion.[3] Moreover, it can be used to derive Kane's equations, which are particularly suited for describing the motion of complex spacecraft.[4] Appell's formulation is an application of Gauss' principle of least constraint.[5]

Derivation

The change in the particle positions rk for an infinitesimal change in the D generalized coordinates is

drk=

D
\sum
r=1

dqr

\partialrk
\partialqr

Taking two derivatives with respect to time yields an equivalent equation for the accelerations

\partialak
\partial\alphar

=

\partialrk
\partialqr

The work done by an infinitesimal change dqr in the generalized coordinates is

dW=

D
\sum
r=1

Qrdqr=

N
\sum
k=1

Fkdrk=

N
\sum
k=1

mkakdrk

where Newton's second law for the kth particle

Fk=mkak

has been used. Substituting the formula for drk and swapping the order of the two summations yields the formulae

dW=

D
\sum
r=1

Qrdqr=

N
\sum
k=1

mkak

D
\sum
r=1

dqr\left(

\partialrk
\partialqr

\right)=

D
\sum
r=1

dqr

N
\sum
k=1

mkak\left(

\partialrk
\partialqr

\right)

Therefore, the generalized forces are

Qr=

N
\sum
k=1

mkak\left(

\partialrk
\partialqr

\right)

N
= \sum
k=1

mkak\left(

\partialak
\partial\alphar

\right)

This equals the derivative of S with respect to the generalized accelerations

\partialS
\partial\alphar

=

\partial
\partial\alphar
1
2
N
\sum
k=1

mk\left|ak\right|2=

N
\sum
k=1

mkak\left(

\partialak
\partial\alphar

\right)

yielding Appell's equation of motion

\partialS
\partial\alphar

=Qr.

Examples

Euler's equations of rigid body dynamics

\boldsymbol\omega

, and the corresponding angular acceleration vector

\boldsymbol\alpha=

d\boldsymbol\omega
dt

The generalized force for a rotation is the torque

bf{N}

, since the work done for an infinitesimal rotation

\delta\boldsymbol\phi

is

dW=N\delta\boldsymbol\phi

. The velocity of the

k

-th particle is given by

vk=\boldsymbol\omega x rk

where

rk

is the particle's position in Cartesian coordinates; its corresponding acceleration is

ak=

dvk
dt

=\boldsymbol\alpha x rk+\boldsymbol\omega x vk

Therefore, the function

S

may be written as

S=

1
2
N
\sum
k=1

mk\left(akak\right) =

1
2
N
\sum
k=1

mk\left\{\left(\boldsymbol\alpha x rk\right)2+\left(\boldsymbol\omega x vk\right)2+2\left(\boldsymbol\alpha x rk\right)\left(\boldsymbol\omega x vk\right)\right\}

Setting the derivative of S with respect to

\boldsymbol\alpha

equal to the torque yields Euler's equations

Ixx\alphax-\left(Iyy-Izz\right)\omegay\omegaz=Nx

Iyy\alphay-\left(Izz-Ixx\right)\omegaz\omegax=Ny

Izz\alphaz-\left(Ixx-Iyy\right)\omegax\omegay=Nz

See also

Further reading

. E. T. Whittaker. 1937 . . 4th . Dover Publications . New York . ISBN .

Notes and References

  1. Gibbs. JW. 1879. On the Fundamental Formulae of Dynamics.. American Journal of Mathematics. 2. 1. 49 - 64. 10.2307/2369196. 2369196.
  2. Appell . P . 1900 . Sur une forme générale des équations de la dynamique. . Journal für die reine und angewandte Mathematik . 121 . 310 - ? .
  3. Deslodge. Edward A.. 1988. The Gibbs–Appell equations of motion. American Journal of Physics. 56. 9. 841–46. 10.1119/1.15463. 1988AmJPh..56..841D . 123074999 .
  4. Deslodge. Edward A.. 1987. Relationship between Kane's equations and the Gibbs-Appell equations. Journal of Guidance, Control, and Dynamics. American Institute of Aeronautics and Astronautics. 10. 1. 120–22. 10.2514/3.20192. 1987JGCD...10..120D .
  5. Lewis. Andrew D.. August 1996. The geometry of the Gibbs-Appell equations and Gauss' principle of least constraint. Reports on Mathematical Physics. 38. 1. 11–28. 10.1016/0034-4877(96)87675-0. 1996RpMP...38...11L .