Ecuaciones Diferenciales
Ecuaciones Diferenciales
Ecuaciones Diferenciales
) p = 0
Ecuacin lineal de donde obtenemos p. A continuacin se calcula u en
funcin de p, con lo cual y
2
(x)= u(x) y
1
(x) es una solucin de la
ecuacin diferencial original, siendo y
1
, y
2
linealmente independientes ya
que u(x) no es constante. Por tanto y
1
, y
2
forman un conjunto fundamental
de soluciones de la ecuacin original. La solucin general es de la forma y =
c
1
y
1
+ c
2
y
2.
Ecuaciones Diferenciales De Segundo Orden Con Coeficientes
Constantes
Esta ecuacin es de la forma
y'' + ay'+ by = q(x)
Supongamos q (x) = 0. Esta ecuacin tiene siempre soluciones del tipo
exponencial de la forma y= e
r x
. La sustitucin de esta solucin en la ecuacin
diferencial nos da
(r
2
+ar+ b)e
r x
= 0
5
Como e
r x
0, x, se tiene que r
2
+ar +b =0. Esta ecuacin se llama ecuacin
caracterstica de la ecuacin diferencial. Las races dan valores de r para los
cuales e
r x
es una solucin de la ecuacin.
Estas races son
=
Las races pueden ser: Dos races reales distintas. Una raz real doble. Races
complejas conjugadas
1. La ecuacin caracterstica tiene dos races reales distintas. Sean r
1
y
r
2
estas races, por tanto y = e
r
1
x
, y
2
= e
r
2
x
son soluciones de la ecuacin
diferencial. Estas soluciones son linealmente independientes en cualquier
intervalo por ser el wronskiano 0, luego y
1
,y
2
forman un conjunto fundamental
de soluciones, por tanto la solucin general de la ecuacin homognea es
y =c
2
e
r
1
x
+ c
2
e
r
2
x
2. La ecuacin caracterstica tiene dos races reales iguales. En este caso: a
2
-4
b=0. Una solucin de esta ecuacin diferencial es: y =
.
La solucin general de la ecuacin diferencial ser:
y() = (
3. La ecuacin caracterstica tiene dos races complejas conjugadas. En este caso:
a
2
- 4b < 0 . Sean:
=
, =
6
Las races son: r
1
= m +in, r
2
= m in . La solucin general de esta ecuacin
diferencial ser:
y= e
m x
(C
1
cos n x+ C
2
sen n x)
Ecuacin Diferencial De Euler
Esta ecuacin es de la forma:
x
2
y'' + axy'+ by = 0 , x >0
0
2 1
2
2
2
= + + y a
dx
dy
x a
dx
x d
x
Donde
2 1
, a a son constantes reales y
0
0
= a
.
Para resolver esta ecuacin se hace el cambio |x| = e
t
, con lo cual la ecuacin
diferencial dada se transforma en la ecuacin de coeficientes constantes.
En efecto:
Suponiendo , 0 > x y tomando
t
e x = t = ln x.
Entonces:
;
1
dt
dy
x dx
dt
dt
dy
dx
dy
= =
y
;
1 1
2 2 2
2
2
2
dt
dy
x x dt
y d
dx
y d
=
Sustituyendo en (3): 0
2 1 2
2
2
0
= + + y a
dx
dy
x a
dx
x d
x a obtenemos:
0
2 1
2
2
0
= + +
(
y a
dt
dy
a
dt
dy
dt
y d
a
Es decir: 0 ) (
2 0 1 2
2
0
= + + y a
dt
dy
a a
dt
y d
a ecuacin diferencial lineal con
coeficientes constantes.
Finalmente, resuelta esta ecuacin, se deshace el cambio y por sustitucin se
obtiene la solucin del problema dado.
7
y'' + (a 1)y' + by = 0
Aplicaciones Fsicas
El estudio analtico del comportamiento de muchos sistemas en Fsica
conduce a ecuaciones lineales.
Un ejemplo tpico en Mecnica son los osciladores armnicos:
1. El oscilador armnico simple: su ecuacin es
my
+ ky = 0
2. El oscilador armnico amortiguado: su ecuacin es
my
+ y + ky = 0
3. El oscilador armnico (amortiguado) forzado: su ecuacin es
my
+ y
+ ky = F (t)
En todos los casos se trata de ecuaciones lineales con coeficientes
constantes. (Su resolucin ha sido tratada en los cursos de Fsica).
El estudio analtico del comportamiento de muchos sistemas elctricos
conduce tambin a ecuaciones lineales. Un ejemplo tpico son los circuitos
elctricos RCL, cuya ecuacin es:
que es, tambin, una ecuacin lineal con coeficientes constantes.
8
Ecuaciones Diferenciales No Lineales
Homogneas de Segundo Orden
Una ecuacin diferencial de segundo orden es una ecuacin que implica la funcin
desconocida y , sus derivadas y', y'' de y, y la variable x.
Consideraremos solamente las ecuaciones diferenciales explicitas de la forma,
= (, ,
)
En general, poco se sabe sobre ecuaciones diferenciales de segundo orden no
lineales
= (, ,
= (,
)
Hgase el cambio de variable v = y. Entonces la nueva ecuacin satisfecha por v
es
= (, )
Esto es una ecuacin diferencial de primer orden. Una vez que se encuentre v su
integracin da la funcin y.
Ejemplo 1: Encuentre la solucin de
Solucin: Puesto que y falta, fije v = y '. Entonces, tenemos
9
Esto es una ecuacin diferencial linear de primer orden. Su resolucin da
Como v (1) = 1, conseguimos . Por lo tanto, tenemos
Teniendo en cuenta que y ' = v , obtenemos la ecuacin siguiente despus de la
integracin
La condicin y (1) = 2 da . Por lo tanto, tenemos
Observe que esta solucin est definida para x > 0.
Ecuaciones Sin La Variable X
= (,
)
Cambiamos v = y '. Entonces
10
Conseguimos
Esto es otra vez una ecuacin diferencial de primer orden. Una vez que v
entonces se encuentre podemos conseguir y a travs
Cual es una ecuacin en variables separadas.
Ejemplo 2: Encuentre la solucin general de la ecuacin
Solucin: Puesto que la variable v falta, fije v = y. Los frmulas arriba conducen a
Esto una ecuacin diferencial separable. Su resolucin da
11
Puesto que , conseguimos y ' = 0 o
Puesto que esto es una ecuacin diferencial separable, conseguimos, despus
de la resolucin,
,
Con las constantes y C. Todas las soluciones de nuestra ecuacin inicial son
Observe que debemos prestar la atencin especial a las soluciones constantes al
solucionar cualquier ecuacin separable. Esta puede ser fuente de errores...
12
APLICACIN A LA ECUACIN DE BALANCE DEL BUQUE
En el estudio del movimiento de balance del buque no vamos a disponer de una
funcin potencial de ninguno de los dos tipos por lo que nuestro acercamiento al
conjunto de bifurcacin se producir por medio de una analoga, es decir se
resolver la ecuacin de balance del buque por mtodos aproximados analticos
en las proximidades de la frecuencia de resonancia y se encontrar que la
amplitud de balance del movimiento estacionario satisface una ecuacin similar a
una de las dos superficies de equilibrio, segn los parmetros considerados en la
ecuacin. En estas condiciones se asimilar el comportamiento cualitativo del
buque en la zona de resonancia al de una de las dos catstrofes elementales
estudiadas.
La ecuacin de balance del buque una vez adimensionalizada (proceso recogido
en el Anexo I) es una ecuacin diferencial no lineal de segundo grado del tipo:
Ecuacin no autnoma (aparece explcitamente el tiempo t), con
La solucin de la ecuacin antes mencionada cuando "f" es sinusoidal de
frecuencia w w^ (es decir en la zona de resonancia) presenta un estado
estacionario peridico de frecuencia w.
Las tcnicas de medianizacin de Krylov y Bogoliouvov son aplicables a
sistemas de la forma X = xg(x,t); xeR" y x << 1 donde "g" es de periodo T en t. La
ecuacin (1.0) se puede convertir en una del tipo (2.9) por la transformacin de
Van der Pol y resolverse analticamente haciendo uso del teorema de
medianizacin.
1.0
13
La amplitud y fase de la solucin, funcin de los parmetros en juego, presentan
particularidades diferentes segn el nmero de alinealidades consideradas e
incluso segn el signo de alguna de ellas. El comportamiento de dichas soluciones
es anlogo al de alguna de las catstrofes elementales estudiadas por Thom.
ECUACIN DE BALANCE NO LINEAL
La ecuacin que rige el movimiento de balance independiente es una ecuacin
diferencial de segundo orden que se puede escribir:
Donde, es el ngulo de balance del buque con respecto a la superficie en reposo
del mar I.d ^/dt es el momento inercial e I representa la suma del momento de
inercia del buque ms el momento de inercia de la masa aadida que, en primera
aproximacin, se puede suponer invariante con el ngulo de escora, D(d/dt) es el
momento de las fuerzas amortiguadoras o resistentes que ser funcin de la
velocidad del movimiento, R() es el momento de las fuerzas restauradoras o
adrizantes, funcin del ngulo de balance y, por ltimo, M(t) ser el momento de
las fuerzas excitadoras, dependiente del tiempo.
El hecho de elegir un balance independiente, indica que no se tendrn en cuenta
posibles acoplamientos con los otros modos de movimiento, es decir, el
movimiento es balance puro. Esta simplificacin es solamente posible en el caso
de considerar mares transversales al buque, aquellos cuya direccin de
propagacin de las olas sea perpendicular a la direccin de la eslora del buque.
Cualquier otra direccin de propagacin del oleaje provocar otros movimientos
simultneos. No obstante, existe la posibilidad de aproximar esta ecuacin para
ngulos de incidencia de las olas prximos a 90 grados.
La expresin ms general del momento de las fuerzas amortiguadoras D(d/dt) en
la ecuacin (3.1) deber incluir los efectos de absorcin de energa producida por
las fuerzas de friccin en el casco del buque, generacin de olas en el movimiento
1.2
14
y formacin de vrtices, aparte de otros fenmenos de menor importancia como
calor generado y tensin superficial. Una buena aproximacin ser:
D(d/dt) = kd/dt + k1 d/dt.d/dt
donde el trmino k1 d/dt aparece asociado con formacin de olas y
amortiguamiento viscoso y kd/dt representa la resistencia friccional y
Torbellinos. Dalzell encontr que el modelo de Amortiguamiento dado por la
ecuacin (3.2) se puede aproximar con un elevado grado de fiabilidad mediante la
expresin:
ms sencilla desde el punto de vista del tratamiento matemtico.
Cuando en el estudio del movimiento de balance se considera un momento
restaurador lineal R() = .GM. , se est dando el mismo tratamiento a buques
cuyas caractersticas restauradoras pueden ser sensiblemente diferentes, este
serla R()= .GM. el caso de los dos buques cuyas curvas de estabilidad
esttica se representan en la figura 3.1 (a) y (b). En el caso (a), en un tratamiento
lineal se estarla subvalorando su capacidad de respuesta en una banda amplia de
ngulo de escora, mientras en (b) ocurrira todo lo contrario. La aproximacin
lineal solamente sera vlida en el tramo recto de la curva de estabilidad esttica,
es decir, para ngulos menores de, aproximadamente, 8. De una forma ms
general el momento de las fuerzas restauradoras R() se debe aproximar por una
expresin polinomial de potencias impares (con objeto de conservar la simetra
respecto al origen) del ngulo de balance multiplicada por el desplazamiento del
buque:
15
donde los coeficientes c^ se deben obtener de un ajuste adecuado de la curva de
estabilidad esttica del buque en la situacin de carga estudiada. En detrimento de
lo que podra ser un mejor ajuste de la curva, se considera adecuado tomar ci =
GM, altura metacntrica transversal, con lo que se conserva la pendiente en el
origen de la curva, al mismo tiempo que no se introducen parmetros de escaso
significado en la ingeniera naval.
El momento de excitacin en la hiptesis de Froude de longitud de olas grande en
comparacin con la manga del buque, vendr dado por:
16
17
que es la que se asume rige el movimiento de balance del buque. La ecuacin
(3.6) se puede adimensionalizar para una mayor facilidad de clculo, para ello se
introducen como parmetros adimensionales
donde representa el ngulo de escora que presenta el brazo mximo en la
curva de estabilidad esttica, de esta manera para valores de x < 1 el buque
estar oscilando con ngulos de escora con suficiente reserva de estabilidad y
para valores de x > 1 se encontrar en la zona en que el brazo de palanca de
adrizado decrece al aumentar el ngulo, con el peligro que ello conlleva.
Por otro lado se escoge la escala de tiempos adecuada para que la frecuencia
propia "corregida" w sea igual a la unidad y la nueva frecuencia adimensional w
coincida con el primitivo factor de sintonizacin /.
La ecuacin de balance obtenida quedar:
18
Conclusin
Hemos visto varios mtodos de resolucin de problemas matemticos utilizando
ecuaciones diferenciales lineales y no lineales, como tambin las aplicaciones en
diversas ramas de la ciencia que por otros mtodos resultara sumamente tedioso.
Ya que lo que se busca, ms que nada es la eficiencia de un sistema determinado,
esta puede ser lograda mediante la Matemtica Aplica, ahorrando tiempo valioso
en los procesos y desarrollo de clculos; y obteniendo precisin en los valores
finales de estos, as, reafirmando la importancia del conocimiento de estos
mtodos para poder estar siempre a la cabeza en un mercado competitivo.