赫兹接触应力计算工具

赫兹接触理论的核心是计算接触区域的半宽（或半径）a和中心最大压力p₀。这首先需要将两个接触物体等效为一个弹性半空间与一个刚性曲面的接触，关键参数是等效弹性模量E*和等效曲率半径R*。

$$ E^* = \left(\frac{1-\nu_1^2}{E_1}+ \frac{1-\nu_2^2}{E_2}\right)^{-1}, \quad R^* = \left(\frac{1}{R_1}+ \frac{1}{R_2}\right)^{-1}$$

E₁, E₂为两物体的杨氏模量；ν₁, ν₂为泊松比；R₁, R₂为曲率半径（凸面为正，凹面为负）。E*综合了材料的软硬，R*综合了几何的弯曲程度。

对于球与球（或球与平面）的接触，接触半径a和最大接触压力p₀由以下公式给出：

$$ a = \left( \frac{3PR^{*}}{4E^{*}}\right)^{1/3}, \quad p_0 = \frac{3P}{2\pi a^2}= \left( \frac{6PE^{*2}}{\pi^3 R^{*2}}\right)^{1/3} $$

P为法向载荷。接触面上的压力呈半椭球分布：p(r) = p₀√(1 - (r/a)²)。而真正危险的是次表面的最大剪切应力τ_max，对于球接触，它约为0.31p₀，位于表面以下约0.48a深处。

滚动轴承：滚珠或滚子与内外圈滚道的接触是典型的赫兹接触。工程师使用此理论计算接触应力，以选择材料、热处理工艺和润滑方案，防止过早出现接触疲劳剥落（点蚀）。

齿轮传动：齿轮齿面在啮合过程中是线接触（近似圆柱接触）。赫兹应力是齿面强度（如接触疲劳强度）计算的基础，直接关系到齿轮的承载能力和寿命预测。

轮轨系统：火车车轮与钢轨的接触是赫兹接触理论的重要应用场景。通过计算接触斑（接触区域）的大小和应力，可以分析钢轨的磨损、塑性变形和滚动接触疲劳裂纹的萌生。

人工关节（如髋关节）：金属或陶瓷股骨头与髋臼衬垫之间的接触也适用赫兹理论。分析接触应力有助于优化假体设计，减少磨损颗粒的产生，延长人工关节的使用寿命。

开始使用此模拟器时，存在一些CAE初学者容易陷入的误区。首先一个常见的误解是认为“只要增强材料（采用高强度钢），就能承受任何高载荷”。材料强度固然重要，但赫兹接触应力本质上是“面压力”。例如，当载荷增加至2倍时，接触半径仅按 $a \propto P^{1/3}$ 关系增至约1.26倍，而最大接触压力 $p_0$ 却会按 $P^{2/3}$ 比例跃升至约1.59倍。单纯提高材料强度往往难以应对这种情况，通常需要通过增大形状（曲率半径）来增加接触面积，进行根本性的设计调整。

其次是参数设置的注意事项。模拟器中将“平面”视为曲率半径∞，但实际工程中几乎不存在完全理想的平面。例如在球轴承的滚珠与滚道接触中，滚道沟槽的曲率半径 $R_2$ 应设置为略大于滚珠半径 $R_1$ 的值（如1.05倍）。若此处简化为“平面”，接触压力的计算结果会远高于实际情况（即过于保守），可能导致过度设计。

最后需要理解本工具的基本前提。赫兹理论建立在完全弹性体、光滑表面、弹性极限内三大基础之上。这意味着当载荷引起塑性变形、表面存在粗糙度、或主要关注疲劳寿命时，计算结果仅能作为“初步近似”。例如淬火钢表面存在微小凹痕时，该处会产生应力集中，破坏可能从完全不同于理论最大剪切应力位置处开始。模拟结果应始终视为理想化的“基准值”，必须通过实物测试或更高级的CAE分析进行验证。