Numerical Study of the Elastic Pendulum on the Rotating Earth

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.


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 ω sin λ, where ω is earth's angular velocity, and λ is the geographical latitude of the site.The precession period would therefore be T 24/ sin λ 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 24/ sin λ 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.

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 S, 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 z axis of S coincide with earth's rotation axis.In the present work, the equations for the elastic pendulum are solved numerically in the noninertial reference system S , 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 S Figure 1 .The system S is thus circulating, but it is not rotating.The differential equations are solved in S , and the solutions are then transformed into the system of the local observer S , where the z axis points vertically upward, the x axis points southward, and the y axis points eastward Figure 1 .Since the inertial forces in the circulating system S are less complicated than those in the circulating and rotating system S , one might expect that also the trajectories are simpler in S than in S see below .A famous historical example is given by the simple heliocentric planetary orbits as compared to their geocentric counterparts.
Neglecting dissipative forces, we have Newton's law for the motion of the pendulum bob where F g is the force of earth's gravity, F s is the force of the suspension string, m is the mass of the bob, and a is the bob acceleration in the inertial system S. Since the coordinate axes of systems S and S remain parallel, the accelerations a and a , measured by observers in S and S , respectively, are simply related by a a 0 a , where a 0 is the acceleration of the origin of S as measured by the observer in S. Therefore, a 0 is the centripetal acceleration perpendicular to earth's rotation axis, that is, to the z axis of S.
We further assume that F g mg 0 is of constant magnitude and has a fixed direction in the meridional plane containing the origin of S .As both g 0 and a 0 lie in the meridional plane, by appropriately choosing t 0, one obtains the effective gravity in S as g ef g 0 − a 0 −g ef cos λ cos ωt , cos λ sin ωt , sin λ .

2.2
Here, g ef |g 0 − a 0 |, and, as already mentioned, ω is the angular velocity of the earth, and λ is the angle of g ef 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 S , where r t is the vector of the instantaneous bob position in S , and l 0 is the unstressed length of the string.The elastic constant is then k mg ef / l e − l 0 , where l e is the equilibrium length of the string, subject to the weight mg ef 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 r t of the bob in S , where g 0 − a 0 is given by 2.2 above.The solution of the differential equation in S , obtained for a given set of parameters, is then transformed into the system S of the local observer.The transformation consists of a rotation around the z axis for ωt followed by a rotation around the y axis for π/2 − λ 4 .One has to transform also the initial conditions, because they are usually given in S , but the equation is solved in S .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 l 0 0. 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 l 0 and the equilibrium length l e .The present calculations have been performed with a fixed equilibrium length l e 10 m, but with the unstressed length l 0 varying between l 0 0 ideally elastic string and l 0 l e 10 m inextensible string .For a given mass of the bob, this corresponds to variation of the suspension string elastic constant.The two parameters, l 0 and l e , also determine the angular frequencies ω p and ω s of pendulum-mode and of spring-mode oscillations, respectively.For a fixed-string length l e , and for small amplitudes, we have the pendulum-mode angular frequency ω p g ef l e , 2.6 while the spring-mode angular frequency is given by ω s k m g ef l e − l 0 .

2.7
It follows that the ratio ω s /ω p of frequencies is determined by the ratio l 0 /l e , ω s ω p l e l e − l 0 .

2.8
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 S of the local observer r 0 l e sin θ 0 cos φ 0 , sin θ 0 sin φ 0 , − cos θ 0 , v 0 0, 0, 0 .

2.9
Here, θ 0 is the initial angle between the pendulum string and the negative z axis, and φ 0 is the initial angle between the oscillation plane and the y 0 plane.All the calculations have been performed with θ 0 5 • and with φ 0 0 • displacement towards the south .The oscillation plane is defined with the z axis and with the vector of bob position r 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.

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 l 0 , while keeping the equilibrium length l e , due to the bob weight, at a constant value of 10 meters.
In the limit of a very stiff string, that is, for l 0 l e , we obtain the Foucault pendulum with a constant precession rate, which is represented by the straight line in Figure 2.For λ 55 • , the corresponding precession velocity is 12.3 • per hour, so the precession period is 29.3 hours.For softer strings l 0 < l e , 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 l 0 ≈ 10 −4 l e 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 l 0 is reduced towards zero, the precession time tends towards 24 hours.
In the limit of an ideally elastic string l 0 0 , 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 S , 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 S , 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 S with z axis vertical, y axis eastward, and x axis southward .Consider, for example, an elliptical orbit, initially in the horizontal plane at latitude 45 • .For the local observer in S , 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 S may be attributed to the inertial forces centrifugal and Coriolis present in the noninertial system S .However, it seems easier to understand the 24-hour period by considering the periodic transformation of the elliptical orbit from the system S to the system S .The addition of a small perturbation to the ideally elastic pendulum, in the form of a small nonzero unstressed string length l 0 , 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 l 0 ≈ 10 −4 l e .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 l 0 /l e ≈ 10 −4 , our calculations show a structure at the elasticity given by l 0 /l e 3/4 see Figure 3 .This structure corresponds to the well-known parametric resonance for which ω s /ω p 2 2.8 .

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 l 0 0 , the precession is periodic with a period of 24 hours.Introduction of a small perturbation l 0 10 −4 l e causes a discontinuous transition from the 24-hour precession time to the latitude-dependent precession time of the Foucault pendulum T 24/ sin λ h .In the broad range of l 0 between 10 −4 l e and l e , the precession time mainly has the value as for the inextensible Foucault

Figure 1 :
Figure 1:Views along the y, y axes a and z, z axes b of reference systems S and S , which have parallel axes, that are fixed relative to the distant stars.The origin of S is at earth's center, and the origins of S and S are at the pendulum suspension point.The local lab system S has the z axis vertical, the x axis southward, and the y axis eastward.

Figure 2 :
Figure 2: The dependence of the azimuth angle of the oscillation plane on time for latitude 55 • .Different curves represent different string elasticities see text .

Figure 3 :
Figure 3:The dependence of pendulum precession time on the suspension string elasticity, that is, on the length l 0 of the unstressed string at fixed equilibrium length l e 10 m.The curves connect points corresponding to a given geographical latitude.