공기 저항이 속력에 비례하는 경우는 물체의 궤적은 closed form이 있다. 그러나 저항력이 속력의 제곱에 비례하게 주어지는 경우는 수치적으로 해결해야 한다. 움직이는 방향의 단면적이 $A=\frac{1}{4}\pi D^2$인 물체가 밀도가 $\rho$인 공기 속에서 $\vec{v}$의 속도로 움직일 때 저항력은
$$ {\vec F}_D = \frac{1}{4} \rho A v \vec{v}=\frac{1}{16}\pi \rho D^2 v\vec{v}=c v \vec{v}$$
로 표현할 수 있다. 따라서 물체의 운동방정식은
$$ m \ddot{\vec r} = m\vec{g}- c v \vec{v},$$
또는 성분으로 쓰면
$$ m \ddot{x} = - c \sqrt{ \dot{x}^2 + \dot{y}^2} \dot{x}, $$
$$m \ddot{y} = -mg - c \sqrt{\dot{x}^2+ \dot{y}^2} \dot{y}$$
로 주어진다. 아래의 mathematica 코드는 구체적인 수치(SI-단위 기준, 발사각 $\theta_0$, 발사속력 $v_0$)를 대입해서 공기저항이 있을 때와 없을 때 물체의 궤적을 보여준다.
$(1,1)$에서 $(0,0)$까지를 연결하는 일차 곡선 ($y=x$), 이차 곡선 ($y=x^2$), 6차 곡선 ($y=x^6$), 사분원 ($y=1-\sqrt{1-x^2}$) 그리고 cycloid ($x=1-R(\theta-\sin\theta), y=1-R(1-\cos\theta)$) 위를 움직이는 물체의 운동을 mathematica를 사용해서 animation으로 표현하는 방법을 알아보자. 중력가속도는 $g=1$로 놓는다. 적절한 generalized coordinate를 선택해서 Euler-Lagrange equation을 이용하면 운동방정식을 쉽게 구할 수 있다.
이 곡선이 $(0,0)$을 지나야 하는 조건에서 $R$과 그 때의 $\theta$ 값을 구할 수 있다. 먼저 $\theta$을 소거하면,
$$ R \cos \Big( \frac{1+\sqrt{2R-1}}{R} \Big) + 1 = R $$
을 얻고, 이 식의 근을 구하면 $R$ 값이 정해진다. 또한, $$1=R(\theta- \sin \theta)$$을 풀어서 $(0,0)$에 도달할 때 $\theta=\theta_0$ 값을 얻을 수 있다. 그리고 시간과 $\theta$ 변수의 관계는 위의 표현을 적분에 대입하면
$$t = \sqrt{ \frac{R}{g} } \theta $$
로 주어짐을 알 수 있다. 이는 사이클로이드가 일정하게 굴러가는 바퀴의 한 점이 그리는 자취이기 때문이다.
아래는 5가지 경우의 곡선 각각에서 운동을 보여주는 mathematica 코드다. 곡선의 모양에 따라 바닥에 도착하는 시간이 다름을 알 수 있다. 출발 높이가 $H$일 때 걸리는 시간은 $\sqrt{H/g}$ 단위로
직선: $t_1=\int_0^1{ \frac{dy}{\sqrt{1-y}}} = 2$
2차 곡선: $t_2=\int_0^1\sqrt{ \frac{1+4y}{8y(1-y)}}dy= 1.86336$