円運動と遠心力 — ω²r加速度からタービン・回転機械CAEまで
1. 円運動と回転機械CAE — タービン翼の遠心荷重
ジェットエンジン、蒸気タービン、ターボチャージャー、モーター、遠心ポンプ——現代の機械システムの多くは高速回転する部品を含む。回転体の設計では、回転によって生じる遠心荷重が支配的な機械荷重となり、これに熱荷重・ガス荷重・振動荷重が重なる。
円運動の物理は高校で学ぶが、それをジェットエンジンのタービンディスクや風力発電ロータの設計に接続するためには、座標系の概念・FEM荷重設定・回転振動解析（キャンベル線図）まで理解が必要だ。
2. 向心加速度と向心力
等速円運動する質点は速さが一定でも方向が常に変化するため、加速度を持つ。この加速度は常に円の中心方向（向心方向）を向き、大きさは:
$$a_c = \frac{v^2}{r} = \omega^2 r$$ここで $v$ は速さ [m/s]、$r$ は半径 [m]、$\omega$ は角速度 [rad/s]。対応する向心力は:
$$F_c = \frac{mv^2}{r} = m\omega^2 r$$向心力は「中心に向かって引っ張る力」——ひもにつないだ球の場合はひもの張力、地球を回る衛星は重力、カーブ走行車両は路面との摩擦力が向心力の役割を担う。
カーブ走行と横G
自動車が半径 $r$ のカーブを速度 $v$ で走るとき、必要な向心力（摩擦力）は $F = mv^2/r$。横加速度（横G）は $a = v^2/r$。乾燥路面の摩擦係数を $\mu = 0.8$ とすると、限界速度は $v_{max} = \sqrt{\mu g r}$。半径 50 m のコーナーでは約 72 km/h が限界だ。これを超えると横滑りするから、タイヤ–路面連成の FEM シミュレーション（Abaqus + タイヤモデル）で走行安定性を評価する。
3. 角速度・角加速度・線速度の関係
角速度 $\omega$ [rad/s] は角度の時間変化率:
$$\omega = \frac{d\theta}{dt}$$角加速度 $\alpha$ [rad/s²] はその時間変化率:
$$\alpha = \frac{d\omega}{dt}$$半径 $r$ の点での線速度・接線加速度:
$$v = r\omega, \quad a_t = r\alpha$$RPMとrad/sの換算
$$\omega \,[\text{rad/s}] = \frac{2\pi n}{60} \quad (n\text{ : 回転数 [rpm]})$$| 機器 | 典型回転数 [rpm] | ω [rad/s] |
|---|---|---|
| 洗濯機（脱水） | 1,200 | ~126 |
| 自動車エンジン（常用） | 2,000〜4,000 | 210〜420 |
| 電動モーター（産業用） | 1,500〜3,600 | 157〜377 |
| ターボチャージャー | 100,000〜200,000 | 10,500〜21,000 |
| ジェットエンジン高圧タービン | 12,000〜16,000 | 1,260〜1,680 |
| 歯科用ドリル | 300,000〜500,000 | 31,400〜52,400 |
4. 遠心力（非慣性系の擬似力）とCAEでの取り扱い
慣性系（絶対座標系）では遠心力は存在しない——向心力がなければ質点は直進するだけだ。しかし回転座標系（ターボ機械と一緒に回る座標系）で運動を記述すると、座標系の非慣性性の補正として「遠心力」と「コリオリ力」の2つの擬似力が現れる。
$$\vec{F}_{cf} = -m\vec{\omega}\times(\vec{\omega}\times\vec{r}) = m\omega^2 r\,\hat{r} \quad \text{（遠心力：外向き）}$$ $$\vec{F}_{Cor} = -2m\vec{\omega}\times\vec{v} \quad \text{（コリオリ力：速度に垂直）}$$*DLOAD, TYPE=CENTRIF に角速度と回転軸を指定するだけで自動計算してくれる。
🧑🎓 学生
コリオリ力は回転体FEMでも考慮する必要がありますか？
🎓 博士
静解析では静止状態を扱うので不要だが、回転体の固有振動数を求める「回転振動解析」ではコリオリ力を考慮しないと進行波・後退波の分岐（コリオリ効果による振動数分裂）を正確に評価できない。Abaqusの *FREQUENCY に CORIOLIS=YES を付けるのが重要だ。
5. タービン翼の遠心応力
均一断面・等密度のバー（長さ方向に均一な丸棒）が回転軸から距離 $r_h$〜$r_t$ の範囲に存在するとき、付け根（$r = r_h$）での遠心応力は:
$$\sigma_{centrifugal} = \frac{1}{2}\rho\omega^2(r_t^2 - r_h^2)$$ここで $\rho$ は材料密度 [kg/m³]。例えばニッケル超合金（$\rho = 8,200\,\text{kg/m}^3$）のタービン翼、$\omega = 1,500\,\text{rad/s}$、$r_t = 0.25\,\text{m}$、$r_h = 0.10\,\text{m}$ の場合:
$$\sigma = \frac{1}{2}\times 8200 \times 1500^2 \times (0.25^2 - 0.10^2) \approx 570\,\text{MPa}$$これに熱応力・ガス荷重が加わる。ニッケル超合金の1,000℃での降伏応力は~200〜600 MPaまで下がるため、材料選択と冷却設計が極めて重要だ。
テーパー翼での応力最適化
実際のタービン翼は先端ほど断面積が小さい（テーパー形状）。先端の質量を減らすことで付け根への遠心荷重を低減できる。Frenet積分を数値的に解いて最適テーパー比を求め、FEMで確認するのが実務の流れだ。
6. 回転体の軸対称FEM
タービンディスク、フライホイール、ロールなど軸対称形状の回転体は、軸対称解析（Axisymmetric Analysis）を使うことで3D問題を2D（断面形状のみ）に削減できる。計算コストを100分の1以下に抑えられる。
軸対称要素での遠心力
軸対称体に遠心荷重 $\rho\omega^2 r$（径方向体積力）を入力すると、軸対称2D解析で径方向応力 $\sigma_r$、周方向応力 $\sigma_\theta$（フープ応力）、軸方向応力 $\sigma_z$ の3成分が計算される。フープ応力 $\sigma_\theta$ が最大になり、ディスクの降伏を支配することが多い。
$$\sigma_\theta \text{（フープ応力）} = \rho\omega^2\frac{3+\nu}{8}\left(r_o^2 + r_i^2 + \frac{r_i^2 r_o^2}{r^2} - \frac{1+3\nu}{3+\nu}r^2\right) \quad \text{（均一ディスク、解析解）}$$3Dモデルが必要な場合
- 翼（ブレード）のような非軸対称構造が付いている場合
- 部分的な欠陥・損傷の評価
- キー溝・ボルト穴など周方向不均一性が重要な場合
7. ジャイロ効果と歳差運動
高速回転する物体は角運動量 $\vec{L} = I\vec{\omega}$ を持つ。外部トルク（力のモーメント）が加わると、角運動量ベクトルの向きが変化し歳差運動（Precession）が起きる。
$$\vec{\tau} = \frac{d\vec{L}}{dt} = I\vec{\omega} \times \vec{\Omega}_{prec}$$直感に反するのがジャイロ効果の特徴で、コマを指で押しても倒れずに向きが変わるのはその典型だ。
自動車のホイールとジャイロモーメント
走行中の自動車のホイールは高速回転している。コーナリング（ステアリング操作）でホイール軸の向きが変わると、ジャイロモーメントがサスペンションに作用する。高回転・大径ホイール（スポーツカー）ではこのモーメントが無視できず、ハンドリング設計に影響する。FEM多体動力学解析（MBD）でジャイロ効果を含めてシミュレーションする。
8. 回転機械の振動解析とキャンベル線図
回転機械では、回転速度に比例した周波数成分が構造の固有振動数と一致したとき共振（Resonance）が生じる。この現象を体系的に管理するのがキャンベル線図（Campbell Diagram）だ。
キャンベル線図の構成
横軸に回転速度 [rpm または rad/s]、縦軸に振動数 [Hz または rad/s] をとる。
- 水平線：構造の固有振動数（回転速度によらず一定、ただし回転が速いとコリオリ効果で分岐）
- 右上がりの直線（次数線）：回転速度 $n$ の $k$ 倍の励振周波数（$kn/60$ [Hz]）。$k$ は整数（1, 2, 3...）で「次数（Engine Order）」と呼ぶ
- 交点：固有振動数と励振周波数が一致する共振点。運転域内にあると危険
臨界速度（Critical Speed）
回転軸の横振動の固有振動数と一致する回転速度を臨界速度と呼ぶ。これを超えて安定に運転する（超臨界運転）設計も多い——蒸気タービン軸や高速遠心分離機では、臨界速度を素早く通過して超臨界域で運転する。
9. 実践：ターボチャージャーの破裂解析
ターボチャージャーのコンプレッサーホイールは、過速度（規定の最大回転数の120〜150%）での破裂強度を試験・解析で確認することが安全規格（ISO 22728 等）で要求される。
破裂の物理
コンプレッサーホイールの最大遠心応力（フープ応力）が材料の引張強さ $\sigma_B$ に達すると、ホイールは急激に崩壊する。アルミ合金ホイール（A2024-T6, $\sigma_B \approx 480\,\text{MPa}$）の場合、遠心応力が $\sigma_B$ を超える回転数が「破裂回転数」だ。
FEM解析手順
- 静解析：遠心荷重 + 温度分布（タービン側からの熱伝導）を入力
- 弾塑性解析：大変形・温度依存材料特性を考慮した非線形FEM
- 破裂回転数の推定：フープ応力が $\sigma_B$ に達する回転数を求める
- 破片の運動解析：破裂後の破片がハウジングを貫通しないかの評価（Explicit FEM）
熱応力との複合荷重
排気ガスで加熱されるタービン側ホイール（最高 1,000℃超）は熱膨張と遠心力の複合荷重を受ける。温度依存の降伏応力（高温では著しく低下）と熱膨張による寸法変化を同時考慮する「熱構造連成解析」が不可欠だ。