赫兹接触力学计算器
参考资料： 接触应力计算公式
赫兹接触力学计算器
※ 赫兹接触压力分布为半椭球形（球接触）或半椭圆柱形（圆柱接触）。最大压力 p₀ 出现在接触中心。
理论与主要公式
等效弹性模量：
\[\frac{1}{E^*}= \frac{1-\nu_1^2}{E_1}+ \frac{1-\nu_2^2}{E_2}\]
球接触：
\[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}\]
什么是赫兹接触力学
物理模型与关键公式
首先，两个不同材料的弹性性质被综合成一个等效弹性模量，这是赫兹理论的核心简化：
\[\frac{1}{E^*}= \frac{1-\nu_1^2}{E_1}+ \frac{1-\nu_2^2}{E_2}\]
\(E_1, E_2\) 是两物体的杨氏模量，\(\nu_1, \nu_2\) 是泊松比。\(E^*\) 越小，表示整体接触刚度越软。
对于两个球体（或球与平面）的接触，接触区域是一个半径为 a 的圆，其大小和最大接触压力 \(p_0\) 由以下公式给出：
\[a = \left(\frac{3FR^{*}}{4E^{*}}\right)^{1/3},\quad p_0 = \frac{3F}{2\pi a^2}\]
其中，\(F\) 是法向载荷，\(R^*\) 是等效曲率半径（\(1/R^* = 1/R_1 + 1/R_2\)）。a 增大时，压力 \(p_0\) 会以平方反比的关系迅速减小。
现实世界中的应用
滚动轴承设计：滚珠与内外圈滚道的接触是典型的赫兹接触。工程师使用此理论计算接触应力，确保在数百万次的循环载荷下材料不发生疲劳剥落（点蚀），从而确定轴承的额定寿命。
齿轮传动：齿轮齿面在啮合时是线接触（近似圆柱接触）。赫兹接触压力是评估齿面接触疲劳强度（如点蚀）的关键指标，直接影响齿轮的承载能力和设计尺寸。
轮轨关系：火车车轮与钢轨的接触分析至关重要。通过赫兹理论计算接触斑（接触区域）的大小和应力，可以研究钢轨的磨损、塑性变形以及滚动接触疲劳裂纹的萌生。
密封与连接件：许多机械密封（如 O 型圈）和过盈配合的螺栓连接，其接触压力分布都基于赫兹接触原理进行初步估算，以确保密封有效或连接可靠。
常见误解与注意事项
首先，你是否认为"接触半径与载荷成正比"？实际上，在球体接触中，关系式为 \(a \propto F^{1/3}\)，即使载荷加倍，接触半径也仅增加约 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 分析进行验证，这是铁律。
使用指南
- 输入两个接触体的半径 R₁、R₂（单位 mm），球体输入正值，平面输入极大值
- 分别设定材料 1 和材料 2 的弹性模量 E₁、E₂（单位 GPa）及泊松比 ν₁、ν₂
- 输入法向载荷 F（单位 N）和圆柱接触长度 L（单位 mm，球体接触可忽略）
- 计算器自动求解等效曲率半径、等效弹性模量、接触半径 a、最大 Hertz 压力 pmax 及表面沉降量 δ
- 查看压力分布曲线 p(r) 和次表面最大剪应力深度 zc ≈ 0.78a 处的应力
具体计算示例
轴承钢球（R₁=10 mm，E₁=210 GPa，ν₁=0.3）与淬硬钢套（R₂=50 mm，E₂=210 GPa，ν₂=0.3）接触，载荷 F=5000 N。计算得等效半径 Req≈8.33 mm，等效模量 E*≈115 GPa，接触半径 a≈0.65 mm，最大 Hertz 压力 pmax≈5.70 GPa，表面沉降 δ≈0.050 mm。次表面最大剪应力 τmax≈1.77 GPa 出现在深度 z≈0.48a≈0.31 mm 处，该应力水平远超轴承钢的许用接触应力，可能导致疲劳剥落。
实务注意事项
- 接触应力高度集中在小面积内，pmax 与 F2/3 成正比；轴承失效常因次表面疲劳而非表面剥落，需监控 z=0.78a 深度应力
- 泊松比误差会影响沉降量计算精度 ±5%；精密机械应采用实测材料参数而非手册标准值
- 圆柱接触（如齿轮啮合）的压力分布为椭圆形，短轴半径 b 远小于长轴 a，需确保接触长度 L 充分覆盖实际接触带宽度
- 温度升高时弹性模量下降，高速轴承应考虑热软化导致接触应力上升 15-25%
- 接触压力超过材料屈服强度的 80% 时，Hertz 理论失效，需引入弹塑性修正或 FEM 非线性分析
参考资料
您可以在数据资料：接触应力计算公式 找到详细的赫兹接触应力计算方法。理论参考 Johnson《Contact Mechanics》(1985)。
