梁の自由振動解析
梁の自由振動の理論基礎
梁の振動 — 動的解析の原点
先生、梁の自由振動は固有振動数解析の最も基本的な問題ですか？
そう。梁の自由振動は理論解が存在する数少ない問題の一つであり、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年にせん断変形と回転慣性を加えた改良版を発表するまで、これが唯一の梁振動理論だった。
梁の自由振動の数値計算手法
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%以内に収まる。
梁の自由振動の実務適用
梁振動の実務適用
梁の自由振動解析は実務でどう使いますか？
配管振動の例
ポンプに接続された配管のスパン $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から外し、問題を解決した。
梁の自由振動のソフトウェア比較
梁振動のツール
梁振動の解析にはどんなツールが使えますか？
手計算から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の条件を満たすメッシュ設定が精度確保の鍵だ。
梁の自由振動の先端研究
梁振動の先端トピック
梁振動に先端研究はありますか？
梁自体は古典的だが、新しい応用がある。
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の使用が推奨される。
関連トピック
なった
詳しく
報告