4. 刚性弹体侵彻金属薄靶的剩余速度计算
侵彻问题建模要点
刚性弹体侵彻金属薄靶的核心是：弹体保持刚性（不设失效、不设塑性），靶板用 Johnson-Cook 本构描述应变率与温升软化，用 Gruneisen 状态方程描述高压下的体积压缩行为。薄靶侵彻的典型特征是冲塞（plugging）与花瓣开裂，靶板厚度远小于弹径，因此网格在厚度方向至少布置 4~6 层实体单元，才能捕捉剪切冲塞带。
- 靶板 *MAT_JOHNSON_COOK：ρ=7850 kg/m³，G=77 GPa，A=350 MPa，B=275 MPa，n=0.36，C=0.022，m=1.0，Tm=1811 K，Tr=293 K
- 靶板 *EOS_GRUNEISEN：C=4569 m/s，S1=1.49，γ0=2.17，a=0.46
- 弹体 *MAT_RIGID：ρ=7850 kg/m³，E=210 GPa，ν=0.3，约束全部自由度
- 接触：*CONTACT_ERODING_SURFACE_TO_SURFACE，FS=0.2，FD=0.2
模型搭建与求解设置
采用 1/4 对称模型降低计算量。弹体为圆柱平头弹，直径 10 mm，长 30 mm；靶板为 100 mm × 100 mm 方板，厚度 3 mm。对称面施加法向约束，靶板四周固定。弹体初速度沿 Z 向，分别取 500、700、900 m/s 三组进行对比。
*KEYWORD
*CONTROL_TERMINATION
$ TERMINATION TIME
200.0
*CONTROL_TIMESTEP
$ DTINIT TSSFAC ISDO
0.0 0.9 0
*DATABASE_BINARY_D3PLOT
$ DT/CYCL
2.0
*DATABASE_HISTORY_NODE
$ 节点集用于提取弹体速度
NODE_BULLET
*MAT_JOHNSON_COOK
$ MID RO G E PR DT VP
2 7.85 0.77 2.10 0.30 0.0 0.0
$ A B N C M TM TR EPS0
0.0035 0.00275 0.36 0.022 1.0 1811 293 1.0
*EOS_GRUNEISEN
$ EOSID C S1 S2 S3 GAMAO A E0
2 0.4569 1.49 0.0 0.0 2.17 0.46 0.0
*MAT_RIGID
$ MID RO E PR N COUPLE M ALIAS
1 7.85 2.10 0.30 0.0 0.0 0.0
$ CMO CON1 CON2
1.0 7.0 7.0
*CONTACT_ERODING_SURFACE_TO_SURFACE
$ SSID MSID SSTYP MSTYP SBOXID MBOXID SPR MPR
2 1 0 0 0 0 0 0
$ FS FD DC VC VDC PENCHK BT DT
0.2 0.2 0.0 0.0 0.0 0 0.0 0.0
*INITIAL_VELOCITY_GENERATION
$ NSID STYP OMEGA VX VY VZ
1 0 0.0 0.0 0.0 -0.5
*SET_NODE_LIST
$ NSID
1
$ NID1 NID2 ...
1 2 3
*END注意单位制统一为 cm-g-μs，速度 500 m/s 对应 0.05 cm/μs，代码中 VZ 取 -0.5 表示沿 Z 负方向运动。*DATABASE_HISTORY_NODE 指定弹体尾部节点，用于输出速度时程曲线。
剩余速度提取与经验公式对比
求解完成后，从 d3plot 中提取弹体节点速度时程，取弹体穿透靶板后速度趋于稳定的值作为剩余速度 Vr。经典经验公式采用 Recht-Ipson 模型：
Vr = (V0 - Vbl) / (1 + m/M) （V0 > Vbl 时）
Vbl = 靶板临界穿透速度其中 m 为弹体质量，M 为冲塞块质量。对于 3 mm 薄靶，冲塞块质量约为 π/4 × d² × h × ρ = π/4 × 1.0² × 0.3 × 7.85 ≈ 1.85 g，弹体质量约 18.5 g，质量比 m/M ≈ 10。
结果判读与验证
从速度时程曲线可以看到三个阶段：初始恒速段（弹体未接触靶板）、减速段（侵彻过程中弹体受靶板阻力减速）、稳定段（弹体穿透后以恒定剩余速度飞行）。取稳定段平均值即为 Vr。
与 Recht-Ipson 公式对比时需注意：公式中的 Vbl 需通过实验或仿真标定，本例中 3 mm 薄靶的 Vbl 约在 180~220 m/s 量级。若仿真得到的 Vr 与公式偏差超过 10%，应检查以下方面：
- 靶板网格厚度方向层数是否足够（建议 ≥ 4 层）
- 接触刚度是否过大导致弹体异常减速
- Johnson-Cook 参数是否适用于当前应变率范围（薄靶侵彻应变率约 10³~10⁴ s⁻¹）
- 是否启用了侵蚀接触以处理单元删除后的接触面更新