2. 薄壁结构几何与材料建模准备

薄壁吸能构件的仿真精度，往往在网格划分之前就已经被决定。几何简化是否保留了真实的变形路径、单元类型是否匹配壁厚与截面尺度、材料曲线是否覆盖了压溃过程中的大应变与应变率效应，这三件事共同决定了后续显式动力学求解能否复现试验中的渐进折叠。本节围绕薄壁管、帽形梁和多胞截面三类典型构件，说明几何简化、单元选择与网格划分的具体做法。

一、典型吸能构件的几何特征与简化原则

薄壁管（圆管、方管）的典型特征是壁厚远小于截面尺寸，例如外径 60 mm、壁厚 1.5 mm 的圆管，径厚比约 40；帽形梁由帽形截面与平板点焊或激光焊连接而成，其闭合腔体是吸能主体；多胞截面则在单一外轮廓内布置多个胞元，通过增加棱线数量提高单位质量的吸能量。

几何简化的核心原则是保留影响变形路径的特征，删除不影响刚度与质量分布的特征。具体而言：

• 保留：诱导槽、翻边、加强筋、焊点位置、截面棱线圆角（圆角半径影响局部屈曲波长，一般不宜删除）。
• 删除：工艺小孔、倒角小于 0.5 mm 的过渡、非承载的安装凸台、装饰性特征。
• 简化：将焊接简化为刚性连接或等效梁单元，将螺栓连接简化为刚性单元或梁单元加接触。

需要特别注意的是，诱导槽的深度与间距直接决定折叠半波长。若在简化中抹平诱导槽，仿真会得到随机屈曲甚至整体失稳，与试验中规则的轴对称或非轴对称折叠不符。因此诱导槽必须保留，且网格需沿槽方向加密。

二、壳单元与实体单元的选择依据

单元类型的选择取决于壁厚与截面特征尺寸之比，以及是否需要捕捉厚度方向的应力梯度。工程上常用的判断依据如下：

对于薄壁吸能构件，壳单元是首选。LS-DYNA 中常用 *SECTION_SHELL 的 ELFORM=2（Belytschko-Tsay）或 ELFORM=16（全积分壳），前者计算效率高，后者在翘曲敏感区域更稳定。壳单元厚度方向积分点通常取 5 个（NIP=5），以捕捉弯曲与膜力的耦合。若壁厚超过 3 mm 且截面尺寸较小，可考虑实体单元，但需注意沙漏控制与接触厚度定义。

三、中面抽取与网格划分策略

壳单元建模的前提是抽取中面。对于等壁厚薄壁件，中面即壁厚中心面；对于变壁厚或带加强筋的结构，需分别抽取各板的中面并保证连接处共节点或使用刚性连接。中面抽取后应检查：面法向是否一致、自由边是否闭合、是否存在重叠面。

网格划分策略需兼顾精度与效率：

• 沿壁厚方向不划分网格（壳单元已包含厚度），但沿折叠波长方向至少布置 4～6 个单元。折叠半波长约为 λ ≈ 1.5√(r·t)，其中 r 为截面圆角半径，t 为壁厚。
• 诱导槽附近网格尺寸取 0.5～1.0 mm，远离诱导槽区域可放宽至 2～3 mm，形成梯度网格。
• 焊点位置需布置节点，焊点间距通常为 20～40 mm，焊点用 *CONSTRAINED_SPOTWELD 或 *MAT_SPOTWELD 模拟。
• 多胞截面在胞元连接处需保证共节点，避免使用接触代替连接，否则会引入非物理滑移。

四、材料建模准备

薄壁吸能构件常用材料为低碳钢、铝合金和复合材料。金属材料在压溃过程中经历大应变（等效塑性应变可达 0.5～1.0）和应变率效应（应变率 10²～10³ s⁻¹）。LS-DYNA 中推荐使用 *MAT_PIECEWISE_LINEAR_PLASTICITY（MAT_24）或 *MAT_PLASTIC_KINEMATIC（MAT_3），并输入真实应力-塑性应变曲线。应变率效应可通过 Cowper-Symonds 模型参数 C 和 P 描述：

*MAT_PIECEWISE_LINEAR_PLASTICITY
$ MID  RO  E   PR  SIGY  ETAN  FAIL  TDEL
   1  7.85e-9 210000 0.3 220  1000  0.8  0.0
$ C    P
  40   5
$ EPS  ES
  0.0  220
  0.05 280
  0.10 320
  0.20 360
  0.40 400
  0.80 430

其中 RO 为密度（单位：t/mm³），E 为弹性模量（MPa），SIGY 为屈服强度，ETAN 为切线模量，C 和 P 为应变率参数。真实应力-塑性应变曲线应通过单轴拉伸试验获取，并转换为真实应力与塑性应变。对于铝合金，需注意其应变率敏感性较低，C 值通常取 10～20。

完成上述准备后，建议在提交计算前做一次快速检查：壳单元法向是否一致、自由边是否闭合、材料曲线是否覆盖预期最大塑性应变、单位制是否统一为 t-mm-s-MPa。这些检查能避免大部分因建模准备不足导致的求解失败或结果异常。