Lévano Huamaccto If 2017

UNIVERSIDAD NACIONAL DEL CALLAO
FACULTAD DE CIENCIAS NATURALES Y MATEMATICA
UNIDAD DE INVESTIGACION DE LA FACULTAD DE CIENCIAS
NATURALES Y MATEMATICA
11:: U 030£l'lcVEE':I£S 030I£AnNACl0NALDELCALLA0;

\§gCi0}g42\_\ ? n} mm$mW :
SQ B 1 2013 f
{IE 6 030; 1 030)
:2 E 034»\§ )2
'l\> 030D:A\
55:5 (03
0 //$50
EM; 2018}
INFORME FINAL DEL PROYECTO DE INVESTIGACION
034CAOS
EN LA DINAMICA DE UN OSCILADOR NO LINEAL
AMORTIGUADO Y FORZADO 035
AUTOR: Lic. CARLOS ALBERTO LEVANO HUAMACCTO
(Periodo de Ejec}401ciénz
De 01 de enero del 2017 al 12 de Diciembre del 2017)
(Resolucién de aprobacién N°110-2017R)
' Callao, 2017

INDICE
IJNDICE 1
1.1.|ndice de tablas 3
1.2.Indice de }401guras 4
1.3.1ndice de gra}401cas 5
lI.RESUMEN Y ABSTRACT
2.1. Resumen 6
2.2. Abstract ' 7
III. INTRODUCCION
3.1. Exposicién del problema de investigacién 8
3.2. Importancia y justi}401cacién

de la investigacién 1 10
IV. MARCO TEORICO
4.1. Antecedentes o estudios brevios 12
4.2. Bases teéricas
4.2.1.La Teoria de| Caos 12
4.2.2. Indicadores de| Caos 14
a) Seccién de Poincaré 15
b) Espectros de Potencia 15
c) Exponentes de Lyapunov 16
- 4.3. Ecuacion de Duf}401ng 19
4.4. Marco conceptual 28
a) Caos 28
b) Mecénica no lineal y Caos 28
c)Linealidad y no linealidad 29
V. MATERIALES Y METODOS
5.1. Materiales 30
5.1.1. De ejecucién 30
5.1.2. De impresién 30
5.2. La poblacion de la investigacién y la muestra 31
030n\
5.3. Las Técnicas de anélisis o Ios medios estadisticos 31
5.4.Las técnicas de anélisis o Ios métodos estadisticos aplicados 31
5.5. Metodologiaa 31
VI. RESULTADOS 34
VII. DISCUSION 41
VIII. REFERENCIALES 42
IX. APENDICES 43
X. ANEXOS 51
Matriz de Consistencia
2
TABLA DE CONTENIDO
Valores para los diferentes parémetros de|

Tabla N°6.1 _ _ _ Pag.34
Sistema no Lmeal Amomguado y Forzado
3
TABLA DE FIGURAS
Figura 4.1 Super}401cie

de la Seccién de Poincare Pag.15
Figura 4.2 Espectro de potencia
Figura 4.3 Oscilador Forzado lineal y espacio fase
Figura 4.4 Oscilador Forzado No lineal y espacio fase Pag.2D
Figura 4.5 Sistema No Lineal Fisico Pag.22
Figura 4.6 Sistema No Lineal Fisico moviendo
Figura 4.7 Sistema No Lineal Fisico desplazado
4 .
TABLA DE GRAFICAS
Grafica 6.1 Espacio de fase sin Caos frecuencia 0.01 Pag.35
Gra}401ca
6.2 Espacio de fase sin Caos frecuencia 0.05 Pag.36
Gra}401ca
6.4 Espacio de fase con Caos frecuencia 0.1 Pag.38
Gra}401ca
6.5 Espacio de fase oon Caos frecuencia 0.3 Pag.39
6.6
Gra}401ca Espacio de fase oon Caos frecuencia 0.6 Pag.40
5
RESUMEN
En este trabajo de investigacién determinamos el estarid de Caos de un sistema
fisico Oscilador no lineal amortiguado y forzado con una fuerza externa de tipo
senoidal a través del indicador Ilamado exponentes de Lyapunov; para caracterizar
la dinémica de| Oscilador no lineal determinas primero el espacio de fase donde
L se observa si presenta comportamiento desordena para una frecuencia relativa
' determinada, luego calculamos su oorrespondiente exponente de Lyapunov a
' través de un programa desarrollado por Wolf et al. Se obtuvo valores del
exponentes de Lyapunov positivos para dinémicas desordenadas del Oscilador no
Lineal Amortiguado y Forzado y, valores de| exponente de Lyapunov negativos
para dinémicas ordenadas del Oscilador corroborando Ia hipétesis de este
trabajo de investigacion.
W
' 6
ABSTRACT
We determine the state of Chaos of a non-linear oscillator damped and forced
physical system with an external force of type sinusoidal through the indicator
called exponents of Lyapunov; characterize the dynamics of the non-linear
oscillator. }401rst
determine the phase space where it is observed if it presents
disordered behavior for a given relative frequency, then we calculate its
corresponding exponent of Lyapunov through a program developed by Wolf et al.
Positive values of the Lyapunov exponents were obtained for disordered dynamics
of the Negligible and Foroed Non-linear Oscillator and, negative Lyapunov
exponent values for ordered dynamics of the Oscillator corroborating the
hypothesis of this research work.
-% 7
Ill. INTRODUCCION
3.1. La exposicién del problema de investigacién
Con el avance de la Ciencia se sabe que son mas sistemas no lineaies que los
sistemas Iineales que se puede encontrar en la naturaleza ,y cuando se
estudian los primeros, se analiza su dinamica en aproximacién lineal, 0 sea las
variables contenidas en la ecuacién diferencial 0 de movimiento son de primer
orden, en esa aproximacion la ecuacién de movimiento del sistema mencionado
Iogra cumplir con el principio de superposicién, importante principio de las
ecuaciones Iineales; sin embargo, Ia aproximacion lineal Ileva a una descripcién
de la dinémica de los sistemas mencionados restringida, como en ei caso del
sistema oscilatorio periédico Ilamado Pénduio Fisico, que se ie estudia
menos de cinco grados a pesar que este

generalmente para éngulos peque}401os,
muestra movimiento en la practica para anguios mayores(Tay|or, 2013).
En |as Oltimas décadas se ha asociado el comportamiento desordenado que
presenta Ia dinémica de los sistemas no Iineales con el problema Ilamado Caos,
caracterizado por una fuerte sensibilidad a las condiciones iniciales de las
ecuaciones no Iineales; se sabe de los sistemas no Iineales que su dinémica
puede mostrar comportamiento ordenado como comportamiento desordenado
y, las condiciones para que se dé el segundo caso es materia de estudio en
este trabajo de investigacién.
Los Sistemas Oscilatorios Iineales periédicos han sido ampliamente estudiados
experimental y teéricamente en la fisica y sirven de modelos para estudiar las
interaociones entre sistemas, en las }401ltimasdécadas |os Sistemas
amortiguados y forzados son usados como modelos con frecuencia ; sin
embargo, el interés por conocer mas sobre los Sistemas no Lineales
Amortiguado y Forzado , y en especial saber las condiciones cuando muestran
031 comportamiento desordenado o Caos es lo que se busca(Douglas G., 2006).
\
031 Es importante determinar cuéndo un sistema presenta comportamiento de Caos
031 porque permite conocer esos comportamientos indeseables de los sistemas y,
de los sistemas flsicos que siwen de modelo.
El problema identi}401cado
es la detenninacién de| estado de Coas en Ios
030 Sistemas no Lineales que muestran comportamiento desordenado. Se plantea
la determlnacién del estado del Caos usando los Exponentes de Lyapunov,
propuesto por el ruso matemético del mismo nombre que propone que el
cambio en la dinémica se puede medir a través de la divergencia de las
variables dinémicas, que es cuanti}401cada

por unos parémetros llamados
Exponentes de Lyapunov (WolfA.,1984).
La formulacién y planteamiento del Problema son:
a) Problema General
Problema general que se propuso en este 'trabajo dé investigacién es el
siguiente,
g,Es posible determinar el Caos en la dinémica de un Oscilador no Lineal
Amortiguado y Forzado con los Exponentes de Lyapunov?
b) Problemas Especificos:
' - (,Cuél es la técnica numérica apropiada para calcular los Exponentes
de Lyapunov en la determinacién de| estado dinémico del Caos de un
Oscilador no Lineal Amortiguado y Forzado?
- ¢;Es posible determinar un modelo fisico para estudiar el Caos en la
dinémica de un Oscilador no lineal Amortiguado y Forzado?
N 9
c) Objetivo General y Alcances de la investigacion
-Objetivo General
Los objetivos planteados en el presente trabajo fueron:
o Desarrollar la técnica numérica apropiada para calcular los
Exponentes de Lyapunov y determinar el estado de| Caos en la dinémlca
de un Oscilador no Lineal Amortiguado y Forzado.
o Desarrollar y calcular los Exponentes de Lyapunov de un Oscilador
no lineal Amortiguado y Forzado en Resonancia.
-Objetivos Especificos
Exponentes de Lyapunov y determinar el estado del Caos en la dinémica
o Desarrollar y calcular Ios "Exponentes de Lyapunov de un Oscilador
3.2. lmportancia y la Justificacion de la lnvestigacién
Sabemos que el estudio de los Sistemas Fisicos no Lineales se ha
_ lncrementado y el problema de Caos se observa en varios sistemas, no
solamente encontramos Sistemas no Iineales en Fisica sino también en otras
areas como Quimica, Biologia, Economla, Medicina, lngenieria y por ello se
dice que el Caos en ubicuo. Los trabajos de modelamiento que realizamos, |os
fisicos, son interdisciplinarios y por ello surge la necesidad de estudiar y
comprender |as herramientas para afrontar el fenémeno del Caos. En ese
sentido es importante que los modelos no Iineales, en especial los mecénicos,
_ 030 que proponemos para estudiar otros fenémenos de la naturaleza debe ser
bien estudiado, saber cuando pueden funcionar en estado de Caos
X 10
031 controlando Ios parémetros propios del sistema para saber |as condiciones
que Ilevan al Caos.
o La investigacién que deseamos realizar corresponde a un valor te6rico 024
024 préctico, desde el punto de vista de una Investigacién Bésica, y constituiré un
' aporte importante dentro de la }402sica

computacional, Ia cual, es una de las
Iineas de especializacién de Fisica.
030 030 11
IV. MARCO TEORICO
_ 4.1. Antecedentes o Estudios Previos
' Zhang e Harne (2016) propone el estudio del Oscilador con efecto
débilmente no lineal, Ilamado Oscilador de Duf}401n,

muestra que este
Oscilador presenta comportamiento de Caos para determinados parémetros.
Cabana e Laroze (2015) muestra el anélisis de las ecuaciones de Landau-
Lifshitz o Landau-Lifshitz-Gilbert, que describe Ia dinémica de una particula
en un campo magnético oscilante, se analiza la dinémica hacienda variar
parémetros del Sistema y se estudia el Caos en su espacio de fase.
Mouronte M (2004) , muestra él estudia el Caos del sistema, a través del
espacio de fase donde |as bifurcaciones que son desdoblamiento de la
dinémica de ,es una forma cualitativa de describir el Caos en la dinémica de
un sistema. 030
Kapitaniak T. (1994), muestra que estudio del comportamiento de los
sistemas cuasi periédicos forzados a través de los Atractores, que son
regiones en el espacio de fase o oonfiguracién hacia donde tiende la
dinémica, es una forma de estudio de los sistemas cuasi periédicos.
La revisién bibliogréfica realizada muestra que en los trabajos no usan los
V exponentes de Lyapunov especi}401camente

para el estudio de| estado el Caos
de un sistema de un Oscilador Amortiguado y Forzado. El antecedente
descrito ha motivado a desarrollar Ia presente investigacién como una
contribucién dirigida a dotar de los instrumentos necesarios de consulta,
actualizacién y orientacién bésica para aquellos que tienen interés de esta
especializacién dentro de la Fisica no Lineal.
12
i
4.2. Bases Teoricas
4.2.1. La Teoria de| Caos
Usualmente en }401sica
aprendimos que las leyes son descritas por ecuaciones
de movimiento (ecuaciones diferenciales), como por ejemplo la Iey de
Newton entre otros, y la solucién analitica o numérica de estas ecuaciones
se hacen teniendo |as condiciones iniciales y de contorno; podemos conocer
como es su dinémica en el tiempo, esto es el punto de vista determinista,
entonces, sistemas fisicos conservativos o disipativos son deterministas
porque elios obedecen ecuaciones diferenciales deterministas (Taylor J.,
2013).
La teoria de Caos, es una teoria que estudia el comportamiento desordenado
de los sistemas no Iineales , como |os fisicos, a los que se conoce su
031 ecuacién dinémica resulto de la aplicacién de una iey o principio, en ei caso
de la Mecanica por ejemplo las Leyes de Newton; por ello se dice también
Caos determinista, |os sistemas no Iineales en régimen pertubativos
usualmente exhiben comportamiento complicado Ilamado Caos(José J. e
Saletan E.,1998).
Si |os estados iniciales de sistemas deterministas son exactamente
conocidos, futuros estados dinamicos son predicidos teéricamente; esta
visién determinista sobrevivio para ser cuestionada después de las famosas
obras visionarias de Henri Poincaré sobre la mecénica celeste realizadas a
}401nales
del siglo XIX que eran contran'o al determinismo(Ra}401ada
A.,1990).
Poicaré proponia que supiéramos exactamente la Iey de la naturaleza y la
situacion del universo en el momento inicial, se podria predecir exactamente
Ia situacion de ese mismo universo en un momento sucesivo, pero incluso si
fuera ei caso de que las leyes naturales ya no tuvieran ningiin secreto para
. nosotros, pero no es asi en los resultados numéricos que nos describen que
& 13
solo .podn'amos conocer Ia situacién inicial aproximadamente. Si nos
permitiera predecir la situacion siguiente con la misma aproximacién, eso es
todo lo que necesitamos, de esta manera, las mediciones efectuadas sobre
el estado de un sistema en un momento dado pueden no permitirnos predecir
la situacion futura incluso moderadamente lejos, a pesar de que las
ecuaciones de gobierno son exactamente conocidas (Douglas R., 2006).
No ocurrio hasta |os a}401os

70, después de casi un siglo de investigaciones,
aprendieron que los sistemas caéticos solo pueden resolverse
numéricamente, y no hay formas sencillas y generales de predecir cua'ndo un
sistema exhibiré caos. También hemos aprendido que el caos determinista
esté siempre asociado con sistemas no Iineales, Ia no Iinealidad es una
condicion necesaria para el caos, pero no una su}401ciente(José

J. e Saletan
E.,1998).
Solo en los }401ltimos

a}401os,
con el advenimiento de las computadoras que se
permitio estudiar el caos ya que ahora es posible realizar célculos numéricos
de la evolucién temporal de las propiedades de los sistemas sensibles a las
condiciones iniciales, comenzamos a comprender la existencia de| caos
cuando |as computadoras estaban fécilmente disponibles para calcular Ias
historias de largo tiempo neoesarias para explicar el comportamiento
discutido, el modelo deterministico siempre produciré la misma salida de una
condicion inicial o estado inicial determinado(Doug|as R., 2006).
4.2.2. indicadores de| Caos
La teoria de| caos, es una teoria que trata sobre el fenomeno del
comportamiento desordenado que los sistemas fisicos no Iineales presentan
, en especial, muestran una fuene sensibilidad a las condiciones iniciales; Ia
no Iinealidad de las ecuaciones (deterministas) hacen que no se puedan
V resolver analiticamente porque no cumplen con el principio de superposicion,
es por ello que en esta teoria se han desarrollado indicadores de| Caos para
\ 14
determinar si el sistema muestra el caos , entre los indicadores tenemos :
_ Espectros de Potencias y Super}401cies

de Seocién de Poincaré y Exponentes
' de Lyapunov(Doug|as R., 2006).
a) Seccién de Poincaré
Este indicador es bésicamente el mapa que se muestra al interoeptar las
1 curvas de las trayectorias, tomadas periédicamente. este procedimiento es
més cualitativo para averiguar si el sistemaVpresenta Caos(Douglas R.,
2006).
FIGURA N° 4.1
A SUPERFICIE DE LA SECCION DE POINCARE
(3) H 030?/" 034"| 0343 030v
3&1 ~
-}401?:§
mo» 2
53$
. _ _ _ $9; r
:<-on
Fuente: Ra}401ada.
Dinémica CIésica.1989
b) Espectros de Potencias
Aqui lo que se hace es determinar el médulo de la transformada de Fourier
en funcién de las frecuencias. una sistema oscilatorio lineal muestra un
Respectro de potencia de un solo pico, de un sistema oscilatorio cuasi

» 15
periodico dosomés pioo,pero siempre discretos latransformada de Fourier; I
sin embargo, cuando un sistema }401sico

esté en estado de caos debe mostrar
una transformada de Fourier continua, para ello se tiene que hacer variar
todos los parémetros y variables del sistema y no siempre se puede
determinar el limite a partir de| cual el sistema presenta caos (Ra}401ada
A.,1990)
FIGURA N°4.2
I ESPECTRO DE POTENCIA DE UN SISTEMA PERIODICO, CUASI
PERIODICO Y EN CAOS
030 ; I. q -_ .; ;:',-;2;": 030_ 031r::-_;'- .
o . .. - .- ,. 3 o _v ». '
Fuente: Ra}401ada,
Dinémica 030CIésica,
1989
c) Exponentes de Lyapunov
Los términos no Iineales de las ecuaciones diferenciales ampli}401can
diferencias en las condiciones iniciales, de esta

exponencialmente peque}401as
manera, Ias leyes de evolucion deterministas crean comportamientos
caéticos, incluso en ausencia de nuevas }402uctuaciones

internas, en el
régimen caotico no es posible predecir exactamente Ia evolucién del estado
del sistema durante un tiempo arbitrariamente, esta es la caracteristica
impredecible del caos.
\ b 16
La evolucién temporal se rige por un espectro continuo de frecuencias
responsables de un comportamiento aperiédico, para este movimiento se
presentan patrones estacionarios, es decir. patrones que se repiten de forma
no periédica. Lyapunov creé un método conocido como exponente
caracteristico de Lyapunov para cuanti}401car

la sensible dependencia de las
condiciones iniciales de| comportamiento caético,A la estabilidad de los
sistemas dinémicos; con este método es posible determinar |os requisitos
minimos de las ecuaciones diferenciales que son necesarias para crear el
caos. Estudiemos el caso de sistemas con una sola variable que asumen dos
estados iniciales x., y x., + s, di}401riendo

por una peque}401a
cantidad £(Ra}401ada
A.,1990)
Queremos investigar los posibles valores de x,, después de n iteraciones de
los dos valores iniciales, Ia diferencia d,. entre los dos valores xn después de
n iteraciones (omitiendo por simplicidad el subindice) se da
aproximadamente.
d,, = f(x',, + 5) 024-

f(x,,) = sen) 030 (1)
Donde, A, es el exponente de Lyapunov que representa el coe}401ciente

del
crecimiento exponencial promedio por unidad de tiempo entre los dos
estados. De lo que vemos. si A es negativo,'|as dos orbitas eventualmente
convergerén, pero si es positivo |as trayectorias préximas divergen. La
diferencia entre los dos estados iniciales se escribe como
«:1 = 034X0
+ s) ~ r(x.,) 2 5(: 024:)KO (2)
Ahora, para evitar Ia confusién que a veoes se encuentra en la literatura
caética, recordamos que,
\ 17
Xn+1 = f(Xn)=-f(f(Xn-1)): f(f(f(Xn-2)))= = f(f(f(-»-(f(Xo)))))- (3)
Donde el superindice n indica Ia n-ésima iteracién de| mapa. Después de
un gran numero de iteraciones. Ia diferencia entre los estados cercanos
es. 4
I d = f(x,1 + 2)
034 f(x,,) = f(x0 + 2)
024 f}401xo)
024 = sen" (4)
Dividiendo Ia ecuacién por s y tomando el Iogaritmo a ambos lados,
obtenemos I
In{[f n (X0 + 8) _
f S034
031 030_O
=
In[e 21}.] =
nix
-
(5)
_ Mxo) é ilnldf (6)
Donde,
«u M, {T}{-- 024m(,,
}-..{ 024m
-nio_ 024dF1(X;1_1) dFl{Xn_7) dfutxo)
& 18
Reempiazando la ecuacién 7 en la ecuacién 6 para calcular el exponente
de Lyapunov, que puede tener valores positive 0 negativo, el valor positivo
que el sistema presenta comportamiento caético.

signi}401ca
I = lim (5)1 N dfil 8

(X0) s 024nan H" dxo ()
4.3. Ecuacién de Duf}401ng
' Un caso ilustrativo disipativo es el movimiento de una particula con masa m
sometida a un potencial y a una fuerza disipativa es proporcional a, pig
respectivamente; se considera débilmente no lineal cuando, 6<<1.
dzx dz 3
024,+B 024 024x+8x
=FCos(mt) (9)
dt dt
Donde (3. 5 , on son parémetros.
19
FiGURA N° 4.3
OSCILADOR AMORTIGUADO FORZADO LINEAL Y SU ESPACIO DE FASE
(3) (b)
6.0 6.0
030 = * if u 031
U
""lt:uJ ; 5 030 ' _ / 03403 - » - 030vtcm

»-
030(Swan
031 4035.0 u)6u.o
031 ""94.a 0.0 4,0
I x
Fuente: Mouronte M., Caracterizacién y Anélisis de Sistemas Dinémicos no
Lineales Mediante el Estudio del Mapa de frecuencias y Espacio de fases, 2004
FIGURA N 4.4
034 '
OSCILADOR AMORTIGUADO FORZADO NO LINEAL EN CAOS Y SU
ESPACIO DE FASE
(a) (b) '
is .- J ' I xi}
'»<o.q _ ; I _. 024v;o.f;._ 151.: .

031 I ' § E _V 031 . _ 030_ 030/j-.'.~.¢?
030 §_,5_
030;~ 030_x..
034an 1
4 _ I '7 O :.%
l'8l0.0 1055.0 no0.o 24,0 my .31)

t x
Fuente: Mouronte M., Caracterizacién y Anélisis de Sistemas Dinémioos no
Lineales Mediante el Estudio de| Mapa de frecuencias y Espacio de fases,
2004
20
La Dinémica esté compuesta por dos campos de estudio, Ia Dinémica Lineal
y la Dinémica no Lineal, esta Liltima trata del comportamiento desordenado
dinémico de los Sistemas no Lineales, es aqui donde Ios sistemas presentan
el Caos, a mejor dicho el Caos determinista, porque cuentan con ecuaciones
dinémicas con fuerte sensibilidad a las condiciones iniciales. Esto signi}401ca
que una peque}401a

variacién en las condiciones iniciales puede alterar
signi}401cativamente
el comportamiento dinémico del sistema convirtiendo los
valores que va a tomar en ei tiempo (Saletan & José, 1998).
Otros sistemas también se puede Ilegar a la forma de| Oscilador de Duf}401ng,
si consideramos un cuerpo de masa ,M, ai cual le llamaremos la particuia(
sin rotar), suspendido entre dos muelles iguales, tal como se muestra en la
figura 5. Si ambos muelles estén sin de deformar(o sea, sin tensién y, por
tanto; sin energia potencial los dos lados) cuando Ia panicula se encuentre
en su posicién de equilibrio y, si despreciamos Ias fuerzas gravitatorias, cada
uno de los muelles ejerceré una fuerza,
, 024x(s
l)
024
sobre la particula (x es la constante de recuperacién de cada mueile), 1 , es
la Iongitud sin deformar dei resorte. (Saletan & José, 1998).
\ 21
_ FIGURA 4.5
SISTEMA NO LINEAL COMPUESTO DE UNA PARTlC_ULA DE MASA M
" ~
Fuente: Marién J. Dinémica Clésica de las partfculas y Sistemas, 1981
FIGURA 4.6
SISTEMA NO LINEAL COMPUESTO DE UNA PARTICULA DE
MASA M MOVIENDOSE LONIGITUDINALMENTE
K 030 E \\: 030
______""_ .._'f.__\;,:M
K 030er
/ , 031
s
' Fuente: Marion J., Dinémica Clésica de las particulas y Sistemas, 1981
- 22
La fuerza que acttia sobre la particula de masa ,M, es ,F, igual a,
F ' 0242z:(.-:l)$'crz(9),
024
como : = wxz + I3
de forma que,
_ 9 x I
Jen = 024
= 024 024 024 024 024
5 \/x3 + :2
Resulta la fuerza en funcién de la posicién,
F= 1 024.'+
\;A 035
(v 030x2
030
e :2 L) = zzcc 1 1
F (10)
H031031 V1 +G} 031
Si suponemos que, . es peque}401o

y desarrollamos en serie el radical,
tendremos(Marion J.,1981),
12' _,cl( 031£)3

1_"L(§.)2 (11
_ 030
1 4 .1 )
23
Y cuando despreciamos todos los términos, salvo el primero, seré
k
F 5-! (11))?
Si alargamos inicialmente ambos muelles una distancia, d ,para sujetar a
ellos Ia masa , m , en posicién de equilibrio, la fuerza seré (Marion J.,1981),
' _ vcd N 3 3
» F(X)§ 0242 x030
024[rc(l
024l- 024d,/I
]x
De donde podemos determinar Ia energia potencial considerando la
expresién,
F _ dU(x)
_ dx
_ 55 _ xii-d) 3 _ __.s'u-{ac}
A 2(1)x [ F J) 030
024 .2; ' (12)
Obtenemos la energia potencial siguiente,
_E E 030
U(x)-z(2(l'))x i[j""F 034"
2 +4 x.
4 031] (13)
24
Luego considerando que el sistema se mueve solamente en forma horizontal
la energia cinética viene dado por , '1 030=}M:*

y conjuntamente con la
energia potencial , U(x) , podemos expresar el Lagrangiano del sistema
L=5M72'2 0245(2(§))x2 0245F3:£]x 030*

Z Z I
(14, 4- A 030
Cuando ademés estén presente otras fuerzas, disipativas y no disipativas en
Mecénica Clésica se usa la siguiente ecuacién,
d ''at.\ (al.) 1 2
030 " T 024 024
= 030
- + ~ 15
d: (.6471) aa, 0 Q3
031 ( )
Donde:
Q} _ representa a las fuerzas disipativas,
Q} , representa |as fuerzas que no son disipativas ni conservativas.
En nuestro caso usaremos consideraremos que nuestra particula esté
sometida a una fuerza disipativa de la forma, Q5 =

030
024b( 030:y 024
una
030:)
fuerza
externa dependiente del tiempo, Q)-2 to (.'os(ut)

024 ambas en direccion
horizontal.
\ 25
030 ' Reemplazando el Lagragiano en la ecuacién dinémica y las fuerzas
respectivas se obtiene Ia siguiente expresién,
3 1 1 -,2 030
_:(_; .2 _1li((§-t}402j
.4) =_ 030g
¢:(a2(zM? 034 any 2(2 J1 4 :3 1 b(a't)+
' I-'0 Cos(wt)
(16)
030 desarrollando Ia expresién de la manera que se muestra para obtener una
ecuacion dinémica de un Sistema no Lineal Amortiguado y Forzado,
+2 .1 +
030 = 024b¢_ +034§)
F0 Cos(wt)
(17)
A la ecuacién,
3 ii
}402 030
E . L-#3 .a_5 .
1:1 + M (dc) + 2 (in) X + [ 131:: J1 " M COS(wt) (18)
_ 26
Se hacen |as siguientes representaciones,
E ,= 55 _. [ 031 =030 030 030 031' 030')I

M 2/3» 2(m) 024wa ~
2
6-
' nos permite expresar de la siguiente forma Ia ecuacién dinémica,
- i . F
.1: + 23:: + (A)(j.22. +
0316x3 = B 024:_ 031Cas(_wt) (19)
De donde observamos que s << 1

030 es bien peque}401o
por lo tanto podemos
a}401rmar
que la ecuacién dinémica corresponde a un Sistema débilmente no
lineal Amortiguado y Forzado.
27
4.4. Marco conceptual
a) Caos
Es un fenomeno, que provoca Ia imposibilidad de realizar cualquier prediccion
}401able
de su comportamiento a pesar de de tener un carécter determinista, aun
en ausencia de las perturbaciones extemas al sistema, los temas de esta
clase desarrollan una comportamiento caétlca (Alemén R., 2010).
b) Mecénlca no lineal y Caos
La mayoria de los sistemas cuyos ecuaciones de movimiento son no Iineales
pueden exhibir Caos. Este fenémeno que se mani}401esta

en muchos campos
diferentes como los sistemas mecénicos oscilanles, reacciones quimicas, }402ujo
de }402uidos,
rayos léser, crecimiento de poblaclones, transmislén de
enfermedades , y muchos otros, supone que aunque un sistema obedezca a
ecuaciones deterministas de movlmiento(|eyes de Newton) su comportamiento
detallado futuro puede, en la préctica, ser impredecible(Taylor, J., 2013).
El comportamiento de un sistema mecénico caético puede ser muy
complicado, y la necesidad de describir este comportamiento ha generado una
enorme variedad de maneras nuevas de visualizar el movimiento de los
sistemas como érbitas del espacio de estados, secclones de Poincare,
diagramas de bifurcacién.
La teoria del Caos esté dividida en dos amplias éreas de estudios, sistemas
disipativos, tales como el péndulo amortiguado, y no disipativos, o sistemas
Hamiltonianos (Taylor, J., 2013).
\ 28
' c) Linealidad y no Iineaiidad
La ecuacion diferencial es lineal si relaciona Ia variable 0 variables
dependientes y sus derivadas solo Iinealmente. El francés Poincare observo
que la tierra, la Iuna y el sol. sus ecuaciones exhiben sensibilidad a las
condiciones iniciales. Caracteristica de| movimiento caético.
La no linealidad es esencial para el caos, pero la no linealidad no garantiza el
caos ademés el sistema tiene que ser 034oomplicado

como
035 |os sistemas
disipativos, sus ecuaciones de movimiento cuando se ponen en su forma
esténdar autonoma deben ser no Iineales y tener N variables con N>=3(Tay|or,
J., 2013).
\ 29
5.1. Materiales
5.1.1. De ejecucién
' o El trabajo de investigacién pertenece a la investigacién teérica,
especi}401camente
a la Fisica computacional, por ello se usé ordenadores
Core I-3 para hacer las simulaciones y célculos, programas en fortran
90 para Runge Kutta de 410 orden (Press W., Teukolsky 3., Vettering
W., Flannery B.,1992) para estudiar el comportamiento de la dinémica
de| espacio de fase.
- Internet . para la busqueda de informacién relacionado al trabajo de
investigacién.
o Fotocopias, de documentos que ayudan en el desarrollo de| trabajo de
investigacién.
o Dispositivos de almacenamiento (CD, USB) . para guardar Ias corridas
con el programa.
5412! De la impresibn
Para Ia parte }401nal

del trabajo de investigacién se uso los siguientes
materiales:
o Papel Bond en formato A-4
o Elementos de escritorio
o Tintas para impresién
X 30
5.2. La poblacién y muestra
Como es un trabajo de investigacién de Fisica Teérica no requiere de
poblacién y muestra.
5.3 Las
030 técnicas, procedimientos e instrumentos de recoleccién de datos
No se usan técnicas estadisticas por tratarse de un trabajo de la Fisica
Teérica.
031 5.4. Las Técnicas de anélisis o métodos estadisticos aplicados
No aplica para este trabajo.
5.5. La Metodologia
-En este trabajo se consideré el siguiente sistema, }401gura

4.7,
FIGURA 4.7
SISTEMA N0 LINEAL COMPUESTO DE UNA PARTICULA DE
MASA M y RESORTES ESTIRADO X
K .: \\\s
M V L X \ \«
030
__.____ ________., 030M
I. 030 9 030I;
¥ Fuente: Marién J., Dinémica Clésica de las particulas y Sistemas, 1981
31
024La
recopilacién de datos se hace a través de un programa mencionado
Iineas arriba de la ecuacién de la dinamica del sistema,
11 ax 024.
024f
+ 2;? (- 024)
+ W 030):
+ 6.1.63 = F(.'ns(mt) , (20)
dz 030 dt
y el programa para obtener |os exponentes de Lyapunov fue adaptado de
(\NoIfA., Swift J., Swinney H. and Vastano J., 1985).
Para ello a la ecuacién ultima se convierte en un sistema dinémico de dos
variables, donde Ia segunda ecuacién contiene el término que contiene al
tiempo, I-'Cas(wt),
dz: _ '
dL 3
iii _ _9 _ 2 __ 3 ,
at ..,B at w 3: 6x + I'Co.s(cot)
(21)
V El Sistema. se hizo discreto y se Ileva al programa de Runge Kutta de 4to
orden con intervalo de tiempo 0,015, para analizar si se observa
comportamiento desordenado o Caos en el espacio de fase para con los
parametros de la tabla N° 6.1.Para esos valores se calculan |os exponentes
de Lyapunov méximo positivo, que mostré Caos.
K 32 .
l
i
Luego haciendo un cambio de variable:
Del tiempo, "L = w,,L, ~
de la sustitucién de la frecuencia, fl 024

3 .La ecuacién adimensional
Wu
resulta,
2 2 .
i+2_(3wo +wD2x+.5:c =iCvs

031 (£1) (22)
:12 -1: M me
dia: 23 das _ 5 _3 Po
- 024;
+ + It + 024zx
= :2-C0S(QT) (23)
03 it ma dr mg mg, M
, . . . 2:3 031 2
-En el analisis se mantuvo constante |as cantidades , :3
024 ,
OJ.) 0
i y solamente se hace variar Ia frecuencia relativa Q .

0
Se determiné Ios espacios de fases para las diferentes frecuencias desde,
m = }402,}402
hasta
030lma
m = 3,4 mu con los datos que se muestra en el cuadro
. . . _ Z}402 s Fa
sugulente, mantenlendo constante las cantidades - 024,024 024 024g-, 030T
.
can an we M
\ 33
VI. RESULTADOS
TABLA N 6.1
034 .
VALORES PARA LOS DIFERENTES PARAMETROS DEL SISTEMA NO
LINEAL AMORTIGUADO Y FORZADO
1 9031 9 = I;
034O2 1:932!" - W0
N0 linealidad Fuerza externa Frecuencia
034E111
Ejij
E
u
j
j
j
j
ET
Fuente: Autoria propia
\ 34
GRAFICA N 6.1
035
ESPACIO DE FASE SIN CAOS DEL OSCILADOR NO LINEAL AMORTIGUADO
Y FORZADO, FRECUENCIA 0,01 Y EXPONENTE LYAPUNOV NEGATIVO -
0.215798E+O1
- 20 _ ""_' 034 034

'sa!ida35,dat' u 23 034
- 024- 024
o
\ -1.. 030 I, f
'5 \$-.\ U 030
\ 034\1 030:,-
-10 V ,
in
now
030 2 4 5 s in 12 14 16 13
Fuente: Autoria propia obtenido con gra}401cador

Gnuplot.
\ 35
GRAFICA N° 62
20
M an 024 024
0305 1" "\\\

030
g /5" ,,
1/ 03 035 'z\ \_r "#7 031 «
5 030f 03
1: I r 030. )
.» 1 « 2 c, 030
030S % . . \ /I
\ \ V. . 030
.10 . V/4,4, *" 0307' M
\ 030/'
,;2a_-2_oM -15 -10 -5 0 5 to 15 20

.324538E+02

Gnuplot.
\ 36
GRAFICA N° 6.3
Y FORZADO, FRECUENCIA 0,08 Y EXPONENTE LYAPUNOV NEGATNO
20 /\
15 H _'__,A _a_:_ Rt 034 -1if\* ;*;-\:\\ 030
,4- 030? 030. 031 *;§\\_\'_\ 030030'7///

030 030
AA\\.~v~'v*~~A
f W
° 1: I)))))! 031) ) 030l)
Va 1* 3"! 031 034031
034'
. ' A \
-212 H
_ _-1s__A_ -1o -5 0 5 1a 15 20

Gnuplot.
X 37
n I// ., h 034N§Cn:gRT|
5 0301 030T;33:::9P»:~§
\;:\ 031_;.~.MW§ 030
1 034 .
Q sra}401h}401q}401r
035 5 %
agbtenidoks 030
2° }402
. "5 um»
Ill 030-I 031
\\V"77>
I K
1 030
\,/
V 030
4% 030:-I
"
X 39
g / %%
05$ 035?
g9F 031:Fof 030 /1
\;
4] 030? 20
A9
%
VII. DISCUSION
El espacio de fase del Sistema no Lineal Oscilatorio Amortiguado y Forzado que
observa de las }401guras

6.1,6.2,6.3 muestran una dinérnica desordenada, las Iineas
se cruzan varias veces , la dinémica se realiza alrededor de dos centros de
oscilacion, que se les llama atractores, |as frecuencias relativas corresponde a
Q=0.1, 0.3 , 0.6 respectivamente.
El comportamiento observado aqui del oscilador no lineal, también Io podemos
observar en el Oscilador no lineal Amortiguado y Forzado de un péndulo donde Ia
forma de la ecuacion dinémica se aproxima al sistema usado en nuestro trabajo de
investigacion.
Los exponentes de Lyapunov determinados para esos tres casos resultaron
teniendo con valores positivos, resultado que esta en la direccién de demostrar de
la hipotesis de este trabajo de investigacion, que los exponentes de Lyapunov son
también indicadores de| estado de Caos del Sistema no Lineal Oscilatorio
Amortiguado y Forzado.
Los resultados obtenidos en este trabajo va en la misma direccion con otros
trabajos de investigacién, donde usan también al exponente de Lyapunov como
indicador, pero para otros sistemas no Iineales (Cabanas A.,Laroze D.,2014).
\ 41
r
SUBROUTINE derivs(x,y.dydx)
REAL x,y(").dydx('),A,C,E,D
A=0.01
c=o.o2 030 .
D=3O.0
E=0.08
dydx(1)=y(2)
! dydx(2)=-((w1"2) 030y(1))
dydX(2)=-A"Y(2)-Y(1)-C*((Y(1))**3)+D 030C0S(E*X)
I dydx(3)=y(4)
I dydx(4)= 024((w2**2)*v(3))
return
END
! rkdumb.for
SUBROUTINE rkdumb(vstart,nvar,x1,x2,nstep,derivs)
INTEGER nstep,nvar,NMAX,NSTPMX
PARAMETER (NMAX=50,NSTPMX=100000)
REAL x1 ,x2,vstart(nvar),xx(NSTPMX),y(NMAX, NSTPMX)
EXTERNAL derivs
COMMON lpathl xx,y
!CU USES rk4
INTEGER i,k
REAL h,x,dv(NMAX),v(NMAX)
do 11 i=1,nvar
v(i)=vstar1(i)
V(i.1)=V(i)
11 continue
xx(1)=x1
44
031 x=x1 030
h=(x2-x1)/nstep
do 13 k=1,nstep
call derivs(x,v,dv)
call rk4(v,dv,nvar,x,h.v.derivs)
if(x+h.eq.x)pause 'stepsizé not signi}401cant

in rkdumb'
x=x+h
xx(k+1)=x
do 12 i=1,nvar
y(i,k+1)=v(i)
12 continue
13 continue
return
END
! rk4.for
SUBROUTINE rk4(y,dydx,n,x,h,yout,derivs)
INTEGER n,NMAX
REAL h,x,dydx(n).y(n),yout(n) '
EXTERNAL derivs
PARAMETER (NMAX=50)
INTEGER i
REAL h6,hh,xh,dym(NMAX).dyt(NMAX),yt(NMAX)
hh=h"0.5 -
V h6=h/6.
xh=x+hh
do 11 i=1,n
Yt(i)=Y(i)+hh*dYdX(i) '
11 continue
cail derivs(xh,yt,dyt) .
do 12 i=1,n -
Vt(i)=y(i)+hh"dyt(i)
45
12 continue
call derivs(xh,yt,dym)
do 13 i=1,n
Yt(i)=y(i)+h*dYm(i)
dym(i)=dyt(i)+dym(i)
13 continue
call derivs(x+h,yt,dyt)
do 14 i=1,n
yout(i)=y(i)th6 031(dydx(i)+dyt(i)+2. 030dym(i))
14 continue
return
END
9.2.PROGRAMA DE EXPONENTES DE LYAPUNOV ADAPTADO DE (WOLF A.
1985)
program ode
nn=n"(n+1)=
nn1=nn+1
parameter (n=3, nn=12, nn1=13)
implicit double precision(a-h,o-z)
external fcn
dimension y(nn1),yprime(nn1),v(nn1),A(nn1),B(nn1),C(nn1)
dimension D(nn1),cum(n),znorm(n),gsc(n)
open(1,}401|e="Lout.txt")
111 format(1x,'nstep, stpsze, irate, io:') .
read(*,*) nstep,stpsze,irate,io
write(1,*) nstep,stpsze,irate.io
write(1,*)
v(1)=10.0
v(2)=1.0
v(3)=1.0
K 46
tme=0.0
do 10 i=n+1,nn
v(i)=0.0
10 continue
do 20 i=1,n
v((n+1)*i)=1.0
cum(i)=0.0
20 continue
do 100 m=1,nstep
if ((m/100)*100.EQ.m) wn'te(*,*)m
do 25 j=1,irate
do 26 i=1,nn
030/(i)=V(i)
26 continue
t=tme
call fcn(t, y, yprime)
do 27 i=1,nn
A(i)=yprime(i)
27 continue
do 28 i=1,nn
y(i)=v(i)+(stpsze*A(i))/2.0
28 continue
t=tme+stpszel2.0
do 29 i=1,nn
B(i)=yprime(i)
29 continue
do 30 i=1,nn
y(i)=v(i)+(stpsze'B(i))/2.0
30 continue
t=tme+stpsze/2.0
X 47
do 31 i=1,nn
C(i)=yprime(i)
31 continue
031 do 32 i=1,nn
y(i)=V(i)+(StPSZe*C(i))
32 continue
t=tme+stpsze
call fcn(t, y, yprime) 030
do 33 i=1,nn
D(i)=yprime(i)
33 continue
do 34 i=1,nn
v(i)=v(i)+stpsze*(A(i)+D(i)+2.0*(B(i)+C(i)))/6.0
34 continue
tme=tme+stpsze
25 continue
znorm(1)=0.0
do 38 j=1,n
znorm(1)=znorm(1)+v(n*j+1) 030*2
38 continue
znorm(1)=sqrt(znorm(1))
do 40 j=1,n
v(n*j+1)=v(n*j+1)/znorm(1)
40 continue
do 80 j=2,n '
do 50 k=1 ,1-1
gsc(k)=0.0
do so |=1,n _
gsc(k)=gsc(k)+v(n*I+j)*v(n*|+k)
50 continue
\ 48
do 60 k=1,n
do 60 l=1,j-1
v(n*k+j)=v(n*k+j)-gsc(|)"v(n*k+|)
60 continue
znorm(j)=0.0
do 70 k=1,n
znorm(j)=znorm(j)+v(n*k+j)**2
70 continue
znorm(j)=sqrt(znorm(j))
do 80 k=1,n
v(n*k+j)=v(n*k+j)/znorm(j)
80 continue
do 90 k=1,n
cum(k)=cum(k)+d|og(znorm(k))/a|og(2.0)
90 continue
if(mod(m,io).nei0) goto 100
write(1,334) tme,(cum(k)/tme,k=1_n)
334 format(1x,f12.6,2x,e12.6,2x,e12.6,2x,e12.6)
_ 100 continue
end
*m*m**..m*....Mmm...m..........m..............*........*
subroutine fcn(t,y,yprime)
b=4.0
sg=16.0
r=45.92
VPrime(1)=S9*(V(2)-Y(1))
yprime(2)=-y(1)*y(3)+r*y(1)-y(2)
VPrime(3)=Y(1)*Y(2)-b*Y(3)
do 10 i=0,2
vprime(4+i)=S9"(V(7+i)-Y(4+i))
vprime(7+i)=(r-y(3))*v(4+i)-v(7+i)-v(1) 030v(10+i)
\ 49
' 030 vr>rime(10+i)=y(2)*y(4+i)+y(1)*y(7+i)-b*v(10+i)
10 continue
. return
_ . end
031 V i X _ 50
/ X. ANEXOS
Matnz de Consistencia
FORMULACI 0 N DEL PROBLEMA OBJETIVOS HIPOTESIS VARIABLES INDICADORES METODOLOGIA
I I Detenninacion del Caos

030 (IE5 posible determinar el Caos I I Los Exponentes de
I I I en la dnnémnca de un I
en la dlnémlca de un Oscilador I I Lyapunov determnnan el
I I Oscilador no Lineal I I
General no Llneal Amortiguado y Forzado I estado dlnémlco del Caos
Amomguado y Forzado I I I I I I
con los Exponentes de en un Oscilador no Lineal Esta lnvestlgaclon es de tipo:
con los Exponentes do I I I
Lyapunov? L AmomguadoyForzado. Célculos aplloada. Con un dlse}401o de
ya unov. _ I I
p numéncos ejecuclon transversal. El
I I lndependiente: iguales o Sistema do anélisls es el

I I Desarrollar la tecnlca I I . _ I
¢;Cuél es la técmca numénca II I La técnica numénca Exponentes de mayores a Osculador no Lineal
I numenca apropnada para I
apropiada para calcular Ios _ basada en la Lyapunov cero de los Amomguado y Forzado; a
' calcular los Exponentes I I I I
Exponentes de Lyapunov en la d L ortonormahzaclén de exponentes este Oscilador no Lineal se
e a unov
determinacion del estado I y P y Gram-Schmidt local de Lyapunov formularé el Lagrangiano de|
I I determinar el estado del I I I I
dlnémlco del Caos de un I I perrmte determinar los sistema y luego la ecuaclén
I _ Caos en la dlnémlca de I I
I Oscllador no LinealAmomguado I _ Exponentes de Lyapunov dlnémloa de Lagrange, con
Especu}401c un OSCIIBUOF no Lnneal I I
y Forzado? ' I an un Oscalador no Llneal ello Implemenlaremos un
a Amortlguado y Forzado. I I
Amomguadoy Forzado, programa basado La técnica
I _ I I numérica basada en 031

la
.. (,Es posible detennlnar un Desarrollary calcularlos I I Dependlentez I
I I I El modelo }401sloo
Osculador I I I orlonormallzaclon de Gram-
modelo flslco para estudiar el Exponentes de Lyapunov I I Estado dunémuco Componamne I
I I I I no lineal Amortiguado y I I Schmidt local.
Caos en la dunémuca de un en un Oscllador no Inneal I I del Oscllador no mo Caolloo
_ I _ Forzado permite estudiar I I
Oscllador no lineal Amomguado Amortiguado y Forzado I Lnneal de la vanable
I la teoria del Caos con los I I I
yFotzado? en Resonancia. Amortiguado y dlnémlca.
Exponentes de Lyapunov.
Forzado.
51
UNIVERSIDAD NACIONAL DEL CALLAO
FACULTAD DE CIENCIAS NATURALES Y MATEMATICA
UNIDAD DE INVESTIGACION DE LA FACULTAD DE CIENCIAS
NATURALES Y MATEMATICA
11:: U 030£l'lcVEE':I£S 030I£AnNACl0NALDELCALLA0;

\§gCi0}g42\_\ ? n} mm$mW :
SQ B 1 2013 f
{IE 6 030; 1 030)
:2 E 034»\§ )2
'l\> 030D:A\
55:5 (03
0 //$50
EM; 2018}
INFORME FINAL DEL PROYECTO DE INVESTIGACION
034CAOS
EN LA DINAMICA DE UN OSCILADOR NO LINEAL
AMORTIGUADO Y FORZADO 035
AUTOR: Lic. CARLOS ALBERTO LEVANO HUAMACCTO
(Periodo de Ejec}401ciénz
De 01 de enero del 2017 al 12 de Diciembre del 2017)
(Resolucién de aprobacién N°110-2017R)
' Callao, 2017

INDICE
IJNDICE 1
1.1.|ndice de tablas 3
1.2.Indice de }401guras 4
1.3.1ndice de gra}401cas 5
lI.RESUMEN Y ABSTRACT
2.1. Resumen 6
2.2. Abstract ' 7
III. INTRODUCCION
3.1. Exposicién del problema de investigacién 8
3.2. Importancia y justi}401cacién

de la investigacién 1 10
IV. MARCO TEORICO
4.1. Antecedentes o estudios brevios 12
4.2. Bases teéricas
4.2.1.La Teoria de| Caos 12
4.2.2. Indicadores de| Caos 14
a) Seccién de Poincaré 15
b) Espectros de Potencia 15
c) Exponentes de Lyapunov 16
- 4.3. Ecuacion de Duf}401ng 19
4.4. Marco conceptual 28
a) Caos 28
b) Mecénica no lineal y Caos 28
c)Linealidad y no linealidad 29
5.1. Materiales 30
5.1.1. De ejecucién 30
5.1.2. De impresién 30
5.2. La poblacion de la investigacién y la muestra 31
030n\
5.3. Las Técnicas de anélisis o Ios medios estadisticos 31
5.4.Las técnicas de anélisis o Ios métodos estadisticos aplicados 31
5.5. Metodologiaa 31
VI. RESULTADOS 34
VII. DISCUSION 41
VIII. REFERENCIALES 42
IX. APENDICES 43
X. ANEXOS 51
Matriz de Consistencia
2
TABLA DE CONTENIDO
Valores para los diferentes parémetros de|

Tabla N°6.1 _ _ _ Pag.34
Sistema no Lmeal Amomguado y Forzado
3
TABLA DE FIGURAS
Figura 4.1 Super}401cie

de la Seccién de Poincare Pag.15
Figura 4.2 Espectro de potencia
Figura 4.3 Oscilador Forzado lineal y espacio fase
Figura 4.4 Oscilador Forzado No lineal y espacio fase Pag.2D
Figura 4.5 Sistema No Lineal Fisico Pag.22
Figura 4.6 Sistema No Lineal Fisico moviendo
Figura 4.7 Sistema No Lineal Fisico desplazado
4 .
TABLA DE GRAFICAS
Gra}401ca
6.2 Espacio de fase sin Caos frecuencia 0.05 Pag.36
Gra}401ca
6.4 Espacio de fase con Caos frecuencia 0.1 Pag.38
Gra}401ca
6.5 Espacio de fase oon Caos frecuencia 0.3 Pag.39
6.6
Gra}401ca Espacio de fase oon Caos frecuencia 0.6 Pag.40
5
RESUMEN
En este trabajo de investigacién determinamos el estarid de Caos de un sistema
fisico Oscilador no lineal amortiguado y forzado con una fuerza externa de tipo
senoidal a través del indicador Ilamado exponentes de Lyapunov; para caracterizar
la dinémica de| Oscilador no lineal determinas primero el espacio de fase donde
L se observa si presenta comportamiento desordena para una frecuencia relativa
' determinada, luego calculamos su oorrespondiente exponente de Lyapunov a
' través de un programa desarrollado por Wolf et al. Se obtuvo valores del
exponentes de Lyapunov positivos para dinémicas desordenadas del Oscilador no
Lineal Amortiguado y Forzado y, valores de| exponente de Lyapunov negativos
para dinémicas ordenadas del Oscilador corroborando Ia hipétesis de este
trabajo de investigacion.
W
' 6
ABSTRACT
We determine the state of Chaos of a non-linear oscillator damped and forced
physical system with an external force of type sinusoidal through the indicator
called exponents of Lyapunov; characterize the dynamics of the non-linear
oscillator. }401rst
determine the phase space where it is observed if it presents
disordered behavior for a given relative frequency, then we calculate its
corresponding exponent of Lyapunov through a program developed by Wolf et al.
Positive values of the Lyapunov exponents were obtained for disordered dynamics
of the Negligible and Foroed Non-linear Oscillator and, negative Lyapunov
exponent values for ordered dynamics of the Oscillator corroborating the
hypothesis of this research work.
-% 7
Ill. INTRODUCCION
3.1. La exposicién del problema de investigacién
Con el avance de la Ciencia se sabe que son mas sistemas no lineaies que los
sistemas Iineales que se puede encontrar en la naturaleza ,y cuando se
estudian los primeros, se analiza su dinamica en aproximacién lineal, 0 sea las
variables contenidas en la ecuacién diferencial 0 de movimiento son de primer
orden, en esa aproximacion la ecuacién de movimiento del sistema mencionado
Iogra cumplir con el principio de superposicién, importante principio de las
ecuaciones Iineales; sin embargo, Ia aproximacion lineal Ileva a una descripcién
de la dinémica de los sistemas mencionados restringida, como en ei caso del
sistema oscilatorio periédico Ilamado Pénduio Fisico, que se ie estudia
menos de cinco grados a pesar que este

generalmente para éngulos peque}401os,
muestra movimiento en la practica para anguios mayores(Tay|or, 2013).
En |as Oltimas décadas se ha asociado el comportamiento desordenado que
presenta Ia dinémica de los sistemas no Iineales con el problema Ilamado Caos,
caracterizado por una fuerte sensibilidad a las condiciones iniciales de las
ecuaciones no Iineales; se sabe de los sistemas no Iineales que su dinémica
puede mostrar comportamiento ordenado como comportamiento desordenado
y, las condiciones para que se dé el segundo caso es materia de estudio en
este trabajo de investigacién.
Los Sistemas Oscilatorios Iineales periédicos han sido ampliamente estudiados
experimental y teéricamente en la fisica y sirven de modelos para estudiar las
interaociones entre sistemas, en las }401ltimasdécadas |os Sistemas
amortiguados y forzados son usados como modelos con frecuencia ; sin
embargo, el interés por conocer mas sobre los Sistemas no Lineales
Amortiguado y Forzado , y en especial saber las condiciones cuando muestran
031 comportamiento desordenado o Caos es lo que se busca(Douglas G., 2006).
\
031 Es importante determinar cuéndo un sistema presenta comportamiento de Caos
031 porque permite conocer esos comportamientos indeseables de los sistemas y,
de los sistemas flsicos que siwen de modelo.
El problema identi}401cado
es la detenninacién de| estado de Coas en Ios
030 Sistemas no Lineales que muestran comportamiento desordenado. Se plantea
la determlnacién del estado del Caos usando los Exponentes de Lyapunov,
propuesto por el ruso matemético del mismo nombre que propone que el
cambio en la dinémica se puede medir a través de la divergencia de las
variables dinémicas, que es cuanti}401cada

por unos parémetros llamados
Exponentes de Lyapunov (WolfA.,1984).
La formulacién y planteamiento del Problema son:
a) Problema General
Problema general que se propuso en este 'trabajo dé investigacién es el
siguiente,
g,Es posible determinar el Caos en la dinémica de un Oscilador no Lineal
Amortiguado y Forzado con los Exponentes de Lyapunov?
b) Problemas Especificos:
' - (,Cuél es la técnica numérica apropiada para calcular los Exponentes
de Lyapunov en la determinacién de| estado dinémico del Caos de un
Oscilador no Lineal Amortiguado y Forzado?
- ¢;Es posible determinar un modelo fisico para estudiar el Caos en la
dinémica de un Oscilador no lineal Amortiguado y Forzado?
N 9
c) Objetivo General y Alcances de la investigacion
-Objetivo General
Los objetivos planteados en el presente trabajo fueron:
Exponentes de Lyapunov y determinar el estado de| Caos en la dinémlca
o Desarrollar y calcular los Exponentes de Lyapunov de un Oscilador
-Objetivos Especificos
Exponentes de Lyapunov y determinar el estado del Caos en la dinémica
o Desarrollar y calcular Ios "Exponentes de Lyapunov de un Oscilador
3.2. lmportancia y la Justificacion de la lnvestigacién
Sabemos que el estudio de los Sistemas Fisicos no Lineales se ha
_ lncrementado y el problema de Caos se observa en varios sistemas, no
solamente encontramos Sistemas no Iineales en Fisica sino también en otras
areas como Quimica, Biologia, Economla, Medicina, lngenieria y por ello se
dice que el Caos en ubicuo. Los trabajos de modelamiento que realizamos, |os
fisicos, son interdisciplinarios y por ello surge la necesidad de estudiar y
comprender |as herramientas para afrontar el fenémeno del Caos. En ese
sentido es importante que los modelos no Iineales, en especial los mecénicos,
_ 030 que proponemos para estudiar otros fenémenos de la naturaleza debe ser
bien estudiado, saber cuando pueden funcionar en estado de Caos
X 10
031 controlando Ios parémetros propios del sistema para saber |as condiciones
que Ilevan al Caos.
o La investigacién que deseamos realizar corresponde a un valor te6rico 024
024 préctico, desde el punto de vista de una Investigacién Bésica, y constituiré un
' aporte importante dentro de la }402sica

computacional, Ia cual, es una de las
Iineas de especializacién de Fisica.
030 030 11
IV. MARCO TEORICO
_ 4.1. Antecedentes o Estudios Previos
' Zhang e Harne (2016) propone el estudio del Oscilador con efecto
débilmente no lineal, Ilamado Oscilador de Duf}401n,

muestra que este
Oscilador presenta comportamiento de Caos para determinados parémetros.
Cabana e Laroze (2015) muestra el anélisis de las ecuaciones de Landau-
Lifshitz o Landau-Lifshitz-Gilbert, que describe Ia dinémica de una particula
en un campo magnético oscilante, se analiza la dinémica hacienda variar
parémetros del Sistema y se estudia el Caos en su espacio de fase.
Mouronte M (2004) , muestra él estudia el Caos del sistema, a través del
espacio de fase donde |as bifurcaciones que son desdoblamiento de la
dinémica de ,es una forma cualitativa de describir el Caos en la dinémica de
un sistema. 030
Kapitaniak T. (1994), muestra que estudio del comportamiento de los
sistemas cuasi periédicos forzados a través de los Atractores, que son
regiones en el espacio de fase o oonfiguracién hacia donde tiende la
dinémica, es una forma de estudio de los sistemas cuasi periédicos.
La revisién bibliogréfica realizada muestra que en los trabajos no usan los
V exponentes de Lyapunov especi}401camente

para el estudio de| estado el Caos
de un sistema de un Oscilador Amortiguado y Forzado. El antecedente
descrito ha motivado a desarrollar Ia presente investigacién como una
contribucién dirigida a dotar de los instrumentos necesarios de consulta,
actualizacién y orientacién bésica para aquellos que tienen interés de esta
especializacién dentro de la Fisica no Lineal.
12
i
4.2. Bases Teoricas
4.2.1. La Teoria de| Caos
Usualmente en }401sica
aprendimos que las leyes son descritas por ecuaciones
de movimiento (ecuaciones diferenciales), como por ejemplo la Iey de
Newton entre otros, y la solucién analitica o numérica de estas ecuaciones
se hacen teniendo |as condiciones iniciales y de contorno; podemos conocer
como es su dinémica en el tiempo, esto es el punto de vista determinista,
entonces, sistemas fisicos conservativos o disipativos son deterministas
porque elios obedecen ecuaciones diferenciales deterministas (Taylor J.,
2013).
La teoria de Caos, es una teoria que estudia el comportamiento desordenado
de los sistemas no Iineales , como |os fisicos, a los que se conoce su
031 ecuacién dinémica resulto de la aplicacién de una iey o principio, en ei caso
de la Mecanica por ejemplo las Leyes de Newton; por ello se dice también
Caos determinista, |os sistemas no Iineales en régimen pertubativos
usualmente exhiben comportamiento complicado Ilamado Caos(José J. e
Saletan E.,1998).
Si |os estados iniciales de sistemas deterministas son exactamente
conocidos, futuros estados dinamicos son predicidos teéricamente; esta
visién determinista sobrevivio para ser cuestionada después de las famosas
obras visionarias de Henri Poincaré sobre la mecénica celeste realizadas a
}401nales
del siglo XIX que eran contran'o al determinismo(Ra}401ada
A.,1990).
Poicaré proponia que supiéramos exactamente la Iey de la naturaleza y la
situacion del universo en el momento inicial, se podria predecir exactamente
Ia situacion de ese mismo universo en un momento sucesivo, pero incluso si
fuera ei caso de que las leyes naturales ya no tuvieran ningiin secreto para
. nosotros, pero no es asi en los resultados numéricos que nos describen que
& 13
solo .podn'amos conocer Ia situacién inicial aproximadamente. Si nos
permitiera predecir la situacion siguiente con la misma aproximacién, eso es
todo lo que necesitamos, de esta manera, las mediciones efectuadas sobre
el estado de un sistema en un momento dado pueden no permitirnos predecir
la situacion futura incluso moderadamente lejos, a pesar de que las
ecuaciones de gobierno son exactamente conocidas (Douglas R., 2006).
No ocurrio hasta |os a}401os

70, después de casi un siglo de investigaciones,
aprendieron que los sistemas caéticos solo pueden resolverse
numéricamente, y no hay formas sencillas y generales de predecir cua'ndo un
sistema exhibiré caos. También hemos aprendido que el caos determinista
esté siempre asociado con sistemas no Iineales, Ia no Iinealidad es una
condicion necesaria para el caos, pero no una su}401ciente(José

J. e Saletan
E.,1998).
Solo en los }401ltimos

a}401os,
con el advenimiento de las computadoras que se
permitio estudiar el caos ya que ahora es posible realizar célculos numéricos
de la evolucién temporal de las propiedades de los sistemas sensibles a las
condiciones iniciales, comenzamos a comprender la existencia de| caos
cuando |as computadoras estaban fécilmente disponibles para calcular Ias
historias de largo tiempo neoesarias para explicar el comportamiento
discutido, el modelo deterministico siempre produciré la misma salida de una
condicion inicial o estado inicial determinado(Doug|as R., 2006).
4.2.2. indicadores de| Caos
La teoria de| caos, es una teoria que trata sobre el fenomeno del
comportamiento desordenado que los sistemas fisicos no Iineales presentan
, en especial, muestran una fuene sensibilidad a las condiciones iniciales; Ia
no Iinealidad de las ecuaciones (deterministas) hacen que no se puedan
V resolver analiticamente porque no cumplen con el principio de superposicion,
es por ello que en esta teoria se han desarrollado indicadores de| Caos para
\ 14
determinar si el sistema muestra el caos , entre los indicadores tenemos :
_ Espectros de Potencias y Super}401cies

de Seocién de Poincaré y Exponentes
' de Lyapunov(Doug|as R., 2006).
a) Seccién de Poincaré
Este indicador es bésicamente el mapa que se muestra al interoeptar las
1 curvas de las trayectorias, tomadas periédicamente. este procedimiento es
més cualitativo para averiguar si el sistemaVpresenta Caos(Douglas R.,
2006).
FIGURA N° 4.1
A SUPERFICIE DE LA SECCION DE POINCARE
(3) H 030?/" 034"| 0343 030v
3&1 ~
-}401?:§
mo» 2
53$
. _ _ _ $9; r
:<-on
Fuente: Ra}401ada.
Dinémica CIésica.1989
b) Espectros de Potencias
Aqui lo que se hace es determinar el médulo de la transformada de Fourier
en funcién de las frecuencias. una sistema oscilatorio lineal muestra un
Respectro de potencia de un solo pico, de un sistema oscilatorio cuasi

» 15
periodico dosomés pioo,pero siempre discretos latransformada de Fourier; I
sin embargo, cuando un sistema }401sico

esté en estado de caos debe mostrar
una transformada de Fourier continua, para ello se tiene que hacer variar
todos los parémetros y variables del sistema y no siempre se puede
determinar el limite a partir de| cual el sistema presenta caos (Ra}401ada
A.,1990)
FIGURA N°4.2
I ESPECTRO DE POTENCIA DE UN SISTEMA PERIODICO, CUASI
PERIODICO Y EN CAOS
030 ; I. q -_ .; ;:',-;2;": 030_ 031r::-_;'- .
o . .. - .- ,. 3 o _v ». '
Fuente: Ra}401ada,
Dinémica 030CIésica,
1989
c) Exponentes de Lyapunov
Los términos no Iineales de las ecuaciones diferenciales ampli}401can
diferencias en las condiciones iniciales, de esta

exponencialmente peque}401as
manera, Ias leyes de evolucion deterministas crean comportamientos
caéticos, incluso en ausencia de nuevas }402uctuaciones

internas, en el
régimen caotico no es posible predecir exactamente Ia evolucién del estado
del sistema durante un tiempo arbitrariamente, esta es la caracteristica
impredecible del caos.
\ b 16
La evolucién temporal se rige por un espectro continuo de frecuencias
responsables de un comportamiento aperiédico, para este movimiento se
presentan patrones estacionarios, es decir. patrones que se repiten de forma
no periédica. Lyapunov creé un método conocido como exponente
caracteristico de Lyapunov para cuanti}401car

la sensible dependencia de las
condiciones iniciales de| comportamiento caético,A la estabilidad de los
sistemas dinémicos; con este método es posible determinar |os requisitos
minimos de las ecuaciones diferenciales que son necesarias para crear el
caos. Estudiemos el caso de sistemas con una sola variable que asumen dos
estados iniciales x., y x., + s, di}401riendo

por una peque}401a
cantidad £(Ra}401ada
A.,1990)
Queremos investigar los posibles valores de x,, después de n iteraciones de
los dos valores iniciales, Ia diferencia d,. entre los dos valores xn después de
n iteraciones (omitiendo por simplicidad el subindice) se da
aproximadamente.
d,, = f(x',, + 5) 024-

f(x,,) = sen) 030 (1)
Donde, A, es el exponente de Lyapunov que representa el coe}401ciente

del
crecimiento exponencial promedio por unidad de tiempo entre los dos
estados. De lo que vemos. si A es negativo,'|as dos orbitas eventualmente
convergerén, pero si es positivo |as trayectorias préximas divergen. La
diferencia entre los dos estados iniciales se escribe como
«:1 = 034X0
+ s) ~ r(x.,) 2 5(: 024:)KO (2)
Ahora, para evitar Ia confusién que a veoes se encuentra en la literatura
caética, recordamos que,
\ 17
Xn+1 = f(Xn)=-f(f(Xn-1)): f(f(f(Xn-2)))= = f(f(f(-»-(f(Xo)))))- (3)
Donde el superindice n indica Ia n-ésima iteracién de| mapa. Después de
un gran numero de iteraciones. Ia diferencia entre los estados cercanos
es. 4
I d = f(x,1 + 2)
034 f(x,,) = f(x0 + 2)
024 f}401xo)
024 = sen" (4)
Dividiendo Ia ecuacién por s y tomando el Iogaritmo a ambos lados,
obtenemos I
In{[f n (X0 + 8) _
f S034
031 030_O
=
In[e 21}.] =
nix
-
(5)
_ Mxo) é ilnldf (6)
Donde,
«u M, {T}{-- 024m(,,
}-..{ 024m
-nio_ 024dF1(X;1_1) dFl{Xn_7) dfutxo)
& 18
Reempiazando la ecuacién 7 en la ecuacién 6 para calcular el exponente
de Lyapunov, que puede tener valores positive 0 negativo, el valor positivo
que el sistema presenta comportamiento caético.

signi}401ca
I = lim (5)1 N dfil 8

(X0) s 024nan H" dxo ()
4.3. Ecuacién de Duf}401ng
' Un caso ilustrativo disipativo es el movimiento de una particula con masa m
sometida a un potencial y a una fuerza disipativa es proporcional a, pig
respectivamente; se considera débilmente no lineal cuando, 6<<1.
dzx dz 3
024,+B 024 024x+8x
=FCos(mt) (9)
dt dt
Donde (3. 5 , on son parémetros.
19
FiGURA N° 4.3
OSCILADOR AMORTIGUADO FORZADO LINEAL Y SU ESPACIO DE FASE
(3) (b)
6.0 6.0
030 = * if u 031
U
""lt:uJ ; 5 030 ' _ / 03403 - » - 030vtcm

»-
030(Swan
031 4035.0 u)6u.o
031 ""94.a 0.0 4,0
I x
Fuente: Mouronte M., Caracterizacién y Anélisis de Sistemas Dinémicos no
Lineales Mediante el Estudio del Mapa de frecuencias y Espacio de fases, 2004
FIGURA N 4.4
034 '
OSCILADOR AMORTIGUADO FORZADO NO LINEAL EN CAOS Y SU
ESPACIO DE FASE
(a) (b) '
is .- J ' I xi}
'»<o.q _ ; I _. 024v;o.f;._ 151.: .

031 I ' § E _V 031 . _ 030_ 030/j-.'.~.¢?
030 §_,5_
030;~ 030_x..
034an 1
4 _ I '7 O :.%
l'8l0.0 1055.0 no0.o 24,0 my .31)

t x
Fuente: Mouronte M., Caracterizacién y Anélisis de Sistemas Dinémioos no
Lineales Mediante el Estudio de| Mapa de frecuencias y Espacio de fases,
2004
20
La Dinémica esté compuesta por dos campos de estudio, Ia Dinémica Lineal
y la Dinémica no Lineal, esta Liltima trata del comportamiento desordenado
dinémico de los Sistemas no Lineales, es aqui donde Ios sistemas presentan
el Caos, a mejor dicho el Caos determinista, porque cuentan con ecuaciones
dinémicas con fuerte sensibilidad a las condiciones iniciales. Esto signi}401ca
que una peque}401a

variacién en las condiciones iniciales puede alterar
signi}401cativamente
el comportamiento dinémico del sistema convirtiendo los
valores que va a tomar en ei tiempo (Saletan & José, 1998).
Otros sistemas también se puede Ilegar a la forma de| Oscilador de Duf}401ng,
si consideramos un cuerpo de masa ,M, ai cual le llamaremos la particuia(
sin rotar), suspendido entre dos muelles iguales, tal como se muestra en la
figura 5. Si ambos muelles estén sin de deformar(o sea, sin tensién y, por
tanto; sin energia potencial los dos lados) cuando Ia panicula se encuentre
en su posicién de equilibrio y, si despreciamos Ias fuerzas gravitatorias, cada
uno de los muelles ejerceré una fuerza,
, 024x(s
l)
024
sobre la particula (x es la constante de recuperacién de cada mueile), 1 , es
la Iongitud sin deformar dei resorte. (Saletan & José, 1998).
\ 21
_ FIGURA 4.5
SISTEMA NO LINEAL COMPUESTO DE UNA PARTlC_ULA DE MASA M
" ~
Fuente: Marién J. Dinémica Clésica de las partfculas y Sistemas, 1981
FIGURA 4.6
SISTEMA NO LINEAL COMPUESTO DE UNA PARTICULA DE
MASA M MOVIENDOSE LONIGITUDINALMENTE
K 030 E \\: 030
______""_ .._'f.__\;,:M
K 030er
/ , 031
s
' Fuente: Marion J., Dinémica Clésica de las particulas y Sistemas, 1981
- 22
La fuerza que acttia sobre la particula de masa ,M, es ,F, igual a,
F ' 0242z:(.-:l)$'crz(9),
024
como : = wxz + I3
de forma que,
_ 9 x I
Jen = 024
= 024 024 024 024 024
5 \/x3 + :2
Resulta la fuerza en funcién de la posicién,
F= 1 024.'+
\;A 035
(v 030x2
030
e :2 L) = zzcc 1 1
F (10)
H031031 V1 +G} 031
Si suponemos que, . es peque}401o

y desarrollamos en serie el radical,
tendremos(Marion J.,1981),
12' _,cl( 031£)3

1_"L(§.)2 (11
_ 030
1 4 .1 )
23
Y cuando despreciamos todos los términos, salvo el primero, seré
k
F 5-! (11))?
Si alargamos inicialmente ambos muelles una distancia, d ,para sujetar a
ellos Ia masa , m , en posicién de equilibrio, la fuerza seré (Marion J.,1981),
' _ vcd N 3 3
» F(X)§ 0242 x030
024[rc(l
024l- 024d,/I
]x
De donde podemos determinar Ia energia potencial considerando la
expresién,
F _ dU(x)
_ dx
_ 55 _ xii-d) 3 _ __.s'u-{ac}
A 2(1)x [ F J) 030
024 .2; ' (12)
Obtenemos la energia potencial siguiente,
_E E 030
U(x)-z(2(l'))x i[j""F 034"
2 +4 x.
4 031] (13)
24
Luego considerando que el sistema se mueve solamente en forma horizontal
la energia cinética viene dado por , '1 030=}M:*

y conjuntamente con la
energia potencial , U(x) , podemos expresar el Lagrangiano del sistema
L=5M72'2 0245(2(§))x2 0245F3:£]x 030*

Z Z I
(14, 4- A 030
Cuando ademés estén presente otras fuerzas, disipativas y no disipativas en
Mecénica Clésica se usa la siguiente ecuacién,
d ''at.\ (al.) 1 2
030 " T 024 024
= 030
- + ~ 15
d: (.6471) aa, 0 Q3
031 ( )
Donde:
Q} _ representa a las fuerzas disipativas,
Q} , representa |as fuerzas que no son disipativas ni conservativas.
En nuestro caso usaremos consideraremos que nuestra particula esté
sometida a una fuerza disipativa de la forma, Q5 =

030
024b( 030:y 024
una
030:)
fuerza
externa dependiente del tiempo, Q)-2 to (.'os(ut)

024 ambas en direccion
horizontal.
\ 25
030 ' Reemplazando el Lagragiano en la ecuacién dinémica y las fuerzas
respectivas se obtiene Ia siguiente expresién,
3 1 1 -,2 030
_:(_; .2 _1li((§-t}402j
.4) =_ 030g
¢:(a2(zM? 034 any 2(2 J1 4 :3 1 b(a't)+
' I-'0 Cos(wt)
(16)
030 desarrollando Ia expresién de la manera que se muestra para obtener una
ecuacion dinémica de un Sistema no Lineal Amortiguado y Forzado,
+2 .1 +
030 = 024b¢_ +034§)
F0 Cos(wt)
(17)
A la ecuacién,
3 ii
}402 030
E . L-#3 .a_5 .
1:1 + M (dc) + 2 (in) X + [ 131:: J1 " M COS(wt) (18)
_ 26
Se hacen |as siguientes representaciones,
E ,= 55 _. [ 031 =030 030 030 031' 030')I

M 2/3» 2(m) 024wa ~
2
6-
' nos permite expresar de la siguiente forma Ia ecuacién dinémica,
- i . F
.1: + 23:: + (A)(j.22. +
0316x3 = B 024:_ 031Cas(_wt) (19)
De donde observamos que s << 1

030 es bien peque}401o
por lo tanto podemos
a}401rmar
que la ecuacién dinémica corresponde a un Sistema débilmente no
lineal Amortiguado y Forzado.
27
4.4. Marco conceptual
a) Caos
Es un fenomeno, que provoca Ia imposibilidad de realizar cualquier prediccion
}401able
de su comportamiento a pesar de de tener un carécter determinista, aun
en ausencia de las perturbaciones extemas al sistema, los temas de esta
clase desarrollan una comportamiento caétlca (Alemén R., 2010).
b) Mecénlca no lineal y Caos
La mayoria de los sistemas cuyos ecuaciones de movimiento son no Iineales
pueden exhibir Caos. Este fenémeno que se mani}401esta

en muchos campos
diferentes como los sistemas mecénicos oscilanles, reacciones quimicas, }402ujo
de }402uidos,
rayos léser, crecimiento de poblaclones, transmislén de
enfermedades , y muchos otros, supone que aunque un sistema obedezca a
ecuaciones deterministas de movlmiento(|eyes de Newton) su comportamiento
detallado futuro puede, en la préctica, ser impredecible(Taylor, J., 2013).
El comportamiento de un sistema mecénico caético puede ser muy
complicado, y la necesidad de describir este comportamiento ha generado una
enorme variedad de maneras nuevas de visualizar el movimiento de los
sistemas como érbitas del espacio de estados, secclones de Poincare,
diagramas de bifurcacién.
La teoria del Caos esté dividida en dos amplias éreas de estudios, sistemas
disipativos, tales como el péndulo amortiguado, y no disipativos, o sistemas
Hamiltonianos (Taylor, J., 2013).
\ 28
' c) Linealidad y no Iineaiidad
La ecuacion diferencial es lineal si relaciona Ia variable 0 variables
dependientes y sus derivadas solo Iinealmente. El francés Poincare observo
que la tierra, la Iuna y el sol. sus ecuaciones exhiben sensibilidad a las
condiciones iniciales. Caracteristica de| movimiento caético.
La no linealidad es esencial para el caos, pero la no linealidad no garantiza el
caos ademés el sistema tiene que ser 034oomplicado

como
035 |os sistemas
disipativos, sus ecuaciones de movimiento cuando se ponen en su forma
esténdar autonoma deben ser no Iineales y tener N variables con N>=3(Tay|or,
J., 2013).
\ 29
5.1. Materiales
5.1.1. De ejecucién
' o El trabajo de investigacién pertenece a la investigacién teérica,
especi}401camente
a la Fisica computacional, por ello se usé ordenadores
Core I-3 para hacer las simulaciones y célculos, programas en fortran
90 para Runge Kutta de 410 orden (Press W., Teukolsky 3., Vettering
W., Flannery B.,1992) para estudiar el comportamiento de la dinémica
de| espacio de fase.
- Internet . para la busqueda de informacién relacionado al trabajo de
investigacién.
o Fotocopias, de documentos que ayudan en el desarrollo de| trabajo de
investigacién.
o Dispositivos de almacenamiento (CD, USB) . para guardar Ias corridas
con el programa.
5412! De la impresibn
Para Ia parte }401nal

del trabajo de investigacién se uso los siguientes
materiales:
o Papel Bond en formato A-4
o Elementos de escritorio
o Tintas para impresién
X 30
5.2. La poblacién y muestra
Como es un trabajo de investigacién de Fisica Teérica no requiere de
poblacién y muestra.
5.3 Las
030 técnicas, procedimientos e instrumentos de recoleccién de datos
No se usan técnicas estadisticas por tratarse de un trabajo de la Fisica
Teérica.
031 5.4. Las Técnicas de anélisis o métodos estadisticos aplicados
No aplica para este trabajo.
5.5. La Metodologia
-En este trabajo se consideré el siguiente sistema, }401gura

4.7,
FIGURA 4.7
SISTEMA N0 LINEAL COMPUESTO DE UNA PARTICULA DE
MASA M y RESORTES ESTIRADO X
K .: \\\s
M V L X \ \«
030
__.____ ________., 030M
I. 030 9 030I;
¥ Fuente: Marién J., Dinémica Clésica de las particulas y Sistemas, 1981
31
024La
recopilacién de datos se hace a través de un programa mencionado
Iineas arriba de la ecuacién de la dinamica del sistema,
11 ax 024.
024f
+ 2;? (- 024)
+ W 030):
+ 6.1.63 = F(.'ns(mt) , (20)
dz 030 dt
y el programa para obtener |os exponentes de Lyapunov fue adaptado de
(\NoIfA., Swift J., Swinney H. and Vastano J., 1985).
Para ello a la ecuacién ultima se convierte en un sistema dinémico de dos
variables, donde Ia segunda ecuacién contiene el término que contiene al
tiempo, I-'Cas(wt),
dz: _ '
dL 3
iii _ _9 _ 2 __ 3 ,
at ..,B at w 3: 6x + I'Co.s(cot)
(21)
V El Sistema. se hizo discreto y se Ileva al programa de Runge Kutta de 4to
orden con intervalo de tiempo 0,015, para analizar si se observa
comportamiento desordenado o Caos en el espacio de fase para con los
parametros de la tabla N° 6.1.Para esos valores se calculan |os exponentes
de Lyapunov méximo positivo, que mostré Caos.
K 32 .
l
i
Luego haciendo un cambio de variable:
Del tiempo, "L = w,,L, ~
de la sustitucién de la frecuencia, fl 024

3 .La ecuacién adimensional
Wu
resulta,
2 2 .
i+2_(3wo +wD2x+.5:c =iCvs

031 (£1) (22)
:12 -1: M me
dia: 23 das _ 5 _3 Po
- 024;
+ + It + 024zx
= :2-C0S(QT) (23)
03 it ma dr mg mg, M
, . . . 2:3 031 2
-En el analisis se mantuvo constante |as cantidades , :3
024 ,
OJ.) 0
i y solamente se hace variar Ia frecuencia relativa Q .

0
Se determiné Ios espacios de fases para las diferentes frecuencias desde,
m = }402,}402
hasta
030lma
m = 3,4 mu con los datos que se muestra en el cuadro
. . . _ Z}402 s Fa
sugulente, mantenlendo constante las cantidades - 024,024 024 024g-, 030T
.
can an we M
\ 33
VI. RESULTADOS
TABLA N 6.1
034 .
VALORES PARA LOS DIFERENTES PARAMETROS DEL SISTEMA NO
LINEAL AMORTIGUADO Y FORZADO
1 9031 9 = I;
034O2 1:932!" - W0
N0 linealidad Fuerza externa Frecuencia
034E111
Ejij
E
u
j
j
j
j
ET
Fuente: Autoria propia
\ 34
GRAFICA N 6.1
035
0.215798E+O1
- 20 _ ""_' 034 034

'sa!ida35,dat' u 23 034
- 024- 024
o
\ -1.. 030 I, f
'5 \$-.\ U 030
\ 034\1 030:,-
-10 V ,
in
now
030 2 4 5 s in 12 14 16 13

Gnuplot.
\ 35
GRAFICA N° 62
20
M an 024 024
0305 1" "\\\

030
g /5" ,,
1/ 03 035 'z\ \_r "#7 031 «
5 030f 03
1: I r 030. )
.» 1 « 2 c, 030
030S % . . \ /I
\ \ V. . 030
.10 . V/4,4, *" 0307' M
\ 030/'
,;2a_-2_oM -15 -10 -5 0 5 to 15 20

.324538E+02

Gnuplot.
\ 36
GRAFICA N° 6.3
Y FORZADO, FRECUENCIA 0,08 Y EXPONENTE LYAPUNOV NEGATNO
20 /\
15 H _'__,A _a_:_ Rt 034 -1if\* ;*;-\:\\ 030
,4- 030? 030. 031 *;§\\_\'_\ 030030'7///

030 030
AA\\.~v~'v*~~A
f W
° 1: I)))))! 031) ) 030l)
Va 1* 3"! 031 034031
034'
. ' A \
-212 H
_ _-1s__A_ -1o -5 0 5 1a 15 20

Gnuplot.
X 37
n I// ., h 034N§Cn:gRT|
5 0301 030T;33:::9P»:~§
\;:\ 031_;.~.MW§ 030
1 034 .
Q sra}401h}401q}401r
035 5 %
agbtenidoks 030
2° }402
. "5 um»
Ill 030-I 031
\\V"77>
I K
1 030
\,/
V 030
4% 030:-I
"
X 39
g / %%
05$ 035?
g9F 031:Fof 030 /1
\;
4] 030? 20
A9
%
VII. DISCUSION
El espacio de fase del Sistema no Lineal Oscilatorio Amortiguado y Forzado que
observa de las }401guras

6.1,6.2,6.3 muestran una dinérnica desordenada, las Iineas
se cruzan varias veces , la dinémica se realiza alrededor de dos centros de
oscilacion, que se les llama atractores, |as frecuencias relativas corresponde a
Q=0.1, 0.3 , 0.6 respectivamente.
El comportamiento observado aqui del oscilador no lineal, también Io podemos
observar en el Oscilador no lineal Amortiguado y Forzado de un péndulo donde Ia
forma de la ecuacion dinémica se aproxima al sistema usado en nuestro trabajo de
investigacion.
Los exponentes de Lyapunov determinados para esos tres casos resultaron
teniendo con valores positivos, resultado que esta en la direccién de demostrar de
la hipotesis de este trabajo de investigacion, que los exponentes de Lyapunov son
también indicadores de| estado de Caos del Sistema no Lineal Oscilatorio
Amortiguado y Forzado.
Los resultados obtenidos en este trabajo va en la misma direccion con otros
trabajos de investigacién, donde usan también al exponente de Lyapunov como
indicador, pero para otros sistemas no Iineales (Cabanas A.,Laroze D.,2014).
\ 41
r
SUBROUTINE derivs(x,y.dydx)
REAL x,y(").dydx('),A,C,E,D
A=0.01
c=o.o2 030 .
D=3O.0
E=0.08
dydx(1)=y(2)
! dydx(2)=-((w1"2) 030y(1))
dydX(2)=-A"Y(2)-Y(1)-C*((Y(1))**3)+D 030C0S(E*X)
I dydx(3)=y(4)
I dydx(4)= 024((w2**2)*v(3))
return
END
! rkdumb.for
SUBROUTINE rkdumb(vstart,nvar,x1,x2,nstep,derivs)
INTEGER nstep,nvar,NMAX,NSTPMX
PARAMETER (NMAX=50,NSTPMX=100000)
REAL x1 ,x2,vstart(nvar),xx(NSTPMX),y(NMAX, NSTPMX)
EXTERNAL derivs
COMMON lpathl xx,y
!CU USES rk4
INTEGER i,k
REAL h,x,dv(NMAX),v(NMAX)
do 11 i=1,nvar
v(i)=vstar1(i)
V(i.1)=V(i)
11 continue
xx(1)=x1
44
031 x=x1 030
h=(x2-x1)/nstep
do 13 k=1,nstep
call derivs(x,v,dv)
call rk4(v,dv,nvar,x,h.v.derivs)
if(x+h.eq.x)pause 'stepsizé not signi}401cant

in rkdumb'
x=x+h
xx(k+1)=x
do 12 i=1,nvar
y(i,k+1)=v(i)
12 continue
13 continue
return
END
! rk4.for
SUBROUTINE rk4(y,dydx,n,x,h,yout,derivs)
INTEGER n,NMAX
REAL h,x,dydx(n).y(n),yout(n) '
EXTERNAL derivs
PARAMETER (NMAX=50)
INTEGER i
REAL h6,hh,xh,dym(NMAX).dyt(NMAX),yt(NMAX)
hh=h"0.5 -
V h6=h/6.
xh=x+hh
do 11 i=1,n
Yt(i)=Y(i)+hh*dYdX(i) '
11 continue
cail derivs(xh,yt,dyt) .
do 12 i=1,n -
Vt(i)=y(i)+hh"dyt(i)
45
12 continue
call derivs(xh,yt,dym)
do 13 i=1,n
Yt(i)=y(i)+h*dYm(i)
dym(i)=dyt(i)+dym(i)
13 continue
call derivs(x+h,yt,dyt)
do 14 i=1,n
yout(i)=y(i)th6 031(dydx(i)+dyt(i)+2. 030dym(i))
14 continue
return
END
9.2.PROGRAMA DE EXPONENTES DE LYAPUNOV ADAPTADO DE (WOLF A.
1985)
program ode
nn=n"(n+1)=
nn1=nn+1
parameter (n=3, nn=12, nn1=13)
implicit double precision(a-h,o-z)
external fcn
dimension y(nn1),yprime(nn1),v(nn1),A(nn1),B(nn1),C(nn1)
dimension D(nn1),cum(n),znorm(n),gsc(n)
open(1,}401|e="Lout.txt")
111 format(1x,'nstep, stpsze, irate, io:') .
read(*,*) nstep,stpsze,irate,io
write(1,*) nstep,stpsze,irate.io
write(1,*)
v(1)=10.0
v(2)=1.0
v(3)=1.0
K 46
tme=0.0
do 10 i=n+1,nn
v(i)=0.0
10 continue
do 20 i=1,n
v((n+1)*i)=1.0
cum(i)=0.0
20 continue
do 100 m=1,nstep
if ((m/100)*100.EQ.m) wn'te(*,*)m
do 25 j=1,irate
do 26 i=1,nn
030/(i)=V(i)
26 continue
t=tme
do 27 i=1,nn
A(i)=yprime(i)
27 continue
do 28 i=1,nn
y(i)=v(i)+(stpsze*A(i))/2.0
28 continue
t=tme+stpszel2.0
do 29 i=1,nn
B(i)=yprime(i)
29 continue
do 30 i=1,nn
y(i)=v(i)+(stpsze'B(i))/2.0
30 continue
t=tme+stpsze/2.0
X 47
do 31 i=1,nn
C(i)=yprime(i)
31 continue
031 do 32 i=1,nn
y(i)=V(i)+(StPSZe*C(i))
32 continue
t=tme+stpsze
call fcn(t, y, yprime) 030
do 33 i=1,nn
D(i)=yprime(i)
33 continue
do 34 i=1,nn
v(i)=v(i)+stpsze*(A(i)+D(i)+2.0*(B(i)+C(i)))/6.0
34 continue
tme=tme+stpsze
25 continue
znorm(1)=0.0
do 38 j=1,n
znorm(1)=znorm(1)+v(n*j+1) 030*2
38 continue
znorm(1)=sqrt(znorm(1))
do 40 j=1,n
v(n*j+1)=v(n*j+1)/znorm(1)
40 continue
do 80 j=2,n '
do 50 k=1 ,1-1
gsc(k)=0.0
do so |=1,n _
gsc(k)=gsc(k)+v(n*I+j)*v(n*|+k)
50 continue
\ 48
do 60 k=1,n
do 60 l=1,j-1
v(n*k+j)=v(n*k+j)-gsc(|)"v(n*k+|)
60 continue
znorm(j)=0.0
do 70 k=1,n
znorm(j)=znorm(j)+v(n*k+j)**2
70 continue
znorm(j)=sqrt(znorm(j))
do 80 k=1,n
v(n*k+j)=v(n*k+j)/znorm(j)
80 continue
do 90 k=1,n
cum(k)=cum(k)+d|og(znorm(k))/a|og(2.0)
90 continue
if(mod(m,io).nei0) goto 100
write(1,334) tme,(cum(k)/tme,k=1_n)
334 format(1x,f12.6,2x,e12.6,2x,e12.6,2x,e12.6)
_ 100 continue
end
*m*m**..m*....Mmm...m..........m..............*........*
subroutine fcn(t,y,yprime)
b=4.0
sg=16.0
r=45.92
VPrime(1)=S9*(V(2)-Y(1))
yprime(2)=-y(1)*y(3)+r*y(1)-y(2)
VPrime(3)=Y(1)*Y(2)-b*Y(3)
do 10 i=0,2
vprime(4+i)=S9"(V(7+i)-Y(4+i))
vprime(7+i)=(r-y(3))*v(4+i)-v(7+i)-v(1) 030v(10+i)
\ 49
' 030 vr>rime(10+i)=y(2)*y(4+i)+y(1)*y(7+i)-b*v(10+i)
10 continue
. return
_ . end
031 V i X _ 50
/ X. ANEXOS
Matnz de Consistencia
FORMULACI 0 N DEL PROBLEMA OBJETIVOS HIPOTESIS VARIABLES INDICADORES METODOLOGIA
I I Detenninacion del Caos

030 (IE5 posible determinar el Caos I I Los Exponentes de
I I I en la dnnémnca de un I
en la dlnémlca de un Oscilador I I Lyapunov determnnan el
I I Oscilador no Lineal I I
General no Llneal Amortiguado y Forzado I estado dlnémlco del Caos
Amomguado y Forzado I I I I I I
con los Exponentes de en un Oscilador no Lineal Esta lnvestlgaclon es de tipo:
con los Exponentes do I I I
Lyapunov? L AmomguadoyForzado. Célculos aplloada. Con un dlse}401o de
ya unov. _ I I
p numéncos ejecuclon transversal. El
I I lndependiente: iguales o Sistema do anélisls es el

I I Desarrollar la tecnlca I I . _ I
¢;Cuél es la técmca numénca II I La técnica numénca Exponentes de mayores a Osculador no Lineal
I numenca apropnada para I
apropiada para calcular Ios _ basada en la Lyapunov cero de los Amomguado y Forzado; a
' calcular los Exponentes I I I I
Exponentes de Lyapunov en la d L ortonormahzaclén de exponentes este Oscilador no Lineal se
e a unov
determinacion del estado I y P y Gram-Schmidt local de Lyapunov formularé el Lagrangiano de|
I I determinar el estado del I I I I
dlnémlco del Caos de un I I perrmte determinar los sistema y luego la ecuaclén
I _ Caos en la dlnémlca de I I
I Oscllador no LinealAmomguado I _ Exponentes de Lyapunov dlnémloa de Lagrange, con
Especu}401c un OSCIIBUOF no Lnneal I I
y Forzado? ' I an un Oscalador no Llneal ello Implemenlaremos un
a Amortlguado y Forzado. I I
Amomguadoy Forzado, programa basado La técnica
I _ I I numérica basada en 031

la
.. (,Es posible detennlnar un Desarrollary calcularlos I I Dependlentez I
I I I El modelo }401sloo
Osculador I I I orlonormallzaclon de Gram-
modelo flslco para estudiar el Exponentes de Lyapunov I I Estado dunémuco Componamne I
I I I I no lineal Amortiguado y I I Schmidt local.
Caos en la dunémuca de un en un Oscllador no Inneal I I del Oscllador no mo Caolloo
_ I _ Forzado permite estudiar I I
Oscllador no lineal Amomguado Amortiguado y Forzado I Lnneal de la vanable
I la teoria del Caos con los I I I
yFotzado? en Resonancia. Amortiguado y dlnémlca.
Exponentes de Lyapunov.
Forzado.
51

Lévano Huamaccto If 2017

Cargado por

Copyright:

Formatos disponibles

Lévano Huamaccto If 2017

Cargado por

Información del documento

Descripción original:

Título 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

Lévano Huamaccto If 2017

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE CIENCIAS NATURALES Y MATEMATICA

UNIDAD DE INVESTIGACION DE LA FACULTAD DE CIENCIAS

11:: U 030£l'lcVEE':I£S 030I£AnNACl0NALDELCALLA0;

AMORTIGUADO Y FORZADO 035

AUTOR: Lic. CARLOS ALBERTO LEVANO HUAMACCTO

(Resolucién de aprobacién N°110-2017R)

' Callao, 2017

2.2. Abstract ' 7

3.1. Exposicién del problema de investigacién 8

3.2. Importancia y justi}401cacién

IV. MARCO TEORICO

4.1. Antecedentes o estudios brevios 12

4.2. Bases teéricas

4.2.1.La Teoria de| Caos 12

4.2.2. Indicadores de| Caos 14

- 4.3. Ecuacion de Duf}401ng 19

4.4. Marco conceptual 28

b) Mecénica no lineal y Caos 28

5.2. La poblacion de la investigacién y la muestra 31

5.4.Las técnicas de anélisis o Ios métodos estadisticos aplicados 31

Valores para los diferentes parémetros de|

Figura 4.1 Super}401cie

Figura 4.2 Espectro de potencia

Figura 4.3 Oscilador Forzado lineal y espacio fase

Figura 4.4 Oscilador Forzado No lineal y espacio fase Pag.2D

Figura 4.5 Sistema No Lineal Fisico Pag.22

Figura 4.6 Sistema No Lineal Fisico moviendo

Figura 4.7 Sistema No Lineal Fisico desplazado

Grafica 6.1 Espacio de fase sin Caos frecuencia 0.01 Pag.35

Grafica 6.3 Espacio de fase sin Caos frecuencia 0.08 Pag.37

En este trabajo de investigacién determinamos el estarid de Caos de un sistema

senoidal a través del indicador Ilamado exponentes de Lyapunov; para caracterizar

la dinémica de| Oscilador no lineal determinas primero el espacio de fase donde

L se observa si presenta comportamiento desordena para una frecuencia relativa

' determinada, luego calculamos su oorrespondiente exponente de Lyapunov a

exponentes de Lyapunov positivos para dinémicas desordenadas del Oscilador no

Lineal Amortiguado y Forzado y, valores de| exponente de Lyapunov negativos

para dinémicas ordenadas del Oscilador corroborando Ia hipétesis de este

We determine the state of Chaos of a non-linear oscillator damped and forced

called exponents of Lyapunov; characterize the dynamics of the non-linear

disordered behavior for a given relative frequency, then we calculate its

corresponding exponent of Lyapunov through a program developed by Wolf et al.

of the Negligible and Foroed Non-linear Oscillator and, negative Lyapunov

exponent values for ordered dynamics of the Oscillator corroborating the

hypothesis of this research work.

3.1. La exposicién del problema de investigacién

sistemas Iineales que se puede encontrar en la naturaleza ,y cuando se

estudian los primeros, se analiza su dinamica en aproximacién lineal, 0 sea las

variables contenidas en la ecuacién diferencial 0 de movimiento son de primer

orden, en esa aproximacion la ecuacién de movimiento del sistema mencionado

Iogra cumplir con el principio de superposicién, importante principio de las

ecuaciones Iineales; sin embargo, Ia aproximacion lineal Ileva a una descripcién

de la dinémica de los sistemas mencionados restringida, como en ei caso del

sistema oscilatorio periédico Ilamado Pénduio Fisico, que se ie estudia

menos de cinco grados a pesar que este

muestra movimiento en la practica para anguios mayores(Tay|or, 2013).

En |as Oltimas décadas se ha asociado el comportamiento desordenado que