User:Prof McCarthy/generalized coordinates

In analytical mechanics, specifically the study of the rigid body dynamics of multibody systems, the term generalized coordinates refers to the parameters that describe the configuration of the system relative to some reference configuration. These parameters must uniquely define the configuration of the system relative to the reference configuration.^[1]

An example of a generalized coordinate is the angle that locates a point moving on a circle, in contrast to its x and y coordinates. The adjective "generalized" distinguishes these parameters from the traditional use of the term coordinate as measured against a specified line, such as Cartesian coordinates .

Parameters that are convenient for the specification of the configuration of a system are selected to be generalized coordinates. If these parameters are independent of one another, then number of independent generalized coordinates is defined by the number of degrees of freedom of the system.^[2] ^[3]

The generalized velocities are the time derivatives of the generalized coordinates of the system.

Simple pendulum

The relationship between the use of generalized coordinates and Cartesian coordinates to characterize the movement of a mechanical system can be illustrated by considering the simple pendulum.

Coordinates

A simple pendulum consists of a mass M hanging from a pivot point so that it is constrained to move on a circle of radius L. The position of the mass is defined by the coordinate vector r=(x, y} measured in the plane of the circle such that y is in the vertical direction. The coordinates x and y are related by the equation of the circle

f (x, y) = x^{2} + y^{2} - L^{2} = 0,

that constrains the movement of M. This equation also provides a constraint on the velocity components,

\frac{d f}{d t} (x, y) = x \dot{x} + y \dot{y} = 0.

Now introduce the parameter θ, that defines the angular position of M from the vertical direction. It can be used to define the coordinates x and y, such that

𝐫 = (x, y) = (L \sin θ, - L \cos θ) .

The use of θ to define the configuration of this system avoids the constraint provided by the equation of the circle.

Virtual work

Notice that the force of gravity acting on the mass M is formulated in the usual Cartesian coordinates,

𝐅 = (0, - M g),

where g is the acceleration of gravity.

The virtual work of gravity on the mass M as it follows the trajectory r is given by

δ W = 𝐅 \cdot δ 𝐫 .

The variation δr can be computed in terms of the coordinates x and y, or in terms of the parameter θ,

δ 𝐫 = (δ x, δ y) = (L \cos θ, L \sin θ) δ θ .

Thus, the virtual work is given by

δ W = - M g δ y = - M g L \sin θ δ θ .

Notice that the coefficient of δy is the y-component of the applied force. In the same way, the coefficient of δθ is viewed as the generalized force along generalized coordinate θ, given by

F_{θ} = - M g L \sin θ .

Kinetic energy

To complete the analysis consider the kinetic energy T of the mass, using the velocity,

𝐯 = (\dot{x}, \dot{y}) = (L \cos θ, L \sin θ) \dot{θ},

so,

T = \frac{1}{2} M 𝐯 \cdot 𝐯 = \frac{1}{2} M ({\dot{x}}^{2} + {\dot{y}}^{2}) = \frac{1}{2} M L^{2} {\dot{θ}}^{2} .

Lagrange's equations

Lagrange's equations for the pendulum in terms of the coordinates x and y are given by,

\frac{d}{d t} \frac{\partial T}{\partial {\dot{q}}_{i}} - \frac{\partial T}{\partial q_{i}} = F_{q_{i}} + λ \frac{\partial f}{\partial q_{i}} .

This yields the three equations

M \ddot{x} = λ (2 x), M \ddot{y} = - M g + λ (2 y), x^{2} + y^{2} - L^{2} = 0,

in the three unknowns, x, y and λ.

Compute two derivatives of the constraint equation and assemble the equations into matrix form

[\begin{matrix} M & 0 & - 2 x \\ 0 & M & - 2 y \\ - 2 x & - 2 y & 0 \end{matrix}] {\begin{matrix} \ddot{x} \\ \ddot{y} \\ λ \end{matrix}} = {\begin{matrix} 0 \\ - M g \\ 2 {\dot{x}}^{2} + 2 {\dot{y}}^{2} \end{matrix}}

Using the parameter θ, Lagrange's equations take the form

M L^{2} \ddot{θ} = - M g L \sin θ,

or

\ddot{θ} + \frac{g}{L} \sin θ = 0.

This shows that the parameter θ is a generalized coordinate that can be used in the same way as the Cartesian coordinates x and y to analyze the pendulum.

Double pendulum

The benefits of generalized coordinates become apparent with the analysis of a double pendulum. For the two masses m_i, i=1, 2, let r=(x_i, y_i), i=1, 2 define their two trajectories. These vectors satisfy the two constraint equations,

f_{1} (x_{1}, y_{1}, x_{2}, y_{2}) = 𝐫_{1} \cdot 𝐫_{1} - L_{1}^{2} = 0, f_{2} (x_{1}, y_{1}, x_{2}, y_{2}) = (𝐫_{2} - 𝐫_{1}) \cdot (𝐫_{2} - 𝐫_{1}) - L_{2}^{2} = 0.

The formulation of Lagrange's equations for this system yields six equations in the four Cartesian coordinates x_i, y_i i=1, 2 and the two Lagrange multipliers λ_i, i=1, 2 that arise from the two constraint equations.

Coordinates

Now introduce the generalized coordinates θ_i i=1,2 that define the angular position of each mass of the double pendulum from the vertical direction. In this case, we have

𝐫_{1} = (L_{1} \sin θ_{1}, - L_{1} \cos θ_{1}), 𝐫_{2} = (L_{1} \sin θ_{1}, - L_{1} \cos θ_{1}) + (L_{2} \sin θ_{2}, - L_{2} \cos θ_{2}) .

The force of gravity acting on the masses is given by,

𝐅_{1} = (0, - m_{1} g), 𝐅_{2} = (0, - m_{2} g)

where g is the acceleration of gravity. Therefore, the virtual work of gravity on the two masses as they follow the trajectories r_i, i=1,2 is given by

δ W = 𝐅_{1} \cdot δ 𝐫_{1} + 𝐅_{2} \cdot δ 𝐫_{2} .

The variations δr_i i=1, 2 can be computed to be

δ 𝐫_{1} = (L_{1} \cos θ_{1}, L_{1} \sin θ_{1}) δ θ_{1}, δ 𝐫_{2} = (L_{1} \cos θ_{1}, L_{1} \sin θ_{1}) δ θ_{1} + (L_{2} \cos θ_{2}, L_{2} \sin θ_{2}) δ θ_{2}

Virtual work

Thus, the virtual work is given by

δ W = - (m_{1} + m_{2}) g L_{1} \sin θ_{1} δ θ_{1} - m_{2} g L_{2} \sin θ_{2} δ θ_{2},

and the generalized forces are

F_{θ_{1}} = - (m_{1} + m_{2}) g L \sin θ_{1}, F_{θ_{2}} = - m_{2} g L \sin θ_{2} .

Kinetic energy

Compute the kinetic energy of this system to be

T = \frac{1}{2} m_{1} 𝐯_{1} \cdot 𝐯_{1} + \frac{1}{2} m_{2} 𝐯_{2} \cdot 𝐯_{2} = \frac{1}{2} (m_{1} + m_{2}) {\dot{θ}}_{1}^{2} + \frac{1}{2} m_{2} {\dot{θ}}_{2}^{2} + m_{2} L_{1} L_{2} \cos (θ_{2} - θ_{1}) {\dot{θ}}_{1} {\dot{θ}}_{2} .

Lagrange's equations

Lagrange's equations yield two equations in the unknown generalized coordinates θ_i i=1, 2, given by

(m_{1} + m_{2}) L_{1} {\ddot{θ}}_{1} + m_{2} L_{1} L_{2} {\ddot{θ}}_{2} \cos (θ_{2} - θ_{1}) + m_{2} L_{1} L_{2} \sin (θ_{2} - θ_{1}) = - (m_{1} + m_{2}) g L \sin θ_{1},

and

m_{2} L_{2} {\ddot{θ}}^{2} + m_{2} L_{1} L_{2} {\ddot{θ}}_{1} \cos (θ_{2} - θ_{1}) - m_{2} L_{1} L_{2} \sin (θ_{2} - θ_{1}) = - m_{2} g L \sin θ_{2} .

The use of the generalized coordinates θ_i i=1, 2 provides an alternative to the Cartesian formulation of the dynamics of the double pendulum.

Constraint equations

Generalized coordinates may be independent (or unconstrained), in which case they are equal in number to the degrees of freedom of the system, or they may be dependent (or constrained), related by constraints on and among the coordinates. The number of dependent coordinates is the sum of the number of degrees of freedom and the number of constraints. For example, the constraints might take the form of a set of configuration constraint equations:^[1]

f_{i} ({q_{n}}, t) = 0,

where q_n is the n-th generalized coordinate and i denotes one of a set of constraint equations, taken here to vary with time t. The constraint equations limit the values available to the set of q_n, and thereby exclude certain configurations of the system.

It can be advantageous to choose independent generalized coordinates, as is done in Lagrangian mechanics, because this eliminates the need for constraint equations. However, in some situations, it is not possible to identify an unconstrained set. For example, when dealing with nonholonomic constraints or when trying to find the force due to any constraint—holonomic or not, dependent generalized coordinates must be employed. Sometimes independent generalized coordinates are called internal coordinates because they are mutually independent, otherwise unconstrained, and together give the position of the system.

A system with $m$ degrees of freedom and n particles whose positions are designated with three dimensional vectors, ${𝐫_{i}}$ , implies the existence of $3 n - m$ scalar constraint equations on those position variables. Such a system can be fully described by the scalar generalized coordinates, ${q_{1}, q_{2}, ..., q_{m}}$ , and the time, $t$ , if and only if all $m$ ${q_{j}}$ are independent coordinates. For the system, the transformation from old coordinates to generalized coordinates may be represented as follows:^[4]^: 260

𝐫_{1} = 𝐫_{1} (q_{1}, q_{2}, ..., q_{m}, t)

,

𝐫_{2} = 𝐫_{2} (q_{1}, q_{2}, ..., q_{m}, t)

, ...

𝐫_{n} = 𝐫_{n} (q_{1}, q_{2}, ..., q_{m}, t)

.

This transformation affords the flexibility in dealing with complex systems to use the most convenient and not necessarily inertial coordinates. These equations are used to construct differentials when considering virtual displacements and generalized forces.

Examples

Double pendulum

A double pendulum constrained to move in a plane may be described by the four Cartesian coordinates {x₁, y₁, x₂, y₂}, but the system only has two degrees of freedom, and a more efficient system would be to use

{q_{1}, q_{2}} = {θ_{1}, θ_{2}}

,

which are defined via the following relations:

{x_{1}, y_{1}} = {L_{1} \sin θ_{1}, L_{1} \cos θ_{1}}

{x_{2}, y_{2}} = {L_{1} \sin θ_{1} + L_{2} \sin θ_{2}, L_{1} \cos θ_{1} + L_{2} \cos θ_{2}}

Bead on a wire

A bead constrained to move on a wire has only one degree of freedom, and the generalized coordinate used to describe its motion is often

q_{1} = l

,

where l is the distance along the wire from some reference point on the wire. Notice that a motion embedded in three dimensions has been reduced to only one dimension.

Motion on a surface

A point mass constrained to a surface has two degrees of freedom, even though its motion is embedded in three dimensions. If the surface is a sphere, a good choice of coordinates would be:

{q_{1}, q_{2}} = {θ, ϕ}

,

where θ and φ are the angle coordinates familiar from spherical coordinates. The r coordinate has been effectively dropped, as a particle moving on a sphere maintains a constant radius.

Generalized velocities and kinetic energy

Each generalized coordinate $q_{i}$ is associated with a generalized velocity ${\dot{q}}_{i}$ , defined as:

{\dot{q}}_{i} = \frac{d q_{i}}{d t}

The kinetic energy of a particle is

T = \frac{m}{2} ({\dot{x}}^{2} + {\dot{y}}^{2} + {\dot{z}}^{2})

.

In more general terms, for a system of $p$ particles with $n$ degrees of freedom, this may be written

T = \sum_{i = 1}^{p} \frac{m_{i}}{2} ({\dot{x}}_{i}^{2} + {\dot{y}}_{i}^{2} + {\dot{z}}_{i}^{2})

.

If the transformation equations between the Cartesian and generalized coordinates

x_{i} = x_{i} (q_{1}, q_{2}, ..., q_{n}, t)

y_{i} = y_{i} (q_{1}, q_{2}, ..., q_{n}, t)

z_{i} = z_{i} (q_{1}, q_{2}, ..., q_{n}, t)

are known, then these equations may be differentiated to provide the time-derivatives to use in the above kinetic energy equation:

{\dot{x}}_{i} = \frac{d}{d t} x_{i} (q_{1}, q_{2}, ..., q_{n}, t) .

It is important to remember that the kinetic energy must be measured relative to inertial coordinates. If the above method is used, it means only that the Cartesian coordinates need to be inertial, even though the generalized coordinates need not be. This is another considerable convenience of the use of generalized coordinates.

Applications of generalized coordinates

Such coordinates are helpful principally in Lagrangian mechanics, where the forms of the principal equations describing the motion of the system are unchanged by a shift to generalized coordinates from any other coordinate system. The amount of virtual work done along any coordinate $q_{i}$ is given by:

δ W_{q_{i}} = F_{q_{i}} \cdot δ q_{i}

,

where $F_{q_{i}}$ is the generalized force in the $q_{i}$ direction. While the generalized force is difficult to construct 'a priori', it may be quickly derived by determining the amount of work that would be done by all non-constraint forces if the system underwent a virtual displacement of $δ q_{i}$ , with all other generalized coordinates and time held fixed. This will take the form:

δ W_{q_{i}} = f (q_{1}, q_{2}, ..., q_{n}) \cdot δ q_{i}

,

and the generalized force may then be calculated:

F_{q_{i}} = \frac{δ W_{q_{i}}}{δ q_{i}} = f (q_{1}, q_{2}, ..., q_{n})

.

Generalized coordinates and virtual work

The principle of virtual work states that the virtual work of the applied forces is zero for all virtual movements of the system from static equilibrium, that is, δW=0 for any variation δr.^[4] When formulated in terms of generalized coordinates, this is equivalent to the requirement that the generalized forces for any virtual displacement are zero, that is F_i=0.

Let the forces on the system be F_j, j=1, ..., m and let the virtual displacement of each point of application of these forces be δr_j, j=1, ..., m, then the virtual work generated by a virtual displacement of these forces from the equilibrium position is given by

δ W = \sum_{j = 1}^{m} 𝐅_{j} \cdot δ 𝐫_{j} .

Now assume that each δr_j depends on the generalized coordinates q_i, i=1, ..., n, then

δ 𝐫_{j} = \frac{\partial 𝐫_{j}}{\partial q_{1}} δ q_{1} + \dots + \frac{\partial 𝐫_{j}}{\partial q_{n}} δ q_{n},

and

δ W = (\sum_{j = 1}^{m} 𝐅_{j} \cdot \frac{\partial 𝐫_{j}}{\partial q_{1}}) δ q_{1} + \dots + (\sum_{j = 1}^{m} 𝐅_{j} \cdot \frac{\partial 𝐫_{j}}{\partial q_{n}}) δ q_{n} .

The n terms

F_{i} = \sum_{j = 1}^{m} 𝐅_{j} \cdot \frac{\partial 𝐫_{j}}{\partial q_{i}}, i = 1, \dots, n,

are the generalized forces acting on the system. Kane^[5] shows that these generalized forces can also be formulated in terms of the ratio of time derivatives,

F_{i} = \sum_{j = 1}^{m} 𝐅_{j} \cdot \frac{\partial 𝐯_{j}}{\partial {\dot{q}}_{i}}, i = 1, \dots, n,

where v_j is the velocity of the point of application of the force F_j.

In order for the virtual work to be zero for an arbitrary virtual displacement, each of the generalized forces must be zero, that is

δ W = 0 \Rightarrow F_{i} = 0, i = 1, \dots, n .

References

1 2 Jerry H. Ginsberg (2008). "§7.2.1 Selection of generalized coordinates". Engineering dynamics, Volume 10 (3rd ed.). Cambridge University Press. p. 397. ISBN 0-521-88303-2.
↑ Farid M. L. Amirouche (2006). "§2.4: Generalized coordinates". Fundamentals of multibody dynamics: theory and applications. Springer. p. 46. ISBN 0-8176-4236-6.
↑ Florian Scheck (2010). "§5.1 Manifolds of generalized coordinates". Mechanics: From Newton's Laws to Deterministic Chaos (5th ed.). Springer. p. 286. ISBN 3-642-05369-6.
1 2 Torby, Bruce (1984). "Energy Methods". Advanced Dynamics for Engineers. HRW Series in Mechanical Engineering. United States of America: CBS College Publishing. ISBN 0-03-063366-4.
↑ T. R. Kane and D. A. Levinson, Dynamics: theory and applications, McGraw-Hill, New York, 1985

Greenwood, Donald T. (1987). Principles of Dynamics (2nd edition ed.). Prentice Hall. ISBN 0-13-709981-9. {{cite book}}: |edition= has extra text (help)
Wells, D. A. (1967). Schaum's Outline of Lagrangian Dynamics. New York: McGraw-Hill.

[Ginsberg-1] 1 2 Jerry H. Ginsberg (2008). "§7.2.1 Selection of generalized coordinates". Engineering dynamics, Volume 10 (3rd ed.). Cambridge University Press. p. 397. ISBN 0-521-88303-2.

[Amirouche-2] Farid M. L. Amirouche (2006). "§2.4: Generalized coordinates". Fundamentals of multibody dynamics: theory and applications. Springer. p. 46. ISBN 0-8176-4236-6.

[Scheck-3] Florian Scheck (2010). "§5.1 Manifolds of generalized coordinates". Mechanics: From Newton's Laws to Deterministic Chaos (5th ed.). Springer. p. 286. ISBN 3-642-05369-6.

[Torby1984-4] 1 2 Torby, Bruce (1984). "Energy Methods". Advanced Dynamics for Engineers. HRW Series in Mechanical Engineering. United States of America: CBS College Publishing. ISBN 0-03-063366-4.

[5] T. R. Kane and D. A. Levinson, Dynamics: theory and applications, McGraw-Hill, New York, 1985

[1]

[2]

[3]

[4]

[5]

Simple pendulum

Coordinates

Virtual work

Kinetic energy

Lagrange's equations

Double pendulum

Coordinates

Virtual work

Kinetic energy

Lagrange's equations

Constraint equations

Examples

Double pendulum

Bead on a wire

Motion on a surface

Generalized velocities and kinetic energy

Applications of generalized coordinates

Generalized coordinates and virtual work

See also

References