🙋
空気抵抗があると、ボールの飛び方はどう変わるんですか？教科書の放物運動と全然違うって聞いたけど。
🎓
ざっくり言うと、飛距離が大幅に短くなって、軌道が非対称になるんだ。真空だと綺麗な放物線だけど、空気抵抗があると、上がり方はゆるやかで、下り方は急な軌道になるよ。このシミュレーターで、初速と角度はそのままに「抗力係数」のスライダーを0から0.5くらいに動かしてみて。飛距離がどんどん縮むのがわかるはずだ。
🙋
え、すごく縮みました！抗力係数って何ですか？ボールの形で決まるんですか？
🎓
その通り。抗力係数 $C_d$ は物体の形で決まる抵抗の大きさを表す無次元数だよ。例えば、丸い球は $C_d \approx 0.47$、流線形は0.04くらい。ゴルフボールの表面にあるディンプル（小さなくぼみ）は、この $C_d$ を約0.25まで下げて飛距離を伸ばすための工夫なんだ。シミュレーターで「直径」も変えてみると、大きいボールほど抵抗が大きくなってすぐ落ちるのが体感できるよ。
🙋
なるほど！じゃあ、このシミュレーターはどうやって計算してるんですか？空気抵抗があると数式が解けないって聞いたけど。
🎓
鋭いね。空気抵抗が速度の2乗に比例すると、運動方程式が非線形になって、紙と鉛筆では解けない（解析解がない）んだ。そこで使うのが「RK4（ルンゲ・クッタ4次法）」という数値積分。CAEソフトでも時間経過を追う計算の基本だよ。このシミュレーターも、微小な時間刻み（0.01秒）ごとに速度と位置を更新して、軌道を描いている。パラメータを変えると、リアルタイムで飛距離や最高高度が計算される仕組みだ。
空気抵抗（抗力）は速度の2乗に比例して大きくなり、常に進行方向と逆向きに働きます。そのため、特に高速域で減速が大きく、到達距離が短くなったり、頂点が低くなったりします。真空軌道は放物線ですが、空気抵抗があると軌道が非対称になります。
抗力係数や直径を大きくすると、空気抵抗が増加し、軌道がより急激に減速・落下します。逆に小さくすると真空軌道に近づきます。質量を大きくすると慣性が増し、空気抵抗の影響が相対的に小さくなります。
RK4（4次ルンゲ・クッタ法）は、運動方程式を高精度で近似計算する手法です。1ステップごとに4回の評価を行い、誤差が非常に小さいため、このシミュレーションでは実用的に正確な軌道が得られます。刻み幅は内部で自動調整されています。
基本的な物理モデルは正しいですが、実際の弾道では風、回転（マグヌス効果）、空気密度の変化なども影響します。本ツールはCAE数値解析の入門用であり、厳密な設計にはより高度なモデルが必要です。ただし、傾向を掴むには十分な精度です。
スポーツ工学：ゴルフボールのディンプルやサッカーボールの縫い目は、抗力係数 $C_d$ を最適化して飛距離や軌道の安定性を向上させるために設計されています。シミュレーションでパラメータを変えると、その効果が定量的に理解できます。
砲弾・ロケットの弾道計算：実戦での砲弾やロケットの軌道は、空気抵抗を無視できません。発射角度や初速だけでなく、弾頭形状（$C_d$）や質量を考慮した精密な数値シミュレーションが命中精度を決定します。
ドローン・パラシュートの降下制御：ドローンの飛行やパラシュートによる降下は、空気抵抗と重力のバランスで成り立っています。安全な着地速度や飛行時間を予測するために、同様の運動方程式が用いられます。
CAE数値解析の入門教育：この「空気抵抗つき放物運動」は、非線形常微分方程式を数値的に解くという、構造動解析や流体動解析など本格的なCAEの根幹にあるプロセスを学ぶ最適な題材です。
このシミュレーターを使い始める際、特にCAE初心者が陥りがちなポイントがいくつかあります。まず「抗力係数は固定値ではない」という点。ツールでは定数として設定していますが、実は物体の形状や速度（より正確にはレイノルズ数）によって変化します。例えば、球の抗力係数は低レイノルズ数域では高く、ある領域を超えると急降下し、ほぼ一定（約0.47）になります。つまり、初速が大きく変わると、設定した$C_d$値自体が現実からずれる可能性があるんです。
次に、パラメータの無次元化を意識しよう。直径を2倍にすると、投影面積は4倍になりますが、質量は体積に比例するので8倍になります（密度が同じ場合）。つまり、直径だけを2倍にすると「質量の増加による慣性の効果」と「面積増加による抵抗の効果」が競合し、直感に反する結果が出ることも。例えば、直径2cmと4cmの鉄球を同じ条件で投げると、大きい方が意外と遠くまで飛ぶかもしれません。パラメータを変える時は、個々の数値ではなく、現象を支配する無次元数（例えば、抗力と重力の比）がどう変わるかを考えるクセをつけましょう。
最後に、「RK4は万能ではない」という落とし穴。このツールで使っているRK4法は高精度ですが、時間刻み$\Delta t$を粗くしすぎると計算が発散したり、誤差が大きくなったりします。逆に細かすぎると計算コストが無駄に増えます。実務では、現象の時間スケール（例えばボールが頂点に達する時間）に対して、$\Delta t$を適切に設定することが重要です。このツールの0.01秒は、多くのケースでバランスが取れていますが、極端に速いまたは遅い現象をシミュレートする際は要注意です。