接触类型
几何参数
材料（两体共用）
荷载
赫兹接触（点接触）
$$a = \left(\frac{3FR^*}{4E^*}\right)^{1/3}, \quad p_0 = \frac{3F}{2\pi a^2}$$
$$\tau_\text{max}\approx 0.31\,p_0 \quad \text{at depth}\approx 0.48\,a$$
$\frac{1}{R^*}=\frac{1}{R_1}+\frac{1}{R_2}$, $\frac{1}{E^*}=\frac{1-\nu_1^2}{E_1}+\frac{1-\nu_2^2}{E_2}$
什么是滚动接触疲劳与点蚀
🧑🎓
“滚动接触疲劳”和“点蚀”是什么？听起来好专业。
🎓
简单来说，就是轴承或齿轮滚着滚着，表面突然出现小坑（点蚀）坏掉了。这就像你反复弯折一根铁丝，它最终会断一样，只不过这是发生在材料表面下。在实际工程中，比如汽车变速箱里的齿轮，如果设计不好，运行几十万公里后表面就可能出现麻点状的剥落，导致噪音增大甚至失效。
🧑🎓
诶，真的吗？那这个小坑是从哪里开始出现的呢？是表面吗？
🎓
有趣的是，裂纹往往不是从最表面开始的！根据赫兹理论，两个滚动的球或圆柱接触时，压力会把材料往下压，在表面下方某个深度会产生最大的剪切应力。疲劳裂纹就在这个“最累”的地方萌生，然后慢慢扩展到表面，形成剥落。你可以在模拟器里改变“法向力 F”这个滑块，会看到下方的应力云图中，那个红色的最大剪应力区域的位置和大小会跟着变化。
🧑🎓
原来藏在下面啊！那工程师怎么知道一个轴承能用多久呢？有办法算吗？
🎓
当然有，这就是“L10寿命”的概念。简单理解，就是一批轴承里90%都能无故障运行的寿命。它和接触压力强相关。你可以试试看：在模拟器里把载荷F加大一倍，你会发现计算出的L10寿命会急剧缩短，可能只剩原来的十分之一！这就是为什么工程师要精打细算每一个参数，比如把滚子半径（R1, R2）设计得大一点，接触压力就会降下来，寿命就能大幅提升。
物理模型与关键公式
赫兹接触理论是分析两个弹性体在压力下接触状态的基础。它首先计算出一个“当量半径”和“当量弹性模量”，用来简化复杂的接触问题。
$$\frac{1}{R^*}= \frac{1}{R_1}+ \frac{1}{R_2}, \quad \frac{1}{E^*}= \frac{1-\nu_1^2}{E_1}+ \frac{1-\nu_2^2}{E_2}$$
$R_1, R_2$是两个接触体的曲率半径（凸为正，凹为负），$E_1, E_2$是弹性模量，$\nu$是泊松比。$R^*$和$E^*$是综合了双方属性的当量参数。
基于当量参数，可以计算出接触区域的半径$a$和最大接触压力$p_0$。压力在接触区内呈半球形分布。
$$a = \left(\frac{3FR^*}{4E^*}\right)^{1/3}, \quad p_0 = \frac{3F}{2\pi a^2}$$
$F$是法向载荷。$a$是接触圆半径，$p_0$是中心最大压力。压力越大，材料越容易疲劳。
点蚀疲劳的关键在于次表面的剪切应力。最大剪切应力$\tau_{max}$与最大接触压力$p_0$成正比，并发生在接触中心正下方的特定深度。
$$\tau_\text{max}\approx 0.31\,p_0 \quad \text{at depth} \approx 0.48\,a$$
这个公式告诉我们，导致材料“内伤”的最大剪切应力值约为表面最大压力的31%，且发生在约0.48倍接触半径的深度处。这就是疲劳裂纹的起源点。
现实世界中的应用
汽车变速箱齿轮设计：工程师使用此理论计算齿轮齿面在啮合时的接触压力，优化齿形和热处理工艺，确保在数十万公里的生命周期内不发生点蚀失效，避免换挡顿挫和异响。
风力发电机主轴轴承：主轴轴承承受巨大的交变载荷。通过精确计算接触应力与L10寿命，可以选择合适的轴承型号和游隙，保障在恶劣工况下20年以上的可靠运行，减少昂贵的停机维护。
高铁轮对与轨道接触分析：分析车轮踏面与钢轨的接触斑（接触区域）应力，预测滚动接触疲劳裂纹（如“蛇行磨损”）的萌生，用于优化车轮廓形和制定打磨维护周期。
医疗器械（人工关节）：在人工髋关节球头与髋臼的接触分析中，应用赫兹理论计算接触应力，用于选择材料（如陶瓷对陶瓷）和设计关节面，以最小化磨损和延长植入寿命。
常见误解与注意事项
这类计算中最初容易陷入的误区有几个。首先是“接触半宽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 $ 这类帕里斯公式。
最终可迈向将计算多维化的阶段。本工具模型为二维（平面应变状态）的圆柱接触。实际的球轴承属于三维点接触。此时接触区域变为圆形，最大剪应力的发生深度也会改变（从约0.48a变为约0.78a）。此外，处理像齿轮这样载荷位置移动的“行进接触”，需要追踪随时间变化的应力场。首先完全理解当前的简单模型，是应对这些复杂问题的最佳基石。