単純支持・片持ち・固定梁の固有モードをアニメーション表示。複数モードの重ね合わせと固有振動数スペクトルをリアルタイムで確認できます。
梁の横振動を記述する基本方程式は、オイラー-ベルヌーイの梁理論に基づく次の4階微分方程式です。
$$EI \frac{\partial^4 w(x, t)}{\partial x^4}+ \rho A \frac{\partial^2 w(x, t)}{\partial t^2}= 0$$ここで、$w(x,t)$は位置$x$、時間$t$における梁の横たわみ、$EI$は曲げ剛性(材質と断面形状で決まる)、$\rho A$は単位長さあたりの質量です。この式は「曲げ剛性による復元力」と「質量による慣性力」のバランスを表しています。
上記の方程式から導かれる固有振動数 $f_n$ は、境界条件に依存するパラメータ $\beta_n L$ を用いて次のように計算されます。
$$f_n = \frac{(\beta_n L)^2}{2\pi L^2}\sqrt{\frac{EI}{\rho A}}$$$L$は梁の長さ、$n$はモード次数(1,2,3...)です。$\beta_n L$は特性値で、境界条件によって決まる数値です(例:単純支持梁では$n\pi$、片持ち梁1次モードでは約1.875)。この式から、梁が長いほど、または剛性($EI$)が高いほど固有振動数が高くなることがわかります。
自動車・航空機の構造設計:ボディや翼のフレームは梁構造とみなせます。走行風やエンジン振動による共振を避けるため、固有振動数を計算し、危険な周波数帯からずらす設計が行われます。シミュレーターで学ぶモード形状は、どこを補強すれば剛性が上がるかを直感的に理解する助けになります。
建築・橋梁の耐震設計:高層ビルの骨組や橋桁は、地震動による横揺れ(曲げ振動)を受けます。固有振動モードを把握することで、地震エネルギーがどのように建物に吸収されるかを予測し、制振装置を設置する最適な位置を決定します。
機械の振動トラブルシューティング:回転機械(ポンプ、ファン、タービン)から特定の周波数で異常振動が発生した場合、その周波数が支持構造(架台や配管)の何次モードの固有振動数と一致していないか調査します。一致していれば、共振が原因と特定し、剛性追加や質量変更などの対策を講じます。
電子部品の信頼性評価:基板上の大型ICやヒートシンクは片持ち梁とみなせます。輸送中の振動や衝撃で破損しないよう、その固有振動数を評価します。シミュレーターの片持ち梁モードは、はんだ接合部に最もストレスがかかる形状を可視化してくれます。
このツールで遊んでいると、いくつか勘違いしやすいポイントがあるよ。まず、「モード次数が高いほど振動が激しい」と思いがちだけど、実は違うんだ。モードが上がると固有振動数(周波数)は確かに高くなるけど、実際にそのモードでどれだけ激しく振動するかは、加わる外力の周波数とエネルギー次第。例えば、モード5の周波数が100Hzだとして、外力が10Hzならほとんど励起されない。逆に、基本モード(モード1)の共振が起きると、大きな振幅で壊れることもあるから、基本モードこそ最重要なんだ。
次に、「重ね合わせ」は単なる足し算じゃないってこと。シミュレーターでは各モードの振幅を足しているように見えるけど、実際の動的応答では位相(振動のタイミング)が超重要。例えばモード1とモード2が同じ振幅でも、位相が逆なら打ち消し合って振動は小さくなる。CAEのモード重ね合わせ法では、この位相関係をきちんと考慮して計算するんだ。
最後に、ツールのパラメータで気をつけてほしいのは、「曲げ剛性(EI)」は断面形状の影響が大きいってこと。材質のヤング率Eだけ変えても、断面二次モーメントIは変わらないよね。例えば、幅10mm、高さ20mmの長方形断面の梁のIは、高さ方向で計算すると約6667 mm^4。これを「幅20mm、高さ10mm」に寝かせると、Iは約1667 mm^4に激減する。つまり、梁の向きを変えるだけで剛性が1/4になり、固有振動数は半分になってしまうんだ。設計ではこの「断面の向き」が命取りになることもあるから要注意だね。
梁の振動解析は、まさに構造力学の基礎中の基礎。これがわかると、実に様々な工学分野の扉が開けるんだ。まず直接繋がるのは「モード解析」や「実験モード解析(EMA)」。CAEで固有値解析をした結果を、実際の構造物に加速度センサーを貼ってハンマーで叩き(インパクト試験)、計測した周波数応答関数(FRF)と比較して検証する、という流れだ。このツールで学んだ「節」の位置は、センサーを貼ってはいけない場所(振動が検出できない)でもあるんだよ。
また、自動車や航空機の「NVH(騒音・振動・ハーシュネス)」評価にも直結する。車のドアミラーの「ブルブル」という振動や、高速走行時の「ブーン」という唸り音は、まさに梁や板の特定のモードが励起されて起こる現象。エンジニアはこのツールで見ているようなモード形状を頭に描きながら、対策(リブ追加や質量付加)を考えるんだ。
さらに発展すると、「波動伝播」や「音響学」にも繋がる。音そのものが空気の縦波だが、固体中を伝わる弾性波(特に曲げ波)の理解は、機械の騒音低減や、建物の防音設計に不可欠。梁の振動方程式を少し拡張して損失項を入れれば、振動がどのように減衰していくかもモデル化できる。こうして、一つの梁のモデルから、構造、音響、制御まで幅広い分野が統合的に理解できるようになるんだ。
このシミュレーターに慣れたら、次は「連続体としての梁」から「離散化モデル」への頭の切り替えをやってみよう。CAEソフトは実際、梁をたくさんの小さな要素(有限要素)に分割して計算している。例えば、このツールの片持ち梁を10個の要素に分けて考え、各節点の質量と剛性をマトリックス(行列)で表すと、運動方程式は $M\ddot{u} + Ku = 0$ という形になる。ここで$M$は質量行列、$K$は剛性行列、$u$は変位ベクトルだ。この方程式の固有値と固有ベクトルを解くことが、そのままCAEソフトの「固有値解析」の核心なんだ。
数学的には、オイラー-ベルヌーイの微分方程式を「変数分離法」で解く過程を追ってみることをお勧めする。解が $w(x,t) = W(x) \cdot T(t)$ と分離できると仮定すると、空間項$W(x)$はサイン、コサイン、双曲線関数の組み合わせになり、そこに境界条件を代入することであの$\beta_n L$(特性値)が求まる。この一連の流れを紙に書いて追うことで、ツールの背後にある「魔法」が確かな「論理」に変わるはずだ。
次の具体的なトピックとしては、「梁の強制振動と周波数応答」を学ぶのが良いね。今は自由振動(自然に揺れる様子)を見ているが、実際にはエンジンなどから特定の周波数で力が加わり続ける(強制振動)。その時に振幅がどうなるか、特に共振点付近でどのくらい大きくなるかは、減衰(ダンピング)の大きさで決まる。減衰を考慮したモデルを学べば、振動解析の実務への第一歩はほぼ完璧だと言っていいだろう。