パラメータ
プリセット
理論・主要公式
\(m\ddot{x} = mg - k(x - L_0) - b\dot{x}\)
コード弛緩時（\(x \le L_0\)）
\(m\ddot{x} = mg - b\dot{x}\)
最低点近似（\(b=0\)）:
\(x_{max} = L_0 + \dfrac{mg}{k} + \sqrt{\left(\dfrac{mg}{k}\right)^2 + \dfrac{2mgL_0}{k}}\)
バンジーコードの物理 — 弾性ポテンシャルと自由落下
🙋
先生、バンジージャンプって跳んでからコードが張るまでは普通に自由落下してるんですよね？コードが伸び始めたら何が起きるんですか？
🎓
そう、最初はただの自由落下だ。コードが自然長 \(L_0\) を超えて伸び始めると、バネの力 \(F = k(x - L_0)\) が上向きに働く。運動方程式は \(m\ddot{x} = mg - k(x - L_0) - b\dot{x}\)（下向きを正）になる。\(b\dot{x}\) は空気抵抗や内部損失による減衰力だ。
🙋
なるほど。最低点ではどうなってるんですか？速度がゼロになる瞬間ですよね？
🎓
その通り。最低点では速度ゼロだから、エネルギー保存則で求めると \(x_{max} = L_0 + mg/k + \sqrt{(mg/k)^2 + 2mgL_0/k}\) になる。例えばm=70kg、k=40 N/m、L₀=20mなら、\(mg/k = 70 \times 9.81 / 40 \approx 17.2\ \text{m}\)、\(x_{max} \approx 20 + 17.2 + 21.3 \approx 58.5\ \text{m}\)。コードが約3倍に伸びるんだよ。
🙋
3倍も！最低点での引っ張りGってどのくらいかかるんですか？
🎓
最低点での加速度はコード力が最大の瞬間。\(a = k(x_{max} - L_0)/m - g\) で、上の例では約\((40 \times 38.5)/70 - 9.81 \approx 22 - 9.81 \approx 12.2\ \text{m/s}^2 \approx 1.2G\)程度。バネが硬いほどGが増え、一般的なバンジーでは2〜3Gが普通だ。乗用車のブレーキが約0.8G、戦闘機のパイロットが9Gまで耐えることを考えると、バンジーのGは実はそれほど大きくないよ。
🙋
思ったより低いんですね。減衰がゼロだと永遠に揺れ続けるんですか？実際はどこでエネルギーが失われるんでしょう？
🎓
理想的には永遠に振動する。実際には①コードの内部摩擦（ゴム・ラテックスの粘弾性）、②空気抵抗、③コードを結んでいる点での微小エネルギー散逸でダンピングされる。シミュレーターで減衰係数bをゼロにすると振動が続くのを確認できるよ。また「エネルギー変換タブ」で位置エネルギー・運動エネルギー・弾性エネルギーの移り変わりも見られる。
よくある質問
バンジーの最低到達点はどうやって計算するの？
減衰なしの場合、エネルギー保存則を使います。ジャンプ台（基準点）から最低点 \(x_{max}\) まで落下したとき、重力ポテンシャルエネルギーの減少 = コードの弾性ポテンシャルエネルギーの増加：
\(mgx_{max} = \frac{1}{2}k(x_{max} - L_0)^2\)
この二次方程式を解くと \(x_{max} = L_0 + mg/k + \sqrt{(mg/k)^2 + 2mgL_0/k}\) が得られます。
バネ定数kが大きいとどうなるの？
コードが「硬い」ほど、①最低到達点が浅くなる（コードがあまり伸びない）、②振動周期 \(T \approx 2\pi\sqrt{m/k}\) が短くなる（速く往復する）、③最低点でのGが増える。逆にkが小さいほどコードはよく伸び、最低点が深くなって地面への衝突リスクが上がります。バンジー設計ではジャンプする人の体重に合わせてkを調整することが安全の要です。
減衰係数の効果は？
減衰係数 \(b\)（単位 N·s/m）が大きいほど振動が早く収まります。\(b = 0\) は完全弾性（永続振動）、\(b = 2\sqrt{mk}\) で臨界制動（一度も振り返らずゆっくり平衡点に落ち着く）、それ以上は過制動です。実際のバンジーコード（天然ゴム系）の等価減衰係数は数 N·s/m 程度です。空気抵抗の寄与は速度の2乗に比例しますが、低速域での近似として \(b\dot{x}\) で近似できます。
実際のバンジーコードはどんな素材？
天然ラテックスゴムや合成ゴムのロープを束ねたもので、自然長の2〜3.5倍まで伸縮できます。フックの法則（線形バネ）から外れた非線形弾性を持ちますが、安全設計では線形近似でマージンを確保します。安全係数は通常3〜5倍（コードの破断強度 / 設計最大荷重）。バネ定数は体重65kgの人でおおよそ30〜60 N/m程度ですが、ジャンプ台の高さと最低許容高さによって設計されます。
ルンゲ-クッタ法とは？なぜEuler法より良いの？
ルンゲ-クッタ4次法（RK4）は1ステップにつき4回の傾き評価を行い、その加重平均で次の値を推定します。精度はステップ幅 \(h\) の5乗で向上（誤差 \(O(h^5)\)）。単純なEuler法では誤差 \(O(h^2)\) で、バンジーのような非線形振動ではエネルギーが増え続ける不安定解が出ることがあります。RK4は同じステップ数でEulerより数桁精度が高く、物理シミュレーションの標準的な手法です。
バンジージャンプ物理シミュレーターとは
バンジージャンプ物理シミュレーターの物理モデルでは、ジャンパーを質量\(m\)の質点とみなし、コードをフックの法則に従う線形ばねと減衰要素の並列結合としてモデル化する。自然長\(L_0\)のコードが伸びた際の復元力は\(F_{\text{spring}} = -k(x - L_0)\)、減衰力は\(F_{\text{damping}} = -c v\)で与えられる。ここで\(x\)は基準点からの変位、\(v\)は速度、\(k\)はばね定数、\(c\)は減衰係数である。重力\(mg\)を加味した運動方程式は\(m \frac{d^2x}{dt^2} = mg - k(x - L_0) - c v\)となり、この2階常微分方程式を1階連立系に変換し、4次のルンゲ-クッタ法で数値積分する。時間刻み\(\Delta t\)を十分小さく設定することで、コードがたるんでいる自由落下区間（\(x < L_0\)では\(k=0\)）と伸張区間を統一的に扱う。シミュレーションは最低到達点\(x_{\min}\)、振動周期\(T\)、位置エネルギー\(U = mgx\)と運動エネルギー\(K = \frac{1}{2}mv^2\)および弾性エネルギー\(E_{\text{elastic}} = \frac{1}{2}k(x-L_0)^2\)の総和が減衰により減少する様子をリアルタイムで可視化する。
実世界での応用
産業での実際の使用例
本シミュレーターは、アミューズメント業界のバンジージャンプ施設設計において、安全マージン検証に活用されています。例えば、大手遊具メーカー「サノヤス・ライド」が、新規タワー型アトラクションの開発時に、参加者の体重分布やコードの経年劣化を想定した動的負荷解析を実施。ルンゲ-クッタ法による高精度な数値積分で、最低到達点での最大張力や減衰特性を予測し、安全基準（ISO 10356）への適合確認に利用しています。
研究・教育での活用
大学の物理工学科では、非線形振動の教材として採用。例えば東京工業大学の「機械力学」講義では、学生が質量や減衰係数を変化させてエネルギー保存則の破綻を観察し、現実の空気抵抗やコードの非弾性散逸を理解します。また、制御工学の研究室では、アクティブダンパー設計の前段階として、本シミュレーターで得た振動周期データをフィードバック制御のモデル検証に使用しています。
CAE解析との連携や実務での位置付け
実務では、本ツールは簡易動的解析の「一次スクリーニング」として位置づけられます。例えば、建設会社が橋梁点検用ゴンドラのワイヤー安全設計を行う際、まず本シミュレーターでパラメータスタディを実施し、危険領域を特定。その後、詳細な有限要素法（FEM）解析ソフト「ANSYS Mechanical」に受け渡し、ワイヤー内部応力分布や疲労寿命を精査するワークフローが一般的です。これにより、計算コストを削減しつつ、設計初期段階でのリスク低減を実現しています。
よくある誤解と注意点
「減衰係数を大きくすれば落下が穏やかになる」と思いがちですが、実際は減衰が強すぎるとバンジーコードが伸び切る前に急激に減速され、逆にジャンパーに大きな衝撃（加速度）が生じる可能性があります。適切な減衰は振動を抑えつつ、急激な制動を避けるバランスが重要です。
「バネ定数が大きいほど最低到達点が低くなる」と考えられがちですが、実際にはバネ定数が大きいコードは伸びにくく、エネルギーを吸収しづらいため、むしろ最低到達点が高くなることがあります。最低到達点は質量・コード長・バネ定数の相互作用で決まるため、単純な比例関係ではない点に注意が必要です。
「ルンゲ-クッタ法を使えば常に正確な解が得られる」と思いがちですが、実際には時間刻み幅が大きすぎると数値誤差が蓄積し、特に減衰が小さい場合に振動の周期や振幅がずれることがあります。シミュレーションの精度を保つには、適切な刻み幅の設定と結果の物理的な妥当性の確認が欠かせません。