パラメータ設定
Duffing方程式：
$\ddot{x}+ 2\zeta\omega_0\dot{x} + \omega_0^2 x + \varepsilon x^3 = F\cos(\Omega t)$
▲ 非線形周波数応答曲線（FRC）と骨格曲線（Backbone）
▲ 非線形ポテンシャル V(x) = ½ω₀²x² + ¼εx⁴
理論・主要公式
一次近似の振幅方程式（$x \approx A\cos\Omega t$）：
$$\left[(\omega_0^2 + \tfrac{3}{4}\varepsilon A^2 - \Omega^2)^2 + (2\zeta\omega_0\Omega)^2\right]A^2 = F^2$$
骨格曲線（減衰ゼロの自由振動）：
$$\Omega_{bb}= \omega_0\sqrt{1 + \frac{3\varepsilon A^2}{4\omega_0^2}}$$
ジャンプ現象は3値解の折り畳み点（Fold）で発生
非線形振動（Duffing振動子）とは
🙋
「非線形振動」って、普通のバネの振動と何が違うんですか？
🎓
大まかに言うと、バネの硬さが伸び縮みに応じて変わる振動だね。普通の線形バネは「力 = 剛性 × 変位」で一定だけど、非線形バネは変位が大きくなるほど硬くなったり（ハードニング）、逆に柔らかくなったり（ソフトニング）するんだ。このシミュレーターの「非線形係数 ε」をプラスやマイナスに変えてグラフを見てみると、その違いが一目瞭然だよ。
🙋
え、バネの硬さが変わるって、そんなことあるんですか？例えばどんなものに使われてるんですか？
🎓
実務では結構あるよ。例えば自動車のエンジンマウントや建物の免震ゴムベアリング。大きな力がかかるとゴムが硬化して、それ以上変形しにくくなる性質（ハードニング）を利用してるんだ。シミュレーターで「ε」を大きくして「励振力振幅 F」も大きくすると、応答曲線が右に曲がるのがわかる。これがハードニング効果だね。
🙋
「ジャンプ現象」って聞きました。周波数をゆっくり変えるだけで、急に振幅が跳ね上がるって本当ですか？
🎓
本当だよ。これが非線形振動の面白いところで、設計では注意が必要な現象なんだ。上のスライダーで「周波数比 Ω/ω₀」をゆっくり増加させていってみて。あるところで、グラフ上の点が突然、上の枝にジャンプするのが観測できるはず。逆に周波数を下げていくと、今度は下にジャンプする。この不連続な変化をCAEで予測するのは結構難しいんだ。
よくある質問
これはDuffing振動子の非線形性（εx³項）によるジャンプ現象です。ハードニングばね（ε>0）では高周波側に、ソフトニングばね（ε<0）では低周波側に曲線が傾き、同一周波数で3つの振幅解が生じます。中央の解は不安定で実際には観測されず、周波数を掃引すると振幅が不連続に跳びます。
骨格曲線は減衰と励振がない場合（ζ=0, F=0）の自由振動における振幅と周波数の関係です。非線形性により振幅が大きくなると固有振動数が変化する様子を示し、ハードニングでは高周波側、ソフトニングでは低周波側に曲がります。周波数応答曲線のピークの軌跡と一致します。
振動応答が励振と同じ角周波数Ωの単一の余弦波（x ≈ A cosΩt）で近似できると仮定しています。高調波成分（3Ω, 5Ωなど）を無視するため、非線形性が強い場合や共振付近では誤差が生じることがありますが、基本的なジャンプ現象や骨格曲線の傾向を捉えるには十分実用的です。
減衰比ζを大きくすると共振ピークが低くなり、ジャンプ現象が起こりにくくなります。非線形係数εの絶対値を大きくすると曲線の傾きが強くなり、ジャンプの周波数幅が拡大します。ε=0では線形の共振曲線に戻ります。これらのパラメータを調整して実際の振動系の特性を再現できます。
実世界での応用
自動車・機械の防振：エンジンマウントやパワートレーンマウントにはゴムが使われ、大振幅時に硬化（ハードニング）する特性を持ちます。これをDuffingモデルで近似し、特定の回転数で振幅が急増する「ジャンプ」を回避する設計が行われます。
建築・土木の免震構造：建物と基礎の間に設置する積層ゴムベアリングは、大地震時の大きな変位で剛性が増す（ハードニング）ように設計され、変位を抑制します。非線形応答解析は安全設計に不可欠です。
MEMS/NEMSデバイス：微細な機械素子はしばしば大きな変形を示し、強い幾何学的非線形性（Duffing型）を発現します。共振周波数の振幅依存性を理解し、センサーの感度や動作範囲を最適化するために利用されます。
CAE解析（FEM）：ANSYSやAbaqusなどの非線形過渡応答解析では、材料や幾何学の非線形性からDuffing型の応答が現れます。特にゴム材料の超弾性モデルと組み合わせた振動解析は、実機試験前の重要な検証ステップです。
よくある誤解と注意点
まず、「非線形係数εが大きいほど必ず危険」というのは誤解です。確かにεがプラス（ハードニング）で大きいと、共振点が高周波側にシフトし、想定外の回転数で大きな振動が起こる可能性があります。しかし、この特性を逆手に取れば、ある周波数帯域の振動を意図的に抑制する設計も可能です。例えば、ε=0.5、F=0.3の設定で周波数を掃引すると、振幅のピークが右にずれて幅も狭まるのが確認できます。これは、ある意味で「振動しにくい領域」を広げているとも言えます。
次に、減衰比ζの設定を軽視しがちです。線形振動ではζが小さいと共振ピークが鋭く高くなるだけですが、非線形ではジャンプ現象の発生範囲や、安定解が存在する周波数幅に直接影響します。ζを0.01から0.05に少し増やすだけで、ジャンプが起こる周波数領域が大きく変わったり、場合によってはジャンプそのものが消失したりします。実務では、この減衰の値が材料や構造でばらつくため、シミュレーション結果を過信せず安全マージンを考えることが重要です。
最後に、「骨格曲線」と「周波数応答曲線」を混同しないこと。骨格曲線（$$ \omega_0^2 + \tfrac{3}{4}\varepsilon A^2 $$で表されるピークの軌跡）は、あくまで減衰ゼロの理想的な共振周波数の変化を示します。実際の応答曲線は、この骨格曲線を中心に、減衰によって幅を持ち、かつジャンプ現象で「折れ曲がった」形になります。設計では骨格曲線の傾向をまず掴み、その後に減衰を加えた実際の応答を評価する、という二段階の思考が役立ちます。