J. Peraire 16.07 Dynamics Fall 2004 Version 1.1 Lecture D28 - Central Force Motion: Kepler’s Laws When the only force acting on a particle is always directed to- wards a fixed point, the motion is called central force motion. This type of motion is particularly relevant when studying the orbital movement of planets and satellites. The laws which gov- ern this motion were first postulated by Kepler and deduced from observation. In this lecture, we will see that these laws are a con- sequence of Newton’s second law. An understanding of central force motion is necessary for the design of satellites and space vehicles. Kepler’s Problem We consider the motion of a particle of mass m, in an inertial reference frame, under the influence of a force, F, directed towards the origin. We will be particularly interested in the case when the force is inversely proportional to the square of the distance between the particle and the origin, such as the gravitational force. In this case, F = ?μr2mer, where μ is the gravitational parameter, r is the modulus of the position vector, r, and er = r/r. It can be shown that, in general, Kepler’s problem is equivalent to the two-body problem, in which two masses, M andm, move solely due to the influence of their mutual gravitational attraction. This equivalence is obvious when M >> m, since, in this case, the center of mass of the system can be taken to be at M. However, even in the more general case when the two masses are of similar size, the problem can be reduced to a Kepler problem (see note below). 1 Although most problems in celestial mechanics involve more than two bodies, many problems of practical interest can be accurately solved by just looking at two bodies at a time. When more than two bodies are involved, the problem is considerably more complicated, and, in this case, no general solutions are known. This problem was studied by Kepler (1571-1630) who lived before Newton was born. His interest was in describing the motion of planets around the sun. He postulated the following laws: 1.- The orbits of the planets are ellipses with the Sun at one focus 2.- The line joining a planet to the Sun sweeps out equal areas in equal intervals of time 3.- The square of the period of a planet is proportional to the cube of the major axis of its elliptical orbit In this lecture, we will start from Newton’s laws and verify that the above three laws can indeed be derived from Newtonian mechanics. Note Equivalence between the two-body problem and Kepler’s problem Here we consider the problem of two isolated bodies of masses M and m which interact though gravitational attraction. Let rM and rm denote the position vectors of the two bodies relative to a fixed origin O. Thus, the position of the center of gravity, G, of the two bodies will be rG = MrM +mrmM +m . Since the two bodies are isolated, we will have, from momentum conservation, that ˙rG =constant, and ¨rG = 0. Therefore, the position of the center of gravity, at all times, can be found trivially from the initial conditions. If the position vector of m as observed by M, r = rm ?rM, is known, then the position vectors of M and m could be computed as rM = rG ? mM +mr, rm = rG + MM +mr . Therefore, the problem of determining rM and rm is equivalent to that of determining r. In order to derive an equation for r, we first consider the equations of motion for m and M independently. Thus, we have M¨rM = GMmr2 er, m¨rm = ?GMmr2 er , 2 where r = |r|, er = r/r, and G is the gravitational constant. Subtracting these two expressions, we obtain, ¨r = ¨rm ? ¨rM = ?GM +mr2 er , or, Mm M +m¨r = ?G Mm r2 er . The above expression shows that the motion of m relative to M is in fact a Kepler problem in which the force is given by ?GMmer/r2 (this is indeed the real force), but the mass of the orbiting body (m in this case), has been replaced by the reduced mass, Mm/(M +m). Note that when M >>m, the reduced mass becomes m. However, the above expression is general and applies to general masses M and m. Alternatively, the above expression can be written as m¨r = ?G(M +m)mr2 er , which is again a Kepler problem for an orbiting body of mass M, in which the gravitational parameter μ is given by μ = G(M +m). Equations of Motion The equation of motion (F = ma), is ?μmr2 er =m¨r. Since the only force in the system is directed towards point O, the angular momentum of m with respect to the origin will be constant. Therefore, the position and velocity vectors, r and ˙r, will be in a plane orthogonal to the angular momentum vector, and, as a consequence, the motion will be planar. Using cylindrical coordinates, with ez being parallel to the angular momentum vector, we have, ?μr2er = (¨r?r˙θ2)er + (r¨θ+ 2˙r˙θ)eθ. Now, we consider the radial and circumferential components of this equation separately. Circumferential component We have, 0 =r¨θ+ 2˙r˙θ . Using the following identity, 1 r parenleftbigg d dt(r 2 ˙θ) parenrightbigg = r¨θ+ 2˙r˙θ, 3 the above equation implies that r2 ˙θ = h≡ constant. (1) We note that the constant of integration, h, that will be determined by the initial conditions, is precisely the magnitude of the specific angular momentum vector, i.e. h = |r×v|. In a time dt, the area, dA, swept by r will be dA =rrdθ/2. Therefore, dA dt = 1 2r 2 ˙θ = h 2, which proves Kepler’s second law. Radial component The radial component of the equation of motion reads, ?μr2 = ¨r?r˙θ2. (2) Since ?r2 ddt parenleftbig1rparenrightbig = ˙r, and r2 =h/˙θ, from equation 1, we can write ˙r = ?h˙θ ddt parenleftbigg1 r parenrightbigg = ?h ddθ parenleftbigg1 r parenrightbigg . Differentiating with respect to time, ¨r = ?h d 2 dθ2 parenleftbigg1 r parenrightbigg ˙θ = ?h2 r2 d2 dθ2 parenleftbigg1 r parenrightbigg . Inserting this expression into equation 2, and using equation 1, we obtain the following differential equation for 1/r as a function of θ. d2 dθ2 parenleftbigg1 r parenrightbigg + 1r = μh2. This is a linear second order ordinary differential equation which has a general solution of the form, 1 r = μ h2(1 +ecos(θ+ψ)) , where e and ψ are two constants of integration. If we choose θ to be zero when r is minimum, then e will be positive, and ψ = 0. The equation describing the trajectory will be r = h 2/μ 1 +ecosθ . (3) We shall see below that this is the equation of a conic section in polar coordinates. 4 Conic Sections Conic sections are planar curves that are defined as follows: given a line, or directrix, and a point, or focus O, a conic section is the locus of points, P, such that the ratio of the distance between the point and the focus, PO, to the distance between the point and the directrix, PA, is a constant e. That is, e = PO/PA. Since, PO = r, and PA = p/e?rcosθ, we have r = p1 +ecosθ . (4) Here, p is the so called parameter of the conic and is equal to r when θ = ±90o. The constant e ≥ 0 is called the eccentricity, and, depending on its value, the conic surface will be either an open or closed curve. In particular, we have that when e= 0 the curve is a circle e< 1 the curve is an ellipse e= 1 the curve is a parabola e> 1 the curve is a hyperbola. Comparing equations 4 and 3, we see that the trajectory of a mass under the influence of a central force will be a conic curve with parameter p = h2/μ. When e< 1, the trajectory is an ellipse, thus proving Kepler’s first law. The point in the trajectory which is closest to the focus is called the periapsis and is denoted by pi. For elliptical orbits, the point in the trajectory which is farthest away from the focus is called the apoapsis and is denoted by α. When considering orbits around the earth, these points are called the perigee and apogee, whereas for orbits around the sun, these points are called the perihelion and aphelion, respectively. Elliptical Trajectories If a is the semi-major axis of the ellipse, then 2a = rpi + rα, using equation 4 to evaluate rpi and rα, we obtain a = p/(1?e2). Thus, rpi = p1 +e = a(1?e), rα = p1?e = a(1 +e). 5 Also, the distance between O and the center of the ellipse will be a?rpi =ae. The distance between point D and the directrix will be equal to DO/e, which in turn will be equal to the sum of the distance between the focus and the center of the ellipse, plus the distance between the focus and the directrix. That is, DO/e =ae+p/e. Therefore, DO = ae2 +p =a. Hence, using Pythagoras’ theorem, b2 + (ae)2 =a2, the semi-minor axis of the ellipse will be b = a√1?e2. The area of the ellipse is given by A = piab. Also, since dA/dt =h/2 is a constant, we have A = hT/2, where T is the period of the orbit. Equating these two expressions and expressing h in terms of the semi-major axis as h2 = μp= μa(1?e2), we have μ = parenleftbigg2pi T parenrightbigg2 a3, (5) which proves Kepler’s third law. ADDITIONAL READING J.L. Meriam and L.G. Kraige, Engineering Mechanics, DYNAMICS, 5th Edition 3/13 (except energy analysis) 6