Research Article | Open Access

Aleš Stanovnik, Borut Jurcic-Zlobec, "Numerical Study of the Elastic Pendulum on the Rotating Earth", *International Scholarly Research Notices*, vol. 2012, Article ID 806231, 7 pages, 2012. https://doi.org/10.5402/2012/806231

# Numerical Study of the Elastic Pendulum on the Rotating Earth

**Academic Editor:**G. Cleaver

#### Abstract

The elastic pendulum is a simple physical system represented by nonlinear differential equations. Analytical solutions for the bob trajectories on the rotating earth may be obtained in two limiting cases: for the ideally elastic pendulum with zero unstressed string length and for the Foucault pendulum with an inextensible string. The precession period of the oscillation plane, as seen by the local observer on the rotating earth, is 24 hours in the first case and has a well-known latitude dependence in the second case. In the present work, we have obtained numerical solutions of the nonlinear equations for different string elasticities in order to study the transition from one precession period to the other. It is found that the transition is abrupt and that it occurs for a quite small perturbation of the ideally elastic pendulum, that is, for the unstressed string length equal to about 10^{−4} of the equilibrium length due to the
weight of the bob.

#### 1. Introduction

The pendulum is such a fundamental physical system that all the details and various aspects of its motion are of interest. In the limit of inextensible suspension strings and small amplitudes, the behaviour of the system is described by simple, linear equations. However, a pendulum with an elastic suspension string is described by nonlinear differential equations, which couple horizontal and vertical oscillations [1]. This coupling has been mainly investigated in the vicinity of the so-called autoparametric resonance, where the string elasticity is such that the frequency of vertical, that is, spring-mode oscillations is double the frequency of horizontal, that is, pendulum-mode oscillations [1, 2].

It is perhaps surprising that no published work about the influence of earth's rotation on the behaviour of the nonlinear elastic pendulum could be found in the available literature. In the limit of an inextensible suspension string, one obtains the well-known and popular Foucault pendulum, for which the oscillation plane rotates, that is, precesses at a constant rate [3]. This rate is given by , where is earth's angular velocity, and is the geographical latitude of the site. The precession period would therefore be h. In the other limit of extreme elasticity, where the force of the suspension string is proportional to the total string length (i.e., the unstressed length is equal to zero), the equation of motion becomes linear and leads to a precession period equal to earth's rotation period, that is, to 24 hours [4]. The precession period is therefore an interesting parameter of the elastic pendulum, which undergoes a transition from a latitude-independent value of 24 h in the case of an extremely soft string to the latitude-dependent value of h in the case of a very stiff suspension string. In an attempt to learn more about this transition, we have performed numerical studies for various elasticities of the suspension string and for different latitudes of the pendulum site. Some results and basic findings are presented below.

#### 2. Method of Analysis

Instead of directly solving the equations of motion in the noninertial laboratory system by including the centrifugal and the Coriolis inertial forces, we have calculated the pendulum trajectories using a slightly different procedure. Such a procedure was followed also in deriving analytical solutions for the ideal elastic pendulum [4].

By neglecting the tidal gravitational forces of external bodies, such as the sun and moon, we have a quasi-inertial system , with origin at the earth's centre and with the orientation of coordinate axes fixed with respect to the distant stars (Figure 1). Let the axis of coincide with earth's rotation axis. In the present work, the equations for the elastic pendulum are solved numerically in the noninertial reference system , the origin of which is fixed at the point of pendulum suspension, and its coordinate axes remain parallel to the axes of the inertial system (Figure 1). The system is thus circulating, but it is not rotating. The differential equations are solved in , and the solutions are then transformed into the system of the local observer , where the axis points vertically upward, the axis points southward, and the axis points eastward (Figure 1). Since the inertial forces in the circulating system are less complicated than those in the circulating and rotating system , one might expect that also the trajectories are simpler in than in (see below). A famous historical example is given by the simple heliocentric planetary orbits as compared to their geocentric counterparts.

**(a)**

**(b)**

Neglecting dissipative forces, we have Newton's law for the motion of the pendulum bob where is the force of earth's gravity, is the force of the suspension string, is the mass of the bob, and is the bob acceleration in the inertial system . Since the coordinate axes of systems and remain parallel, the accelerations and , measured by observers in and , respectively, are simply related by , where is the acceleration of the origin of as measured by the observer in . Therefore, is the centripetal acceleration perpendicular to earth's rotation axis, that is, to the axis of .

We further assume that is of constant magnitude and has a fixed direction in the meridional plane containing the origin of . As both and lie in the meridional plane, by appropriately choosing , one obtains the effective gravity in as Here, , and, as already mentioned, is the angular velocity of the earth, and is the angle of with respect to the equatorial plane.

The force of the string is proportional to its dilatation and points towards the suspension point, that is, towards the origin of , where is the vector of the instantaneous bob position in , and is the unstressed length of the string. The elastic constant is then , where is the equilibrium length of the string, subject to the weight of the bob. The string force may therefore be written as By inserting all this into Newtons law (2.1) and rearranging the terms, we obtain a nonlinear differential equation for the position vector of the bob in , where is given by (2.2) above. The solution of the differential equation in , obtained for a given set of parameters, is then transformed into the system of the local observer. The transformation consists of a rotation around the axis for followed by a rotation around the axis for [4]. One has to transform also the initial conditions, because they are usually given in , but the equation is solved in .

Let us observe that the above nonlinear differential equation (2.5) reduces to a linear differential equation for the case of an ideal elastic suspension string, that is, for . The analytical solutions for such a case have been derived in [4].

From the differential equation (2.5), we see that the suspension string is determined by two parameters: the unstressed string length and the equilibrium length . The present calculations have been performed with a fixed equilibrium length m, but with the unstressed length varying between (ideally elastic string) and m (inextensible string). For a given mass of the bob, this corresponds to variation of the suspension string elastic constant. The two parameters, and , also determine the angular frequencies and of pendulum-mode and of spring-mode oscillations, respectively. For a fixed-string length , and for small amplitudes, we have the pendulum-mode angular frequency while the spring-mode angular frequency is given by It follows that the ratio of frequencies is determined by the ratio , The spring-mode frequency is thus always greater than the pendulum mode frequency and therefore determines the computation time. On the PC with an Intel Pentium, dual core, 3 GHz, and 64-bit processor, the Mathematica programme package requires between a few minutes and a few hours to calculate the bob trajectory during the first 48 hours of oscillation.

The motion of the bob has been studied for the usual initial conditions with which a Foucault pendulum is started in the system of the local observer Here, is the initial angle between the pendulum string and the negative axis, and is the initial angle between the oscillation plane and the plane. All the calculations have been performed with and with (displacement towards the south).

The oscillation plane is defined with the axis and with the vector of bob position at moments of maximum displacement. This amplitude vector is obtained with a special fitting procedure. The current maximum displacement vector therefore defines the oscillation plane, which is time dependent, that is, it precesses.

#### 3. Results

For various values of the string elasticity, we show in Figure 2 the azimuth angle of the oscillation plane (defined in the previous section) as a function of time at 55° geographical latitude. The string elasticity varied by varying the unstressed string length (), while keeping the equilibrium length (), due to the bob weight, at a constant value of 10 meters.

In the limit of a very stiff string, that is, for , we obtain the Foucault pendulum with a constant precession rate, which is represented by the straight line in Figure 2. For , the corresponding precession velocity is 12.3° per hour, so the precession period is 29.3 hours. For softer strings (), the oscillation plane initially keeps up with the Foucault pendulum, but after about 15 hours, it falls behind, resulting in a longer precession time. At string elasticity corresponding to about (see below), a discontinuity occurs. Instead of slowing down after 15 hours, the precession speeds up, completing a full circle in less time. As the string elasticity parameter is reduced towards zero, the precession time tends towards 24 hours.

In the limit of an ideally elastic string (), the equation of motion of the pendulum bob (2.5) becomes linear and may be solved analytically. As has been shown in a previous publication [4], the solution, as seen by an observer in , is quite simple. The trajectory is an ellipse with fixed orientation relative to the distant stars, but with the center of the ellipse moving on a circle. However, the local observer , with whom we are mainly concerned, will see a more complicated bob trajectory, which is obtained by a time-dependent transformation of the ellipse into the local system (with axis vertical, axis eastward, and axis southward). Consider, for example, an elliptical orbit, initially in the horizontal plane at latitude 45°. For the local observer in , the trajectory would be such that after 12 hours the ellipse will appear to be in a vertical plane. Alternatively, the relatively complicated bob trajectory as seen by the local observer may be attributed to the inertial forces (centrifugal and Coriolis) present in the noninertial system . However, it seems easier to understand the 24-hour period by considering the periodic transformation of the elliptical orbit from the system to the system . The addition of a small perturbation to the ideally elastic pendulum, in the form of a small nonzero unstressed string length , might be expected to contribute some modifications to the trajectory. Such trajectories therefore represent the faster branch in Figure 2.

In order to learn more about the transition from the 24-hour precession time to the precession time of the Foucault pendulum, we plot in Figure 3 the time required for the oscillation plane to return to its initial position and precession velocity as a function of string elasticity at different geographical latitudes. In Figure 2, for example, one obtains these precession times at the interception of the curves with a horizontal line going through 360°. (We note that for latitudes below 45° the oscillation plane does not necessarily complete a full rotation, but may oscillate between positive and negative azimuth angles). We see in Figure 3 that the transition between the two limiting values of precession time is abrupt and occurs at surprisingly low values of the elasticity parameter . The exact value of the elasticity at which the transition occurs seems to depend on latitude (Figure 3). This, as well as a dependence of the transition on initial conditions and on earth's angular velocity, remains to be the subject of further studies.

Finally, we should point out that in addition to the discontinuity in precession time at , our calculations show a structure at the elasticity given by (see Figure 3). This structure corresponds to the well-known parametric resonance for which (2.8).

#### 4. Conclusions

We have performed calculations of the precession of the elastic pendulum and have found a most interesting behaviour of the precession time as a function of suspension string elasticity. In the limit of an ideally elastic string with zero unstressed string length (), the precession is periodic with a period of 24 hours. Introduction of a small perturbation () causes a discontinuous transition from the 24-hour precession time to the latitude-dependent precession time of the Foucault pendulum ( h). In the broad range of between and , the precession time mainly has the value as for the inextensible Foucault pendulum. An exception is the structure observed at , which is due to the parametric resonance . We may therefore conclude that the agreement of our numerical results with limiting cases where analytical solutions exist, as well as with the well-known parametric resonance, allows for confidence in the correctness of the procedure and the results.

Surely, the reader is wondering why should there be a discontinuity of the pendulum precession time and why should it appear at that particular value of string elasticity. The authors do not have answers to these questions, but they do believe that answers will be found by further investigation of the long-term behaviour of the elastic pendulum. Maybe, this simple system can teach us something about more complex systems.

#### References

- E. Breitenberger and R. D. Mueller, “The elastic pendulum: a nonlinear paradigm,”
*Journal of Mathematical Physics*, vol. 22, no. 6, pp. 1196–1210, 1980. View at: Google Scholar - H. N. Núñez-Yépez, A. L. Salas-Brito, C. A. Vargas, and L. Vicente, “Onset of chaos in an extensible pendulum,”
*Physics Letters A*, vol. 145, no. 2-3, pp. 101–105, 1990. View at: Google Scholar - J. B. Marion,
*Classical Dynamics*, Harcourt Brace Jovanovich, New York, NY, USA, 1970. - A. Stanovnik, “The Foucault pendulum with an ideal elastic suspension string,”
*European Journal of Physics*, vol. 27, no. 2, pp. 205–213, 2006. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2012 Aleš Stanovnik and Borut Jurčič-Zlobec. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.