1 概念与研究对象
接触力学是研究两物体在相互接触时、随载荷与边界条件变化的力—位移关系及其内部力学响应的学科。其核心目标在于描述接触界面处的压力分布与接触面积如何随法向与切向载荷演化，并进一步推导由此产生的应力、应变与应变能等结果。由于实际工程材料往往存在弹性与非弹性特征、表面存在微观粗糙度、界面还伴随摩擦与粘着行为，接触问题通常带有明显的非线性。
1.1 接触问题的基本描述
1.1.1 接触几何与接触状态（分离/贴合/滑移）
“接触几何”不仅指两接触体在空间中的相对位置，也包括接触区边界的形状与面积。随着载荷变化，接触状态可能发生变化，例如：在法向载荷不足时出现分离；在载荷增大后形成贴合接触；在切向载荷超过摩擦能力时出现滑移。对工程分析而言，区分这些状态并建立其判别条件，是求解接触力学的前提。
1.1.2 载荷—位移关系的关键量
接触问题常用若干关键量刻画输入与响应关系。法向方向上，通常以法向位移（或压入量）与法向载荷之间的关系描述“压多少承受多少”；切向方向上，关注切向力与相对滑移或微滑（partial slip）的对应关系。除力与位移外，接触面积、接触压力峰值位置、平均压力、以及接触区边缘的奇异性等量，往往决定疲劳、磨损与密封等后续失效评估的结论。
1.2 理论框架的核心假设
1.2.1 连续介质与材料本构
接触力学通常采用连续介质理论，将材料视作在空间尺度上具有连续力学性质。材料的本构关系（如弹性或粘弹性）决定应力—应变如何关联，并进而影响接触压力分布与接触刚度。对于工程常见场景，常将材料简化为各向同性或横向各向同性等形式，以获得可解析或便于数值求解的模型。
2.2.2 弹性、塑性与粘弹性建模思路
在弹性范围内，接触问题可较容易建立力—位移关系与应力场。若载荷足够使局部达到屈服，则需引入塑性，处理残余变形与接触历史的影响。若材料表现出显著应力松弛或滞回特征，还可采用粘弹性框架，将时间或频率尺度引入本构，从而描述速率效应与能量耗散。工程上往往需在“模型复杂度”与“预测精度”之间做权衡。
1.3 接触力学的工程意义
1.3.1 可靠性与寿命预测
接触区域的应力与剪应力循环是疲劳失效的重要来源。接触力学提供从载荷与几何到应力场的桥梁，使得裂纹萌生位置、循环应力幅值及其演化趋势可以被估计。对可靠性设计而言，这意味着能够在不完全依赖试验的情况下进行参数敏感性分析与寿命筛选。
1.3.2 摩擦磨损与密封性能
摩擦决定切向力与发热、微滑状态以及磨损速率；接触压力分布影响密封界面是否持续处于有效贴合。接触力学因此常与表面与界面过程模型耦合，用以解释“为什么压力提高能改善密封但可能加速磨损”等具有工程权衡性质的现象。
1.3.3 轮—路与齿—轮等典型应用
轮胎与路面之间、齿轮啮合面之间均存在显著接触变形与摩擦效应。接触力学的经典模型与扩展模型可用于理解变形—载荷耦合、接触椭圆特征、以及由界面应力导致的疲劳与磨耗。对传动系统，合理的接触分析可帮助优化齿面修形与润滑策略。
2 接触条件与边界要求
接触问题之所以难以直接求解，关键在于接触区域未知且需满足“无穿透”“力平衡”“摩擦边界条件”等约束。良好的边界描述会显著影响数值算法的稳定性与收敛性，也影响物理上接触状态判别的合理性。
2.1 接触约束方程
2.1.1 法向间隙与无穿透条件
法向方向上，两接触面在几何上可以用间隙量表达。无穿透条件要求接触体不能相互穿过：当间隙为正时应无法向接触力；当间隙逼近零时可产生接触压力。由此形成一种互补关系（接触状态与压力共同满足），使系统呈现典型的非线性约束特征。
2.1.2 接触压力的有效性条件
除无穿透外，接触压力还必须遵循物理有效性，例如通常要求接触压力非负（只允许“顶住”而不允许“拉住”）。当切换为分离状态时，压力应降为零。该条件与位移或间隙的变化共同决定接触域边界的移动方式。
2.2 切向接触与摩擦模型
2.2.1 库仑摩擦与静/动摩擦
库仑摩擦是工程中最常见的摩擦描述之一：切向接触力受法向压力调控，并与摩擦系数相关。通常需要引入静摩擦与动摩擦概念，用以描述微滑与完全滑移之间的转变：当切向载荷未超过摩擦阈值时可处于贴合的“静状态”；超过阈值则进入滑动。
2.2.2 罚函数与粘着/滑移判别
在数值求解中常采用罚函数等思路将接触约束“软化”，通过引入惩罚刚度来逼近互补条件。与此同时，粘着—滑移判别需要确定接触界面的相对运动状态：是保持粘着（或等效静摩擦）还是发生滑移。具体判别往往依赖局部法向压力与切向牵引之间的比较。
2.3 粗糙度与真实接触
2.3.1 名义面积与真实接触面积
实际工程表面通常具有微观起伏，导致名义接触面积与真实承压面积差异显著。真实接触面积一般远小于名义面积，它决定了压力集中程度与界面强度。粗糙度还会引入接触点的统计分布，从而影响整体摩擦系数与磨损趋势。
2.3.2 表面形貌对压力分布的影响
表面形貌会改变局部峰谷几何，使得压力从“相对平滑”的分布转变为具有离散承载点或小区域的分布。尤其在高压力或较软材料情况下，粗糙度导致的局部变形与材料流动会显著改变接触刚度与能量耗散路径。因而，在需要高精度时，需引入粗糙度统计模型或直接进行表面级建模。
3 经典接触模型
经典模型的价值在于提供可解析或半解析的关系，用于理解关键尺度规律，并为复杂数值计算提供校验基准。模型通常基于简化几何与材料假设，例如理想化的半空间、线弹性与光滑界面。
3.1 刚体—弹性体接触
3.1.1 点接触：Hertz 理论（基本推导思路）
Hertz 理论描述光滑刚体与弹性体之间的点接触。其基本思路包括：将接触区域内的法向位移与压力分布建立联系；利用弹性半空间的解与对称性推导接触区尺寸；并得到载荷—压入量的幂律关系。该理论还能给出接触压力的分布特征，例如在中心附近达到峰值并向外衰减。
3.1.2 线接触：相关弹性解的适用场景
线接触对应接触长度远大于接触半径的情形，常用于描述滚动体在局部区域的接触。与点接触相比，压力分布在横向方向表现出不同的衰减规律，同时载荷与接触宽度之间的关系也随之改变。该类模型常与工程几何参数（曲率）相联系，便于估算滚动接触应力水平。
3.2 线性弹性与半空间假设
3.2.1 弹性常数与等效模量
经典接触模型通常需要将两材料的弹性参数合并为等效模量。等效模量体现了两接触体共同的“柔度”，从而决定接触刚度与应力响应。其计算依赖材料的弹性模量与泊松比，能够把“双体问题”简化为“单等效弹性体”处理。
3.2.2 接触刚度与等效曲率
接触几何以曲率参数表征。等效曲率将两接触体在法向截面上的弯曲程度合并，进而影响接触区尺寸与压力峰值。接触刚度则可以理解为法向载荷随压入量变化的斜率，它与等效模量、等效曲率以及接触形状直接相关。
3.3 接触面积与应力特征
3.3.1 最大接触压力位置
在典型光滑弹性接触中，最大接触压力通常出现在接触区内部的特定位置，常与对称性一致，且随几何与载荷变化而改变其量级。最大值对于疲劳与材料失效评估尤为关键，因为它决定局部是否可能超过材料屈服或疲劳阈值。
3.3.2 应力场随距离的衰减规律
接触引起的应力在空间中并非均匀分布。经典弹性解表明，离开接触表面或远离接触区边缘后，应力会逐渐衰减，且不同应力分量可能具有不同衰减率。理解这种衰减规律有助于确定次表面风险点的位置，并用于寿命预测模型的参数化。
4 多因素扩展与改进模型
现实工程往往超出理想光滑弹性接触范畴。扩展模型的重点在于引入非线性本构、速率相关效应、粗糙度统计与多层结构界面等因素，以更贴近实际的接触演化过程。
4.1 材料非线性与塑性接触
4.1.1 弹塑性接触的判据与求解思路
弹塑性接触需要判断局部是否达到屈服并进入塑性区。常见思路是：在法向载荷逐渐增大的过程中，先以弹性解估计应力水平；当应力达到屈服条件后，将塑性区引入并更新接触压力分布与接触面积。与纯弹性相比，载荷—位移关系通常呈现更强的非线性。
4.1.2 残余变形与接触历史效应
塑性变形不会完全可逆。加载—卸载循环后，接触体会保留残余凹陷或塑性流动，导致后续接触的初始几何状态改变。于是接触结果不仅与当前载荷相关，还与历史加载路径有关，这对评估疲劳累计损伤与磨损进程尤为重要。
4.2 粘弹性接触与速率效应
4.2.1 频率/时间尺度的影响
粘弹性材料具有与时间或频率相关的刚度与阻尼特性。相同法向载荷下，不同加载速率可能对应不同的等效模量，从而改变接触面积、压力水平以及界面温升趋势。工程上常需要根据工作频率或循环周期选取合适的粘弹参数。
4.2.2 能量耗散与滞回特性
由于粘弹性材料的滞回响应，接触过程会产生能量耗散。该耗散与界面黏性剪切、材料内部损失相关，可能进一步影响润滑性能、热-力耦合与磨损发展。因而，粘弹性模型不仅用于几何与应力预测，也用于评估热生与寿命衰减。
4.3 表面形貌与粗糙接触统计模型
4.3.1 统计分布与等效压力
粗糙表面可用统计方法描述其高度分布与空间相关性。将局部峰谷视为承压单元后，可以得到真实接触面积随载荷增长的规律，并将离散压力合成为等效的平均压力与压力分散度。这类模型通常适用于不追求点级细节、而强调整体趋势的工程分析。
4.3.2 功率谱/高度相关的建模思路
表面形貌的空间尺度可通过功率谱等量表征，从而反映不同波长的起伏对接触的贡献。通过引入高度相关函数，可以估计在某一特征尺度上的接触高度超越概率。该过程使得模型能够在多尺度粗糙表面中给出更合理的压力分布统计特性。
4.4 复合材料与层状结构接触
4.4.1 涂层—基体的等效响应
涂层覆盖在基体上时，接触变形受涂层厚度、弹性模量与粘结层性能影响。涂层较薄时可能更接近基体主导的响应；当涂层厚度或模量比发生变化时，接触区的应力集中可能显著改变。等效响应方法常用于把层状问题简化为具有等效材料参数的连续体。
4.4.2 界面强度与脱粘风险
界面强度决定涂层或层状材料是否发生剥离。接触载荷会引入法向与切向应力分量，界面处的拉应力、剪应力以及应力梯度共同影响脱粘风险。对于涂层耐磨或防护设计，接触力学模型用于评估“承载能力”与“界面可靠性”的平衡。
5 接触疲劳、磨损与失效相关
接触失效往往不是单一机理导致，而是疲劳、磨损与环境耦合作用的结果。接触力学提供应力场与界面状态的基础输入，使得机理模型能够被参数化并用于工程预测。
5.1 接触疲劳机理
5.1.1 裂纹萌生与应力循环
在循环载荷下，接触区的反复加载会引起材料内部的应力循环，最终导致裂纹萌生。裂纹往往在应力幅值较大且缺陷更敏感的位置产生。接触力学通过给出次表面应力分布与循环应力幅值，帮助判断萌生的可能位置与方向性趋势。
5.1.2 次表面最大剪应力观点
很多接触疲劳分析会强调次表面的剪应力集中现象。由于接触压力在表面附近形成特定的应力结构，剪应力最大值可能出现在表层下方，从而解释了“表面看似完好但内部已发生损伤”的现象。以最大剪应力为代表的准则可用于初步判别疲劳风险。
5.2 磨损模型与接触演化
5.2.1 干摩擦与磨粒磨损的基本分类
磨损可按界面机理划分为不同类型，例如干摩擦引起的材料转移与磨粒造成的磨削/犁削等。不同磨损形式对压力、相对速度、润滑状态和粗糙度敏感性不同。接触力学给出压力与微滑等关键变量后，磨损模型才能进一步估计质量损失或表面形貌变化。
5.2.2 表面形貌的逐步变化（“越磨越光”或反之）
磨损会改变表面起伏与粗糙度，从而反馈到下一阶段的接触压力分布与摩擦水平。某些条件下表面峰值逐渐被削弱，可能出现“越磨越光”；在另一侧，磨粒嵌入或材料剥落可能导致粗糙度反向增长，表现为“越磨越糟”。因此需把接触演化与表面演化作为耦合过程来理解。
5.3 密封与泄漏风险
5.3.1 压紧力与接触状态敏感性
密封效果与压紧力密切相关。过小的压紧力可能导致局部未能形成有效贴合，出现通道并引发泄漏；过大的压紧力则可能加速材料老化与磨损，反而降低长期稳定性。接触力学用来分析压紧力变化如何影响接触面积与压力保持能力，从而评估泄漏敏感性。
5.3.2 失效形式：泄漏、磨穿与老化
密封失效常见形式包括泄漏通道形成、材料磨穿以及随时间累积的老化或性能衰退。压力分布、剪切作用与摩擦发热会共同促进界面退化。通过把接触分析嵌入寿命框架，可以更系统地预测何时进入不可靠区间。
6 工程数值方法与建模流程
当几何复杂、材料非线性或接触条件难以满足解析解假设时，工程上常依赖数值方法（尤其是有限元）进行求解。接触算法的选择直接影响效率与精度，因此需要建立清晰的建模流程与验证策略。
6.1 有限元中的接触算法概览
6.1.1 刚度法/迭代求解与收敛问题
有限元求解接触问题通常需要迭代，因为接触状态会随位移更新。刚度法等迭代流程在每个步长内更新接触压力与接触域，并通过平衡方程迭代到收敛。收敛困难可能源于接触状态频繁切换、网格质量不足或罚参数设置不合理，因此需要经验性调参与策略选择。
6.1.2 接触面离散与网格敏感性
接触面附近的网格分辨率会影响压力峰值与接触边界位置的预测。网格过粗可能导致压力被“平均化”而低估峰值；网格过细则可能提高计算成本并产生数值噪声。良好的网格策略通常包括对接触区进行局部加密并进行网格无关性验证。
6.2 约束处理方式
6.2.1 拉格朗日乘子法与罚函数法
拉格朗日乘子法通过引入额外未知量严格满足约束，通常具有较好的约束精度；罚函数法则通过惩罚项近似无穿透条件，结构实现更简便，但需要合理选择罚参数以平衡精度与数值稳定性。两者各有适用场景：前者可能更严格，后者可能更易部署。
6.2.2 无穿透与接触面法向处理
数值实现中通常对法向间隙进行监控，并在计算中施加无穿透处理，避免出现非物理穿透。接触面法向定义、主从面选择、以及接触法向一致性都会影响结果稳定性。对于强非线性接触，必须确保这些实现细节与几何尺度匹配。
6.3 后处理与验证
6.3.1 压力分布与接触面积提取
求解得到的接触压力分布和接触面积是后处理重点。接触压力云图可用于检查压力峰值位置是否符合理论或经验预期；接触面积可用于评估接触刚度变化与密封有效性等指标。若需要与实验对比，通常还需将数值结果转化为与测量方式一致的量。
6.3.2 网格收敛与与解析解对比
验证通常包含网格收敛性检查与对解析/半解析解的对比。对于满足经典假设的简化模型，可用解析解建立基准曲线，例如载荷—位移关系与压力分布形式。通过对比可以定位误差来源，如网格、材料参数、摩擦与约束实现等。
6.4 多尺度与耦合建模思路
6.4.1 粗糙度尺度与结构尺度的衔接
多尺度方法用于连接表面粗糙细节与结构整体响应。常见策略是：在结构尺度上进行有限元计算，同时用统计粗糙度模型提供等效的接触参数或边界条件；或在局部区域进行更细致的表面级分析以校准等效参数。这样可以在不计算全部微观细节的前提下，提高预测的物理一致性。
6.4.2 热-力-接触耦合的简述
摩擦接触会产生热量，热量进一步影响材料性能与接触状态，从而形成热-力-接触耦合问题。工程上通常通过在接触区计算热源，再将温度场反馈到材料本构与接触边界条件中。该方向对高速、重载或强摩擦场景尤其重要。
7 典型应用场景（工程化视角）
接触力学在工程中的意义体现在：把复杂界面现象转化为可计算、可优化的指标，从而服务于设计与失效分析。下面列举若干高频应用场景，展示经典模型与扩展模型的落地思路。
7.1 轴承与滚动接触
7.1.1 载荷传递与接触椭圆特征
轴承中滚动体与滚道之间的接触通常表现为沿某方向延展的接触区。载荷由接触压力分布承载，接触区往往具有椭圆特征（由几何曲率和相对运动共同决定）。根据接触力学的关系，可以估算最大压力与接触尺寸，为材料疲劳与寿命评估提供输入。
7.1.2 微滑与效率/磨损影响
即使在滚动体“名义纯滚动”条件下，也可能因为弹性变形、载荷变化与润滑状态等导致微滑。微滑与摩擦共同决定能量损失与磨损发展。接触模型给出微滑范围与压力水平，从而帮助评估效率与寿命之间的权衡。
7.2 齿轮接触
7.2.1 点/线近似与啮合运动
齿轮啮合可在局部近似为点或线接触，具体取决于几何与运动状态。通过把啮合瞬间映射到等效接触问题，可以估算接触应力峰值以及随啮合位置变化的载荷分布。该分析对于齿面修形与材料选择具有直接指导意义。
7.2.2 胶合与疲劳的工程考虑
在润滑不足或滑移较大时，齿面可能发生胶合风险；在循环应力作用下则可能出现疲劳裂纹。接触力学提供应力场、温升相关的输入变量与接触压力分布，从而与润滑与材料疲劳模型联用进行风险评估。
7.3 轮胎与路面接触
7.3.1 变形与力学耦合的直观理解
轮胎接触路面时，轮胎材料的变形决定法向反力与接触面积分布。接触力学帮助解释：为什么外观上接触区大小的变化并不总与载荷简单成比例，而与轮胎结构与材料本构有关。变形也影响滚阻与操控指标。
7.3.2 抓地与磨耗的相关性
抓地能力与界面摩擦以及有效接触状态相关。与此同时，压力与滑移状态决定磨耗速率与表面形貌演化。通过将摩擦与接触压力耦合，可以从机理上解释磨耗“为什么会在某些路况与工况下更快”。
7.4 密封与法兰连接
7.4.1 压紧—变形—泄漏的链路
法兰连接中，螺栓施加的压紧力导致密封件产生接触变形与压力建立。接触状态（贴合与否、压力是否维持）共同决定是否形成泄漏通道。接触力学将压紧力转化为接触压力分布与接触面积，从而支持密封设计与装配参数选择。
7.4.2 材料选型与表面加工的作用
材料的弹性、粘弹性与界面强度决定密封件在载荷循环下的承压与保持能力；表面加工（如粗糙度与形貌）影响真实接触面积与压力集中程度。结合接触分析可以对材料等级、硬度、涂层方案以及加工精度进行取舍。
7.5 轻度“梗”式理解：接触即妥协
7.5.1 “接触面积不是越大越好”的工程直觉
在直觉层面，工程师常发现：增加压紧或追求更大的名义接触面积，未必带来更好的综合性能。因为更大的接触面积可能伴随更高的摩擦功与磨损风险，或引起更强的材料退化。接触力学提供了从压力分布与摩擦边界到“代价”的量化视角，使这种直觉更可验证。
7.5.2 摩擦与磨损：免费的却不便宜
很多接触系统中摩擦看似“不可避免”，但它的能量耗散与材料损伤成本并不便宜。尤其当载荷、粗糙度或润滑状态变化时，磨损可能加速并导致性能劣化。接触力学与磨损模型的组合，使这种“免费的却不便宜”的经验感受能够转化为可预测指标。