3. Taylor杆撞击刚性墙的塑性变形分析

Taylor杆撞击刚性墙：从几何到塑性变形的显式动力学建模

Taylor杆试验是研究材料动态塑性行为最经典的实验之一：一根圆柱形金属杆以一定初速度垂直撞击刚性墙，撞击端迅速镦粗、杆体整体缩短，通过测量撞击前后的长度变化与端面直径，可以反推材料的动态屈服强度。本例使用LS-DYNA建立二维轴对称模型，采用*MAT_PLASTIC_KINEMATIC描述杆体弹塑性行为，通过对比不同撞击速度下的变形模式，理解应变率与惯性效应对塑性变形的影响。

几何模型与网格

Taylor杆为圆柱体，利用轴对称性只需建立1/2截面模型。取杆长L₀ = 5.0 cm，半径R₀ = 0.5 cm，撞击端位于z = 0平面，刚性墙用*RIGIDWALL_PLANAR定义在z = 0处。网格采用二维轴对称实体单元（*SECTION_SHELL配合轴对称选项，或直接使用*SECTION_2D），沿轴向划分100个单元、径向划分20个单元，撞击端附近适当加密以捕捉大变形梯度。

*KEYWORD
*TITLE
Taylor Rod Impact - Plastic Deformation
*CONTROL_TERMINATION
  50.0
*CONTROL_TIMESTEP
  0.0  0.9  0.0  0.0  0.0
*DATABASE_BINARY_D3PLOT
  1.0
*DATABASE_GLSTAT
  0.5
*DATABASE_MATSUM
  0.5
*PART
Rod
         1         1         1
*SECTION_SHELL
         1         1  0.100000       4.0       0.0       0.0       0.0
*MAT_PLASTIC_KINEMATIC
         1  7.85e-3  2.10e-2  0.300000  2.50e-3  1.00e-2  0.0
  0.0       0.0       0.0
*EOS_GRUNEISEN
         1  0.456900       1.490000       0.0       0.0       2.0       0.0
  0.0
*INITIAL_VELOCITY
         1       0.0       0.0     -0.05
*RIGIDWALL_PLANAR
  0.0  0.0  0.0  0.0  0.0  1.0  0.0
  0.0  0.0  0.0
*SET_NODE_LIST
         1
         1         2         3         4         5
*BOUNDARY_SPC_SET
         1         0         0         0         1         1         1
*END

上述关键字中，*MAT_PLASTIC_KINEMATIC的卡片依次为：密度ρ = 7.85×10⁻³ g/cm³、弹性模量E = 2.10×10⁻² Mbar（即210 GPa）、泊松比ν = 0.3、屈服应力σy = 2.50×10⁻³ Mbar（250 MPa）、切线模量Et = 1.00×10⁻² Mbar、硬化参数β = 0（随动硬化）。*INITIAL_VELOCITY中z方向速度为-0.05 cm/μs，即500 m/s。*BOUNDARY_SPC_SET约束对称轴节点的径向位移。

不同撞击速度下的变形模式

Taylor杆的变形模式随撞击速度变化显著。低速撞击时（约100~200 m/s），杆体以均匀镦粗为主，长度缩短较小，端面直径增大有限；中速撞击时（约300~500 m/s），撞击端出现明显鼓胀，形成典型的“蘑菇头”形状，塑性波从撞击端向自由端传播；高速撞击时（约800 m/s以上），撞击端材料发生剧烈塑性流动甚至局部失效，杆体可能出现剪切带或断裂。

结果提取与变形判读

仿真结束后，通过LS-PrePost提取以下关键量：

• 杆体最终长度L_f：取自由端节点z坐标与撞击端节点z坐标之差，与初始长度L₀对比得到长度缩短率(L₀-L_f)/L₀。
• 撞击端直径D_f：取撞击端附近径向最大变形节点的r坐标乘以2，与初始直径D₀对比得到镦粗比D_f/D₀。
• 塑性应变分布：通过*DATABASE_ELOUT输出单元有效塑性应变，观察应变沿轴向的梯度——撞击端应变最大，向自由端递减。
• 速度时程：通过*DATABASE_NODOUT输出撞击端节点速度，观察速度衰减过程，判断塑性波传播时间。

参数影响与模型调试

Taylor杆仿真的关键参数包括：屈服应力σy、切线模量Et、硬化参数β以及撞击速度v₀。其中σy决定塑性变形的起始阈值，Et决定硬化段的斜率，β控制随动硬化与各向同性硬化的比例。对于大多数金属材料，β取0~0.5之间。撞击速度v₀直接影响变形模式，建议以100 m/s为步长进行参数扫描，观察长度缩短率与镦粗比的变化趋势。

学习卡片：Taylor杆仿真快速检查清单