パラメータ
\(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}}\)
質量・コード長・バネ定数・減衰係数を動かすたびに、バンジーの自由落下と減衰振動がルンゲ-クッタ法で解き直されます。最低到達点・振動周期に加え、運動/位置/弾性エネルギーの推移までグラフで追えます。
バンジージャンプ物理シミュレーターの物理モデルでは、ジャンパーを質量\(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」に受け渡し、ワイヤー内部応力分布や疲労寿命を精査するワークフローが一般的です。これにより、計算コストを削減しつつ、設計初期段階でのリスク低減を実現しています。
「減衰係数を大きくすれば落下が穏やかになる」と思いがちですが、実際は減衰が強すぎるとバンジーコードが伸び切る前に急激に減速され、逆にジャンパーに大きな衝撃(加速度)が生じる可能性があります。適切な減衰は振動を抑えつつ、急激な制動を避けるバランスが重要です。
「バネ定数が大きいほど最低到達点が低くなる」と考えられがちですが、実際にはバネ定数が大きいコードは伸びにくく、エネルギーを吸収しづらいため、むしろ最低到達点が高くなることがあります。最低到達点は質量・コード長・バネ定数の相互作用で決まるため、単純な比例関係ではない点に注意が必要です。
「ルンゲ-クッタ法を使えば常に正確な解が得られる」と思いがちですが、実際には時間刻み幅が大きすぎると数値誤差が蓄積し、特に減衰が小さい場合に振動の周期や振幅がずれることがあります。シミュレーションの精度を保つには、適切な刻み幅の設定と結果の物理的な妥当性の確認が欠かせません。
体重80kg、コード長25m、バネ定数60N/mの設定例:自由落下区間(0~25m)では加速度9.81m/s²で、自然長25m到達時の速度は22.1m/sに達します。減衰なしの場合、最低到達点は x_max = 25 + 80·9.81/60 + √((80·9.81/60)² + 2·80·9.81·25/60) ≈ 66.8m で、コードの伸びは約41.8mになります。最低点での最大加速度は約2.2gです。最低点では運動エネルギーが弾性ポテンシャル E = ½·60·41.8² ≈ 52.4kJ に変換され、これは重力位置エネルギー mg·x_max ≈ 52.4kJ と一致します(エネルギー保存)。減衰係数を加えると振動は数周期で減衰し、最低到達点も浅くなります。