L2 Degrees off Freedom and Constraints, Rectilinear Motion

L2 Degrees off Freedom and Constraints, Rectilinear Motion...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: S. Widnall 16.07 Dynamics Fall 2009 Version 1.0 Lecture L2 - Degrees of Freedom and Constraints, Rectilinear Motion Degrees of Freedom Degrees of freedom refers to the number of independent spatial coordinates that must be specified to determine the position of a body. If the body is a point mass, only three coordinates are required to determine its position. On the other hand, if the body is extended, such as an aircraft, three position coordinates and three angular coordinates are required to completely specify its position and orientation in space. Kinematic Constraints In many situations the number of independent coordinates will be reduced below this number, either because the number of spacial dimensions is reduced or because there are relationships specified among the spatial coordinates. When setting up problems for solution it is useful to think of these relationships as constraints. For example, if a point mass is constrained to move in a plane (two dimensions) the number of spatial coordinates necessary to describe its motion is two. If instead of being a point mass, this body has extended dimensions, such as a flat plate confined to a plane, it requires three coordinates to specify its position and orientation: two position coordinates and one angular coordinate. 1 If a particle is confined to move on a curve in either two or three dimensions, such as a bead moving on a wire, the number of independent coordinates necessary to describe its motion is one. Another source of constraints on the motion of particles is connections between them. For example, the two particle connected by a cable passing over a pulley are constrained to move in equal and opposite directions. More complex arrangements are possible and can be analyzed using these ideas. Two gears in contact are constrain to move together according to their individual geometry. A cylinder rolling on a plane is constrained in two ways. Contact with the plane reduces the two-dimensional motion to one spatial coordinate along the plane, and the constraint of rolling provides a relationship between the angular coordinates and the spatial position, resulting in a single degree of freedom system. 2 Internal Force-Balance Constraints Another type of constraint occurs when we consider the of a system of particles and the necessary force balance that occurs between the parts. These constraints follow directly from Newton's third law: the force of action and reaction between two bodies are equal in magnitude and opposite in direction. We will pursue these ideas in greater depth later in the course. For now, we will give a simple example to illustrate the principle. Consider the systems shown in a) and b). System a) consists of two masses m in contact resting on a frictionless plane in the presence of gravity. A force F is applied to mass 1, and it is obvious that the two masses will accelerate at a = F/(2m). If we look at the two masses separately, we can determine what internal force must exist between them to cause the motion. It is clear that each mass feels a net force of F/2, since its acceleration is a = F/(2m). This net force arises because between the two masses there is an equal and opposite force F/2 acting across the interface. Another way to look at this is that the interface between the bodies is a "body" of zero mass, and therefore can have no net force acting upon it otherwise its acceleration would be infinite. System b) is a bit more complex, primarily because the forces between one mass and the mass above it are shear forces and must be supplied by friction. Assuming that the friction coefficient is large enough to accelerate the three masses an equal amount given by a = F/(3m), by the reasoning we have discussed, the force balance is as sketched in b): equal and opposite normal forces F 2/3 on the vertical surfaces, and equal and opposite shear forces F 1/3 on the horizontal surfaces. Rectilinear Motion In many case we can get an exact expression for the position of a particle as a function of time. We start by considering the simple motion of a particle along a straight line. The position of particle A at any instant can be specified by the coordinate s with origin at some fixed point O. 3 The instantaneous velocity is v= ds =v . dt (1) We will be using the "dot" notation, to indicate time derivative, e.g. () d/dt. Here, a positive v means is moving in the opposite direction. The acceleration is a= dv d2 s =v= 2 =s. dt dt that the particle is moving in the direction of increasing s, whereas a negative v, indicates that the particle (2) The above expression allows us to calculate the speed and the acceleration if s and/or v are given as a function of t, i.e. s(t) and v(t). In most cases however, we will know the acceleration and then, the velocity and the position will have to be determined from the above expressions by integration. Determining the velocity from the acceleration From a(t) If the acceleration is given as a function of t, a(t), then the velocity can be determined by simple integration of equation (2), v(t) = v0 + t a(t) dt . (3) t0 Here, v0 is the velocity at time t0 , which is determined by the initial conditions. From a(v) If the acceleration is given as a function of velocity a(v), then, we can still use equation (2), but in this case we will solve for the time as a function of velocity, t(v) = t0 + v v0 dv . a(v) (4) Once the relationship t(v) has been obtained, we can, in principle, solve for the velocity to obtain v(t). A typical example in which the acceleration is known as a function of velocity is when aerodynamic drag forces are present. Drag forces cause an acceleration which opposes the motion and is typically of the form velocity). a(v) v 2 (the sign "" means proportional to, that is, a(v) = v 2 for some , which is not a function of From a(s) When the acceleration is given as a function of s then, we need to use a combination of equations (1) and (2), to solve the problem. From a= dv dv ds dv = =v dt ds dt ds 4 (5) we can write a ds = v dv . This equation can now be used to determine v as a function of s, s 2 v 2 (s) = v0 + 2 a(s) ds . s0 (6) (7) where, v0 , is the velocity of the particle at point s0 . Here, we have used the fact that, v v v2 v2 v2 - 0 . v dv = d( ) = 2 2 2 v0 v0 A classical example of an acceleration dependent on the spatial coordinate s, is that induced by a deformed linear spring. In this case, the acceleration is of the form a(s) s. Of course, when the acceleration is constant, any of the above expressions (3, 4, 7), can be employed. In this case we obtain, v = v0 + a(t - t0 ), If a = g, this reduces to the familiar v = v0 + g(t - t0 ), or 2 v 2 = v0 + 2g(s - s0 ) . or 2 v 2 = v0 + 2a(s - s0 ) . Determining the position from the velocity Once we know the velocity, the position can be found by integrating ds = vdt from equation (1). Thus, when the velocity is known as a function of time we have, t s = s0 + v(t) dt . t0 (8) If the velocity is known as a function of position, then t = t0 + Here, s0 is the position at time t0 . s s0 ds . v(s) (9) It is worth pointing out that equation (6), can also be used to derive an expression for v(s), given a(v), s v v s - s0 = ds = dv . (10) a(v) s0 v0 This equation can be used whenever equation (4) is applicable and gives v(s) instead of t(v). For the case of constant acceleration, either of equations (8, 9), can be used to obtain, 1 s = s0 + v0 (t - t0 ) + a(t - t0 )2 . 2 5 In many practical situations, it may not be possible to carry out the above integrations analytically in which case, numerical integration is required. Usually, numerical integration will also be required when either the velocity or the acceleration depend on more than one variable, i.e. v(s, t), or, a(s, v). Example Terminal Velocity Terminal velocity occurs when the acceleration becomes zero and the velocity Consider an air-dropped payload starting from rest. The force on the body is a combination of gravity and air drag and has the form 1 F = mg - v 2 CD A 2 Applying Newton's law and solving for the acceleration a we obtain 1 CD A a = g - v 2 m 2 The quantity m CD A Reentry, Ballistic Coefficient, Terminal Velocity (11) (12) characterizes the combined effect of body shape and mass on the acceleration; it is an m CD A . important parameter in the study of reentry; it is called the Ballistic Coefficient. It is defined as = Unlike many coefficients that appear in aerospace problems, the Ballistic Coefficient is not non-dimensional, but has units of mass/length2 or kg/meters2 in mks units. Also, in some applications, the ballistic coefficient is defined as the inverse, B = CD A m , so it pays to be careful in its application. Equation 13 then becomes 1 a = g - v 2 / 2 (13) Terminal velocity occurs when the force of gravity equals the drag on the object resulting in zero acceleration. This balance gives the terminal velocity as vterminal = 2gm = CD A 2g (14) For the Earth, atmospheric density at sea level is = 1.225kg/m3 ; we shall deal with the variation of atmospheric density with altitude when we consider atmospheric reentry of space vehicles. Typical value of range from = 1 (Assuming CD = .5, a tennis ball has = 35.) to = 1000 for a reentry vehicle. As an example, consider a typical case of = 225, where the various parameters then give the following expression for the acceleration a = g - 0.002725v 2 m/s . Here g = 9.81m/s2 , is the acceleration due to gravity and v is the downward velocity. It is clear from this expression that initially the acceleration will be g. Therefore, the velocity will start to increase and keep on increasing until a = 0, at which point the velocity will stay constant. The terminal velocity is then given by, 2 0 = g - 0.002725 vf or vf = 60m/s . 6 To determine the velocity as a function of time, the acceleration can be re-written introducing the terminal velocity as, a = g(1 - (v/vf )2 ) . We then use expression (4), and write 1 v dv 1 v 1/2 1/2 vf vf + v t= = ( + )d v = ln . g 0 1 - (v/vf )2 g 0 1 + (v/vf ) 1 - (v/vf ) 2g vf - v Solving for v we obtain, v = vf e2gt/vf - 1 m/s . e2gt/vf + 1 We can easily verify that for large t, v = vf . We can also find out how long does it take for the payload to reach, say, 95% of the terminal velocity, t= vf 1.95 ln = 11.21s . 2g 0.05 To obtain an expression for the velocity as a function of the traveled distance we can use expression (10) and write 1 s= g v 2 vf v dv = - ln(1 - (v/vf )2 ) . 1 - (v/vf )2 2g 0 Solving for v we obtain v = vf We see that for, say, v = 0.95vf , s = 427.57m. This is the distance traveled by the payload in 11.21s, which can be compared with the distance that would be traveled in the same time if we were to neglect air resistance, sno drag 2 1 - e-2gs/vf m/s . = gt2 /2 = 615.75m. Example Spring-mass system Here, we consider a mass allowed to move without friction on a horizontal slider and subject to the force exerted by a linear spring. Initially the system is in equilibrium (no force on the spring) at s = 0. Suddenly, the mass is given a velocity v0 and then the system is left free to oscillate. We know that the effect of the spring is to cause an acceleration to the body, opposing the motion, of the form a = -s, where > 0 is a constant. Using equation (7), we have 2 v 2 = v0 - s2 . The displacement can now be obtained using expression (9), s s ds 1 t= = arcsin , 2 - s2 v0 v0 0 7 which gives, v0 s = sin t . Finally, the velocity as a function of time is simply, v = v0 cos undamped harmonic oscillator. t. We recognize this motion as that of an References [1] R. Dugas, A History of Mechanics, Dover, 1988. ADDITIONAL READING J.L. Meriam and L.G. Kraige, Engineering Mechanics, DYNAMICS, 5th Edition 1/1, 1/2, 1/3, 1/4, 1/5 (Effect or Altitude only), 1/6, 1/7, 3/2 8 MIT OpenCourseWare http://ocw.mit.edu 16.07 Dynamics Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. ...
View Full Document

Ask a homework question - tutors are online