Exercícios Resolvidos de Equação da Onda

Ver Teoria

Enunciado

Dada a seguinte equação do da onda:

y t t = α 2 y x x ,     0 < x < L ,     t > 0 y 0 , t = 0 ,     y L , t = 0 y x , 0 = 0 ,     y t x , 0 = g x

  1. O método de seperação de variáveis supõe que a solução da equação da onda é y x , t = X x T t . Usando esse método mostre como se chega ao sistema abaixo onde λ é a constante de separação. Analise os casos λ = 0 , λ < 0 e λ > 0 .
  2. X ' ' + λ X = 0 X 0 = X L = 0 T ' ' + α 2 λ T = 0 T 0 = 0

  3. Substituindo os valores L = 5 , α = 3 e g x = x 2 escolha a extensão par ou ímpar de período 2 L que resolve o problema em (a) encontrando a solução u x , t . Especifique os coeficientes dessa solução SEM calcular a integral.

Passo 1

Vamos lá a mais uma questão, e essa daqui vai dar um trabalho...

Mas calma! Não se assustem! A questão em si não é difícil, ela só é bem trabalhosa... Isto acontece pois temos que desenvolver toda a solução da equação da onda, e no final ainda temos que indicar como calcularíamos os coeficientes da solução (ainda bem que não precisamos calcular, né)

Sendo assim, vamos começar pela letra a). Acompanhem o raciocínio com calma pra não se perderem! Assim, será só sucesso!!

A questão pede para que a gente trabalhe com a seguinte equação da onda:

y t t = α 2 y x x ,     0 < x < L ,     t > 0 y 0 , t = 0 ,     y L , t = 0 y x , 0 = 0 ,     y t x , 0 = g x

O legál é que ele já nos diz como deve ser o formato da solução da nossa EDP. Vamos assumir que ele será dado por y x , t = X x T t . Sendo assim, vamos considerar a solução como o produto de duas funções, uma de x e a outra de t .

Essa consideração vai ajudar e tanto a nossa vida... Mas bom, continuando, vamos trabalhar com o método em si. Inicialmente, precisamos calcular as derivadas parciais da nossa solução u x , t . Isto é importante para que a gente consiga substituir na nossa EDP e fazer a análise.

Sendo assim, calculando as derivadas parciais (se formos derivar em x , consideramos a função T t uma constante, e vice versa), temos:

y x = X ' T y t = X T ' ⟹ y x x = X ' ' T y t t = X T ' '

Passo 2

Bom, brincadeiras de lado, agora que encontramos as derivadas, vamos substituir na EDP que rege o nosso problema. Com isto, temos:

y t = α 2 y x x

X T ' ' = α 2 X ' ' T

Separando tudo que tem T no primeiro membro de tudo que tem X no segundo, temos:

1 α 2 T ' ' T = X ' ' X

A igualdade acima é bem interessante. No primeiro membro temos uma função apenas de t . Já no segundo membro temos uma função de x (lembrando que estamos considerando o coeficiente k como uma constante!).

Isto implica que a equação acima só é verdade se ambos os membros forem iguais a uma constante λ (funções de variáveis distintas iguais entre si, só podem ser iguais a uma constante). Portanto, temos:

X ' ' X = 1 α 2 T ' ' T = λ ⟹ X ' ' - λ X = 0 T ' ' - α 2 λ T = 0

E agora simplificamos o nosso problema. Transformarmos uma EDP em duas EDO’s homogêneas.

Estamos quase chegando no sistema que a questão pede para demonstrarmos. Só nos resta determinarmos as condições de contorno. Para isto, temos:

y 0 , t = X 0 T t = 0 y L , t = X L T t = 0 y x , 0 = X x T 0 = 0

Observem as equações acima, as quais consistem em um produto de funções igual a zero. Nelas, ou o termo T é zero ou o termo X será igual a zero.

Como queremos que a função y x , t não seja igual a zero para qualquer valor de t ou de x , temos que considerar os termos T t e X x diferentes de zero. Portanto, obtemos que:

X 0 = X L = 0 T 0 = 0

Logo, essas são as nossas condições de contorno, e temos o seguinte sistema:

X ' ' - λ X = 0 X 0 = X L = 0 T ' - α 2 λ T = 0 T 0 = 0

A única diferença deste sistema para o da questão é o sinal negativo nas duas EDO’s. Como não sabemos nada sobre o sinal de λ (é até a próxima discussão nossa), podemos dizer que λ = - λ (só trocar o sinal nas equações), e portanto,

X ' ' + λ X = 0 X 0 = X L = 0 T ' + α 2 λ T = 0 T 0 = 0

E assim demonstramos como se chega no sistema da questão!

Passo 3

Nosso trabalho agora é resolvê-las utilizando as condições de contorno do enunciado do problema. Vamos começar resolvendo a EDO em relação a X . Somente depois vamos partir para o caso da EDO em função de T .

Com isto, temos o seguinte problema:

X ' ' + λ X = 0 X 0 = X L = 0

E quem é λ ? Ele é positivo? Ou nulo? Será se é negativo??? Está claro que nada sabemos sobre ele. Consequentemente, temos que testar todas as possibilidades!

OBS: podemos pensar que a questão está nos obrigando a verificar essas possibilidades, talvez este seja o maior estímulo...

Sendo assim, vamos as nossas hipóteses! A começar, temos λ < 0 λ = - μ 2 .

Para este caso, a nossa equação característica se transforma em:

r 2 - μ 2 = 0 ⟹ r = ± μ

Como Δ > 0 , a solução desta EDO possui a seguinte forma: X x = c 1 e μ x + c 2 e - μ x .

Para descobrirmos os valores dos coeficientes c 1 e c 2 vamos utilizar as condições X 0 = X L = 0 . Substituindo as condições de contorno, temos o seguinte sistema:

X 0 = c 1 + c 2 = 0 X L = c 1 e μ L + c 2 e - μ L = 0

Ao resolvermos o sistema acima, vamos obter que c 1 = c 2 = 0 . Esta solução é inviável, pois nos retornaria X x = 0 , independente do valor de x , o qual foi o mesmo problema que tivemos no Passo 2, só que em relação a função T t .

Portanto, a EDP não tem solução para λ < 0 !

Passo 4

Agora vamos para a segunda! Vamos assumir que λ = 0 . E com isto, vamos repetir o mesmo procedimento da hipótese anterior. Só que, neste caso, a equação característica transforma em:

r 2 = 0

Como Δ = 0 , a solução desta EDO possui a seguinte forma: X x = c 3 + c 4 x . Substituindo as condições de contorno, temos o seguinte sistema:

X 0 = c 3 = 0 X L = c 4 L = 0

Cuja solução é c 3 = c 4 = 0 . Da mesma forma que ocorreu no Passo 3, esta solução é considerada inviável, pois nos retornaria X x = 0 .

Logo, a EDP não tem solução para λ = 0 .

Passo 5

Por fim, vamos para a última possibilidade. Vamos assumir que λ > 0 λ = μ 2 . E com isto, vamos repetir o mesmo procedimento. Neste caso, a equação característica transforma em:

r 2 + μ 2 = 0 ⟹ r = ± μ i

Como Δ < 0 , a solução desta EDO possui a seguinte forma: X x = c 5 sen ⁡ ( μ x ) + c 6 cos ⁡ μ x . Substituindo as condições de contorno, temos o seguinte sistema:

X 0 = c 6 = 0 X L = c 5 sen ⁡ μ L + c 6 cos ⁡ μ L = 0

Da primeira equação verificamos que c 6 = 0 . Com isto, a nossa segunda equação se torna:

X L = c 5 sen ⁡ μ L = 0

E isto nos dá duas possibilidades: c 5 = 0 ou sen ⁡ μ L = 0 . É claro que a primeira possibilidade não nos interessa, pois teríamos X x = 0 , independente do valor de x . Portanto, vamos verificar quais os valores de μ que satisfazem sen ⁡ μ L = 0 . Portanto, temos:

sen ⁡ μ L = 0 ⟹ μ L = n π

μ = n π L

Sendo n qualquer número inteiro. Com isto, o valor de λ é:

λ = μ 2 = n 2 π 2 L 2

Guardem o valor de λ que vamos utilizar para solucionar a função T t ! Para a nossa função X vamos finalizando por aqui. O que obtivemos foi:

X x = c 5 sen ⁡ n π x L

E quem é essa constante c 5 ? Bom, na real, podemos desconsiderar ela por enquanto. Vamos trabalhar com as soluções fundamentais, as quais são da forma:

X n x = sen ⁡ n π x L

Passo 6

Bom, só nos resta agora encontrar a solução para a função T (t). Para isto, vamos resolver a EDO de segunda ordem associada a ela, a qual é:

T ' ' + α 2 λ T = 0 ⟹ T ' ' + n 2 π 2 α 2 L 2 T = 0

Cuja equação característica é:

r 2 + n 2 π 2 α 2 L 2 = 0 ⟹ r = ± n π α L i

E, como já vimos anteriormente, a sua solução é:

T t = c 7 sen ⁡ n π α L t + c 8 cos ⁡ n π α L t

Precisamos aplicar agora as sua condição de contorno, T 0 = 0 . Isto nos da:

T 0 = c 8 = 0

Portanto, a solução é:

T t = c 7 sen ⁡ n π α L t

Novamente, vamos trabalhar com as soluções fundamentais aqui. Portanto,

T n t = sen ⁡ n π α L t

Então, para cada n pertencente aos naturais, quando λ > 0 , temos a seguinte solução fundamental:

y n x , t = sen ⁡ n π x L sen ⁡ n π α L t

Por fim, a solução geral é obtida através de uma combinação linear das duas soluções fundamentais obtidas anteriormente. Portanto, a solução para a nossa EDP será:

y x , t = ∑ n = 1 ∞ c n u n x , t = ∑ n = 1 ∞ c n sen ⁡ n π x L sen ⁡ n π α L t

E a letra a) está resolvida! Foi uma resolução bem grande, com muitos procedimentos. Mas no final, ela acabou sendo bastante repetitiva. Só devemos tomar muito cuidado para não errarmos nenhuma conta!

Passo 7

E para finalizarmos a questão, devemos resolver a letra b). Agora, só devemos indicar como que calcularemos o coeficiente c n que apareceu na solução da nossa EDP.

Aaah! Não podemos esquecer de substituir L = 5 , k = 3 e g x = x 2 viu!!! Primeiro, vamos determinar a derivada y t x , t e calcularmos y t x , 0 . Logo,

y x , t = ∑ n = 1 ∞ c n sen ⁡ n π x L sen ⁡ n π α L t

y t x , t = ∑ n = 1 ∞ c n n π α L sen ⁡ n π x L cos ⁡ n π α L t

y t x , 0 = ∑ n = 1 ∞ c n n π α L sen ⁡ n π x L

Igualando y t x , 0 = g x , temos:

∑ n = 1 ∞ c n n π α L sen ⁡ n π x L = x 2

Então, os coeficientes c n são os coeficientes da série de Fourier em senos de f x (extensão ímpar). Portanto, ele será:

c n n π α L = 2 L ∫ 0 L f x sen ⁡ n π x L d x

c n = 2 n π α ∫ 0 L f x sen ⁡ n π x L d x

c n = 2 3 n π ∫ 0 5 x 2 sen ⁡ n π x 5 d x

O que nos daria um trabalho calcular... Mas, para a nossa alegria, não precisamos calcular nada! Basta indicarmos como calcular ambos os coeficientes! Logo... finalizamos a questão!!!

Resposta

  1. Para λ < 0 e para λ = 0 , não temos nenhuma solução fundamental. Já para λ > 0 , obtivemos y n x , t = sen ⁡ n π x L sen ⁡ n π α L t
  2. Os coeficientes da solução u x , t são calculados através das seguintes integrais:

c n = 2 3 n π ∫ 0 5 x 2 sen ⁡ n π x 5 d x

Exercícios de Livros Relacionados

Os Problemas 19 e 20 indicam como é possível mostrar que a solução formal (20), (22) das Eqs. (1), (3) e (9) é, de fato, a solução desse problema.19- Usando a identidade trigonométrica s e nAc o sB=½
Seja h ( ξ ) o deslocamento inicial em [ 0 ,L ] estendido a ( - L ,0 ) como uma função ímpar e estendido ao resto da reta como uma função periódica de período 2 L . Supondo a continuidade de h e suas
Combinando os resultados dos Problemas 16 e 17, mostre que a solução do problemaé dada por
O movimento de uma membrana circular elástica, como a membrana de um tambor, é determinado pela equação de onda bidimensional em coordenadas polaresSupondo que u ( r ,θ ,t )=R r Θ ( θ ) T ( t ) , enco
Considere a equação de onda a 2 u x x = u t t em um meio unidimensional infinito, sujeita às condições iniciaisUsando a forma da solução obtida no Problema 13, mostre queUse a primeira das equações no