理論メモ
中実ディスク(r=0)の応力:$$\sigma_r = \frac{3+\nu}{8}\rho\omega^2(R^2 - r^2)$$ $$\sigma_\theta = \frac{3+\nu}{8}\rho\omega^2 R^2 - \frac{1+3\nu}{8}\rho\omega^2 r^2$$ 環状ディスクは Lamé 定数 A, B を内外境界条件で決定。
外径・内径・回転速度・材料を変えて、ディスクの半径方向応力 σr と円周方向応力 σθ の分布をリアルタイムに可視化。破裂速度と安全率も自動算出します。
回転する薄肉円盤(円板)の半径方向と円周方向の応力分布は、Laméの方程式で記述されます。中実円盤(内径r=0)の場合の解は以下の通りです。
$$\sigma_r = \frac{3+\nu}{8}\rho\omega^2(R^2 - r^2)$$ $$\sigma_\theta = \frac{3+\nu}{8}\rho\omega^2 R^2 - \frac{1+3\nu}{8}\rho\omega^2 r^2$$ここで、$\sigma_r$: 半径方向応力 [Pa], $\sigma_\theta$: 円周方向応力(フープ応力)[Pa], $\rho$: 材料密度 [kg/m³], $\omega$: 角速度 [rad/s], $R$: 円盤外半径 [m], $r$: 円盤内任意点の半径 [m], $\nu$: ポアソン比。中心(r=0)では両応力が等しく最大値となります。
内孔を持つ環状ディスクの場合は、一般解の積分定数A, Bを、内径と外径で半径方向応力が0という境界条件(自由端)から決定します。
$$\sigma_r = A - \frac{B}{r^2}- \frac{3+\nu}{8}\rho\omega^2 r^2$$ $$\sigma_\theta = A + \frac{B}{r^2}- \frac{1+3\nu}{8}\rho\omega^2 r^2$$境界条件 $\sigma_r = 0$ at $r=r_i, R$ を代入してA, Bを求めます。この解では内孔縁($r=r_i$)で円周方向応力$\sigma_\theta$が最大となり、中実円盤より高い応力集中が生じることが特徴です。
タービン・発電機ディスク:ガスタービンや蒸気タービンの回転ディスクは、高温高応力下で作動します。フープ応力を正確に評価し、クリープ破壊や低サイクル疲労を防ぐために、この解析が設計の基本となります。
遠心分離機・遠心鋳造機:高速回転するロータ(ボウル)は、中実または小さな内孔を持つ形状が一般的です。材料の降伏応力と破裂速度を比較することで、運転速度の安全限界を設定します。
自動車のブレーキディスク・クラッチ板:回転部品として遠心力による応力が発生します。特にベンチレーテッドディスクなど複雑な形状でも、基本的な円盤応力の理解が熱応力との組み合わせ解析の基礎になります。
フライホイールエネルギー蓄積装置:高速回転により運動エネルギーを蓄えるフライホイールは、軽量化と高強度が求められます。繊維強化複合材料を用いる場合も、円周方向の引張応力が設計上のキーとなります。
このツールを使い始める際、特にCAE初心者が陥りがちな落とし穴がいくつかあります。まず、「材料の降伏応力さえ超えなければ安全」という誤解です。確かに破裂速度は降伏応力から計算しますが、実機では疲労破壊が大きな問題になります。例えば、毎分10万回転を繰り返す遠心分離機のロータは、最大応力が降伏応力の50%以下でも、繰り返し荷重により微小なき裂が成長して破断することがあります。ツールで安全率を確認したら、次は疲労寿命の検討が必要です。
次に、パラメータ入力の単位ミス。これは本当に多いです。ツールは内部でSI単位(m, kg, rad/s)で計算していますが、設計図面ではmmで、回転速度はrpmで記載されていることがほとんど。例えば、外径100mmを入力する際、うっかり「100」のまま入力すると(正しくは0.1)、応力は実は1万分の1で計算され、とんでもなく安全そうに見える結果が出て大事故のもとになります。入力後は、常識的なオーダーの応力値(例えば鋼なら数十~数百MPa)が出ているか、必ずサニティチェックを。
最後に、「薄肉円盤」理論の限界を理解すること。このツールの基礎式はディスクの厚さが一定で薄いことが前提です。しかし、実際のタービンディスクはボス部やブレード取付部で厚みが大きく変化します。そのような複雑形状では、このツールの結果はあくまで一次評価。ツールで内孔縁の応力集中を確認したら、その値を参考に、より詳細な3次元FEM解析で形状最適化を行うのが実務の流れです。
この回転体応力解析の考え方は、遠心機やタービン以外の様々な工学分野でも応用されています。まず、高速回転機械全般です。例えば、自動車のターボチャージャーインペラーや、CNC工作機械のスピンドル、さらにはハードディスクドライブ(HDD)内のプラッタも高速回転体です。HDDの場合、プラッタの変形が微小でもヘッドの浮上量に影響するため、変形解析が極めて重要になります。
また、慣性力を利用した装置の設計にも直結します。遠心鋳造機では、溶融金属を型の中で回転させて遠心力で押し付けるため、型自体の強度計算に同じ理論が使えます。同じく、回転式の遠心分離・遠心沈降装置では、ロータの強度が最大処理G(遠心加速度)を決定し、分離性能の上限を左右します。
少し視点を変えると、材料科学や実験力学との接点もあります。遠心機は材料に体積力を加える優れた実験装置です。例えば、コンクリートや地盤材料の遠心模型実験では、このツールで計算するような応力が模型内に再現され、実物規模の挙動を調べられます。つまり、回転体の設計者だけでなく、回転そのものを「道具」として使う研究者にも必要な基礎知識なのです。
このツールの背後にある理論を深め、実務に活かすための学習ステップを紹介します。まず第一歩は、基礎式の「意味」を図解で理解することです。「ラメの方程式」という名前でビビらず、力のつり合い(半径方向)と変形の適合条件(円周方向)からどのように導かれるか、を追ってみてください。多くの機械力学の教科書に、微小要素に働く力を描いた図とともに解説があります。
次に、「厚肉円筒」の理論に進むことをお勧めします。これは、このツールの対象である「薄肉円盤」に軸方向の厚みがある場合や、内圧・外圧が同時にかかるシリンダーの計算に発展します。基本式は非常によく似ており、$$\sigma_\theta = \frac{p_i r_i^2 - p_o r_o^2}{r_o^2 - r_i^2} + \frac{(p_i - p_o) r_i^2 r_o^2}{(r_o^2 - r_i^2) r^2}$$ のような形で表されます(p: 圧力, ri, ro: 内径・外径)。回転と内圧が組み合わされた、ポンプのインペラーなどの解析に必要です。
最終的には、有限要素法(FEM)による回転体解析を自分の手で行えるようになるのが目標です。FEMソフトでは「スピニング」や「遠心力」荷重を簡単に設定できます。このツールで得た「内孔縁で応力最大」といった直感を、複雑な3D形状でどこまで適用できるか、FEM結果と比較検証してみると、理解が一気に深まります。CAEエンジニアとしての強力な武器になる分野です。