Departamento de Matemáticas y Física

ITESO

PERÍODO DE OSCILACIÓN DE UN PÉNDULO SIMPLE

Fís. Jorge Eduardo Aguilar Rosas

Introducción.

Generalmente, al estudiar sistemas oscilatorios nos gustaría que el movimiento fuera un

Movimiento Armónico Simple, como en el caso de un sistema Masa-Resorte; sin embargo, existen
sistemas en donde las oscilaciones no tienen un Movimiento Armónico Simple. Como un ejemplo
de esta situación consideraremos en esta ocasión al sistema conocido como Péndulo Simple, y
realizaremos el análisis para determinar el período de oscilación correspondiente.

1. Ecuación de Movimiento.

El Péndulo Simple está compuesto de una masa “pequeña” m que está unida a un extremo de una
cuerda “ligera”, de longitud l, que tiene el otro extremo fijo, como se indica en la figura 1. Al decir
que la masa es “pequeña” se considera que sus dimensiones son pequeñas de tal manera que no
es necesario tomar en cuenta su momento de inercia respecto a su centro de masa, y así poder
tratarla como una masa “puntual”; mientras que al decir que la cuerda es “ligera” se refiere a que su
masa es despreciable comparada con la masa suspendida, de tal forma que se puede omitir sus
efectos en el movimiento del Péndulo.

θ
l

m
mg

Figura 1. Fuerzas que actúan sobre la masa de un péndulo simple.

La posición de la masa en la parte más baja corresponde a la posición de equilibrio, de tal manera
que si es llevada fuera de esta posición, las fuerzas que actúan sobre la masa tratarán de llevarla
hacia la posición de equilibrio. Sea θ0 el ángulo de amplitud máxima, y θ el ángulo para cualquier
tiempo. Las fuerzas que actúan sobre la masa son: el peso, mg; y, la tensión, T. La componente
tangencial de la ecuación de movimiento, de acuerdo a la Segunda Ley de Newton, es:

−mgsen θ = ma T ,
Periodo de Oscilación de un Jorge Eduardo Aguilar Rosas
Péndulo Simple Departamento de Matemáticas y Física
ITESO

donde aT es la aceleración tangencial, dada por el producto del radio de giro l y la aceleración
angular, de tal manera que la ecuación de movimiento se puede escribir como:

d2θ g
2
+ senθ = 0. (1)
dt l

La ecuación de movimiento no corresponde a la de un Movimiento Armónico Simple, que tiene una

forma general dada por:

d2 ξ 2
2
+ ω0 ξ = 0, (2)
dt

donde ξ es la variable de posición, y ω0 es la frecuencia angular de oscilación.

2. Pequeñas Oscilaciones.

Generalmente en los movimientos oscilatorios que no son armónicos simples se busca determinar
las condiciones bajo las cuales el movimiento se pudiera considerar armónico simple. En este caso
al comparar la ecuación de movimiento (ec. 1), con la ecuación general de un Movimiento
Armónico Simple, es fácil ver que si

senθ ≈ θ, (3)

la ecuación de movimiento queda en la forma de un Movimiento Armónico Simple,

d2 θ g
2
+ θ = 0. (4)
dt l

La condición 3 indica se refiere a ángulos pequeños, por lo que la condición señala que se trata de
“pequeñas” oscilaciones. Comparando las ecuaciones 2 y 4 se puede identificar a la frecuencia
angular

1/ 2
g
ω0 = . (5)
l

Como la frecuencia angular es igual al cociente de 2π entre el período, T0, entonces tenemos que
el período para pequeñas oscilaciones es:

1/ 2
l
T0 = 2π . (6)
g

3. Período de Oscilación en General.

Para determinar el período de oscilación del péndulo simple en general comenzaremos por escribir
la ecuación de movimiento (ec. 1) en la forma:

2
Jorge Eduardo Aguilar Rosas Periodo de Oscilación de un
Departamento de Matemáticas y Física Péndulo Simple
ITESO

dω 2
= −ω 0 senθ.
dt

siendo ω la velocidad angular. Multiplicando la ecuación por dθ, tenemos

dω 2
dθ = −ω0 senθdθ
dt
2
ωdω = −ω0 senθdθ.

Si realizamos la integración en el intervalo de tiempo desde t0 hasta t, en donde la posición angular

cambia desde θ0 hasta θ, y la velocidad angular pasa de cero a ω, tenemos:

ω θ
2
ωdω = −ω0 senθdθ
0 θ0

ω2
= ω 0 (cos θ − cos θ 0 ),
2
2

de donde se tiene a velocidad angular

dθ
= 2ω0 (cos θ − cos θ 0 ) .
1/ 2
ω=
dt

Esta ecuación es una relación de la posición angular y el tiempo. Despejando las variables en la
igualdad, e integrando de nuevo en el tiempo desde t1 hasta t2, en donde la posición angular
cambia desde θ1 hasta θ2, obtenemos:

t2 θ2
dθ
dt = .
2ω0 (cos θ − cos θ 0 )
1/ 2
t1 θ1

Si se considera que θ1 = 0 y que θ2 = θ0, con t1 = 0 se tiene que t2 = T/4, siendo T el período,
entonces se tiene que

θ0
4 dθ
T= . (7)
2ω0 0 (cos θ − cos θ 0 )1/ 2
La integral a resolver es una Integral Elíptica de primera clase, cuya solución se puede obtener en
forma numérica o en series. Para obtener la solución de la integral en series consideremos el
cambio de variable de ángulo medio:

θ θ
sen = sen 0 senβ. (8)
2 2

Cuando θ es igual a 0, β es 0 también; mientras que cuando θ es igual a θ0, β es π/2.

Utilizando identidades trigonométricas y el cambio de variable obtenemos que

3
Periodo de Oscilación de un Jorge Eduardo Aguilar Rosas
Péndulo Simple Departamento de Matemáticas y Física
ITESO

θ0
cos θ = 1 − 2sen 2 sen 2 β; (9)
2

al evaluar en θ0 nos da que

θ0
cos θ 0 = 1 − 2sen 2 . (10)
2

Con las expresiones 9 y 10 el denominador del integrando resulta:

θ0
[cos θ − cos θ 0 ]1/ 2 = 2sen cos β. (11)
2

Por otra parte, tomando el diferencial de la relación del cambio de variable, llegamos a la expresión
para el diferencial del ángulo,

θ0
2sen cos β
2
dθ = 1/ 2
dβ (12)
θ0
1 − sen 2 sen 2 β
2

Así la integral de la ecuación 7 se transforma de la siguiente manera:

θ0
θ0 π/2 2sen cos βdβ
dθ 2
=
0 [cos θ − cos θ 0 ]
1/ 2 1/ 2
0 2 θ0 θ0
1 − sen sen 2 β 2sen cos β
2 2

θ0 π/2
dθ 2dβ
= (13)
0 [cos θ − cos θ 0 ] 1/ 2
0 2 θ0
1/ 2
1 − sen sen 2 β
2

El desarrollo en serie binomial del término del denominador del integrando es:

[1 − q sen β]
2 −1 / 2
= 1+
−1
2
(
− q sen 2 β +
−1
2
) −3
2
1
2
(
− q sen 2 β
2
)
+
−1
2
−3
2
−5
2
1
(2)(3 )
( 3
− q sen 2 β + ... )

[1 − qsen β] 2 −1/ 2
=
∞
(2n − 1)!! qn sen 2nβ, (14)
n=0 (2n )!!

4
Jorge Eduardo Aguilar Rosas Periodo de Oscilación de un
Departamento de Matemáticas y Física Péndulo Simple
ITESO

2 1
siendo q = sen (θ0/2); y, el símbolo “!!”, la manera de denotar al doble factorial . Sustituyendo el
desarrollo en la integral de la ecuación 13, tenemos:

θ0
dθ
= 2
∞
(2n + 1)!!qn π / 2sen 2nβdβ. (15)
0 [cos θ − cos θ 0 ]1/ 2 n=0 (2n )!! 0

2
La integral sobre β que aparece en cada uno de los términos de la suma es :

1
El doble factorial de (2n-1) es el producto de (2n-1) y los números impares menores a (2n-1); mientras que el doble
factorial de (2n) es el producto de (2n) y los números pares menores a 2n. En particular, (0)!! = 1, y (-1)!! = 1.
2
Podemos mostrar que el resultado de la integral es el indicado en la ecuación 20 mediante inducción.

(i) Si n = 0, la integral resulta:

π/2
π
1dβ = ,
0
2

Por lo que la ecuación 16 es válida para n = 0.

(ii) Supongamos que es válida para n = m-1, es decir que

π/2
sen 2(m−1)β dβ =
π (2(m − 1) − 1)!! ,
0
2 (2(m − 1))!!

y mostremos que también es válida para n = m. La integral para n = m es:

π/2 π/2
sen 2m βdβ = sen 2(m−1)β sen 2 βdβ
0 0
π/2 π/2
= sen 2(m−1)βdβ − sen 2(m−1)β cos 2 βdβ.
0 0

La primera de las integrales corresponde a la integral para n = m-1, y en la segunda integral realizamos una integración
por partes con u = cosβ, y dv = sen2(m-1)βcosβdβ, quedando

π/2
sen 2m βdβ =
π (2(m − 1) − 1)!! − 1 π / 2sen 2mβdβ.
0
2 (2(m − 1))!! (2m − 1) 0

De ahí, factorizando la integral tenemos:

1+
1
π/2
sen 2m β dβ =
π (2(m − 1) − 1)!! .
(2m − 1) 0
2 (2(m − 1))!!

Finalmente, despejando la integral obtenemos:

π/2
sen 2m β dβ =
π (2(m − 1) − 1)!! (2m − 1)
0
2 (2(m − 1))!! (2m)
π/2
sen 2m β dβ =
π (2m − 1)!! ,
0
2 (2m))!!

con lo queda mostrado que la ecuación 16 también es valida para n = m.

5
Periodo de Oscilación de un Jorge Eduardo Aguilar Rosas
Péndulo Simple Departamento de Matemáticas y Física
ITESO

π/2
sen 2n βdβ =
π (2n − 1)!! . (16)
0
2 (2n)!!

Sustituyendo el resultado 16, y la equivalencia de q en la ecuación 15, para luego ponerlo en la

ecuación 7, tenemos que el periodo de oscilación queda escrito de la forma:

T=
4
2
π ∞
(2n + 1)!! 2
sen 2n
θ0
2ω0 2 n= 0 (2n)!! 2

T = T0
∞
(2n + 1)!! 2
sen 2n
θ0
, (17)
n= 0 (2n)!! 2

en donde se ha considerado a T0 como el período que se obtuvo para pequeñas oscilaciones (ec.
6). Considerando los primeros términos de las series tenemos:

1 θ 9 θ 225 θ
T = T0 1 + sen 2 0 + sen 4 0 + sen 6 0 + ... (18)
4 2 64 2 2304 2

Si sólo se considera el primer término del desarrollo se tiene la aproximación para “pequeñas”
oscilaciones; pero, conforme se toman amplitudes “grandes”, la presencia de términos de mayor
orden debe de ser considerada. Si suponemos que el segundo término es un 1.0% respecto al
primero, tenemos que
1 θ
sen 2 0 = 0.01,
4 2

0
de donde se obtiene que la amplitud máxima es de θ0 = 23.1 , y el tercer término para este ángulo
es de 0.000225 (es decir, de 0.0225% respecto al primero). Entonces se puede considerar que
0
para amplitudes máximas menores a 23.0 se tiene un Movimiento Armónico Simple, con una
precisión del orden del 1.0%. La figura 2 muestra la gráfica del porcentaje de la desviación del
periodo respecto al periodo para “pequeñas” oscilaciones en función del ángulo de amplitud
máxima.

5
P o rc. d e De sv . T T o T o %

0
0 0.2 0.4 0.6 0.8 1
Ángulo rad
Figura 2. Gráfica del período en función del ángulo de amplitud máxima.

6
Jorge Eduardo Aguilar Rosas Periodo de Oscilación de un
Departamento de Matemáticas y Física Péndulo Simple
ITESO

4. Posición en Función del Tiempo.

La solución de la ecuación de movimiento para obtener la posición angular en términos del tiempo
se puede obtener numéricamente y compararla con la función θ0cos(2πt/T), en donde T es el
período, para distintos valores de la amplitud máxima θ0. Las figuras 3 a 4 muestran las gráficas de
las funciones para los valores de amplitud máxima de π/2, y π/4. En la gráfica en donde la amplitud
máxima es de π/2 es clara la diferencia entre la solución numérica y la función θ0cos(2πt/T);
mientras que para la amplitud máxima de π/4 la diferencia no es notoria.

0.8
1.5

1.25
0.6
P o si ci o n A n g u l a r ra d

P o si ci o n A n g u l a r ra d
1

0.75 0.4

0.5
0.2
0.25

0
0
0 0.05 0.1 0.15 0.2 0.25 0.3
Tiempo t s 0 0.05 0.1 0.15 0.2 0.25
Tiempo t s
Figura 3. Gráfica de la posición angular en función del Figura 4. Gráfica de la posición angular en función del
tiempo, y de la función θ0cos(2πt/T) (línea punteada), en tiempo, y de la función θ0cos(2πt/T) (línea punteada), en
el intervalo de tiempo cero a T/4, con la amplitud máxima el intervalo de tiempo cero a T/4, con la amplitud máxima
de π/2. de π/4.