空気抵抗ありの斜方投射を計算する方法｜kvとkv²モデルの違いと近似計算の考え方

空気抵抗を考慮した斜方投射の計算では、抵抗力を速度に比例するモデル（kv）で扱うか、速度の二乗に比例するモデル（kv²）で扱うかによって計算方法が大きく変わります。特に現実の空気抵抗は一般的にkv²モデルが近いため、解析解が得られない場合にどのような近似を使えばよいのかが重要になります。

空気抵抗のモデルはなぜkvではなくkv²が一般的なのか

物体が空気中を運動するときに受ける抵抗力は、低速では速度に比例する場合がありますが、日常的な投射運動や落下運動では速度の二乗に比例する抵抗がよく用いられます。

空気抵抗をFとして表すと、速度比例モデルではF=-kv、速度二乗モデルではF=-kv²となります。ここでkは物体の形状、大きさ、空気密度などによって決まる係数です。

例えば小さな粒子がゆっくり流体中を動く場合は粘性抵抗が支配的になりkvモデルが適しています。一方、ボールや弾丸などが空気中を飛ぶ場合は空気の流れによる圧力抵抗が大きくなるため、kv²モデルの方が現実に近くなります。

速度に比例する抵抗の場合、運動方程式は線形微分方程式になります。そのため数学的に解きやすく、斜方投射でもx方向とy方向の運動を比較的簡単に求めることができます。

例えば水平方向では、抵抗力によって速度が指数関数的に減少する形になります。鉛直方向でも重力と抵抗の影響を組み合わせることで解析的な式を作ることができます。

しかしkv²モデルでは速度の大きさそのものが関係するため、水平速度と鉛直速度を合わせた速度ベクトルが非線形に影響します。そのため一般的な初等関数による解析解を求めることは困難になります。

速度比例モデルから二乗抵抗モデルに近づけるために、平均速度を利用して抵抗係数を調整する考え方は、近似手法としては意味があります。

例えばkv²の抵抗を一時的にkvの形に置き換える場合、ある代表速度v_aveを決めて、kv²=k’vとなるようにk’=kv_aveと置くことができます。このようにすれば、その速度付近では二乗抵抗の影響を近似できます。

ただし問題点として、斜方投射では速度は常に変化します。打ち出し直後は速度が大きく、頂点付近では鉛直速度が小さくなり、着地直前では再び速度が変化します。そのため一つの平均速度だけで全運動を表すことには限界があります。

kv²モデルをできるだけ正確に扱いたい場合は、微分方程式を数値的に解く方法が一般的です。

具体的には、時間を小さな区間に分割して、その瞬間の速度から空気抵抗を計算し、次の位置と速度を求める方法を使います。代表的な方法としてオイラー法やルンゲ・クッタ法があります。

例えば0.001秒ごとに計算すると、各時刻で速度ベクトルを求め、空気抵抗F=-kv²を計算できます。その結果、解析解がなくても実際の飛行軌道にかなり近い結果を得ることができます。

どうしても解析式に近い形で扱いたい場合は、速度区間ごとに抵抗係数を変える方法があります。

例えば飛行開始直後、中間、終盤で代表速度を設定し、それぞれ異なるk’を使う区分近似を行えば、単純な平均速度による方法より精度を上げることができます。

また、初期条件から試行計算を行い、得られた平均速度を使って再度係数を修正する反復計算を行う方法もあります。これは数値計算ほどではありませんが、手計算による近似としては有効です。

実際の工学分野では、空気抵抗が関係する飛翔体の解析では数値シミュレーションが広く利用されています。ゴルフボール、野球ボール、ロケットなども基本的には非線形な運動方程式をコンピューターで解いています。

そのため、kv²モデルに対して解析解がないことは大きな問題ではありません。重要なのは適切な物理モデルを選び、そのモデルを正しく数値的に解くことです。

平均速度を使ってkvモデルへ変換する方法は、簡単な近似として利用できますが、精密な軌道を求める目的なら時間刻みを細かくした数値計算の方が信頼性があります。

空気抵抗ありの斜方投射では、kvモデルは計算しやすい一方で、現実の空気抵抗を完全には表せません。kv²モデルは現実に近いものの解析解が難しいため、通常は数値計算によって解きます。

平均速度から新しい抵抗係数を作る方法は近似としては成立しますが、速度変化が大きい運動では誤差が発生します。より厳密な結果を求めるなら、ルンゲ・クッタ法などによる数値解析を利用するのが最も一般的な方法です。