2. 接触对与接触属性基础
在 Abaqus/Standard 或 Abaqus/Explicit 中建立接触,第一步不是直接给两个零件加“接触”两个字,而是要把“谁和谁可能接触”明确成一对可被求解器识别的对象,这就是接触对(contact pair)。接触对由两个面(或一组面)组成:一个被指定为主面(master surface),另一个被指定为从面(slave surface)。求解器在每一个增量步中检查从面节点是否穿透主面,并据此施加法向约束和切向摩擦约束。因此,接触对的定义方式、主从面的选择、离散化方法以及法向/切向参数,直接决定了接触计算的稳定性与精度。
一、接触对的定义方式
在 Abaqus/CAE 中,接触对通过 Interaction 模块创建,常用入口是 Create Interaction,类型选择 Surface-to-surface contact (Standard) 或 General contact (Explicit)。对于 Standard 分析,典型流程是:先分别定义两个 Surface(可以是基于几何的面,也可以是基于网格节点的面),再在接触对中指定主面和从面,最后赋予一个 接触属性(contact property)。接触属性与接触对是分离的:同一个接触属性可以被多个接触对引用,这样便于统一管理摩擦系数、法向刚度等参数。
接触属性本身由两部分组成:法向行为(Normal Behavior)和切向行为(Tangential Behavior)。法向行为决定“压不进去”的约束方式,切向行为决定“能不能滑、滑多少”的摩擦规律。两者必须同时考虑,缺一不可。
二、主从面选择原则
主从面的选择不是随意的,它直接影响接触搜索效率和结果精度。一般遵循以下原则:
- 刚度大、网格粗的一侧作主面:主面可以穿透从面,但从面节点不允许穿透主面。若把软材料或细网格面设为主面,容易产生明显穿透。
- 凸面作主面,凹面作从面:凸面与凹面接触时,凹面更容易被搜索到,作为从面可减少漏检。
- 大面作主面,小面作从面:小面作为从面时,节点数量少,搜索代价低。
- 刚体作主面:刚体不参与变形,作为主面可避免从面节点穿透刚体。
当两个面刚度、网格相近时,可以互为从面,但 Abaqus 中一个接触对只能指定一个从面。此时可考虑使用通用接触(General Contact),它自动处理面与面之间的主从关系,适合复杂装配体。
三、接触离散化方法
接触离散化决定接触约束施加在哪些点上。Abaqus 提供两种主要方法:
- 点对面(Node-to-Surface):从面节点与主面片之间建立约束。计算量小,适合主面光滑、从面网格较粗的情况,但容易产生“节点滑过主面”的误差。
- 面对面(Surface-to-Surface):在从面片和主面片之间建立积分点约束,精度更高,对网格不匹配不敏感,是 Standard 中的默认推荐方法。
在 Abaqus/Explicit 中,通用接触默认采用有限滑移(Finite Sliding)与面对面离散,适合大滑动、大变形问题。若接触面之间相对滑动很小,可改用小滑移(Small Sliding),提高计算效率。
四、法向行为参数设置
法向行为常用两种模型:
- 硬接触(Hard Contact):默认模型。当间隙为零时,接触压力可以任意大;当间隙大于零时,压力为零。它不允许穿透,但容易引起收敛困难。
- 软接触(Soft Contact):允许一定穿透,压力随穿透量线性或指数增长。适合需要缓和接触冲击、改善收敛的场合。
在硬接触中,还可设置允许穿透容差(Allowable Interference)和接触刚度(Contact Stiffness)。默认刚度由 Abaqus 根据材料弹性模量和网格尺寸自动估算,若收敛困难,可适当降低刚度,但会增大穿透量。
五、切向行为参数设置
切向行为主要描述摩擦。常用模型是库仑摩擦(Penalty Friction),需要输入摩擦系数 μ。Abaqus 默认使用罚函数法,允许弹性滑移。关键参数包括:
- 摩擦系数 μ:例如钢-钢干燥表面可取 0.1~0.3,具体需根据实验或手册确定。
- 弹性滑移容差(Elastic Slip Tolerance):默认 0.005,表示允许的弹性滑移量占特征长度的比例。值越小,摩擦约束越刚,收敛越难。
- 剪切应力上限(Shear Stress Limit):当接触压力很大时,可设置最大剪应力,避免摩擦应力无限增长。
若接触面之间不允许相对滑动,可使用粗糙摩擦(Rough Friction),此时 μ 为无穷大,从面节点一旦接触即被完全约束。
六、参数设置示例与判断依据
以一个钢制滑块在钢板上滑动的二维模型为例。滑块底面为从面,网格尺寸 0.5 mm;钢板顶面为主面,网格尺寸 2 mm。接触属性设置如下:
*Surface Interaction, name=IntProp-1
*Friction
0.15,
*Surface Behavior, pressure-overclosure=HARD
其中 0.15 为摩擦系数,采用默认弹性滑移容差 0.005。若计算中出现“接触压力振荡”或“穿透过大”,可先检查主从面是否选反,再尝试将法向行为改为软接触,或把弹性滑移容差从 0.005 提高到 0.01。若摩擦应力远大于材料屈服强度,应设置剪切应力上限,例如 *Friction, shear stress limit=100。