46
Cálculo III Aula 26 – Separação de Variáveis e a Equação da Onda. Marcos Eduardo Valle Depart. Matemática Aplicada IMECC – Unicamp Marcos Eduardo Valle MA311 – Cálculo III 1 / 16

Cálculo III - Aula 26 Separação de Variáveis e a Equação ...valle/Teaching/MA311/Aula26.pdfCálculo III Aula 26 – Separação de Variáveis e a Equação da Onda. Marcos Eduardo

  • Upload
    others

  • View
    8

  • Download
    0

Embed Size (px)

Citation preview

Cálculo IIIAula 26 – Separação de Variáveis e a Equação da Onda.

Marcos Eduardo ValleDepart. Matemática Aplicada

IMECC – Unicamp

Marcos Eduardo Valle MA311 – Cálculo III 1 / 16

Equação da Onda

Considere uma corda elástica de comprimento L presa nasextremidades em suportes de mesmo nível horizontal.

Vamos denotar por u(x, t) o deslocamento vertical da corda no ponto0 ≤ x ≤ L no instante t ≥ 0.

Desprezando efeitos de amortecimento e supondo que a amplitude domovimento não é grande, u satisfaz a equação diferencial parcial

a2uxx = utt ,

em que a é a velocidade de propagação de ondas ao longo da corda(depende da tensão e da massa por unidade de comprimento).

Para descrever o movimento da corda, precisamos também dascondições iniciais e de contorno.

Marcos Eduardo Valle MA311 – Cálculo III 2 / 16

Condições de Iniciais e de Contorno

Como as extremidades da corda permanecem fixas, as condições decontorno são:

u(0, t) = 0 e u(L , t) = 0.

As condições iniciais são:• Posição inicial:

u(x, 0) = f(x), 0 ≤ x ≤ L .

• Velocidade inicial:

ut(x, 0) = g(x), 0 ≤ x ≤ L .

em que f e g são funções tais que

f(0) = f(L) = 0 e g(0) = g(L) = 0.

Marcos Eduardo Valle MA311 – Cálculo III 3 / 16

Corda Elástica com Deslocamento Não-Nulo

Iniciaremos o estudo do problema de vibrações de uma corda elásticaadmitindo que a velocidade inicial é nula, ou seja,

g(x) = 0, ∀0 ≤ x ≤ L .

Em outras palavras, considere o problemaa2uxx = utt ,

u(0, t) = 0 e u(L , t) = 0,

u(x, 0) = f(x), 0 ≤ x ≤ L ,

ut(x, 0) = 0, 0 ≤ x ≤ L , .

em que f(0) = f(L) = 0 descreve a configuração inicial da corda.

Marcos Eduardo Valle MA311 – Cálculo III 4 / 16

Separação de Variáveis

Vamos admitir que u pode ser escrita como

u(x, t) = X(x)T(t),

em que X depende apenas de x e T depende somente de t .

Derivando e substituindo na equação diferencial parcial, obtemos

X ′′

X=

1a2

T ′′

T= −λ,

em que λ é uma constante de separação.

Equivalentemente, temos as equações diferenciais ordinárias

X ′′ + λX = 0 e T ′′ + a2λT = 0.

Marcos Eduardo Valle MA311 – Cálculo III 5 / 16

Usando a condição de contorno, encontramos o problema

X ′′ + λX = 0, X(0) = X(L) = 0,

cuja solução é

Xm(x) = sen(mπx

L

)e λm =

(mπL

)2, m = 1, 2, . . . .

Com as constantes de separação acima, obtemos a EDO

T ′′ + ω2T = 0, com ω =mπa

L,

cujas soluções são

T(t) = k1 cos(ωt) + k2 sen(ωt).

Marcos Eduardo Valle MA311 – Cálculo III 6 / 16

Como a velocidade inicial é nula, deduzimos

ut(x, 0) = X(x)T ′(0) = 0,∀ 0 ≤ x ≤ L =⇒ T ′(0) = 0.

ComoT ′(t) = −ωk1 sen(ωt) + ωk2 cos(ωt),

temosT ′(0) = 0 =⇒ k2 = 0.

Assim, as soluções fundamentais da equação da onda, com ascondições de contorno e a segunda condição inicial, são

um(x, t) = sen(mπx

L

)cos(ωt) = sen

(mπxL

)cos

(mπatL

),

para m = 1, 2, . . ..

Note que um é periódica no tempo com período 2L/ma.Marcos Eduardo Valle MA311 – Cálculo III 7 / 16

A superposição das soluções fundamentais fornece

u(x, t) =∞∑

m=1

cm sen(mπx

L

)cos

(mπatL

).

Finalmente, a condição inicial u(x, 0) = f(x), fornece

u(x, 0) =∞∑

m=1

cm sen(mπx

L

)= f(x).

Portanto, admitindo que f é uma função ímpar com período T = 2L ,concluímos que os coeficiente satisfazem

cm =2L

∫ L

0f(x) sen

(mπxL

)dx, m = 1, 2, . . . .

Marcos Eduardo Valle MA311 – Cálculo III 8 / 16

Concluindo, a solução do problemaa2uxx = utt ,

u(0, t) = 0 e u(L , t) = 0,

u(x, 0) = f(x), 0 ≤ x ≤ L ,

ut(x, 0) = 0, 0 ≤ x ≤ L , .

é

u(x, t) =∞∑

m=1

cm sen(mπx

L

)cos

(mπatL

),

em que

cm =2L

∫ L

0f(x) sen

(mπxL

)dx, m = 1, 2, . . . .

Marcos Eduardo Valle MA311 – Cálculo III 9 / 16

Observações

• A solução é a superposição de funções periódicas no tempo comperíodo 2L/ma.

• As quantidades mπa/L são chamadas frequências naturais dacorda.

• O fator sen(

mπxL

)é chamado modo natural de vibração.

• O período do modo natural de vibração 2L/m é chamadocomprimento da onda.

Marcos Eduardo Valle MA311 – Cálculo III 10 / 16

Exemplo 1

Considere uma corda vibrante de comprimento L = 30 que satisfaz aequação da onda

4uxx = utt , 0 < x < 30 e t > 0.

Suponha que as extremidades da corda estão fixas e que a corda écolocada em movimento sem velocidade inicial da posição inicial

u(x, 0) =

x10 , 0 ≤ x ≤ 10,(30−x)

20 , 10 < x ≤ 30.

Encontre o deslocamento u(x, t) da corda.

Marcos Eduardo Valle MA311 – Cálculo III 11 / 16

Exemplo 1

Considere uma corda vibrante de comprimento L = 30 que satisfaz aequação da onda

4uxx = utt , 0 < x < 30 e t > 0.

Suponha que as extremidades da corda estão fixas e que a corda écolocada em movimento sem velocidade inicial da posição inicial

u(x, 0) =

x10 , 0 ≤ x ≤ 10,(30−x)

20 , 10 < x ≤ 30.

Encontre o deslocamento u(x, t) da corda.

Resposta: A solução é

u(x, t) =9π2

∞∑m=1

1m2

sen(mπ

3

)sen

(mπx30

)cos

(2mπt

30

).

Marcos Eduardo Valle MA311 – Cálculo III 11 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=0.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=1.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=2.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=3.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=4.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=5.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=6.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=7.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=8.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=9.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=10.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=11.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=12.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=13.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=14.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=15.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=16.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=17.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=18.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=19.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=20.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=21.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=22.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=23.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=24.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=25.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=26.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=27.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=28.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=29.00

Solução do Exemplo 1 para diferentes t .

Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

-1

-0.5

0

0.5

1

0 5 10 15 20 25 30

t=30.00

Solução do Exemplo 1 para diferentes t .Marcos Eduardo Valle MA311 – Cálculo III 12 / 16

Corda Elástica com Velocidade Não-Nula

Um problema semelhante ao discutido anteriormente, consiste noestudo das vibrações de uma corda que é colocada em movimento apartir do repouso com uma velocidade dada.

Formalmente, temos o problemaa2uxx = utt ,

u(0, t) = 0 e u(L , t) = 0,

u(x, 0) = 0, 0 ≤ x ≤ L ,

ut(x, 0) = g(x), 0 ≤ x ≤ L ,

em que g(x) é a velocidade inicial da corda no ponto x.

Marcos Eduardo Valle MA311 – Cálculo III 13 / 16

Procedendo de forma semelhante, concluímos que a solução é

u(x, t) =∞∑

m=1

km sen(mπx

L

)sen

(mπatL

),

em que

km =2

mπa

∫ L

0g(x) sen

(mπxL

)dx, m = 1, 2, . . . .

Esclarecemos que o fator 2/(mπa) multiplicando a integral acimaaparece porque precisamos identificar a série de Fourier em senos daderivada

ut(x, t) =∞∑

m=1

mπaL

km sen(mπx

L

)cos

(mπatL

),

com a série de Fourier em senos de g.

Marcos Eduardo Valle MA311 – Cálculo III 14 / 16

Problema Geral para a Corda Elástica

Finalmente, a solução do problema gerala2uxx = utt ,

u(0, t) = 0 e u(L , t) = 0,

u(x, 0) = f(x), 0 ≤ x ≤ L ,

ut(x, 0) = g(x), 0 ≤ x ≤ L ,

em que f(x) e g(x) descrevem, respectivamente, a posição e avelocidade inicial da corda no ponto x, é obtido superpondo assoluções dos problemas anteriores, ou seja,

u(x, t) = v(x, t) + w(x, t),

em que v e w são as soluções da corda elástica com deslocamentonão-nulo e velocidade não-nula, respectivamente.

Muito grato pela atenção!

Marcos Eduardo Valle MA311 – Cálculo III 15 / 16