4. 碰撞仿真中的连接与接触建模

在整车碰撞仿真中,白车身由数百个冲压件通过点焊、结构胶、螺栓和铰链连接而成。碰撞载荷下,这些连接部位往往是应力集中、失效起始和能量传递路径改变的关键位置。连接建模是否合理,直接决定碰撞力传递、变形模式和吸能结果的可靠性。本节围绕四类典型连接的简化与精细建模方法,以及通用接触、面面接触、自接触的设置与稳定性调试展开。

一、四类连接的建模策略选择

连接建模的核心判断依据是:该连接在碰撞过程中是否可能失效、是否需要传递弯矩、以及局部细节是否影响整体变形模式。据此可将建模方法分为简化与精细两个层次。

焊点(点焊):简化方法采用 *FASTENER 或 CONNECTOR 单元,在两层壳单元对应节点间建立刚性或弹塑性连接,计算成本低,适合焊点数量大、非关注区域的快速建模。精细方法采用实体单元或 *BEAM 单元配合失效准则,可模拟焊点拉剪失效。常用失效判据为等效塑性应变或力-位移曲线,例如设定焊点直径 6 mm、失效拉力约 5 kN 作为示例参数。当分析涉及焊点撕裂、母材撕裂或接头脱开时,应使用精细模型。

粘接(结构胶):简化方法用 *COHESIVE 单元或 *SURFACE-BASED COHESIVE 行为,在胶层位置建立零厚度或薄层内聚力单元,通过双线性牵引-分离法则描述胶层损伤。关键参数包括法向刚度、切向刚度、界面强度和断裂能。示例:某结构胶法向强度 30 MPa、切向强度 20 MPa、断裂能 2 kJ/m²。精细方法则用实体单元离散胶层厚度(通常 0.2~1 mm),可捕捉胶层厚度方向的应力梯度,但单元尺寸需小于胶层厚度,计算量显著增加。

螺栓:简化方法用 *BEAM 或 *CONNECTOR 连接被连接件,忽略螺纹细节,仅传递力和弯矩。精细方法用实体单元建立螺栓杆、螺母和垫圈,并设置螺栓预紧力(*PRE-TENSION SECTION),可模拟预紧力松弛和螺栓剪切断裂。碰撞中若螺栓可能被剪断,应使用精细模型并定义剪切失效。

铰接:简化方法用 *CONNECTOR 中的 HINGE 类型,释放转动自由度,仅传递力和两个方向的弯矩。精细方法用实体单元建立铰链销轴和耳板,设置接触和间隙,可模拟销轴磨损、耳板孔变形和铰链卡滞。车门铰链、发动机罩铰链在碰撞中可能发生大转动,建议采用精细模型。

四类连接建模:简化 vs 精细 决策与失效路径 焊点 简化:FASTENER / CONNECTOR 刚性或弹塑性,无失效 精细:BEAM/实体 + 失效 拉剪失效,d=6mm,F≈5kN 失效路径:焊点撕裂→母材撕裂 →接头脱开 粘接(结构胶) 简化:COHESIVE 单元 双线性牵引-分离法则 精细:实体单元离散胶层 厚度 0.2~1mm,应力梯度 失效路径:胶层内聚失效→ 界面脱粘→母材撕裂 螺栓 简化:BEAM / CONNECTOR 忽略螺纹,传递力与弯矩 精细:实体杆+螺母+垫圈 预紧力+剪切失效 失效路径:预紧力松弛→ 螺栓剪断→孔壁挤压 铰接 简化:CONNECTOR HINGE 释放转动,传递力与弯矩 精细:实体销轴+耳板 接触+间隙,大转动 失效路径:销轴磨损→ 耳板孔变形→铰链卡滞 接触建模:通用接触 / 面面接触 / 自接触 设置与稳定性 通用接触(General Contact) *CONTACT, GENERAL 自动识别所有外表面 适合整车多部件碰撞 摩擦系数 μ=0.1~0.3 注意:初始穿透检查 面面接触(Surface-to-Surface) *CONTACT PAIR 主面/从面明确指定 适合已知接触区域 可控制接触刚度与穿透 注意:主从面网格匹配 自接触(Self-Contact) *CONTACT, SELF 同一部件自身折叠接触 适合薄壁管压溃、折叠 需设置合理的接触厚度 注意:避免自锁与振荡

二、接触设置的关键参数与操作

接触定义是碰撞仿真中最容易引发计算失败和结果失真的环节。Abaqus/Explicit 中常用三种接触类型,其设置要点如下。

通用接触:通过 *CONTACT, GENERAL 或 Interaction 模块的 General Contact (Explicit) 创建。它会自动搜索模型中所有外表面并建立接触关系,适合整车碰撞中大量部件之间的随机接触。关键设置包括:摩擦系数(钢-钢取 0.1~0.15,钢-铝取 0.2~0.3)、接触厚度(通常取壳单元厚度的 1.0~1.5 倍)、初始穿透检查与修正。通用接触的优点是建模效率高,缺点是接触对数量大时计算成本上升。

面面接触:通过 *CONTACT PAIR 指定主面和从面。主面应选择刚度较大、网格较粗的一侧,从面选择刚度较小、网格较细的一侧。接触属性中需定义法向硬接触和切向摩擦。对于碰撞中可能分离的部件,应设置 *SURFACE BEHAVIOR 允许分离。面面接触适合已知接触区域且需要精确控制接触刚度的场景,如保险杠与壁障、车门与门槛。

自接触:通过 *CONTACT, SELF 定义,用于同一部件或同一集合内部的接触。薄壁吸能盒、前纵梁在压溃过程中会发生自身折叠和贴合,必须开启自接触。设置时需注意接触厚度不能过大,否则会人为增加结构刚度;也不能过小,否则会穿透。建议取壳单元厚度的 1.0~1.2 倍。

三、接触稳定性调试方法

接触引发的数值不稳定主要表现为:接触穿透、接触力振荡、能量不守恒和计算终止。调试时按以下顺序排查。

第一步:检查初始穿透。在 Abaqus/CAE 中使用 Tools → Query → Contact Penetration 查看初始穿透量。若穿透量超过接触厚度的 10%,应调整部件位置或使用 *CONTACT CONTROLS, RESET 消除初始穿透。

第二步:调整接触刚度。默认接触刚度可能过大导致振荡,或过小导致穿透。可通过 *SURFACE BEHAVIOR, PENALTY 中的刚度比例因子调整,建议从 0.1 开始逐步增大至 1.0,观察接触力和穿透量的变化。

第三步:检查能量曲线。在 History Output 中输出 ALLKE(动能)、ALLIE(内能)、ALLSE(弹性应变能)和 ALLPD(塑性耗散能)。正常碰撞中,动能应单调下降,内能单调上升,总能量守恒误差应小于 5%。若动能出现异常反弹或总能量不守恒,说明接触定义有问题。

第四步:使用接触稳定化。对于严重的不稳定,可在 *CONTACT CONTROLS 中设置 STABILIZE 参数,添加粘性阻尼。阻尼系数建议从 1e-4 开始,过大会掩盖真实物理响应。

调试经验:碰撞仿真中,接触问题占计算失败原因的 60% 以上。建议在提交完整模型前,先用 5~10 ms 的短时计算检查接触状态和能量曲线,确认稳定后再进行全时程计算。

四、连接与接触的协同建模示例

以某轿车前纵梁与吸能盒的连接为例:纵梁与吸能盒之间采用 12 个焊点连接,焊点直径 6 mm,采用 *FASTENER 简化建模,同时定义焊点失效拉力 5 kN。吸能盒自身在压溃过程中开启自接触,接触厚度取 1.2 倍壳厚。吸能盒前端与刚性壁障之间采用面面接触,摩擦系数 0.15。计算结果显示:碰撞前 20 ms 吸能盒有序压溃,焊点未失效;20~35 ms 焊点开始逐个失效,纵梁参与吸能;总吸能量约 45 kJ,满足设计目标。若焊点过早失效,应检查焊点失效参数是否偏小,或改用精细焊点模型重新评估。

连接与接触建模没有唯一正确的方案,关键在于根据分析目标、计算资源和精度要求做出合理取舍。简化模型用于快速迭代和趋势分析,精细模型用于失效机理研究和最终验证。两者结合,才能在保证计算效率的同时获得可靠的碰撞安全性能预测结果。