🧑🎓
このシミュレーターで計算できる「転がり接触疲労」って何ですか？軸受が壊れる原因ですか？
🎓
その通り！ざっくり言うと、ボールベアリングや歯車が「転がりながら押し合う」ことで、表面のすぐ下で疲労が溜まって、小さな穴（ピッティング）がポコポコできる現象だ。例えば、自動車のホイールベアリングが「カラカラ」と異音を出すのは、これが原因のことが多いんだ。このツールでは、接触面の圧力や、疲労の元になる最大せん断応力がどこで起こるかを計算できるよ。
🧑🎓
え、表面の下で壊れるんですか？じゃあ、左のグラフで赤く表示されている山が接触圧力で、その下の青い線がせん断応力ってことですか？
🎓
鋭いね！その解釈で正解。表面の圧力はヘルツ接触で放物線状に分布する。でも、材料をずらそうとする「せん断応力」は表面直下、深さ0.48aのところで最大になる。シミュレーターで「法線力F」のスライダーをグイッと大きくしてみて。接触幅aが広がると同時に、最大せん断応力の位置（グラフの青丸）も深く移動するのがわかるだろう？これが設計で超重要なポイントなんだ。
🧑🎓
なるほど！で、右側に出てくる「L10寿命」って、どうやって求めてるんですか？信頼性90%の寿命ってことですよね？
🎓
そう、90%が壊れない寿命だ。これは実務で非常によく使われる基準で、軸受メーカーのカタログ寿命計算の基礎になってる。計算式は、最大接触圧力$p_0$と材料の耐久限界から求めるんだ。ツールでは「弾性係数E₁, E₂」や「半径R₁, R₂」を変えると、接触状態が変わり、寿命がガラッと変わるのが体感できる。硬い材料を使ったり、曲率を大きくすると寿命が伸びる理由がわかるよ。
二つの円柱（または球）が接触するときの、接触幅と最大接触圧力を求めるヘルツ接触の基本式です。等価半径$R^*$と等価弾性係数$E^*$を使って計算を簡略化します。
$$a = \left(\frac{3FR^*}{4E^*}\right)^{1/3}, \quad p_0 = \frac{3F}{2\pi a^2}$$
$a$: 接触半幅 [mm]
$p_0$: 最大接触圧力 [MPa]
$F$: 法線方向の荷重 [N]
$R^*$: 等価半径 ($1/R^* = 1/R_1 + 1/R_2$) [mm]
$E^*$: 等価ヤング率 ($1/E^* = (1-\nu_1^2)/E_1 + (1-\nu_2^2)/E_2$) [GPa]
転がり接触疲労の起点となる、接触面直下に発生する最大せん断応力$\tau_{max}$とその発生深度$z_{max}$を求める式です。この応力が材料の耐久限界を超えると、微小き裂が発生します。
$$\tau_\text{max}\approx 0.31\,p_0 \quad \text{at depth}\quad z_\text{max}\approx 0.48\,a$$
$\tau_{max}$: 最大せん断応力 [MPa]。圧力分布の内部で発生し、材料の降伏や疲労き裂の起点となる。
$z_{max}$: 最大せん断応力の発生深度 [mm]。表面ではなく、少し下がった位置で最大となるのがポイント。
よくある誤解と注意点
この手の計算で最初にハマるポイントをいくつか挙げておくよ。まず、「接触幅aが大きければ大きいほど安全」という誤解。確かに接触圧力p0は下がるけど、最大せん断応力τ_maxの発生深度z_maxも深くなるんだ。例えば、表面硬化処理をした歯車の場合、硬化層の深さがこのz_maxより浅いと、一番応力が高い場所が軟らかい芯材部分になってしまい、かえって早期破損の原因になる。ツールでFを大きくしてz_maxの動きを確認してみると、この関係がよくわかるはずだ。
次に、材料パラメータの入力ミス。特にポアソン比νは0.3前後で適当に入れがちだけど、これが等価ヤング率E*を計算する式 $1/E^* = (1-\nu_1^2)/E_1 + (1-\nu_2^2)/E_2$ に入っている。νを0.25と0.33で変えただけでも、E*は数%変わり、寿命L10にはさらに大きな影響が出る。実務では材料証明書の値を必ず確認しよう。
最後に、この計算は「理想的な状態」の基礎値であることを忘れないで。実際の軸受や歯車には、潤滑油の影響、表面粗さ、残留応力、組み付け誤差など、計算に入っていない要素が山ほどある。例えば、ツールで計算した寿命が10,000時間でも、潤滑不良があればその1/10以下になることも珍しくない。このシミュレーション結果は「比較のためのベンチマーク」と捉え、安全率をしっかり見込んだ設計を心がけよう。
関連する工学分野
このヘルツ接触と疲労寿命の計算は、実は思っているよりずっと広い分野で顔を出すんだ。まず真っ先に挙がるのはトライボロジーだね。転がり接触面には必ず潤滑油膜が介在する。この油膜圧力が接触圧力分布を変え、最大せん断応力の値を軽減したり、逆に表面ピッチングを促進したりする。このツールで求めた接触圧力は、弾流体潤滑（EHL）解析の重要な入力条件になるんだ。
もう一つは材料表面工学との深い関係。さっきも少し触れたけど、ショットピーニングや浸炭焼入れなどの表面改質技術は、まさに最大せん断応力が発生する深さ（z_max≈0.48a）をターゲットにしている。表面に圧縮の残留応力を作り込むことで、実際に働くせん断応力を打ち消し、疲労寿命を飛躍的に伸ばすことができる。この計算は、最適な硬化層の深さを決めるための第一歩なんだ。
さらに視野を広げると、生体力学の分野でも同じ式が活躍している。人工関節、特に股関節のボールとソケットの接触は、まさにヘルツ接触そのもの。ここでは金属やセラミックの耐久性だけでなく、接触圧力が骨の吸収や溶解を引き起こさないかどうかも評価される。一見、機械設計と無縁に見える分野が、実は同じ物理原理で繋がっているんだ。
発展的な学習のために
このツールの計算に慣れてきたら、次は「なぜそうなるのか」の背景を学ぶと一気に世界が広がるよ。まずは応力テンソルと主せん断応力の概念を押さえよう。ツールで表示されているτ_maxは、接触体内の無数の点の中で最も大きなせん断応力値だ。この値は、三つの主応力（σ1, σ2, σ3）から $ \tau_{max} = ( \sigma_1 - \sigma_3 ) / 2 $ で計算される。ヘルツ接触場では、これが深さ0.48aで発生する、というのが理論的に導かれるんだ。
次におすすめなのは、疲労寿命計算の別アプローチを調べてみること。このツールで使っているのは、最大接触圧力p0に基づく比較的簡便な方法（Lundberg-Palmgren理論の基礎）。より詳細に評価するには、「応力-寿命（S-N）曲線」を材料ごとに用意し、累積損傷則（マイナー則）で変動荷重を考慮したり、あるいは「き裂進展力学」に基づいて、初期き裂から破壊に至るまでのサイクル数を計算する方法がある。例えば、$ da/dN = C(\Delta K)^m $ のようなパリス則だね。
最終的には、この計算を多次元化するステップに進もう。このツールは2次元（平面ひずみ状態）の円柱接触がモデルだ。実際のボールベアリングは3次元の点接触。その場合、接触領域は円形になり、最大せん断応力の発生深度も変わる（約0.48aから約0.78aへ）。また、歯車のように荷重位置が移動する「走行接触」を扱うには、時間変化する応力場を追跡する必要がある。まずは今のシンプルなモデルを完全に理解することが、それら複雑な問題への最高の足がかりになるんだ。