Abstracts
This paper deals with the nonlinear oscillation of a simple pendulum and presents not only the exact formula for the period but also the exact expression of the angular displacement as a function of the time, the amplitude of oscillations and the angular frequency for small oscillations. This angular displacement is written in terms of the Jacobi elliptic function sn(u;m) using the following initial conditions: the initial angular displacement is different from zero while the initial angular velocity is zero. The angular displacements are plotted using Mathematica, an available symbolic computer program that allows us to plot easily the function obtained. As we will see, even for amplitudes as high as 0.75pi (135º) it is possible to use the expression for the angular displacement, but considering the exact expression for the angular frequency w in terms of the complete elliptic integral of the first kind. We can conclude that for amplitudes lower than 135º the periodic motion exhibited by a simple pendulum is practically harmonic but its oscillations are not isochronous (the period is a function of the initial amplitude). We believe that present study may be a suitable and fruitful exercise for teaching and better understanding the behavior of the nonlinear pendulum in advanced undergraduate courses on classical mechanics.
simple pendulum; large-angle period; angular displacement
Este artigo aborda a oscilação não-linear de um pêndulo simples e apresenta não apenas a fórmula exata do período mas também a dependencia temporal do deslocamento angular para amplitudes das oscilações e a freqüência angular para pequenas oscilações. O deslocamento angular é escrito em termos da função elíptica de Jacobi sn(u;m) usando as seguintes condições iniciais: o deslocamento angular inicial é diferente de zero enquanto que a velocidade angular inicial é zero. Os deslocamentos angulares são plotados usando Mathematica, um disponível programa simbólico de computador que nos permite plotar facilmente a função obtida. Como veremos, mesmo para amplitudes tão grandes quanto 0,75pi (135º) é possível usar a expressão para o deslocamento angular mas considerando a expressão exata para a freqüência angular w em termos da integral elíptica completa de primeira espécie. Concluímos que, para amplitudes menores que 135º, o movimento periódico exibido por um pêndulo simples é praticamente harmônico, mas suas oscilações não são isócronas (o período é uma função da amplitude inicial). Acreditamos que o presente estudo possa tornar-se um exercício conveniente e frutífero para o ensino e para uma melhor compreensão do pêndulo não-linear em cursos avançados de mecânica clássica na graduação.
pêndulo simples; período a grandes ângulos; deslocamento angular
NOTAS E DISCUSSÕES
Exact solution for the nonlinear pendulum
Solução exata do pêndulo não linear
A. Beléndez1 1 E-mail: a.belendez@ua.es. ; C. Pascual; D.I. Méndez; T. Beléndez; C. Neipp
Departamento de Física, Ingeniería de Sistemas y Teoría de la Señal, Universidad de Alicante, Alicante, Spain
ABSTRACT
This paper deals with the nonlinear oscillation of a simple pendulum and presents not only the exact formula for the period but also the exact expression of the angular displacement as a function of the time, the amplitude of oscillations and the angular frequency for small oscillations. This angular displacement is written in terms of the Jacobi elliptic function sn(u;m) using the following initial conditions: the initial angular displacement is different from zero while the initial angular velocity is zero. The angular displacements are plotted using Mathematica, an available symbolic computer program that allows us to plot easily the function obtained. As we will see, even for amplitudes as high as 0.75p (135º) it is possible to use the expression for the angular displacement, but considering the exact expression for the angular frequency w in terms of the complete elliptic integral of the first kind. We can conclude that for amplitudes lower than 135º the periodic motion exhibited by a simple pendulum is practically harmonic but its oscillations are not isochronous (the period is a function of the initial amplitude). We believe that present study may be a suitable and fruitful exercise for teaching and better understanding the behavior of the nonlinear pendulum in advanced undergraduate courses on classical mechanics.
Keywords: simple pendulum, large-angle period, angular displacement.
RESUMO
Este artigo aborda a oscilação não-linear de um pêndulo simples e apresenta não apenas a fórmula exata do período mas também a dependencia temporal do deslocamento angular para amplitudes das oscilações e a freqüência angular para pequenas oscilações. O deslocamento angular é escrito em termos da função elíptica de Jacobi sn(u;m) usando as seguintes condições iniciais: o deslocamento angular inicial é diferente de zero enquanto que a velocidade angular inicial é zero. Os deslocamentos angulares são plotados usando Mathematica, um disponível programa simbólico de computador que nos permite plotar facilmente a função obtida. Como veremos, mesmo para amplitudes tão grandes quanto 0,75p (135º) é possível usar a expressão para o deslocamento angular mas considerando a expressão exata para a freqüência angular w em termos da integral elíptica completa de primeira espécie. Concluímos que, para amplitudes menores que 135º, o movimento periódico exibido por um pêndulo simples é praticamente harmônico, mas suas oscilações não são isócronas (o período é uma função da amplitude inicial). Acreditamos que o presente estudo possa tornar-se um exercício conveniente e frutífero para o ensino e para uma melhor compreensão do pêndulo não-linear em cursos avançados de mecânica clássica na graduação.
Palavras-chave: pêndulo simples, período a grandes ângulos, deslocamento angular.
Perhaps one of the nonlinear systems most studied and analyzed is the simple pendulum [1-12], which is the most popular textbook example of a nonlinear system and is studied not only in advanced but also in introductory university courses of classical mechanics. The periodic motion exhibited by a simple pendulum is harmonic only for small angle oscillations [1]. Beyond this limit, the equation of motion is nonlinear: the simple harmonic motion is unsatisfactory to model the oscillation motion for large amplitudes and in such cases the period depends on amplitude. Application of Newton's second law to this physical system gives a differential equation with a nonlinear term (the sine of an angle). It is possible to find the integral expression for the period of the pendulum and to express it in terms of elliptic functions. Although it is possible in many cases to replace the nonlinear differential equation by a corresponding linear differential equation that approximates the original equation, such linearization is not always feasible. In such cases, the actual nonlinear differential equation must be directly dealt with.
In this note we obtain analytical exact formulas for the period of a simple pendulum and we present some figures in which the angular displacement is plotted as a function of the time for different initial amplitudes. Moreover, we stress that one may show to the students the advantage of using available symbolic computer programs such as Mathematica.
As we can point out previously, one of the simplest nonlinear oscillating systems is the simple pendulum. This system consists of a particle of mass m attached to the end of a light inextensible rod, with the motion taking place in a vertical plane. The differential equation modelling the free undamped simple pendulum is
where q is the angular displacement, t is the time and w0 is defined as
Here l is the length of the pendulum and g is the acceleration due to gravity. Because of the presence of the trigonometric function sinq, Eq. (1) is a nonlinear differential equation.
We consider that the oscillations of the pendulum are subjected to the initial conditions
where q0 is the amplitude of the oscillation. The system oscillates between symmetric limits [-q0,+q0]. The periodic solution q(t) of Eq. (1) and the angular frequency w (also with the period T = 2p/w) depends on the amplitude q0.
Equation (1), although straightforward in appearance, is in fact rather difficult to solve because of the nonlinearity of the term sinq. In order to obtain the exact solution of Eq. (1), this equation is multiplied by dq / dt, so that it becomes
which can be written as
Equation (5), which corresponds to the conservation of the mechanical energy, is immediately integrable, taking into account initial conditions in Eq. (3). From Eq. (5) we can obtain
which can be written as follows
where we use the trigonometric relation
Now let
and
From Eqs. (3), (9) and (10) we have
It is easy to obtain the value of dq / dt as a function of dy / dt as follows. Firstly, from Eq. (10), we have
and secondly
Then, we have
Substituting Eqs. (9), (10) and (14) into Eq. (7), one get
which can be rewritten as follows
We do define new variables t and z as
Then Eq. (16) becomes
where 0 < k < 1, and
Solving Eq. (18) for dt gives
The time t to go from the point (1,0) to the point (z,dz/dt) in the lower half plane of the graph of dz/dt as a function of z is
Equation (21) can be rewritten as follows
which allows us to obtain t as a function of z and k as
where K(m) and F(j;m) are the complete and the incomplete elliptical integral of the first kind, defined as follows [13]
and z = sinj.
The period of oscillation T is four times the time taken by the pendulum to swing from q = 0 (z = 0) to q = q0 (z = 1). Therefore
where
is the period of the pendulum for small oscillations. Equation (23) can be written as follows
which can be written in terms of the Jacobi elliptic function sn(u;m) [14]
In terms of the original variables and taking into account Eqs. (9), (10), (17) and (26), Eq. (29) becomes
which allows us to express q as a function of t as
In Figs. 1-6 we have plotted the exact angular displacement q as a function of w0t (Eq. (31)) for different values of the initial amplitude q0. Plots have been obtained using the Mathematica program. In these figures we have also included the angular displacement q(t) = 2arcsin . As we can see, for amplitudes q0 < 0.75p (135º) it would be possible to use the approximate expression
for the angular displacement of the pendulum and Eq. (32) corresponds to a simple harmonic oscillator, but with the following exact expression for the angular frequency
However, for larger initial amplitudes it is easy to see that the angular displacement does not correspond to a simple harmonic oscillator and the exact expression must be considered.
We can conclude that for amplitudes as high as 135º the effect of the nonlinearity is seen only in the fact that the angular frequency of the oscillation w depends on the amplitude q0 of the motion. For these amplitudes the harmonic function provides an excellent approximation to the periodic solution of Eq. (1) and the periodic motion exhibited by a simple pendulum is practically harmonic but its oscillations are not isochronous (the period is a function of the amplitude of oscillations). Finally, we point out that in the process of plotting the exact displacement one may show to the students the advantage of using available symbolic computer programs such as Mathematica.
Acknowledgments
This work was supported by the "Ministerio de Educación y Ciencia", Spain, under project FIS2005-05881-C02-02, and by the "Generalitat Valenciana", Spain, under project ACOMP/2007/020.
Recebido em 30/7/2007; Aceito em 28/8/2007
- [1] J.B. Marion, Classical Dynamics of Particles and Systems (Harcourt Brace Jovanovich, San Diego, 1970).
- [2] L.P. Fulcher and B.F. Davis, Am. J. Phys. 44, 51 (1976).
- [3] W.P. Ganley, Am. J. Phys. 53, 73 (1985).
- [4] L.H. Cadwell and E.R. Boyco, Am. J. Phys. 59, 979 (1991).
- [5] M.I. Molina, Phys. Teach. 35, 489 (1997).
- [6] R.B. Kidd and S.L. Fogg, Phys. Teach. 40, 81 (2002).
- [7] L.E. Millet, Phys. Teach. 41, 162-3 (2003).
- [8] R.R. Parwani, Eur. J. Phys. 25, 37 (2004).
- [9] G.E. Hite Phys. Teach. 43, 290 (2005).
- [10] A. Beléndez, A. Hernández, A. Márquez, T. Beléndez and C. Neipp, Eur. J. Phys. 27, 539 (2006).
- [11] A. Beléndez, A. Hernández, T. Beléndez, C. Neipp and A. Márquez, Eur. J. Phys. 28, 93 (2007).
- [12] R.E. Mickens, Oscillations in Planar Dynamics Systems (World Scientific, Singapore,1996).
- [13] L.M. Milne-Thomson, in M. Abramowitz and I.A. Stegun (eds) Handbook of Mathematical Functions (Dover Publications, Inc., Nova Iorque, 1972).
- [14] L.M. Milne-Thomson in M. Abramowitz and I.A. Stegun (eds) Handbook of Mathematical Functions (Dover Publications, Inc., Nova Iorque, 1972).
Publication Dates
-
Publication in this collection
18 Mar 2008 -
Date of issue
2007
History
-
Accepted
28 Aug 2007 -
Received
30 July 2007