軸受寸法
運転条件
Sommerfeld S
—
偏心率 ε
—
h_min (μm)
—
P_max (MPa)
—
摩擦係数 f
—
動力損失 (W)
—
軸受断面図
周方向油膜圧力分布
ゾンマーフェルト数 S vs 最小油膜厚さ h_min/C
短軸受近似レイノルズ方程式でジャーナル軸受の油膜圧力分布・負荷容量・最小油膜厚さをリアルタイム計算。ゾンマーフェルト数と軸受性能の関係を視覚的に理解できます。
このシミュレーターの基礎は、短軸受近似（無限幅近似）を仮定したレイノルズ方程式です。軸受円周方向の油膜圧力分布を計算します。
$$ \frac{d}{d\theta}\left( h^3 \frac{dp}{d\theta}\right) = 6 \mu U R \frac{dh}{d\theta}$$ここで、$h = C(1 + \epsilon \cos\theta)$ は油膜厚さ、$p$は油膜圧力、$\mu$は油の粘度、$U$は軸表面速度、$R$は軸半径、$\theta$は円周角度、$C$はラジアルクリアランス、$\epsilon$は偏心率です。右辺の $dh/d\theta$ が「くさび効果」を表し、これが圧力を発生させる源です。
レイノルズ方程式を解いて得られる圧力分布を積分すると、軸受が支えられる荷重（負荷容量）と摩擦係数が求まります。これが設計の直接的な指標となります。
$$ W = L \int_{0}^{2\pi}p R \cos\theta \, d\theta, \quad f = \frac{F_f}{W} $$$W$は負荷容量、$L$は軸受幅、$F_f$は粘性抵抗による摩擦力、$f$は摩擦係数です。ゾンマーフェルト数 $S$ は、これらの結果を無次元化して整理する際に現れる重要なパラメータです。
大型発電用タービン・発電機：高速回転する巨大な回転子を支えるためにすべり軸受が用いられます。シミュレーションでゾンマーフェルト数や最小油膜厚さを確認し、起動時・停止時を含む全ての運転条件で安定した流体潤滑状態を確保することが、長期信頼性に直結します。
船舶用プロペラシャフト軸受：海水による腐食環境かつ大きな推力を支えるため、白金属などのすべり軸受が使われます。軸受にかかる荷重とプロペラの回転数から油膜圧力分布を予測し、軸受の幅やクリアランスを最適化するのに本ツールのような解析が役立ちます。
産業用モーター・ポンプ：工場で多用される中規模モーターの軸受設計では、コストと性能のバランスが重要です。異なる粘度の潤滑油を使用した場合の摩擦係数と負荷容量を比較し、エネルギー損失を抑えつつ十分な耐久性を持つ設計を探ります。
自動車エンジンのクランクシャフト軸受：燃焼圧力による衝撃荷重と高速回転が同時に作用する過酷な環境です。軸受の幅やクリアランスを変えた時の油膜圧力のピーク値と最小油膜厚さを評価し、疲労破壊や焼き付きを防ぐ設計検証に活用されます。
このシミュレーターを使い始めるとき、いくつかハマりやすい落とし穴があるよ。まず第一に、「ゾンマーフェルト数Sが大きければ常に安全」と思いがちだけど、そう単純じゃない。確かにSが大きいと油膜は厚くなるけど、発熱が無視できなくなるんだ。高速・高粘度でSを大きくしすぎると、油のせん断による発熱で粘度が下がり（サーマルランアウェイ）、かえって油膜が薄くなる悪循環に陥ることもある。例えば、回転速度を1000rpmから10000rpmに上げると、発熱はほぼ10倍になる。冷却設計は必須だね。
次に、パラメータの入力順序。いきなり「偏心率」を大きく設定して「計算」を押すと、すぐに接触状態を示す結果が出ちゃう。実務的な設計手順は、まず目標荷重と回転数で必要な粘度やクリアランスを決め、その状態での偏心率や最小油膜厚さを結果として確認する流れだ。「偏心率」は原因ではなく、荷重と剛性の結果として生じる現象と捉えるのがコツだよ。
最後に、このシミュレーターの大きな前提「無限幅近似」を見落とさないで。現実の軸受は側面から油が漏れる「端洩れ」があるから、計算された負荷容量は実際より過大評価されがち。特に幅が狭い（L/D < 0.5くらい）軸受では要注意。シミュレーターでL/Dを0.3と1.0で比べてみて。圧力分布の山の形が、幅広だと台形、幅狭だと尖った三角形になるのがわかるはず。これが端洩れの影響のイメージだ。
このすべり軸受の流体潤滑シミュレーションは、実はいろんな分野の基礎技術と直結しているんだ。まず真っ先に挙がるのはトライボロジーだね。摩擦・摩耗・潤滑の科学で、軸受はその代表的な応用先。ここで学んだ油膜形成の考え方は、歯車の歯面やカム・フォロアなど、全ての「動く機械要素」の接触面設計に通じる。
もう一つは回転機械の振動工学。すべり軸受の油膜は、軸を支えるだけでなく、ばねとダンパーの役割も果たしている。油膜の剛性と減衰特性は、軸の危険速度（共振回転数）や不安定振動（オイルホイップ）の発生に直結する。シミュレーターで「回転速度」を上げていくと、軸心軌道が変わっていく様子を見られるが、あれが振動特性の変化の一端を表しているんだ。
さらにCFD（数値流体力学）との関わりも深い。このツールで解いているレイノルズ方程式は、ナビエ-ストークス方程式を油膜という薄い領域に適用して大幅に簡略化したもの。CFDの専門家から見れば「すごく特殊化されたCFDソルバー」なんだ。だから、この計算の裏側を理解することは、より複雑な流れの数値シミュレーションを学ぶ第一歩にもなるよ。
もしこのシミュレーションに興味が湧いて、もっと深く知りたくなったら、次のステップを踏んでみるといい。まずは「有限長軸受」のモデルを学ぼう。このツールの基礎である「無限幅」の仮定を外し、軸方向の圧力変化も考慮するんだ。支配方程式は $$ \frac{\partial}{\partial x}\left(\frac{h^3}{\mu}\frac{\partial p}{\partial x}\right) + \frac{\partial}{\partial z}\left(\frac{h^3}{\mu}\frac{\partial p}{\partial z}\right) = 6U \frac{\partial h}{\partial x} $$ という2次元の偏微分方程式になる。これを解くには数値解析（有限差分法など）の知識が必要で、一気に難易度が上がるが、現実にぐっと近づく。
数学的な背景としては、偏微分方程式の数値解法と無次元解析の考え方を押さえたい。なぜ無次元数（ゾンマーフェルト数など）を使うのか？それは、物理的に相似な現象をグループ化し、実験結果や設計データを整理するためだ。例えば、模型実験の結果を実機にスケールアップするときの必須技術だね。
次の具体的なトピックとしては、「動圧スラスト軸受」や「静圧軸受」を調べてみよう。ここで学んだ動圧効果（くさび効果）は、回転する円盤の斜面を押し上げる力（スラスト軸受）にも応用される。一方、外部から高圧油を供給して軸を持ち上げる静圧軸受は、発熱が少なく超精密位置決めが可能で、まったく異なる原理と面白さがある。両方知ることで、軸受設計の選択肢が広がるはずだ。