赫兹接触理论（球−平面）
赫兹接触理论（球−平面）的理论基础
概要
老师，我听说Hertz接触在FEA接触分析验证中最常用，为什么呢？
两个弹性体接触时，接触半径、最大面压、接触变位的所有值都存在闭形式的解析解。1882年由Hertz推导。球与平面的接触最为基本，是非线性接触分析Code Verification的入口。
非线性的情况下也有解析解吗？
Hertz理论的巧妙之处就在这里。虽然包含接触面积随荷载变化的几何非线性，但通过弹性半空间的Boussinesq解的叠加和几何适合条件的联立，可得到闭形式。但有前提条件：接触面与平面相比足够小、在弹性范围内、无摩擦。
支配方程
请给出具体的公式。
半径为 $R$ 的弹性球以荷载 $P$ 作用在弹性平面上。等效弹性系数 $E^* = [(1-\nu_1^2)/E_1 + (1-\nu_2^2)/E_2]^{-1}$，等效半径 $R^* = R$（平面 $R_2 = \infty$）
接触半径: $a = \left(\frac{3PR^*}{4E^*}\right)^{1/3}$
最大面压: $p_0 = \frac{2aE^*}{\pi R^*} = \frac{1}{\pi}\left(\frac{6PE^{*2}}{R^{*2}}\right)^{1/3}$
接近量: $\delta = \frac{a^2}{R^*} = \left(\frac{9P^2}{16R^*E^{*2}}\right)^{1/3}$
面压分布如何？
接触面内的面压为半椭圆分布。
中心最大值为 $p_0$，接触边缘 $r = a$ 处为零。该分布可通过Boussinesq解的叠加再现，接触面下的应力场也可精确计算。最大von Mises应力出现在表面正下方 $z \approx 0.48a$ 处，为滚动疲劳的起始点。
基准数值例
请给出具体数值的验证例。
钢球（$R = 10$ mm、$E = 200$ GPa、$\nu = 0.3$）对钢平面（同材料）。$P = 100$ N。
$E^* = 200/(2 \times (1-0.09)) = 109.9$ GPa
$a = (3 \times 100 \times 0.01 / (4 \times 109.9 \times 10^9))^{1/3} = 0.0727$ mm
$p_0 = 3P/(2\pi a^2) = 9013$ MPa
$\delta = a^2/R = 0.000529$ mm
| 物理量 | 理论值 | Abaqus (C3D20R) | 误差 |
|---|---|---|---|
| 接触半径 $a$ [mm] | 0.0727 | 0.0731 | 0.55% |
| 最大面压 $p_0$ [MPa] | 9013 | 8940 | 0.81% |
| 接近量 $\delta$ [mm] | 0.000529 | 0.000527 | 0.38% |
精度在1%以内呢。网格细度如何？
接触面内布置了20个以上的单元。由于接触半径 $a$ 很小，接触部分的单元尺寸需要 $a/10 \approx 7$ μm。远离接触区域可以粗化，因此偏置网格是必需的。
赫兹接触理论（球−平面）的数值计算方法
接触算法选择
FEA中的接触算法有哪些种类？
主要有三种。
| 方法 | 原理 | 精度 | 计算成本 |
|---|---|---|---|
| 惩罚法 | 穿透量与反力成正比 | 取决于惩罚系数 | 低 |
| 拉格朗日乘数法 | 严格满足穿透量=0 | 高精度 | 高（自由度增加） |
| 增强拉格朗日法 | 惩罚法+拉格朗日法混合 | 高精度 | 中等 |
对Hertz验证推荐使用拉格朗日乘数法或增强拉格朗日法。惩罚法会因惩罚系数值而改变面压，不适合验证。
各求解器的默认设置如何？
Abaqus默认采用增强拉格朗日法（*SURFACE INTERACTION + *CONTACT PAIR）。Nastran的SOL 101 GLUE/SLIDE接触采用惩罚制。Ansys默认采用增强拉格朗日法。接触问题的V&V中，使用的算法必须明确说明。
网格设计
接触问题的网格设计中最重要的是什么？
接触面内的单元尺寸要小于等于接触半径 $a$ 的1/10。粗网格会使接触面呈阶梯状，面压分布不准确。
推荐配置：
- 接触面附近：$h_{elem} = a/10$。采用HEX20结构网格
- 离接触面远的区域：用偏置逐渐粗化，比例约2.0
- 球的1/4对称模型（两个对称面）
- 平面侧模型尺寸至少为 $10a$ 以上半径
接触面应该使用哪一个面作为主面/从面？
Abaqus约定凸面（球）为从面，平面为主面。从面这一侧应该网格更细。Nastran和Ansys也遵循类似原则，将刚性更高的面设为主面。当两个面网格密度差异大时，会影响面压精度，因此接触面两侧的网格密度应保持一致。
收敛要点
接触分析经常不收敛，有什么排查办法？
对于Hertz这样光滑接触的问题，通常收敛无问题。不收敛时按以下顺序检查：
1. 接触初始状态：检查球与平面初始是否略有重叠（过度），或间距是否太大。初始间隙过大会导致接触检测失败
2. 荷载步分割：从最终荷载的1/10开始，逐步增加。Abaqus中设置Step内的初始/最大增量为较小值
3. 接触稳定化：*CONTACT CONTROLS, STABILIZE 添加微小粘性阻尼来帮助收敛。最终结果需确认稳定化无影响
4. 单元类型：HEX8在接触面易出现振荡（接触/非接触状态反复）。改用HEX20
什么是振荡？
接触边缘的节点在迭代计算中反复出现open/closed状态，这是收敛最大的障碍。采用面-面离散化（而非节点-面）可以改善。Abaqus中*CONTACT PAIR, INTERACTION=... 指定 TYPE=SURFACE TO SURFACE。
赫兹接触理论（球−平面）的实际应用
验证步骤
如何系统地进行Hertz接触的V&V验证？
分阶段进行。
1. Boussinesq问题验证（前置步骤）：确认集中荷载下半空间的应力场与理论值一致
2. 赫兹接触（同材·球-平面）：最基本设置中验证接触半径、面压、接近量
3. 赫兹接触（异材·球-球）：检验等效弹性系数和等效半径的计算
4. 含摩擦的Hertz接触：与Mindlin理论比较，验证摩擦力分布
5. 弹塑性接触：与硬度关系（Meyer定律）比对。无理论解，需要交叉验证
各阶段的评价指标是什么？
第1阶段：应力与理论值的一致程度。第2-3阶段：接触半径（误差<5%）、最大面压（误差<5%）、荷载-变位关系的 $P \propto \delta^{3/2}$ 指数一致性。第4阶段：滑移起始荷载。第5阶段：残留变形和面压分布在求解器间的一致性。
滚动疲劳应用
Hertz接触与滚动轴承设计的关系是什么？
轴承的接触应力分布用Hertz理论算出，接触面下的交变剪应力支配滚动疲劳寿命。ISO 281的动额定载荷基于Hertz理论的Lundberg-Palmgren理论计算。
用FEA求解Hertz接触后，验证最大剪应力 $\tau_{max}$ 的位置和数值，可确认轴承设计工具的合理性。理论值为 $\tau_{max} = 0.31 p_0$（$z = 0.48a$）。
椭圆接触（接近线接触）如何处理？
Hertz一般理论涵盖椭圆接触，用包含椭圆积分的公式表示。实用中可用Hamrock-Dowson近似式。滚子轴承和齿轮齿面接触属于此类。FEA用三维模型直接求解，但通过理论值比较来保证精度。
结果评估方法
接触面积的准确评估方法是什么？
FEA中用接触状态（open/closed/sliding）标志判断接触区域。Abaqus中，COPEN（接触张开量）≤0的节点处于接触状态。该节点群的外边界对应接触半径。
更准确地，对面压>0的区域积分求面积。CPRESS > 0 的单元面积合计应与 $\pi a^2$ 一致。网格粗时会因接触边界的阶梯状离散化而误差大，需用细网格。
赫兹接触理论（球−平面）的软件比较
Abaqus中的实现
Abaqus中如何设置Hertz接触？
用*CONTACT PAIR将球面设为从表面，平面设为主表面。无摩擦情况下，设置*SURFACE INTERACTION + *SURFACE BEHAVIOR, PRESSURE-OVERCLOSURE=HARD。
在NLGEOM=YES的*STATIC步中以Ramp方式加载。用自动增量（*CONTROLS, PARAMETERS=FIELD控制接触稳定性），直到接触在最终荷载处稳定。
输出：CPRESS（面压）、COPEN（间隙）、U（位移）作为场输出。反力RF和位移U作为历史输出记录，绘制 $P-\delta$ 曲线。
General Contact 和 Contact Pair 哪个更适合？
V&V推荐用Contact Pair。General Contact算法自动选择，透明度较低。Contact Pair可明确指定Surface-to-Surface离散化+增强拉格朗日法，便于在报告中说明设置。
Nastran中的实现
Nastran中的接触设置如何？
SOL 400（非线性）+ BCTABLE/BCONPRG定义接触。BCBODY1为主体，BCBODY2为从体。用BCTPARM控制接触算法参数。
Nastran的接触设置比Abaqus复杂，选项众多，但可以明确切换惩罚法和拉格朗日法。
Abaqus和Nastran结果是否一致？
在相同网格密度和接触算法下，接触半径和最大面压误差在2～3%以内。荷载-变位曲线的 $P \propto \delta^{3/2}$ 指数两者都接近理论值1.5（约1.48～1.52）。
开源软件验证
CalculiX和Code_Aster也能做接触分析吗？
CalculiX用*CONTACT PAIR实现Abaqus兼容性，支持HEX20单元的面-面接触。Hertz问题上与商用求解器不相上下。
Code_Aster用 DEFI_CONTACT + FORMULATION='DISCRETE' 或 'CONTINUE'定义接触。CONTINUE对应增强拉格朗日法，适合Hertz问题。
开源软件接触分析的可信度如何评估？
Hertz问题是最好的试金石。若开源结果与理论值误差≤5%，网格收敛可验证，则该接触算法的基本可信度有保障。再与商用求解器交叉验证，就构成双重独立验证。
赫兹接触理论（球−平面）的前沿研究
弹塑性接触
荷载很大导致塑性变形时怎样？
最大von Mises应力达到屈服应力时，塑性变形开始。屈服起始荷载为 $P_y = (\pi \sigma_y)^3 R^{*2} / (6E^{*2}) \times C_y^3$，其中 $C_y \approx 1.6$ 是理论推导值（Tabor, 1951）。
塑性进展时，面压分布从Hertz的半椭圆逐渐变为均匀分布，接触半径比弹性预测更大。此域无解析解，FEA交叉验证或实验对比成为Validation手段。
FEA求解弹塑性接触有什么注意事项？
材料模型选择影响结果很大。等向强化和移动强化会改变卸载后的残余变形。von Mises（J2塑性）是标准，但需考虑接触面下三轴应力状态下的模型合理性。网格需要比弹性情况细2～3倍，塑性区的单元尺寸≤接触半径/20。
动态接触（碰撞）
球撞击平面时怎样？
准静态应用Hertz理论可预测碰撞时间 $T$ 和最大变形 $\delta_{max}$：
$m$ 为球质量，$v$ 为碰撞速度。用显式法（显式法、LS-DYNA、Abaqus/Explicit）可验证。反发系数 $e = 1$（完全弹性碰撞）的再现是验证指标。
CFL条件如何设置？
接触部单元尺寸小时，CFL条件 $\Delta t \leq h_{min}/c_L$（$c_L$为纵波速度）限制时间步极小。质量缩放（人工增加密度）可增大 $\Delta t$，但碰撞分析中惯性很重要，应慎用。
考虑粗糙度的接触
实际表面有粗糙度，与理论偏离多大？
Greenwood-Williamson（1966）模型最常用。将表面粗糙度视为多个微小突起的Hertz接触，统计处理。粗糙度影响下，见面积比名义值小得多，实际接触面约为名义面积的1～10%。
FEA可直接建立分形表面模型，或用包含粗糙参数的表面相互作用模型（Abaqus的*SURFACE BEHAVIOR with roughness等）。计算成本剧增，多尺度方法和ROM是研究方向。
赫兹接触理论（球−平面）的故障排查
不收敛的情况
接触分析不收敛时的排查步骤是什么？
按以下顺序检查。
1. 接触初始状态：确认球与平面初始是否略有重叠（穿透），或间距过大。初始间隙太大会导致接触检测失败
2. 荷载步分割：从最终荷载的1/10开始，逐步增加。Abaqus中用Step内的初始/最大增量控制
3. 接触稳定化：*CONTACT CONTROLS, STABILIZE 添加微小粘性阻尼助力收敛。最终结果需验证稳定化无影响
4. 单元类型：HEX8在接触面易出现颤振（接触/非接触反复）。改用HEX20
什么是颤振？
接触边缘节点在迭代中反复出现open/closed状态，这是收敛最大障碍。采用面-面离散化（非节点-面）可改善。Abaqus中*CONTACT PAIR, INTERACTION=... 指定TYPE=SURFACE TO SURFACE。
面压与理论值不符的情况
面压分布不成半椭圆是怎回事？
可能原因有多个。
- 网格不足：接触面内最少需10个单元。5个以下会导致半椭圆分布呈锯齿状离散化
- 惩罚法影响：惩罚系数小时穿透量大，面压过小。改用拉格朗日法
- 球面几何精度：CAD球面被分面化时接触不均。用二次单元（HEX20）的中间节点正确放置在球面上
中间节点位置真的那么重要？
非常重要。HEX20中间节点在辺正中点（直线边）与放在球面正确位置（曲线边）时，接触面几何精度相差很大。Abaqus的*SURFACE, TYPE=ELEMENT启用curved elements，或网格生成时把中间节点精确放在曲面上，可改善。
荷载-变位关系验证
如何验证 $P \propto \delta^{3/2}$？
用History输出记录反力 $P$ 和接近量 $\delta$，绘制双对数图，检查斜率是否为1.5。最小二乘法拟合得指数 $n$。理论值$n = 1.5$，可接受范围为 $n = 1.48～1.52$。
该指数检验比单个值（$a$、$p_0$）比较更鲁棒。受网格粗度影响小，是接触算法总体合理性的指标。
相关主题
更详细
错误