Oscilaciones libres

Cuando una partícula se desplaza x de la posición de equilibrio, actúa sobre ella una fuerza que es proporcional al desplazamiento x y de sentido contrario a éste, tal como se muestra en la figura.

La ecuación del movimiento se escribe

ma=-kx

Teniendo en cuenta que la aceleración es la derivada segunda de la posición x, expresamos la ecuación del movimiento como ecuación diferencial de segundo orden.

d 2 x d t 2 + ω 0 2 x=0   ω 0 2 = k m

ω0 se denomina frecuencia propia o natural del oscilador armónico.

La ventaja de expresar las oscilaciones en términos de una ecuación diferencial es que podemos establecer analogías entre sistemas físicos oscilantes completamente diferentes: mecánicos, eléctricos, hidráulicos, etc.

La solución de esta ecuación diferencial es

x=Asin( ω 0 t)+Bcos( ω 0 t) v= dx dt = ω 0 ( Acos( ω 0 t)−Bsin( ω 0 t) )

Las constantes A y B se determinan a partir de la posición inicial x0 y la velocidad inicial v0 en el instante t=0.

t=0 { x 0 =B v 0 = ω 0 A x= v 0 ω 0 sin( ω 0 t)+ x 0 cos( ω 0 t)

Sea un oscilador de frecuencia angular ω0=100 rad/s. Sabiendo que la partícula parte de la posición x0=5 con velocidad inicial nula, v0=0, obtener y representar su posición x en función del tiempo t

O bien, de forma equivalente, x=A·sin(ω0t+φ), donde A es la amplitud y φ la fase inicial

A= v 0 2 ω 0 2 + x 0 2  tanφ= x 0 ω 0 v 0

>> syms w0 x0 v0;
>> x=dsolve('D2x+w0^2*x=0','x(0)=x0','Dx(0)=v0');
>> x=simplify(x)
x =x0*cos(t*w0) + (v0*sin(t*w0))/w0
>> xx=subs(x,{w0 x0 v0},{100 5 0})
xx =5*cos(100*t)
>> ezplot(xx,[0 0.1*pi])
>> xlabel('t')
>> ylabel('x')
>> title('oscilaciones libres')
>> grid on

Trayectoria en el espacio de las fases

El espacio de las fases nos muestra otra perspectiva del comportamiento de un oscilador. Se representa el momento lineal (o la velocidad) v en el eje vertical, y la posición del móvil x en el eje horizontal.

x=A·sin(ω0t+φ)
v=A·ω0·cos(ω0t+φ)

Eliminando el tiempo t en estas dos ecuaciones, obtenemos la ecuación de la trayectoria, una elipse.

x 2 A 2 + y 2 (A ω 0 ) 2 =1

Energía del oscilador

La fuerza es proporcional al desplazamiento x y de sentido contrario a éste.

F=−m ω 0 2 x

Como la fuerza F es conservativa. El trabajo de dicha fuerza es igual a la diferencia entre el valor inicial y el final de la energía potencial Ep.

∫ x 2 x 1 F·dx = E p1 − E p2 ∫ x 2 x 1 F·dx = ∫ x 2 x 1 −m ω 0 2 x·dx =− 1 2 m ω 0 2 x 2 | x 1 x 2 = 1 2 m ω 0 2 x 1 2 − 1 2 m ω 0 2 x 2 2

La expresión de la energía potencial es

E p (x)= 1 2 m ω 0 2 x 2 +c

Donde c es cualquier constante. Se toma como nivel cero de la energía potencial Ep=0 cuando el móvil está en el origen, x=0, por lo que c=0

La función Ep(x) representa una parábola cuyo vértice está en el origen, que tiene un mínimo en x=0 cuyo valor es Ep=0.

fplot(@(x) x.^2/2,[-1,1])
grid on
xlabel('x')
ylabel('E_p(x)')
title('Energía potencial')

La característica esencial de una oscilación libre es que la amplitud A se mantiene constante y por tanto, la energía total se mantiene constante.

E= 1 2 m v 2 + 1 2 k x 2 = 1 2 m v 2 + 1 2 m ω 0 2 x 2 = 1 2 m( v 0 2 + x 0 2 ω 0 2 )= 1 2 m ω 0 2 A 2

Actividades

Se introduce

Se pulsa en el botón Nuevo.


Fuerza proporcional al desplazamiento

En el apartado anterior, hemos estudiado el caso en el que la fuerza F es proporcional al desplazamiento x y de sentido contario a éste, F=-kx. En este apartado, vamos a considerar el caso hipotético en el que la fuerza es proporcional al desplazamiento, F=kx

La ecuación diferencial del movimiento es

d 2 x d t 2 − ω 0 2 x=0,   ω 0 2 = k m

La solución de esta ecuación diferencial es

x=A·exp( ω 0 t)+B·exp(− ω 0 t) v= dx dt = ω 0 ( A·exp( ω 0 t)−B·exp( ω 0 t) )

Las constantes A y B se determinan a partir de la posición inicial x0 y la velocidad inicial v0 en el instante t=0.

t=0 { x 0 =A+B v 0 = ω 0 ( A−B ) A= 1 2 ( x 0 + v 0 ω 0 ), B= 1 2 ( x 0 − v 0 ω 0 ) x= 1 2 ( x 0 + v 0 ω 0 )exp( ω 0 t )+ 1 2 ( x 0 − v 0 ω 0 )exp( − ω 0 t ) x= x 0 cosh( ω 0 t )+ v 0 ω 0 sinh( ω 0 t )

Eliminando el signo menos de la fuerza F, cambia el comportamiento de la solución de la ecuación diferencial, de oscilatorio a exponencial

Vamos a examinar distintos casos:

Referencias

Orion Ciftja. Reversing the sign of Hooke’s law. Phys. Educ. 61 (2026) 055005