1. 接触刚度的重要性
接触刚度定义了接触面上力与相对位移（或穿透）之间的关系。在abaqus 中，合理设置接触刚度能够准确模拟接触界面的力学行为，避免过度穿透或不合理的接触力传递，从而提高分析结果的准确性。不合适的接触刚度可能导致计算不收敛、结果偏差较大等问题。
2. abaqus 中的接触类型与刚度设置方式
abaqus 提供了多种接触类型，如面面接触（surface - to - surface contact）和点面接触（node - to - surface contact），不同接触类型的刚度设置方法略有差异。
l面面接触刚度设置：
（1）通用接触（general contact）：在 abaqus/explicit 中，通用接触可通过 “scale factor (general contact)” 来缩放默认接触刚度，以指定分段线性的接触压力 - 穿透关系。用户需在 “overclosure” 栏输入正值的穿透量，有 “measure”（直接输入穿透量）和 “factor”（输入最小单元尺寸的分数）两种方式。
（2）接触对（contact pair）：对于接触对定义的面面接触，在创建接触属性时设置刚度。执行 “interaction→property→create” 命令，选择 “contact” 类型。在 “edit contact property” 对话框的 “normal behavior” 选项中，选择定义接触刚度等法向接触属性。对于硬接触（hard contact），可在 “contact stiffness” 栏输入接触刚度，即接触压力 - 穿透关系曲线的斜率，该值必须为正值。在 abaqus/explicit 中，还可通过 “maximum stiffiness” 栏设置模型的最大接触刚度，默认选项为 “infinite (no slip)”，表示在运动接触中最大接触刚度为无限大，在罚函数接触中最大接触刚度为默认罚函数刚度；若选择 “specify”，则可由用户指定最大接触刚度。
l点面接触刚度设置：与面面接触类似，先定义接触对，在接触属性设置中确定点面接触的刚度。在 “normal behavior” 中设置合适的接触刚度值。点面接触常用于模拟一个表面上的离散点与另一个面之间的接触，如模拟铆钉连接等情况。
3. 接触刚度取值的确定方法
接触刚度的取值并非固定，需根据具体问题进行调试和确定。
l参考经验值：对于一些常见的材料组合和接触情况，可参考相关文献或工程经验获取初始接触刚度值。例如，在金属接触问题中，根据材料的弹性模量和接触表面的特性，可大致估算接触刚度的范围。
l试算法：先设置一个较小的接触刚度值，进行初步计算。观察计算结果中接触区域是否存在过度穿透现象。若有穿透，逐步增大接触刚度值，重新计算，直至穿透量在可接受范围内且计算结果收敛。反之，如果设置的接触刚度过大导致计算不收敛，则适当减小刚度值，并调整其他相关计算参数，如时间步长、收敛容差等。
l理论计算法：基于接触力学理论，结合材料的弹性常数、接触表面的几何形状等因素，通过公式计算接触刚度的理论值。例如，对于两个弹性半空间体的接触问题，可利用 hertz 接触理论计算接触刚度的近似值，作为 abaqus 中设置接触刚度的参考。
4. 案例演示：螺栓连接结构的接触刚度设置
以螺栓连接两块钢板的结构为例，展示在 abaqus 中设置接触刚度的具体步骤。
l模型建立：在 abaqus/cae 中创建钢板和螺栓的几何模型，并进行合理的网格划分。对于螺栓连接部分，网格应足够细密，以准确捕捉接触应力分布。例如，对于钢板采用六面体网格，网格尺寸设置为 5mm；对于螺栓和螺母，由于结构相对复杂，采用四面体网格，网格尺寸设置为 2mm。
l接触对定义：将螺母底面定义为主面，钢板的上表面定义为从面，创建面面接触对。在 “interaction” 模块中，选择 “create interaction”，选择 “surface - to - surface contact” 类型。
l接触属性设置：执行 “interaction→property→create”，创建接触属性。在 “edit contact property” 对话框中，“normal behavior” 选择硬接触（hard contact）。根据经验或初步估算，在 “contact stiffness” 栏先输入一个较小的值，如 1000mpa。同时，根据实际连接材料的特性，设置摩擦系数，假设为钢 - 钢连接，摩擦系数设置为 0.15。
l 加载与约束：在一块钢板的一端施加固定约束，限制其所有自由度。在螺栓的头部施加预紧力，如 5000n，在另一块钢板的一端施加横向力，如 1000n。
l计算与结果分析：提交计算任务，观察计算结果。若接触面上出现穿透现象，返回接触属性设置，增大接触刚度值，如调整为 2000mpa，重新计算，直至接触穿透量满足工程要求且计算收敛。分析计算结果中螺栓连接结构的接触应力分布、螺栓的受力情况以及两块钢板之间的相对位移等，评估接触刚度设置对结果的影响。
5. 注意事项与常见问题解决
l主从面选择：在定义接触对时，主面应选择刚度较大或网格较粗的面，从面选择相对较软或网格较细的面，以提高计算的收敛性和准确性。例如，在螺栓连接中，螺栓头或螺母的底面通常作为主面，被连接件的接触表面作为从面。
l网格质量：接触区域的网格质量对接触分析结果影响较大。应确保接触表面的网格划分精细且均匀，避免主从面网格尺寸差异过大。不均匀的网格可能导致接触刚度在局部区域的计算不准确，进而影响整体分析结果。
l收敛问题：接触刚度设置过高可能导致计算不收敛。此时，可适当减小接触刚度值，并结合调整其他收敛控制参数，如迭代次数、收敛容差等。同时，检查模型的约束设置是否合理，避免出现过约束或欠约束情况。若计算过程中出现接触状态频繁变化导致的收敛困难，可尝试调整接触算法，如从罚函数法切换为增强拉格朗日法（仅适用于 abaqus/standard）。
l穿透问题：接触刚度设置过低会导致接触面上出现穿透现象。逐步增大接触刚度值，同时观察计算结果，直到穿透量在可接受范围内。在某些情况下，即使接触刚度设置合理，仍可能出现少量穿透，这可能是由于模型的复杂性或数值计算误差引起的。此时，可通过调整 “adjust” 参数（abaqus 提供的自动调整接触刚度机制）来进一步减少穿透，同时确保不引入新的数值问题。
