23/08/2010

EDO: Queda dos Corpos com Resistência do Ar

Geralmente tratamos problemas de Física desprezando a resistência do ar para facilitar nossos cálculos. No entanto, essa resistência existe e, aqui, vamos ver uma equação diferencial para a queda de corpos considerando a resistência do ar.

Consideremos um corpo de massa $m$ em queda livre vertical, onde atuam somente a força da gravidade $g$ e a resistência do ar proporcional à velocidade do corpo. Precisamos admitir que tanto a gravidade como a massa deste corpo permaneçam constantes durante a queda.

Já vimos em outro artigo deste blog que, para corpos em queda livre, devemos adotar um sentido (para baixo ou para cima) como positivo. Convenientemente, vamos adotar o sentido para baixo como sendo sentido positivo.

edo-queda-dos-corpos-em-queda-livre-com-resistencia-do-ar

Segundo a Segunda Lei de Newton: A força resultante que atua sobre um corpo é igual à taxa de variação da quantidade de movimento (momentum) do corpo, ou, para uma massa constante.
$$
F= m \frac{dv}{dt} \tag{1}
$$
onde $F$ é a força resultante que atua sobre o corpo e $v$ é a velocidade do corpo, consideradas no instante $t$.

Em um corpo em queda livre, há duas forças atuando sobre ele: a primeira é a força da gravidade $g$, dada pelo peso $p$ do corpo, que é igual a $mg$:
$$
p = mg \tag{2}
$$
A segunda força é devida à resistência do ar, dada por:
$$
-kv \tag{3}
$$
onde $k \geqslant 0$ é uma constante de proporcionalidade. O sinal negativo se dá pelo fato desta força estar atuando no sentido contrário à queda do corpo, se opondo à velocidade, atuando no sentido para cima.
edo-queda-dos-corpos-em-queda-livre-com-resistencia-do-ar
A força resultante será:
$$
F = mg - kv \tag{4}
$$
Se substituirmos $(4)$ em $(1)$, obteremos:
$$
mg - kv = m \frac{dv}{dt}\\
\ \\
mg = m \frac{dv}{dt} + kv
$$
Isolando a força da gravidade:
$$
g = \frac{dv}{dt} + v \frac{k}{m} \tag{5}
$$
Agora, $(5)$ é a equação de movimento do corpo. Para resolver esta equação, usamos a técnica de fator integrante dado por:
$$
\mu = e^{kt/m} \tag{6}
$$
Multiplicando $(6)$ por $(5)$, obtemos:
$$
e^{kt/m} \frac{dv}{dt} + \frac{k}{m}e^{kt/m}v = ge^{kt/m} \tag{7}
$$
Observe que o lado esquerdo de $(7)$ é derivada da função:
$$
w(t) = e^{kt/m}v \tag{8}
$$
Pois:
$$
\frac{dw}{dt} = \frac{k}{m} e^{kt/m}v + e^{kt/m} \frac{dv}{dt} \tag{9}
$$
Assim, comparando $(7)$ com $(9)$, segue que:
$$
\frac{dw}{dt} = g \Rightarrow w(t) = \int dw = \int ge^{kt/m} dt = \frac{mg}{k}e^{kt/m} + C
$$
ou seja:
$$
e^{kt/m} v = \frac{mg}{k} e^{kt/m}+C\\
\ \\
v(t) = \frac{mg}{k} + Ce^{-kt/m}
$$
Usando o fato que $v(0)=v_0$, segue que:
$$
v(t) = \frac{mg}{k} + \left(v_0 - \frac{mg}{k}\right) e^{-kt/m}
$$
que é a velocidade em cada instante.

Vejam que, se $k=0$, temos que a resistência do ar é desprezível e a equação $(5)$ se reduz a:
$$
g = \frac{dv}{dt} \tag{10}
$$
Se $k > 0$, a velocidade limite $v_1$ é obtida fazendo $\cfrac{dv}{dt} = 0$, pois na há um equilíbrio da força da gravidade com a força devida ao atrito do ar. Assim:
$$
v_1 = \frac{mg}{k} \tag{11}
$$
Uma observação importante é que as equações $(5)$ e $(11)$ são válidas somente se as condições dadas forem satisfeitas. Se a resistência do ar não for proporcional à velocidade e sim ao quadrado da velocidade.

Exemplo:

Deixa-se cair um corpo de massa de $5\ kg$ de uma altura de $100\ m$, com velocidade inicial zero. Supondo que não haja resistência do ar, determine:

  1. A expressão da velocidade do corpo no instante $t$;
  2. A expressão da posição do corpo no instante $t$;
  3. O tempo necessário para o corpo atingir o solo.

Primeiramente adotamos o sistema de coordenada como na figura abaixo, sendo positivo o sentido para baixo.
exermplo-edo-queda-dos-corpos-em-queda-livre-com-resistencia-do-ar

Como nao há resistência do ar, usamos a equação $(10)$:
$$
g = \frac{dv}{dt}
$$
Essa é uma equação linear separável. Assim:
$$
dv = g\ dt\\
\ \\
\int dv = \int g\ dt\\
\ \\
v = gt + C
$$
1. Como $v(0) = 0$, segue que $v(t) = gt$.

2. Para determinar a expressão da posição $x$ no instante $t$, fazemos:
$$
\frac{dx}{dt} = gt\\
\ \\
dx = gt\ dt\\
\ \\
\int dx = g\int t\ dt\\
\ \\
x(t) = \frac{gt^2}{2}+C
$$
Sendo $x(0)=0$, segue que:
$$
x(t) = \frac{gt^2}{2}
$$
3. Para $x(t)=100$, fazemos:
$$
100 = \frac{gt^2}{2}
$$
Se adotarmos a aproximação $g=9,8\ m/s^2$, teremos:
$$
100 = \frac{9,8t^2}{2}\\
\ \\
t = \sqrt{20,4}\\
\ \\
t \approx 4,5\ s
$$

Veja mais:

COMO REFERENCIAR ESSE ARTIGO: Título: EDO: Queda dos Corpos com Resistência do Ar. Publicado por Kleber Kilhian em 23/08/2010. URL: . Leia os Termos de uso.


Siga também o blog pelo canal no Telegram.
Achou algum link quebrado? Por favor, entre em contato para reportar o erro.
Para escrever em $\LaTeX$ nos comentários, saiba mais em latex.obaricentrodamente.com.

35 comentários:

  1. A versão final ficou muito boa. Obrigado por citar o blog e o post sobre a velocidade terminal de um pára-quedas. Abraços!

    ResponderExcluir
  2. Que bom que gostou. Sem sua ajuda acho que não sairia dessa forma.

    Um abraço!

    ResponderExcluir
  3. Interessante que você utilizou -kv para a força de arrasto. No livro Dinâmica Clássica de Partículas, Marion e Thornton pág. 52, é utilizado -kmv, a massa então é cancelada em ambos os lados da equação e a solução não depende da massa.

    ResponderExcluir
  4. humm, vou pesquisar este livro para ver o desenvolvimento.

    Agradeço seu comentário e observação.

    Um abraço.

    ResponderExcluir
  5. Uhuuulll...
    Estavamos queimando nossos neurônios aqui... Valeu meeeesmo!!!!
    Agora vamos formar... \o

    ResponderExcluir
  6. Oi!. Se não houvesse o atrito a equação seria a=g certo? Mas isto não resultaria em uma equação horária onde y cresceria indefinidamente $y=y_{0}+v_{0}t+\frac{g}{2}t^2$?

    ResponderExcluir
  7. Sim, mas se adotamos altura inicial sendo zero e a velocidade inicial igual a zero, caímos na mesma equação: $x(t)=\dfrac{gt^2}{2}$.

    ResponderExcluir
  8. Oi! Você poderia escrever para mim os vetores aceleração, velocidade e posição?

    ResponderExcluir
  9. É que estou encafifado com os sinais. Me parece que na letra c) do exemplo você trocou pois se a pedra caiu de uma altura de 100 m então temos $x(0)=100$ e $x(t)=0$....

    ResponderExcluir
  10. Adriano, veja que o eixo vertical está orientado para baixo, sendo positivo para baixo. para $x(t)=0$ temos $t=0$, que é no instante em que o corpo é largado de uma altura de 100m. Para $x(t)=100$, ou seja, a distância total até o solo, gasta-se 4,5s.

    Tem uma regrinha que usava, vou tentar lembrar agora:

    Orientação Vertical de Baixo para Cima
    1) Na descida, a velocidade é negativa;
    2) $g$ é sempre negativa;
    3) $y_0$ está no ponto de origem da queda do corpo;
    4) $y$ é o ponto final do deslocamento do corpo.

    Lançamento Vertical
    1) Na subida, a velocidade é positiva;
    2) $g$ é sempre negativa;
    3) $y_0$ é o ponto de origem do lançamento;
    4) $y$ é o ponto máximo que o corpo atinge.

    Com isso podemos montar uma figura representativa que facilita a visualização.

    ResponderExcluir
    Respostas
    1. Veja as imagnes neste link:

      http://obaricentrodamente.blogspot.com.br/2009/05/equacao-do-movimento.html

      Excluir
  11. Isso é coisa do tinhoso...
    Eu não acredito que vou ter que aprender tudo isso!
    Me salva, me diz que só vou ter aprender isso beeeem mais tarde... (sou nono ano (fundamental II)

    ResponderExcluir
    Respostas
    1. Você só terá que aprender isso se for fazer exatas na faculdade.

      Excluir
  12. Eduardo Muller24/11/12 17:45

    Profesor, uma dúvida... Vi em um livro de física que a força de arrasto é proporcional ao quadrado da velocidade (mais precisamente F = C . r . A . v^2 onde C é o coeficiente de arasto, r é a densidade do ar e A a área da seção transversal do corpo) mas aqui está dizendo que a força é apenas proporcional à velocidade... Como definimos o coeficiente k então?

    ResponderExcluir
  13. Eduardo, este artigo não mostra como determinar o coeficiente de arrasto, já que este é determinado experimentalmente. Veja os links do Wikipédia:

    Equação do arrasto: http://pt.wikipedia.org/wiki/Equa%C3%A7%C3%A3o_do_arrasto

    Coeficiente do arrasto: http://pt.wikipedia.org/wiki/Coeficiente_de_arrasto

    Um abraço.

    ResponderExcluir
  14. Num exercicio ele me da g=dv/dt + (kv^2)/m, qu método de integração eu utilizo, ja que tenho v^2?

    ResponderExcluir
    Respostas
    1. Olá amigo, neste caso, faça a separação das variáveis. Veja esta EDO resolvida neste artigo:

      http://fatosmatematicos.blogspot.com.br/2010/04/velocidade-terminal-de-um-para-quedas.html

      Abraços!

      Excluir
  15. se eu mudasse a forma geometrica do objeto que está em queda livre, por exemplo, uma pessoa saltando e depois abrindo o para-quedas... consequentemente a restistencia do ar aumentaria, como eu resolveria esse problema??

    obrigado!

    ResponderExcluir
  16. Veja o problema do para-quedas no link abaixo:

    http://fatosmatematicos.blogspot.com.br/2010/04/velocidade-terminal-de-um-para-quedas.html

    Abraços!

    ResponderExcluir
    Respostas
    1. precisa ser convidado para olhar.
      meu email: pcaetaninha@gmail.com

      Excluir
  17. Belo texto. Parabéns. A algum tempo atrás, eu escrevi um pequeno texto para o Portal Píon da Sociedade Brasileira de Física reportando o mesmo problema, porém especificamente, perguntando se é mais rápido para um corpo subir ou descer.

    (http://www.sbfisica.org.br/v1/pion/index.php/publicacoes/artigos/54-subir-ou-descer)

    Disponibilizei-o recentemente no meu blog (http://jefsrodrigues.blogspot.com.br/2012/11/o-que-e-mais-rapido-subir-ou-descer.html). São temas relacionados, fantástico!

    Att.


    ResponderExcluir
  18. o que vem a ser o w na equação?? e se eu no caso for proporcional a velocidade ao quadrado (kv²) é so substituir v por v² nas equaçoes?

    ResponderExcluir
  19. Veja que para resolver a equação (5), usamos o fator integrante dado em (6), de modo que $w$ é uma mudança de variável.

    Para $v^2$, neste caso, faça a separação das variáveis. Veja esta EDO resolvida neste artigo:

    http://fatosmatematicos.blogspot.com.br/2010/04/velocidade-terminal-de-um-para-quedas.html

    Abraços.

    ResponderExcluir
    Respostas
    1. mas teria algum problema nas equações se substituir diretamente v por v²?

      desde já , obrigado.

      Excluir
    2. Sim. Para cada caso, tem um método específico de resolução. PAra $v^2$ faz-se a separação de variáveis e para $v$ usa-se o fator integrante.

      Excluir
    3. Obrigado. Você sabe me dizer quando se utiliza v ou v². Obrigado de novo!!

      Excluir
  20. Caro Professor,tem um exercicio semelhante a este que não estou conseguindo resolver (traduzido do meirovitch-Methods of Analitical Dynamics):

    Uma partícula de massa m é lançada verticalmente com velocidade V num campo gravitacional uniforme g. Considerando o atrito viscoso do ar diretamente
    proporcional à velocidade ( f = -cV ) determine o tempo necessário para que a partícula retorne ao ponto de partida. Com que velocidade essa partícula retorna?

    O Sr poderia, por favor, dar uma ajuda?

    ResponderExcluir
  21. Interessante desenvolvimento. Em modelos mais avançados, tem-se que a força de arrasto é proporcional ao quadrado da velocidade, mas aí tem-se o inconveniente da EDO do sistema não ser mais linear.

    ResponderExcluir
  22. Anônimo1/1/15 00:46

    Tem um livro, acho, Kazunoru Watari, da USP, que fala mais dessa exponencial em queda livre, bem como correlatos. É mais complicado que o $$h=gt^2/2$$ do colegial. Excelente.

    ResponderExcluir
  23. OLA GOSTARIA DE SABER QUAL A INFLUENCIA DA RESISTENCIA DO AR

    ResponderExcluir
  24. oi pessoal queria uma ajuda,como isolar a massa na formula da velocidade do paraquedista em queda live v(t) = g * m * ( 1 – e – (c / m) * t))/ c

    ResponderExcluir
  25. Professor, não entendi a parte que você considerou o fato que v(0) = v0, sei que usou isso porque se trata de uma queda livre e portanto a velocidade inicial é igual a zero, mas não entendi porque você substituiu a constante c por $v(0) - (mg/k)$. Fazendo as contas eu achei que a constante c deveria ser substituída por $(-mg/k) * e^(k/m)$ quando V(0)=0. Como você conseguiu substituir dessa forma? Obrigado desde já pela atenção!

    ResponderExcluir
  26. Poderia ser mais claro(Passo a Passo), para leigos.

    ResponderExcluir
  27. Eu pergunteiao meu professor ele me explicou q a substituiçao de v^2 por v e pq v^2 e a velocidade relativa e ela depende da velocidade do movel q tbm varia sendo assim podemos. Substituir v^2por v(t) e utilizar o fator integrante e resolver a equaçao

    ResponderExcluir

Whatsapp Button works on Mobile Device only

Pesquise no blog