钽靶板在冲击下层裂过程的数值模拟

王云天，曾祥国，陈华燕，杨 鑫，王 放，祁忠鹏

（1. 四川大学建筑与环境学院深地科学与工程教育部重点实验室，四川 成都 610065；

2. 成都理工大学环境与土木工程学院地质灾害防治与地质环境保护国家重点实验室，四川 成都 610059；3. 西南大学材料与能源学院，重庆 400715）

延性金属在强动载荷下的动态破坏是冲击波物理中重要的基础科学问题之一，在武器制造和装甲防护等特殊环境中具有广泛应用。层裂是延性金属在强动载荷下的一种典型动态破坏形式，不仅包含了金属材料在动态破坏过程中的丰富信息，而且层裂实验原理简单、理论分析成熟，因此一直作为研究材料动态破坏过程的主要方法[1]。Curran 等[2]通过对大量实验结果进行观察，提出延性金属的层裂是由于内部的孔洞（或微损伤）成核、长大以及聚集，最后导致材料内部发生断裂，产生层裂飞片的动力学过程，这也是目前得到广泛认可的延性金属层裂的基本物理内涵。

平板撞击实验具有能够在加载过程中保持一维、平面应变载荷条件以及便于测量和分析回收的特点，是研究材料层裂的常用方式，在层裂实验中，自由面速度曲线是研究层裂演化过程最重要的参考数据之一[3]。图1 为平板撞击实验原理示意图。飞片在t0时刻以一定速度撞击样片，二者相撞时会分别传入冲击波，当冲击波到达样片自由面和飞片后表面时会形成反射稀疏波，由于飞片厚度比样片厚度小，两束稀疏波会在样片内部相遇并形成拉伸区域，如果拉伸应力超过材料的动态承载极限，材料便会发生层裂[2]。通过VISAR（Velocity interferometer system for any reflector）、DISAR（Displacement interferometer system for any reflector）等诊断技术，可以得到样片的自由面速度曲线，从中可以分析得到延性金属材料动态损伤宏观响应的一系列重要信息，如加载应力幅值、拉伸应变率、层裂强度和层裂片厚度等[4]。

Meyers 等[5]系统地总结和评述了20 世纪80 年代之前关于层裂的主要研究结果。2003 年，Antoun 等[3]在其著作中对21 世纪以来的层裂研究历史和现状作了权威性总结。近年来，国内外学者通过实验、理论和数值模拟方法对层裂开展了深入研究。陈永涛等[6]采用等厚对称碰撞实验，对纯铁在冲击压力超过相变阈值下的层裂特性进行了研究，发现样品出现2 次层裂现象。Thomas 等[7]采用对称碰撞的方式对铁在低于冲击相变阈值下的层裂特性进行了研究，碰撞速度在0.2～0.4 km/s 范围内，获得了铁在低速冲击下的Hugoniot 弹性极限与层裂强度数据。翟少栋等[8]利用轻气炮与激光加载等方式对铝的层裂特性进行了研究，分析表明层裂强度具有明显的应变率效应。Koller 等[9]研究了无氧铜在不同冲击波作用下的层裂特性，结果表明冲击波波形与峰值应力以及拉伸脉冲的持续时间对层裂损伤有显著影响。张林等[10]开展了钽在30 GPa 以内冲击波加载下的层裂实验，并对不同的层裂强度计算模型进行了分析讨论。实验研究由于其自身的特点，难以直接获得样品内部物理量随时间的变化数据，只能通过自由面速度等间接获得样品内部层裂信息。而通过数值模拟方法可以提供完整的层裂演化信息，是实验研究的有力补充。目前关于层裂现象的数值仿真技术主要有有限元方法和分子动力学方法。张凤国等[11]基于应力瞬时断裂判据和Tuler-Butcher 判据对延性金属在三角波加载下的多次层裂破坏问题进行了研究，分析表明延性金属层裂片厚度主要由材料本身的物理性质和加载条件决定。种涛等[12]使用一维流体动力学软件模拟了铁的冲击相变与层裂损伤之间的相互影响。贺年丰等[13]采用光滑粒子流体动力学（Smoothed particle hydrodynamics, SPH）方法对爆轰加载下锡的层裂进行了数值模拟，得到的结果与实验结果基本吻合。席涛等[14]通过分子动力学方法研究了超高应变率下铜的层裂特性，获得了铜在109～1010s-1时的层裂强度。Glam 等[15]采用数值模拟方法研究了铝在不同温度下的层裂强度与Hugonoit 弹性极限，结果表明预加热的铝有更高的层裂强度。

图1 平板撞击实验原理及自由面速度曲线示意图Fig. 1 Theoretical illustration of the flat plate impact test and the schematic diagram of the free surface velocity curves

实验研究成本高、周期长、得到的实验数据有限，数值模拟的加载范围与条件更加灵活，并且计算快速、结果直观，已经逐渐成为实验研究的重要补充。基于网格的计算方法对弹塑性材料在中低速冲击下比较可靠，但金属在高速冲击产生的高温高压环境下会表现出流体的性质，在利用传统网格方法进行数值模拟时会遇到一些困难，网格可能会发生严重扭曲，甚至会导致计算无法正常进行。SPH 方法[16]是一种无网格粒子算法，相较于传统网格方法，具有精确追踪自由表面、变形边界等优势。SPH 方法不依赖网格，而是将材料离散为带有质量、能量、动量的运动粒子进行计算，非常适合模拟层裂、剥落及碎片等现象。Libersky 等[17]率先将材料强度效应引入SPH 方法，成功地进行了高速碰撞的数值模拟计算。徐志宏等[18]利用改进的SPH 方法研究了超高速碰撞问题，分析了靶板厚度、弹丸速度和形状对碎片云的影响。Zhou 等[19]则利用SPH 方法开展了铝球超高速撞击铝板的数值模拟，并与实验结果进行了对比。综上所述，SPH 方法在模拟高速冲击中是一种有效、可行的数值方法。数值模拟结果受状态方程、本构模型和失效模型等因素的影响很大，准确合理的模型选择至关重要。

利用平板撞击实验研究金属材料在高速冲击下的动态损伤时，对冲击波的分析一般采用流体力学近似的方法，忽略固体材料抗剪切能力。要想准确地描述金属材料的复杂力学行为，就必须考虑冲击过程中非各向同性应力状态下的应力、应变、应变率和温度的变化，以及它们之间的相互关系。迄今为止，已经有一些材料模型用来描述延性金属在强动载荷下的动态破坏现象，其中被商用有限元软件广泛采用的主要有3 种：Johnson-Cook（JC）模型[20]、Steinberg-Cochran-Guinan（SCG）模型[21]以及Zerilli-Armstrong（ZA）模型[22]。ZA 模型基于不同微观结构（FCC、BCC 结构）特点建立，参数的物理意义明确。JC 模型是一种获得广泛应用的动态本构模型。在JC 模型中假定材料的屈服强度与应变率相关而与压力无关，其模型参数较少，通过一维应力实验可快速获得，一般认为JC 模型在应变率低于104s-1时有较好的适用性。SCG 模型则考虑了材料在强动载荷下压力和温度对剪切模量与屈服强度的影响，忽略材料在高压下的应变率效应，常用于描述高压状态下材料的动态破坏，例如铝在10 GPa 以上时的应变率效应可以忽略，屈服强度和剪切模量只与温度、压力相关[21]，除此之外，SCG 模型在延性金属层裂的数值模拟中也得到广泛应用[4,13,23-24]。

钽具有高密度、高熔点、良好的延展性和热传导性以及耐腐蚀、可焊接性好、高冲击阻抗与高断裂韧性等特点，在国防领域具有广泛应用。本研究选取钽作为研究对象，利用AUTODYN 软件中的Lagrange 模块与SPH 模块，对钽在高速冲击下的层裂行为进行研究，对比JC、ZA、SCG 模型在层裂数值模拟中的结果以及SPH 方法与Lagrange 方法在层裂模拟中的差异，分析自由面速度曲线的物理意义与获得的物理量之间的关系。得到的结果可以为延性金属层裂的数值模拟方法中材料模型的选择和实际工程应用提供参考。

1 计算方法和模型

1.1 本构模型、状态方程与失效模型

1.1.1 JC 模型

JC 模型[20]考虑了延性金属在强动载荷下的应变、应变率硬化和温升软化等因素，假设材料为各向同性，将应变、应变率硬化效应和温升软化效应解耦，广泛应用于表征材料在冲击作用下的动态力学行为[25-26]。JC 模型的表达式为

表1 JC 模型参数[27]Table 1 Parameters of Johnson-Cook model[27]

1.1.2 ZA 模型

Zerilli 等[22]基于热激活和位错动力学理论，建立了面心立方金属（FCC）与体心立方金属（BCC）的本构模型，ZA 模型具有明确的物理意义，综合考虑了应变率与温度的影响，能够准确地描述多种延性金属的动力学行为[28-29]。对于FCC 金属，ZA 模型的表达式为

对于BCC 金属，ZA 模型的表达式为

式中：C1、k1、C2、C3、C4、C5为材料常数，T 为温度， λ为晶粒尺寸。本研究中计算所采用的ZA 模型参数如 表2[22]所示。

表2 ZA 模型参数[22]Table 2 Parameters of Zerilli-Armstrong model[22]

1.1.3 SCG 模型

金属材料的屈服强度Y 与剪切模量G 之间存在比例关系[21,30]，Steinberg 等[21]通过研究发现，对于应变率不敏感的材料，屈服强度Y 与剪切模量G 成正比，即Y 与G 的比值为常数。基于此，Steinberg、Cochran 和Guinan 提出了一种适用于金属材料在高压高应变率下的本构模型，考虑了静水压力对剪切模量和屈服强度的强化效应以及温度软化效应，同时忽略了应变率效应，并假设金属在熔化状态下可以忽略材料强度。SCG 模型中屈服强度Y 和剪切模量G 的关系表示为

式中： η为压缩比， εp为等效塑性应变，T 为温度，p 为压强，G0和Y0分别为初始状态下的剪切模量和屈服强度（T = 300 K，p = 0，ε = 0）， β和n 为材料常数。

屈服强度满足

式中：Ymax为最大屈服强度。本研究中计算所采用的SCG 模型参数如表3[21]所示。JC 模型与ZA 模型中采用的剪切模量、屈服强度、熔点等材料参数与SCG 模型相同。

表3 SCG 模型参数[21]Table 3 Parameters of Steinberg-Cochran-Guinan model[21]

1.1.4 Mie-Grüneisen 状态方程

Mie-Grüneisen 状态方程广泛应用于冲击动力学相关问题的数值模拟研究中，内嵌在AUTODYN与LS-DYNA 等商业软件中，可以较为准确地描述绝大部分金属固体在冲击载荷下的动力学行为。常用于有限元软件中的形式为[13]

式中： ρ为冲击后材料密度， ρ0为初始材料密度，C0为初始波速，S1为材料常数， γ为Grüneisen 系数，e 为质量内能。本研究中计算所采用的Mie-Grüneisen 状态方程参数如表4[31]所示。在计算中，材料密度均选用表4 中给出的数值。

表4 Mie-Grüneisen 状态方程参数[31]Table 4 Parameters of Mie-Grüneisen equation of state[31]

1.1.5 Grady 失效模型

在层裂的数值模拟中，常用最大主应力失效模型与Grady 失效模型。Rinehart[32]将材料静态拉伸断裂模型推广到动态断裂中，提出了最大拉应力失效模型：当材料内部拉应力 σ达到层裂强度 σspall时，材料发生层裂。最大主应力失效模型需预设层裂强度，而层裂强度是一个与冲击过程相关的物理量。Grady[33]基于断裂力学理论提出了层裂的能量平衡破碎模型，将层裂强度 σspall与层裂碎片平均尺寸和应变 率联系起来，层裂强度随冲击加载过程而改变。对于延性金属，层裂强度 σspall的表达式为[34]

式 中：Cb为材料体积声速， εc为材料临界失效应变，金属材料一般取0.15[33]。

1.1.6 层裂强度、拉伸应变率与层裂片厚度

层裂强度是表征材料在强动载荷下性能的一个重要指标，在工程应用中具有重要参考意义，但其在实验中难以直接测量。目前关于如何计算层裂强度还没有统一的方法，常用的方法是通过样片自由面速度曲线来进行计算。Novikov[35]基于声学近似的方法给出了计算层裂强度的公式

式中：Δus为自由面速度曲线最大值与第一次速度拉回时的差值。

通过自由面速度曲线，可以得到样片在层裂过程中的平均拉伸应变率计算公式为[36]

式中：Δts分别为自由面上速度最大值和第一个极小值之间的时间差。

为了能够更系统地了解自由面速度曲线，计算了自由面速度回跳速率，即极小值到峰值速度之间的 斜率。回跳速率的计算公式为

式中：Δur与Δtr分别为自由面上速度第一个极小值与其后峰值之间的速度差和时间差。

通过自由面速度曲线可以获得层裂片厚度dsp，计算公式为

式中：Δt 为第一个层裂振荡周期持续时间，CL为弹性纵波声速。

1.2 研究方法

数值模拟采用在冲击动力学研究中广泛应用的有限元软件AUTODYN[34,37-41]。AUTODYN 包含欧拉（Euler）求解器、拉格朗日(Lagrange)求解器、任意拉格朗日欧拉（ALE）求解器以及SPH 求解器。相对于其他同类软件，AUTODYN 中内嵌了丰富的材料模型及对应参数，便于快速部署应用。本研究中采用Lagrange 与SPH 两种求解器对延性金属在平板撞击下的层裂问题进行分析。Lagrange 求解器对于侵彻、碰撞类低速、中速以下变形问题的模拟具有很好的适应性和精确性。SPH 求解器是近年来快速发展的一种粒子类方法，没有网格束缚，可以有效地避免网格畸变问题，能方便地处理断裂、层裂、破碎等材料破坏问题[13,34,42]。本研究中算例的模型尺寸参数如表5 所示。其中：df为飞片厚度，ds为样片厚度，Ds为样片直径，v 为撞击速度。采用二维轴对称方式建模，有限元网格尺寸为0.05 mm，SPH 粒子尺寸为0.10 mm。平板撞击模型与306 m/s 撞击速度下的结果如图2 所示。由图2 可以看出，层裂区域内的材料已经完全分离，分离界面如图2 中白色区域所示，这表明钽靶板在306 m/s 的撞击速度下发生了完全层裂。

表5 验证模型参数设置Table 5 Parameter settings of simulation validation cases

图2 306 m/s 撞击速度下平板撞击模型与模拟结果Fig. 2 Configuration of the plate impact simulations and simulation results at 306 m/s

2 模拟结果与分析

2.1 本构模型验证

为了分析不同本构模型与求解器对模拟结果的影响，分别采用JC、ZA、SCG 模型与Lagrange 和SPH 求解器对撞击速度为306 m/s、飞片厚度为3.00 mm 的情况进行了模拟计算，得到的自由面速度曲线与实验数据[43]的对比如图3 所示。从图3 中可以看出，（1）在0～2 μs 内，即碰撞初始到自由面速度达到最大值阶段，Lagrange 求解器的模拟结果与实验数据吻合较好，SPH 求解器由于其粒子特性在速度上升段与实验数据稍有不同；二者的区别还体现在，Lagrange 求解器的结果可以观察到明显的Hugoniot 弹性极限信号，而SPH 求解器的结果中Hugoniot 弹性极限信号不是特别明显，可见SCG 与ZA 本构模型对钽的弹塑性行为有较好的描述。除SPH 求解器结合ZA 本构模型得到的自由面速度曲线最大值稍低之外，其他模拟结果都与实验数据非常接近。（2）对比层裂信号出现时的自由面速度，Lagrange 求解器结合JC 模型的偏差最大，差值达到65%，Lagrange 求解器结合ZA、SCG 模型的偏差在25%左右。SPH 求解器得到的结果中，JC 模型偏差最大，约为55%，而ZA、SCG 模型得到的结果与实验数据基本相同。（3）对比层裂信号出现后速度回拉曲线与振荡周期，Lagrange 求解器得到的模拟结果中速度回跳速率快于实验数据。具体来看，SCG 模型的振荡周期与实验数据较为一致，但回跳速率最大值过高，JC、ZA 模型的振荡周期均小于实验数据。SPH 求解器的结果中，JC 模型偏差较大，ZA 模型速度回跳的最大值低于实验数据，SCG 模型与实验数据比较吻合。

图3 不同模型的自由面速度曲线与实验数据[43]的对比（撞击速度306 m/s）Fig. 3 Comparison of free surface velocity profiles between different simulations and experiment data[43] (impact velocity 306 m/s)

为了进一步验证SPH 方法结合SCG 本构模型的准确性，分别对比了撞击速度为212 与412 m/s 下的实验数据，结果如图4 所示。从图4 中可以看出，SPH 求解器结合SCG 本构模型得到的模拟结果与实验数据在各个阶段都有较好的一致性。因此，本研究采用SPH 求解器与SCG 模型开展钽在平板撞击下的层裂行为模拟。为了探究不同拉伸应变率下钽的层裂行为，保持样片厚度4.95 mm、直径50.00 mm不变，通过改变飞片的厚度来改变飞片自由面反射的稀疏波到达样片自由面的距离，改变样片中的冲击波卸载时间，以此达到不同的拉伸应变率，模拟参数和结果如表6 所示。其中：ps为冲击压力，umax为峰值层裂速度。

按照模型参数开展了不同拉伸应变率下的层裂行为模拟，模拟结果如图5 所示，拉伸应变率范围为2.13 × 104～5.40 × 104s-1。从图5 中可以看出，在相同的加载速度下，样片自由面速度曲线峰值平台持续时间逐渐延长。这是由于随着飞片厚度的增大，其自由面反射稀疏波到达样片自由面的距离增大，导致自由面速度衰减的稀疏波波头到达自由面的时间延长，自由面速度的衰减时间逐渐延后；在飞片厚度相同时，自由面速度曲线随加载速度增大而升高，峰值速度平台持续时间基本相同。

图4 不同速度下SCG-SPH 模型模拟结果与实验数据对比Fig. 4 Comparison of SCG-SPH model simulation results and experimental data with different velocities

表6 层裂模型参数与结果Table 6 Parameters of plate impact simulations and results

图5 数值模拟得到的不同加载情况下的自由面速度曲线Fig. 5 Numerical simulation of free surface velocity profiles with different loading conditions

2.2 层裂特性分析

层裂强度是金属材料动态损伤行为的重要参数之一，本研究采用式(11)来计算钽在不同拉应变率下的层裂强度。现有研究表明，在不同加载条件下，金属材料会表现出不同的强度性质[44]。图6 给出了钽的层裂强度随拉伸应变率的变化情况，从图6 中可以看出，层裂强度与加载速度之间没有明确的对应关系，表明层裂强度不是一个仅受加载速度影响的物理量。进一步分析可以发现，层裂强度随拉伸应变率的增长而增加，结合式(12)，可以推断出层裂强度受冲击加载速度和应力脉冲持续时间等因素的综合作用。这表明层裂强度不是一个特定的材料常数，而是随加载情况的不同而改变。

文献[3]给出了某铜材料的层裂强度与拉伸应变率的拟合关系

图6 层裂强度与拉伸应变率的关系Fig. 6 Relationship between spall strength and tensile strain rate

通过数值分析方法，对数值模拟得到的数据进行拟合，得到钽的层裂强度与拉伸应变率之间的关系为

图7 给出了对数坐标下层裂强度与拉伸应变率之间的关系。从图7 可以看出，在对数坐标下钽的层裂强度随拉伸应变率的变化呈线性增长趋势。由于本研究中数值模拟的最高拉伸应变率为5.4 × 104s-1，在图7 中还对比了更高拉伸应变率下的实验数据[45]。更高拉伸应变率下层裂强度与模拟结果呈现相同的增长趋势，并且位于模拟结果的线性拟合曲线上，这再次验证了数值模拟结果的准确性。

除分析层裂强度与拉伸应变率的关系外，图8 给出了层裂强度与回跳速率之间的关系。由图8可知，层裂强度与回跳速率之间近似呈线性增长关系，表明回跳速率与层裂强度存在一定的内在联系。

图7 对数坐标下层裂强度与拉伸应变率的关系Fig. 7 Relationship between spall strength and tensile strain rate in logarithmic coordinates

图8 层裂强度与回跳速率的关系Fig. 8 Relationship between spall strength and bounce rate

本研究中的层裂强度根据式(11)进行计算，这也是目前研究中广泛采用的声学近似计算方法。式(11)进行了一定的简化，并未考虑塑性波和层裂片厚度等因素的影响。Stepanov 等[46]指出，对于弹塑性材料，层裂冲击应以弹性纵波速度CL传播。而在它前面的入射稀疏塑性波应以体积声速Cb传播，当弹性波重要时，可以采用有效声速Ce

将式(11)修正为

在式(11)基础上，Kanel[47]指出层裂强度的计算应考虑层裂片厚度的影响，并给出了一种考虑层裂片厚度影响的修正公式

3 种方法得到的层裂强度计算结果如表7 所示。从表7 可以发现，第2 种方法对第1 种方法的修正可达8%，第3 种方法的计算结果与第2 种相差不大。这表明通过自由面速度曲线来计算层裂强度与选用的计算方法相关，而且不同计算模型之间存在一定差异，差异对结果的影响应该得到关注。

表7 不同计算层裂强度的公式得到的数据对比Table 7 Comparison of data obtained by different formulas for calculating the fracture strength

通过式(14)计算得到样片中的层裂片厚度dsp，结果如表8 所示，其中 δ为误差。由表8 可以看出，通过计算得到的层裂片厚度与飞片厚度基本相同，表明样片发生了完全层裂，这与模拟得到的结果相同，同时也验证了所选模型的准确性。样片发生完全层裂时，自由面速度曲线表现出规则的振荡，如图5中所示，自由面速度曲线振幅逐渐减小，表明层裂片中的脉冲在往返过程中逐渐衰减。完全层裂的应力 波传播过程如图1 所示。

表8 不同撞击速度下层裂片厚度Table 8 Spall scab thickness at different impact velocities

2.3 Hugoniot 弹性极限

Hugoniot 弹性极限描述了在一维应变压缩条件下，金属材料发生非弹性变形时的纵向应力值，是应力-体积状态转折点，也是重要的评价金属材料抵抗动态破坏能力的指标[48]。从自由面速度曲线可以得到钽的Hugoniot 弹性极限

式中：uHEL为Hugoniot 弹性极限出现时对应的自由面速度，对应自由面速度曲线中出现的第一个短平台。表9 给出了由式(21)计算得到的Hugoniot 弹性极限。从表9 中可以看出，Hugoniot 弹性极限随加载条件的改变而改变， σHEL的范围为1.77～2.09 GPa，张林[45]通过实验数据得到钽在冲击压力范围8.9～22.6 GPa 下的 σHEL为1.8～2.6 GPa，表明通过数值模拟得到的Hugoniot 弹性极限是合理的，当然在更高的冲击压力下，钽有更高的Hugoniot 弹性极限。

表9 不同撞击速度下的Hugoniot 弹性极限Table 9 Hugoniot elastic limit at different impact velocities

Hugoniot 弹性极限与冲击压力的关系如图9 所示。通过数值分析，采用指数形式拟合σHEL与冲击压力数据，得到

将式(21)扩展到17 GPa 后，与实验数据进行对比，发现有较好的吻合性。这表明冲击压力越大，钽发生非弹性变形时的纵向应力值越高，即动态屈服强度越高。

2.4 回跳速率分析

图9 不同冲击压力下Hugoniot 弹性极限Fig. 9 Hugoniot elastic limit under different impact pressures

以自由面速度曲线中第一个极小值为起点，至其后第一个峰值速度平台，分析自由面速度曲线的变化情况，如图10 所示。从图10 中可以看出，自由面速度回跳曲线随加载速度和飞片厚度的变化没有明显的变化趋势，结合表6 中计算的加载应变率信息，发现随着加载应变率的提高，曲线逐渐变得陡峭，斜率逐渐增大。

利用式(13) 计算了钽的自由面速度回跳速率，图11 给出了回跳速率与拉伸应变率的关系。从图11可知，回跳速率随着拉伸应变率的升高而增长，二者近似呈线性变化趋势。Kanel 等[49]在对铝和镁的层裂特性分析中，发现了自由面速度回跳速率与损伤演化之间存在正相关的特点，二者呈线性关系

式中：D 为样片中孔隙率的增长速率。结合本研究中钽的数值模拟数据，可以推断钽在层裂过程中样片内部的损伤演化速率同样随着拉伸应变率的增长而增加，与其他延性金属具有相似的动态损伤演化特点。

图10 自由面速度曲线第一个极小值之后的速度曲线Fig. 10 Free surface velocity profiles after the first minimal speed

图11 应变率对钽样品回跳速率的影响Fig. 11 Influence of strain rate on bounce rate of tantalum samples

3 结 论

以延性金属钽为研究对象，通过对不同本构模型结合Lagrange 以及SPH 求解器，对钽在平面冲击下的层裂行为进行了数值模拟，分析了JC、ZA、SCG 本构模型以及Lagrange 和SPH 求解器在数值模拟中的差异，通过改变飞片厚度与加载速度来改变冲击加载应变率，对不同加载条件下的自由面速度曲线特性进行了分析，得到以下主要结论。

（1）JC 模型在模拟中与实验数据相差较大，ZA 与SCG 模型在模拟弹塑性转变过程中有较好的表现，综合对比，SCG 模型得到的模拟结果与实验数据有更好的一致性。

（2）Lagrange 求解器在模拟钽的弹塑性转变过程有较好的表现，能够观测到明显的Hugoniot 弹性极限信号。SPH 求解器得到的结果与实验数据更符合，能够更全面地描述层裂全过程中的自由面速度曲线变化情况，表明利用SPH 方法研究钽的层裂行为是可行的。

（3）钽的层裂强度与加载速度没有明确的对应关系，层裂强度随拉伸应变率的增加而增大，二者近似呈指数关系。进一步对比层裂强度随自由面速度回跳速率的关系，结果表明自由面速度回跳速率也是层裂行为的一种反映，应该得到关注。通过观察发现自由面速度回跳速率与拉伸应变率之间近似为线性关系，随拉伸应变率的增加而增大。