E = 200 GPa
ρ = 7850 kg/m³
L = 1.00 m
I = — m⁴
アニメーション中...
固有振動モード形シミュレーターとは
🧑🎓
固有振動の「モード形」って何ですか？教科書の図は静止してるけど、実際はどう動くんですか？
🎓
ざっくり言うと、構造物が「その周波数だけで揺れた時」の変形のカタチだよ。例えばギターの弦をはじくと、1次モードなら真ん中が一番大きく振れる半波の形で揺れるよね。このシミュレーターでは、上の「n」や「m」のスライダーを動かすと、そのモードのアニメーションがリアルタイムで見られるんだ。触ってみると、2次モードは節（動かない点）が1つ、3次は2つ…と増えていくのがわかるよ。
🧑🎓
え、そうなんですか！「単純支持」と「片持ち」で結果が全然違うみたいです。境界条件ってそんなに効くんですか？
🎓
めちゃくちゃ効くんだ。実務で多い片持ち梁は、固定端が動かないから、1次モードでも先端が一番大きく揺れる。シミュレーターで「境界条件」を切り替えて、同じ「n=1」でも形が全く違うことを確かめてみて。あと、材料の「E（ヤング率）」や寸法「L」を変えると、下に表示される「固有周波数」の値がガラッと変わる。剛い（Eが大きい）ほど高周波で振動するんだ。
🧑🎓
矩形板の振動も見られますね。梁と比べて何が難しいんですか？
🎓
良いところに気づいたね。板はx方向とy方向の振動が組み合わさるから、モード形が複雑になる。例えば「m=1, n=1」は全体がドーム状に膨らむけど、「m=2, n=1」はx方向に節が1本入る。シミュレーターで「m」と「n」をバラバラに変えて、どんなリッチな模様（モード）が現れるか試してみよう。FEM（有限要素法）で解析するときも、この2次元の振動をどう離散化するかが鍵になるんだ。
物理モデルと主要な数式
このシミュレーターの基礎は、オイラー・ベルヌーイ梁理論です。梁の曲げ振動を記述する4階の偏微分方程式は以下の通りです。
$$ \frac{\partial^2}{\partial x^2}\left( EI \frac{\partial^2 w}{\partial x^2}\right) + \rho A \frac{\partial^2 w}{\partial t^2}= 0 $$
ここで、$w(x,t)$は梁のたわみ、$E$はヤング率、$I$は断面二次モーメント、$\rho$は密度、$A$は断面積です。この式は、「曲げ剛性による復元力」と「質量による慣性力」のバランスを表しています。
上記の方程式を、変数分離法と境界条件（支持条件）を用いて解くことで、固有角周波数$\omega_n$とモード形$\varphi_n(x)$が得られます。例えば単純支持梁の場合、解は以下の美しい形になります。
$$ \varphi_n(x) = \sin\!\left(\dfrac{n\pi x}{L}\right), \quad \omega_n = \left(\dfrac{n\pi}{L}\right)^{\!2}\!\sqrt{\dfrac{EI}{\rho A}} $$
$n$はモード次数（1,2,3...）、$L$は梁の長さです。$\omega_n$は$n$の2乗で増加し、長さ$L$の2乗に反比例することがわかります。シミュレーターはこの式を用いて周波数を瞬時に計算しています。
実世界での応用
建築・土木構造物の耐震設計：高層ビルや橋梁は、地震動の周波数と構造物の固有振動数が一致すると共振し、大きく揺れます。設計ではFEMを用いて固有値解析を行い、主要なモード形と周波数を把握し、それらが地震の卓越周波数と重ならないように調整します。
自動車・航空機のNVH対策：車のボディや航空機の翼の振動・騒音（Noise, Vibration, Harshness）は、エンジンや気流による強制振動が固有振動モードで増幅されることで発生します。CAEシミュレーションでモード形を事前に把握し、リブを追加するなどして剛性を上げ、問題周波数を設計範囲外にシフトさせます。
精密機械・半導体製造装置：微細な加工や計測を行う装置は、外部振動や内部モーターの振動によるわずかな変位が精度を大きく損ないます。装置フレームの固有振動モードを解析し、防振材を設置する位置や構造の弱点を特定して改良します。
楽器の音響設計：ギターのボディや太鼓の膜、ビルの壁など、あらゆるものは固有振動数で音を発します。楽器製作者は、木材の選定や形状を調整することで、望ましいモード形（音色）を実現し、望ましくないモード（雑音）を抑制しています。
よくある誤解と注意点
このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず、「n=1が一番低い周波数だから一番危険」と思いがちだけど、必ずしもそうとは限らない。確かに基本モードは最も起こりやすいけど、例えば回転機械の場合は運転速度（回転数）に応じて高次モードが励振されることもある。実務では、対象がどんな外力を受けるかを考えて、どのモードに注目すべきかを判断するんだ。
次に、「節（動かない点）は完全に固定されている」という誤解。シミュレーター上の節は、あくまで「そのモードだけ」で振動した時の話。実際の構造物は全てのモードが混ざって振動するので、節の位置は時間とともに変わったり、完全に止まったりしない。あくまで理論上の目安として捉えよう。
パラメータ設定で気をつけるのは、寸法の単位系の統一だ。例えば、長さLを[mm]、ヤング率Eを[GPa]、密度ρを[kg/m³]で入力すると、計算される固有周波数はとんでもない値になる。シミュレーター内部では単位を調整しているかもしれないが、自分で計算する時は必ずSI単位系（m, Pa, kg/m³）などに揃える癖をつけよう。E=200 GPaは200×10^9 Paだね。
関連する工学分野
固有振動の考え方は、構造力学以外の幅広い分野で顔を出す、いわば「工学の共通言語」だよ。例えば、音響工学では、スピーカーのコーン紙や楽器の胴体の振動モードが、音色や歪みに直結する。このシミュレーターで板の複雑なモード形を見たよね。あれがスピーカーで起こると、特定の周波数で音が濁る原因になるんだ。
もう一つはマイクロエレクトロメカニカルシステム（MEMS）の世界。微小なシリコンの梁や膜を振動させてセンサーやフィルターとして使うんだけど、その設計にはこの梁・板の振動理論がそのまま使える。寸法が微小だから固有周波数は非常に高く（MHz帯）、Eや密度ρの製造ばらつきが性能に直撃する。シミュレーターでパラメータを少し変えた時の周波数変化の激しさは、MEMS設計の繊細さを体感できる例だと言えるね。
さらに流体構造連成（FSI）解析にも発展する。飛行機の翼のフラッター（自励振動）や、化学プラントで配管を流れる流体による振動は、構造の固有振動と流体の力が絡み合って起こる。まずは流体を無視した「構造だけ」の固有振動を把握することが、複雑な連成現象への第一歩になるんだ。
発展的な学習のために
このシミュレーターに慣れたら、次は「連続体の理論」と「離散化されたFEM」の橋渡しを意識して学ぶのがおすすめだ。まず、板の振動の基礎式である「薄板の曲げ理論（Kirchhoff-Loveの仮定）」に触れてみよう。支配方程式は梁のをさらに発展させたもので、次の4階の偏微分方程式になるよ。
$$ D \nabla^4 w + \rho h \frac{\partial^2 w}{\partial t^2} = 0 $$
ここで、$D=Eh^3/(12(1-\nu^2))$は曲げ剛性、$h$は板厚、$\nu$はポアソン比だ。この式を様々な境界条件で解くことで、シミュレーターで見た複雑な矩形板のモード形の数式表現が得られるんだ。
そして、この連続体の微分方程式をどうやってコンピュータが解くのかを知るために、FEMの固有値解析の基礎を学ぼう。キーワードは「質量マトリクス」と「剛性マトリクス」。構造を小さな要素（例えば三角形や四角形）に分割し、それぞれの変形と質量を行列で表現する。最終的には $(\mathbf{K} - \omega^2 \mathbf{M})\mathbf{u}=0$ という一般化固有値問題に帰着する。シミュレーターの背後でも、おそらく同様の計算が高速で行われているはずだ。この先は、より現実的な3次元構造のモード解析や、減衰を考慮した複素固有値解析へとステップアップしていけるよ。