梁の自由振動解析
理論と物理
梁の振動 — 動的解析の原点
先生、梁の自由振動は固有振動数解析の最も基本的な問題ですか？
そう。梁の自由振動は理論解が存在する数少ない問題の一つであり、FEMの精度検証に不可欠だ。
支配方程式
オイラー・ベルヌーイ梁の自由振動方程式：
分離変数法 $w(x,t) = W(x) e^{i\omega t}$ で空間部分：
4階の常微分方程式ですね。解はどうなりますか？
一般解：
ここで $\beta^4 = \rho A \omega^2 / (EI)$。
4つの定数 $C_1 \sim C_4$ は4つの境界条件で決まる。非自明解の条件（振動数方程式）から固有振動数が得られる。
各境界条件の振動数
$n$ 次の固有振動数 $f_n = (\beta_n L)^2 / (2\pi L^2) \sqrt{EI/(\rho A)}$：
| 境界条件 | $(\beta_1 L)^2$ | $(\beta_2 L)^2$ | $(\beta_3 L)^2$ | モード形状 |
|---|---|---|---|---|
| 片持ち | 3.516 | 22.03 | 61.70 | 先端が最大 |
| 単純支持 | $\pi^2 = 9.870$ | $4\pi^2 = 39.48$ | $9\pi^2 = 88.83$ | $\sin(n\pi x/L)$ |
| 両端固定 | 22.37 | 61.67 | 120.9 | 両端ゼロ |
| 自由-自由 | 22.37 | 61.67 | 120.9 | 両端が自由端 |
片持ちと両端固定で1次の $(\beta L)^2$ が3.516 vs. 22.37…6倍以上の差。振動数は $\sqrt{6} \approx 2.5$ 倍ですか？
$f \propto (\beta L)^2$ だから22.37/3.516 = 6.36倍。境界条件で振動数が6倍以上変わる。FEMで理論値と合わない場合、まず境界条件を疑う。
ティモシェンコ梁の振動
ティモシェンコ梁の振動は？
せん断変形と回転慣性を含むため、高次モードでEB梁と差が出る：
常にティモシェンコ梁の方が低い振動数。差は$f \cdot h / c_s$（$c_s$ = せん断波速度）が大きいとき顕著。
高次モードほど差が大きい？
そう。1次モードでは差がわずか（1〜2%）だが、10次モードでは10〜30%の差になることがある。高次モードまで正確に求めたい場合はティモシェンコ梁が必要。
FEMの検証
梁の自由振動はFEMのベンチマーク問題として最適：
1. 理論解が正確にわかっている
2. 要素数を変えてメッシュ収束を確認できる
3. EB梁要素とティモシェンコ梁要素の差を体験できる
4. 梁要素とシェル/ソリッド要素の比較ができる
FEMを学ぶ人にとって最良の練習問題ですね。
片持ち梁の1次固有振動数をFEMで計算し、理論値 $f_1 = 3.516/(2\pi L^2) \sqrt{EI/\rho A}$ と比較する。これだけで要素の精度、境界条件、質量の設定が正しいか全て確認できる。
まとめ
梁の自由振動を整理します。
要点：
- $EI W'''' = \rho A \omega^2 W$ — 4階ODEの固有値問題
- 理論解が存在 — FEMの最良のベンチマーク
- 境界条件で振動数が6倍以上変わる — 最初にチェック
- 高次モードではティモシェンコ梁が必要 — EB梁は高次で不正確
- $(\beta_n L)^2$ の値を暗記 — 片持ち: 3.516, 単純梁: $\pi^2$
Euler-Bernoulliはりの歴史
「オイラー-ベルヌーイはり理論」はレオンハルト・オイラー（1744年）とヤコブ・ベルヌーイによる協作だ。当初ベルヌーイが横たわりから見た曲げを研究し、オイラーがそれを一般化した。Timoshenkoが1921年にせん断変形と回転慣性を加えた改良版を発表するまで、これが唯一の梁振動理論だった。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
FEMによる梁の振動解析
梁要素で振動解析するとき、要素数はどの程度必要ですか？
EB梁要素（エルミート補間）は曲げの4次多項式を正確に表現するから、静解析では1要素で十分だった。しかし振動解析ではモード形状を表現するために複数要素が必要。
$n$ 次モードの半波長を表現するには：
- 1次モード: 最低4要素
- 5次モード: 最低20要素
- $n$ 次モード: 最低 $4n$ 要素（目安）
静解析より多くの要素が必要なんですね。
振動のモード形状は正弦/余弦の組み合わせで、高次モードほど波長が短い。短い波長を解像するには細かいメッシュが必要。これは板やソリッドでも同じだ。
質量マトリクスの影響
一貫質量と集中質量で結果が変わりますか？
EB梁要素の場合：
- 一貫質量 — $f$ を過大評価する傾向（剛性の上界定理）
- 集中質量 — $f$ を過小評価する傾向
- 真値はこの間
どちらが正確ですか？
一般に一貫質量のほうが正確。ただし要素数が十分多ければ差は小さい。EB梁要素4個で片持ち梁の1次固有振動数を求めると、一貫質量で0.1%以内の精度が出る。
ソルバー別の設定
```
SOL 103
CEND
METHOD = 10
BEGIN BULK
CBAR, ...
EIGRL, 10, , , 10
```
```
*BEAM SECTION, SECTION=RECT, ELSET=beam
0.01, 0.001
*STEP
*FREQUENCY, EIGENSOLVER=LANCZOS
10, ,
*END STEP
```
梁要素の固有振動数解析はシンプルですね。
梁の振動は「設定がシンプルで理論解と直接比較できる」理想的なFEM演習だ。新しいソルバーを使い始めるとき、まず梁の振動解析を行って設定の確認をすることを推奨する。
まとめ
梁の振動の数値手法、整理します。
要点：
- $n$ 次モードに最低 $4n$ 要素 — 静解析より多くの要素が必要
- 一貫質量のほうが正確 — $f$ の過大評価傾向だが精度は高い
- EB梁要素4個で1次モード0.1%精度 — 非常に効率的
- FEMの最初の検証に最適 — 新ソルバーの設定確認に使う
固有振動数の解析解と境界条件
両端固定はりの固有振動数はfn=（βnL)²/(2πL²)√(EI/ρA）で与えられ、βnLは超越方程式cos(βnL)cosh(βnL)=1の解（第1次=4.730）だ。両端単純支持ならβnL=nπで明示的に求まる。FEM結果との比較では、1次モードで要素数10個以上あれば誤差0.1%以内に収まる。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
梁振動の実務適用
梁の自由振動解析は実務でどう使いますか？
配管振動の例
ポンプに接続された配管のスパン $L = 3$ m、外径 $D = 100$ mm、肉厚 $t = 5$ mm、鋼管。
ポンプの回転数が1750 rpmなら振動数は29.2 Hz…共振！
まさにポンプの振動数と一致してしまう！
対策：スパンを短くする（支持点を追加）or 配管に質量を追加して振動数を下げる。手計算でこの判断ができるのが梁振動理論の強みだ。
実務チェックリスト
梁振動解析のチェックリストは？
「流体の付加質量」は見落としがちですね。
水が入った配管の実効密度は鋼管だけの2倍近くになることがある。配管振動では管内流体の質量を忘れると固有振動数が $\sqrt{2} \approx 1.4$ 倍ずれる。
橋梁の固有振動数と歩行者共振
歩行者が橋を歩く際の歩行周波数は1.6〜2.4Hz（歩行）、2.4〜4.0Hz（走行）だ。橋の1次固有振動数がこの範囲に入ると共振が起き、2000年にロンドンのミレニアムブリッジが開通2日目に横揺れで歩行者が立てなくなった事件が有名だ。質量ダンパー（TMD）を設置して0.8Hzから外し、問題を解決した。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
梁振動のツール
梁振動の解析にはどんなツールが使えますか？
手計算からFEMまで幅広い選択肢がある。
| ツール | 特徴 |
|---|---|
| 手計算（$f_n$ の公式） | 最速。概算に最適 |
| Excel/Python | 公式のパラメトリック計算 |
| CAESAR II | 配管スパンの振動評価 |
| 汎用FEM | 複雑な境界条件、集中質量、変断面に対応 |
手計算で十分な場合も多いんですね。
一様断面の梁なら手計算で十分。変断面、集中質量、非一様な境界条件がある場合はFEMが必要。
選定ガイド
梁振動は「手計算が最初のツール」。FEMは複雑な場合のみ。
手計算の理論解がFEMの結果を検証する「物差し」になる。理論解を知らずにFEMだけで振動解析をするのは、地図なしで旅行するようなものだ。
COMSOLのはり固有値解析
COMSOLのBeam物理インターフェースはEuler-BernoulliとTimoshenkoを選択でき、断面特性（I・A・J等）を自動計算してから固有値を求める。MEMSセンサー設計では共振周波数を0.1%以内の精度で予測する必要があり、Timoshenko補正とアスペクト比L/h>20の条件を満たすメッシュ設定が精度確保の鍵だ。
選定で最も重要な3つの問い
- 「何を解くか」：梁の自由振動解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
梁振動の先端トピック
梁振動に先端研究はありますか？
梁自体は古典的だが、新しい応用がある。
MEMS共振器
MEMS（微小電気機械システム）のカンチレバーは梁の自由振動そのもの。質量センサー（バイオセンサー）、周波数標準（共振器）に使われる。固有振動数の変化でピコグラムレベルの質量変化を検出する。
ナノスケールの振動
カーボンナノチューブやグラフェンの振動は非局所弾性理論（Eringen理論）に基づくEB/ティモシェンコ梁として記述される。分子間力の長距離効果で古典的梁理論とは異なる振動数が得られる。
エネルギーハーベスティング
圧電カンチレバー梁の振動で環境振動から発電する。梁の固有振動数を環境振動の周波数に合わせて共振させることでエネルギー変換効率を最大化する。
まとめ
梁振動の先端研究、まとめます。
18世紀の梁理論が21世紀のナノテクノロジーとエネルギー技術を支えている。
非線形大変形はりの振動周波数シフト
大変形になるとはりの曲げ剛性が変形に依存し、振動周波数が振幅に依存する「非線形振動」になる。振幅が長さの10%以上になると固有振動数が最大20%上昇するhardening型非線形が生じる。MEMS（微小電気機械）のナノはりではこの効果が支配的で、共振周波数が変位センサとして機能する。
トラブルシューティング
梁振動のトラブル
梁の振動解析でよくあるトラブルは？
1次固有振動数が理論値と合わない
確認項目（優先度順）：
1. 密度 — $f \propto 1/\sqrt{\rho}$。単位を確認
2. 境界条件 — ピン(回転自由) vs. 固定(回転拘束)で振動数が大幅に変わる
3. 断面諸元（$I$） — 強軸/弱軸の間違い
4. 梁理論の種類 — EB梁 vs. ティモシェンコ梁で差が出る（特に太短い梁）
5. 要素数 — 高次モードでは要素数不足が精度低下
ティモシェンコ梁の結果がEB理論と一致しない
正常な結果。ティモシェンコ梁はせん断変形を含むから、EB梁より低い振動数が出る。差を確認するには $L/h$ を大きくしてEB梁の理論値に収束するか見る。
高次モードが不正確
高次モードの波長が短くなると、メッシュが不十分になりやすい。$n$ 次モードに最低 $4n$ 要素。要素数を増やして収束を確認する。
まとめ
梁振動のトラブル対処、整理します。
FEMはり解析で高次モードが不正確な場合
FEMで高次固有モードの精度が落ちる主因はメッシュ不足と数値積分次数の不適切な選択だ。欲しい最高モード次数をNとすると、要素数は最低N×4以上必要という経験則がある。また縮減積分要素（Abaqus B21R等）は低次モードには良いが高次では誤差が大きくなるため、完全積分要素B22の使用が推奨される。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——梁の自由振動解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告