FUB-HEP/94-5 Euler Equations for Rigid-Body -- a Case for Autoparallel Trajectories in Spaces with Torsion next up previous
Next: References

Institut für Theoretische Physik,
Freie Universität Berlin
Arnimallee 14, D - 14195 Berlin

  FUB-HEP/94-5
Euler Equations for Rigid-Body --
a Case for Autoparallel Trajectories in Spaces with Torsion

P. Fizievgif and H. Kleinertgif

July 10, 1998

Abstract:

A recently discovered action principle for the motion of particles in spaces with curvature and torsion states that the classical trajectories of spinless point particles are straightest paths (autoparallels), not shortest paths (geodesics), as commonly believed. To illustrate the correctness of the new action principle, it is applied to the motion of a rigid body in the body-fixed coordinate system where it yields the correct Euler equations of motion as autoparallels rather than geodesics.

03.20.+i
04.20.Fy
02.40.+m

1) Introduction
In the literature on gravity with curvature and torsion [1], there is a widespread belief that in spaces with torsion (for the geometry of such spaces see [2]), spinless particles move on shortest paths [3]. However, it was discovered in Ref. [4] (in a completely different context, when solving the path integral of the hydrogen atom) that the correct trajectories are the straightest paths in a given geometry, thus calling for a revision of the field-theoretic arguments in [3].

In Ref. [5], a classical action principle was found to comply with this physical fact. The key to the new action principle is the observation, that a space tex2html_wrap_inline893 with nonzero Riemann curvature and torsion may be mapped locally into a euclidean space tex2html_wrap_inline895 by an anholonomic transformation tex2html_wrap_inline897 [6], [4]. For the particle with mass m in tex2html_wrap_inline895, the equations of motion tex2html_wrap_inline903 yield straight-line trajectories. Under an anholonomic transformation, these lines go over into autoparallels satisfying the equation
 equation37
with tex2html_wrap_inline905 being the Cartan connection (tex2html_wrap_inline907 is the inverse of tex2html_wrap_inline909, i.e., tex2html_wrap_inline911). The cartan connection can be written as tex2html_wrap_inline913, where tex2html_wrap_inline915 is the Levi-Cevita connection, the Christoffel symbol tex2html_wrap_inline917 which depends only on the metric tensor tex2html_wrap_inline919, and tex2html_wrap_inline921 is a combination of torsion tensors, the so-called contortion tensor tex2html_wrap_inline923 Although the torsion tensor tex2html_wrap_inline925 is the antisymmetric part of the connection, tex2html_wrap_inline927, it does contribute to the symmetric part of the affine connection:
 equation94
Thus it contributes the equation for the autoparallel (1).

It is a well-known fact of classical mechanics, that equations of motion can be transformed nonholonomically without changing the physical trajectories.

This is in contrast to the behavior of the action principle under nonholonomic transformations. In a flat space, the orbits in the euclidean space tex2html_wrap_inline895 may be derived from Hamilton's principle tex2html_wrap_inline931 applied to the classical action tex2html_wrap_inline933. Under the anholonomic transformation tex2html_wrap_inline935, this action goes into the curvilinear form tex2html_wrap_inline937. In the space tex2html_wrap_inline893, a naive application of Hamilton's principle tex2html_wrap_inline941 would produce wrong equations of motion in the space tex2html_wrap_inline893. One would find tex2html_wrap_inline945 which are the equations for the geodesics rather than the autoparallels; they lack torsion force tex2html_wrap_inline947 contained in tex2html_wrap_inline949 via (2).

The problem of describing the dynamics of a rigid body within the body-fixed frame may be used as a case to find out which of the trajectories is correct. In the usual parameter space of Euler angles, there is no torsion and a constant Riemann curvature [7]. Thus the trajectories are geodesics. This changes when describing the system within the body-fixed frame. The transformation between the coordinates is non-holonomic and the new space possesses torsion. We shall apply the variational principle of Ref. [5] and derive, within the body-fixed reference system both the Euler equations for the angular momentum and the equations for the translational motion and show, that the rigid body moves along autoparallel trajectories in the body-fixed reference system.

2) The two coordinate frames of the Rigid Body
Since their discovery more than two centuries ago, the Euler equations
 equation139
have played a key role in understanding the rotations of a rigid body around a fixed point O. (see for example [8] - [12] and the references therein). The vectors refer to the body-fixed system, tex2html_wrap_inline951 being the angular momentum, tex2html_wrap_inline953 the instantaneous angular velocity, and tex2html_wrap_inline955 the moments of the external forces. The dot denotes differentiation with respect to the time t. The unit basis vectors tex2html_wrap_inline959 in B may be assumed to point along the eigenvectors of the inertia tensor with respect to O. Then the angular moment tex2html_wrap_inline951 has the components tex2html_wrap_inline963 (no sum over i).

To describe also the translational motion of the rigid body one usually chooses the point O to coincide with the center of mass moving through space. The motion satisfies the additional equations:
 equation152
where tex2html_wrap_inline967 is the linear momentum, M the body's total mass, tex2html_wrap_inline971 the velocity of the center of mass, and tex2html_wrap_inline973 the driving force. Presently, the equations (3) and (4) play an important role in missile dynamics analysis [12].

Usually, the two sets of equations are derived as the Newton equations of motion in the stationary reference system S:
  eqnarray167
Let tex2html_wrap_inline977 and tex2html_wrap_inline979 be the angular and linear velocities in S and tex2html_wrap_inline983 those in B. They are related by a linear transformation (see [8] - [10]):
  eqnarray177
The tex2html_wrap_inline987 matrix elements tex2html_wrap_inline989 and tex2html_wrap_inline991 depend on the coordinates tex2html_wrap_inline993 in the system S. They satisfy tex2html_wrap_inline997, tex2html_wrap_inline999, making the transformation (7) and (8) anholonomic coordinate transformations (see [8]-[14] for the general theory of the dynamics in anholonomic coordinates).

Under the transformation (7), the equations (5) go over into (1). The term tex2html_wrap_inline1001 describes the moment of the gyroscopic force. The additional term tex2html_wrap_inline1003 in equation (4) arise similarly. Both terms are a consequence of the anholonomy of the transformation (7), (8).

Within the stationary system S, the equations of motion (5) and (6) can be derived via Hamilton's action principle. If one transforms the classical action to the system B, however, the coordinates become nonholonomic and a naive application of Hamilton's principle would produce wrong equations of motion lacking the additional gyroscopic moments.

In 1901, Poincaré showed [13] how to vary an action expressed in terms of nonholonomic coordinates. Following his treatment, one certainly recovers the gyroscopic forces. His treatment is reviewed in the textbook [14].

We shall now first demonstrate that Poincaré's treatment may be viewed geometrically as an application of a recently proposed action principle in spaces with curvature and torsion.

Let us rederive the equations of motion in the appropriate geometric language.

3) The tex2html_wrap_inline1009 geometry in a rotating body-fixed system produced by an anholonomic transformation from the rest system
Consider first a rigid body rotating around a fixed point. Within the body-fixed system B, we introduce anholonomic coordinates tex2html_wrap_inline1011 corresponding to the transformation (7). They are the components of the body's rotation vector in the axis-angle parametrization. The anholonomic transformation (7) defines their infinitesimal increments tex2html_wrap_inline1013. For a precise specification, let us go to the system S where the standard Euler angles tex2html_wrap_inline1017 [10] parametrize (holonomically) the body's configuration space tex2html_wrap_inline1019. The components of the angular velocity in the system B are tex2html_wrap_inline1021 (with the short notation tex2html_wrap_inline1023). The basic relation
 equation219
has coefficients tex2html_wrap_inline989 which form the matrix
 equation224
The symbol tex2html_wrap_inline1027 stands for increments do not belong to an integrable function; they do not satisfy the Schwarz integrability condition which would read, with the coefficients of (7), tex2html_wrap_inline1029.

The metric tex2html_wrap_inline1031 and its inverse tex2html_wrap_inline1033 have the matrices
 equation235
where tex2html_wrap_inline1035 is the determinant of g. The space of nonholonomic coordinates is flat. The space of Euler angles has a constant Riemann curvature. From g we find the Christoffel symbol:
displaymath1041
The associated Riemannian Ricci tensor is tex2html_wrap_inline1043 with the constant contraction tex2html_wrap_inline1045.

On the other hand, the space of Euler angles is free of torsion. Since the mapping is nonholonomic, the space tex2html_wrap_inline1047 can carry torsion. To calculate this we observe that the only nonzero components of the affine-flat Cartan connection are
 equation275
Using the metric tex2html_wrap_inline1049 we find the antisymmetric part is the torsion tensor tex2html_wrap_inline1051, where tex2html_wrap_inline1053 denotes the Levi-Civita antisymmetric tensor tex2html_wrap_inline1055 [2].

Since torsion involves a single derivative of the transformation matrices between two sets of coordinates, the chain rule of differentiation makes it an additive quantity. With the tex2html_wrap_inline1057-space being free of torsion, the tex2html_wrap_inline1047-space must have a torsion which is the negative of tex2html_wrap_inline1061 after transforming the indices appropriately:
 equation306
Indeed, inverting (7) to tex2html_wrap_inline1063 the torsion arising when going from tex2html_wrap_inline1065- to tex2html_wrap_inline1047-coordinates is by definition
 equation315
which due to tex2html_wrap_inline1069 and tex2html_wrap_inline1071 coincides with (13). Explicitly, we find from the above tex2html_wrap_inline1061 using the metric tex2html_wrap_inline1075:
 equation336

4) Action principle for the rotational motion in the body-fixed tex2html_wrap_inline1077
Pure rotations around an arbitrary fixed point of the body are governed by the Euler equations (3). The kinetic energy is tex2html_wrap_inline1079 ( with an extended Einstein's summation convention to include also diagonal components of the inertia tensor).

An application of Hamilton's action principle tex2html_wrap_inline1081 to the classical action tex2html_wrap_inline1083 in the system S certainly yields the correct equations of motion. Under the transformation (7), they become the Euler equations (3).

Let us now transform the action tex2html_wrap_inline1087 via (7) to the body-fixed system B. The result is very simple:
 equation353
By applying naively Hamilton's action principle we would find the equations tex2html_wrap_inline1089 for each i which are not the correct Euler equations (3) -- the gyroscopic moments being missed. The contradiction is caused by the fact that the variations tex2html_wrap_inline1093 in the space of Euler angles tex2html_wrap_inline1095 and the variations tex2html_wrap_inline1097 in the space of anholonomic coordinates tex2html_wrap_inline1099, are related with each other in a path dependent way (``nonlocal" on the time axis) [5]. Explicitly, there is the functional equation
 equation367
Its most important consequence is that closed paths in the space tex2html_wrap_inline1095 are not mapped into closed paths in the space tex2html_wrap_inline1099, due to a nonvanishing Burgers vector tex2html_wrap_inline1105.

The usual Hamilton action principle in the system S proceeds by considering a variation tex2html_wrap_inline1109 between two paths tex2html_wrap_inline1111 tex2html_wrap_inline1113 with common ends: tex2html_wrap_inline1115. Together, they form a closed path in the space tex2html_wrap_inline1095. The above naive application of Hamilton's action principle in the system B employed analogous closed-path variations tex2html_wrap_inline1119 between two paths tex2html_wrap_inline1121 tex2html_wrap_inline1113, with also common ends: tex2html_wrap_inline1125. This, however, cannot be correct. Under a transformation (7), variations tex2html_wrap_inline1127 do not produce closed-path variations tex2html_wrap_inline1129 in the space tex2html_wrap_inline1099. The nonzero Burgers vector tex2html_wrap_inline1133 of the anholonomy causes a closure failure. This has to be accounted for in a correct derivation of the equations of motion.

The anholonomy of the transformation (7) requires distinguishing two types of variations of the paths in the space tex2html_wrap_inline1099: closed path variations tex2html_wrap_inline1097 of the tpye described above, and anholonomic variations tex2html_wrap_inline1139, which are images of the tex2html_wrap_inline1093 variations in the space tex2html_wrap_inline1095. Only the latter will produce the correct equations of motion from an action principle in the system B.

To calculate the variation of the classical action in the system B we derive the following simple formula for the anholonomic variation tex2html_wrap_inline1145 of the angular velocity tex2html_wrap_inline1147:
 equation390
(corresponding to formula (9) in Ref. [5]).

Indeed, the transformation (7) and the definition (9) lead to the equation
 equation399
which by (13) becomes
 equation411
leading to (18) via (15). The variations tex2html_wrap_inline1127 are still performed in the system S. They can be mapped locally into holonomic closed-path variations in the system B, tex2html_wrap_inline1151, which have the usual fixed-endpoint property tex2html_wrap_inline1153 reflecting the closed-path condition tex2html_wrap_inline1155 in the system S.

Let us emphasize that our nonholonomic variations (15) are completely intrinsic to the system B. By applying these to the action tex2html_wrap_inline1157, an integration by parts with the boundary conditions tex2html_wrap_inline1153 leads to
 equation426
Since tex2html_wrap_inline1129 are now ordinary (holonomic) variations, we find the Euler equations (1). Thus the anholonomic action principle
 equation436
produces the correct equations of motion in the body-fixed system B.
5) Autoparallels
Let us now verify that these body-fixed equations of motion correspond to the equations for autoparalles rather than geodesics. The Lagrangian of a symmetric rigid body in the body-fixed coordinate system has the general form tex2html_wrap_inline1163 with the generalized coordinates tex2html_wrap_inline1165 and the metric of the kinetic energy tex2html_wrap_inline1167. The equations of motion for autoparallels in this space are independent of I and read tex2html_wrap_inline1171. Now, since tex2html_wrap_inline1173 is a constant, the space has obviously a vanishing Riemann connection. Because of (15), it has a nonvanishing torsion. Using the decomposition (2), the equation for the autoparallels reduce to tex2html_wrap_inline1175. The antisymmetry of the torsion tensor (15) causes the torsion force to vanish, i.e., the autoparallels happen to coincide with the geodesics.

Let us now go over to the asymmetric rigid body with the metric tex2html_wrap_inline1177 (no sum over i). This is done by an additional holonomic transformation tex2html_wrap_inline1181. Since the transformation coefficients are independent of the coordinates, this does not change the torsion of the space. Since the equations of motion contain the torsion tensor with the first index raised by tex2html_wrap_inline1183 and the second lowered by tex2html_wrap_inline1185 [see Eq. (2)], the torsion tensor tex2html_wrap_inline1187 changes into tex2html_wrap_inline1189 and the equation of motion becomes tex2html_wrap_inline1191. The autoparallel is now different from the geodesic by containing the effect of the gyroscopic force.
6) Action principle for general motion in the body-fixed tex2html_wrap_inline1077
To extend the action principle to the general motion of a rigid body one usually chooses the point O in B to be the center of mass of the body and studies the movements of O and rotations around O. Then tex2html_wrap_inline1195 are the body's principal moments of inertia. The configuration space is a 6-dimensional manifold tex2html_wrap_inline1197 and the kinetic energy tex2html_wrap_inline1199 consists of two terms: The first is the translational kinetic energy tex2html_wrap_inline1201 in the system B the second is the rotational kinetic energy tex2html_wrap_inline1203 in the system B.

In the system S, the position and orientation of the body are parametrized by the (holonomic) coordinates tex2html_wrap_inline1207, with tex2html_wrap_inline1209 being the mass center coordinates and tex2html_wrap_inline1057 the Euler angles.

The body-fixed basis vectors in B (see Fig. 1) are some function of the Euler angles tex2html_wrap_inline1213. Their components in the cartesian basis of the system S are tex2html_wrap_inline1217. They form a tex2html_wrap_inline987 orthogonal matrix tex2html_wrap_inline1221. The components of the velocity of the center of mass in the system S are tex2html_wrap_inline1225. The components of tex2html_wrap_inline1227 with respect to the basis tex2html_wrap_inline1229 are tex2html_wrap_inline1231.

The last relation permits us to introduce a new set of anholonomic coordinates tex2html_wrap_inline1233 describing the center of mass motion in the system S as seen from the system B. By analogy with formula (9), we define the infinitesimal increments tex2html_wrap_inline1235 as:
 equation498
The set tex2html_wrap_inline1237 gives a complete set of an anholonomic coordinates for the configuration space tex2html_wrap_inline1239 of the moving body.

Let us write down the anholonomic geometry of this space. Introducing tex2html_wrap_inline1241 matrices, whose rows and columns are labeled by capital Latin and Greek letters,
 equation507
the geometry possesses the affine connection tex2html_wrap_inline1243 with the torsion tensor tex2html_wrap_inline1245 and a vanishing Cartan curvature tensor tex2html_wrap_inline1247, implying an nonvanishing Riemann curvature tensor tex2html_wrap_inline1249. This geometry is a simple extension of the 3-dimensional anholonomic geometry on tex2html_wrap_inline1009 described in Section 2. It reflects the structure of the present configuration space tex2html_wrap_inline1253, i.e., the space of configurations of the rigid body comprising translations of the center of mass and rotations around it.

As a direct consequence of the definition (23), we have the relation tex2html_wrap_inline1255 which, together with (9), leads to the following simple form of the rigid body's Lagrangian in the system B:


 equation544

As before, a naive application of Hamiltons principle in the system B would produce a wrong equation tex2html_wrap_inline1257 for the motion of the center-of-mass. The correct anholonomic variation of the velocity
 equation551
Indeed, the transformation (8) and the definition (23) lead to
 equation560
The variations tex2html_wrap_inline1259 may be mapped locally into holonomic closed-path variations tex2html_wrap_inline1261 in the system B with the property tex2html_wrap_inline1263 reflecting the closed path condition tex2html_wrap_inline1265. Formula (23) follows using the equations tex2html_wrap_inline1267, and the relation tex2html_wrap_inline1269.

By applying the anholonomic variation (23) to the translational energy, we obtain
 equation582
An integration by parts using of the fixed-end conditions tex2html_wrap_inline1271 leads to the following expression for the total anholonomic variation of the action tex2html_wrap_inline1273:
 eqnarray593
From the second term, we find the equations (2) for the translational motion.

By an extension of the geometric arguments in Sec. 2, the extended equations of motion describe again autoparallels, not geodesics.
5) Conclusion
By subjecting the action tex2html_wrap_inline1157 in the body-fixed system to the new anholonomic variations with respect to translational and rotational degrees of freedom of the rigid body we have been able to derive both the correct Euler equations (3) and the equations (4) completely within the body-fixed system without reference to the stationary systems. Geometrically, the equations of motion in the body-fixed coordinates correspond to autoparralels rather than geodesics.

This has important implications on the classical field theory of gravity in spaces with curvature and torsion. It calls for a revision of the customary variational derivations of the field equations [1]. We have just learned that these variations require the knowledge of the torsion in the "superspace" inhabited by the metric tensor and the torsion field. DeWitt has shown us how to construct the "supermetric" in a the space of metrics with Riemannian geometry. His construction will require an extension to spaces with torsion, before it will be possible to derive correct field equations for metric and the torsion fields from a given field action.

Apart from that, the above action principle for the rigid body which is intrinsic to the body-fixed system may have practical consequences for variational calculations of the mechanical motion of rigid bodies in accelerated frames of references, for instance encounter maneuvers in satellite physics.  
 
 
Acknowledgement:
We are greateful to Dr. S.V. Shabanov for useful discussions.




next up previous
Next: References

Hagen Kleinert
Fri Jul 10 11:41:47 MET DST 1998