Exercícios Resolvidos de Equações Diferenciais Ordinárias (EDO)

Ver Teoria

Enunciado

Uma avançada turbina chamada de ventilador não canalizado (UDF) é uma das novas tecnologias mais promissoras que tem sido desenvolvida para o futuro transporte de aeronaves.

Turbinas, que têm sido usadas por décadas, combinam o poder e a confiabilidade dos motores a jato com a eficiência dos propulsores. Eles constituem uma importante melhoria dos antigos propulsores movidos a pistão. Suas aplicações têm sido limitadas a pequenas aeronaves do tipo comutador, isto porque eles não são tão rápidos, nem poderosos quanto as turbinas usadas em grandes aeronaves. Esse tipo de turbina (UDF) implica em avanços significantes na tecnologia de propulsão. Novos materiais, aerodinâmica e velocidades de alta rotação habilitam esta turbina a voar tão rápido quanto as turbinas a jato, e com grande aproveitamento de combustível. A UDF é também menos barulhenta que o sistema convencional de turbinas.

Durante um teste de vôo de uma turbina UDF de uma aeronave, o motor é levado para um nível de 40   N , o que significa os 20   k g da aeronave tendo alcançado uma velocidade de 180   m / s . As válvulas de regulação do motor são então levadas para atingir um nível de 60   N , e a aeronave começa a acelerar.

A equação diferencial que determina a aceleração da aeronave é:

a = T m - 0,000062 v 2

Onde:

a = d v d t

T   =   nível atingido em N

m   =   massa em k g

v   =   velocidade em m / s

Escreva no MATLAB um programa para determinar a nova velocidade depois de uma mudança no nível do motor através do gráfico da solução para a equação diferencial. Considere um intervalo de tempo de 6 min.

Passo 1

O enunciado nos dá a seguinte EDO de primeira ordem:

a = d v d t = v ' = T m - 0,000062 v 2

Substituindo m por 20   k g , que é o valor da massa da aeronave, temos a EDO:

v ' = g ( t , v ) = T 20 - 0,000062 v 2

Onde T é o nível do motor em N , o tempo t em s é a variável independente e v ( t ) , a velocidade da aeronave em m / s , é a função que queremos descobrir resolvendo a EDO.

O enunciado diz que quando o nível do motor é T =   40   N , a aeronave atinge uma velocidade de 180   m / s . Em seguida, o nível do motor sobe para T =   60   N . Ao final, nos é perguntado qual é a velocidade da aeronave com o motor atinge um nível de T =   60   N . Para encontrar esta velocidade, temos que considerar T = 60 na EDO acima. Logo, a EDO que temos que encontrar uam solução é:

v ' = g ( t , v ) = 60 20 - 0,000062 v 2

De maneira indireta, o enunciado nos dá o intervalo de integração e a condição inicial para resolver a EDO.

O enunciado nos diz pra considerarmos um intervalo de 6 min. Considerando o instante inicial de t 0 = 0 e o final t f = 6 * 60 = 360   s , temos que o intervalo de integração é t i n t = t 0 t f = 0 360 .

O enunciado diz que a aeronave tem uma velocidade de 180   m / s quando o nível do motor vai para T = 60   N . Logo, a condição inicial é v t 0 = v 0 = 18 0 .

Portanto, temos todas as entradas necessárias para usar a função o d e 45 do MATLAB.

Além disso, o exercício pede para mostrarmos no gráfico da solução da EDO a velocidade que a aeronave atinge quando o nível do motor é T =   60   N . Por isso, no final, também teremos que gerar um gráfico de v ( t ) considerando T = 60   N .

Passo 2

Primeiramente, vamos relembrar a função o d e 45 , para vermos o que entraremos no MATLAB. A função o d e 45 possui a seguinte estrutura:

[ t , v ]   =   o d e 45 ( E D O _ f u n , t _ i n t , v 0 )

Ela exige três entradas, fornecidas pelo enunciado:

  • E D O _ f u n : uma função representando a EDO de primeira ordem:
  • v ' = 60 20 - 0,000062 v 2

  • t _ i n t : intervalo de integração, normalmente no formato [ t 0   t f ] : o intervalo começa em t 0 e termina em t f . Neste exercício, t 0 = 0 e t f = 3 60 , logo, t _ i n t = 0 3 60 ;
  • v 0 : condição inicial de v ( t ) em t 0 = 0 . Neste exercício, v 0 = 18 0 .

E ela nos dá duas saídas: um vetor com valores de t (começando em t 0 = 0 e terminando em t f = 3 60 , correspondendo ao nosso intervalo de integração) e outro com os valores de v (iniciando em v 0 = 18 0 e terminando no valor de v correspondente a t f ), representando a solução da EDO. Em outras palavras, 𝑜𝑑𝑒45 nos retorna os pontos que representam a função v = f ( t ) para o intervalo de integração t _ i n t . Usaremos os vetores t e v para gerar o gráfico pedido.

Depois de relembrar a função o d e 45 , vamos abrir um script novo no MATLAB para criarmos uma função com a EDO do enunciado. Nesse script, digite as linhas de código:

Após criar a função E D O _ e x 3 , abra um novo script para implementar a função o d e 45 . Neste script, digite as linhas de código:

Na primeira linha, entramos com o intervalo de integração t _ i n t = 0 360 .

Na segunda linha, entramos com o valor inicial v 0 = 18 0 .

Na terceira linha, entramos com a função o d e 45 para encontrar uma solução da EDO. Repare que, no lugar de 𝐸𝐷𝑂_𝑓𝑢𝑛, coloquei um string com o nome da função que representa a EDO: ' E D O _ e x 3 ’. Não esqueça de pôr o nome entre aspas simples, como no print acima.

A solução da EDO é dada pelos vetores t e v , saídas de o d e 45 .

Passo 3

Com os códigos do passo 2, encontramos a solução da EDO. Mas falta ainda uma coisinha: gerar o gráfico de v = f ( t ) . Para isto, usaremos a saída de o d e 45 os vetores t e v .

No mesmo script onde você implementou o d e 45 , digite as linhas de código:

A primeira linha gera o gráfico da velocidade em função do tempo.

A segunda linha dá um nome para o eixo horizontal do gráfico, usando a função x l a b e l . Foi dado um nome de ‘ T e m p o   [ s ] ’ .

A terceira linha dá um nome para o eixo vertical do gráfico, usando a função y l a b e l . Foi dado um nome de ‘ V e l o c i d a d e   [ m / s ] ’ .

Ao rodar o programa, o gráfico de tempo versus velocidade é gerado e nele podemos encontrar a resposta do exercício:

Repare que a curva se estabiliza em torno de 220   m / s . Isto quer dizer que a nova velocidade da aeronave depois de uma mudança no nível do motor de 40   N para 60   N é 220   m / s .

Resposta

A nova velocidade da aeronave depois de uma mudança no nível do motor de 40   N para 60   N é 220   m / s .