ヘルツ接触力学シミュレーター 戻る
接触力学

ヘルツ接触力学シミュレーター

球・円筒・平板の組み合わせを選択し、荷重・材料・寸法を入力するだけで接触半径・最大圧力・沈込み量・表面下応力分布を瞬時に計算。

パラメータ入力


一時停止中はスライダーを動かすと結果が即座に更新されます。

ライブ計算結果(荷重スイープ中)
法線荷重 P [N]
接触半径 a [mm]
最大接触圧 p₀ [MPa]
沈込み量 δ [μm]
降伏判定 p₀ / σy(1.60超で初期降伏)
球-平板ヘルツ接触の可視化(荷重 P をスイープ)
理論・主要公式

等価弾性係数:

$$\frac{1}{E^*}= \frac{1-\nu_1^2}{E_1}+ \frac{1-\nu_2^2}{E_2}$$

球接触 (R* = R等価曲率半径):

$$a = \left(\frac{3FR^{*}}{4E^{*}}\right)^{1/3},\quad p_0 = \frac{3F}{2\pi a^2}$$

円筒接触 (単位長さ当たり):

$$b = \sqrt{\frac{4F R^{*}}{\pi E^{*} L}},\quad p_0 = \frac{2F}{\pi b L}$$

ヘルツ接触力学とは

🙋
ヘルツ接操作して何ですか?ベアリングのカタログで見たことあるけど、具体的に何がわかるんですか?
🎓
大まかに言うと、球や円筒が押し合った時に、どれくらいの面積で接触して、どれくらいの圧力がかかるかを計算する理論だよ。例えば、ボールベアリングの玉とレースの間の圧力分布を予測して、寿命を計算するのに使われるんだ。このシミュレーターで「接触形式」を「球-球」に、法線力Fを変えてみると、接触半径aがどう変わるかすぐにわかるよ。
🙋
え、接触圧力って表面で一番大きいのではないんですか?でも「最大せん断応力」って項目がありますよね?
🎓
そこがヘルツ接触の面白いところ!実は最大のせん断応力は表面ではなく、材料の内部で発生するんだ。球が接触する場合、表面下だいたい接触半径の半分の深さ(z≈0.48a)で最大になる。これが材料の降伏(塑性変形)の始まりになることが多いんだよ。シミュレーターで「降伏応力σ_y」を設定して、最大接触圧力p0と比べてみると、安全率がわかる。
🙋
なるほど!じゃあ、材料の組み合わせで結果は大きく変わるんですか?「等価弾性係数E*」って何を表しているんですか?
🎓
その通り!鋼の玉とゴムの板では変形が大きく異なるよね。E*は、2つの異なる材料の接触を、1つの仮想的な材料と1つの曲率半径を持つ接触に置き換えるための「代表的な硬さ」みたいなものだ。シミュレーターのE₁, E₂, ν₁, ν₂を動かしてみると、E*が変わり、接触半径aや沈み込み量δがどう影響されるか体感できるよ。実務では、硬い材料と柔らかい材料の組み合わせを最適化するのに使うんだ。

よくある質問

形状の選択肢に「球 — 平板」があります。これを選択すると平板側の曲率半径は自動的に無限大として扱われ、R2の入力は不要になります。円筒と平面の場合も同様に「円筒 — 平板」を選択してください。
表面下の最大せん断応力の発生深さが分かります。この位置が材料の降伏点を超えると、表面下から亀裂が発生する可能性があります。軸受や歯車の疲労寿命評価に役立ててください。
荷重の単位はニュートン(N)です。本シミュレーターは静的な垂直荷重のみを対象としており、動的荷重や摩擦を伴う接線力は考慮していません。純粋な押し付け接触の解析にご利用ください。
剛体として扱いたい材料のヤング率に非常に大きな値(例:1e15 MPa)、ポアソン比に0.3を入力してください。等価弾性係数の計算でその材料の項が実質的にゼロとなり、剛体接触として正しく計算されます。

実世界での応用

転がり軸受(ベアリング):ボールやローラーとレース道との接触圧力を計算し、疲労寿命(L10寿命)を予測するために不可欠です。材料や曲率半径を最適化して、高い荷重能力と長寿命を実現します。

歯車の歯面接触:かみ合う歯面は瞬間的に線接触(円筒接触に近似)します。ここに発生するヘルツ応力から、歯面のピッチング(表面疲労)や焼き付きのリスクを評価します。

鉄道の車輪とレール:円筒と平面の接触に近い状態です。巨大な荷重が一点に集中するため、接触圧力とそれに伴う内部せん断応力が、レールのスケーリング(鱗状剥離)やき裂進展の起点となります。

プローブピンやコンタクトレンズ:微小な領域での接触を評価します。例えば半導体テストのプローブカードでは、ピン先端とパッドの接触圧力が電気的接信性とパッド損傷に直結するため、設計パラメータの重要な指標です。

よくある誤解と注意点

まず、「接触半径が荷重に比例する」と思っていませんか? 実は、球接触では $a \propto F^{1/3}$ の関係で、荷重を2倍にしても接触半径は約1.26倍にしかなりません。逆に、最大接触圧力 $p_0$ は $F^{1/3}$ に比例して増加します。つまり、荷重を増やすと接触面積は少ししか広がらず、圧力が急激に高くなるのがポイントです。例えば、荷重を8倍にすると接触半径は2倍ですが、最大接触圧力は2倍になります。

次に、材料定数の入力ミス。特にポアソン比 $\nu$ の単位は無次元ですが、0.3を0.03と入力してしまうと、等価弾性係数 $E^*$ が大きく変わって計算結果が大きく狂います。また、鋼のヤング率は約206 GPaですが、単位を「MPa」のまま「206000」と入力するか、「GPa」で「206」と入力するか、ツールの仕様を必ず確認しましょう。単位系の不一致は最も多いエラーの元です。

最後に、ヘルツ理論の適用限界。この理論は「完全な弾性体」「接触面が滑らか」「摩擦がない」という理想条件が前提です。実際には、表面粗さや潤滑油膜、材料の塑性変形が発生すれば、計算結果は現実からズレます。例えば、計算上の最大接触圧力 $p_0$ が材料の降伏応力 $\sigma_y$ を超えていても、局所的な塑性変形で応力が緩和される場合があります。シミュレーターの結果は「第一近似」として捉え、重要な設計ではより詳細なCAE解析で検証するのが鉄則です。

使い方ガイド

  1. R1、R2に接触する2つの部材の曲率半径(mm)を入力。平板との接触は「球 — 平板」「円筒 — 平板」を選択(平板側は自動的に∞として扱われ、R2の入力は不要)
  2. E1、E2に各部材のヤング率(GPa)を入力。鋼は200、アルミは70、銅は110が目安
  3. nu1、nu2にポアソン比を入力。鋼は0.3、アルミは0.33、ゴムは0.5が一般的
  4. Fに接触荷重(N)、Llenに接触長さまたは幾何学的長さ(mm)を入力して計算実行

具体的な計算例

鋼球同士(R1=R2=10mm、E1=E2=200GPa、nu1=nu2=0.3)に荷重F=1000Nを加えた場合、接触半径a≈0.32mm、最大ヘルツ接触圧力≈4540MPa、押込み量δ≈0.021mmと計算されます。点接触では高荷重時に接触圧が数GPaに達するため、材料の許容接触圧との照合が設計の要点です。円筒ローラと平板の接触(「円筒 — 平板」、R1=5mm)にF=5000N、L=50mmを加えると接触半幅b≈0.076mm、最大圧力≈836MPa(線荷重100N/mm)が得られ、圧延機やコンベアベアリング設計に活用できます。

準拠規格・前提条件

準拠/参考: ヘルツ接触理論(Johnson『Contact Mechanics』1985)。等価弾性 \(1/E^* = (1-\nu_1^2)/E_1 + (1-\nu_2^2)/E_2\)。球:\(a = (3FR^*/4E^*)^{1/3},\; p_0 = 3F/(2\pi a^2),\; \delta = a^2/R^*\)。円筒(単位長):\(b = \sqrt{4FR^*/(\pi E^* L)},\; p_0 = 2F/(\pi b L)\)。

モデルの前提: 線形弾性・等方・滑面・無摩擦・静的法線荷重のみ。微小変形で接触域が代表寸法より十分小さい。最初の降伏は表面下(球 z≈0.48a)で発生し、von Mises基準で \(p_0/\sigma_y \approx 1.60\) が初期降伏の目安。\(\tau_{max}\approx 0.31\,p_0\)(球)。

適用範囲・限界: 弾性域の第一近似。塑性・粘着・接線力(摩擦)・潤滑膜・表面粗さ・動荷重は対象外。EHL条件下では圧力分布が理論から数十%ずれうる。重要設計では詳細CAEで検証のこと。

実務での注意点

  1. 接触圧力が材料の降伏応力を超えないか確認が必須。焼き入れ鋼は1500MPa、靱性鋼は800MPa程度が上限目安
  2. 繰り返し接触による疲労を考慮し、許容接触応力は静的圧力の60~70%に設定する慣例がある
  3. 潤滑膜厚さと表面粗さが接触挙動に影響。EHL条件下では圧力分布が理論値から30%程度変動する可能性
  4. 温度上昇によるヤング率低下。接触部温度が100℃上昇するとE値は1~3%減少し、接触圧力は3~5%増加する傾向