📚 课程目录第4章材料模型与接触碰撞的耦合关系

4. 材料模型与接触碰撞的耦合关系

在 LS-DYNA 接触碰撞建模中，接触算法只负责“力如何从主面传到从面”，而真正决定接触力大小、接触面积演化、能量吸收比例和失效触发时刻的，是材料模型。同一套接触参数（如 SOFT=1、FS=0.1、VDC=20）配在不同材料上，接触响应可能相差数倍。本节把材料本构与接触响应之间的耦合关系拆开讲清楚，重点说明金属塑性、泡沫、橡胶、复合材料四类材料在接触界面上表现出的不同行为，以及应变率效应和失效准则如何改变接触力曲线。

一、耦合的本质：接触力由“接触刚度×材料刚度”共同决定

接触力可近似理解为 F = K_contact · δ，其中穿透量 δ 由接触算法控制，而等效刚度 K 同时受罚函数刚度和材料本身刚度影响。当材料进入塑性、压溃或软化阶段，接触界面上的反力不再随穿透线性增长，接触压力分布也会从“边缘集中”转向“中心均布”。因此材料模型的每一个阶段（弹性、屈服、硬化、损伤、失效）都会在接触力—位移曲线上留下特征拐点。

判断材料与接触是否耦合正确，最直接的方法是提取“接触力—时间”曲线和“接触面积—时间”曲线，与材料应力—应变曲线对照：屈服点应对应接触力增速变缓，失效点应对应接触力骤降或接触面分离。

二、四类材料在接触界面上的典型行为

金属塑性：常用 *MAT_PLASTIC_KINEMATIC（MAT_003）或 *MAT_JOHNSON_COOK（MAT_015）。金属在碰撞中先弹性、后屈服、再硬化，接触力曲线呈“先陡后缓”。若使用随动硬化，包辛格效应会让卸载后再加载的接触力偏低。金属薄壁件压溃时，接触面会随褶皱形成而反复接触—分离，此时接触厚度 SST/MST 和 SOFT 选项对结果影响显著。

泡沫：常用 *MAT_CRUSHABLE_FOAM（MAT_063）或 *MAT_LOW_DENSITY_FOAM（MAT_057）。泡沫的应力—应变曲线有长平台段，接触力在压溃后基本保持恒定，是理想的吸能材料。泡沫接触建模的关键是 LCID 定义的应力—应变曲线必须与接触面积匹配，否则会出现“接触力虚高”。泡沫还容易发生单元负体积，需配合 SOFT=1 和较小的 SLSFAC。

橡胶：常用 *MAT_MOONEY-RIVLIN_RUBBER（MAT_027）或 *MAT_OGDEN_RUBBER（MAT_077）。橡胶几乎不可压缩，接触时体积刚度极大，容易产生接触振荡。建议使用 SOFT=2 或 SOFT=1 配合 PENALTY 缩放，并适当增大 VDC 阻尼。橡胶的接触力—位移曲线呈非线性上翘，接触面积随压缩量快速增大。

复合材料：常用 *MAT_COMPOSITE_DAMAGE（MAT_022）或 *MAT_ENHANCED_COMPOSITE_DAMAGE（MAT_054/055）。复合材料接触响应具有方向性：纤维方向刚度高、基体方向刚度低，接触力随铺层角度变化。层间失效（分层）会改变接触面的局部刚度，需配合 *CONTACT_AUTOMATIC_SINGLE_SURFACE 和 *MAT_COHESIVE 使用。

三、应变率效应：碰撞速度越高，接触力越“硬”

金属和泡沫在高速碰撞下表现出明显的应变率强化。LS-DYNA 中常用 Cowper-Symonds 模型描述：σ_d / σ_s = 1 + (ε̇ / C)^(1/P)。以低碳钢为例，取 C=40 s⁻¹、P=5，当应变率达到 1000 s⁻¹ 时，动态屈服强度约为静态的 1.6 倍。这意味着同样接触参数下，高速碰撞的接触力峰值会显著高于准静态，若仍按静态曲线标定接触刚度，会出现“接触力偏低、穿透偏大”的假象。

泡沫的应变率效应相对较弱，但高密度泡沫在 ε̇ > 100 s⁻¹ 时平台应力可提高 20%～40%。橡胶的应变率效应则表现为刚度随速率上升，需在 *MAT_OGDEN_RUBBER 中通过 LCID 提供不同速率下的曲线。

四、失效准则：接触力骤降与接触面分离的触发点

失效准则决定材料何时“退出”接触计算。金属常用失效应变 FS 或 *MAT_ADD_EROSION；泡沫常用体积应变或压实应变；复合材料常用 Hashin 准则或 Chang-Chang 准则。失效一旦触发，单元被删除，接触面出现“空洞”，接触力骤降，同时可能引发相邻单元的连锁失效。

失效准则设置过松，接触力曲线会“虚高”且不出现下降段；设置过严，则接触力过早跌落，吸能量被低估。建议先用单轴压缩仿真标定失效应变，再放入整车或部件模型。

五、实操示例：泡沫吸能盒接触建模参数

以某泡沫吸能盒为例，尺寸 200 mm × 100 mm × 80 mm，密度 0.06 g/cm³，采用 *MAT_CRUSHABLE_FOAM。关键卡片如下：

*MAT_CRUSHABLE_FOAM
$#   mid        ro         e        pr        lc
         1   6.0e-8    50.0      0.01         1
$#   tssf      damp
       0.5       0.1

*DEFINE_CURVE
$#   lcid      sidr       sfa       sfo
         1         0       1.0       1.0
$#                 x                   y
                 0.0                 0.0
                0.05                 1.2
                0.30                 2.5
                0.55                 4.8
                0.70                12.0
                0.75                25.0

*CONTACT_AUTOMATIC_SINGLE_SURFACE
$#   ssid      msid     sstyp     mstyp
         1         0         0         0
$#     fs        fd        dc        vc
       0.15      0.15       0.0       0.0
$#   soft    sofscl    lcidab    maxpar
         1       0.1         0       0.0

其中 LCID=1 的应力—应变曲线在 0.70 应变处应力快速上升，对应泡沫压实。接触采用 SOFT=1，SOFS=0.1 降低接触刚度，避免泡沫单元负体积。运行后提取接触力曲线，若在 0.70 应变附近出现明显上升拐点，说明压实段被正确捕捉；若接触力在 0.3 应变前就异常升高，应检查 SST 和 SOFS 是否过小。

六、耦合检查清单

检查项	金属	泡沫	橡胶	复合材料
接触刚度是否匹配材料模量	高 E，需较大 SOFS	低 E，SOFS 取 0.05～0.2	不可压缩，SOFT=2 优先	方向性，分区设置
应变率是否启用	必须启用 C、P	高密度泡沫建议启用	按速率曲线	按加载速率
失效准则是否标定	FS 或 ADD_EROSION	压实应变	通常不失效	Hashin / Chang-Chang
接触力曲线特征	屈服后增速变缓	长平台段	非线性上翘	方向相关台阶

材料模型与接触不是两个独立模块，而是通过接触界面上的力—位移关系耦合在一起。建模时应先明确材料在碰撞中的主导阶段（屈服、压溃、大变形还是分层），再据此选择接触算法和参数，最后用接触力曲线与材料曲线对照验证。只有两者匹配，仿真结果才具备工程参考价值。