自由落下運動 - 空気抵抗有無の比較1

これまでに導いた空気抵抗無しと有りの2つの自由落下運動を比較してみよう。導いた一般解を再度書き出す。

 

ここでは落下物の速度、は位置、は重力加速度、は時間、は空気抵抗係数、は落下物の質量、は任意定数である。いずれも上方向(重力と逆方向)をとしていることに注意。

 

空気抵抗なしの自由落下

一般解

前回求めた一般解を再掲する。

 

速度

\begin{eqnarray}
v=-gt+v_0 \tag{1}
\end{eqnarray}

 

位置

\begin{eqnarray}
x=-\frac{1}{2}gt^2+v_0t+x_0 \tag{2}
\end{eqnarray}

 

 

特殊解

初期値を代入して特殊解を作る。

 

今回はにおいて、基準点から静かに落下させたとする。すなわち初期値はでになる。

 

空気抵抗なしの場合、式(1)より、式(2)よりである。

 

速度

\begin{eqnarray}
v=-gt
\end{eqnarray}

 

位置

\begin{eqnarray}
x=-\frac{1}{2}gt^2
\end{eqnarray}

 

 

空気抵抗ありの自由落下

一般解

前回求めた一般解を再掲する。

 

速度

\begin{equation}
v=C_1\exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k} \tag{3}
\end{equation}

 

位置

\begin{eqnarray}
x&=&-\frac{mC_1}{k}\exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k}t+C_2 \tag{4}
\end{eqnarray}

 

特殊解

同様に、初期値をでとして任意定数を定める。

 

任意定数の計算

式(3)にを代入して、

\begin{eqnarray}
v&=&C_1\exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k} \tag{3}\\
0&=&C_1\exp(0)-\frac{mg}{k}\\
C_1&=&\frac{mg}{k}
\end{eqnarray}

 

式(4)にを代入して、

\begin{eqnarray}
x&=&-\frac{mC_1}{k}\exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k}t+C_2 \tag{4}\\
0&=&-\frac{m}{k}\frac{mg}{k}\exp(0)-0+C_2\\
C_2&=&\frac{m^2g}{k^2}
\end{eqnarray}

 

これらの定数を一般解に代入すると特殊解になる。

 

速度

\begin{eqnarray}
v&=&\frac{mg}{k} \exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k}\\
&=&\frac{mg}{k} \left[ \exp{\left(-\frac{kt}{m}\right)}-1 \right]
\end{eqnarray}

 

位置

\begin{eqnarray}
x&=&-\frac{m}{k} \frac{mg}{k} \exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k}t+\frac{m^2g}{k^2}\\
&=&-\frac{m^2g}{k^2} \exp{\left(-\frac{kt}{m}\right)}-\frac{mg}{k}t+\frac{m^2g}{k^2}
\end{eqnarray}

 

何が言いたいか

同じような微分方程式から生み出されたにもかかわらず、導かれた速度や位置の式は空気抵抗の有無で大きく異なっている。次回これらの運動をグラフにプロットし、本当に似た運動を示しているのか確かめる。