🙋
転がり接触応力って、軸受とか歯車でよく聞くけど、具体的に何が計算できるんですか？
🎓
大まかに言うと、ボールとレースが押し合った時に、どれくらいの小さな面積で力を受け持っているか、その部分の圧力はどれくらいか、を計算するんだ。このシミュレーターでは、右側のパラメータで「物体の半径」や「荷重」を変えると、リアルタイムで接触の広さ（接触半径a）と最大圧力（p₀）が計算されてグラフが変わるよ。まずは荷重Pのスライダーを動かしてみて。
🙋
え、荷重を増やすと接触面積が広がるけど、最大圧力も上がるんですね。でも、表面が一番圧力が高いのではないの？
🎓
そこが転がり接触の面白いところなんだ。表面の直下、だいたい接触半径の半分くらいの深さ（z=0.47a）で、せん断応力が最大になるんだよ。下の「表面下応力分布」グラフで、z方向の深さを表すスライダーを動かすと、その深さでの応力の大きさがわかる。疲労き裂はこの最大せん断応力の位置から始まることが多いんだ。
🙋
材料を「鋼 vs 鋼」から「鋼 vs ゴム」に変えると、結果がすごく変わります！応力が大きく異なる。これはどういうことですか？
🎓
いいところに気づいたね。材料の組み合わせで「等価弾性率E*」が変わるからなんだ。ゴムみたいに柔らかい材料が相手だと、同じ荷重でも接触面積が大きく広がって、圧力が分散される。だから最大接触圧p₀が劇的に下がるんだ。実務では、クッション材を挟んで衝撃を和らげる設計に応用されているよ。
本ツールは2つの弾性球体の接触を前提としています。円筒や平面との接触には対応していませんが、等価半径R*の計算式を適宜変更することで、球-平面接触（R₂=∞）など一部のケースを近似できます。
画面上の材料パラメータ入力欄で、各球体のヤング率E₁、E₂とポアソン比ν₁、ν₂を直接数値入力できます。入力後、等価弾性率E*が自動計算され、接触半径や最大接触圧に即座に反映されます。
接触面下の最大せん断応力は深度約0.48a（aは接触半径）で発生します。本ツールでは深度0～2aの範囲で応力分布をプロットし、最大せん断応力の値とその深度を数値で表示します。
本ツールはヘルツ理論に基づく弾性接触の厳密解を計算するため、材料が等方性で弾性範囲内、かつ表面が滑らかであれば高精度です。ただし、塑性変形や摩擦、表面粗さの影響は考慮していないため、実機評価の参考値としてご利用ください。
転がり軸受：ボールやローラーとレース道の接触応力を計算し、基本定格寿命（L10寿命）を予測する基礎となります。荷重条件や材料の組み合わせを変えて、最適な軸受サイズを選択するために使われます。
歯車：歯面同士の接触（歯面接触）において、発生する接触圧力（ヘルツ応力）を計算します。これが歯面の疲労強度（ピッチング強度）を評価するための重要な指標となり、歯車のモジュールや歯幅を決定する設計パラメータになります。
レールと車輪：鉄道車輪がレールヘッドに接触する部分の応力を解析します。非常に大きな荷重が繰り返されるため、接触面下のせん断応力分布を評価することで、レールのき裂進展や剥離のメカニズムを理解し、保守間隔を計画するのに役立ちます。
カム・フォロワ：エンジンのカムシャフトとロッカーアームやフォロワーの接触部の応力を計算します。高速で繰り返される接触荷重に対して、十分な耐久性を持つように材料や表面処理（浸炭焼入れなど）を選定する設計プロセスで活用されます。
このツールを使い始める際、いくつかつまずきやすいポイントがあるから気をつけてね。まず第一に、「接触半径が大きいから安全」と単純に考えないこと。確かに接触面積は広がるけど、最大接触圧p₀は荷重Pの増加に対して、$a$の2乗に反比例して急激に上がるんだ。例えば、荷重を8倍にすると接触半径は2倍になるけど、最大接触圧は4倍になる。面積が4倍になっても、単位面積あたりの負担は増えてしまうわけ。だから、荷重増加の影響は思ったより深刻だということだね。
次に、材料定数の入力だ。特にポアソン比νを軽視しがちだけど、これが等価弾性率E*に効いてくる。例えば、鋼（ν=0.3）とゴム（ν≈0.5）では、νの2乗の項$(1-\nu^2)$の値が大きく変わる。実務でよくあるのが、材料メーカーのカタログ値と実際の材質の微妙な違い。特に樹脂や複合材料はバッチ差があるから、安全側を見て少し大きめの弾性率を使うなどの工夫が必要だよ。
最後に、この計算は完全な弾性体かつ表面が滑らかという理想条件が前提だということを忘れてはいけない。実際の部品には表面粗さや潤滑油の影響が大きく、それらは応力分布を変える。計算結果は「目安」として捉え、特に疲労寿命を評価する時は、ここで求めたτ_maxに安全率をかけるか、実機テストで検証するのが鉄則だ。