World Library  
Flag as Inappropriate
Email this Article

Generalized coordinates

Article Id: WHEBN0000473514
Reproduction Date:

Title: Generalized coordinates  
Author: World Heritage Encyclopedia
Language: English
Subject: Lagrangian Mechanics, Virtual work, Analytical mechanics, Equations of motion, Maupertuis' principle
Collection: Dynamical Systems, Lagrangian Mechanics, Rigid Bodies
Publisher: World Heritage Encyclopedia
Publication
Date:
 

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] The generalized velocities are the time derivatives of the generalized coordinates of the system.

An example of a generalized coordinate is the angle that locates a point moving on a circle. The adjective "generalized" distinguishes these parameters from the traditional use of the term coordinate to refer to Cartesian coordinates: for example, describing the location of the point on the circle using x and y coordinates.

Although there may be many choices for generalized coordinates for a physical system, parameters which are convenient are usually selected for the specification of the configuration of the system and which make the solution of its equations of motion easier. If these parameters are independent of one another, the number of independent generalized coordinates is defined by the number of degrees of freedom of the system.[2] [3]

Contents

  • Constraints and degrees of freedom 1
    • Holonomic constraints 1.1
    • Non-holonomic constraints 1.2
  • Physical quantities in generalized coordinates 2
    • Kinetic energy 2.1
    • Generalized momentum 2.2
  • Examples 3
    • Simple pendulum 3.1
    • Double pendulum 3.2
  • Generalized coordinates and virtual work 4
  • See also 5
  • Notes 6
  • References 7
  • Bibliography of cited references 8

Constraints and degrees of freedom

Open curved path F(x, y) = 0
Closed curved path C(x, y) = 0
One generalized coordinate, one degree of freedom, on paths in 2d. Only one number is needed to uniquely specify positions on the curve, the examples shown are the arc length s or angle θ. Both of the Cartesian coordinates (x, y) are unnecessary since either x or y is related to the other by the equations of the curves. They can also be parameterized by s or θ.
Open curved path F(x, y) = 0
Closed curved path C(x, y) = 0
The arc length s along the curve is a legitimate generalized coordinate since the position is uniquely determined, but the angle θ is not since there are multiple positions owing to the self-intersections of the curves.

Generalized coordinates are usually selected to provide the minimum number of independent coordinates that define the configuration of a system, which simplifies the formulation of Lagrange's equations of motion. However, it can also occur that a useful set of generalized coordinates may be dependent, which means that they are related by one or more constraint equations.

Holonomic constraints

Top: one degree of freedom, bottom: two degrees of freedom, left: an open curve F (parameterized by t) and surface F, right: a closed curve C and closed surface S. The equations shown are the constraint equations. Generalized coordinates are chosen and defined with respect to these curves (one per degree of freedom), and simplify the analysis since even complicated curves are described by the minimum number of coordinates required.

For a system of N particles in 3d real coordinate space, the position vector of each particle can be written as a 3-tuple in Cartesian coordinates;

\mathbf{r}_1 = (x_1,y_1,z_1) \,, \quad \mathbf{r}_2 = (x_2,y_2,z_2) \,, \ldots \,, \mathbf{r}_N = (x_N,y_N,z_N)\,.

Any of the position vectors can be denoted rk where k = 1, 2, ..., N labels the particles. A holonomic constraint is a constraint equation of the form for particle k[4][nb 1]

f(\mathbf{r}_k, t) = 0

which connects all the 3 spatial coordinates of that particle together, so they are not independent. The constraint may change with time, so time t will appear explicitly in the constraint equations. At any instant of time, when t is a constant, any one coordinate will be determined from the other coordinates, e.g. if xk and zk are given, then so is yk. One constraint equation counts as one constraint. If there are C constraints, each has an equation, so there will be C constraint equations. There is not necessarily one constraint equation for each particle, and if there are no constraints on the system then there are no constraint equations.

So far, the configuration of the system is defined by 3N quantities, but C coordinates can be eliminated, one coordinate from each constraint equation. The number of independent coordinates is n = 3NC. (In D dimensions, the original configuration would need ND coordinates, and the reduction by constraints means n = NDC). It is ideal to use the minimum number of coordinates needed to define the configuration of the entire system, while taking advantage of the constraints on the system. These quantities are known as generalized coordinates in this context, denoted qj(t). It is convenient to collect them into an n-tuple

\mathbf{q}(t) = (q_1(t), q_2(t), \ldots, q_n(t))

which is a point in the configuration space of the system. They are all independent of one other, and each is a function of time. Geometrically they can be lengths along straight lines, or arc lengths along curves, or angles; not necessarily Cartesian coordinates or other standard orthogonal coordinates. There is one for each degree of freedom, so the number of generalized coordinates equals the number of degrees of freedom, n. A degree of freedom corresponds to one quantity that changes the configuration of the system, for example the angle of a pendulum, or the arc length traversed by a bead along a wire.

If it is possible to find from the constraints as many independent variables as there are degrees of freedom, these can be used as generalized coordinates[5] The position vector rk of particle k is a function of all the n generalized coordinates and time,[6][7][8][5][nb 2]

\mathbf{r}_k = \mathbf{r}_k(\mathbf{q}(t),t) \,,

and the generalized coordinates can be thought of as parameters associated with the constraint.

The corresponding time derivatives of q are the generalized velocities,

\dot{\mathbf{q}} = \frac{d\mathbf{q}}{dt} = (\dot{q}_1(t), \dot{q}_2(t), \ldots, \dot{q}_n(t))

(each dot over a quantity indicates one time derivative). The velocity vector vk is the total derivative of rk with respect to time

\mathbf{v}_k = \dot{\mathbf{r}}_k = \frac{d\mathbf{r}_k}{dt} = \sum_{j=1}^n \frac{\partial \mathbf{r}_k}{\partial q_j}\dot{q}_j +\frac{\partial \mathbf{r}_k}{\partial t}\,.

and so generally depends on the generalized velocities and coordinates. Since we are free to specify the initial values of the generalized coordinates and velocities separately, the generalized coordinates and velocities can be treated as independent variables. The generalized coordinates qj and velocities dqj/dt are treated as independent variables.

Non-holonomic constraints

A mechanical system can involve constraints on both the generalized coordinates and their derivatives. Constraints of this type are known as non-holonomic. First-order non-holonomic constraints have the form

g(\mathbf{q}, \dot{\mathbf{q}}, t) = 0\,,

An example of such a constraint is a rolling wheel or knife-edge that constrains the direction of the velocity vector. Non-holonomic constraints can also involve next-order derivatives such as generalized accelerations.

Physical quantities in generalized coordinates

Kinetic energy

The total kinetic energy of the system is the energy of the system's motion, defined as[9]

T = \frac {1}{2} \sum_{k=1}^N m_k \dot{\mathbf{r}}_k \cdot \dot{\mathbf{r}}_k\,,

in which · is the dot product. The kinetic energy is a function only of the velocities vk, not the coordinates rk themselves. By contrast an important observation is[10]

\dot{\mathbf{r}}_k \cdot \dot{\mathbf{r}}_k = \sum_{i,j=1}^n \left(\frac{\partial \mathbf{r}_k}{\partial q_i}\cdot\frac{\partial \mathbf{r}_k}{\partial q_j}\right)\dot{q}_i\dot{q}_j + \sum_{i=1}^n \left(2\frac{\partial \mathbf{r}_k}{\partial q_i}\cdot\frac{\partial \mathbf{r}_k}{\partial t}\right) \dot{q}_i + \left(\frac{\partial \mathbf{r}_k}{\partial t}\cdot\frac{\partial \mathbf{r}_k}{\partial t}\right) \,,

which illustrates the kinetic energy is in general a function of the generalized velocities, coordinates, and time if the constraint also varies with time, so T = T(q, dq/dt, t).

In the case the constraint on the particle is time-independent, then all partial derivatives with respect to time are zero, and the kinetic energy has no time-dependence and is a homogeneous function of degree 2 in the generalized velocities;

\dot{\mathbf{r}}_k\cdot \dot{\mathbf{r}}_k = \sum_{i,j=1}^n \left(\frac{\partial \mathbf{r}_k}{\partial q_i}\cdot\frac{\partial \mathbf{r}_k}{\partial q_j}\right) \dot{q}_i \dot{q}_j \,.

Still for the time-independent case, this expression is equivalent to taking the line element squared of the trajectory for particle k,

ds_k^2 = d\mathbf{r}_k\cdot d\mathbf{r}_k = \sum_{i,j=1}^n \left(\frac{\partial \mathbf{r}_k}{\partial q_i}\cdot\frac{\partial \mathbf{r}_k}{\partial q_j}\right) dq_i dq_j \,,

and dividing by the square differential in time, dt2, to obtain the velocity squared of particle k. Thus for time-independent constraints it is sufficient to know the line element to quickly obtain the kinetic energy of particles and hence the Lagrangian.[11]

It is instructive to see the various cases of polar coordinates in 2d and 3d, owing to their frequent appearance. In 2d polar coordinates (r, θ),

\left(\frac{ds}{dt}\right)^2 = \dot{r}^2 + r^2\dot{\theta}^2 \,,

in 3d cylindrical coordinates (r, θ, z),

\left(\frac{ds}{dt}\right)^2 = \dot{r}^2 + r^2\dot{\theta}^2 + \dot{z}^2 \,,

in 3d spherical coordinates (r, θ, φ),

\left(\frac{ds}{dt}\right)^2 = \dot{r}^2+r^2\dot{\theta}^2 +r^2\sin^2\theta \, \dot{\varphi}^2 \,.

Generalized momentum

The generalized momentum "canonically conjugate to" the coordinate qi is defined by

p_i =\frac{\partial L}{\partial\dot q_i}.

If the Lagrangian L does not depend on some coordinate qi, then it follows from the Euler–Lagrange equations that the corresponding generalized momentum will be a conserved quantity, because its time derivative is zero so the momentum must be a constant of the motion;

\dot{p}_i = \frac{d}{dt}\frac{\partial L}{\partial\dot q_i} = \frac{\partial L}{\partial q_i}=0\,.

Examples

Simple pendulum

Dynamic model of a 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 constrained dynamics of a simple pendulum.[12][13]

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,

\dot{f}(x, y)=2x\dot{x} + 2y\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

\mathbf{r}=(x, y) = (L\sin\theta, -L\cos\theta).

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

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

\mathbf{F}=(0,-mg),

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

\delta W = \mathbf{F}\cdot\delta \mathbf{r}.

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

\delta \mathbf{r} =(\delta x, \delta y) = (L\cos\theta, L\sin\theta)\delta\theta.

Thus, the virtual work is given by

\delta W = -mg\delta y = -mgL\sin\theta\delta\theta.

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

F_{\theta} = -mgL\sin\theta.

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

\mathbf{v}=(\dot{x}, \dot{y}) = (L\cos\theta, L\sin\theta)\dot{\theta},

so,

T= \frac{1}{2} m\mathbf{v}\cdot\mathbf{v} = \frac{1}{2} m (\dot{x}^2+\dot{y}^2) = \frac{1}{2} m L^2\dot{\theta}^2.

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

\frac{d}{dt}\frac{\partial T}{\partial \dot{x}} - \frac{\partial T}{\partial x} = F_{x} + \lambda \frac{\partial f}{\partial x},\quad \frac{d}{dt}\frac{\partial T}{\partial \dot{y}} - \frac{\partial T}{\partial y} = F_{y} + \lambda \frac{\partial f}{\partial y}.

This yields the three equations

m\ddot{x} = \lambda(2x),\quad m\ddot{y} = -mg + \lambda(2y),\quad x^2+y^2 - L^2=0,

in the three unknowns, x, y and λ.

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

\frac{d}{dt}\frac{\partial T}{\partial \dot{\theta}} - \frac{\partial T}{\partial \theta} = F_{\theta},

which becomes,

mL^2\ddot{\theta} = -mgL\sin\theta,

or

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

This formulation yields one equation because there is a single parameter and no constraint equation.

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

A double pendulum

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

f_1 (x_1, y_1, x_2, y_2) = \mathbf{r}_1\cdot \mathbf{r}_1 - L_1^2 = 0, \quad f_2 (x_1, y_1, x_2, y_2) = (\mathbf{r}_2-\mathbf{r}_1) \cdot (\mathbf{r}_2-\mathbf{r}_1) - L_2^2 = 0.

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

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

\mathbf{r}_1 = (L_1\sin\theta_1, -L_1\cos\theta_1), \quad \mathbf{r}_2 = (L_1\sin\theta_1, -L_1\cos\theta_1) + (L_2\sin\theta_2, -L_2\cos\theta_2).

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

\mathbf{F}_1=(0,-m_1 g),\quad \mathbf{F}_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 ri, i=1,2 is given by

\delta W = \mathbf{F}_1\cdot\delta \mathbf{r}_1 + \mathbf{F}_2\cdot\delta \mathbf{r}_2.

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

\delta \mathbf{r}_1 = (L_1\cos\theta_1, L_1\sin\theta_1)\delta\theta_1, \quad \delta \mathbf{r}_2 = (L_1\cos\theta_1, L_1\sin\theta_1)\delta\theta_1 +(L_2\cos\theta_2, L_2\sin\theta_2)\delta\theta_2

Thus, the virtual work is given by

\delta W = -(m_1+m_2)gL_1\sin\theta_1\delta\theta_1 - m_2gL_2\sin\theta_2\delta\theta_2,

and the generalized forces are

F_{\theta_1} = -(m_1+m_2)gL_1\sin\theta_1,\quad F_{\theta_2} = -m_2gL_2\sin\theta_2.

Compute the kinetic energy of this system to be

T= \frac{1}{2}m_1 \mathbf{v}_1\cdot\mathbf{v}_1 + \frac{1}{2}m_2 \mathbf{v}_2\cdot\mathbf{v}_2 = \frac{1}{2}(m_1+m_2)L_1^2\dot{\theta}_1^2 + \frac{1}{2}m_2L_2^2\dot{\theta}_2^2 + m_2L_1L_2 \cos(\theta_2-\theta_1)\dot{\theta}_1\dot{\theta}_2.

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

(m_1+m_2)L_1^2\ddot{\theta}_1+m_2L_1L_2\ddot{\theta}_2\cos(\theta_2-\theta_1) + m_2L_1L_2\ddot{\theta_2}^2\sin(\theta_1-\theta_2) = -(m_1+m_2)gL_1\sin\theta_1,

and

m_2L_2^2\ddot{\theta}_2+m_2L_1L_2\ddot{\theta}_1\cos(\theta_2-\theta_1) + m_2L_1L_2\ddot{\theta_1}^2\sin(\theta_2-\theta_1)=-m_2gL_2\sin\theta_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.

Generalized coordinates and virtual work

The principle of virtual work states that if a system is in static equilibrium, the virtual work of the applied forces is zero for all virtual movements of the system from this state, that is, δW=0 for any variation δr.[15] 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 Fi=0.

Let the forces on the system be Fj, j=1, ..., m be applied to points with Cartesian coordinates rj, j=1,..., m, then the virtual work generated by a virtual displacement from the equilibrium position is given by

\delta W = \sum_{j=1}^m \mathbf{F}_j\cdot \delta\mathbf{r}_j.

where δrj, j=1, ..., m denote the virtual displacements of each point in the body.

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

\delta \mathbf{r}_j = \frac{\partial \mathbf{r}_j}{\partial q_1} \delta{q}_1 + \ldots + \frac{\partial \mathbf{r}_j}{\partial q_n} \delta{q}_n,

and

\delta W = \left(\sum_{j=1}^m \mathbf{F}_j\cdot \frac{\partial \mathbf{r}_j}{\partial q_1}\right) \delta{q}_1 + \ldots + \left(\sum_{j=1}^m \mathbf{F}_j\cdot \frac{\partial \mathbf{r}_j}{\partial q_n}\right) \delta{q}_n.

The n terms

F_i = \sum_{j=1}^m \mathbf{F}_j\cdot \frac{\partial \mathbf{r}_j}{\partial q_i},\quad i=1,\ldots, n,

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

F_i = \sum_{j=1}^m \mathbf{F}_j\cdot \frac{\partial \mathbf{v}_j}{\partial \dot{q}_i},\quad i=1,\ldots, n,

where vj is the velocity of the point of application of the force Fj.

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

\delta W = 0 \quad \Rightarrow \quad F_i =0, i=1,\ldots, n.

See also

Notes

  1. ^ Some authors set the constraint equations to a constant for convenience with some constraint equations (e.g. pendulums), others set it to zero. It makes no difference because the constant can be subtracted to give zero on one side of the equation. Also, in Lagrange's equations of the first kind, only the derivatives are needed.
  2. ^ Some authors e.g. Hand & Finch take the form of the position vector for particle k, as shown here, as the condition for the constraint on that particle to be holonomic.

References

  1. ^ Ginsberg 2008, p. 397,  §7.2.1 Selection of generalized coordinates
  2. ^ Farid M. L. Amirouche (2006). "§2.4: Generalized coordinates". Fundamentals of multibody dynamics: theory and applications. Springer. p. 46.  
  3. ^ Florian Scheck (2010). "§5.1 Manifolds of generalized coordinates". Mechanics: From Newton's Laws to Deterministic Chaos (5th ed.). Springer. p. 286.  
  4. ^ Goldstein 1980, p. 12
  5. ^ a b Kibble & Berkshire 2004, p. 232
  6. ^ Torby 1984, p. 260
  7. ^ Goldstein 1980, p. 13
  8. ^ Hand & Finch 2008, p. 15
  9. ^ Torby 1984, p. 269
  10. ^ Goldstein 1980, p. 25
  11. ^ Landau & Lifshitz 1976, p. 8
  12. ^ Greenwood, Donald T. (1987). Principles of Dynamics (2nd ed.). Prentice Hall.  
  13. ^ Richard Fitzpatrick, Newtonian Dynamics, http://farside.ph.utexas.edu/teaching/336k/Newton/Newtonhtml.html.
  14. ^ Eric W. Weisstein, Double Pendulum, scienceworld.wolfram.com. 2007
  15. ^ Torby, Bruce (1984). "Energy Methods". Advanced Dynamics for Engineers. HRW Series in Mechanical Engineering. United States of America: CBS College Publishing.  
  16. ^ T. R. Kane and D. A. Levinson, Dynamics: theory and applications, McGraw-Hill, New York, 1985

Bibliography of cited references

  • Ginsberg, Jerry H. (2008). Engineering dynamics (3rd ed.). Cambridge UK:  
  •  
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
 
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
 
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.
 


Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.