単振動・減衰振動・強制振動とは
🧑🎓
単振動と減衰振動、強制振動って何が違うんですか？グラフを見ると、減衰振動だけ振幅がだんだん小さくなってますね。
🎓
ざっくり言うと、エネルギーが保存されるか、失われるか、外から与えられるかの違いだね。単振動は摩擦や抵抗がない理想的な振動で、バネに重りをつけて宇宙空間で振らせたイメージ。減衰振動は現実世界の振動で、空気抵抗や摩擦でエネルギーが失われて振幅が減衰する。シミュレーターの「減衰係数γ」のスライダーを大きくすると、減衰が早くなるのがわかるよ。
🧑🎓
え、そうなんですか！じゃあ強制振動は？「外力振幅F₀」ってパラメータがありますが、これは何を表してるんですか？
🎓
強制振動は、ブランコを外から定期的に押して動かし続けるようなものだ。F₀はその「押す力の強さ」を表しているんだ。面白いのは「強制角周波数ω_d」を変えたときの挙動だ。上のスライダーでω_dを固有角周波数ω₀に近づけてみて。振幅が急に大きくなる「共鳴」という現象が起こるよ。実務では、この共鳴を避ける設計がすごく重要になるんだ。
🧑🎓
共鳴って、橋が壊れるあの現象ですよね！「Q値」ってのも出てきますが、これは何の指標なんですか？
🎓
その通り！Q値は振動系の「鋭さ」や「持続性」を表す品質係数だ。数式では$Q = \frac{\omega_0}{2\gamma}$で、減衰係数γが小さい（=摩擦が少ない）ほどQ値は大きくなる。シミュレーターでγを小さくして強制振動の共鳴ピークを見てごらん。Q値が大きいと、共鳴時のピークが鋭く高くなるし、減衰振動も長く続く。例えば高音質スピーカーはQ値を調整して、特定の周波数をきれいに響かせているんだ。
物理モデルと主要な数式
このシミュレーターの根幹をなすのは、減衰と外力を含む運動方程式（強制振動の式）です。単振動と減衰振動はこの式の特別な場合として導かれます。
$$m\frac{d^2x}{dt^2}+ c\frac{dx}{dt}+ kx = F_0 \cos(\omega_d t)$$
ここで、$m$: 質量, $c$: 減衰係数, $k$: ばね定数, $F_0$: 外力振幅, $\omega_d$: 強制角周波数です。両辺を$m$で割り、$\omega_0 = \sqrt{k/m}$（固有角周波数）, $\gamma = c/(2m)$（減衰係数）と置き換えると、より見慣れた形になります。
強制振動の定常状態における振幅$X$は、以下の式で与えられます。これが共鳴曲線を決定する式です。
$$X=\frac{F_0/m}{\sqrt{(\omega_0^2-\omega_d^2)^2+(2\gamma\omega_d)^2}}$$
$\omega_d = \omega_0$のとき分母の第一項が0になり、振幅$X$は最大値$X_{max}= \frac{F_0/m}{2\gamma\omega_0}$に達します（共鳴）。減衰係数$\gamma$が小さいほど、このピークは鋭く高くなります。
実世界での応用
機械・構造物の設計：橋梁、ビル、航空機の翼、自動車のシャシなど、あらゆる構造物には固有振動数があります。エンジン振動や風による周期的な力（外力）がこの周波数に一致すると共鳴が起き、大振幅振動による疲労破壊や騒音の原因となります。CAEシミュレーションで固有値を事前に解析し、共鳴を避ける設計が必須です。
電子回路・通信機器：LC共振回路や水晶振動子は、電気的な振動系です。特定の周波数の信号だけを選択的に増幅・発振させるために用いられ、その選択性の鋭さはQ値で評価されます。ラジオのチューニングや時計の基準信号生成はこの原理の応用です。
音響機器：スピーカーや楽器は、機械的振動を音に変換する装置です。スピーカーのエンクロージャー（箱）やダイアフラム（振動板）の設計では、目的の周波数帯域で望ましい振動特性（Q値）を得ることが重要です。低音を力強く響かせるためには、適度な減衰（ダンピング）が必要です。
計測・センシング：原子間力顕微鏡（AFM）のカンチレバーやMEMS加速度センサーは、微小な質量とばねからなる振動子です。外部からの微小な力や質量変化を、振動子の共振周波数のシフトとして検出することで、ナノスケールの計測を実現しています。
よくある誤解と注意点
まず、「減衰係数γ」と「Q値」の関係を逆に覚えてしまうケースが多いよ。γは「振動を減らす係数」だから、大きくすればするほど振動は早く収まる。一方、Q値はその逆数みたいなもので、γが小さい（減衰が少ない）ほどQ値は大きくなる。例えば、γを0.1から0.01に小さくすると、Q値は約5倍になり、共鳴ピークは鋭く高くなるんだ。実務で「Q値が高いシステム」と言われたら、「減衰が小さくて、特定の周波数に敏感に反応するシステム」と読み替えよう。
次に、「固有角周波数ω₀」はパラメータを変えても変わらないと思いがちな点。このシミュレーターでは、ω₀は「質量m」と「ばね定数k」から決まる値だ。だから、強制振動の周波数ω_dを変えてもω₀自体は変わらない。しかし実世界では、大きな振幅が生じると材料が変形してkが変わったり（非線形性）、質量の分布が変わって実効的なmが変化したりする。つまり、共鳴点そのものが動くことがあるんだ。シミュレーターは線形モデルだから、そこは頭の片隅に入れておこう。
最後に、強制振動の「定常状態」をすぐに見ようとしないこと。外力を加え始めた直後は「過渡状態」で、振動が不安定だ。特に減衰が小さい（Q値が高い）系では、定常状態に落ち着くまでに非常に長い時間がかかる。例えば、γ=0.05、ω_d=ω₀の設定でシミュレーションをスタートさせると、振幅が最大値に達するまでに数十周期かかることもある。実験データと比較する時は、必ず十分な時間が経過した「定常状態」の波形を見るようにしよう。
関連する工学分野
このシミュレーターの背後にある「線形振動理論」は、音響工学の基礎そのものだ。スピーカーのコーン紙やヘッドホンのダイアフラムは、まさに質量（m）とバネ（復元力k）、ダンパー（減衰c）で構成される振動系だ。Q値の調整は音質設計の肝で、低音用スピーカーはQ値を適度に下げてダンピングを効かせ、ブーミーな音を防いでいる。
制御工学とも深く結びついている。サーボモータの位置決め制御では、指令値への追従応答がまさに減衰振動の形を示す。減衰係数γが小さすぎると（=Q値が高すぎると）、目標位置を行き過ぎたり戻ったりを繰り返す（オーバーシュート）。逆に大きすぎると、目標に到達するまでに時間がかかりすぎる。最適な応答を得るために、PID制御のパラメータ調整は、この振動系の減衰特性を理想的な形に整える作業なんだ。
さらに地震工学では、建物を一種の振動系（特に減衰振動系）とみなす。地盤から入力される地震動が「外力」に相当する。建物の固有周期（=2π/ω₀）と地震動の卓越周期が一致すると共鳴（共振）し、大きな被害をもたらす。現代の免震構造は、この固有周期を意図的に長くしたり、減衰係数を大幅に増加（ダンパー設置）させたりして、地震エネルギーを吸収するように設計されているんだ。
発展的な学習のために
まず次のステップは、「2自由度振動」を学ぶことだ。このシミュレーターは質量が1つの「1自由度系」。現実の多くの構造物は複数の質量とバネが連なった多自由度系だ。2自由度系を学ぶと、新しい重要な概念「連成振動」と「モード」が登場する。例えば、車のサスペンション（車体とタイアという2つの質量）では、2つの異なる固有振動数が生じ、振動のエネルギーが片方からもう片方に移動する現象を観察できる。
数学的には、微分方程式の解法をしっかり理解することが全ての基礎になる。特に「過渡解」と「定常解」の和としての一般解、そしてその中で時間とともに消えていく過渡解の部分の意味を掴もう。また、共鳴曲線の式 $X(\omega_d)$ を導出する過程で出てくる複素数を用いた解法（$e^{i\omega t}$を使う方法）は、電気回路の交流理論や波動工学にまで通じる強力なツールなので、ぜひマスターしたい。
最後に、この線形モデルを超えた「非線形振動」の世界に挑戦してみよう。現実のバネは大きく伸びるとフックの法則から外れ、減衰も速度に単純比例しないことが多い。非線形の世界では、振幅が跳び上がる「ジャンプ現象」や、外力周波数を変えていった時と戻した時で振る舞いが変わる「ヒステリシス」、さらには予測不能な「カオス振動」といった、豊かで複雑な現象が待っている。このシミュレーターで基礎を固めたら、ぜひその先の広大な振動学の世界を探検してみてほしい。