2. SHPB分离式霍普金森压杆仿真与波形校验
为什么SHPB适合用二维轴对称显式模型
SHPB（分离式霍普金森压杆）的核心物理过程是：子弹撞击入射杆产生一个短时应力脉冲，脉冲沿入射杆传播到试样，一部分反射回入射杆形成反射波，另一部分穿过试样进入透射杆形成透射波。由于压杆和试样都是圆柱体，且载荷、几何、边界条件均绕轴线对称，因此可以用二维轴对称模型（LS-DYNA 中 *SECTION_SHELL 配合轴对称选项，或直接用 *SECTION_SOLID 的轴对称单元）把三维问题降为二维，显著减少单元数和计算时间，同时保留波传播的一维特征。
一维应力波理论给出三个基本关系：入射波应变 εI、反射波应变 εR、透射波应变 εT 满足 εT = εI + εR（界面连续条件），试样两端应力差决定其变形。仿真中提取这三个波形并检查它们是否满足该关系，是校验模型是否正确的关键步骤。
建模要点与关键参数
建立二维轴对称模型时，压杆和试样均采用轴对称实体单元。压杆通常用线弹性材料（*MAT_ELASTIC），试样根据研究目的可选用弹性或塑性材料。子弹以初速度 *INITIAL_VELOCITY 施加。接触方面，子弹与入射杆、入射杆与试样、试样与透射杆之间均需定义面面接触（*CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE）。
网格尺寸必须足够细以分辨应力波的高频成分。一般要求沿传播方向单元尺寸 Δx ≤ c/(10·fmax)，其中 c 为杆中纵波波速，fmax 为脉冲最高频率。对于典型钢杆（c≈5100 m/s）和脉冲上升时间约 20 μs，单元尺寸取 1~2 mm 较为合适。
| 参数 | 典型值 | 说明 |
|---|---|---|
| 入射杆/透射杆直径 | 14.5 mm | 常用 SHPB 杆径 |
| 杆材料 | 钢，E=210 GPa，ρ=7850 kg/m³ | 线弹性 |
| 试样直径 | 10 mm | 略小于杆径 |
| 试样长度 | 5 mm | L/D≈0.5 |
| 子弹速度 | 10~20 m/s | 控制应变率 |
| 单元尺寸 | 1 mm | 保证波形分辨率 |
波形提取与校验方法
在入射杆中部和透射杆中部各设置一个历史变量输出点（*DATABASE_HISTORY_SOLID），记录轴向应变 εxx 随时间的变化。仿真结束后，从 d3plot 或 binout 中读取这两个点的应变时程曲线。
校验步骤如下：
- 从入射杆信号中分离出入射波 εI 和反射波 εR（两者时间上分离，入射波先到达）。
- 从透射杆信号中提取透射波 εT。
- 检查 εT 与 εI + εR 是否吻合。若偏差较大，需检查接触定义、网格密度或边界条件。
- 由 εI 和 εR 计算试样两端位移和应变，进而得到应力-应变曲线。
关键代码片段
*KEYWORD
*CONTROL_TERMINATION
$ 终止时间根据波传播时间设定
500.0
*CONTROL_TIMESTEP
$ 时间步长缩放因子
0.9
*SECTION_SOLID
$ 轴对称单元选项
1, 1
*MAT_ELASTIC
$ 钢杆材料
1, 7850.0, 2.1e11, 0.3
*INITIAL_VELOCITY
$ 子弹初速度
2, 0, 15.0, 0, 0, 0
*CONTACT_2D_AUTOMATIC_SURFACE_TO_SURFACE
$ 子弹-入射杆接触
1, 2, 0, 0, 0, 0, 0, 0
*DATABASE_HISTORY_SOLID
$ 入射杆应变片位置
101
$ 透射杆应变片位置
201
*DATABASE_BINARY_D3PLOT
$ 输出间隔
5.0
结果判读与常见问题
仿真完成后，用 LS-PrePost 打开 d3plot，通过 ASCII 菜单读取历史变量曲线。正常波形应满足：入射波为矩形脉冲，反射波与入射波形状相似但幅值较小（弹性试样时反射波很小），透射波幅值介于两者之间。若出现以下情况需排查：
- 透射波幅值异常大：检查试样与透射杆接触是否过刚，或试样材料参数是否合理。
- 波形振荡严重：网格过粗或时间步长过大，需细化网格或减小时间步缩放因子。
- 入射波与反射波无法分离：子弹过长或试样过短，调整几何尺寸使脉冲宽度小于杆长。