什么是滚动接触应力
🙋
滚动接触应力是什么？听起来好复杂，是轴承里用的那个吗？
🎓
简单来说，就是两个滚动的圆球或者圆柱体，它们接触时内部产生的“挤压力”。比如你用手捏一个橡皮球，表面会凹下去，内部也会被挤压。在工程中，轴承里的滚珠和轨道之间、齿轮的牙齿咬合时，都会产生这种应力。试着在模拟器里把“载荷 P”调大，你会看到接触区域的压力瞬间飙升，这就是为什么重载下零件容易损坏的原因。
🙋
诶，真的吗？那为什么工程师总说零件是从“里面”开始坏的？表面不是压力最大吗？
🎓
这是个好问题！这正是赫兹接触理论最反直觉的地方。表面确实承受最大的压应力 $p_0$，但材料失效（比如疲劳裂纹）往往始于表面下方。你可以试着改变模拟器里两个物体的半径 $R_1$ 和 $R_2$，观察“最大剪应力”的位置。在实际工程中，比如汽车变速箱的齿轮，裂纹就是从表面下约0.47倍接触半径的深度萌生，然后扩展到表面，最终导致点蚀或剥落。
🙋
哦！所以设计时不能只看表面？那这个“等效弹性模量 $E^*$”又是干嘛的？材料参数 $E_1$ 和 $ν_1$ 我都调了，感觉结果变化好大。
🎓
没错，要“透视”到材料内部！$E^*$ 是个超级有用的合成参数，它把两种不同材料的“软硬”程度（弹性模量 $E$）和横向变形能力（泊松比 $ν$）打包在一起了。你刚才调整 $E_1$ 和 $ν_1$ 时，接触半径 $a$ 的变化，本质上就是 $E^*$ 在起作用。工程现场常见的是钢对钢的轴承，但如果是陶瓷球轴承（$E$ 很大）或橡胶轮（$E$ 很小），接触面积和压力分布会天差地别，这个模拟器能让你直观地看到这种差异。
物理模型与关键公式
赫兹接触理论的核心是计算两个弹性曲面在法向载荷下的接触半径和最大压力。它将问题简化为一个等效的弹性球体与刚性平面的接触。
$$a = \left(\dfrac{3PR^*}{4E^*}\right)^{1/3}, \quad p_0 = \dfrac{3P}{2\pi a^2}$$
$a$ ：接触圆半径。 $p_0$ ：接触区中心的最大压力。 $P$ ：法向载荷。 $R^*$ ：等效曲率半径，$1/R^*= 1/R_1 + 1/R_2$。 $E^*$ ：等效弹性模量，$1/E^*= (1-ν_1^2)/E_1 + (1-ν_2^2)/E_2$。
为了预测疲劳寿命，我们需要找到材料内部的最大剪应力及其位置，这通常是滚动接触疲劳裂纹的起源点。
$$\tau_{max}\approx 0.31\,p_0 \quad \text{at}\quad z \approx 0.47\,a$$
$\tau_{max}$ ：次表面最大剪应力。 $z$：从接触表面向下的深度。公式表明，最危险的应力点不在表面，而在表面之下约0.47倍接触半径的深度。
现实世界中的应用
滚动轴承设计：这是赫兹理论最经典的应用。工程师通过计算滚珠与内外圈滚道之间的接触应力，来选择合适的材料、热处理工艺，并预测轴承的疲劳寿命（如L10寿命），确保机床主轴或风力发电机齿轮箱能可靠运行数万小时。
齿轮传动系统：在汽车变速箱或工业减速机中，齿轮齿面在啮合时承受循环接触应力。运用赫兹理论可以优化齿形（如修形），使接触压力分布更均匀，防止早期点蚀和胶合失效，降低噪音和振动。
轮轨关系研究：高铁或地铁的车轮与钢轨的接触是一个典型的赫兹接触问题。分析接触斑（接触区域）的大小和应力，对于研究钢轨的磨损、滚动接触疲劳（如轨头龟裂）以及列车运行稳定性至关重要。
人工关节（如髋关节）仿真：在生物医学工程中，金属或陶瓷股骨头与聚乙烯髋臼杯之间的接触也适用赫兹理论。通过模拟接触应力，可以评估假体的磨损速率，优化设计以延长其使用寿命，减少患者二次手术的风险。
常见误解与注意事项
开始使用此工具时，有几个容易陷入的误区需要注意。首先，切勿简单认为“接触半径大就安全”。虽然接触面积确实会扩大，但最大接触压力p₀会随载荷P的增加，按接触半径a的平方反比急剧上升。例如，当载荷增至8倍时，接触半径变为2倍，但最大接触压力会升至4倍。即使面积扩大为4倍，单位面积承受的负荷反而增加。因此，载荷增加的影响比想象中更为严重。
其次，关于材料参数的输入。尤其容易轻视泊松比ν，但它会直接影响等效弹性模量E*。例如，钢材（ν=0.3）与橡胶（ν≈0.5）的$(1-\nu^2)$项值差异显著。实际工程中常见的情况是，材料制造商提供的目录值与实际材质存在细微差别。特别是树脂和复合材料可能存在批次差异，需要采取保守策略，例如采用稍大的弹性模量值。
最后，切勿忘记本计算基于完全弹性体且表面光滑的理想条件。实际零件会受到表面粗糙度和润滑油的显著影响，这些因素会改变应力分布。计算结果应视为“参考值”，尤其在评估疲劳寿命时，必须将此处求得的τ_max乘以安全系数，或通过实际设备测试进行验证。