16.07 Dynamics Fall 2004 Lecture D2- Curvilinear motion Cartesian Coordinates We will start by studying the motion of a particle. We think of a particle as a body which has mass. but has negligible dimensions. Treating bodies as particles is, of course, an idealization which involves an approximation. This approximation may be perfectly acceptable in some situations and not adequate in some other cases. For instance, if we want to study the motion of planets, it is common to consider each planet as a particle. This same idealization is totally meaningless if we want to study, for example, the motion of a rover on the surface of the planet Kinematics of curvilinear motion In dymamics we study the motion and the forces that cause, or are generated as a result of, the motion. Before we can explore these connections we will look first at the description of motion irrespective of the forces that produce them. This is the domain of kinematics. On the other hand, the connection between forces and motions is the domain of kinetics and will be the subject of the next lecture Position vector and path We consider the general situation of a particle moving in a three dimensional space. To locate the position of a particle in space we need to set up an origin point, O, whose location is known. The position of a particle A, at time t, can then be described in terms of the position vector, r, joining points O and A. In general this particle will not be still, but its position will change in time. Thus, the position vector will be a function of time, i.e. r(t). The curve in space described by the particle is called the path, or trajectory s(t) e introduce the path or arc length coordinate, s, which measures the distance travelled by the particle along the curved path. Note that for the particular case of rectilinear motion(considered in the review notes)the arc length coordinate and the coordinate. s. are the same
J. Peraire 16.07 Dynamics Fall 2004 Version 1.1 Lecture D2 - Curvilinear Motion. Cartesian Coordinates We will start by studying the motion of a particle. We think of a particle as a body which has mass, but has negligible dimensions. Treating bodies as particles is, of course, an idealization which involves an approximation. This approximation may be perfectly acceptable in some situations and not adequate in some other cases. For instance, if we want to study the motion of planets, it is common to consider each planet as a particle. This same idealization is totally meaningless if we want to study, for example, the motion of a rover on the surface of the planet. Kinematics of curvilinear motion In dynamics we study the motion and the forces that cause, or are generated as a result of, the motion. Before we can explore these connections we will look first at the description of motion irrespective of the forces that produce them. This is the domain of kinematics. On the other hand, the connection between forces and motions is the domain of kinetics and will be the subject of the next lecture. Position vector and Path We consider the general situation of a particle moving in a three dimensional space. To locate the position of a particle in space we need to set up an origin point, O, whose location is known. The position of a particle A, at time t, can then be described in terms of the position vector, r, joining points O and A. In general, this particle will not be still, but its position will change in time. Thus, the position vector will be a function of time, i.e. r(t). The curve in space described by the particle is called the path, or trajectory. We introduce the path or arc length coordinate, s, which measures the distance travelled by the particle along the curved path. Note that for the particular case of rectilinear motion (considered in the review notes) the arc length coordinate and the coordinate, s, are the same. 1
Using the path coordinate we can obtain an alternative representation of the the motion of the particle Consider that we know r as a function of s, i. e. r(s), and that, in addition we know the value of the path coordinate as a function of time t, i. e. s(t). We can then calculate the speed at which the particle moves on the path simply as u=s= ds/dt. We also compute the rate of change of speed as at =s=d-s/dt We consider below some motion examples in which the position vector is referred to a fixed cartesian ordinate system. Example Motion along a straight line in 2D Consider for illustration purposes two particles that move along a line defined by a point P and a unit vector m. We further assume that at t=0, both particles are at point P. The position vector of the first particle is given by ri(t)=rp+mt=(rPr+m2t)i+(rPy+m,t)3, whereas the position vector of the second particle is given by r2(t)=rp+mt2=(rPr+mt2)i+(rPy+m,t2)j r() Clearly the path for these two particles is the same, but the speed at which each particle moves along the path is different. This is seen clearly if we parametrize the path with the path coordinate, s. That is, we write r(s)=rp+ms=(rPr + mxs)i+(rPy+mys)3. It is straightforward to verify that s is indee the path coordinate i.e. the distance between two points r(s) and r(s+As) is equal to As. The twe motions introduced earlier simply correspond to two particles moving according to s1(t)=t and s2(t)=t2, respectively. Thus, ri(t=r(si(t)) and r2(t)=r(s2(t)) out that, in many situations, we will not have an expression for the path as a function of s. It is in fact possible to obtain the speed directly from r(t) without the need for an arc length parametrization of the trajectory Velocity Vector We consider the positions of the particle at two different times t and t+ At, where At is a small increment of time. Let Ar=r(r+ At)-r(t), be the displacement vector as shown in the diagram
Using the path coordinate we can obtain an alternative representation of the the motion of the particle. Consider that we know r as a function of s, i.e. r(s), and that, in addition we know the value of the path coordinate as a function of time t, i.e. s(t). We can then calculate the speed at which the particle moves on the path simply as v = ˙s ≡ ds/dt. We also compute the rate of change of speed as at = ¨s = d 2 s/dt2 . We consider below some motion examples in which the position vector is referred to a fixed cartesian coordinate system. Example Motion along a straight line in 2D Consider for illustration purposes two particles that move along a line defined by a point P and a unit vector m. We further assume that at t = 0, both particles are at point P. The position vector of the first particle is given by r1(t) = rP + mt = (rPx + mxt)i + (rPy + myt)j, whereas the position vector of the second particle is given by r2(t) = rP + mt 2 = (rPx + mxt 2 )i + (rPy + myt 2 )j. Clearly the path for these two particles is the same, but the speed at which each particle moves along the path is different. This is seen clearly if we parametrize the path with the path coordinate, s. That is, we write r(s) = rP + ms = (rPx + mxs)i + (rPy + mys)j. It is straightforward to verify that s is indeed the path coordinate i.e. the distance between two points r(s) and r(s + ∆s) is equal to ∆s. The two motions introduced earlier simply correspond to two particles moving according to s1(t) = t and s2(t) = t 2 , respectively. Thus, r1(t) = r(s1(t)) and r2(t) = r(s2(t)). It turns out that, in many situations, we will not have an expression for the path as a function of s. It is in fact possible to obtain the speed directly from r(t) without the need for an arc length parametrization of the trajectory. Velocity Vector We consider the positions of the particle at two different times t and t + ∆t, where ∆t is a small increment of time. Let ∆r = r(r + ∆t) − r(t), be the displacement vector as shown in the diagram. 2
8(t) (t+△t) r(t) t) The average velocity of the particle over this small increment of time is Thich is a vector whose direction is that of Ar and whose magnitude is the length of Ar divided by At. If At is small, then Ar will become tangent to the path, and the modulus of Ar will be equal to the distance the particle has moved on the curve As The instantaneous velocity vector is given by △r_dr(t) and is always tangent to the path. The magnitude, or speed, is given by △sds Acceleration vector In an analogous manner, we can define the acceleration vector. Particle A at time t, occupies position r(t), and has a velocity u(t), and, at time t+ At, it has position r(t +At)=r(t)+Ar, and velocity u(t+At)=v(t)+Av. Considering an infinitesimal time increment, we define the acceleration vector as the derivative of the velocity vector with respect to time, a= lim △Udvd2r △t-0△ t dt dt2 r+△r We note that the acceleration vector will reflect the changes of velocity in both magnitude and direction. The acceleration vector will, in general, not be tangent to the trajectory (in fact it is only tangent when the velocity vector does not change direction). A sometimes useful way to visualize the acceleration vector is to
The average velocity of the particle over this small increment of time is vave = ∆r ∆t , which is a vector whose direction is that of ∆r and whose magnitude is the length of ∆r divided by ∆t. If ∆t is small, then ∆r will become tangent to the path, and the modulus of ∆r will be equal to the distance the particle has moved on the curve ∆s. The instantaneous velocity vector is given by v = lim ∆t→0 ∆r ∆t ≡ dr(t) dt ≡ r˙ , (1) and is always tangent to the path. The magnitude, or speed, is given by v = |v| = lim ∆t→0 ∆s ∆t ≡ ds dt ≡ s . ˙ Acceleration Vector In an analogous manner, we can define the acceleration vector. Particle A at time t, occupies position r(t), and has a velocity v(t), and, at time t + ∆t, it has position r(t + ∆t) = r(t) + ∆r, and velocity v(t + ∆t) = v(t) + ∆v. Considering an infinitesimal time increment, we define the acceleration vector as the derivative of the velocity vector with respect to time, a = lim ∆t→0 ∆v ∆t ≡ dv dt = d 2r dt2 . (2) We note that the acceleration vector will reflect the changes of velocity in both magnitude and direction. The acceleration vector will, in general, not be tangent to the trajectory (in fact it is only tangent when the velocity vector does not change direction). A sometimes useful way to visualize the acceleration vector is to 3
translate the velocity vectors, at different times, such that they all have a common origin, say, O. Then the heads of the velocity vector will change in time and describe a curve in space called the hodograph. We then see that the acceleration vector is, in fact, tangent to the hodograph at every point V+△v Expressions(1)and (2)introduce the concept of derivative of a vector. Because a vector has both magnitude and direction, the derivative will be non-zero when either of them changes(see the review notes on s). In general, the derivative of a vector will have a component which is parallel to the vector itself. due to the magnitude change; and a component which is orthogonal to it, and is due to the direction Note Unit tanget and arc-length parametrization The unit tangent vector to the curve can be simply calculated as u/v It is clear that the tangent vector depends solely on the geometry of the trajectory and not on the speed at which the particle moves along the trajectory. That is, the geometry of the trajectory determines the tangent vector, and hence the direction of the velocity vector. How fast the particle moves along the trajectory determines the magnitude of the velocity vector. This is clearly seen if we consider the arc-length parametrization of the trajectory r(s). Then, applying the chain rule for differentiation, we have that dr dr ds dt here,s=v, and we observe that dr/ds= et. The fact that the modulus of dr/ds is always unity indicates that the distance travelled, along the path, by r(s), (recall that this distance is measured by the coordinate s), per unit of s is, in fact, unity!. This is not surprising since by definition the distance between two neighboring points is ds, i.e. dr= ds Cartesian Coordinates When working with fixed cartesian coordinates, vector differentiation takes a particularly simple form. Sin the vectors i, j, and k do not change, the derivative of a vector A(t)=A(t)i+Ay(t)3+A2(t)k, is simply A(t)=A,(t)i+ A,(t)j+A2(t)k. That is, the components of the derivative vector are simply the derivatives of the compa
translate the velocity vectors, at different times, such that they all have a common origin, say, O′ . Then, the heads of the velocity vector will change in time and describe a curve in space called the hodograph. We then see that the acceleration vector is, in fact, tangent to the hodograph at every point. Expressions (1) and (2) introduce the concept of derivative of a vector. Because a vector has both magnitude and direction, the derivative will be non-zero when either of them changes (see the review notes on vectors). In general, the derivative of a vector will have a component which is parallel to the vector itself, and is due to the magnitude change; and a component which is orthogonal to it, and is due to the direction change. Note Unit tanget and arc-length parametrization The unit tangent vector to the curve can be simply calculated as et = v/v. It is clear that the tangent vector depends solely on the geometry of the trajectory and not on the speed at which the particle moves along the trajectory. That is, the geometry of the trajectory determines the tangent vector, and hence the direction of the velocity vector. How fast the particle moves along the trajectory determines the magnitude of the velocity vector. This is clearly seen if we consider the arc-length parametrization of the trajectory r(s). Then, applying the chain rule for differentiation, we have that, v = dr dt = dr ds ds dt = etv , where, ˙s = v, and we observe that dr/ds = et. The fact that the modulus of dr/ds is always unity indicates that the distance travelled, along the path, by r(s), (recall that this distance is measured by the coordinate s), per unit of s is, in fact, unity!. This is not surprising since by definition the distance between two neighboring points is ds, i.e. |dr| = ds. Cartesian Coordinates When working with fixed cartesian coordinates, vector differentiation takes a particularly simple form. Since the vectors i, j, and k do not change, the derivative of a vector A(t) = Ax(t)i + Ay(t)j + Az(t)k, is simply A˙ (t) = A˙ x(t)i+A˙ y(t)j +A˙ z(t)k. That is, the components of the derivative vector are simply the derivatives of the components. 4
Thus, if we refer the position, velocity, and acceleration vectors to a fixed cartesian coordinate syst have r(t)=r(t)i+y(t)j+2(t)k v(t)=tr(t1)i+vy(t)j+v(t)k=x(t)i+y(t)j+i(t)k=r(t) a(t)=ar(t)i+ay(t)j+a2(t)k=ix(t)i+iy(t)j+i2(t)k=v(t) (5) Here, the speed is given by=V吗+吗+, and the magnitude of the acceleration is a=V"+吗+吗 The advantages of cartesian coordinate systems is that they are simple to use, and that if a is constant, or a function of time only, we can integrate each component of the acceleration and velocity independently as hown in the ballistic motion example Example Circular Motion We consider motion of a particle along a circle of radius R at a constant speed vo. The parametrization of circle in terms of the arc length is r(s)=RcosRi+Rsin(R)j Since we have a constant speed vo, we have s= vot. Thus, r(t)=Rcos(p-i+ Rsin(p )3 The velocity 2=-(m+mc方
Thus, if we refer the position, velocity, and acceleration vectors to a fixed cartesian coordinate system, we have, r(t) = x(t)i + y(t)j + z(t)k (3) v(t) = vx(t)i + vy(t)j + vz(t)k = ˙x(t)i + ˙y(t)j + ˙z(t)k = r˙(t) (4) a(t) = ax(t)i + ay(t)j + az(t)k = ˙vx(t)i + ˙vy(t)j + ˙vz(t)k = v˙(t) (5) Here, the speed is given by v = q v 2 x + v 2 y + v 2 z , and the magnitude of the acceleration is a = q a 2 x + a 2 y + a 2 z . The advantages of cartesian coordinate systems is that they are simple to use, and that if a is constant, or a function of time only, we can integrate each component of the acceleration and velocity independently as shown in the ballistic motion example. Example Circular Motion We consider motion of a particle along a circle of radius R at a constant speed v0. The parametrization of a circle in terms of the arc length is r(s) = R cos( s R )i + R sin( s R )j . Since we have a constant speed v0, we have s = v0t. Thus, r(t) = R cos(v0t R )i + R sin(v0t R )j . The velocity is v(t) = dr(t) dt = −v0 sin(v0t R )i + v0 cos(v0t R )j , 5
which, clearly, has a constant magnitude u= vo. The acceleration is (t) dr(t) CO in(-n)3 Note that, the acceleration is perpendicular to the path(in this case it is parallel to r), since the velocity vector changes direction, but not magnit We can also verify that, from r(s), the unit tangent vector, et, could be computed directly as dr(s) L cOS Example Motion along a helix The equation r(t)=Rcos ti+Rsintj+htk, defines the motion of a particle moving on a helix of radius R, and pitch 2Th, at a constant speed. The velocity vector is given by dr rsin ti+ rcos ti+ hk and the acceleration vector is given by, du Rcos ti+-Rsinti In order to determine the speed at which the particle moves we simply compute the modulus of the velocit v=Ju=VR Sin t+R2 cost+h2=VR2 +h2 If we want to obtain the equation of the path in terms of the arc-length coordinate we simply write drl=udt=VR2 +h2dt Integrating, we obtain s =So +VR2+h2t, where so corresponds to the path coordinate of the particle t time zero. Substituting t in terms of s, we obtain the expression for the position vector in terms of the arc-length coordinate In this case, r(s)=Rcos(s/VR2+h2)i+R sin(s/vR2+h2)j+hs/vR2+h2k.The figure below shows the particle trajectory for R= l and h=0.1
which, clearly, has a constant magnitude |v| = v0. The acceleration is, a(t) = dr(t) dt = − v 2 0 R cos(v0t R )i − v 2 0 R sin(v0t R )j . Note that, the acceleration is perpendicular to the path (in this case it is parallel to r), since the velocity vector changes direction, but not magnitude. We can also verify that, from r(s), the unit tangent vector, et, could be computed directly as et = dr(s) ds = cos(v0t R )i + sin(v0t R )j . Example Motion along a helix The equation r(t) = R costi + R sin tj + htk, defines the motion of a particle moving on a helix of radius R, and pitch 2πh, at a constant speed. The velocity vector is given by v = dr dt = −R sin ti + R costj + hk , and the acceleration vector is given by, a = dv dt = −R costi + −R sin tj . In order to determine the speed at which the particle moves we simply compute the modulus of the velocity vector, v = |v| = p R2 sin2 t + R2 cos2 t + h 2 = p R2 + h 2 . If we want to obtain the equation of the path in terms of the arc-length coordinate we simply write, ds = |dr| = vdt = p R2 + h 2 dt . Integrating, we obtain s = s0 + √ R2 + h 2 t, where s0 corresponds to the path coordinate of the particle at time zero. Substituting t in terms of s, we obtain the expression for the position vector in terms of the arc-length coordinate. In this case, r(s) = R cos(s/√ R2 + h 2)i+R sin(s/√ R2 + h 2)j +hs/√ R2 + h 2k. The figure below shows the particle trajectory for R = 1 and h = 0.1. 1 0.5 0 0.5 1 1 0.5 0 0.5 1 0 0.5 1 1.5 2 6
Example Ballistic motion Consider the free-Hlight motion of a projectile which is initially launched with a velocity vo vo cos oi to sinoj. If we neglect air resistance, the only force on the projectile is the weight,which causes the projectile to have a constant acceleration a=-g3. In component form this equation can be written as dv / dt= 0 and doy/dt Integrating and imposing initial conditions, we get o, Uy=to sin- gt where we note that the horizontal velocity is constant. A further integration yields the trajectory (uo cos o)t, y=yo+(vo sin o niCi we recog nize as the equation of a parabola The maximum height, ymh, occurs when vy(tmh)=0, which gives tmh=(vo/g) sin o, or, The range, Ir, can be obtained by setting y= yo, which gives tr=(2vo/g)sin o, or sin o cos o sin(2o) We see that if we want to maximize the range ar, for a given velocity vo, then sin(2)=l, or =450 Finally, we note that if we want to model a more realistic situation and include aerodynamic drag forces of he form, say, -ku", then we would not be able to solve for a and y independently, and this would make the problem considerably more complicated(usually requiring numerical integration) ADDITIONAL READING J. L. Meriam and L.G. Kraige, Engineering Mechanics, DYNAMICS, 5th Edition 2/1,2/3,2/4
Example Ballistic Motion Consider the free-flight motion of a projectile which is initially launched with a velocity v0 = v0 cos φi + v0 sin φj. If we neglect air resistance, the only force on the projectile is the weight, which causes the projectile to have a constant acceleration a = −gj. In component form this equation can be written as dvx/dt = 0 and dvy/dt = −g. Integrating and imposing initial conditions, we get vx = v0 cos φ, vy = v0 sin φ − gt , where we note that the horizontal velocity is constant. A further integration yields the trajectory x = x0 + (v0 cos φ) t, y = y0 + (v0 sin φ) t − 1 2 gt2 , which we recognize as the equation of a parabola. The maximum height, ymh, occurs when vy(tmh) = 0, which gives tmh = (v0/g) sin φ, or, ymh = y0 + v 2 0 sin2 φ 2g . The range, xr, can be obtained by setting y = y0, which gives tr = (2v0/g) sin φ, or, xr = x0 + 2v 2 0 sin φ cos φ g = x0 + v 2 0 sin(2φ) g . We see that if we want to maximize the range xr, for a given velocity v0, then sin(2φ) = 1, or φ = 45o . Finally, we note that if we want to model a more realistic situation and include aerodynamic drag forces of the form, say, −κv2 , then we would not be able to solve for x and y independently, and this would make the problem considerably more complicated (usually requiring numerical integration). ADDITIONAL READING J.L. Meriam and L.G. Kraige, Engineering Mechanics, DYNAMICS, 5th Edition 2/1, 2/3, 2/4 7