0313320 - Métodos Numéricos Aplicados à Engenharia Civil

MÉTODO DE
RUNGE-KUTTA
Carolina Mendes Esposito - N°USP: 11734660
Felipe Bisognin - N°USP: 11821556
Paola Jossara Domingos - N°USP: 11858071
Thiago Fugii Moraes - NºUSP: 10774312
Sumário
1. Histórico
2. Explicação do Método
3. Exemplificação do Método
4. Método de Runge-Kutta vs Método de Euler
5. Aplicação Computacional
6. Aplicações na Engenharia
7. Tarefa para Casa
8. Referências Bibliográficas
Histórico
O método de Runge-Kutta é uma técnica numérica amplamente utilizada para resolver equações diferenciais
ordinárias (EDOs). Desenvolvido por Carl Runge e Martin Kutta no início do século XX, este método se tornou
uma ferramenta fundamental em diversas áreas da ciência e engenharia.

O contexto histórico que levou ao desenvolvimento do método de Runge-Kutta está intrinsecamente ligado à
necessidade de resolver equações diferenciais de forma eficiente. No final do século XIX e início do século XX,
os cientistas e engenheiros enfrentavam desafios significativos ao tentar modelar e compreender fenômenos
complexos que eram regidos por EDOs. Na ausência de soluções analíticas para muitos desses problemas, a
busca por métodos numéricos precisos e eficientes tornou-se imperativa.

O método de Runge-Kutta é uma abordagem iterativa que calcula a solução de uma EDO em incrementos
sucessivos. Sua popularidade cresceu rapidamente devido à sua precisão e eficiência em comparação com
outros métodos numéricos da época. Ao longo dos anos, vários aprimoramentos e variações do método foram
propostos, resultando em diferentes formas e ordens do método de Runge-Kutta.
Histórico
Carl Runge Martin Kutta
Carl Runge (1871–1946) foi um matemático Martin Kutta (1867–1944), também matemático
alemão cujas contribuições significativas para a alemão, colaborou com Runge no desenvolvi-
teoria numérica o destacam como um pioneiro mento do método de Runge-Kutta. Nascido em
no campo. Nascido em Bremen, Alemanha, Schlawe, Prússia (atualmente Polônia), Kutta
Runge dedicou-se ao estudo de equações refinou e generalizou os métodos propostos por
diferenciais e métodos numéricos. Em 1895, Runge em 1901. Sua contribuição essencial foi a
publicou um trabalho descrevendo um método formulação de um esquema mais abrangente e a
iterativo para a resolução de equações introdução de uma notação sistemática que
diferenciais ordinárias (EDOs), pavimentando o facilitou a compreensão e a implementação do
caminho para o desenvolvimento do que agora é método. A colaboração entre Runge e Kutta
conhecido como o método de Runge-Kutta. Seu resultou em uma abordagem numérica robusta e
trabalho estabeleceu as bases para abordagens versátil, que se tornou uma ferramenta funda-
numéricas que se tornariam cruciais em diversas mental na resolução de EDOs e teve um impacto
disciplinas científicas e de engenharia. duradouro no campo da análise numérica.
Explicação do Método
Entre os métodos empregados na solução de
problemas de valor inicial, os Métodos de
Runge-Kutta destacam-se pela precisão e
simplicidade.

Cada método de Runge-Kutta utiliza a

comparação com um polinômio de Taylor
apropriado para evitar o cálculo direto de
derivadas.
Explicação do Método
Método geral de Runge-Kutta

As constantes são obtidas na comparação

a expansão da função com a função
Explicação do Método
Métodos de Runge-Kutta de ordem 2
Desenvolvendo k2 em série de Taylor em torno
do ponto (x,y), obtemos:
Explicação do Método
Métodos de Runge-Kutta de ordem 2
Explicação do Método
Métodos de Runge-Kutta de ordem 2

Método de Euler modificado Método de Euler melhorado

Explicação do Método
Métodos de Runge-Kutta de ordem 3
Sistema de 4 equações e 6 incógnitas

Desde que

k1 e k2 possuem as mesmas expressões do

método de ordem 2
Explicação do Método
Métodos de Runge-Kutta de ordem 3

Método de Nystrom
Explicação do Método
Métodos de Runge-Kutta de ordem 4

Os dois métodos de Runge-Kutta de ordem 4 mais utilizados

Exemplificação do Método
Exemplo: Para a equação diferencial y’ = -x + y + 2 e condição inicial y(0) = 2, determine y(0,2). Considere
h = 0,1 e faça arredondamento com 4 casas decimais

Exercício para Casa: resolver a EDO do

f(x,y) = -x + y + 2
enunciado pelo método de Runge-Kutta de
x0 = 0
ordem 2 (utilizando apenas K1 e K2), com-
y0 = 2
parando com o resultado analítico e gráfico.
h = 0,1 (espaçamento, ou passos)

n
0 0 2
1
2
Exemplificação do Método
Exemplo: Para a equação diferencial y’ = -x + y + 2 e condição inicial y(0) = 2, determine y(0,2). Considere
h = 0,1 e faça arredondamento com 4 casas decimais

f(x,y) = -x + y + 2
x0 = 0
y0 = 2
h = 0,1 (espaçamento, ou passos)

Primeiro passo: n = 0
n=0
xn = 0
n
yn = 2
K1 = f(0,2) = -0 + 2 + 2 = 4 0 0 2
K2 = f(0,05, 2,2) = 4,15
K3 = f(0,05, 2,2075) = 4,1575 1
K4 = f(0,1, 2,4157) = 4,3157
2
Exemplificação do Método
Exemplo: Para a equação diferencial y’ = -x + y + 2 e condição inicial y(0) = 2, determine y(0,2). Considere
h = 0,1 e faça arredondamento com 4 casas decimais

f(x,y) = -x + y + 2
x0 = 0
y0 = 2
h = 0,1 (espaçamento, ou passos)

Primeiro passo: n = 0
n=0
xn = 0
n
yn = 2
K1 = f(0,2) = -0 + 2 + 2 = 4 0 0 2 2,4155
K2 = f(0,05, 2,2) = 4,15
K3 = f(0,05, 2,2075) = 4,1575 1
K4 = f(0,1, 2,4157) = 4,3157
2
Exemplificação do Método
Exemplo: Para a equação diferencial y’ = -x + y + 2 e condição inicial y(0) = 2, determine y(0,2). Considere
h = 0,1 e faça arredondamento com 4 casas decimais

f(x,y) = -x + y + 2
x0 = 0
y0 = 2
h = 0,1 (espaçamento, ou passos)

Segundo passo: n = 1
n=1
xn = 0,1
n
yn = 2,4155
K1 = 4,3155 0 0 2 2,4155
K2 = 4,4813
K3 = 4,4896 1 0,1 2,4155
K4 = 4,6645
2
Exemplificação do Método
Exemplo: Para a equação diferencial y’ = -x + y + 2 e condição inicial y(0) = 2, determine y(0,2). Considere
h = 0,1 e faça arredondamento com 4 casas decimais

f(x,y) = -x + y + 2
x0 = 0
y0 = 2
h = 0,1 (espaçamento, ou passos)

Segundo passo: n = 1
n=1
xn = 0,1
n
yn = 2,4155
K1 = 4,3155 0 0 2 2,4155
K2 = 4,4813
K3 = 4,4896 1 0,1 2,4155 2,8642
K4 = 4,6645
2 0,2 2,8642
Exemplificação do Método
Exercício: Usando:

Resolver:

Exercício para Casa: resolver a EDO do

enunciado pelo método de Runge-Kutta
de ordem 4 (utilizando de K1 a K4),
comparando com o resultado analítico.
Exemplificação do Método
Exercício. Usando:

Solução:
Passo 1:

Passo 2: n = 0
Exemplificação do Método
Exercício. Usando:

Solução:
Passo 3: n = 1
Exemplificação do Método
Exercício. Usando:

Solução:
Passo 4: n = 2
Exemplificação do Método
Exercício. Usando:

Solução:

K Xk Yk Y(Xk)

0 0,000000 2,000000 2,000000

1 0,100000 2,005000 2,004837

2 0,200000 2,019025 2,018731

3 0,300000 2,041218 2,040818

Método Runge-Kutta vs Método de Euler
Precisão
O método de Euler é um método de primeira ordem, o que significa que a precisão da solução é proporcional ao tamanho
do passo utilizado. Quanto menor o passo, mais precisa será a solução, mas isso implica em mais cálculos. O método de
Runge-Kutta, por outro lado, pode ser implementado em diferentes ordens, como RK2, RK3, RK4, etc. O RK4 é
especialmente usado, por ser um método de quarta ordem - o que proporciona uma precisão muito maior em comparação
com o método de Euler.

Complexidade Computacional
O método de Euler é relativamente simples de implementar, exigindo apenas uma única avaliação da derivada por
iteração. O método de Runge-Kutta, especialmente em suas versões de ordem superior, envolve mais cálculos por
iteração, tornando-o mais complexo computacionalmente. No entanto, isso é muitas vezes compensado pela maior
precisão, permitindo que o método de Runge-Kutta use passos maiores sem sacrificar a precisão.

Estabilidade
Em alguns casos, o método de Euler pode apresentar instabilidade numérica, especialmente em sistemas de equações
diferenciais que envolvem fenômenos oscilatórios ou não-lineares fortes. O método de Runge-Kutta, devido à sua
precisão superior, tende a ser mais estável em uma variedade mais ampla de situações.
 Aplicação Computacional
Código em Python que resolve uma Parâmetros:
equação diferencial ordinária de f: Função que representa a EDO (dy/dt = f(t, y))
y0: Condição inicial
primeira ordem (EDO) utilizando o
t0: Tempo inicial
Método Runge-Kutta de tn: Tempo final
quarta ordem (RK4). h: Tamanho do passo

Dentro do loop, calcula quatro coeficientes ‘k1’,

‘k2’, ‘k3’ e ‘k4’ usando o método de Runge-Kutta.
Esses coeficientes são baseados em avaliações da
função diferencial ‘f’ em diferentes pontos ao longo A lista ‘resultado’ armazena os pontos da solução,
do intervalo de tempo. começando com a condição inicial ‘(t0, y0)’. Em
seguida, usa-se um loop ‘while’ que continua até
Usamos esses coeficientes para atualizar o valor que o tempo ‘t’ alcance o tempo final ‘tn’.
da solução ‘y’ usando a fórmula do método de
Runge-Kutta. Também atualizamos o tempo ‘t’. Os O exemplo de uso define uma função ‘exemplo_edo’
novos valores (t, y) são adicionados à lista que representa a EDO:
‘resultado’.
Em seguida, são definidos os parâmetros iniciais.

A função ‘runge_kutta’ é chamada para resolver a

EDO. A solução é então impressa para cada ponto
calculado.
Aplicações na Engenharia
Engenharia Civil Engenharia Elétrica
Análise de Estruturas Circuitos Elétricos
Na análise dinâmica de estruturas, como pontes, edifícios e túneis, Podem ser usados para simular a resposta transitória do sistema a
os métodos de Runge-Kutta podem ser usados. diferentes condições iniciais ou entradas

Dinâmica de Solos Comunicações

Em estudos de interação solo-estrutura e análise de assentamento Modelar a propagação de sinais em meios de transmissão e simular
de fundações a resposta dinâmica de sistemas de comunicação

Modelagem Hidráulica Redes Elétricas

Em simulações de escoamento de água em canais, rios e redes de Simular transientes e eventos transitórios, fornecendo insights
distribuição, os métodos de Runge-Kutta podem ser úteis para sobre o comportamento dinâmico da rede
resolver equações diferenciais que representam os fenômenos
hidráulicos, como transientes em sistemas de tubulação

Engenharia Química Engenharia Aeroespacial

Dinâmica de Fluidos Dinâmica de Voo
Modelagem da dinâmica de fluidos em reatores químicos, ajudando Resolver as equações diferenciais que descrevem a dinâmica de
a entender e otimizar a distribuição de temperatura voo de aeronaves e espaçonaves

Processo de Separação Aeroelasticidade

Podem ser empregados para simular o comportamento dinâmico e Modelar as interações entre as forças aerodinâmicas, as
otimizar as condições operacionais para obter os produtos deformações estruturais e os movimentos do veículo
desejados
Aplicações na Engenharia
Pesquisas Acadêmicas
1 Método de Runge-Kutta Aplicado à Deflexão de Viga - UNIJUI 2017 3 Implementação Numérica de Problemas de Viscoelasticidade - 2013

O método numérico, especificamente o de Runge-Kutta, Neste trabalho, foi apresentada uma teoria unificada que assegura
apresentou um erro percentual muito baixo (0,0669%), indicando uma ordem de convergência nominal p>=2 para métodos de
alta precisão em comparação com o método analítico (0,0668%). Runge-Kutta (RK) utilizados como integradores temporais em
O estudo sugere que, mesmo em casos mais complexos, como modelos viscoelásticos analisados pelo método dos elementos
vigas com geometria não uniforme, o método numérico pode finitos. A validade desse conceito foi demonstrada em um modelo
fornecer resultados precisos a um custo reduzido em comparação viscoelástico de grandes deformações de Hartmann, com foco
com softwares de simulação. especial nos métodos RK diagonalmente implícitos (DIRK). O ponto
central desta pesquisa destaca que, para garantir a ordem de
convergência nominal do método de integração, a aproximação da

2 Método RK Aplicado em Sistemas Mecânicos em Fase - IFPE 2021 deformação ao longo do tempo deve ser consistente com a ordem
de convergência do integrador temporal.
Este trabalho destaca a significativa aplicação dos métodos
computacionais na resolução de equações diferenciais que
descrevem sistemas mecânicos, utilizando o software Wolfram
Mathematica® 11. A abordagem se concentra na mecânica
hamiltoniana no espaço de fase e utiliza o método numérico de
Runge-Kutta para resolver numericamente sistemas como o
pêndulo simples e osciladores harmônicos unidimensional,
bidimensional e anarmônico. Os resultados numéricos são
comparados com soluções analíticas, evidenciando a eficácia do
método. Além das implicações físicas, ressalta-se a amplitude dos
métodos numéricos, destacando sua importância na engenharia e
em aplicações simulativas no cotidiano
Tarefa para Casa
Referências Bibliográficas
Blog de Engenharia Química, Santos SP. Métodos de Euler e Runge-Kutta: Exemplo de Aplicação. Disponível em:
https://www.engquimicasantossp.com.br/2016/03/metodos-de-euler-e-runge-kutta.html.

Departamento de Matemática da Universidade Federal de Santa Catarina. Métodos de Euler e Runge-Kutta. Disponível
em: http://www.mtm.ufsc.br/~fermin/runge.

E-Disciplinas - MAP3121 (2022). Notas de Aula: Métodos Numéricos para Equações Diferenciais Ordinárias. Disponível
em: https://edisciplinas.usp.br/course/view.php?id=97397.

Stumpf, Felipe Tempel. Implementação numérica de problemas de viscoelasticidade finita utilizando métodos de
Runge-Kutta de altas ordens e interpolação consistente entre as discretizações temporal e espacial. Disponível em:
https://lume.ufrgs.br/handle/10183/75757.

De Lima, Claudia Claudiane. Aplicação do Método de Runge-Kutta em Sistemas Mecânicos no Espaço de Fase.
Disponível em: https://releia.ifsertao-pe.edu.br/jspui.

Dambroz, Giovani Prates. Método de Runge-Kutta Aplicado à Deflexão de Viga. Disponível em:
file:///C:/Users/Thiago/Downloads/8213-Texto%20do%20artigo-35043-1-10-20170921%20(2).pdf.