$\varphi_n(x) = \sin\!\left(\dfrac{n\pi x}{L}\right)$
$\omega_n = \left(\dfrac{n\pi}{L}\right)^{\!2}\!\sqrt{\dfrac{EI}{\rho A}}$
単純支持・片持ち・両端固定梁の固有振動モード形をリアルタイムアニメーション。固有周波数を自動計算。矩形板の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だね。
長さL=1.5mの鋼製片持ち梁(E=200GPa、ρ=7850kg/m³、断面二次モーメントI=8.33×10⁻⁶m⁴)の場合、1次固有周波数f₁≈18.2Hzで先端が最大変位を示します。幅0.5m、厚さ0.01mのアルミ板(E=73GPa、ρ=2780kg/m³)を両端固定で解析すると、2次モード周波数は約156Hzとなり、板面内に複数の腹が形成される様子をアニメーションで確認できます。