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.
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=mdvdtonde 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=mgA segunda força é devida à resistência do ar, dada por:
−kvonde k⩾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.
A força resultante será:
F=mg−kvSe substituirmos (4) em (1), obteremos:
mg−kv=mdvdt mg=mdvdt+kvIsolando a força da gravidade:
g=dvdt+vkmAgora, (5) é a equação de movimento do corpo. Para resolver esta equação, usamos a técnica de fator integrante dado por:
μ=ekt/mMultiplicando (6) por (5), obtemos:
ekt/mdvdt+kmekt/mv=gekt/mObserve que o lado esquerdo de (7) é derivada da função:
w(t)=ekt/mvPois:
dwdt=kmekt/mv+ekt/mdvdtAssim, comparando (7) com (9), segue que:
dwdt=g⇒w(t)=∫dw=∫gekt/mdt=mgkekt/m+Cou seja:
ekt/mv=mgkekt/m+C v(t)=mgk+Ce−kt/mUsando o fato que v(0)=v0, segue que:
v(t)=mgk+(v0−mgk)e−kt/mque é 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=dvdtSe k>0, a velocidade limite v1 é obtida fazendo dvdt=0, pois na há um equilíbrio da força da gravidade com a força devida ao atrito do ar. Assim:
v1=mgkUma 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:
- A expressão da velocidade do corpo no instante t;
- A expressão da posição do corpo no instante t;
- 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.
Como nao há resistência do ar, usamos a equação (10):
g=dvdtEssa é uma equação linear separável. Assim:
dv=g dt ∫dv=∫g dt v=gt+C1. Como v(0)=0, segue que v(t)=gt.
2. Para determinar a expressão da posição x no instante t, fazemos:
dxdt=gt dx=gt dt ∫dx=g∫t dt x(t)=gt22+CSendo x(0)=0, segue que:
x(t)=gt223. Para x(t)=100, fazemos:
100=gt22Se adotarmos a aproximação g=9,8 m/s2, teremos:
100=9,8t22 t=√20,4 t≈4,5 s
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!
ResponderExcluirQue bom que gostou. Sem sua ajuda acho que não sairia dessa forma.
ResponderExcluirUm abraço!
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.
ResponderExcluirhumm, vou pesquisar este livro para ver o desenvolvimento.
ResponderExcluirAgradeço seu comentário e observação.
Um abraço.
valeu, muito bom
ResponderExcluirUhuuulll...
ResponderExcluirEstavamos queimando nossos neurônios aqui... Valeu meeeesmo!!!!
Agora vamos formar... \o
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=y0+v0t+g2t2?
ResponderExcluirSim, mas se adotamos altura inicial sendo zero e a velocidade inicial igual a zero, caímos na mesma equação: x(t)=gt22.
ResponderExcluirOi! Você poderia escrever para mim os vetores aceleração, velocidade e posição?
ResponderExcluirÉ 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....
ResponderExcluirAdriano, 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.
ResponderExcluirTem 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) y0 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) y0 é 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.
Veja as imagnes neste link:
Excluirhttp://obaricentrodamente.blogspot.com.br/2009/05/equacao-do-movimento.html
Isso é coisa do tinhoso...
ResponderExcluirEu 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)
Você só terá que aprender isso se for fazer exatas na faculdade.
ExcluirProfesor, 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?
ResponderExcluirEduardo, este artigo não mostra como determinar o coeficiente de arrasto, já que este é determinado experimentalmente. Veja os links do Wikipédia:
ResponderExcluirEquaçã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.
Num exercicio ele me da g=dv/dt + (kv^2)/m, qu método de integração eu utilizo, ja que tenho v^2?
ResponderExcluirOlá amigo, neste caso, faça a separação das variáveis. Veja esta EDO resolvida neste artigo:
Excluirhttp://fatosmatematicos.blogspot.com.br/2010/04/velocidade-terminal-de-um-para-quedas.html
Abraços!
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??
ResponderExcluirobrigado!
Veja o problema do para-quedas no link abaixo:
ResponderExcluirhttp://fatosmatematicos.blogspot.com.br/2010/04/velocidade-terminal-de-um-para-quedas.html
Abraços!
precisa ser convidado para olhar.
Excluirmeu email: pcaetaninha@gmail.com
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.
ResponderExcluir(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.
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?
ResponderExcluirVeja que para resolver a equação (5), usamos o fator integrante dado em (6), de modo que w é uma mudança de variável.
ResponderExcluirPara v2, 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.
mas teria algum problema nas equações se substituir diretamente v por v²?
Excluirdesde já , obrigado.
Sim. Para cada caso, tem um método específico de resolução. PAra v2 faz-se a separação de variáveis e para v usa-se o fator integrante.
ExcluirObrigado. Você sabe me dizer quando se utiliza v ou v². Obrigado de novo!!
ExcluirCaro Professor,tem um exercicio semelhante a este que não estou conseguindo resolver (traduzido do meirovitch-Methods of Analitical Dynamics):
ResponderExcluirUma 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?
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.
ResponderExcluirTem um livro, acho, Kazunoru Watari, da USP, que fala mais dessa exponencial em queda livre, bem como correlatos. É mais complicado que o h=gt2/2
OLA GOSTARIA DE SABER QUAL A INFLUENCIA DA RESISTENCIA DO AR
ResponderExcluiroi 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
ResponderExcluirProfessor, 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!
ResponderExcluirPoderia ser mais claro(Passo a Passo), para leigos.
ResponderExcluirEu 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