プリセット
機体パラメータ
発射・空力パラメータ
発射速度の推定
$$v_0 = \eta\sqrt{\frac{2 M_{cw} g h}{m}}$$
$$F_d = \frac{1}{2}\rho C_d A v^2$$
η: 機械効率, A: 断面積(球体)
投射物直径: 0.25m と仮定
CAE応用：弾道解析はロケット・砲弾・スポーツ球の飛行設計に直結。CFD連成で抗力を精密計算し、MBD（マルチボディダイナミクス）でてこ機構を解析します。
点線：理想放物線（空気抵抗なし）/ 実線：実際の弾道（空気抵抗あり）
投石機（トレビュシェ）弾道シミュレーターとは
🧑🎓
投石機って、どうやって石を遠くまで飛ばすんですか？単にてこで投げてるだけに見えるけど…
🎓
ざっくり言うと、重いカウンターウェイトが落ちる「位置エネルギー」を、軽い石を投げる「運動エネルギー」に一気に変換してるんだ。シミュレーターで「カウンターウェイト M」のスライダーを大きくしてみて。重くするほど、石の初速が上がって遠くに飛ぶのが体感できるよ。
🧑🎓
え、そうなんですか！でも、てこの長さも関係ありますよね？「長腕 L_long」と「短腕 L_short」って何が違うの？
🎓
良いところに気づいたね。短い方の腕に重り、長い方の先端に石をつるす。長腕を長くすると、重りの小さな動きが石の大きな動きに変わる「速度比」が上がるんだ。実務では、石を速く飛ばすために長腕を長く設計するけど、その分構造が弱くなるトレードオフがある。実際にパラメータをいじって、射程がどう変わるか確かめてみよう。
🧑🎓
なるほど！でも、空気抵抗の「抗力係数 Cd」を変えると、飛び方が全然変わりますね。教科書だと45度で一番遠くに飛ぶって習ったけど…
🎓
その通り。空気抵抗がゼロ（真空中）なら45度が最適だ。でも現実の空気中では、速く飛ぶほど抵抗が増えるから、少し低い角度の方が飛距離が出るんだ。例えば、Cdを0.1から0.5に上げて、発射角を45度と35度で比べてみて。空気抵抗の影響が大きいと、最適角度が下がるのがよくわかるよ。
物理モデルと主要な数式
カウンターウェイトの位置エネルギーが、投射物の運動エネルギーに変換されると仮定して、発射速度を求めます。機械効率ηでエネルギーロスを考慮しています。
$$v_0 = \eta\sqrt{\frac{2 M_{cw} g h}{m}}$$
$v_0$: 投射物の発射速度 [m/s], $\eta$: 機械効率 (0~1), $M_{cw}$: カウンターウェイト質量 [kg], $g$: 重力加速度 [m/s²], $h$: カウンターウェイトの落下高さ [m], $m$: 投射物質量 [kg]
飛翔中の投射物には、速度の2乗に比例する空気抵抗（抗力）が働きます。これが弾道を理想的な放物線からずらし、飛距離を減少させます。
$$F_d = \frac{1}{2}\rho C_d A v^2$$
$F_d$: 空気抵抗（抗力） [N], $\rho$: 空気密度 [kg/m³], $C_d$: 抗力係数 (無次元), $A$: 投射物の進行方向への断面積 [m²], $v$: 投射物の速度 [m/s]
実世界での応用
兵器・防衛分野：砲弾やロケット弾の弾道計算に直接応用されます。発射条件と空気抵抗を精密にシミュレーションすることで、命中精度を向上させます。
スポーツ工学：野球のボールやゴルフボールの飛翔解析に使われます。ボールの回転による揚力（マグヌス効果）を加えた拡張モデルもあり、変化球の研究などに活用されています。
宇宙工学：ロケットの打ち上げや、大気圏再突入カプセルの軌道計算の基礎となります。高高度では空気密度が変わるため、高度に依存するρのモデルが使われます。
CAE（コンピュータ支援工学）：このシミュレーターの原理は、MBD（マルチボディダイナミクス）で投石機の機構解析をし、得られた発射条件をCFD（数値流体力学）で詳細な空力解析に連成させる、といった本格的なCAEワークフローの入門として位置付けられます。
よくある誤解と注意点
まず、「カウンターウェイトを重くすればするほど飛距離が伸びる」と思いがちですが、そう単純ではありません。確かに初速は上がりますが、同時に投石機の構造にかかる負荷が増大します。例えば、カウンターウェイトを100kgから200kgに倍増させると、発射時の軸受け部分のトルクは単純計算で約2倍に。実物では部材の破損や、腕の振り始めに大きな慣性抵抗が生まれ、逆に効率（η）が低下する可能性があります。シミュレーターでは「機械効率η」を一定と仮定していますが、実機設計ではパラメータ変更に伴う効率の変化を常に考慮する必要があります。
次に、長腕を極端に長く設定すると射程が最大化されると考える点。理論上は速度比が上がりますが、腕自体の質量と慣性モーメントが無視できなくなります。重い腕を加速するのにエネルギーを消費するため、石に伝わるエネルギーが減少します。例えば、長腕を5mから10mにすると、速度比は2倍になりますが、腕の質量を考慮したシミュレーションでは最適射程を与える長さが存在することがわかります。これは「等価質量」の概念で理解できます。
最後に、抗力係数Cdを現実の値に設定する際の落とし穴。Cdは形状に強く依存します。このシミュレーターでは「球」を仮定していますが、実際の中世の石弾は不規則な形状のため、Cdは0.4〜0.6程度と、滑らかな球（〜0.1）より大幅に高くなります。また、飛翔中に回転やゆらぎが生じるとCdは一定ではなくなります。実務的なCFD解析では、この「一つの値で代表させることの難しさ」を認識することが第一歩です。
関連する工学分野
このシミュレーターの核となる計算は、MBD（マルチボディダイナミクス）そのものです。投石機の腕、支点、カウンターウェイト、石を剛体リンクと見なし、それらの間の相対運動と力を計算します。発射瞬間の複雑な力の伝達（ショックロード）を解析するには、こちらの手法が必須です。例えば、自動車のサスペンションやロボットアームの動作解析も、同じMBDの枠組みで行われています。
弾道計算部分は、飛行力学や誘導制御工学の基礎と直結します。ミサイルやドローンの飛翔軌道予測も、質量、初速、発射角、空気抵抗・揚力のモデルをより精密にした上で、同じ運動方程式を解いています。特に「空気抵抗が最適発射角を45度からずらす」という現象は、長距離砲弾やゴルフボールの弾道設計でも重要な考慮事項です。
また、エネルギー変換の考え方は機構学や機械設計の基本です。カウンターウェイトの位置エネルギーが、どのように損失（摩擦、振動、音、熱）を経て投射物の運動エネルギーになるかを追うことは、あらゆる機械システムの効率改善に通じます。例えば、油圧ショベルのブームの動きや、振り子式発電装置の設計でも、同様のエネルギー収支計算が行われています。
発展的な学習のために
まず次のステップとしておすすめなのは、「パラメータスタディ」を体系化することです。例えば「最大射程を達成する発射角θと、カウンターウェイト質量Mの関係」をグラフにプロットしてみてください。空気抵抗Cdを変えると、その曲線がどう変化するか。これにより、複数パラメータが互いにどう影響し合うか（相互作用）の直感が養えます。これは実験計画法（DOE）の入門的な実践です。
数学的背景を深めたいなら、シミュレーターの裏で解かれている常微分方程式の数値解法（オイラー法やルンゲ・クッタ法）を学びましょう。このツールは、運動方程式 $m \frac{d^2 \vec{r}}{dt^2} = \vec{F}_g + \vec{F}_d$ を微小時間刻みで逐次計算しています。これを理解すれば、自分で簡単な弾道プログラムをスクラッチで書けるようになります。
さらに学ぶべきトピックは、「連成解析」の概念です。現実の設計では、まずMBDツールで発射時の初速と角度を求め、その結果をCFDツールに渡して詳細な空力解析を行い、得られた抗力・揚力を再びMBDにフィードバックする、というループを回します。この投石機シミュレーターは、そのプロセスを「MBD（エネルギー変換）→ 飛翔力学（弾道計算）」という形で非常に簡略化して体験できる、優れた入門教材なのです。