強制振動の振幅倍率
$$X = \frac{F_0/k}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$$r = \omega/\omega_n$(周波数比)
Q値:$Q = 1/(2\zeta)$
固有円振動数:$\omega_n = \sqrt{k/m}$
質量・バネ定数・減衰比を動かして共振曲線の形が変わる様子をリアルタイムで観察。Q値・半値幅・位相遅れも同時表示。CAE振動解析の核心がここにある。
$r = \omega/\omega_n$(周波数比)
Q値:$Q = 1/(2\zeta)$
固有円振動数:$\omega_n = \sqrt{k/m}$
このシミュレーターは、バネ、質量、ダッシュポット(減衰要素)からなる「1自由度減衰強制振動系」をモデル化しています。ここでの支配方程式は運動方程式で、外力 $F_0 \sin(\omega t)$ が加わったときの質量の変位 $x(t)$ を記述します。
$$m\ddot{x}+ c\dot{x}+ kx = F_0 \sin(\omega_d t)$$$m$: 質量 [kg], $c$: 減衰係数 [Ns/m], $k$: バネ定数 [N/m], $F_0$: 外力の振幅 [N], $\omega_d$: 外力の円振動数 [rad/s]
定常状態における振動振幅 $X$ は、以下の「振幅倍率」の式で与えられます。これがシミュレーターのメイングラフ(共振曲線)を描くための式です。
$$X = \frac{F_0/k}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}\quad \text{ここで}\quad r = \frac{\omega_d}{\omega_n}, \quad \omega_n = \sqrt{\frac{k}{m}}, \quad \zeta = \frac{c}{2\sqrt{mk}}$$$X$: 定常振幅 [m], $r$: 周波数比, $\omega_n$: 固有円振動数 [rad/s], $\zeta$: 減衰比 (無次元)。分母が最小、つまり $r \approx 1$(加振周波数が固有振動数に近い)で振幅 $X$ が最大(共振)します。減衰比 $\zeta$ が分母の $(2\zeta r)^2$ 項を通じて振幅を抑制します。
自動車・機械設計:エンジンのマウントや車体のサスペンション設計で、エンジン回転数や路面入力による共振を避けるために使用されます。シミュレーターで質量やバネ定数を変えながら共振点がどう移動するかを確認することは、設計パラメータの最適化に直結します。
建築・土木構造物:高層ビルや橋梁は風や地震による周期的な力を受けます。固有振動数がこれらの外力の周波数と一致しないよう設計する必要があり、その基礎理解にこの振動モデルが役立ちます。建物の減衰比は比較的小さく(0.02〜0.05)、共振が起こりやすいため注意が必要です。
家電・電子機器:洗濯機の脱水槽や冷蔵機のコンプレッサーなど、回転部を持つ機器では、起動・停止時に必ず共振周波数域を通過します。減衰を適切に設計し、共振時の振動振幅を許容範囲内に抑えることが信頼性向上の鍵となります。
音響・電子回路:スピーカーの振動板や、ラジオの同調回路(LC回路)も共振現象を利用しています。こちらは機械系よりもQ値(鋭さ)が非常に高く、特定の周波数だけを選択的に増幅・吸収することが可能です。
このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず、「共振周波数はバネ定数だけで決まる」と思いがちだけど、実際は質量との組み合わせで決まる固有振動数 $\omega_n = \sqrt{k/m}$ が基準だ。例えば、バネを2倍硬くしても、質量も同時に2倍にしたら、共振周波数は変わらないんだ。
次に、「減衰比ζを0にすれば無限に振動する」という理解。理論上はそうなるけど、現実世界では摩擦や空気抵抗が必ず存在するから、ζ=0の系は存在しない。実務で金属構造物の解析をするとき、この「見えない減衰」をどう見積もるかが結構難しい。カタログにない材料を使う場合は、似た材料の値から類推したり、実験で測定したりする必要がある。
あと、グラフの縦軸「振幅倍率」の意味をしっかり押さえよう。これは静的な変位($F_0/k$)を何倍したかを示している。例えば、振幅倍率が10と出ても、外力$F_0$が小さければ実際の変位は微少だ。共振の「危険性」を評価するには、この倍率と実際にかかる外力の大きさの両方を見る必要がある。シミュレーション結果の数値だけを見て一喜一憂しないようにね。
この1自由度振動系の考え方は、CAEの世界だけじゃなく、様々な工学分野の基礎になっている。まず音響工学だ。スピーカーや楽器の設計では、この共振曲線がそのまま周波数特性(どの音が大きく出るか)を決める。減衰比を調整することで、音の「きれいさ」や「鳴りの長さ」をコントロールしているんだ。
制御工学とも深く関係する。ロボットアームや精密ステージの位置制御では、機械的な共振周波数が制御系の帯域幅を制限する。制御の教科書で出てくる「二次遅れ系」のボード線図は、このシミュレーターの位相遅れのグラフと本質的に同じものだ。共振ピークをいかに抑えるか(=減衰を入れるか)が、安定で速い制御を実現する鍵になる。
さらに地震工学や建築構造でも必須の知識だ。建物自体を質量、柱をバネ、ダンパーを減衰とみなしてモデル化する。地震の揺れ(入力周波数)が建物の固有周期(共振周波数)に近づかないように設計したり、逆に免震構造では意図的に共振点をずらすことでエネルギーを吸収する。このツールで遊んだ感覚が、大きな構造物の安全性を考える土台になるよ。
このシミュレーターに慣れたら、次は「多自由度振動系」を学ぶのがおすすめだ。現実の構造物は複数の共振周波数(モード)を持つ。例えば自動車のボディは、グローバルな曲げモード、ねじれモード、ドアの局部振動など、無数の共振点がある。CAEソフトではモード解析として、これら全ての共振周波数と振動の形(モード形状)を一度に計算する。
数学的には、行列と固有値問題の理解が不可欠になる。1自由度の運動方程式 $m\ddot{x}+kx=0$ が、多自由度では $[M]\{\ddot{x}\} + [K]\{x\} = \{0\}$ という行列形式になる。ここで出てくる固有値が各モードの固有振動数$\omega_n^2$に、固有ベクトルがモード形状に対応する。このツールで学んだ「$\omega_n = \sqrt{k/m}$」が、どう行列の世界に拡張されるかを追ってみると、理解がぐっと深まるはずだ。
実務に直結する次のステップとしては、CAEソフト(例えばANSYSやAbaqus)で簡単な梁や板のモデルを作り、モード解析を実際に実行してみること。計算結果として出てくる複数のモードを、このシミュレーターで学んだ「1つ1つの共振曲線の重ね合わせ」としてイメージできるようになれば、振動解析の第一歩は成功だ。