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 v = ˙s ≡ ds/dt. We also compute the rate of change of speed as at = ¨s = d2s/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 +mt2 = (rPx + mxt2)i+ (rPy + myt2)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) = t2, 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 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 = d2r 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. 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 = drdt = drds dsdt = 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) = ˙Ax(t)i+ ˙Ay(t)j+ ˙Az(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 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 = radicalBig v2x + v2y + v2z, and the magnitude of the acceleration is a = radicalBig a2x + a2y + a2z. 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) = Rcos( sR)i+ Rsin( sR)j . Since we have a constant speed v0, we have s = v0t. Thus, r(t) = Rcos(v0tR )i + Rsin(v0tR )j . The velocity is v(t) = dr(t)dt = ?v0 sin(v0tR )i+ v0 cos(v0tR )j , 5 which, clearly, has a constant magnitude |v| = v0. The acceleration is, a(t) = dr(t)dt = ?v 20 R cos( v0t R )i? v20 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(v0tR )i + sin(v0tR )j . Example Motion along a helix The equation r(t) = Rcosti+Rsintj +htk, defines the motion of a particle moving on a helix of radius R, and pitch 2pih, at a constant speed. The velocity vector is given by v = drdt = ?Rsinti+ Rcostj + hk , and the acceleration vector is given by, a = dvdt = ?Rcosti+?Rsintj . In order to determine the speed at which the particle moves we simply compute the modulus of the velocity vector, v = |v| = radicalbig R2 sin2 t + R2 cos2 t+ h2 = radicalbig R2 + h2 . If we want to obtain the equation of the path in terms of the arc-length coordinate we simply write, ds = |dr| = vdt = radicalbig R2 + h2 dt . Integrating, we obtain s = s0 + √R2 + h2 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) = Rcos(s/√R2 + h2)i+Rsin(s/√R2 + h2)j+hs/√R2 + h2k. 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-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? 12gt2 , 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 20 sin2 φ 2g . The range, xr, can be obtained by setting y = y0, which gives tr = (2v0/g)sinφ, or, xr = x0 + 2v 20 sinφcosφ g = x0 + v20 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