Oscilaciones Acopladas

Acoplamiento de dos osciladores
Valerin Reyes Juan Enrique

Ivan Llana Carhuamaca
MARCO TEORICO
Cuando se conectan dos (o ms) osciladores de tal forma que la energa pueda
pasar de uno a otro nos encontramos frente a una situacin que corresponde al
caso de oscilaciones acopladas. Pero siempre es posible describir el movimiento
de un sistema oscilante en funcin de las coordenadas normales.
Acoplamiento de dos osciladores armnicos
Dos osciladores armnicos iguales de constantes
constante
k1
unidos por un muelle de
k . El movimiento est limitado a la recta que une las masas.
Acoplamiento fuerte.- Ocurre cuando la constante del muelle de

acoplamiento es grande en comparacin con la de los muelles de los
osciladores.
Acoplamiento dbil.- Ocurre cuando la constante del muelle de
acoplamiento es pequea en comparacin con la de los muelles de los
osciladores.
DESARROLLO
Se tiene un sistema de dos masas m=1
1 como se muestra en la figura.
k =c
k 1=1
EJERCICIO 1.-
y tres resortes con constantes
k 1=1
1, c y
Desarrolle analticamente (halle
x ( t ) , y (t ) ) el sistema anterior con las
siguientes condiciones iniciales

v x ( t=0 )=v y ( t=0 ) =0
(1)
x ( t=0 )=1, y ( t=0 )=1 , y
(2)
x ( t=0 )=1, y ( t=0 )=1 , y
(3)
x ( t=0 )=1, y ( t=0 )=0 , y
v x ( t=0 )=v y ( t=0 ) =0

v x ( t=0 )=v y ( t=0 ) =0
SOLUCION
De la figura usando la segunda ley de Newton tenemos el siguiente sistema de
ecuaciones diferenciales
m x =k 1 x +k ( yx ) (1)
m y =k ( yx )k 1 y (2)
Reemplazando los valores de la grfica ( m=1 ,
ecuaciones
(1) ,
(2)
k 1=1
k =c ) en las
y reagrupando tenemos:
x + (1+ c ) xcy (3)

y + ( 1+c ) y cx (4)
Luego usando el operador De I reescribimos
( D2 + ( 1+c ) I ) xcy =0
( D2 + ( 1+c ) I ) ycx=0
Aplicando a la primera expresin
(1/c) ( D2 + ( 1+c ) I )
y sumando a la segunda
tenemos
2
(1/c) ( D2 + ( 1+c ) I ) xcx =0
( D2 + ( 1+c ) I )
[ 2c 2 I ] x=0
[( D 2+ I ) (D 2 +(1+ 2 cI ))]x=0
Suponiendo una solucin oscilante
x ( t )=e it
tenemos entonces
(2 +1 ) (2 + ( 1+2 c ) )=0
De donde obtenemos los valores de
1 =1
2 = 1+2 c
y cuya solucin es
x ( t )=a 1 cos ( t ) +a2 sin ( t ) + a3 cos ( t 1+ 2 c ) + a4 sin ( t 1+2 c )

x ( t )=a1 sin ( t ) +a 2 cos (t )a3 1+2 c sin ( t 1+2 c ) +a 4 1+2 c cos ( t 1+2 c )
Y como
x
[+ (1+ c ) x ]/c
y =
entonces
y (t )=a1 cos ( t ) +a 2 sin ( t )a 3 cos ( t 1+2 c )a4 sin ( t 1+2 c )

y (t )=a 1 sen ( t ) +a 2 cos +a3 1+ 2c sen ( t 1+2 c )a4 1+2 c cos ( t 1+2 c )
Usando las condiciones iniciales tenemos
(1)
v x ( t=0 )=v y ( t=0 ) =0
x ( t =0 )=1, y ( t =0 )=1 , y
x ( 0 )=a1 +a 3=1
y ( 0 )=a 1a3 =1
Entonces
a1=1
a3 =0
x ( 0 )=a2 +a 4=0
y ( 0 )=a 2a 4=0
Entonces
a2=a4 =0
Luego las ecuaciones para
(2)
x(t )
y (t)
son
x ( t )=cos (t )
(1a)
y (t )=cos ( t )
(1b)
x ( t=0 )=1, y ( t=0 )=1 , y
v x ( t=0 )=v y ( t=0 ) =0
x ( 0 )=a1 +a 3=1
y ( 0 )=a 1a3 =1
Entonces
a1=0
a3 =1
x ( 0 )=a2 +a 4=0
y ( 0 )=a 2a 4=0
Entonces
a2=a4 =0
(3)
x(t )
y (t)
son
x ( t )=cos ( t 1+ 2c )
(2a)
y (t )=cos ( t 1+2 c )
(2b)
x ( t =0 )=1, y ( t =0 )=0 , y
v x ( t=0 )=v y ( t=0 ) =0
x ( 0 )=a1 +a 3=1
y ( 0 )=a 1a3 =0
Entonces
a1=1 /2
a3 =1/2
x ( 0 )=a2 +a 4=0
y ( 0 )=a 2a 4=0
Entonces
a2=a4 =0
x(t )
y (t)
son
x ( t )=[ cos (t ) +cos ( t 1+2 c ) ] / 2
( t )cos ( t 1+2 c )
cos /2
y ( t )=
(3a)
(3b)
EJERCICIO 2.El siguiente cdigo en Mathematica resuelve la ecuacin del movimiento del
sistema:
Para las condiciones (1) tenemos

sol 1=Flatten [DSolve [{x '' [t ]+(1+ c) x[ t ]cy [t ]= 0, y ''[ t]+(1+c ) y [t ]cx [t ]= 0, x [0]=1, x ' [0]=
1
1
{x [t] t (1+2 t ), y [t] t (1+2 t )}
2
2
Transformando las exponenciales complejas a senos y cosenos obtenemos
1
x (t)= ( t + t )=cos ( t)
2
1
y (t)= ( t + t ) =cos ( t)
2
Con esto comprobamos las ecuaciones (1a) y (1b).
sol 2=Flatten [DSolve [{x '' [t ]+(1+ c) x [t ]cy [t]= 0, y ''[ t]+(1+c ) y [t ]cx [t ]= 0, x [0]=1, x ' [0]=
1
{x [t] 12c t ( 1+ 2 12 c t ),
2
1
y [ t] 12 c t (1+ 2 12 c t )}
2
Transformando las exponenciales complejas dado que
12 c=i 1+2 c a
senos y cosenos obtenemos

1
x (t)= ( i 1 +2 ct + i 1+2 c t )=cos ( t 1+ 2c )
2
y (t)=
1 i 1+ 2c t i 1+2 ct
(
+
)=cos (t 1+2 c)
2

sol 3=Flatten[ DSolve[{x ''[t ]+(1+ c) x [t ]cy [t]= 0, y ''[ t ]+(1+c ) y [t ]cx [t ]= 0, x [0]= 1, x '[0]=
{x [t]
1 t 12 c t t 12 ct 2 t +12 ct t + 2 12 ct
( +
+
+
),
4
1
t
y [t] t 12 c t ( 12 ct 2 t +12 c t + t +2 12 c t )}
4
Transformando las exponenciales complejas dado que
12 c=i 1+2 c a
senos y cosenos obtenemos
1
( i 1+ 2 c t+ i 12c t + it + t )= [cos ( t ) +cos ( t 1+2 c)]
2
1
x( t)=
4
1
( i 1+ 2 c t+ i 12c t it t )= [cos ( t )cos ( t 1+2 c )]
2
1
y (t )=
4
Pregunta De los tres casos anteriores Cules corresponden a modos normales
de oscilacin puros?
Los casos (1) y (2) corresponden a modos normales de oscilacin puros pues solo
dependen a la vez de una sola frecuencia normal de oscilacin, mientras que el
caso (3) no corresponde a un modos normal de oscilacin puro pues depende
ambas frecuencias normales.
EJERCICIO 3.El siguiente cdigo grafica el resultado sol3 para
c=1 .
Plot [{x [t ]/. sol 3/. c 1, y [t ]/. sol 3/. c 1}, {t , 0,30 }]
Grafica del resultado sol3 para
c=0.1
Plot [{x [t ]/. sol 3/. c 0.1, y [t ]/. sol 3/. c 0.1 }, {t , 0,30}]
Grfica (2)
Grafica del resultado sol3 para
c=5
Plot [{x [t ]/. sol 3/. c 5, y [t ]/. sol 3/. c 5 }, {t , 0,30}]
Grfica (3)
Pregunta Qu caractersticas y diferencias nota en la grficas para diferentes

valores de c ?
Podemos observar en la grfica (2) que para un acoplamiento dbil ( c=0.1 )
x
y
las variables
e
varan ms lentamente, lo cual reflejara un
movimiento lento de las masas, mientras que para un acoplamiento fuerte (
c=5 ) vemos en la grfica (3) un cambio mucho ms rpido en las variables
x e y esto originaria un movimiento relativamente rpido en comparacin
con el acoplamiento dbil.
Animacin Haga la siguiente animacin
1) Modos de oscilacin puros
1a) Usando la sol1
t =[x [t ]]/. sol 1
xx
t =[ y [t ]]/. sol 1
yy
t , c 0 :=
plotosc
{
Graphics
, Disk [{N [2.5+ yy [t ]/. c c 0 ],0 }, 0.1]}, Axes True ,
PlotRange {{1.1,5 .1}, {0.14,0 .14 }}
Animate [ plotosc [t ,0.1], {t , 0,100}, AnimationRate 1]
En la animacin podemos observar que las masas se mueven en la misma
direccin de vaivn y con la misma velocidad lo cual confirma los resultados
obtenidos en las ecuaciones (1a) y (1b) las cuales sealan que las masas deben
moverse en sentidos iguales.
caso1_modospuros.nb
1b) Usando la sol2
t =[x [t ]]/. sol 1

xx
t =[ y [t ]]/ . sol 1
yy
t , c 0 :=
plotosc
{
Graphics
, Disk [{N [3.75+ yy [ t ]/. c c 0 ],0 },0.1]}, Axes True ,
PlotRange {{1.1,5 .1}, {0.14,0 .14 }}
En la animacin podemos observar que las masas se mueven en direccin
opuesta y con la misma velocidad lo cual confirma los resultados obtenidos en
las ecuaciones (2a) y (2b) las cuales sealan que las masas deben moverse en
sentidos contrarios.
caso2_modospuros.nb
2) Para la solucin 3 se mostrara la animacin para un acoplamiento dbil y
uno fuerte usando el cdigo
2a) Acoplamiento dbil
t =[x [t ]]/ . sol 3
xx
t =[ y [t ]]/. sol 3
yy
t , c 0 :=
plotosc
{
Graphics
, Disk [{N [ 4 + yy [t ]/. c c 0], 0 },0.1]}, Axes True ,
PlotRange {{0.1,5 .1 }, {0.14,0 .14 }}

Con esto comprobamos que el acoplamiento dbil hace que las masas oscilen
lentamente en el tiempo, esto se puede apreciar claramente en la siguiente
animacin.
Acoplamiento_Debil.nb
2b) Acoplamiento fuerte
t =[x [t ]]/. sol 3
xx
t =[ y [t ]]/ . sol 3
yy
t , c 0 :=
plotosc
{
Graphics
, Disk [{N [ 4 + yy [t ]/. c c 0], 0 },0.1]}, Axes True ,
PlotRange {{0.1,5 .1 }, {0.14,0 .14 }}
Animate [ plotosc [t ,5], {t , 0,100}, AnimationRate 1]
Con esto comprobamos que el acoplamiento fuerte hace que las masas oscilen
rpidamente en el tiempo, esto se puede apreciar claramente en la siguiente
animacin.
Acoplamiento_Fuerte.nb
Otro mtodo de solucin:
a) Calculando las coordenadas normales de oscilacin.
Primero calculamos las matrices T y V , para ellos usamos las expresiones
de la energa cintica y potencial del sistema.
1
1
T = x 2 + y 2
2
2
1
1
1
U= x 2 + c ( yx )2+ y 2
2
2
2
Y como los elementos de las matrices
estn relacionados de la
siguiente forma
2
T=
1
T
2 i=1 j=1 (i , j) i j
2
1
U= V (i , j) i j
2 i=1 j=1
1=x
Donde
2 = y
, comparando las expresiones obtenemos los valores
que corresponden a los elementos de las matrices.
( )
T= 1 0
0 1
V = 1+c c
c 1+c
b) Usando la siguiente expresin para calcular las frecuencias normales de

oscilacin tenemos
det |T2j V |=0
Luego
det
1+ c 2
c
2 =0
c
1+ c
2 2
Entonces
(1+c ) c =0
2
1 =1
Luego
2 =1+ 2 c
( 1+2 c 2) ( 12 ) =0
entonces las frecuencias normales de oscilacin
son
1 =1
2 = 1+2 c
c) Procederemos a calcular los autovectores correspondientes a cada
autovalor.
1+c2j
c
c
1+c 2j
Primero consideraremos la condicin que los vectores estn normalizados,

entonces tenemos
T T =1
Considerando
(s)=
a1 s
a2 s
( )
y usando la matriz
T obtenemos la condicin de
normalizacin:
a21 s + a22 s=1
Para
21 =1
(cc cc )( aa )=0
11
21
Entonces
a11 =a21
Luego normalizando obtenemos que
2
a11 = =a21
2
(1) = 2 1
2 1
()
Para
22 =1+ 2 c
c a =0
(c
c c )( a )
12
22
Entonces
a22=a12
Luego normalizando obtenemos que
2
a12= =a22
2
(2) = 2 1
2 1
( )
A
es
( )
d) Luego la matriz
A=( )= 2
2
2
(1 )
(2)
2 = 2 1 1
2 1 1
2
2
Luego las coordenadas normales

T
C=A
T
)
son
Donde
()
= x
y
y reemplazando los valores obtenemos
C 1 = 2 1 1 1 0 x
C=
2 1 1 0 1 y
C2
( ) (
)( )( )
De donde obtenemos las coordenadas normales

2
C1 = (x+ y )
2
2
C2 = ( x y )
2
Despejando obtenemos las coordenadas
en funcin de las
coordenadas normales
2
x= (C 1+C 2)
2
entonces
2
x = ( C 1+ C 2)
2
2
y= ( C 1C2 )
2
entonces
2
y = ( C 1C2)
2
Reemplazamos en
x + (1+ c ) xcy (3)
y + ( 1+c ) y cx (4)
Tenemos en (3)
2 ( C + C ) + ( 1+ c ) 2 C +C c 2 C C =0
( 1 2)
(
1
2
2)
2
2
2 1
Simplificamos
C 1 + C 2+ ( 1+ 2c ) C2 +C 1=0
Tenemos en (4)
2 ( C C ) + ( 1+c ) 2 C C c 2 C =0
( 1 2)
( 1 )
1
2
2
2
2
Simplificamos
C 1 C 2( 1+2 c ) C2 +C 1=0
Cuyo sistema es
C 1 + C 2+ ( 1+ 2c ) C2 +C 1=0
C 1 C 2( 1+2 c ) C2 +C 1=0
Y cuya solucin
C1 ( t )=a1 cos ( t )+ a2 sin ( t )
C2 ( t )=b1 cos ( 1+2 c t ) +b2 sin ( 1+2 c t)
CONCLUSION.Las oscilaciones acopladas mantienen su movimiento gracias a la transferencia

de energa de un resorte a otro, adems estos movimientos desordenados pueden
ser descritos gracias a las frecuencias normales de oscilacin.
Adems la solucin usando coordenadas normales muestran que estn solo
dependen a la vez de una sola frecuencia normal de oscilacin.
El acoplamiento fuerte genera una gran energa cintica en comparacin al
acoplamiento dbil.

Oscilaciones Acopladas

Cargado por

Copyright:

Formatos disponibles

Oscilaciones Acopladas

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Oscilaciones Acopladas

Cargado por

Copyright:

Formatos disponibles

Acoplamiento de dos osciladores

Valerin Reyes Juan Enrique

unidos por un muelle de

k . El movimiento est limitado a la recta que une las masas.

Acoplamiento fuerte.- Ocurre cuando la constante del muelle de

y tres resortes con constantes

Desarrolle analticamente (halle

x ( t ) , y (t ) ) el sistema anterior con las

siguientes condiciones iniciales

x ( t=0 )=1, y ( t=0 )=1 , y

x ( t=0 )=1, y ( t=0 )=1 , y

x ( t=0 )=1, y ( t=0 )=0 , y

v x ( t=0 )=v y ( t=0 ) =0

x + (1+ c ) xcy (3)

Luego usando el operador De I reescribimos

(1/c) ( D2 + ( 1+c ) I ) xcx =0

x ( t )=a 1 cos ( t ) +a2 sin ( t ) + a3 cos ( t 1+ 2 c ) + a4 sin ( t 1+2 c )

y (t )=a1 cos ( t ) +a 2 sin ( t )a 3 cos ( t 1+2 c )a4 sin ( t 1+2 c )

v x ( t=0 )=v y ( t=0 ) =0

Luego las ecuaciones para

x ( t=0 )=1, y ( t=0 )=1 , y

v x ( t=0 )=v y ( t=0 ) =0

Luego las ecuaciones para

v x ( t=0 )=v y ( t=0 ) =0

Luego las ecuaciones para

x ( t )=[ cos (t ) +cos ( t 1+2 c ) ] / 2

Para las condiciones (1) tenemos

senos y cosenos obtenemos

Con esto comprobamos las ecuaciones (2a) y (2b).

Transformando las exponenciales complejas dado que

senos y cosenos obtenemos

EJERCICIO 3.El siguiente cdigo grafica el resultado sol3 para

Grafica del resultado sol3 para

Grafica del resultado sol3 para

Plot [{x [t ]/. sol 3/. c 5, y [t ]/. sol 3/. c 5 }, {t , 0,30}]

Pregunta Qu caractersticas y diferencias nota en la grficas para diferentes

t =[x [t ]]/. sol 1

Animate [ plotosc [t ,0.1], {t , 0,100}, AnimationRate 1]

Y como los elementos de las matrices

, comparando las expresiones obtenemos los valores

que corresponden a los elementos de las matrices.

b) Usando la siguiente expresin para calcular las frecuencias normales de

entonces las frecuencias normales de oscilacin

Primero consideraremos la condicin que los vectores estn normalizados,

Luego normalizando obtenemos que

Luego normalizando obtenemos que

Luego las coordenadas normales

y reemplazando los valores obtenemos

De donde obtenemos las coordenadas normales

CONCLUSION.Las oscilaciones acopladas mantienen su movimiento gracias a la transferencia

También podría gustarte