오늘은 기본적인 자유 낙하 운동을 기술해보려고 한다.
질량이 m 인 물체를 지구 상공 어딘가 높은 곳에서 떨어뜨린다고 생각해보자. 이 때 공기저항은 속도에 비례한다. 따라서,
F=ma=mg−kv
여기서 k는 상수이고, 우리는 가속도 a가 속도 v의 시간미분임을 잘 알고 있다.
mdvdt=mg−kvdvmg−kv=dtm
양 변 적분해주면,
−1kln(mg−kv)=1mt+Cmg−kv=e−kmt+C∴
그러므로 시간이 충분히 지나 t \to \infty 일때, v = \frac{mg}{k} 로 수렴한다. 이를 종단속도 v_f 라고 부른다. 즉, 물체는 시간이 충분히 지나면 공기의 저항과 지구중력의 크기가 같아져 등속운동을 한다. 따라서 단순히 종단속도v_f만 구하고 싶다면, 두 힘은 같다고 놓고 풀면 쉽게 구해낼 수 있다.
F=ma=mg-kv_f = 0 \\ mg=kv_f \\ v_f = \frac{mg}{k}
만약 물체에 작용하는 저항이 속도의 제곱에 비례한다면 어떨까?
F = ma = mg - kv^2
마찬가지로 미분방정식을 풀어내면 된다.
m \frac{dv}{dt} = mg-kv^2 \\ \frac{dv}{mg-kv^2} = \frac{1}{m} dt \\ \int \frac{dv}{mg-kv^2} = \int \frac{1}{m} dt
이제 적분을 수행해야되는데, 속도 제곱에 비례하기 때문에 적분이 조금 까다롭다. 우리가 잘 아는 적분형태 \int \frac{1}{1+u^2} = arctan(u) 로 바꾸어 주자.
Let \ u = \frac{i\sqrt{k}v}{\sqrt{mg}}, \ du = \frac{i\sqrt{k}}{\sqrt{mg}} dv \\ \int \frac{dv}{mg-kv^2} \rightarrow -\frac{i}{\sqrt{mkg}} \int \frac{1}{u^2+1}
이제 마저 적분을 수행하면,
-\frac{i}{\sqrt{mkg}} \int \frac{1}{u^2+1} = \int \frac{1}{m} dt \\ -\frac{i \ arctan (u)}{\sqrt{mkg}} = \frac{t}{m} +C
이 때
tan(u) = \frac {sin (u)}{cos (u)} = \frac{\frac{e^{ix}-e^{-ix}}{2i}}{\frac{e^{ix}+e^{-ix}}{2}} \\ = \frac{e^{ix}-e^{-ix}}{i(e^{ix}+e^{-ix})}
이고,
tanh (iu) =\frac{sinh(iu)}{cosh{iu}} = \frac{e^{ix}-e^{-ix}}{e^{ix}+e^{-ix}}
이므로,
itan(u) = tanh(iu)
라는 사실을 유도할 수 있다. 따라서, i arctan (u) = tanh^{-1} (iu) 이고, 우리가 원래 풀던 미분방정식은
-i arctan(u) = \frac{t\sqrt{kg}}{\sqrt{m}} + C \\ -tanh^{-1}(iu) = \frac{t\sqrt{kg}}{\sqrt{m}} +C \\ tanh^{-1}(\frac{\sqrt{k}}{\sqrt{mg}}v)=\frac{t\sqrt{kg}}{\sqrt{m}} + C
식을 정리하면, 다음과 같은 결론을 얻을 수 있다.
\therefore v = \sqrt{\frac{mg}{k}} tanh (\sqrt{\frac{kg}{m}}t) \\ if \quad t \ \to \infty, \quad v_f = \sqrt{\frac{mg}{k}}
정리하자면, 물체에 적용되는 항력이 속도에 비례할 때 종단속력은 v_f = \frac{mg}{k}, 속도 제곱에 비례할 때는 \sqrt{\frac{mg}{k}} 라는 것을 구할 수 있었다.
각각에 대해서 시간에 따른 속도변화 그래프는 아래와 같다.