回転体の動特性解析
理論と物理
回転体の振動の特殊性
先生、回転体の振動は静止構造と何が違いますか？
回転体には3つの特殊効果がある：
1. 遠心力による剛性変化（スピンスティフニング） — プレストレスモーダルで扱った
2. コリオリ力 — 回転座標系で見たとき、振動するものに作用する見かけの力
3. ジャイロ効果 — 回転する物体の向きが変わるときに生じる効果
コリオリ力とジャイロ効果が追加されるんですね。
運動方程式：
- $[G]$ — ジャイロ行列（反対称。速度に比例する力）
- $[K_\sigma]$ — 幾何剛性（遠心力によるスピンスティフニング）
- $[K_c]$ — 遠心力軟化（spin softening）
$[G]$ が反対称！ これは減衰マトリクス $[C]$ と同じ位置に入りますが、非対称ですよね。
$[G]$ は反対称（$[G]^T = -[G]$）で、系にエネルギーを加えることも抜くこともしない。しかしモードの分裂（前進波と後退波の振動数が異なる現象）を引き起こす。
キャンベルダイアグラム
回転速度 $\Omega$ vs. 固有振動数 $f$ をプロットしたキャンベルダイアグラムが回転体振動の基本ツール。
特徴：
- 前進波（Forward Whirl） — 回転方向と同じ向きの振動。$\Omega$ が増すと振動数が上昇
- 後退波（Backward Whirl） — 回転方向と逆向き。$\Omega$ が増すと振動数が低下
- 励振線（$f = n\Omega$） — 不釣り合い（$n=1$）や高次励振との交点が共振
前進波と後退波に分裂する…これがジャイロ効果ですか。
そう。静止状態では同じ振動数だった2つのモードが、回転速度の増加とともに分裂する。タービンブレードやロータの設計ではこの分裂を正確に予測する必要がある。
危険速度
危険速度（critical speed）は励振線と固有振動数線の交点に対応する回転速度。特に不釣り合い励振（$n=1$）との交点が最も重要。
危険速度を避けて運転するのが基本ですか？
基本的にはそう。API 617（コンプレッサー）やAPI 612（蒸気タービン）では、危険速度から±15%のマージンを設定する。
まとめ
回転体の振動を整理します。
要点：
- 3つの特殊効果 — 遠心力剛性変化、コリオリ力、ジャイロ効果
- ジャイロ行列 $[G]$ は反対称 — モードの前進波/後退波分裂を引き起こす
- キャンベルダイアグラム — 回転速度 vs. 固有振動数。共振条件の特定
- 危険速度 — 励振線との交点。±15%のマージンで回避
- 複素固有値解析が必要 — ジャイロ効果で固有値が複素になる
ジャイロ効果の発見と応用の歴史
ジャイロ効果は1852年にフランスの物理学者レオン・フーコーが独楽の安定性を研究中に発見した。ジャイロスコープという名前も彼がつけた。回転体の角運動量がモーメントに抵抗する性質は1900年代初頭にラジオコンパスに応用され、現在もINS（慣性航法装置）の核心技術だ。FEMでのジャイロ効果は反対称陀螺行列[G]として表現される。
各項の物理的意味
- 慣性項（質量項）：$\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での回転体解析
回転体の振動解析をFEMでどうやりますか？
Nastran
```
SOL 107 $ 複素固有値（ジャイロ効果含む）
CEND
CMETHOD = 10
BEGIN BULK
RFORCE, 100, 1, , 100., 0., 0., 1. $ 回転速度100 rad/s, z軸回転
```
SOL 107でジャイロ行列を含む複素固有値解析。RFORCEで回転速度と軸を指定。
Abaqus
```
*STEP
*COMPLEX FREQUENCY, CORIOLIS=ON
20, ,
*DLOAD
element_set, CENTRIF, omega_squared, x, y, z
*END STEP
```
CORIOLIS=ONでジャイロ効果を有効化。
Ansys
```
/SOLU
ANTYPE, MODAL
MODOPT, QRDAMP, 20
CORIOLIS, ON, , ON ! コリオリ/ジャイロ効果を有効化
OMEGA, , , 100. ! z軸回転100 rad/s
SOLVE
```
AnsysではCORIOLISコマンドでジャイロ効果を有効化するんですね。
AnsysのQRDAMP法は実モードに投射してから減衰/ジャイロ効果を含む複素固有値を求める。大規模モデルに効率的だ。
キャンベルダイアグラムの作成
手順：
1. 回転速度を10〜20点設定（0から最大運転速度まで）
2. 各速度でプレストレスモーダル＋複素固有値解析
3. 回転速度 vs. 固有振動数をプロット
4. 励振線（$f = \Omega, 2\Omega, 3\Omega, ...$）を重ね書き
5. 交点を危険速度として特定
Ansys Workbenchでは自動化されていますよね。
Ansys Workbenchの「Rotordynamics Analysis」でキャンベルダイアグラムを自動生成できる。回転速度のパラメトリックスイープと結果のプロットが一体化している。
まとめ
回転体の数値手法、整理します。
要点：
- 複素固有値解析が必要 — ジャイロ行列 $[G]$ で複素固有値
- SOL 107（Nastran）, COMPLEX FREQUENCY CORIOLIS=ON（Abaqus）, QRDAMP+CORIOLIS（Ansys）
- キャンベルダイアグラムを自動生成 — Ansysが最も手軽
- 回転速度を変えた多数ケースの解析 — パラメトリックスタディ
Campbell図の作成と危険速度の特定
Campbell図は縦軸を固有振動数、横軸を回転数として各回転数での固有振動数曲線と回転数の整数倍（エンジンオーダー）直線を描いた図だ。交点が「危険速度（クリティカルスピード）」となり、その回転数での運転を避けるか共鳴通過を素早く行う設計とする。ISO 10816ではキャンベル図を用いた危険速度の評価が回転機械の設計標準として要求されている。
線形要素（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要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
回転体振動の実務
回転体の振動解析は実務でどう使われていますか？
回転機械の設計では危険速度の回避が最優先事項。API規格（617, 612, 610等）で規定。
実務チェックリスト
軸受の剛性・減衰が重要なんですね。
回転体の振動特性は軸受の特性に大きく依存する。軸受の剛性と減衰は回転速度に依存するため、キャンベルダイアグラムの各点で正しい軸受特性を使う必要がある。
航空機エンジンファンブレードのスピン試験
GE90エンジン（B777搭載）のファンブレードはチタン製で直径3.3m・総重量280kgに達する。スピン試験では25,000rpmまで段階的に回転数を上げ、各回転数でのFRF計測によりCampbell図を確認する。FEM予測との整合性が±5%以内に収まることが型式証明の要件で、ブレード1枚ごとに個別の質量・重心測定によるバランス調整も必須だ。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
回転体解析のツール
回転体の振動解析にはどんなツールが使えますか？
| ツール | 特徴 |
|---|---|
| Ansys Rotordynamics | WorkbenchでCampbell Diagram自動生成。GUI最良 |
| Nastran SOL 107/110 | 航空宇宙のタービン/コンプレッサー |
| MADYN 2000 | ロータダイナミクス専用。軸受データベース内蔵 |
| DyRoBeS | ロータダイナミクス専用。API準拠 |
| XLTRC2 | テキサスA&M大学。ロータ解析の研究標準 |
選定ガイド
専用ツール（MADYN, DyRoBeS）が軸受データベースを持っているのが強みですね。
回転体の振動は軸受特性が支配する。軸受メーカーのデータベースと連携できる専用ツールが実務では不可欠。
AxSTREAMとRotordynamics専用ツール
SoftInWay社のAxSTREAMはタービン・コンプレッサー設計専用のRotordynamics解析モジュールを持ち、API 684要求のLateral Analysis（横方向回転動力学）とTorsional Analysis（ねじり振動）を1ツールで実行できる。GE Power子会社では発電用蒸気タービンのロータ設計に採用し、クリティカルスピード予測精度を2%以内に達成している。
選定で最も重要な3つの問い
- 「何を解くか」：回転体の動特性解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
回転体振動の先端研究
回転体振動の最前線を教えてください。
非線形ロータダイナミクス
軸受の非線形性（油膜のキャビテーション、接触）や大変形を含む非線形ロータダイナミクス。時刻歴解析で非線形応答を直接シミュレーション。ルーンバーグ解析やHBM（Harmonic Balance Method）で効率化。
磁気浮上軸受（AMB）との連成
能動磁気浮上軸受（AMB: Active Magnetic Bearing）のフィードバック制御と構造の振動を連成解析。制御系と構造のモデルをFEMと制御工学の双方で扱う。
ブリスク（ブレード一体ディスク）の振動
ブリスク（blisk: bladed disk）はブレードとディスクが一体成形された構造。ブレード間のミスチューン（製造誤差による微小な差異）で振動特性が劇的に変化する。確率論的FEMでミスチューンの影響を評価。
まとめ
回転体振動の先端研究、まとめます。
ロータダイナミクスの動的安定性解析
回転機械の不安定振動（ホワール）はジャイロ効果だけでなく流体軸受の油膜力・内部減衰・クロスカップリング剛性によっても生じる。軸剛性と速度域によって「前進ホワール（同方向回転）」と「後退ホワール（逆方向回転）」が現れ、油膜軸受では前進ホワールが不安定になりやすい。API 684はロータダイナミクスの解析要件を規定しており、圧縮機・ポンプの設計でのコンプライアンス確認に使われる。
トラブルシューティング
回転体振動のトラブル
回転体の振動解析でよくあるトラブルは？
前進波と後退波の分裂が出ない
ジャイロ効果が有効になっていない。確認：
- Nastran: SOL 107を使っているか（SOL 103ではジャイロ非対応）
- Abaqus: CORIOLIS=ONが設定されているか
- Ansys: CORIOLIS, ON が設定されているか
危険速度が実験と合わない
確認項目：
- 軸受の剛性が正しいか — 軸受剛性は回転速度に依存。正しい速度の値を使っているか
- 回転体の質量分布 — ディスクの質量、ブレードの質量が正しいか
- 軸の剛性 — 段付き軸の断面諸元
キャンベルダイアグラムの線が交差しない
危険速度が運転範囲外にある可能性。運転速度の範囲を超えた範囲までプロットして確認。
まとめ
回転体振動のトラブル対処、整理します。
回転試験でのCampbell図と予測の不一致
Campbell図で実測固有振動数がFEM予測より低い場合、ディスクハブの接触剛性が過大評価されている可能性がある。ブレードとディスクの接触部（ダブテール）の摩擦スティフネスは接触圧・摩擦係数・面粗さで変わり、FEM設定のデフォルト接触剛性は往々にして過大だ。回転数0での実測FRFとFEMを比較し、接触剛性を実測値にキャリブレートしてからCampbell予測に進むこと。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——回転体の動特性解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告