Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Tema2 Ejercicios

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 3

Tema 2: Ecuaciones diferenciales ordinarias

—Ejercicios—

Ejercicio 1. Resuelve los siguientes PVI por medio de la transformada de Laplace. Representa gráficamente
la solución.
( (
y 0 + y = f (t) 1, 0 ≤ t < 1
a) donde f (t) =
y(0) = 0 0, caso contrario
( (
y 0 + y = f (t) 1, 1 ≤ t < 2
b) donde f (t) =
y(0) = 0 0, caso contrario
( (
y 00 + 4y = f (t) sin(t), 3π2 ≤t
c) 0
donde f (t) =
y(0) = 0, y (0) = 1 0, caso contrario
(
00 0
y + 4y + 3y = 1 − u(t − 1) + u(t − 2) − u(t − 4)
d)
y(0) = 0, y 0 (0) = 1

Ejercicio 2. Resuelve el PVI


(
y 0 − 2y = t sin πt
y(1) = 1
de dos formas diferentes: condición inicial arbitraria y traslación en tiempo de la ecuación/solución.
Ejercicio 3. Resuelve la EDO
y 00 − 2y 0 + 3y = e−t cos 2t
en los siguientes casos.
a) y(0) = 0, y 0 (0) = 0
b) y(0) = y0 , y 0 (0) = y1
c) y(0) = 0, y(2π) = 0
d) y(2π) = 0, y(2π) = 0 (Prueba también con una traslación en tiempo de la ecuación/solución.)
Ejercicio 4. Dada una viga horizontal de longitud L, se tiene que la dislocación transversal (deflexión o
deformación) que sufre ésta al soportar un peso viene determinada por la ecuación diferencial
d4 y
EI = w(x) ,
dx4
donde EI es la resistencia a la flexión (producto de la elasticidad del material y el segundo momento de
inercia), y w es la distribución del peso a lo largo de la viga. Suponiendo que L = 6 m, EI = 1 N/m2 , y
que se distribuye uniformemente una carga de 10 kg en el tercio central. Calcula la dislocación de la viga
en los siguientes casos.
a) La viga está empotrada en una pared en su extremo izquierdo y está libre en el derecho (y(0) =
y 0 (0) = y 00 (0) = y 000 (0) = 0).
b) Ambos extremos de la viga están empotrados en sendas paredes (y(0) = y 0 (0) = y(L) = y 0 (L) = 0).
c) Ambos extremos de la viga están únicamente soportados por sendas columnas o puntos de apoyo
(y(0) = y 00 (0) = y(L) = y 00 (L) = 0).
d) ¿Cuáles serían las condiciones en caso de que estuviese empotrada en el extremo izquierdo y soportada
en el derecho?

28 de diciembre de 2020
Ejercicio 5. Considera la ecuación
y 0 = t2 − y 2
en el intervalo [0, 5]. La utilizaremos para testear algunos métodos.
a) En primer lugar, resuelve esta ecuación utilizando métodos explícitos que ya tengas implementados
como Runge-Kutta IV y Crank-Nicolson, utilizando para ello las condiciones iniciales y(0) = −0.675 y
y(0) = −0.676. ode45 también te puede servir para este experimento. ¿Qué observas? ¿Qué deduces
de ello?
En adelante, se fija para la ecuación anterior la condición inicial y(0) = −0.675, queriéndose determinar el
valor de la solución en t = 5. Para ello, se presentan las siguientes tablas de Butcher.
0
1 1 1 1 1 1
3 3 4 4 0 2 2
3 2 1
I. 2
3 − 13 1 II. 4 4 4 III. 1 − 12 1
2
1 1 −1 1 1 1 1 1
2 2 2 2
1 3 3 1
8 8 8 8

b) Para cada tabla, escribe las ecuaciones correspondientes como métodos de Runge-Kutta. Clasifica
estos métodos como explícitos, diagonalmente implícitos o completamente implícitos.
c) Ejecuta 3 iteraciones completas del primer método, 2 del segundo y 1 del tercero, con paso de
integración h = 0.1. Calcula, cuando sea posible el error local y el global. Puedes apoyarte en
soluciones numéricas de métodos ya implementados con un paso relativamente más pequeño h =
0.001.
d) Implementa los métodos y determina el orden numérico de estos (4, 2 y 2, respectivamente), ana-
lizando el error global en t = 5. Puedes hacerlo frente a una solución considerada «exacta» o con
errores aproximados. Comprueba cuánto se reduce el error cuando se divide el paso por 2 (es decir,
se multiplica por 2 el número de iteraciones, N = 2, 4, 8, 16, . . .).
Ejercicio 6. Utilizando desarrollos de Taylor, calcula el error local del método de Crank-Nicolson. ¿Cuál es,
por tanto, el orden del método?
Ejercicio 7 ([1], 25.16). El oscilador armónico suele representarse como una masa puntual en el extremo
inferior de un muelle fijado por su extremo superior. La ecuación
d2 x dx
m 2
+c + kx = 0
dt dt
obtenida gracias a la ley de Hook, modela el desplazamiento vertical x (m) respecto a la posición de
equilibrio (x = 0) en función del tiempo t (s), siendo m (kg) la masa, c el coeficiente de amortiguamiento
(Ns/m), k la elasticidad (N/m) del muelle.
Utiliza un método numérico para simular posición y velocidad de este sistema cuando x(0) = 1 m, x0 (0) = 0
m/s, m = 20 kg, k = 20 N/m, en los siguientes casos:
a) Infra-amortiguado: c = 5 Ns/m (oscilamiento).
b) Amortiguado crítico: c = 40 Ns/m
c) Sobre-amortiguado: c = 200 Ns/m (sin oscilamiento).
Nota: el tipo amortiguamiento está relacionado con la forma de la solución, siendo el caso crítico c2 −4mk =
0.
Ejercicio 8 ([1], 25.26). Tres «puentistas» están atados entre sí. Cada cuerda-puentista se modela como un
oscilador con respecto a su punto de anclaje. Tomando en cuenta los equilibrios entre las fuerzas ejercidas

28 de diciembre de 2020
(ley de Hook) por los sistemas adyacentes, se define el siguiente sistema de ecuaciones diferenciales
d2 x1
m1 = m1 g − k1 (x1 − 0) + k2 (x2 − x1 )
dt2
d2 x2
m2 2 = m2 g − k2 (x2 − x1 ) + k3 (x3 − x2 )
dt
d2 x3
m3 2 = m1 g − k3 (x3 − x2 )
dt
donde mi es la masa del i-ésimo saltador; ki la constante de la cuerda, xi el desplazamiento (respecto a la
posición de equilibrio) y g = 9.81 m/s2 es la aceleración de la gravedad.
Simula el sistema suponiendo que las posiciones y velocidades iniciales son nulas para los tres saltadores
con lo siguientes datos: m1 = 60 kg, m2 = 70 kg, m3 = 80 kg, k1 = 50 N/m, k2 = 100 N/m, k3 = 50
N/m.

Referencias
[1] Steven C. Chapra y Raymond P. Canale. Métodos numéricos para ingenieros (7ª edición). McGraw-Hill, Education, 2015.

28 de diciembre de 2020

También podría gustarte