1
ÍNDICE
INTRODUCCIÓN .................................................................................................................................... 2
MÉTODO DE EULER ............................................................................................................................. 3
APLICACIONES EN EL CAMPO DE LA INGENIERÍA ................................................................ 4
INGENIERÍA CIVIL: DEFLEXIÓN DE LA VIGA .......................................................................... 4
EJERCICIO DEL MÉTODO DE EULER ......................................................................................... 6
MÉTODO DE RUNGE-KUTTA ............................................................................................................. 7
APLICACIONES DEL MÉTODO DE RUNGE-KUTTA EN EL CAMPO DE LA INGENIERÍA 8
INGENIERÍA MECÁNICA: ANÁLISIS DE VIBRACIONES ......................................................... 8
EJERCICIO DEL MÉTODO DE RUNGE-KUTTA (ORDEN 4) ................................................... 10
CONCLUSIONES .................................................................................................................................. 12
BIBLIOGRAFÍA .................................................................................................................................... 12
2
INTRODUCCIÓN
Este trabajo se refiere a los siguientes temas: “Método de Euler” y “Método de Runge-Kutta”.
Las ecuaciones diferenciales surgen de manera natural al momento del movimiento, situaciones físicas
en los campos de las ciencias naturales, la ingeniería y otros campos de estudio, donde existen envueltas
causas de cambio de una o varias funciones no identificadas en relación a una de una o múltiples
variables autónomas, se los puede resolver mediante el método de Euler y también con el método de
Runge-Kutta, a continuación, se detallará de que trata cada uno de los métodos. El método de Euler se
refiere al proceso iterativo que proporciona una aproximación, se trata de solucionar una ecuación
diferencial de primer orden, con un valor inicial establecido, o sea, encontramos una respuesta numérica
(Cabrera Gómez, y otros, 2024).
Los métodos de Runge-Kutta son vistos como métodos de un solo paso, posiblemente uno de
los métodos numéricos más empleados y uno de los más exactos, empleados para conseguir soluciones
aproximadas a un problema con valor inicial (Segarra Escandón, 2020). Es crucial señalar que los
métodos de Runge- Kutta son alargamientos de la fórmula fundamental de Euler. Por lo tanto, se afirma
que el método de Euler es un método de primer orden de Runge-Kutta.
El método de Euler se aplica en varias ingenierías, pero es más utilizada es en la ingeniería
civil, por ejemplo, para la flexión en vigas homogéneas. El método de Runge-Kutta se lo puede aplicar
en la ingeniería mecánica para los sistemas mecánicos, ingeniería eléctrica para los circuitos eléctricos.
El objetivo de esta investigación es aportar información necesaria para motivar a que se investiguen más
ejemplos de aplicaciones de los dos métodos y resuelvan sus respectivas ecuaciones para comprobar sus
resultados y así aplicarlos en cualquier ingeniería mencionadas anteriormente.
3
MÉTODO DE EULER
Euler en el año 1760 fue el primero en sugerir un enfoque numérico que permitiera la solución
de la ecuación diferencial. “El método de Euler o de las rectas tangentes es una de las técnicas más
simples para aproximar soluciones de una ecuación diferencial” (Córdoba Quiroga, 2004). Considerando
una ecuación diferencial con un valor inicial, se puede determinar un valor inicial en la ecuación
diferencial 𝑦̇ = (𝑡, 𝑦), 𝑦(𝑡0 ) = 𝑦0 , el método de Euler se fundamenta en un concepto tan sencillo como
segmentar el lapso de tiempo para encontrar la solución en subintervalos [𝑡𝑘−1 , 𝑡𝑘 ] de longitud h. El
método de Euler nos indica que el valor de la función al término de un subintervalo es el mismo que al
final de un subintervalo [𝑡𝑘−1 , 𝑡𝑘 ] se define por la siguiente expresión.
𝑦𝑘 = 𝑦𝑘−1 + ℎ𝑓 (𝑡𝑘−1 , 𝑦𝑘−1 ) = 𝑦𝑘−1 + ℎ𝑦̇ (𝑡𝑘−1 , 𝑦𝑘−1 )
Este es un procedimiento explícito, pues emplea valores del paso previo para progresar en un
paso, en el que se determina el punto final 𝑦𝑘 del intervalo [𝑡𝑘−1 , 𝑡𝑘 ] moviéndonos sobre la recta
tangente en el punto (𝑡𝑘−1 , 𝑦𝑘−1 ), que es el punto inicial del subintervalo.
Una de las modificaciones realizadas al método de Euler es el denominado método de Euler
implícito o regresivo, en el que se emplea la derivada del punto final del subintervalo para progresar un
paso. La fórmula es la siguiente: 𝑦𝑘 = 𝑦𝑘−1 + ℎ𝑓 (𝑡𝑘 , 𝑦𝑘 ) = 𝑦𝑘−1 + ℎ𝑦̇ (𝑡𝑘 , 𝑦𝑘 ). Este método consiste em
calcular el punto final (𝑡𝑘 , 𝑦𝑘 ), de tal forma que la recta tangente en el punto (𝑡𝑘 , 𝑦𝑘 ) contenga al punto
inicial del intervalo (𝑡𝑘−1 , 𝑦𝑘−1 ).
4
APLICACIONES EN EL CAMPO DE LA INGENIERÍA
INGENIERÍA CIVIL: DEFLEXIÓN DE LA VIGA
Según (Morales, Hernández, & Gutiérrez, 2015) , “las deflexiones son las alteraciones que
experimenta un elemento debido al impacto de las flexiones internas”. Para calcular la deflexión, se
utilizan las leyes que vinculan fuerzas y desplazamientos mediante dos clases de métodos de cálculo: los
geométricos y los de energía. Según la figura 1, considere una viga horizontal AB, se asume que la viga
tiene una sección transversal uniforme y un material uniforme. El eje de simetría se sitúa en el plano
intermedio, señalando la región sombreada.
Figura 1. Viga sin flectar en su eje de simetría. (Morales, Hernández, & Gutiérrez, 2015)
Cuando la viga está expuesta a fuerzas, asumiendo que se encuentran en un plano que incluye el eje de
simetría, su forma puede verse alterada debido a su elasticidad, tal como se ilustra en la figura 2.
Figura 2. Viga afectada en su eje de simetría (Morales, Hernández, & Gutiérrez, 2015)
Para hallar la deflexión de la viga hay que determinar la ecuación de la curva elástica por lo
cual para poder formular la ecuación debemos saber:
Sea 𝑀(𝑥 ) el momento flector en una sección vertical transversal de la viga en 𝑥. Este momento
flector se establece como la adición algebraica de los momentos de las fuerzas que inciden en un lado de
𝑥, estos momentos se registran en una línea horizontal en la sección transversal de 𝑥. Al determinar los
momentos, seguiremos la regla de que las fuerzas ascendentes generan momentos negativos y las
ascendentes generan momentos positivos, asumiendo que el eje se inclina hacia abajo.
5
No tiene importancia qué lado de 𝑥 se seleccione, ya que los momentos flectores determinados desde
cualquier lado se mantienen constantes.
𝑌̈
𝐸𝐼 = 𝑀(𝑥)
[1 + (𝑦̇ )2 ]32
Donde E representa el módulo de elasticidad de Young y está vinculado al material utilizado en
el diseño de la viga, mientras que I representa el momento de inercia de la sección transversal de la viga
en 𝑥 en relación a una línea horizontal que atraviesa al centro de gravedad de dicha sección transversal.
𝐸𝐿 se llama rigidez y se considerará como una constante.
Si consideramos que la viga solo se dobla ligeramente, lo que es aplicable a numerosos usos
prácticos, la pendiente 𝑦̇ de la curva elástica es tan baja que su cuadrado es insignificante en
comparación con 1, y la ecuación puede ser sustituida por la aproximación apropiada:
𝐸𝐿𝑦̇ = 𝑀(𝑥)
El momento flector 𝑀(𝑥) representa la adición algebraica de los momentos de dichas fuerzas
que operan a un lado del punto P.
Al seleccionar el lado derecho de P, actuarían dos fuerzas:
• La fuerza hacia abajo 𝑤(𝐿 − 𝑥), a una distancia (𝐿 − 𝑥)/2 de P, produciendo un
momento positivo.
• La fuerza hacia arriba 𝑤𝐿/2, a una distancia 𝐿 − 𝑥 de P, produciendo un momento
negativo. En este caso el momento flector es:
(𝐿 − 𝑥 ) 𝑤𝐿 𝑤𝑥 2 𝑤𝐿𝑥
𝑀 (𝑥 ) = 𝑤 (𝐿 − 𝑥 ) − − (𝐿 − 𝑥 ) = −
2 2 2 2
6
Con el valor de 𝑀 (𝑥 ), la ecuación fundamental es:
𝑤𝑥 2 𝑤𝐿𝑥
𝐸𝐿𝑦̈ = −
2 2
Se requieren dos condiciones para determinar 𝑦. Esto significa que 𝑦 = 0 cuando 𝑥 = 0 y en
𝑥 = 𝐿, ya que la viga no presenta deformación en sus extremos o soportes. Como la ecuación
requerida de la curva elástica. Es de interés práctico usar la solución final.
𝑤
𝑦= (𝑥 4 − 2𝐿𝑥 3 + 𝐿3 𝑥)
24𝐸𝐼
Para determinar la deflexión máxima. En términos de simetría o cálculo, el máximo se produce
en 𝑥 = 𝐿/2, donde la flecha máxima será:
𝟓𝒘𝑳𝟒
𝒚𝒎𝒂𝒙 =
𝟑𝟖𝟒𝑬𝑰
Figura 3. Curva de deflexiones para el ejemplo. (Morales, Hernández, & Gutiérrez, 2015)
EJERCICIO DEL MÉTODO DE EULER
Usar el método de Euler para aproximar la solución del P.V.I dado en los puntos x = 0.1, 0.2, 0.3, 0.4,
0.5 usando tamaño de paso h = 0.1.
𝑦𝑛+1 = 𝑦𝑛 + ℎ𝑓(𝑥𝑛 , 𝑦𝑛 )
𝑑𝑦 𝑥
= −𝑦
𝑑𝑥
𝑥0 = 0 𝑦0 = 4
0
𝑥1 = 0,1 𝑦1 = 𝑦0 + ℎ𝑓(𝑥0 , 𝑦0 ) = 4 + 0,1 ∙ (− 4) = 4
0,1
𝑥2 = 0,2 𝑦2 = 𝑦1 + ℎ𝑓(𝑥1 , 𝑦1 ) = 4 + 0,1 ∙ (− ) = 3,9975
4
0,2
𝑥3 = 0,3 𝑦3 = 𝑦2 + ℎ𝑓(𝑥2 , 𝑦2 ) = 3,9975 + 0,1 ∙ (− ) = 3,9925
3,9975
7
0,3
𝑥4 = 0,4 𝑦4 = 𝑦3 + ℎ𝑓(𝑥3 , 𝑦3 ) = 3,9925 + 0,1 ∙ (− ) = 3,2411
3,9925
0,4
𝑥5 = 0,5 𝑦5 = 𝑦4 + ℎ𝑓(𝑥4 , 𝑦4 ) = 3,2411 + 0,1 ∙ (− ) = 3,2288
3,2411
MÉTODO DE RUNGE-KUTTA
Para (Valdillo, 2023), “Los métodos de Runge-Kutta son la clase más importante de los
métodos de un paso para aproximar e.d.o”. El método de Runge-Kutta aproxima la solución de las
series de Taylor, no obstante, en contraposición a la solución formal de las series de Taylor, el método
R-K no demanda una valoración explícita de derivadas superiores a la primera. Algunas evaluaciones
de las primeras derivadas consideran los efectos de las grandes derivadas. Según la cantidad de
términos retenidos efectivamente en las series de Taylor, existen métodos R-K de distinto orden.
Método R-K de segundo orden:
𝑘1 + 𝑘2
𝑥1 = 𝑥0 + ∆𝑥 = 𝑥0 +
2
Donde
𝑘1 = 𝑓(𝑥0 , 𝑡0 )∆𝑡
𝑘2 = 𝑓(𝑥0 + 𝑘1 , 𝑡0 + ∆𝑡)∆𝑡
Método R-K de orden tres:
ℎ
𝑦𝑛+1 = 𝑦𝑛 + (𝑘1 + 4𝑘2 + 𝑘3 ),
6
𝑘1 = 𝑓(𝑡𝑛 , 𝑦𝑛 ),
1 1
𝑘2 = 𝑓 (𝑡𝑛 + ℎ, 𝑦𝑛 + ℎ𝑘1 ),
2 2
𝑘3 = 𝑓(𝑡𝑛 , +ℎ, 𝑦𝑛 − ℎ𝑘1 + 2ℎ𝑘2 ).
Método R-K de orden cuatro:
ℎ
𝑦𝑛+1 = 𝑦𝑛 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 ),
6
𝑘1 = 𝑓(𝑡𝑛 , 𝑦𝑛 ),
8
1 1
𝑘2 = 𝑓 (𝑡𝑛 + ℎ, 𝑦𝑛 + ℎ𝑘1 ),
2 2
1 1
𝑘3 = 𝑓 (𝑡𝑛 , + ℎ, 𝑦𝑛 + ℎ𝑘2 ).
2 2
𝑘4 = 𝑓(𝑡𝑛 + ℎ, 𝑦𝑛 + ℎ𝑘3 )
Método R-K de quinto orden de Butcher (1964):
ℎ
𝑦𝑛+1 = 𝑦𝑛 + (7𝑘1 + 32𝑘3 + 12𝑘4 + 32𝑘5 + 7𝑘6 ),
90
Donde:
𝑘1 = ℎ𝑓(𝑥𝑛 , 𝑦𝑛 )
1 1
𝑘2 = 𝑓 (𝑥𝑛 + ℎ, 𝑦𝑛 + ℎ𝑘1 ),
4 4
1 1 1
𝑘3 = ℎ𝑓 (𝑥𝑛 , + ℎ, 𝑦𝑛 + ℎ𝑘1 + ℎ𝑘2 ).
4 8 8
ℎ 1
𝑘4 = 𝑓(𝑥𝑛 + , 𝑦𝑛 − 𝑘1 ℎ + 𝑘3 ℎ)
2 2
3 3 9
𝑘5 = 𝑓 (𝑥𝑛 , + ℎ, 𝑦𝑛 + ℎ𝑘1 + ℎ𝑘4 )
4 16 16
3 2 12 12 8
𝑘6 = 𝑓 (𝑥𝑛 , +ℎ, 𝑦𝑛 − ℎ𝑘1 + ℎ𝑘2 + 𝑘3 ℎ − 𝑘4 ℎ + 𝑘5 ℎ)
7 7 7 7 7
APLICACIONES DEL MÉTODO DE RUNGE-KUTTA EN EL CAMPO DE LA INGENIERÍA
INGENIERÍA MECÁNICA: ANÁLISIS DE VIBRACIONES
El método de Runge-Kutta es un recurso numérico inestimable para el estudio de vibraciones, en
particular cuando se trata de ecuaciones diferenciales ordinarias (EDO) que representan el
comportamiento dinámico de sistemas vibratorios. Estas EDO, frecuentemente no lineales o con
coeficientes variables, son complicadas de solucionar de manera analítica, lo que hace indispensable la
aproximación numérica.
9
Se considera un sistema de masa-resorte-amortiguador tradicional, que representa uno de los
esquemas de vibración mecánica más sencillos. El movimiento de este sistema es descrito por la ecuación
diferencial siguiente:
𝑚𝑥´´(𝑡) + 𝑐𝑥´´(𝑡) + 𝑘𝑥 (𝑡) = 𝐹(𝑡)
Donde:
• m: masa del sistema.
• c: coeficiente de amortiguamiento.
• k: constante de rigidez del resorte.
• x(t): desplazamiento en función del tiempo.
• x´(t): velocidad en función del tiempo.
• x´´(t): aceleración en función del tiempo.
• F(t): fuerza externa aplicada
Para implementar el procedimiento de Runge-Kutta, es habitual convertir la ecuación de
segundo orden en un sistema de dos ecuaciones de primer orden. Se estipula:
𝑦1 ´ = 𝑦2
(−𝑐𝑦2 − 𝑘𝑦1 + 𝐹 (𝑡))
𝑦2 ´ =
𝑚
Supongamos que deseamos aplicar el método de cuarto orden de Runge-Kutta. Las fórmulas
comunes para este procedimiento son:
ℎ
𝑦𝑛+1 = 𝑦𝑛 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 ),
6
𝑘1 = 𝑓(𝑡𝑛 , 𝑦𝑛 ),
10
1 1
𝑘2 = 𝑓 (𝑡𝑛 + ℎ, 𝑦𝑛 + ℎ𝑘1 ),
2 2
1 1
𝑘3 = 𝑓 (𝑡𝑛 , + ℎ, 𝑦𝑛 + ℎ𝑘2 ).
2 2
𝑘4 = 𝑓(𝑡𝑛 + ℎ, 𝑦𝑛 + ℎ𝑘3 )
Donde:
• h: tamaño del paso de integración.
• 𝑡𝑛 : tiempo en el paso n.
• 𝑦𝑛 : valor de la función en el paso n.
• 𝑓 : función que define el sistema de ecuaciones diferenciales.
Aplicándole la fórmula al ejemplo la función sería:
𝑓 (𝑡, 𝑦) = [𝑦2 ; (−𝑐𝑦2 − 𝑘𝑦1 + 𝐹 (𝑡))/𝑚]
Al reemplazar en las fórmulas de Runge-Kutta, obtendríamos las ecuaciones para determinar
𝑦1(desplazamiento) y 𝑦2 (velocidad) en cada periodo de tiempo.
EJERCICIO DEL MÉTODO DE RUNGE-KUTTA (ORDEN 4)
𝑑𝑦
Dada la ecuación diferencial = 𝑥 √𝑦 𝑦(1) = 4, aplique el método de R-K de cuarto orden para
𝑑𝑥
aproximar 𝑦(1.6) usando ℎ = 0.2
Primera iteración (𝒊 = 𝟎)
𝑥0 = 1, 𝑦0 = 4
𝑘1 = 𝑓(𝑥0 , 𝑦0 ) = 𝑥0 √𝑦0 = 1√4 = 2
1 1
𝑘2 = 𝑓 (𝑥0 + ℎ, 𝑦0 + ℎ𝑘1 ) = 𝑓(1.1, 4.2) = 1.1√4.2 = 2.25433
2 2
1 1
𝑘3 = 𝑓 (𝑥0 , + ℎ, 𝑦0 + ℎ𝑘2 ) = 𝑓(1.1, 4.22543) = 1.1√4.22543 = 2.26114
2 2
𝑘4 = 𝑓(𝑥0 + ℎ, 𝑦0 + ℎ𝑘3 ) = 𝑓(1.2, 4.45222) = 1.2√4.45222 = 2.53204
ℎ
𝑦1 = 𝑦0 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 ) = 4.45210
6
11
Segunda iteración (𝒊 = 𝟏)
𝑥1 = 𝑥0 + ℎ = 1.2, 𝑦1 = 4.45210
𝑘1 = 𝑓(𝑥1 , 𝑦1 ) = 𝑥1 √𝑦1 = 1.2√4.45210 = 2.53200
1 1
𝑘2 = 𝑓 (𝑥1 + ℎ, 𝑦1 + ℎ𝑘1 ) = 𝑓(1.3, 4.7053) = 1.3√4.7053 = 2.81992
2 2
1 1
𝑘3 = 𝑓 (𝑥1 , + ℎ, 𝑦1 + ℎ𝑘2 ) = 𝑓(1.3, 4.73409) = 1.3√4.73409 = 2.82854
2 2
𝑘4 = 𝑓(𝑥1 + ℎ, 𝑦1 + ℎ𝑘3 ) = 𝑓(1.4, 5.01781) = 1.4√5.01781 = 3.13606
ℎ
𝑦2 = 𝑦1 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 ) = 5.01760
6
Tercera iteración (𝒊 = 𝟐)
𝑥2 = 𝑥1 + ℎ = 1.4, 𝑦2 = 5.01760
𝑘1 = 𝑓(𝑥2 , 𝑦2 ) = 𝑥2 √𝑦2 = 1.4√5.01760 = 3.13600
1 1
𝑘2 = 𝑓 (𝑥2 + ℎ, 𝑦2 + ℎ𝑘1 ) = 𝑓(1.5, 5.3312) = 1.5√5.3312 = 3.46341
2 2
1 1
𝑘3 = 𝑓 (𝑥2 , + ℎ, 𝑦2 + ℎ𝑘2 ) = 𝑓(1.5, 5.36394) = 1.5√5.36394 = 3.47403
2 2
𝑘4 = 𝑓(𝑥2 + ℎ, 𝑦2 + ℎ𝑘3 ) = 𝑓(1.6, 5.7124) = 1.6√5.7124 = 3.82410
ℎ
𝑦3 = 𝑦2 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 ) = 5.71210
6
Valores que se ha encontrado:
i 𝒙𝒊 𝒚𝒊 𝒚𝒊+𝟏
0 1 4 4.45210
1 1.2 4.45210 5.01760
2 1.4 5.01760 5.71210
12
CONCLUSIONES
El método de Euler nos permite hallar soluciones aproximadas a un problema de ecuaciones
diferenciales, aunque no es un método muy exacto debido a sus altos porcentajes de error, solo se emplea
para hallar estimaciones aproximadas. La solución de diversos problemas diarios, muchos de ellos de
ingeniería, suelen estar vinculados a resultados numéricos. Con todo lo analizado, se puede concluir que,
para solucionar estos problemas, el Método de Runge-Kutta es la opción más adecuada entre el Método
de Euler. Además, es muy sencillo incrementar la precisión si empleamos pasos más reducidos entre los
puntos, lo que proporciona más sencillez. A pesar de todo, Runge-Kutta continúa siendo el método más
empleado para ofrecer mínimos márgenes de error en relación a la solución auténtica del problema, con el
beneficio de que es más sencillo que desarrollar un software de este método para llevar a cabo las
iteraciones requeridas.
BIBLIOGRAFÍA
Cabrera Gómez, G., Gamboa Pérez, M., Mateo Navas, L. M., Monge Romojaro, M., Pérez Pérez, T., Pons
Bordería, M. J., . . . Taipe Hidalgo, D. P. (2024). Material audiovisual clasificado por asignaturas
o contenidos: recurso para la docencia de Matemáticas y Estadística. Obtenido de
https://docta.ucm.es/bitstreams/c464de7e-7f3b-498b-aaba-8e89e3ef3da9/download
Córdoba Quiroga, I. M. (2004). MÉTODOS NUMÉRICOS . Obtenido de
https://www.academia.edu/download/57871321/METODOS_NUMERICOS.pdf
Morales, D., Hernández, C., & Gutiérrez, J. (2015). Aplicación del Método de Euler de orden superior en
la Deflexión de Vigas Homogéneas . Obtenido de
https://www.academia.edu/download/50259599/Aplicacion_de_Metodos_Numericos_de_Orden.
pdf
Segarra Escandón, J. (2020). MÉTODOS NUMÉRICOS RUNGE-KUTTA Y ADAMS BASHFORTH-
MOULTON EN MATHEMATICA. Obtenido de
https://ojs.urepublicana.edu.co/index.php/ingenieria/article/view/639/507
Valdillo, F. (2023). METODOS DE UN PASO. MÉTODOS RUNGE-KUTTA . Obtenido de
https://www.ehu.eus/~mepvaarf/ficheros/num_ode_un_paso.pdf
