Dyna泡沫负体积如何处理？

Question

在使用LS-DYNA进行泡沫材料仿真时，常遇到“负体积”问题，尤其是在大变形或高应变率工况下。该问题表现为单元体积计算为负值，导致计算中断或结果失真。其主要原因包括网格畸变严重、材料本构模型参数设置不合理（如压碎应力与密实化应变不匹配）、接触算法不稳定或时间步长控制不当。特别是在低密度泡沫中，刚度较低，易发生过度压缩，加剧负体积风险。如何合理设置*MAT_LOW_DENSITY_FOAM*或*MAT_CRUSHABLE_FOAM*参数，并结合*AUTOMATIC_SURFACE_TO_SURFACE*接触与单元删除准则，成为关键。如何在保证计算稳定性的同时有效抑制负体积现象？

Answer 1

一、负体积问题的物理本质与仿真挑战

在使用LS-DYNA进行泡沫材料的大变形动力学仿真过程中，"负体积"（Negative Volume）是一种常见的数值不稳定性现象。其表现为某些有限单元在计算过程中体积变为负值，导致求解器终止运行或输出失真结果。这一问题在低密度泡沫材料中尤为突出，因其刚度低、易压缩，在高应变率加载下极易发生严重网格畸变。

从物理角度看，泡沫材料在受压时经历弹性段、屈服平台段和密实化阶段。若本构模型未能准确描述密实化应变与压碎应力之间的关系，则材料在高度压缩后无法正确响应，进而诱发单元过度扭曲。

负体积的根本原因可归结为：

• 网格质量差或初始尺寸不合理
• 材料模型参数设置不当
• 接触算法引发穿透或非物理反弹
• 显式时间积分步长控制失效

二、材料模型选择与关键参数配置

LS-DYNA提供了多种适用于泡沫类材料的本构模型，其中 *MAT_LOW_DENSITY_FOAM* 和 *MAT_CRUSHABLE_FOAM* 是最常用的选择。二者均基于多线性塑性理论，但适用场景略有不同。

模型类型	适用密度范围	核心参数	推荐使用场景
MAT_LOW_DENSITY_FOAM	< 200 kg/m³	CSS, EPS, LCID	超轻泡沫、缓冲包装
MAT_CRUSHABLE_FOAM	100–600 kg/m³	QEP, PC, GD	汽车吸能结构

以 *MAT_LOW_DENSITY_FOAM* 为例，需特别关注以下参数：

1. CSS (Crush Stress)：必须与实验应力-应变曲线匹配，避免过高估计承载能力
2. EPS (Crush Strain)：定义密实化起始点，通常取0.6~0.8，过小会导致提前硬化
3. LCID：建议通过单轴压缩试验数据拟合生成多段线表
4. NU（泊松比）：低密度泡沫常设为0.0–0.2，防止横向膨胀异常

三、接触算法优化与表面定义策略

接触是诱发负体积的重要外部因素。当使用 *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE* 时，若主从面定义错误或罚函数过大，将导致局部穿透加剧单元畸变。

优化建议包括：

• 确保刚性体作为主面，泡沫作为从面
• 启用 SOFT = 2 使用柔性接触算法降低冲击振荡
• 设置合理的 DEPTH 和 PENMAX 控制穿透量
• 对复杂几何采用 *CONTACT_TIED_SURFACE_TO_SURFACE 预绑定界面

*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$#     SID     SSTYP   MSTYP   SEGID   BCDT    BSTYP   FS      FD
         1       1       1       0       0.0     0       0.7     0.7
$#     COF     BTDT    DCTOL   OFFSET  RCOF    VDC     PENMAX  TSTART
                 0.0     0.0     0.0     0.0     0.0     0.9     0.0

四、时间步长控制与单元删除机制设计

显式动力学求解依赖于稳定的CFL条件。当局部变形剧烈时，自动时间步长可能不足以抑制发散。

可通过如下方式增强稳定性：

• 设置 *CONTROL_TIMESTEP 中的 TSSFAC ≤ 0.8
• 启用 ISDYN = 1 动态缩减因子调整
• 结合单元删除准则，在达到极限状态前移除失效单元

推荐的单元删除逻辑流程如下：

graph TD
            A[开始计算] --> B{单元应变 > EPS_max?}
            B -->|Yes| C[标记为失效]
            B -->|No| D[继续积分]
            C --> E{满足删除条件?}
            E -->|Yes| F[执行单元删除]
            E -->|No| D
            F --> G[更新接触拓扑]
            G --> A

五、综合解决方案与工程实践建议

为系统性抑制负体积现象，应构建多层次的仿真控制框架：

层级	措施	作用
前处理	六面体主导网格 + 尺寸渐变	提升初始网格质量
材料建模	基于实验标定LCID曲线	确保本构真实性
接触设置	AUTOMATIC_S2S + SOFT=2	减少穿透扰动
求解控制	TSSFAC=0.7, ISDYN=1	动态稳定时间步
失效处理	ERODE with MAXSEP	及时清除畸变单元
后处理监控	NSKEW, MINJAC	实时评估网格健康度

此外，建议在初期运行中开启 *DATABASE_BINARY_D3THDT 并监控 VOLUME 场变量变化趋势，定位最早出现负体积的位置，针对性优化局部设置。