回転体の動特性解析
回転体の動特性の理論基礎
回転体の振動の特殊性
先生、回転体の振動は静止構造と何が違いますか？
回転体には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]として表現される。
回転体の動特性の数値計算手法
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ではキャンベル図を用いた危険速度の評価が回転機械の設計標準として要求されている。
回転体の動特性の実務適用
回転体振動の実務
回転体の振動解析は実務でどう使われていますか？
回転機械の設計では危険速度の回避が最優先事項。API規格（617, 612, 610等）で規定。
実務チェックリスト
軸受の剛性・減衰が重要なんですね。
回転体の振動特性は軸受の特性に大きく依存する。軸受の剛性と減衰は回転速度に依存するため、キャンベルダイアグラムの各点で正しい軸受特性を使う必要がある。
航空機エンジンファンブレードのスピン試験
GE90エンジン（B777搭載）のファンブレードはチタン製で直径3.3m・総重量280kgに達する。スピン試験では25,000rpmまで段階的に回転数を上げ、各回転数でのFRF計測によりCampbell図を確認する。FEM予測との整合性が±5%以内に収まることが型式証明の要件で、ブレード1枚ごとに個別の質量・重心測定によるバランス調整も必須だ。
回転体の動特性のソフトウェア比較
回転体解析のツール
回転体の振動解析にはどんなツールが使えますか？
| ツール | 特徴 |
|---|---|
| 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%以内に達成している。
回転体の動特性の先端研究
回転体振動の先端研究
回転体振動の最前線を教えてください。
非線形ロータダイナミクス
軸受の非線形性（油膜のキャビテーション、接触）や大変形を含む非線形ロータダイナミクス。時刻歴解析で非線形応答を直接シミュレーション。ルーンバーグ解析や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予測に進むこと。
関連トピック
なった
詳しく
報告