dynamics

8. Lagrangian dynamics

  1. Lagrangian dynamics
    1. Degrees of freedom and phase space
    2. Lagrange equations
    3. Equilibrium conditions
    4. Dissipative Forces
    5. Forces of constraint

Introduction

ATHLETE Robot

Each arm in a robot has usually 3 joints, with two perpendicular axes allowing two independent rotations, corresponding to two degrees of freedom. Thus, there are 6 degrees of freedom for each arm, enough to be able to reach any point, in any desired direction, within its maximum range. The ATHLETE ( All-Terrain Hex-Legged Extra-Terrestrial Explorer ) robot, used by NASA for lunar exploration, has six 3-joint arms which means 36 degrees of freedom. Plus the 3 degrees of freedom to determine the position of a point in the body of the robot, that makes a total of 39 degrees of freedom. The human arm, not including the hand, has 7 degrees of freedom corresponding to three independent rotations at the shoulder, two rotations at the elbow and two other rotations at the wrist.

8.1. Degrees of freedom and phase space

All mechanical systems considered in the previous chapter have a single degree of freedom (a coordinate or an angle which determines the position) and two state variables: one variable for the position and its derivative with respect to time (linear or angular velocity).

The state of a system with n degrees of freedom is specified by n independent variables dependent on time, called generalized coordinates, which we will denoted by the letters: q 1 , q 2 , …, q n . These variables can be distances, angles or any other quantities. The time derivatives of each of these variables are the generalized velocities: ˙ q i .

The phase space of the system has 2 n dimensions. The 2 n coordinates of a point in phase space are ( q 1 , …, q n , ˙ q 1 , …, ˙ q n ) and the 2 n components of the phase velocity are ( ˙ q 1 , …, ˙ q n , ¨ q 1 , …, ¨ q n ). To determine the expression of the phase velocity at any point in phase space, it is necessary to know n expressions for the generalized accelerations ¨ q i , in terms of the generalized coordinates and velocities, which are called equations of motion.

The equations of motion can be derived from Newton's second law. However, it would be necessary to relate each generalized acceleration ¨ q i to the acceleration of the center of mass of a component of the system and identify all external forces acting on that component. Some of those forces are forces of constraint, for example, the tension on a string or the normal force on a surface. As it was discussed in the previous chapter, the evolution equations can also be obtained from the derivatives of the hamiltonian function. The problem is that in more complicated systems than those considered in the previous chapter, the hamiltonian is not the mechanical energy per unit mass or moment of inertia, but it might be a more complicated expression. This chapter introduces a more general method to derive the equations of motion without the need to identify any forces of constraint.

8.2. Lagrange equations

The total kinetic energy E c of a mechanical system is equal to the sum of all the kinetic energies of translation and rotation of its parts. In general, it is a function that can depend on all the generalized coordinates and velocities and time:

(8)
E c ( q 1 , . . . , q n , ˙ q 1 , . . . , ˙ q n , t )

The motion of a system is usually not completely free, but subject to some constraints; there are forces of constraint resulting from those constraints. For example, in a car moving on a road, the normal force exerted by the road on the tires is the force of constraint that ensures that the trajectory of the car stays over the surface of the road. The static friction on a wheel with traction is also a force of constraint, which ensures that the wheel rotates without sliding on the surface. The fact that the car remains in contact with the road surface as it moves and the friction on the tires prevents them from moving sideways, reduces the three coordinates of position to a single degree of freedom: the displacement along the road. When the wheels rotate without skidding, the angular velocity of the wheels depends on the speed of the car on the road. That dependence also implies a dependence between the angle of rotation of the wheels and the displacement of the car along the road; thus, only one of those two variables is enough to describe the motion of the car and the rotation of the wheels.

Each constraint on a system's motion reduces its number of degrees of freedom. A constraint that can be written as an equation depending on the generalized coordinates of the system is called holonomic. When a system is holonomic, namely all constraints to its motion are holonomic, Newton's second law leads to the following equations (the proof is given in Appendix B ):

(8.2)
d d t E c ˙ q j E c q j = Q j j = 1 , . . . n

where Q j is the j component of the generalized force, defined as:

(8.3)
Q j = i F i · r i q j

where all forces F i (internal or external) are summed and r i is the position of the point force is F i applied. However, some of the forces do not contribute to the generalized force Q j . For example, normal forces and static friction can be neglected because they are applied in a fixed point r i and therefore, F i · d r i = 0 . The tension in a string with constant length can also be ignored because it acts in opposite directions at both ends of the string and the sum of F i · d r i at both ends gives zero.

Among the forces to be included in Q j , some may be conservative and, in these cases, F i · d r i = d U , where U is the potential energy associated to that force. Thus, the contribution of that conservative force to Q j equals U / q j and equations 8.2 can be written as

(8.4)
d d t E c ˙ q j E c q j + U q j = Q j j = 1 , . . . n

where U is the total potential energy of the system and the components Q j of the generalized force include only non-conservative forces. Equations  8.4 are called Lagrange equations, valid for holonomic systems. In the particular case of conservative systems, the right sides of the equations are zero.

Example 8.1

The cart in the figure, with mass m , can move freely on the surface of an incline with mass M . The incline has wheels that let it move freely on the horizontal table. Assuming that the mass of the wheels is much smaller than m and M and that the friction in the axes of the wheels are negligible, find the equations of motion of the system.

Cart on an incline

Resolution. To determine the positions of the cart and the incline at any instant, we need to know position of s of some point of the inclined, on the horizontal table and the position x of a point of the cart on the surface of the incline. The figure above shows how these two variables can be defined. Thus, the system has two degrees of freedom and the generalized velocities are ˙ s and ˙ x .

The derivative ˙ s is also the velocity of the center of mass of the incline. And ˙ x is the velocity of the cart, relative to the incline. Choosing an axis q perpendicular to s and pointing upwards, the vector expressions for the velocity of the incline and the velocity of the cart relative to the incline are:

v p = ˙ s ˆ e s v c/p = ˙ x (cos θ ˆ e s + sin θ ˆ e q )

The velocity of the cart, relative to the table, is the sum of those two vectors:

v c = (˙ s + ˙ x cos θ )ˆ e s + ˙ x sin θ ˆ e q

And the square of its module is,

v 2c = (˙ s + ˙ x cos θ ) 2 + ˙ x 2 sin 2 θ = ˙ s 2 + ˙ x 2 + 2˙ s ˙ x cos θ

The rotational kinetic energy of the wheels is negligible, since their masses are very small. Thus, the total kinetic energy of the system is:

E c = M 2˙ s 2 + m 2 ˙ s 2 + ˙ x 2 + 2˙ s ˙ x cos θ

The gravitational potential energy of the incline can be ignored because it remains constant. Therefore, the potential energy of the system is the gravitational potential energy of the cart:

U = m g x sin θ

Notice that the height of the cart's center of mass of the cart, from the surface of the table, is slightly bigger than x sin θ but the difference is a constant that also be neglected since we only need the derivative of U .

There are no nonconservative forces (or rather, they're effect is negligible). Thus, the right-hand sides of the Lagrange equations  8.4 are zero. In the first Lagrange equation, related to the coordinate x , we must compute the following derivatives:

E c ˙ x = m (˙ x + ˙ s cos θ ) E c x = 0 U x = m g sin θ

And the Lagrange equation becomes,

d d t E c ˙ x E c x + U x = m ¨ x + ¨ s cos θ + g s i n θ = 0

For the coordinate s , the derivatives involved are

E c ˙ s = ( M + m )˙ s + m ˙ x cos θ E c s = 0 U s = 0

And the corresponding Lagrange equation is

d d t E c ˙ s E c s + U s = ( M + m )¨ s + m ¨ x cos θ = 0

Solving the two Lagrange equations for the accelerations ¨ x and ¨ s , we obtain the two equations of motion:

¨ x = ( M + m ) g sin θ M + m sin 2 θ ¨ s = m g sin θ cos θ M + m sin 2 θ

The two accelerations are constant; ¨ x being negative and ¨ s positive. That means that the cart descends the incline while this moves to the right.

Example 8.2

In the system of pulleys and shown in the figure, the pulley in the middle can move vertically while the other two pulleys are fixed to the ceiling. The masses of the two fixed pulleys are m , the mass of the moving pulley is 2 m and the masses of the 3 cylinders are 8 m , 7 m and 5 m ( 7 m for the cylinder in the middle already includes the mass of the support that links it to the moving pulley). The masses of the strings and the friction in the axes of the pulleys are negligible and the string makes the pulleys turn without sliding over them. Find the values of the accelerations of the 3 cylinders.

Cylinders and pulleys

Resolution. This example will also be used to illustrate how to use Maxima to solve Lagrangian mechanics problems. We begin by defining the generalized variables. To determine the vertical position of the cylinders and the moving pulley 3 distances are required, which can be the three variables y 1 , y 2 and y 3 indicated in the figure. y 1 and y 3 are the positions of some points in the two cylinders on the left and the right and y 2 is the position of the moving pulley, which differs from the position of the cylinder in the middle by just a constant.

The constraint that the length of the string remains constant leads to the following equation:

y 1 + 2 y 2 + y 3 = k

where k is a constant. This equation can be used to replace y 3 in terms of y 1 and y 2 . Thus, the system has two degrees of freedom, which can chose as y 1 and y 2 . The two generalized speeds are v 1 = ˙ y 1 and v 2 = ˙ y 2 ; the relation between v 3 and the two generalized speeds is obtained by differentiating the previous equation. The result is trivial, but we will show how it is done in Maxima, would be useful in cases more complicates. Derivatives computed with the Maxima's function diff are partial derivatives. To be find the ordinary derivative of the previous equation with respect to time, it is necessary to declare that the derivative of y 1 is not zero, but the velocity v 1 and similarly for y 2 . As it will be necessary later, we can also declare that the derivatives of v 1 and v 2 , are the accelerations a 1 and a 2 . The command used in Maxima to define the derivative of a variable is gradef. The following commands are used to define set y 3 and v 3 in terms of the degrees of freedom and the corresponding velocities

(%i1) y3: k - y1 - 2*y2$
(%i2) gradef (y1, t, v1)$
(%i3) gradef (y2, t, v2)$
(%i4) gradef (v1, t, a1)$
(%i5) gradef (v2, t, a2)$
(%i6) v3: diff (y3,t);
(%o6)      2 v 2 v 1

Since the string does not slide over the pulleys, the angular velocity of each pulley is ω = V / r , where V is the velocity of the string relative to the center of the pulley and r is the radius of the pulley. Assuming that the disk of each pulley is a uniform cylinder, its moment of inertia relative to its axis is I = M r 2 /2 , where M is the mass of the pulley. Thus, the total rotational kinetic energy is

1 2 I ω 2 = M 4 V 2

The total kinetic energy of the system is:

E c = m 1 v 21 2 + m 2 v 22 2 + m 3 v 23 2 + M 2 v 22 2 + M 1 V 21 4 + M 2 V 22 4 + M 3 V 23 4

Where the indices 1, 2 and 3 refer to the 3 cylinders and the 3 pulleys (from left to right), the masses m i in lowercase letters are the masses of the cylinders and the masses M i in capital letters are the pulley masses. The velocities of the 3 cylinders are v i and the velocities of the strings, relative to the centers of the 3 pulleys are V i . Notice that the pulley in the middle has both translational and rotational kinetic energy.

The gravitational potential energy of the system, up to some constant, is:

U = m 1 g y 1 ( m 2 + M 2 ) g y 2 m 3 g y 3

We now replace the values ​​of the masses in terms of the parameter m and write these energies in terms of the 2 degrees of freedom and the corresponding velocities v 1 = ˙ y 1 and v 2 = ˙ y 2 (notice that V 1 = v 1 , V 2 = v 1 + v 2 and V 3 = v 3 ). That can be done in Maxima as follows:

(%i7) [m1, m2, m3, M1, M2, M3]: [8*m, 7*m, 5*m, m, 2*m, m]$
(%i8) [V1, V2, V3]: [v1, v1+v2, v3]$
(%i9) Ec: expand (m1*v1^2/2 + m2*v2^2/2 + m3*v3^2/2 + M2*v2^2/2
    + M1*V1^2/4 + M2*V2^2/4 + M3*V3^2/4);
(%o9)      16 m v 2 2 + 12 m v 1 v 2 + 15 m v 1 2 2
(%i10) U: expand (-m1*g*y1 - (m2+M2)*g*y2 - m3*g*y3);
(%o10)      g m y 2 3 g m y 1 5 g k m

And the two Lagrange equations are

(%i11) eq1: diff (diff(Ec,v1),t) - diff(Ec,y1) + diff(U,y1) = 0;
(%o11)      3 g m + 12 a 2 m + 15 a 1 m = 0
(%i12) eq2: diff (diff(Ec,v2),t) - diff(Ec,y2) + diff(U,y2) = 0;
(%o12)      g m + 32 a 2 m + 12 a 1 m = 0

Finally, the Lagrange equations are solved to find a 1 and a 2 and the results are used to compute the acceleration a 3

(%i13) solve ([eq1, eq2],[a1,a2]);
(%o13)      a 1 = 9 g 28 , a 2 = 17 g 112
(%i14) subst (%, diff(v3,t));
(%o14)      g 56

It should be noted that the results do not depend on m and the three accelerations are constant. The cylinder on the left side has acceleration equal to 9 g /28 , downwards (because a 1 is positive). The cylinder in the middle and the moving pulley have acceleration 17 g /112 , upwards and the acceleration of the third cylinder is g /56 , upwards. If the 3 cylinders were initially at rest, the cylinder on the left would start to go down and the other two cylinders would go up.

8.3. Equilibrium conditions

In the two examples solved in the previous section, the values ​​obtained for the accelerations were constant. In more general cases, the accelerations will be expressions that depend on the generalized coordinates, velocities and time. The resolution of those systems of differential equations will be undertaken in later chapters.

Even before the equations of motion are solved, it is possible (and convenient) to find the values ​​of the generalized coordinates for which the system will be in equilibrium. The condition for kinetic equilibrium is that the accelerations are zero and if the velocities are also zero, the equilibrium is static.

Recall that in systems with only one degree of freedom, the instability of the equilibrium points is determined from the sign of the derivative of the acceleration, with respect to the generalized coordinate. The equilibrium point is stable when the derivative is negative or unstable when it is positive.

Example 8.3

A motorcyclist moving with velocity v on a curve with radius r tilts his body and the motorcycle at an angle θ , above the horizontal, towards the center of curvature of the curve, to avoid falling to the side. Find the value that θ should have, in terms of v , r and h , which is the distance between the point of contact of the tires with the road, P, and the center of mass, C, of ​​the system.

Biker

Resolution. Due to the inclination of the motorcycle, points P and C are not at the same distance from the center of the curve. Since the distance from point P to the center is r and the speed of point P relative to the road is v , the distance from point C to the center is r h cos θ and the speed of point C is:

v C = r h cos θ r v

in a direction parallel to the velocity of point P. But since the angle θ can vary, point C also has another velocity component h ˙ θ , in the plane perpendicular to the velocity of P. Thus, the kinetic energy of translation is

E c = m 2 h 2 ˙ θ 2 + 1 h r cos θ 2 v 2

There are also rotational kinetic energies, associated with the angular velocity ˙ θ , the angular velocity of the wheels on their axles and the rotation of the whole system in the horizontal plane, since the motorcyclist enters the curve pointing in one direction and leaves pointing in the opposite direction. The calculation of those energies goes beyond the goals of this introductory book; we will consider the case in which the change in those energies can be neglected. The gravitational potential energy of the system is

U = m g h sin θ

The partial derivatives of the energies with respect to θ and ˙ θ are

E c ˙ θ = m h 2 ˙ θ E c ∂θ = m h v 2 r sin θ 1 h r cos θ U ∂θ = m g h cos θ

And the equation of motion is

¨ θ = v 2 h r sin θ 1 h r cos θ g h cos θ

The height of the center of mass, h , is usually much smaller than the radius of the curve. Thus, the expression inside the parentheses is approximately 1 and a good approximation is

¨ θ = v 2 h r sin θ g h cos θ

It the system is in equilibrium, ¨ θ = 0 , and the angle should be:

(8.5)
θ = tan 1 g r v 2

And the derivative of the generalized acceleration with respect to the angle is:

¨ θ ∂θ = v 2 h r cos θ + g h sin θ

Which is positive because 0 θ π /2 . We conclude that the equilibrium is unstable.

Example 8.4

A cart moves on a horizontal table with constant acceleration a . the cart has a vertical pole with a simple pendulum of mass m and length L . Determine the value of the angle θ at which the pendulum remains in equilibrium with respect to the cart. Assume that the mass of the pendulum's string is negligible and that the radius of the sphere is much smaller than L .

Cart with a pendulum

Velocity of the pendulum in the cart

Resolution. The velocity of the cart will always be horizontal and with value a t , where t is the time from when the cart speed was zero. The figure on the right shows the velocity v e/c of the ball, relative to the cart, at an instant when ˙ θ is positive. The relative velocity v e/c equals L ˙ θ and using a coordinates system with axis x in the same direction and sense of the acceleration of the cart and y axis vertical pointing upwards, the components of vector v e/c and the cart's velocity are:

v e/c = L ˙ θ cos θ ˆ ı + sin θ ˆ v c = a t ˆ ı

The velocity of the sphere relative to the table is the sum of those two vectors

v e = a t L ˙ θ cos θ ˆ ı + L ˙ θ sin θ ˆ

In Maxima, if θ is represented by the variable q and ˙ θ by the variable w , the expression of the kinetic energy of the sphere is obtained from its velocity, considering that v 2e = v e · v e

(%i15) ve: [a*t-L*w*cos(q), L*w*sin(q)]$
(%i16) Ec: m*trigsimp(ve.ve)/2;
(%o16)      m w 2 L 2 2 a ˙ θ cos( q ) t w L + a 2 t 2 2

We next introduce the expression of the potential energy of the sphere , compute the derivatives of the generalized coordinate and velocity with respect to time, introduce them into the Laplace equation and solve it to obtain the expression for the angular acceleration ¨ θ which will be denoted by the variable f .

(%i17) U: -m*g*L*cos(q)$
(%i18) gradef (q, t, w)$
(%i19) gradef (w, t, f)$
(%i20) solve (diff(diff(Ec,w),t) - diff(Ec,q) + diff(U,q), f);
(%o20)      f = g sin( q ) a cos( q ) L

Thus, the equation of motion is

(8.6)
¨ θ = a L cos θ g L sin θ

There will be static equilibrium when the angular velocity and acceleration are both zero: ˙ θ = 0 , ¨ θ = 0 , which leads to the condition for the angle in the equilibrium position:

(8.7)
θ = tan 1 a g

And the derivative of the angular acceleration with respect to the angle is

¨ θ ∂θ = a L sin θ g L cos θ

Which is negative because at the point of equilibrium θ is between 0 and π /2 . We conclude that the equilibrium is stable. The pendulum could be swing around the angle θ at the equilibrium position.

Notice that the equation of motion depends on the acceleration of the cart but not on its speed. A measurement of the equilibrium position of the pendulum can be used to find the acceleration of the cart, but not its velocity.

8.4. Dissipative Forces

In all the examples shown in previous sections there were only conservative forces and thus the generalized force was zero. The following examples include non-conservative forces.

Example 8.5

A simple pendulum is formed by a small object of mass m , hanging from a string of length l . The string's mass is negligible compared to m . Determine the equation of motion, including the air resistance.

Simple pendulum

Resolution. The force of resistance of the air is proportional to the square of the velocity of the pendulum, and in the direction opposite to that velocity (see equation  4.14 in chapter 4). Since the velocity of the pendulum has value l ˙ θ , the expression for the air resistance force is:

F r = C l 2 | ˙ θ | ˙ θ

Where C is a constant. Choosing the origin at the point where the string is attached to the support, the position where this resistance force acts is

r = l sin θ ˆ ı cos θ ˆ

And its derivative with respect to θ is

d r d θ = l cos θ ˆ ı + sin θ ˆ = l e θ

where e θ is the tangent unit vector to the circular trajectory of the pendulum, in the sense that θ increases. The generalized force is then

Q θ = F r · d r d θ = C l 2 | ˙ θ | ˙ θ e θ · ( l e θ ) = C l 3 | ˙ θ | ˙ θ

The kinetic and potential energies and their derivatives are similar to the last example in the previous section, setting a = 0

E c = m 2 l 2 ˙ θ 2 U = m g l c o s θ
E c ˙ θ = m l 2 ˙ θ E c ∂θ = 0 U ∂θ = m g l sin θ

The Lagrange equation leads to

(8.8)
¨ θ = g l sin θ C l m | ˙ θ | ˙ θ

8.5. Forces of constraint

One of the advantages of Lagrangian mechanics compared to vector mechanics is not having to identify the forces of constraint, their directions and the points where they are applied. However, in some cases we may want to find those forces. For example, when there is kinetic friction between two surfaces, the frictional force is proportional to the normal force, which is a force of constraint.

There is a method that allows us to calculate the forces of constraint from the Lagrangian equations. It begins by identifying the constraint to which the force of constraint is associated and writing it in the form f ( q 1 , . . . , q n ) = constant . In the case of Example 8.2, the constraint that the length of the string remains constant, y 1 + 2 y 2 + y 3 = k , is responsible for the appearance of the tension force along the string and turns y 3 dependent of y 1 and y 2 . Therefore, to find the tension in the string, the 3 variables ( y 1 , y 2 , y 3 ) are assumed to be independent generalized coordinates, thus increasing the number of Lagrange equations to 3. A function λ called Lagrange multiplier and an additional equation, f ( y 1 , y 2 , y 3 ) = constant , are introduced. In the case of example  8.2 the function f is y 1 + 2 y 2 + y 3 = k and the Lagrange multiplier represents the tension in the string.

The next step is to add a λ f / q j to each Lagrange equation,

(8.9)
d d t E c ˙ q j E c q j + U q j λ f q j = Q j

where j = 1 , . . . n . The following example shows how to calculate the Lagrange multiplier. Each term λ f / q j is the component of the force of constraint acting at q j . In the case of example 8.2, λ f / y 1 , λ f / y 2 and λ f / y 3 are the values ​​of the tension of the string at each of the 3 blocks, which are different.

Example 8.6

A block of mass m slides over an incline of mass M which has wheels that let it move freely on a horizontal table, as shown in the figure. The coefficient of kinetic friction between the block and the incline is µ k . Assuming that the mass of the wheels is much smaller than m and M and that the friction in the wheel axes is negligible, find the equations of motion of the system.

Cart on an incline

Resolution. The figure above shows the two coordinates systems that will be used below. The axes s and q are fixed to the table while the axes x and y move with the incline.

This example is similar to Example 8.1 , but with a non-conservative force: kinetic friction between the block and the incline. Since the kinetic frictional force is equal to µ c N , where N is the normal force between the block and the plane, we need to find that normal force. We then pretend that the block does not keep in contact with the incline and that the two coordinates x and y may change. There are thus 3 degrees of freedom: x , y and s and the equation of the constraint which causes the block to remain in contact with the inclined plane is:

f ( x , y , s ) = y = constant

We introduce a Lagrange multiplier λ and the 3 generalized components of the force of constraint are:

λ f x = 0 λ f y = λλ f s = 0

This shows that the force of constraint points in the direction of axis y and the Lagrange multiplier is the normal force itself between the block and the plane, N .

To determine the components of the velocity in terms of the generalized velocities ( ˙ x , ˙ y , ˙ s ), we will show a method different from the one that was used in the resolution of example  8.1. The position vector of the center of mass of the incline is

r p = s ˆ e s + q ˆ e q

And its derivative is the velocity vector of the incline: v p = ˙ s ˆ e s .

The position of the block relative to the center of mass of the inclined plane is

r b/p = r o + x ˆ ı + y ˆ

where r O is the vector from the center of mass of the incline to the origin of the frame x y . The position of the block relative to the table is r p + r b/p . Since the unit vectors of frame x y with respect to frame s q are

ˆ ı = cos θ ˆ e s + sin θ ˆ e q ˆ = sin θ ˆ e s + cos θ ˆ e q

Then the position of the block with respect to frame s q is,

r b = s + x cos θ y sin θ ˆ e s + q + x sin θ + y cos θ ˆ e q + r o

And the velocity of the block is obtained by differentiation of this last expression

v b = ˙ s + ˙ x cos θ ˙ y sin θ ˆ e s + ˙ x sin θ + ˙ y cos θ ˆ e q

Since the kinetic energy of the wheel rotation is negligible, the total kinetic energy of the system is:

E c = M 2 v 2p + m 2 v 2b = M 2˙ s 2 + m 2 ˙ s 2 + ˙ x 2 + ˙ y 2 + 2˙ s ˙ x cos θ ˙ y sin θ

The height of the block, relative to the table is

h = r b · ˆ e q = q + x sin θ + y cos θ + h o

And, ignoring the constant terms, the gravitational potential energy of the system is

U = m g x sin θ + y cos θ

In this case there is an internal force that performs work: the kinetic frictional force between the block and the incline. To calculate the components Q j of the generalized force, we take into account that in Q j = F · r / q j the vector r is the position of the block relative to the incline r b / p , because the force is internal. Using the expression given above for r b / p , the 3 partial derivatives are r / x = ˆ ı , r / y = ˆ and r / s = 0 . Since the frictional force is µ c N ˆ ı , the three components of the generalized force are then

Q x = µ c R n ˆ ı · ˆ ı = µ c R n Q y = µ c R n ˆ ı · ˆ = 0 Q s = 0

The Lagrange equations 8.9 for the 3 coordinates are

d d t E c ˙ x E c x + U x λ f x = Q x = m ¨ x + ¨ s cos θ + g sin θ = µ c N d d t E c ˙ y E c y + U y λ f y = Q y = m ¨ y ¨ s sin θ + g cos θ N = 0d d t E c ˙ s E c s + U s λ f s = Q s = ( M + m )¨ s + m ¨ x cos θ ¨ y sin θ = 0

These 3 equations can be solved to find the 2 equations of motion for ¨ x and ¨ s in terms of ( x , s , ˙ x , ˙ s ) and the force of constraint N . To substitute y , ˙ y and ¨ y in terms of the generalized coordinates and velocities ( x , s , ˙ x , ˙ s ) we use the constraint equation, f ( x , y , s ) = constant , which in this case is y = constant and, therefore, ¨ y = 0 . Eliminating the terms ¨ y in the Lagrange equations and solving for ¨ x , ¨ s and N we obtain

(8.10)
¨ x = ( M + m ) g β M + m β sin θ ¨ s = m g β cos θ M + m β sin θ N = m M g cos θ M + m β sin θ

where β = sin θ µ c cos θ . In the case when kinetic friction is neglected ( µ c = 0 ), β becomes sin θ and the equations of motion are the same as those obtained in Example 8.1.

Questions

(To check your answer, click on it.)

  1. A very long and homogeneous bar of length L and mass m is falling to the ground. At point A the coefficient of static friction is high enough to prevent point A from moving while the angle θ decreases. Find the expression for the kinetic energy of the bar in terms of the angular velocity ˙ θ Bar leaning against a wall
    1. 1 8 m L 2 ˙ θ 2

    2. 1 6 m L 2 ˙ θ 2

    3. 1 12 m L 2 ˙ θ 2

    4. 1 4 m L 2 ˙ θ 2

    5. 1 2 m L 2 ˙ θ 2
  2. In an Atwood machine, two blocks are hung at the ends of a string passing through a pulley, as shown in the figure. The heavier block descends with constant acceleration and the lighter block rises with the same acceleration value. The masses of the blocks are 3 m and 4 m and the pulley is a homogeneous disk with mass 2 m . Neglecting the friction in the pulley's shaft and the air resistance, determine the value of the acceleration of the blocks. Atwood Machine
    1. g /7
    2. g
    3. 7 g /8
    4. 3 g /4
    5. g /8
  3. The kinetic energy of a particle moving inside a cylinder of radius R is given by m ( R 2 ˙ θ 2 + ˙ z 2 )/2 , where θ and z are the coordinates of the position of the particle in the cylinder, and its potential energy is given by a z 2 /2 + b θ 2 /2 + c z θ , where a , b and c are constant. Find the expression for the acceleration ¨ θ .
    1. b θ + c z m

    2. b θ + c z m R 2

    3. b θ + c z m R

    4. b θ + a z m R

    5. b θ + a z m R 2
  4. The expressions for the kinetic and potential energies of a system with two degrees of freedom, x and θ , are: E c = 5˙ x 2 + 11˙ θ 2 and U = 3 x θ . Find expression for the acceleration ¨ θ .
    1. 3 θ /22
    2. 3 x θ /5
    3. 3 x /22
    4. 3 x θ /22
    5. 3 x /5
  5. The kinetic and gravitational potential energies of a celestial body in orbit around the Sun are given by the expressions:
    E c = m 2( r 2 ˙ θ 2 + ˙ r 2 ) U = 4 π 2 m r
    where m is the mass of the body, r its distance from the Sun, θ is an angle measured on the plane of the orbit with vertex in the Sun and the distances are to be measured in astronomical units and the time in years. Find the equation of motion for ¨ r
    1. r ¨ θ 2 π r 2
    2. r 2 ˙ θ (2 π r ) 2
    3. r ˙ θ 2 2 π r 2
    4. r ˙ θ (2 π r ) 2
    5. r 2 ˙ θ 2 2 π r 2

Problems

  1. In example  8.1, if the masses are m = 0 . 6  kg and M = 2 . 5  kg and the angle is θ = 20 , (a) determine the values of the acceleration of the incline and of the cart relative to the incline. (b) If at an initial moment the incline and the cart are at rest, with x 0 = 20  cm, compute the velocity of the cart, relative to the incline, when it reaches the base of the incline ( x = 0 ) and the time it takes. (c) Find the velocity of the incline when the cart reaches the base of the incline. Cart on an incline
  2. One end of a string is attached to a point P of a pulley, the string is then winded around the pulley and a block is hanged on the other end of the string. The system has a single degree of freedom, which can be the height y of the block. Assume that the pulley is a homogeneous disk with mass equal to the mass of the block and that the mass of the string, the kinetic frictional force on the pulley's axis and the air resistance are negligible. (a) Find the value of the acceleration of the block in as a factor of the acceleration of gravity. (b) If the block starts from rest, determine the value of its velocity after it falls 50 cm. Pulley and Block
  3. A particle with mass m = 2 kg moves along a vertical path with equation y = x 2 , where x is measured horizontally and y vertically (both in meters). Thus, the motion of the particle has only one degree of freedom, which can be chosen as the x coordinate.
    (a) Find the expression of the particle's kinetic energy in terms of x .
    (b) Find the expression of the particle's gravitational potential energy in terms of x (use the value g = 9 . 8  m/s2).
    (c) Assuming that there are no dissipative forces acting on the particle, use the Lagrange equation to find its equation of motion.
    (d) Find the equilibrium points of the system in phase space, and determine whether they are stable or unstable.
  4. Cylinder A in the figure has a mass of 36 grams, cylinder B has a mass of 24 grams and the moment of inertia of the double pulley is 4.43×10−7 kg·m2. The pulley is formed by two discs, of rays 5 cm and 8 cm, glued to each other. Each cylinder is connected to the end of a string with the opposite end attached to the pulley, making the pulley rotate. (a) Neglecting the friction forces in the pulley's shaft and the air resistance, determine the acceleration ​of each cylinder and say whether they are upwards or downwards. (b) Find the tension in each of the two strings. Double pulley
  5. In the system shown in the figure, the mass of the wheels and pulley and the friction in their axles can be neglected. (a) Determine the expressions for the kinetic and potential energies of the system, in terms of the angle θ and the horizontal position x of the cart. (b) Find the expressions of the cart's acceleration and the angular acceleration ¨ θ . (c) Find the value of the angle θ in the equilibrium position of the pendulum and say whether the equilibrium is stable or unstable. (d) Compute the value of the cart's acceleration when the pendulum remains in its equilibrium position.

    Cart with a pendulum
  6. In the system shown in the figure, the fixed pulley has mass m and the moving pulley has mass 2  m (both can be considered uniform disks). The cars mass is 20 m and the mass of the cylinder plus the support connecting it to the moving pulley is 8 m . Assume that the mass of the string and the wheels, the kinetic frictional forces on the pulley and the wheels and the air resistance are all negligible.
    (a) Show that the kinetic and potential energies of the system in terms of the cylinder's height y are given by the expressions:
    E c = 93 2 m ˙ y 2 U = 10 m g y

    (b) Compute the values of the accelerations of the cylinder and the cart. Cart linked to a falling cylinder
  7. A block of mass m slides down an incline which makes an angle θ with the horizontal. The coefficient of kinetic friction between the block and the incline is µ k . Using the Lagrange equation with a multiplier, find the expressions for the normal force of the plane on the block and the acceleration of the block, ¨ x (neglect the air resistance). Block on incline
  8. The bar in the figure is homogeneous, with mass m and length L = 2  m, and is supported on the floor at point A and on a wall at point B. At an initial instant, the bar is placed at rest, with initial angle θ = 30 . If the floor and wall are very smooth, the frictional forces at points A and B are negligible and the bar descends until the angle θ decreases to 0. Assume that points A and B remain always in contact with the floor and wall, that the air resistance is negligible and that the thickness of the bar is much smaller than its length. Bar leaning against a wall
    (a) Prove that at any time the velocity of the bar's center of mass is equal to L ˙ θ /2
    (b) Find the expression of the kinetic energy in terms of the angle θ .
    (c) Find the expression of the gravitational potential energy in terms of the angle θ .
    (d) Find the expression of the angular acceleration.
    (e) Find the expression of the angular velocity.
    (f) The time that the bar takes to fall to the ground is given by the integral 0 π /6 d θ ˙ θ . Using the expression for ˙ θ obtained in the previous item, compute that time. (The integral is improper and can not be calculated analytically, but can be calculated numerically, using the function quad_qags of Maxima.)
  9. In a simple pendulum, with a small object of mass m hung by a string of length l and negligible mass, the point where the string is fixed moves up and down according to the expression a cos( b t ) , where a e b are two constants.
    (a) Ignoring the air resistance, determine the expressions for the kinetic and potential energies in terms of the angle θ between the pendulum and the vertical.
    (b) Find the equation of motion for ¨ θ .
    (c) Explain at what values ​​of the constants a and b the equilibrium point θ = 0 will be stable or unstable.
  10. The jumper in the figure bends grabing his legs at point P, in order to rotate faster, and then extends his body at Q to reduce the rotation as he enters the water. The changes in his angular velocity are a consequence of the change in his moment of inertia.
    (a) If the moment of inertia of the jumper relative to the center of mass is I , which depends on time, write the expressions for their kinetic and potential energies in terms of the position ( x , y ) of the center of mass and the angle of rotation θ .
    (b) Using the Lagrange equation for θ , show that the angular momentum, L = I ˙ θ , remains constant.
    (c) If at the highest point P of the trajectory the moment of inertia is 3.28 kg·m2 and the angular velocity ˙ θ = 4s 1 , and at point Q the moment of inertia is 28.2 kg·m2, determine the angular velocity of the jumper at point Q. Diver
  11. The gravitational potential energy of a celestial body with mass m , in orbit around another body of mass M , is given by the expression (see problem  2 of chapter 6 ):
    U g = G M m r
    where G is the universal gravitational constant and r the distance between the two bodies. It can be shown that the possible orbits of the celestial body are always plane. As such, the orbital motion has two degrees of freedom that can be r and an angle θ measured in the plane of the orbit, with vertex in the body of mass M . In this system of polar coordinates, the square of the velocity of the body with mass m is ( r 2 ˙ θ 2 + ˙ r 2 ).
    (a) Using the Lagrange equation for θ , prove that the angular momentum
    L = m r 2 ˙ θ
    of the body with mass m , relative to the body of mass M , remains constant.
    (b) Find the equation of motion for ¨ r and show that it depends only on r and neither on ˙ r nor on θ or ˙ θ .

Answers

Questions: 1. B. 2. E. 3. B. 4. C. 5. C.

Problems

  1. (a) ¨ x = 4 . 043  m/s2 and ¨ s = 0 . 735  m/s2  (b) ˙ x = 1 . 272  m/s, t = 0 . 315  s.  (c) ˙ s = 0 . 231  m/s.
  2. (a) 2 g /3   (b) 2.56 m/s.
  3. (a) E c = ˙ x 2 4 x 2 + 1   (b) U g = 19 . 6 x 2   (c) ¨ x = x 4˙ x 2 + 19 . 6 4 x 2 + 1   (d) There is a single stable equilibrium point at the origin.
  4. (a) a A = 0 . 2409  m·s−2, upwards, a B = 0 . 3855  m·s−2, downwards.  (b) The tension in the string attached to cylinder A is0.362 N and the tension in the string attached to cylinder B is0.226 N.
  5. (a) In SI units,  E c = 3 . 03˙ x 2 + 0 . 0012˙ θ 2 0 . 012˙ x ˙ θ cos θ U = 9 . 8 x 0 . 1176cos θ

    (b) ¨ x = 2450 147cos θ sin θ 3˙ θ 2 sin θ 15(101 cos 2 θ ) ¨ θ = 2450cos θ 14847sin θ 3˙ θ 2 cos θ sin θ 3(101 cos 2 θ )

    (c) 9.37°, stable.  (d) 1.617 m/s2.
  6. (b) Cylinder: 10 g /93 1 . 05  m/s2. Cart: 20 g /93 2 . 11  m/s2.
  7. N = m g cos θ ,    ¨ x = g sin θ µ c cos θ   
  8. (a) The position of the center of mass is ( x ˆ ı + y ˆ )/2 and the velocity of the center of mass is the derivative of that expression. Replacing x = L cos θ and y = L sin θ the result given is obtained.

    (b) 1 6 m L 2 ˙ θ 2   (c) L 2 m g sin θ   (d ) ¨ θ = 3 g 2 L cos θ   (e) ˙ θ = 3 g L 1 2 sin θ   (f ) 0.3977 s.
  9. (a) Kinetic energy:  m 2 l 2 ˙ θ 2 + a 2 b 2 sin 2 ( b t ) 2 a b l ˙ θ sin θ sin( b t )
    Potential energy:  m g ( a cos( b t ) l cos θ )

    (b) ¨ θ = sin θ l a b 2 cos( ω t ) g

    (c) If a b 2 g , the equilibrium is stable, otherwise the equilibrium is unstable.
  10. (a) E c = m 2 ˙ x 2 + ˙ y 2 + 1 2 I ˙ θ 2 , U = m g y   (b) d( I ˙ θ ) d t = 0 , which implies I ˙ θ = constant ..  (c) 0.465 s−1
  11. (a) The Lagrange equation is:  d d t m r 2 ˙ θ = 0 , which implies m r 2 ˙ θ constant.
    (b) ¨ r = L 2 m 2 r 3 G M r 2 , where L , m , G and M are constants.
Wrong

(Click to continue)

Wrong

(Click to continue)

Wrong

(Click to continue)

Correct

(Click to continue)

Wrong

(Click to continue)

Wrong

(Click to continue)

Wrong

(Click to continue)

Wrong

(Click to continue)

Correct

(Click to continue)

Wrong

(Click to continue)

Wrong

(Click to continue)

Correct

(Click to continue)

Wrong

(Click to continue)

Wrong

(Click to continue)