すべり軸受シミュレーター 戻る
トライボロジー

すべり軸受シミュレーター

短軸受近似レイノルズ方程式でジャーナル軸受の油膜圧力分布・負荷容量・最小油膜厚さをリアルタイム計算。ゾンマーフェルト数と軸受性能の関係を視覚的に理解できます。

軸受寸法

運転条件

計算結果
Sommerfeld S
偏心率 ε
h_min (μm)
P_max (MPa)
摩擦係数 f
動力損失 (W)
軸受断面図
周方向油膜圧力分布
ゾンマーフェルト数 S vs 最小油膜厚さ h_min/C
理論・主要公式

$$h(x) = h_1 + (h_2 - h_1)\frac{x}{L}$$

くさび形油膜厚さ(m):$h_1, h_2$ は両端の膜厚、$L$ はすべり面長さ(m)。

$$\frac{dp}{dx} = 6\mu U \frac{h - \bar{h}}{h^3}$$

レイノルズ方程式(1次元):$\mu$ は粘度(Pa·s)、$U$ は相対速度(m/s)、$\bar{h}$ は最大圧力位置の膜厚。

$$F_L = \int_0^L p(x)\,dx$$

流体浮力(N/m):油膜圧力分布を積分した軸受荷重支持能力。

すべり軸受シミュレーターとは

🙋
このシミュレーターで「ゾンマーフェルト数」ってよく出てきますけど、これって何ですか?設計でどう使うんですか?
🎓
大まかに言うと、軸受が「どれだけうまく浮いているか」を表す無次元の指標だね。数式は $S = \frac{\mu N L D}{W}\left(\frac{R}{C}\right)^2$ で、粘性力と荷重の比を表すんだ。Sが大きいほど(高速・高粘度・低荷重)油膜が厚く安定して、小さいと軸と軸受が接触しやすくなる。実際の設計では、このSの値で潤滑状態を予測するんだよ。上の「回転速度」や「軸受荷重」のスライダーを動かして、Sの値がどう変わるか確かめてみて。
🙋
え、そうなんですか!でも、油膜が薄すぎるとどうなってしまうんですか?「最小油膜厚さ」が画面に出てますけど。
🎓
そこが一番大事なポイントだ!最小油膜厚さ $h_{min}= C(1 - \epsilon)$ が、軸と軸受の最接近距離だ。これが表面粗さより十分に厚くないと、金属同士がぶつかって焼き付きや摩耗の原因になる。実務では、$h_{min}$ を表面粗さの3〜5倍以上に保つ設計が目安だね。シミュレーターで「ラジアルクリアランス」を小さくしたり「偏心率」を大きくすると、$h_{min}$ が一気に小さくなって危険な状態になるのがわかるよ。操作してみて。
🙋
なるほど…。でも、画面のカラフルな圧力分布って、実際の設計ではどう見るんですか?一番赤いところが危ないということ?
🎓
良いところに気が付いたね!その圧力分布こそが軸受が荷重を支える仕組みなんだ。油がくさび形の隙間に絞り込まれて高圧が生まれ、それが軸を持ち上げている。赤い高圧部の位置と大きさで、負荷容量や軸の位置が決まる。現場でトラブルが多いのは、この圧力分布が偏りすぎて一部に極端な高圧がかかる「エッジ負荷」だね。「軸受幅」のパラメータを極端に小さくすると、圧力の山が尖ってしまうのがシミュレーターでも確認できるよ。

物理モデルと主要な数式

このシミュレーターの基礎は、短軸受近似(無限幅近似)を仮定したレイノルズ方程式です。軸受円周方向の油膜圧力分布を計算します。

$$ \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 = (μN/P) × (R/C)^2 で定義されます。μは粘度、Nは回転数、Pは面圧、Rは軸半径、Cはクリアランスです。この値が大きいほど油膜が厚くなりやすく、軸受の負荷容量や安定性を評価する指標として使います。
偏心率ε(0〜1)は軸の偏心の度合いを示します。εを大きくすると最小油膜厚さが減少し、油膜圧力のピークが高くなります。負荷容量は増加しますが、油膜が薄くなるため焼付きリスクが高まります。εが0に近いほど油膜は均一で安定します。
短軸受近似は軸受幅が直径に比べて十分小さい(L/D < 0.5)場合に有効です。円周方向の圧力勾配を無視し、軸方向の流れのみを考慮します。幅が広い軸受では誤差が大きくなるため、その場合は長軸受近似や2次元解析が必要です。
粘度μまたは回転数Nを上げると、レイノルズ方程式の右辺(くさび効果)が大きくなり、油膜圧力が増加します。その結果、同じ偏心率でも負荷容量が向上します。ただし、過度に上げると発熱やキャビテーションのリスクが生じるため、適切な範囲で調整してください。

実世界での応用

大型発電用タービン・発電機:高速回転する巨大な回転子を支えるためにすべり軸受が用いられます。シミュレーションでゾンマーフェルト数や最小油膜厚さを確認し、起動時・停止時を含む全ての運転条件で安定した流体潤滑状態を確保することが、長期信頼性に直結します。

船舶用プロペラシャフト軸受:海水による腐食環境かつ大きな推力を支えるため、白金属などのすべり軸受が使われます。軸受にかかる荷重とプロペラの回転数から油膜圧力分布を予測し、軸受の幅やクリアランスを最適化するのに本ツールのような解析が役立ちます。

産業用モーター・ポンプ:工場で多用される中規模モーターの軸受設計では、コストと性能のバランスが重要です。異なる粘度の潤滑油を使用した場合の摩擦係数と負荷容量を比較し、エネルギー損失を抑えつつ十分な耐久性を持つ設計を探ります。

自動車エンジンのクランクシャフト軸受:燃焼圧力による衝撃荷重と高速回転が同時に作用する過酷な環境です。軸受の幅やクリアランスを変えた時の油膜圧力のピーク値と最小油膜厚さを評価し、疲労破壊や焼き付きを防ぐ設計検証に活用されます。

よくある誤解と注意点

このシミュレーターを使い始めるとき、いくつかつまずきやすい落とし穴があるよ。まず第一に、「ゾンマーフェルト数Sが大きければ常に安全」と思いがちだけど、そう単純じゃない。確かにSが大きいと油膜は厚くなるけど、発熱が無視できなくなるんだ。高速・高粘度でSを大きくしすぎると、油のせん断による発熱で粘度が下がり(サーマルランアウェイ)、かえって油膜が薄くなる悪循環に陥ることもある。例えば、回転速度を1000rpmから10000rpmに上げると、発熱はほぼ10倍になる。冷却設計は必須だね。

次に、パラメータの入力順序。いきなり「偏心率」を大きく設定して「計算」を押すと、すぐに接触状態を示す結果が出てしまう。実務的な設計手順は、まず目標荷重と回転数で必要な粘度やクリアランスを決め、その状態での偏心率や最小油膜厚さを結果として確認する流れだ。「偏心率」は原因ではなく、荷重と剛性の結果として生じる現象と捉えるのがコツだよ。

最後に、このシミュレーターの大きな前提「無限幅近似」を見落とさないで。現実の軸受は側面から油が漏れる「端洩れ」があるから、計算された負荷容量は実際より過大評価されがち。特に幅が狭い(L/D < 0.5くらい)軸受では要注意。シミュレーターでL/Dを0.3と1.0で比べてみて。圧力分布の山の形が、幅広だと台形、幅狭だと尖った三角形になるのがわかるはず。これが端洩れの影響のイメージだ。

使い方ガイド

  1. 軸受半径R(mm)、軸受幅L(mm)、軸受すきまC(μm)を入力。鋼製ジャーナル軸受の場合、典型値はR=25mm、L=50mm、C=50μmです
  2. 回転数N(rpm)を設定。電動機軸受で3000rpm、工作機械主軸で6000rpm以上が一般的。粘度を考慮してゾンマーフェルト数So=(μ・N・R²)/(P・C²)を自動計算
  3. 短軸受レイノルズ方程式を有限差分法で離散化し、油膜圧力分布・最小油膜厚さδmin・負荷容量Wを同時リアルタイム計算。圧力分布の2次元コンター図と軸受性能曲線を表示

具体的な計算例

R=20mm、L=40mm、C=45μm、N=2800rpm、ISO VG 32鉱物油(μ=32mPa・s)のジャーナル軸受を想定。負荷P=5kNのとき、ゾンマーフェルト数So≈0.95。レイノルズ方程式解析により、最小油膜厚さδmin=8.2μm、最大圧力pmax=2.8MPa、負荷容量W=6.2kNが算出されます。回転数を5600rpmに増加させるとSo=1.90となり、δmin=15.1μm、pmax=1.5MPaへ変化。同一負荷で回転速度の上昇が油膜厚さを増加させるメカニズムを可視化

実務での注意点

  1. すきまCは製造誤差と熱膨張を考慮。アルミ軸受座の場合、温度上昇時にすきまが増加してSoが低下し、油膜破壊リスクが急増するため、実装荷重に20%の安全係数を設ける
  2. 粘度入力値は軸受温度における動粘度を使用。常温32cSt油でも80℃では8cSt程度に低下するため、定常運転温度を予測してシミュレーション条件を設定
  3. 最小油膜厚さがすきまの0.2倍未満(C=50μmなら10μm以下)になると、境界潤滑領域に移行して摩耗が急速に進行。警告値を超えたら軸受交換または運転条件見直しが必須