]> Oscilaciones de una partícula bajo la acción de dos muelles elásticos
Anterior

Solución analítica

La ecuación del movimiento es

m d 2 x d t 2 =2k( 1 l 0 x 2 + l 0 2 )x d 2 x d t 2 + 2k m x{ 1 ( 1+ ( x l 0 ) 2 ) 1/2 }

Cuando x<<l0 la ecuación diferencial se aproxima a

d 2 x d t 2 + k m l 0 2 x 3 =0

Las condiciones iniciales son: en el instante t=0, la partícula parte de la posición x=A, con velocidad nula, v=dx/dt=0

d 2 x d t 2 = dv dt = dv dx · dx dt =v dv dx v dv dx = k m l 0 2 x 3

Integramos esta ecuación diferencial

0 v vdv = k m l 0 2 A x x 3 dx v= k 2m l 0 2 ( A 4 x 4 )

es el módulo de la velocidad

Como la partícula parte de x=+A y se dirige hacia el origen la velocidad es negativa en el primer cuarto de periodo

v= k 2m l 0 2 ( A 4 x 4 )

Para obtener la posición x del oscilador en función del tiempo t tenemos que integrar nuevamente

dx dt = k 2m l 0 2 ( A 4 x 4 ) t= 2m l 0 2 k A x dx A 4 x 4

Hacemos el cambio de variable

x=A 1 z 2

La integral se transforma

t= 1 A 2m l 0 2 k 0 z dz 1 z 2 2 z 2 A l 0 k m t= 0 z dz 1 z 2 1 ( 2 2 ) 2 z 2

La integral elíptica primera especie se define

u=F(ϕ,k)= 0 ϕ dθ 1 k 2 sin 2 θ

Se hace el cambio de variable z=sinφ,

u= 0 z dz (1 z 2 )(1 k 2 z 2 )

se designa am u a la función inversa φ(u)

z=senφ=sin am u=sn u

Esta es la definición de la función elíptica sn.

Análogamente, escribiremos cos am u=cn u

Ambas funciones están relacionadas

sn2 u+cn2 u=1

En el caso que estamos analizando

u= A l 0 k m tk= 2 2

Conocido el valor de la variable auxiliar z=sn u, se despeja la posición x

x=A 1 sn 2 u=A·cn( A l 0 k m t )

Periodo del movimiento

El tiempo que tarda en describir un cuarto de periodo P es

P 4 = 2m l 0 2 k A 0 dx A 4 x 4

Haciendo el cambio de variable x=Acosφ

P 4 = l 0 A m k 0 π/2 dϕ 1 ( 2 2 ) 2 sin 2 ϕ

El último término es la integral elíptica completa de primera especie para k= 2 /2  que denominamos F( 2 /2 ) Su valor aproximado es 1.8541. Véase tabla de integrales elípticas de primera especie (Puig Adam P., Cálculo Integral. Editorial Biblioteca Matemática 1972, pág. 72)

P=4·F( 2 /2 ) l 0 A m k 7.4164 l 0 A m k

En la figura, se muestra cómo el periodo P depende de la amplitud A

Actividades

Se introduce

Se pulsa el botón titulado Inicio

Se pulsa el botón titulado Empieza

El programa interactivo calcula la posición x en función del tiempo t, mediante la función elíptica cn

LibresApplet aparecerá en un explorador compatible con JDK 1.1.

 

Referencias

Mohazzabi P. Theory and examples of intrinsically nonlinear oscillators. Am. J. Phys. 72 (4) April 2004, pp. 492-498

Detcheva V., Spassov V., A simple nonlinear oscillator: analytical amd numerical solution.  Phys. Educ. 28 (1993) pp. 39-42

Puig Adam P., Curso teórico práctico de Cálculo Integral aplicado a la Física y Técnica. Editorial Biblioteca Matemática 1972, págs. 71-76

Anterior