📚 课程目录第3章 LS-DYNA显式动力学基础与电池包分析适配

3. LS-DYNA显式动力学基础与电池包分析适配

电池包碰撞安全分析在LS-DYNA中属于典型的非线性瞬态动力学问题：碰撞过程持续几十到几百毫秒，涉及大变形、材料失效、接触状态频繁切换以及螺栓预紧、焊点、胶层等多种连接形式。理解显式时间积分的基本机制，并据此调整时间步长、沙漏控制、质量缩放与接触算法，是保证电池包模型"算得动、算得准"的前提。

一、显式时间积分：为什么适合碰撞问题

LS-DYNA以中心差分法求解运动方程。在时刻 t_n 已知节点加速度 a_n 后，程序按 v_n+1/2 = v_n-1/2 + a_n·Δt 更新速度，再按 d_n+1 = d_n + v_n+1/2·Δt 更新位移，随后由位移计算应变与应力，得到下一步的节点力。整个过程不需要组装和求解整体刚度矩阵，也不存在隐式方法的迭代收敛问题，因此对高度非线性、接触状态剧烈变化的电池包碰撞问题具有天然优势。

代价是条件稳定性：显式积分要求时间步长小于临界值，否则解会发散。临界步长由最小单元尺寸、材料声速和密度共同决定，粗略可写为 Δt_crit ≈ L_e/c，其中 c = √(E/ρ)。这意味着模型中一个尺寸过小的单元，会把整个模型的时间步长拖低，直接决定计算成本。

二、时间步长与质量缩放

电池包模型中常见的"小单元"来源包括：为捕捉薄壁挤压变形而细化的吸能盒网格、模组端板圆角处的过渡单元、焊点与螺栓连接的局部细化区、以及胶层或隔热垫的薄壳单元。若最小单元尺寸为0.5 mm，钢材料声速约5000 m/s，则临界步长约0.1 μs；而碰撞总时长若为100 ms，需要约10⁶ 步，计算量难以接受。

质量缩放通过人为增大密度来放大时间步长。LS-DYNA中常用 *CONTROL_TIMESTEP 的 DT2MS 参数控制：DT2MS 为负值时，仅对时间步长小于 |DT2MS| 的单元进行质量缩放，且不改变整体质量；为正值时则对所有单元统一施加。工程上通常先设定目标步长（如1 μs），再通过 *DATABASE_GLSTAT 输出的"added mass"监控质量增加比例。

经验判据：电池包整体质量增加控制在1%以内通常可接受；若某局部区域质量增加超过5%，应优先检查该区域网格尺寸，而不是继续加大缩放。

下面这张图给出从"最小单元尺寸"到"可接受时间步长"的完整参数传递链，以及质量缩放介入的位置。

三、沙漏控制：单点积分的必要代价

为降低计算量，LS-DYNA的壳单元和实体单元默认采用单点积分，即每个单元只在一个积分点计算应力。这会产生零能变形模式——沙漏模式，表现为网格出现锯齿状、棋盘状的非物理变形，且不消耗能量。电池包中的薄壁护板、模组壳体、底护板等大面积壳单元区域最容易出现沙漏。

控制手段主要有三类：一是使用 *CONTROL_HOURGLASS 设置刚度型沙漏控制（IHQ=4），系数 QH 一般取0.05～0.1，对电池包这类以弯曲和压溃为主的工况较合适；二是对关键区域改用全积分单元（ELFORM=2 或 -16），但会显著增加计算量；三是通过细化网格、避免单元过度扭曲来从源头抑制。判断沙漏是否失控，可查看 glstat 中的"hourglass energy"与总内能之比，通常应低于5%。

四、接触算法选择

电池包内部存在大量接触界面：电芯与模组框架、模组与箱体、箱体与底护板、螺栓与法兰面等。LS-DYNA提供多种接触算法，选择依据是接触面的相对运动特征和穿透风险。

接触类型	适用场景	关键参数
`*CONTACT_AUTOMATIC_SINGLE_SURFACE`	自接触、大变形区域，如吸能盒压溃、护板折叠	SOFT=1，防止薄壳穿透
`*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE`	模组与箱体、箱体与底护板等已知配对界面	FS/FD摩擦系数，SFS/SFM缩放
`*CONTACT_AUTOMATIC_NODES_TO_SURFACE`	螺栓头/垫圈与法兰面、点焊连接	可选 `*TIED` 模拟焊点
`*CONTACT_TIED_SHELL_EDGE_TO_SURFACE`	胶层、焊缝、结构胶连接	失效准则 `OPTION=1`

对于电池包，推荐优先使用自动接触，并开启 SOFT=1 或 SOFT=2 以处理薄壁结构。摩擦系数方面，钢-钢界面可取0.1～0.15，钢-铝界面约0.15～0.2，胶层界面则需根据试验数据标定。接触厚度 SLDTHK 和 SSTHK 应结合壳厚设置，避免初始穿透。

五、电池包多材料多尺度建模的特殊要求

电池包不是单一材料结构，而是钢箱体、铝合金模组框架、电芯（等效均质或分层）、导热胶、隔热垫、线束与BMS支架的复合体。建模时需注意：

材料本构匹配：钢箱体可用 *MAT_PLASTIC_KINEMATIC 或 *MAT_024；铝合金模组框架建议用 *MAT_024 配合应变率效应（C、P 参数）；电芯若关注整体变形，可用等效均质材料，若关注内部短路风险，则需分层建模并引入失效准则。
尺度过渡：电芯尺度（毫米级）与整车尺度（米级）相差三个数量级，直接细化会导致时间步长骤降。可采用子模型技术或壳-实体耦合，在碰撞高风险区局部细化。
连接建模：螺栓用 *CONSTRAINED_NODAL_RIGID_BODY 或梁单元简化；焊点用 *MAT_100 或 *CONSTRAINED_SPOTWELD；胶层用 *MAT_184 或内聚单元，并定义失效位移。
初始穿透检查：多材料装配体在网格划分后极易出现初始穿透，需用 *CONTROL_CONTACT 的 IGNORE=1 或 *CONTACT 的 THKOPT 处理，否则会在计算初期产生虚假接触力。

电池包模型中若同时存在壳单元、实体单元和梁单元，接触定义时务必检查从面节点是否落在主面单元上，避免因单元类型不匹配导致接触失效或穿透。

综合来看，电池包显式分析的核心调参逻辑是：先由最小单元尺寸和材料声速确定临界步长，再通过质量缩放把步长提升到可接受范围，同时用沙漏控制抑制单点积分的零能模式，最后根据接触界面的运动特征选择接触算法并设置合理的摩擦与厚度参数。每一步都会影响后续步骤的可行域，因此建议在正式计算前先用小模型或局部模型验证参数组合，再推广到完整电池包模型。