1
PHYS 705: Classical Mechanics Constraints and Generalized Coordinates 2
Constraints
Problem Statement: In solving mechanical problems, we start with the 2nd law
(e ) (*) FFji i m ii r j
th In principle, one can solve for ri(t) (trajectory) for the i particle by specifying all the external and internal forces acting on it .
However, if constraints are present, these external forces in general are NOT known.
Therefore, we need to understand the various constraints and know how to handle them. 3
Holonomic Constraints
Holonomic constraints can be expressed as a function in terms of the coordinates and time,
rij fr1, r 2 , ; t 0 ri rj 2 2 e.g. (a rigid body) r r c 0 i j ij O
non-holonomic examples: -Gas in a container -Object rolling on a rough surface without slipping… more later
More quantifiers: - Rheonomous: depend on time explicitly - Scleronomous: not explicitly depend on time
e.g. a bead constraints to move on a fixed vs. a moving wire 4
Constraints and Generalized Coordinates
Difficulties involving constraints:
1. Through f r 1 , r 2 , ; t 0 , the individual coordinates ri are no longer independent
eqs of motion (*) for individual particles are now coupled (not independent)
2. Forces of constraints are not known a priori and must be solved as additional unknowns
With holonomic constraints:
Prob #1 can be handled by introducing a set of “proper” (independent) Generalized Coordinates
Prob #2 can be treated with: D’Alembert’s Principle & Lagrange’s Equations (with Lagrange multipliers) 5
Generalized Coordinates
• Without constraints, a system of N particles has 3N dof • With K constraint equations, the # dof reduces to 3N-K
• With holonomic constraints, one can introduce (3N-K) independent
(proper) generalized coordinates q1,,, q 2 q 3N K such that:
r1 r 1qq 1,,,, 2 q 3N K t riandq1 , q 3 NK are related by a point transformation rNN r qq1,,,, 2 q 3 NK t
• The goal is to describe the time evolution of the system in the set of (3N-K) independent (proper) generalized coordinates. 6
Generalized Coordinates
Generalized coordinates can be anything: angles, energy units, momentum
units, or even amplitudes in the Fourier expansion of ri
Each q j is just a number, a scalar
But, they must completely specify the state of a given system
The choice of a particular set of generalized coordinates is not unique.
No specific rule in finding the most “suitable” (resulting in simplest EOM) 7
Generalized Coordinates
Example: In Cartesian coordri :
xyx,,, y we have 4 dof O x 1 1 2 2 2 2 2 l x y l 0 1 2 constraints: 1 1 (x , y ) 2 2 1 1 xx yy l 2 0 m1 2 1 2 1 2 l So, there are only 2 indep dof…
m2 One choice of generalized coordsq j is: y (x2, y2) 1, 2 2 indep dof (Double Plane Pendulum) And, they are linked to the Cart. coord through: 1 (constraints are implicitly 1 tan x 1 y 1 1 encoded in the Pt Trans) 2tan xx 2 1 yy 2 1 8
Non-Holonomic Constraints - can’t use constraint equations to eliminate dependent coordinates - in general, solution is problem specific.
Example in book: Vertical Disk rolling without slipping on a horizontal plane
z y a f Described by 4 coordinates: (x, y) of the contact point v : orientation of disk- angle of disk axis with x-axis f : angle of rotation of the disk O x 9
Non-Holonomic Constraints (rolling disk exp)
Now, consider the constraints:
1. No-slip condition : s af v a f y top view 2. Disk rolling vertically v disk axis see graph v sin
v -v cos x vsin y v cos x 10
Non-Holonomic Constraints (rolling disk exp)
Putting them together, gives the following differential equations of constraint, dx df afsin a sin dt dt dx asin d f 0 or dy df dy acos f d 0 afcos a cos dt dt
The point is that we can’t write this in Holonomic form:
ff xy, , 0 with f being a function! (hw)
Physical intuition Roll the disk in a circle with radius R. Upon completion of the circle, x, y and will have returned to their original values but, f will depend on R (can’t be specified by x ,, y ) 11
How to deal with Constraints? Principle of Virtual Work
Consider the simplest situation, a system in equilibrium first,
(note: i labels the particles) - The net force on each particle vanishes: Fi 0
Consider an arbitrary “virtual” infinitesimal change in the coordinates, ri
- Virtual means that it is done with no change in time during which forces and constraints do not change.
- These virtual displacements are done consistent with the constraints (we will be more specific later). F r 0 Since the net force on each particle, Fi is i i i zero (equilibrium), obviously we have: (virtual work) 12
Principle of Virtual Work
Separating the forces into applied ( a ) and constraint forces , Fi fi
(a ) FFi i f i
(a ) Then, we have Fi r i fr ii 0 i i
Now, what do we mean by the virtual displacements r i being done consistent with the constraint? 13
Principle of Virtual Work
For virtual displacements to be consistent with the constraints means that
the virtual work done by the constraint forces along the virtual displacement must be zero.
Geometric view fi ri >> More on this later for N particles with K constraints, motion is restricted f r or fr 0 on a (3N-K)-D surface and the constraint forces fi i i ii must be to that surface.
For r i to be consistent with constraints means that the net virtual work from the forces of constraints is zero !
fi r i 0 i 14
Principle of Virtual Work
Back to our original equation for a constrained system in equilibrium, 0 (a ) Fi r i fr ii 0 i i
With the virtual displacements satisfying the constraints leaves us with the statement, (a ) Fi r i 0 i
The virtual work of the applied forces must also vanish!
This is called the Principle of Virtual Work. 15
Principle of Virtual Work
(a ) Fi r i 0 i
Note: Since the coordinates (and the virtual variations) are not necessarily independent. They are linked through the constraint equations. The Principle of Virtual Work does not imply,
(a ) Fi 0 for all i independently.
The trick is now to change variables to a set of proper (independent) generalized coordinates. Then, we can rewrite the equation as, ? q 0 j j j With q being independent, we can then claim:? 0 for all j. As we will j j see, this will give us expressions which will lead to the solution of the problem. 16
D’Alembert’s Principle
Now, we consider the more general case when the system is not necessary in equilibrium so that the net force on the particles is NOT zero. We continue to assume the constraints forces to be unknown a priori…
Similar to our discussion on the Principle of Virtual Work, we would like to reformulate the mechanical problem to include the constraint forces such that they “disappear” you solve the “new” problem using only the (given) applied forces.
This is the basis for the D’Alembert’s Principle AND by additionally choosing a set of proper generalized coordinates, the problem can be solved and it will result in the Euler-Lagrange’s Equations. 17
D’Alembert’s Principle
In deriving the Principle of Virtual Work, the system was in equilibrium.
In extending it to include dynamics , we will begin with Newton’s 2nd law,
th Fpiior Fp ii 0 for the i particle in the system.
We again consider a virtual infinitesimal displacement r i consistent with the given constraint. Since we have for all particles, Fi p i 0
We have, Fpi i r i 0 i (a ) Again, we separate out the applied and constraint forces, FFi i f i
(a ) This gives, Fi pr ii fr ii 0 i i 18
D’Alembert’s Principle
Then, following a similar argument for the virtual displacement to be consistent with constraints, i.e, (no virtual work for f ) fi r i 0 i i (a ) We can write down, Fi p i r i 0 i This is the D’Alembert’s Principle.
Again, since the coordinates (and the virtual variations) are not necessary (a ) independent. This does NOT implies, F i p i 0 for the individual i.
We then need to look into changing variables to a set of independent generalized ? q 0 coordinates so that we have j j with the coefficients j in the sum independently equal to zero, i.e., ? 0 j 19
Geometric View of the D’Alembert’s Principle
Consider a particle moving in 3D with one Holonomic constraint,
z (a ) equation of motion: mrF f r (,x yz ,) equation of constraint: g(r , t ) 0
Here, m - F(a) is the known applied force r (t) - And, we model the unknown y g(r , t ) 0 constraint force by the vector f. x Note: r(t) has 3 unknown trajectory (red) is constraint to move in a 3- components + 1 constraint 1=2 dimensional surface (blue g ( r , t ) 0 ). 20
Geometric View of the D’Alembert’s Principle
- There are three unknown components to the constraint force f. A scalar constraint does not specify the vector f completely. - There are multiple choices for f which satisfy g(r, t)=0 BUT there is an additional physical restriction on f that we should consider…
z z
f f
y y
x x Observation: For a given f, adding a component // to the surface will still keep the particle on the surface (satisfying g(r, t)=0) but will result with an additional acceleration along the surface). 21
Geometric View of the D’Alembert’s Principle
A reasonable physical argument is to restrict the choice of f so that:
Constraint Force f needs to lay to the constraint surface
Note that g(r , t ) 0 is the equation for the z constraint surface and
g(r , t ) surface f
So, we can “parametrized” f in term of g (,) r t , y g(r , t ) 0 f g(,) r t where is a parameter x
(a ) mxF gt(,) r 4 unknowns r and This gives, 4 equations g(r , t ) 0 22
Geometric View of the D’Alembert’s Principle
(a ) mrF gt(,) r 4 unknowns r and g(r , t ) 0 4 equations This system is solvable but now we would like to solve the system w/o using the constraint explicitly …
Note that g is to the surface of constraint and we can project the dynamical equation onto the tangent plane of the constraint surface at (r, t).
To do that, take e a and e b as two basis vectors spanning the tangent plane to the constraint surface at (r, t). Dotting the above Eq to e a and e b gives two independent scalar equations, (a ) mrFe a gt( re , ) a 0 (a ) mrFe b gt( re , ) b 0 23
Geometric View of the D’Alembert’s Principle
Together with the constraint equation itself, we then have 3 eqs for the 3 unknown components of r.
(a ) mrF ea, b 0 3 unknowns r 3 equations g(r , t ) 0
So, now, in principle, we can solve for the dynamical equation (EOM), r(t), without knowing the constraint forces f explicitly.
(a ) mrF ea, b 0
This is the D’Alembert’s Principle (for a single particle). 24
Geometric View of the D’Alembert’s Principle
We can generalize the argument to a system of N particles with K constraints (Holonomic):
(a ) Note: The virtual ri miirF i e k 0 displacements consistent i with the constraints are pF(a ) r 0 in the tangent space ii i e i spanned by the basis k
Geometric Interpretation: The K constraints restrict the system to a (3N-K)-D surface within the 3N-D
space. There are (3N-K) ek vectors spanning that tangent plane to the constraint surface so that the above expression gives (3N-K) equations that the problem can be solved without knowing the constraint forces explicitly. 25
D’Alembert’s Principle
(a ) Fi p i r i 0 i
To solve for EOM using the D’Alembert’s Principle …
We still need to look into changing variables to a set of independent generalized coordinates so that we have ? q 0 j j j Then, we can claim the “coefficients” ? in the sum to be independently j equal to zero, i.e., ? 0 j 26
D’Alembert’s Principle
(a ) Fi p i r i 0 i
To solve for EOM using the D’Alembert’s Principle …
We still need to look into changing variables to a set of independent generalized coordinates so that we have ? q 0 j j j The correct “coefficients” allowing us to have ? 0 will give us the Euler- j Lagrange equation and the EL Eq gives an explicit expression for the EOM:
dT T Q 0 j dt q j q j 27
Side Note: Constraint and Work
(a ) Recall that we have from the EOM: mrF gt(,) r
Let F(a) be a conservative force , i.e., F ( a ) U (,) r t so that
mr U g Dotting r into both sides,
d1 2 dT mrr m r U r g r dt2 dt
Consider the full-time derivative of g , we have,
dg g dx g dy g dz g g g r dt x dt y dt z dt t t 28
An Aside: Constraint and Work
As the particle moves, it is constraint to stay on the g=0 surface, d g g So, 0 and, g r dt t dU U Similarly, considering the full-time derivate of U, U r dt t Putting everything together,
mr r U r g r With E=T+U,
dT dU U g dE U g dt dt t t dt t t 29
An Aside: Constraint and Work dE U g dt t t So, either U or g explicitly depends on time, the total energy changes with time.
Since we typically do not consider time-dependent U potential functions, So, we can make the following assertions:
Scleronomous (g not explicitly depends on t) Holonomic Constraints: g g r 0 and constraint force won’t do work! t
Rheonomous (g explicitly depends on t) Holonomic Constraints: g g r 0 and constraint force can do work! t