Physics:Appell's equation of motion

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 and Paul Émile Appell in 1900.

Statement

The Gibbs-Appell equation reads

[math]\displaystyle{ Q_{r} = \frac{\partial S}{\partial \alpha_{r}}, }[/math]

where [math]\displaystyle{ \alpha_r=\ddot{q_r} }[/math] is an arbitrary generalized acceleration, or the second time derivative of the generalized coordinates [math]\displaystyle{ q_r }[/math], and [math]\displaystyle{ Q_r }[/math] is its corresponding generalized force. The generalized force gives the work done

[math]\displaystyle{ dW = \sum_{r=1}^{D} Q_{r} dq_{r}, }[/math]

where the index [math]\displaystyle{ r }[/math] runs over the [math]\displaystyle{ D }[/math] generalized coordinates [math]\displaystyle{ q_r }[/math], which usually correspond to the degrees of freedom of the system. The function [math]\displaystyle{ S }[/math] is defined as the mass-weighted sum of the particle accelerations squared,

[math]\displaystyle{ S = \frac{1}{2} \sum_{k=1}^{N} m_{k} \mathbf{a}_{k}^{2}\,, }[/math]

where the index [math]\displaystyle{ k }[/math] runs over the [math]\displaystyle{ K }[/math] particles, and

[math]\displaystyle{ \mathbf{a}_k = \ddot{\mathbf{r}}_k = \frac{d^2 \mathbf{r}_k}{dt^2} }[/math]

is the acceleration of the [math]\displaystyle{ k }[/math]-th particle, the second time derivative of its position vector [math]\displaystyle{ \mathbf{r}_k }[/math]. Each [math]\displaystyle{ \mathbf{r}_k }[/math] is expressed in terms of generalized coordinates, and [math]\displaystyle{ \mathbf{a}_k }[/math] 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. Moreover, it can be used to derive Kane's equations, which are particularly suited for describing the motion of complex spacecraft. Appell's formulation is an application of Gauss' principle of least constraint.

Derivation

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

[math]\displaystyle{ d\mathbf{r}_{k} = \sum_{r=1}^{D} dq_{r} \frac{\partial \mathbf{r}_{k}}{\partial q_{r}} }[/math]

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

[math]\displaystyle{ \frac{\partial \mathbf{a}_{k}}{\partial \alpha_{r}} = \frac{\partial \mathbf{r}_{k}}{\partial q_{r}} }[/math]

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

[math]\displaystyle{ dW = \sum_{r=1}^{D} Q_{r} dq_{r} = \sum_{k=1}^{N} \mathbf{F}_{k} \cdot d\mathbf{r}_{k} = \sum_{k=1}^{N} m_{k} \mathbf{a}_{k} \cdot d\mathbf{r}_{k} }[/math]

where Newton's second law for the kth particle

[math]\displaystyle{ \mathbf{F}_k = m_k\mathbf{a}_k }[/math]

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

[math]\displaystyle{ dW = \sum_{r=1}^{D} Q_{r} dq_{r} = \sum_{k=1}^{N} m_{k} \mathbf{a}_{k} \cdot \sum_{r=1}^{D} dq_{r} \left( \frac{\partial \mathbf{r}_{k}}{\partial q_{r}} \right) = \sum_{r=1}^{D} dq_{r} \sum_{k=1}^{N} m_{k} \mathbf{a}_{k} \cdot \left( \frac{\partial \mathbf{r}_{k}}{\partial q_{r}} \right) }[/math]

Therefore, the generalized forces are

[math]\displaystyle{ Q_{r} = \sum_{k=1}^{N} m_{k} \mathbf{a}_{k} \cdot \left( \frac{\partial \mathbf{r}_{k}}{\partial q_{r}} \right) = \sum_{k=1}^{N} m_{k} \mathbf{a}_{k} \cdot \left( \frac{\partial \mathbf{a}_{k}}{\partial \alpha_{r}} \right) }[/math]

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

[math]\displaystyle{ \frac{\partial S}{\partial \alpha_{r}} = \frac{\partial}{\partial \alpha_{r}} \frac{1}{2} \sum_{k=1}^{N} m_{k} \left| \mathbf{a}_{k} \right|^{2} = \sum_{k=1}^{N} m_{k} \mathbf{a}_{k} \cdot \left( \frac{\partial \mathbf{a}_{k}}{\partial \alpha_{r}} \right) }[/math]

yielding Appell's equation of motion

[math]\displaystyle{ \frac{\partial S}{\partial \alpha_{r}} = Q_{r}. }[/math]

Examples

Euler's equations of rigid body dynamics

Euler's equations provide an excellent illustration of Appell's formulation.

Consider a rigid body of N particles joined by rigid rods. The rotation of the body may be described by an angular velocity vector [math]\displaystyle{ \boldsymbol\omega }[/math], and the corresponding angular acceleration vector

[math]\displaystyle{ \boldsymbol\alpha = \frac{d\boldsymbol\omega}{dt} }[/math]

The generalized force for a rotation is the torque [math]\displaystyle{ \textbf{N} }[/math], since the work done for an infinitesimal rotation [math]\displaystyle{ \delta \boldsymbol\phi }[/math] is [math]\displaystyle{ dW = \mathbf{N} \cdot \delta \boldsymbol\phi }[/math]. The velocity of the [math]\displaystyle{ k }[/math]-th particle is given by

[math]\displaystyle{ \mathbf{v}_{k} = \boldsymbol\omega \times \mathbf{r}_{k} }[/math]

where [math]\displaystyle{ \mathbf{r}_{k} }[/math] is the particle's position in Cartesian coordinates; its corresponding acceleration is

[math]\displaystyle{ \mathbf{a}_{k} = \frac{d\mathbf{v}_{k}}{dt} = \boldsymbol\alpha \times \mathbf{r}_{k} + \boldsymbol\omega \times \mathbf{v}_{k} }[/math]

Therefore, the function [math]\displaystyle{ S }[/math] may be written as

[math]\displaystyle{ S = \frac{1}{2} \sum_{k=1}^{N} m_{k} \left( \mathbf{a}_{k} \cdot \mathbf{a}_{k} \right) = \frac{1}{2} \sum_{k=1}^{N} m_{k} \left\{ \left(\boldsymbol\alpha \times \mathbf{r}_{k} \right)^{2} + \left( \boldsymbol\omega \times \mathbf{v}_{k} \right)^{2} + 2 \left( \boldsymbol\alpha \times \mathbf{r}_{k} \right) \cdot \left(\boldsymbol\omega \times \mathbf{v}_{k}\right) \right\} }[/math]

Setting the derivative of S with respect to [math]\displaystyle{ \boldsymbol\alpha }[/math] equal to the torque yields Euler's equations

[math]\displaystyle{ I_{xx} \alpha_{x} - \left( I_{yy} - I_{zz} \right)\omega_{y} \omega_{z} = N_{x} }[/math]

[math]\displaystyle{ I_{yy} \alpha_{y} - \left( I_{zz} - I_{xx} \right)\omega_{z} \omega_{x} = N_{y} }[/math]

[math]\displaystyle{ I_{zz} \alpha_{z} - \left( I_{xx} - I_{yy} \right)\omega_{x} \omega_{y} = N_{z} }[/math]

See also

• Principle of stationary action
• Analytical mechanics

Further reading

• Pars, LA (1965). A Treatise on Analytical Dynamics. Woodbridge, Connecticut: Ox Bow Press. pp. 197–227,631–632. 
• Whittaker, ET (1937). A Treatise on the Analytical Dynamics of Particles and Rigid Bodies, with an Introduction to the Problem of Three Bodies (4th ed.). New York: Dover Publications. ISBN. 
• Seeger (1930). "Appell's equations". Journal of the Washington Academy of Sciences 20: 481–484. 
• Brell, H (1913). "Nachweis der Aquivalenz des verallgemeinerten Prinzipes der kleinsten Aktion mit dem Prinzip des kleinsten Zwanges". Wien. Sitz. 122: 933–944.  Connection of Appell's formulation with the principle of least action.
• PDF copy of Appell's article at Goettingen University
• PDF copy of a second article on Appell's equations and Gauss's principle