Dos osciladores acoplados bajo la acción de una fuerza oscilante

Para estudiar un sistema formado por dos osciladores acoplados, vamos a tomar como modelo el sistema formado por dos partículas iguales de masa m situadas en los extremos de dos muelles de idéntica constante elástica k. El acoplamiento se efectúa uniendo las dos partículas mediante un muelle de constante kc, tal como se puede ver en la figura.

Oscilaciones libres

Llamemos x1 y x2 a los desplazamientos de cada una de las partículas a partir de su posición de equilibrio, medidos como positivos cuando están a la derecha. El muelle de la izquierda se ha estirado x1, el de la derecha se ha comprimido x2 y el central se ha deformado x2-x1. Las fuerzas sobre cada una de las partículas se indican en la parte inferior de la figura.

El muelle central ejerce fuerzas iguales y de sentido contrario sobre cada una de las partículas.

Aplicamos la segunda ley de Newton a cada una de las partículas, y escribimos las ecuaciones del movimiento en forma de ecuaciones diferenciales de segundo orden

m d 2 x 1 d t 2 =k x 1 + k c ( x 2 x 1 ) m d 2 x 2 d t 2 =k x 2 k c ( x 2 x 1 )

Sumando y restando las dos ecuaciones diferenciales tenemos, las ecuaciones diferenciales de un Movimiento Armónico Simple.

d 2 ( x 1 + x 2 ) d t 2 + k m ( x 1 + x 2 )=0 d 2 ( x 1 x 2 ) d t 2 + k+2 k c m ( x 1 x 2 )=0

que representan dos movimientos armónicos simples de frecuencias

ω a 2 = k m ω b 2 = k+2 k c m

Las soluciones de estas dos ecuaciones, son respectivamente

Ψa=x1+x2=Ψ0a sin(ωat+φa)

Ψb=x1-x2=Ψ0b sin(ωbt+φb)

Donde las amplitudes Ψ0a y Ψ0b y las fases iniciales φa y φb están determinadas por las condiciones iniciales: posición inicial y velocidad inicial de cada una de las partículas.

Despejando x1 y x2 de las dos ecuaciones anteriores tenemos

x 1 =( Ψ 0a sin( ω a t+ φ a )+ Ψ 0b sin( ω b t+ φ b ) )/2 x 2 =( Ψ 0a sin( ω a t+ φ a ) Ψ 0b sin( ω b t+ φ b ) )/2

El movimiento general de dos osciladores acoplados puede considerarse como la superposición de dos modos normales de oscilación de frecuencias angulares ωa y ωb.

Condiciones iniciales

En el instante t=0, las posiciones iniciales de las partículas son respectivamente x01 y x02. Las velocidades iniciales son cero.

Las ecuaciones se transforman después de algunas operaciones en

x 1 = x 01 + x 02 2 cos( ω a t)+ x 01 x 02 2 cos( ω b t) x 2 = x 01 + x 02 2 cos( ω a t) x 01 x 02 2 cos( ω b t)

Modos normales de vibración

El primer modo normal de vibración de frecuencia ωa se obtiene cuando los dos osciladores se mueven en fase x01 es igual a x02. El muelle central no sufre ninguna deformación y por tanto, no ejerce ninguna fuerza sobre las partículas, las cuales se mueven como si no estuvieran acopladas.

x 1 = x 01 cos( ω a t) x 2 = x 01 cos( ω a t)

El segundo modo normal de frecuencia ωb se obtiene cuando los dos osciladores se mueven en oposición de fase x01 =- x02. Las ecuaciones del movimiento de cada oscilador se reducen a las siguientes.

x 1 = x 01 cos( ω b t) x 2 = x 01 cos( ω b t)

Experiencia en el aula

Una experiencia con osciladores acoplados que se realiza en el aula suele sorprender a los estudiantes. Consiste en una cuerda que se sujeta por sus extremos situados a la misma altura. Se atan dos péndulos iguales, a dos puntos simétricos de la cuerda, tal como se indica en la figura. Se desplaza uno de los péndulos, por ejemplo el de color rojo, de su posición de equilibrio y se suelta.

El péndulo empieza a oscilar pero su amplitud disminuye con el tiempo, el otro péndulo de color azul que estaba inicialmente en reposo, empieza a oscilar con una amplitud que aumenta.

Al cabo de un cierto tiempo, el péndulo rojo se para momentáneamente, y el péndulo azul oscila con la máxima amplitud. Luego, se cambian los papeles, el péndulo azul disminuye su amplitud con el tiempo, y el péndulo rojo va aumentando su amplitud.

Se analiza la situación desde el punto de vista energético, cómo la energía fluye de un péndulo al otro a través del acoplamiento. Si el acoplamiento es débil, como es éste el caso, la suma total de las energías de los dos péndulos debe ser constante.

En la experiencia en el aula, se desplaza x01 la primera partícula y la segunda se se deja en la posición de equilibrio x02 es cero. Las ecuaciones del movimiento de las partículas se pueden escribir de forma más simple usando las relaciones trigonométricas cosA+cosB y cosA-cosB.

x 1 = x 01 cos ω a ω b 2 tcos ω a + ω b 2 t x 2 = x 01 sin ω a ω b 2 tsin ω a + ω b 2 t

Cuando la amplitud de un oscilador varía con el tiempo, se denomina amplitud modulada. La amplitud del primer oscilador x01cos(ωab)/2 es una función coseno que está adelantada π/2 respecto de la amplitud modulada del segundo oscilador, que es una función seno. Debido a la diferencia de fase entre las dos amplitudes modulantes hay un intercambio de energía entre los dos osciladores. Durante un cuarto de periodo modulante, la amplitud de un oscilador disminuye y la del otro aumenta, dando lugar a una transferencia de energía del primero al segundo. Durante el siguiente cuarto de periodo, la situación se invierte y la energía fluye en dirección opuesta. El proceso se repite continuamente.

Tomemos un sistema formado por dos muelles de constante k=10, el muelle que une las dos partículas tiene una constante kc=0.5. La masa de cada una de las partículas es m=1. Las condiciones iniciales se han fijado de modo que se desplaza la primera partícula de la posición de equilibrio y se suelta, x01=1, v01=0, x02=0 y v02=0. Las frecuencias angulares de los dos modos de vibración ωa y ωb, son cercanas, ya que la constante del muelle que une las dos partículas kc es pequeño frente a k.

k=10; kc=0.5; %constantes de los muelles
m=1; %masa
%frecuencias de los modos de vibración
wa=sqrt(k/m);
wb=sqrt((k+2*kc)/m);
x01=1; %se desplaza la primera partícula

x1=@(t) x01*cos((wa-wb)*t/2).*cos((wa+wb)*t/2);
x2=@(t) x01*sin((wa-wb)*t/2).*sin((wa+wb)*t/2);
hold on
fplot(x1,[0,40]);
fplot(x2,[0,40]);
hold off
xlabel('t')
ylabel('x_1,x_2')
title('Osciladores acoplados')

Oscilaciones forzadas

Las fuerzas adicionales que actúan sobre la primera partícula son:

Las fuerza adicional que actúa sobre la segunda partícula es:

En la parte inferior de figura, se muestran las fuerzas que actúan sobre cada una de las partículas. Las ecuaciones del movimiento son:

m d 2 x 1 d t 2 =k x 1 + k c ( x 2 x 1 )λ d x 1 dt + F 0 cos( ω f t ) m d 2 x 2 d t 2 =k x 2 k c ( x 2 x 1 )λ d x 2 dt

Sumanos y restamos miembro a miembro las dos ecuaciones

d 2 ( x 1 + x 2 ) d t 2 + λ m d( x 1 + x 2 ) dt + k m ( x 1 + x 2 )= F 0 m cos( ω f t ) d 2 ( x 1 x 2 ) d t 2 + λ m d( x 1 x 2 ) dt + k+2 k c m ( x 1 x 2 )= F 0 m cos( ω f t )

Llamamos Ψa=x1+x2, Ψb=x1-x2, y 2γ=λ/m

d 2 Ψ a d t 2 +2γ d Ψ a dt + ω a 2 Ψ a = F 0 m cos( ω f t ) d 2 Ψ b d t 2 +2γ d Ψ b dt + ω b 2 Ψ b = F 0 m cos( ω f t )

Estas son las ecuaciones de dos osciladores forzados de frecuencias propias ωa y ωb, cuya solución hemos estudiado en dicha página. La solución completa de cada una de las ecuaciones diferenciales, se compone de dos términos, el término transitorio que depende de las condiciones iniciales y el estacionario que es el que permanece después de un tiempo teóricamente infinito

Ψ a =( C a cos( ω 1 t )+ D a sin( ω 1 t ) )exp( γt )+ A a cos( ω f t )+ B a sin( ω f t ) A a = F 0 m ω a 2 ω f 2 ( ω a 2 ω f 2 ) 2 +4 γ 2 ω f 2 B a = F 0 m 2γ ω f ( ω a 2 ω f 2 ) 2 +4 γ 2 ω f 2 ω 1 2 = ω a 2 γ 2 Ψ b =( C b cos( ω 2 t )+ D b sin( ω 2 t ) )exp( γt )+ A b cos( ω f t )+ B b sin( ω f t ) A b = F 0 m ω b 2 ω f 2 ( ω b 2 ω f 2 ) 2 +4 γ 2 ω f 2 B b = F 0 m 2γ ω f ( ω b 2 ω f 2 ) 2 +4 γ 2 ω f 2 ω 2 2 = ω b 2 γ 2

El desplazamiento x1 y x2 de cada una de las dos partículas es

x 1 = Ψ a + Ψ b 2 x 2 = Ψ a Ψ b 2 x 1 = 1 2 ( C a cos( ω 1 t )+ D a sin( ω 1 t ) + C b cos( ω 2 t )+ D b sin( ω 2 t ) )exp( γt )+ A a + A b 2 cos( ω f t )+ B a + B b 2 sin( ω f t ) x 2 = 1 2 ( C a cos( ω 1 t )+ D a sin( ω 1 t ) C b cos( ω 2 t ) D b sin( ω 2 t ) )exp( γt )+ A a A b 2 cos( ω f t )+ B a B b 2 sin( ω f t )

Derivando con respecto al tiempo t, obtenemos la velocidad de cada una de las partículas

d x 1 dt = 1 2 ( ( ω 1 C a sin( ω 1 t )+ ω 1 D a cos( ω 1 t ) ω 2 C b sin( ω 2 t )+ ω 2 D b cos( ω 2 t ) )γ( C a cos( ω 1 t )+ D a sin( ω 1 t ) + C b cos( ω 2 t )+ D b sin( ω 2 t ) ) )exp( γt ) ω f ( A a + A b 2 sin( ω f t ) B a + B b 2 cos( ω f t ) ) d x 2 dt = 1 2 ( ( ω 1 C a sin( ω 1 t )+ ω 1 D a cos( ω 1 t ) + ω 2 C b sin( ω 2 t ) ω 2 D b cos( ω 2 t ) )γ( C a cos( ω 1 t )+ D a sin( ω 1 t ) C b cos( ω 2 t ) D b sin( ω 2 t ) ) )exp( γt ) ω f ( A a A b 2 sin( ω f t ) B a B b 2 cos( ω f t ) )

Condiciones iniciales

En el instante t=0,los desplazamientos iniciales de las partículas son x01 y x02, sus velocidades iniciales son v01 y v02, respectivamente

x 01 = 1 2 ( C a + C b + A a + A b ) v 01 = 1 2 ( ω 1 D a + ω 2 D b γ C a γ C b + ω f B a + ω f B b ) x 02 = 1 2 ( C a C b + A a A b ) v 02 = 1 2 ( ω 1 D a ω 2 D b γ C a +γ C b + ω f B a ω f B b )

Despejamos los coeficientes Ca, Cb, Da y Db, que determinan el estado transitorio de cada una de las partículas

C a = x 01 + x 02 A a C b = x 01 x 02 A b ω 1 D a = v 01 + v 02 +γ C a ω f B a ω 2 D b = v 01 v 02 +γ C b ω f B b

Tomemos un sistema formado por dos muelles iguales cuya constante k=1, el muelle que une las dos partículas tiene una constante grande kc=10, para que las dos frecuencias angulares de los modos de vibración ωa y ωb estén alejadas una de otra. La masa de cada una de las partículas es m=1. El factor de rozamiento γ=0.15, la amplitud de la fuerza oscilante F0=1. Vamos cambiando la frecuencia de la fuerza oscilante ωf. Las condiciones iniciales se han fijado de modo que se desplaza la primera partícula de la posción de equilibrio y se suelta, x01=1, v01=0, x02=0 y v02=0

k=1; kc=10; %constantes de los muelles
m=1; %masa
wa=sqrt(k/m); %primer modo  
wb=sqrt((k+2*kc)/m; %segundo modo

g=0.15; %coeficiente de rozamiento
F0=1; %amplitud de la fuerza oscilante F0/m
wf=sqrt(3); %frecuencia angular de la fuerza oscilante 
%condiciones iniciales
x01=1; x02=0; %posiciones iniciales
v01=0; v02=0; %velocidades iniciales

%frecuencias 
w1=sqrt(wa^2-g^2);
w2=sqrt(wb^2-g^2);

%amplitudes del estado estacionario
Aa=F0*(wa^2-wf^2)/((wa^2-wf^2)^2+4*g^2*wf^2);
Ab=F0*(wb^2-wf^2)/((wb^2-wf^2)^2+4*g^2*wf^2);
Ba=F0*(2*g*wf)/((wa^2-wf^2)^2+4*g^2*wf^2);
Bb=F0*(2*g*wf)/((wb^2-wf^2)^2+4*g^2*wf^2);

%coeficientes del estado transitorio
Ca=(x01+x02)-Aa;
Cb=(x01-x02)-Ab;
Da=((v01+v02)+g*Ca-wf*Ba)/w1;
Db=((v01-v02)+g*Cb-wf*Bb)/w2;

x1=@(t) (Ca*cos(w1*t)+Da*sin(w1*t)+Cb*cos(w2*t)+Db*sin(w2*t)).*exp(-g*t)/2+
(Aa+Ab)*cos(wf*t)/2+(Ba+Bb)*sin(wf*t)/2;
x2=@(t) (Ca*cos(w1*t)+Da*sin(w1*t)-Cb*cos(w2*t)-Db*sin(w2*t)).*exp(-g*t)/2+
(Aa-Ab)*cos(wf*t)/2+(Ba-Bb)*sin(wf*t)/2;

%representación gráfica
hold on
fplot(x1,[0,40]);
fplot(x2,[0,40]);
hold off
grid on
legend('x_1','x_2')
title('Dos osciladores acoplados, forzados')
ylabel('x_1,x_2')
xlabel('t')

Observamos el estado transitorio y su evolución hacia el estado estacionario, donde las amplitudes de las oscilaciones son constantes

Amplitudes y fases en el estado estacionario

Después de un tiempo teóricamente infinito, el estado transitorio desaparece, se ha establecido el estado estacionario independiente de las condiciones iniciales. La posición de cada una de las partículas viene dada por

x 1 = A a + A b 2 cos( ω f t )+ B a + B b 2 sin( ω f t ) x 2 = A a A b 2 cos( ω f t )+ B a B b 2 sin( ω f t )

Para calcular la amplitud A1 y la fase φ1 de la oscilación de la primera partícula y la amplitud A2 y la fase φ2 de la oscilación de la segunda partícula, resolvemos los sistema de dos ecuaciones

x 1 = A 1 sin( ω f t+ φ 1 )=( A 1 sin φ 1 )cos( ω f t )+( A 1 cos φ 1 )sin( ω f t ) { A 1 sin φ 1 = A a + A b 2 A 1 cos φ 1 = B a + B b 2 x 2 = A 2 sin( ω f t+ φ 2 )=( A 2 sin φ 2 )cos( ω f t )+( A 2 cos φ 2 )sin( ω f t ) { A 2 sin φ 2 = A a A b 2 A 2 cos φ 2 = B a B b 2

Despejamos las aplitudes A1 y A2

A 1 = 1 2 F 0 m ω a 4 + ω b 4 +4 ω f 4 4 ω a 2 ω f 2 4 ω b 2 ω f 2 +2 ω a 2 ω b 2 +16 γ 2 ω f 2 ( ( ω a 2 ω f 2 ) 2 +4 γ 2 ω f 2 )( ( ω b 2 ω f 2 ) 2 +4 γ 2 ω f 2 ) A 2 = 1 2 F 0 m ( ω a 2 ω b 2 ) 2 ( ( ω a 2 ω f 2 ) 2 +4 γ 2 ω f 2 )( ( ω b 2 ω f 2 ) 2 +4 γ 2 ω f 2 )

k=1; kc=10; %constantes de los muelles
m=1; %masa
wa=sqrt(k/m); %primer modo  
wb=sqrt((k+2*kc)/m); %segundo modo
g=0.15; %coeficiente de rozamiento
F0=1; %amplitud de la fuerza oscilante F0/m

%amplitudes
wf=linspace(0,7,100); %frecuencia angular de la fuerza oscilante 
den=((wa^2-wf.^2).^2+4*g^2*wf.^2).*((wb^2-wf.^2).^2+4*g^2*wf.^2);
A1=sqrt((wa^4+wb^4+4*wf.^4-4*wa^2*wf.^2-4*wb^2*wf.^2+2*wa^2*wb^2+
16*g^2*wf.^2)./den)*F0/2;
A2=sqrt((wb^2-wa^2)^2./den)*F0/2;

%representación gráfica
hold on
plot(wf,A1);
plot(wf,A2);
line([wa,wa],[0,1.6], 'LineStyle','--');
line([wb,wb],[0,1.6], 'LineStyle','--');
title('Dos osciladores acoplados, forzados')
legend('A_1','A_2')
ylabel('A_1,A_2')
xlabel('w_f')
grid on
hold off

Las amplitudes A1 y A2 presentan un máximo para las frecuencias angulares próximas a ωa y ωb de los dos modos de vibración, representadas mediante líneas verticales a trazos

Despejamos las fases, φ1 y φ2

tan φ 1 = A a + A b B a + B b tan φ 2 = A a A b B a B b

k=1; kc=10; %constantes de los muelles
m=1; %masa
wa=sqrt(k/m); %primer modo  
wb=sqrt((k+2*kc)/m); %segundo modo
g=0.15; %coeficiente de rozamiento
F0=1; %amplitud de la fuerza oscilante F0/m

wf=linspace(0,7,200); %frecuencia angular de la fuerza oscilante 

%coeficientes del estado estacionario
Aa=F0*(wa^2-wf.^2)./((wa^2-wf.^2).^2+4*g^2*wf.^2);
Ab=F0*(wb^2-wf.^2)./((wb^2-wf.^2).^2+4*g^2*wf.^2);
Ba=F0*(2*g*wf)./((wa^2-wf.^2).^2+4*g^2*wf.^2);
Bb=F0*(2*g*wf)./((wb^2-wf.^2).^2+4*g^2*wf.^2);

%fases
fi1=atan2((Aa+Ab),(Ba+Bb));
fi2=atan2((Aa-Ab),(Ba-Bb));

%representación gráfica
hold on
plot(wf,fi1);
plot(wf,fi2);
line([wa,wa],[-pi,pi], 'LineStyle','--');
line([wb,wb],[-pi,pi], 'LineStyle','--');
title('Dos osciladores acoplados, forzados')
ylabel('\phi_1,\phi_2')
xlabel('w_f')
set(gca,'YTick',-pi:pi/4:pi)
set(gca,'YTickLabel',{'-\pi','-3\pi/4','-\pi/2','-\pi/4','0',
'\pi/4','\pi/2','3\pi/4','\pi'})
legend('\phi_1','\phi_2', 'location','northwest')
grid on
hold off

Las fases φ1 y φ2 presentan un comportamiento más complejo. Las líneas verticales a trazos, señalan las frecuencias de los dos modos de vibración ωa y ωb

Cuando no hay rozamiento γ=0

Poniendo γ=0, los coeficientes valen

A a = F 0 m 1 ω a 2 ω f 2 B a =0 ω 1 = ω a A b = F 0 m 1 ω b 2 ω f 2 B b =0 ω 2 = ω b C a = x 01 + x 02 A a C b = x 01 x 02 A b ω a D a = v 01 + v 02 ω b D b = v 01 v 02

El desplazamiento x1 y x2 de las partículas en función del tiempo t es

x 1 = 1 2 ( C a cos( ω a t )+ D a sin( ω a t ) + C b cos( ω b t )+ D b sin( ω b t ) )+ A a + A b 2 cos( ω f t ) x 2 = 1 2 ( C a cos( ω a t )+ D a sin( ω a t ) C b cos( ω b t ) D b sin( ω b t ) )+ A a A b 2 cos( ω f t )

Parecería que cuando la frecuencia de la fuerza oscilante ωf se hace igual a las frecuencias de los modos de vibración ωa o ωb, los coeficientes Aa o Ab tendieran a infinito. Sin embrago, el cociente

lim ω f ω 0 ( cos( ω f t)cos( ω 0 t) ω 0 2 ω f 2 )= tsin( ω 0 t) 2 ω 0

donde el subíndice '0' representa 'a' o 'b', tal como hemos visto al estudiar el caso de las oscilaciones forzadas con rozamiento nulo γ=0

>> syms w0 wf t;
>> limit((cos(wf*t)-cos(w0*t))/(w0^2-wf^2),wf,w0)
ans =(t*sin(t*w0))/(2*w0)

Supongamos que desplazamos la primera partícula x01 y la soltamos: x02=0, v01=0, v02=0. Los desplazamientos x1 y x2 de las partículas en función del tiempo t son

x 1 = 1 2 ( ( x 01 A a )cos( ω a t )+( x 01 A b )cos( ω b t ) )+ A a + A b 2 cos( ω f t ) x 2 = 1 2 ( ( x 01 A a )cos( ω a t )( x 01 A b )cos( ω b t ) )+ A a A b 2 cos( ω f t ) x 1 = 1 2 x 01 ( cos( ω a t )+cos( ω b t ) )+ F 0 m ( cos( ω f t )cos( ω a t ) ω a 2 ω f 2 + cos( ω f t )cos( ω b t ) ω b 2 ω f 2 ) x 2 = 1 2 x 01 ( cos( ω a t )cos( ω b t ) )+ F 0 m ( cos( ω f t )cos( ω a t ) ω a 2 ω f 2 cos( ω f t )cos( ω b t ) ω b 2 ω f 2 )

Oscilaciones amortiguadas

Las oscilaciones amortiguadas se producen cuando no se aplica una fuerza oscilante a la primera partícula, F0=0, ωf=0. No hay solución particular que describe el estado estacionario, los coeficientes Aa, Ab, Ba, Bb, son nulos. Los desplazamientos x1 y x2 de las partículas en función del tiempo t son

x 1 = 1 2 ( C a cos( ω 1 t )+ D a sin( ω 1 t )+ C b cos( ω 2 t )+ D b sin( ω 2 t ) )exp( γt ) x 2 = 1 2 ( C a cos( ω 1 t )+ D a sin( ω 1 t ) C b cos( ω 2 t ) D b sin( ω 2 t ) )exp( γt )

Donde los coeficientes Ca, Cb, Da, Db, vienen determinados por las condiciones iniciales. Supongamos que desplazamos la primera partícula x01 y la soltamos: x02=0, v01=0, v02=0. Los desplazamientos x1 y x2 de las partículas en función del tiempo t son

x 1 = 1 2 x 01 ( cos( ω 1 t )+cos( ω 2 t )+γ( sin( ω 1 t ) ω 1 + sin( ω 2 t ) ω 2 ) )exp( γt ) x 2 = 1 2 x 01 ( cos( ω 1 t )cos( ω 2 t )+γ( sin( ω 1 t ) ω 1 sin( ω 2 t ) ω 2 ) )exp( γt )

Tomemos de nuevo el sistema estudiado en las oscilaciones libres, formado por dos muelles de constante k=10, el muelle que une las dos partículas tiene una constante kc=0.5. La masa de cada una de las partículas es m=1. Las condiciones iniciales se han fijado de modo que se desplaza la primera partícula de la posción de equilibrio y se suelta, x01=1, v01=0, x02=0 y v02=0. Las frecuencias angulares de los dos modos de vibración ωa y ωb, son cercanas, ya que la constante del muelle que une las dos partículas kc es peueño frente a k.

k=10; kc=0.5; %constantes de los muelles
m=1; %masa
wa=sqrt(k/m); %frecuencias angulares de los modos de vibración
wb=sqrt((k+2*kc)/m);
g=0.15; %coeficiente de rozamiento

%frecuencias 
w1=sqrt(wa^2-g^2);
w2=sqrt(wb^2-g^2);
x01=1; %se desplaza la primera partícula

x1=@(t) x01*(cos(w1*t)+cos(w2*t)+g*(sin(w1*t)/w1+sin(w2*t)/w2))*exp(-g*t)/2;
x2=@(t) x01*(cos(w1*t)-cos(w2*t)+g*(sin(w1*t)/w1-sin(w2*t)/w2))*exp(-g*t)/2;
hold on
fplot(x1,[0,30]);
fplot(x2,[0,30]);
hold off
xlabel('t')
ylabel('x_1,x_2')
legend('x_1','x_2')
grid on
title('Osciladores acoplados')

Actividades

Se introduce

Oscilaciones forzadas

Nota: cuando se establece la constante de amortiguamiento γ=0 y se introduce una frecuencia ωf de la fuerza oscilante igual a una de las frecuencias de los modos de vibración ωa o ωb, el programa genera un error.

Oscilaciones libres o amortiguadas

La flecha de color negro debajo de la primera partícula (de color rojo) indica la fuerza oscilante F0cos(ωft) que se aplica a esta partícula

Las flechas de color rojo y de color azul, debajo de las partículas inican la velocidad v1 de la primera partícula y v2 de la segunda. Se observará el sentido de las flechas, fuerza y de las velocidades de las partículas cuando la frecuencia ωf de la fuerza oscilante sea igual a las frecuencias de los modos de vibración ωa o ωb,

Se representa, el desplazamiento x1 de la primera partícula en función del tiempo t, en color rojo y de la segunda x2 en color azul. Los datos del tiempo t, x1 y x2, se proporcionan en la parte superior izquierda




Referencias

M. Ebrahim Foulaadvand, Davoud Masoumi. Mechanical filtering in forced-oscillation of two coupled pendulums. Department of Physics. Shahid Beheshti University (2009)