• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

修正金属本构模型在超高速撞击模拟中的应用

马坤 李名锐 陈春林 沈子楷 周刚

马坤, 李名锐, 陈春林, 沈子楷, 周刚. 修正金属本构模型在超高速撞击模拟中的应用[J]. 爆炸与冲击, 2022, 42(9): 091406. doi: 10.11883/bzycj-2021-0315
引用本文: 马坤, 李名锐, 陈春林, 沈子楷, 周刚. 修正金属本构模型在超高速撞击模拟中的应用[J]. 爆炸与冲击, 2022, 42(9): 091406. doi: 10.11883/bzycj-2021-0315
MA Kun, LI Mingrui, CHEN Chunlin, SHEN Zikai, ZHOU Gang. The application of a modified constitutive model of metals in the simulation of hypervelocity impact[J]. Explosion And Shock Waves, 2022, 42(9): 091406. doi: 10.11883/bzycj-2021-0315
Citation: MA Kun, LI Mingrui, CHEN Chunlin, SHEN Zikai, ZHOU Gang. The application of a modified constitutive model of metals in the simulation of hypervelocity impact[J]. Explosion And Shock Waves, 2022, 42(9): 091406. doi: 10.11883/bzycj-2021-0315

修正金属本构模型在超高速撞击模拟中的应用

doi: 10.11883/bzycj-2021-0315
基金项目: 国家自然科学基金(11772269,11402213)
详细信息
    作者简介:

    马 坤(1986- ),男,博士研究生,makun@nint.ac.cn

    通讯作者:

    周 刚(1964- ),男,研究员,博士生导师,gzhou@nint.ac.cn

  • 中图分类号: O347.3; TJ012.4

The application of a modified constitutive model of metals in the simulation of hypervelocity impact

  • 摘要: 为更加准确地计算93钨合金弹超高速撞击Q345钢板问题,构建了修正的金属本构模型。引入GRAY三相物态方程描述材料相态变化,采用Johnson-Cook强度模型描述撞击后期材料的力学行为。结合封加波损伤演化模型以及Johnson-Cook失效模型描述不同应力三轴度下材料的拉伸、剪切失效行为;引入曹祥提出的断裂演化模型,描述材料失效后应力归零的过程。通过对比超高速撞击数值模拟结果与实验结果,验证了本构模型的适用性,并进一步分析了典型弹靶撞击条件下破片群的空间分布特征。研究结果表明:基于修正金属本构模型获得的超高速撞击靶板穿孔直径、弹体侵蚀长度、破片群扩展速度结果与实验结果一致;GRAY三相物态方程能够相对准确地给出弹体撞击首层靶板以及剩余弹体、破片群撞击第2层靶板时弹靶材料的熔化情况;封加波损伤演化模型能够准确判断超高速撞击过程中靶板是否产生层裂破坏;综合封加波损伤演化模型、Johnson-Cook失效模型以及曹祥提出的断裂演化模型后,数值模拟获得的破片群撞击后效靶板的穿孔面积与累积数量的统计曲线结果与实验结果一致;获得了典型条件下的柱形93钨弹体超高速撞击Q345靶板破片群空间分布结果,破片群的前端具有较高的质量、轴向动量以及横向动量(绝对值)。
  • 柱形弹体超高速撞击靶板研究是以空间碎片防护[1-2]、反装甲[3]等为背景的研究工作。数值模拟是相关研究手段之一,主要通过离散近似的方式求解带有初边值条件、材料广义应力应变关系的偏微分方程组,进而获得求解域内的力学参量和动力学行为。结果的准确性主要取决于离散求解的算法以及材料的广义应力应变关系(即本构模型)。在算法方面,粒子类无网格法在解决超高速撞击下材料大变形问题优势明显[4-6],目前光滑粒子流体动力学(smooth particle hydrodynamics,SPH)方法应用较为广泛。而在本构模型方面,如何更加准确地反映超高速撞击过程中材料的动力学行为,是力学工作者主要的研究方向,具有重要的学术意义和应用价值。

    金属材料本构模型一般包括物态方程、强度模型和失效模型。在爆炸与冲击动力学问题的数值模拟中,金属材料普遍采用Grüneisen状态方程[7]。对于超高速撞击问题,为了能清楚地获得材料的相态情况,唐蜜[8]构建了GRAY三相完全物态方程并开展了相关研究。在强度模型方面,Johnson-Cook模型[4]和Steinberg模型[9]是超高速撞击数值模拟中普遍使用的强度模型。两个模型各有特点,Johnson-Cook模型基于一维应力实验得到,Steinberg模型基于一维应变实验得到。在超高速撞击后期的靶板惯性扩孔阶段,材料的压力、应变率开始下降,Johnson-Cook强度模型能很好地描述此时的材料动力学行为。在材料失效模型方面,Johnson-Cook失效模型[10-11]在超高速撞击过程中也有较为广泛的应用。但Johnson-Cook失效模型在计算压力远大于屈服强度的层裂问题上,并不能准确计算材料的层裂失效过程。Collé等[12]在材料层裂问题中采用当压力超过失效阈值后偏应力归零而压力不变的方法建立失效模型,处理方法不符合物理实际。Grady失效模型本质上是拉伸失效模型,能够反映材料层裂失效问题,广泛应用于超高速撞击数值模拟中[13-14]。Eftis等[15]基于微孔洞的成核和长大过程建立了损伤演化模型,并模拟了球形弹体超高速撞击问题,获得的靶板的层裂以及断口破坏形貌与实验结果一致。该模型主要考虑的是静水压力产生的层裂问题,没有考虑剪切失效机制。封加波损伤演化模型[16]是具有明确物理意义的描述材料内部微孔洞总体积演化导致层裂问题的材料模型,一般用于模拟飞片撞击、滑移内爆(或外爆)加载柱壳破坏等问题,在超高速撞击问题中的应用较少。Seisson等[17]在数值模拟钢球超高速撞击石墨成坑问题中,钢球材料的失效模型同时考虑了最大拉伸主应力和等效塑性变形两种失效准则。Corbett[18]在计算铝球超高速撞击铝靶问题中,弹体采用了最大拉应力以及最大剪应力失效模型,靶板选用了Johnson-Cook失效模型。Randles等[4]建立了拉伸损伤和Johnson-Cook损伤的联合损伤演化模型,并验证了数值模拟的准确性。因此,同时考虑拉伸和剪切失效是全面且准确的。但还需要注意的是:在大多数的数值模拟中,当材料满足失效准则后,材料立刻失效[19],相应的单元/粒子的屈服强度立刻为零并且不再承受拉伸应力。这将产生大量互相无作用力的单独单元/粒子,这和弹靶撞击形成破片群的物理现象不符。在韧性金属冲击碎裂问题方面,曹祥等[20]以韧性金属长杆动态拉伸碎裂为例,分析了非线性损伤演化过程对韧性碎裂过程的影响规律。该断裂演化模型应用于超高速撞击数值模拟中,将有助于相对准确地计算弹靶材料的碎裂过程。

    在现有的数值模拟研究中,失效模型还不能全面反映超高速撞击下材料在层裂、剪切断裂时的失效机制,并给出失效后的应力归零的断裂演化过程。因此,本文中拟针对超高速撞击问题,在材料本构模型方面,引入GRAY三相物态方程以描述材料相态变化,采用Johnson-Cook强度模型描述靶板后期动态扩孔,结合封加波损伤演化模型[16]以及Johnson-Cook失效模型,描述不同应力三轴度下材料的层裂、剪切失效行为,并引入曹祥断裂演化模型[20],描述材料失效后应力归零的过程。并将该修正金属本构模型应用于柱形93钨超高速撞击Q345钢板数值模拟中。通过与实验结果的对比,验证本构模型的适用性和准确性。初步分析典型弹靶撞击条件下破片群的空间分布特征。

    GRAY物态方程[21]通过初始熔化内能Em1、完全熔化内能Em2、恒定容比热初始内能EGG,将材料的固相-液相区分为固相、熔化相、液相以及热液相。计算中可根据内能E判断材料所处相态,进而给出相应的修正压力pc以及温度T,具体表达式如下。

    EEm1时:

    T=1G{3R+[(3R)2+2G(EE0)]1/122}
    (1)
    pc=12V0(1x)(γeγ0+ax)GT2
    (2)

    式中:x=1V/V0V/V0为比体积,V为当前体积,V0为初始体积;E0为冷能;G=ge/MageMaRγeγ0a均为材料参数。

    Em1EEm2时:

    T=1G{[3R+ν(ΔSα)]+[(3R+ν(ΔSα))2+2G(EE0+ν2δ T(ΔSα))]1/2}
    (3)
    pc=1V0(1x){12(γeγ0+ax)GT2ν(ΔSα)[λTm+(Tνδ T)(γ0ax)]}
    (4)

    式中:Tm为材料熔化温度;δT为温度改变量;v=(EEm1)/(Em2Em1)ΔSα为材料参数。

    Em2EEGG时:

    T=1GαTm{3R+[(3R)2+2(GαTm)×(EE0Tm(ΔS12α))]1/122}
    (5)
    pc=1V0(1x){12(γeγ0+ax)GT2Tm[ΔS12α(1+T2Tm2)](λ+γ0ax)}
    (6)

    式中:λ为与Grüneisen系数相关的参量。

    EEGG时:

    T=1G{(3RΓα)+[(3RΓα)2+2G(EE0Tm(ΔS12α(1Γ2)))]1/122}
    (7)
    pc=1V0(1x){12(γeγ0+ax)GT2Tm[ΔS12α(1Γ2+2ΓTTm)](λ+γ0ax)}
    (8)

    式中:Γ为材料参数。

    在此基础上,结合线性Grüneisen压力p1,以及修正压力项pcpcc,可获得材料的实际压力p

    p=p1+pc+pcc
    (9)
    p1(x,E)=c20xV0(1x)(1sx)2[1(1+12γ0)x+12ax2]+γ0axV0(1x)E
    (10)
    pcc=12V0T20G(γ0γe)
    (11)

    式中:T0为材料初始温度,c0sγ0为材料的冲击雨贡纽参数。

    在液-气相区方面,GRAY物态方程基于Young和Alder的修正气体物态方程描述气相区,在液气共存区采用Maxwell等面积法进行处理。那么在液-气相区需满足条件V>VJVJ为搭接比体积。液-气相区的压力和温度采用下式计算:

    p=RTVbz(1+z+z2z3)(1z)3ayV2bz2+FP(C1+C2T+C3T2)
    (12)
    T=14(C3FeH3){(3R+2H2)[(3R+2H2)216(E+ayVbzC1FeH1)(C3FeH3)]1/2}
    (13)
    FP=[z(θzJ)zJ(θz)]3
    (14)
    Fe=Vb2[(θzJ)3(2z2+θ)z3J(θz)2(θzJ)(2zJ2+θ)z3J]
    (15)

    式中:z=Vb/Vay=ay/MaVbayC1C2C3H1H2H3θzJ均为材料参数。

    Johnson-Cook强度模型主要考虑了材料在冲击动力学问题中面临的大变形、高温、高应变率问题,该模型假设材料的应变硬化、应变率效应以及热软化效应相互独立可解耦。屈服强度Y的表达式为:

    Y=(Q1+Q2εQ3p)(1+Q4ln˙εp˙ε0)[1(TT0TmT0)Q5]
    (16)

    式中:εp为等效塑性应变;˙εp为等效塑性应变率;˙ε0为参考应变率;T为材料温度,T0为参考温度,Tm为材料熔化温度;Q1Q2Q3Q4Q5为材料参数。

    Johnson-Cook失效模型主要包含了应力三轴度、应变率以及温度对材料失效的影响。失效应变的表达式为:

    εf=[D1+D2exp(D3σ)](1+D4ln˙εp˙ε0)[1+D5(TT0TmT0)]
    (17)

    式中:D1D5为材料参数;σ = p/σeq为应力三轴度,其中p为压力,σeq为等效应力。模型假设当材料发生塑性变形时损伤同时产生积累,引入损伤变量DY,单元的损伤演化定义为:

    DY=Δεeq/εf
    (18)

    式中:Δεeq为每一次损伤积累时的等效塑性应变增量,εf取损伤积累时应力状态下的断裂失效应变。DY的初始值为0,当DY达到1时,材料失效。在此之前的损伤发展过程中,材料强度由强度模型决定,不受失效模型影响。

    为了避免材料失效后应力立刻归零,引入曹祥提出的非线性断裂演化模型[20]。该模型以线性内聚力断裂假设Mott-Grady模型为基础,描述了材料在满足失效准则后的断裂(裂纹软化)过程。即材料在达到失效阈值后,材料应力存在一个逐渐卸载至零的过程。此时损伤屈服强度Y与内聚力断裂假设下的损伤因子Dcs、材料初始失效时的等效应力Y0的关系为:

    Y=(1Dcs)Y0
    (19)

    损伤因子假设为:

    Dcs=1eβ(upl/uplufplufpl)1eβ
    (20)

    式中:upl为损伤启动后单元的塑性位移;β为内聚力断裂的非线性系数,uplf为完全断裂时的塑性位移,其计算公式为:

    uplf=WeY0(1β1eβ1)1β0
    (21)

    式中:We为材料断裂过程中消耗的表面能,即材料的断裂能。因此,材料的损伤演化过程主要取决于材料断裂能、非线性系数以及单元的尺寸。

    封加波等[16]将材料宏观微元中的微损伤(微裂纹、微孔洞)总体积与宏观微元体积之比定义为损伤度函数DV。基于Griffith-Orowan提出的裂纹扩展能量观点,假设微孔洞长大过程中因表面积增加所需要的表面能由微孔洞周围介质所存储的应变能提供。此外,提供微孔洞动态长大的周围介质范围由应力波扰动传播的速度即材料声速来确定。损伤度函数的演化方程为:

    ˙DV={DV(ps2σ02)(4/3)B(ω/c0)DVV˙V|ps|σ00|ps|σ0
    (22)

    式中:ps为拉伸状态下基体材料的压力,σ0为孔洞长大的初始应力阈值,B为基体材料的体积模量,ω为表面功(断裂比功)。

    在材料损伤的后期,微孔洞相互间开始连通,材料的损伤快速增长,材料发生分离性断裂,并且所承受的应力在断裂发生后下降到零。为近似描述材料损伤后期的应力状态,仍采用曹祥等[20]提出的断裂演化模型,含损伤的本构关系定义为:

    p=ps[V(1DV),e](1DV)(1Dcs)
    (23)
    G=Gs(1DV15B+20Gs9B+8Gs)(1Dcs)
    (24)
    Y=Ys(14DV)(1Dcs)
    (25)
    Dcs=1eβ(vpl/vplvfplvfpl)1eβ
    (26)

    式中:e为质量内能,G为剪切模量,Gs为基体材料的剪切模量,vpl为损伤启动后当前单元体积相对初始损伤时的膨胀值,vplf为完全断裂时的单元体积相对初始损伤时的膨胀值,其计算公式为:

    vplf=Wep0(1β1eβ1)1β0
    (27)

    式中:p0为单元损伤启动时的静水压力。代表材料成核、长大过程的损伤度DV的初始值为D0,代表材料微孔洞聚集并导致最终破坏的损伤函数Dcs的初始值为0。材料的损伤演化过程是:首先,Dcs为0并且不变,DVD0起始随损伤演化逐渐增加,当达到Dc时,DV值不再变化,Dcs开始随损伤演化而逐渐增加,当Dcs增加到1时表示材料完全层裂破坏。

    选用ANSYS AUTODYN 3D中的光滑粒子流体动力学(SPH)方法模拟超高速撞击过程。根据软件接口编写修正金属本构模型的自定义子程序模块,再通过数值模拟结果和实验结果的对比来验证本构模型的准确性。选取柱形93钨弹体超高速撞击Q345靶板问题进行计算,材料模型参数见表13。GRAY物态方程的其他材料参数可参考Royce[21]的研究工作。

    表  1  93钨合金、Q345钢GRAY物态方程及损伤演化模型材料参数
    Table  1.  The material parameters of the GRAY equation of state and the damage evolution model for 93 tungsten alloy and Q345 steel
    材料ρ0/(kg·cm−3)Sγ0c0/(m·s−1)ageTm0VJ/V0Vb/V0
    93钨合金17.61.23[22]1.67[22]4040[22]1.231 0904.521.340.50
    Q345钢7.831.49[23]2.17[23]4569[23]1.49 9602.381.400.51
    材料ayMa/(g·mol−1)ω/(J·m−2)σ0/GPaD0DcβWe/(kJ·m−2)
    93钨合金10.2183.85 3001.41.0×10−40.0110.024
    Q345钢4.4 55.382 0500.71.0×10−40.0210.041
    下载: 导出CSV 
    | 显示表格
    表  2  93钨合金、Q345钢Johnson-Cook强度模型材料参数
    Table  2.  The material parameters of the Johnson-Cook yield criterion for 93 tungsten alloy and Q345 steel
    材料G/GPaQ1/MPaQ2/MPaQ3Q4Q5Tm/K˙ε0/s−1
    93钨合金132.00600.812000.49440.0590.820336830.001
    Q345钢[24]80.47374795.70.45450.015860.885617950.001
    下载: 导出CSV 
    | 显示表格
    表  3  93钨合金、Q345钢Johnson-Cook失效模型材料参数
    Table  3.  The material parameters of the Johnson-Cook failure model for 93 tungsten alloy and Q345 steel
    材料D1D2D3D4D5Tm/K˙ε0/s−1
    93钨合金0.091000.40900−2.4600−0.218001.61336830.001
    Q345钢0.094361.15059−3.05970.016450.24017950.001
    下载: 导出CSV 
    | 显示表格

    上述Johnson-Cook强度模型、Johnson-Cook失效模型以及封加波损伤演化模型对应的材料参数中,除部分参数引用了相关文献外,其他参数均为实验所得。包括SHPB实验、缺口拉伸实验、飞片撞击实验等。

    针对柱形93钨弹体超高速撞击Q345钢板问题,可通过实验后效以及高速摄像获得靶板穿孔直径、破片群扩展速度、弹体剩余长度的测量结果。因此,根据实验弹靶结构和撞击条件,建立数值模拟模型(光滑长度为0.1 mm)并开展计算。获得的相应参量的数值模拟结果与实验结果一同列于表45中。其中,vx,max为破片群轴向最大速度,vr,max为破片群横向最大速度,ΔL为弹体侵蚀长度。需要说明的是,实验1~6中未使用高速摄像,实验8~9中受弹托材料撞击影响无法获得靶板穿孔直径数据。

    表  4  靶板穿孔直径的数值模拟结果与实验结果对比
    Table  4.  Comparison of numerical simulation results and experimental results of perforation diameter of the target plate
    实验弹直径/mm弹长/mm靶厚/mm撞击速度/(m·s−1)实验穿孔直径/mm数值模拟穿孔直径/mm误差/%
    12.3223.21.528475.3575.4932.54
    22.9214.61.531406.4206.7615.31
    32.3223.21.523005.0084.907−2.02
    42.9214.61.521805.5665.7553.40
    53.4510.51.530006.9177.2775.20
    63.4510.51.530106.8087.2796.92
    73.4510.51.531606.8267.3707.97
    下载: 导出CSV 
    | 显示表格
    表  5  破片群扩展速度、弹体侵蚀长度数值模拟结果与实验结果对比
    Table  5.  Comparison of numerical simulation results and experimental results of fragment group expansion speed and erosion length of projectile
    实验弹直径/
    mm
    弹长/
    mm
    靶厚/
    mm
    初速度/
    (m·s−1)
    实验数值模拟误差
    vx,max/(m·s−1)vr,max/(m·s−1)ΔL/mmvx,max/(m·s−1)vr,max/(m·s−1)ΔL/mmε(vx,max) /%ε(vr,max)/%εL)/%
    73.4510.51.5316031438373.4031168753.60−0.864.545.88
    82.9214.61.5297529537383.2129477802.99−0.205.69−6.85
    92.3223.22.0218021234603.1021494582.841.22−0.44−8.39
    下载: 导出CSV 
    | 显示表格

    对比表45中数据可知,在靶板穿孔直径、破片群扩展速度、弹体剩余长度方面,数值模拟结果和实验结果一致性较好,大部分误差在8%以内。

    实验研究发现,所有靶板厚度(tt)小于0.72倍[25]弹体直径(dp)的实验中,弹靶撞击后将在靶后破片群前部形成“尖端”。主要是因为弹靶撞击产生向靶板内部传播的强冲击波来不及被来自弹体头部侧向自由面、靶板正向自由面的卸载波及时追赶卸载,导致冲击波在靶板背表面反射形成较强的拉伸波,进而使靶板产生拉伸层裂破坏并形成速度较高的“尖端”破片。而在数值模拟中,弹靶材料模型包含反映材料层裂失效机制的损伤演化模型。因此,数值模拟应能够有效预测靶板是否产生与实验相同的层裂结果。根据实验条件建模并开展数值模拟,对比靶后破片群形貌数值模拟结果与实验结果,如表6所示。

    表  6  不同靶厚与弹体直径比的破片群实验结果与数值模拟结果对比
    Table  6.  Comparison of experimental results and numerical simulation results of fragment group with different target thickness and projectile diameter ratio
    实验tt /dp是否有“尖端”高速摄像结果数值模拟结果
    70.435
    80.514
    90.862
    下载: 导出CSV 
    | 显示表格

    表6可以看出,在靶后破片群前部是否形成“尖端”问题上,数值模拟能够获得同实验一致的现象,这一定程度上验证了数值模拟中靶板材料模型所包含的反映材料层裂失效机制的损伤演化模型的准确性。

    为验证修正金属本构模型中GRAY物态方程计算材料相态的准确性,数值模拟了弹体超高速撞击双层靶板问题,分析了首层靶板穿孔及第2层靶板破片撞击坑的材料熔化情况。弹体尺寸为3.5 mm×17.5 mm,撞击速度为2.76 km/s,两层靶板厚度均为1.5 mm,建立三维面对称模型,光滑长度为0.1 mm。

    弹体穿透首层靶板后,破片群形貌如图1所示,图中标记了材料所处相态,蓝色为固相,绿色为固液混合相,红色为液相。弹体撞击首层靶后,有部分粒子发生熔化,但总体上,破片群没有发生大范围的熔化。破片群“尖端”部分的破片均为固相。

    图  1  数值模拟首层靶后破片群形貌及熔化状态
    Figure  1.  The morphology and melting state of the fragment group after the first layer of target obtained by numerical simulation

    对实验结果的分析表明,在该速度段下,首层靶后破片群没有产生大范围材料熔化。而当剩余弹体及破片群撞击第2层靶板后,受二次加载作用影响,撞击区附近将发生大范围的材料熔化。典型实验获得的首层、第2层靶板中心穿孔断口扫描电子显微镜(SEM)结果如图2所示,第1层靶板断口较为光滑,未发现大范围熔化,而第2层靶板断口则有较多熔融后凝固组织。在图1的数值模拟中,能够观察到首层靶板断口、整个破片群尚未发生大范围的材料熔化。

    图  2  实验第1层和第2层靶板断口的SEM结果
    Figure  2.  SEM results of fracture surface of the first and second layers of the target plates

    在数值模拟中,当剩余弹体穿透第2层靶板时,除靶板中心区域受剩余弹体撞击外,其他区域均为破片群撞击影响区。如图3所示,在破片群撞击后的撞击坑处,大量的SPH粒子发生熔化。对比图1中撞击前的破片群的熔化状态,可以判断在图3中靶板破片撞击坑中的SPH粒子的熔化,大部分是破片在撞击第2层靶板过程中发生的。

    图  3  数值模拟弹体穿透第2层靶板时弹靶位置及靶板撞击坑熔化状态
    Figure  3.  Numerical simulation of projectile target position and melting state of impact crater when projectile penetrates the second layer of the target plate

    上述关于材料熔化的数值模拟分析结果和实验中靶板破片撞击坑微观组织分析的结论是一致的,实验通过扫描电镜获得的破片撞击坑的形貌如图4所示。可以看出,撞击坑表面圆润光滑,像液体浇筑后的表面形貌,说明撞击坑中的材料已经发生了熔化。因此,认为数值模拟基本能够反映超高速撞击过程中弹靶材料所处的相态。

    图  4  实验第2层靶板破片撞击坑SEM形貌图
    Figure  4.  SEM morphologies of impact crater of fragments in the second layer of target

    当93钨合金弹体以2~3.5 km/s速度撞击钢板时,撞击后弹靶材料破碎将形成平均尺寸较大的破片群。为了获得破片群特征尺寸的分布规律,可在首层靶板后设置一层薄效应靶。通过统计效应靶穿孔面积与累积数量的关系,间接分析破片群特征尺寸分布规律以及材料的动态破碎行为。因此,为了进一步验证修正金属本构模型的计算结果,参照实验条件建模并开展计算。弹体直径15 mm、长15 mm。首层靶板厚度10 mm,后效靶板厚度1 mm,两层靶板间距1014 mm。实验弹体速度为2.25 km/s,模型光滑长度为0.25 mm。选取数值模拟弹体撞击首层靶板后1 ms时刻后效靶板的穿孔结果进行分析,实验和数值模拟效应靶穿孔形貌对比如图5所示。图6为效应靶穿孔的累积数量与穿孔面积曲线的对比情况。其中,累积数量表示大于横坐标穿孔面积的穿孔总数。需要说明的是,实验中,弹体撞击首层靶板后形成的破片群经历了与空气介质短暂的热量交换才撞击到薄效应靶上,而数值模拟中,则考虑短暂热交换所引起的温度变化对破片穿孔尺寸的影响相对破片密度以及速度对穿孔尺寸的影响是次要因素,因而忽略了空气的影响。

    图  5  实验与数值模拟后效靶穿孔形貌对比
    Figure  5.  Comparison of experiment and numerical simulation on the perforation morphology of aftereffect target
    图  6  实验与数值模拟后效靶穿孔累积数量与穿孔面积统计曲线对比
    Figure  6.  Comparison of experiment and numerical simulation on the statistical curves of cumulative perforation number and perforation area

    在效应靶穿孔形貌上,由于实验难以保证绝对的轴对称撞击,并且材料破碎有一定的随机性,导致实验效应靶穿孔形貌对称性不明显,而数值模拟有明显的对称性。但总的来说,较大破片穿孔更趋近于效应靶中心区域,效应靶边缘区域穿孔较小,这是实验以及数值模拟的共同特征。从图6可以看出,后效靶穿孔累积数量与穿孔面积统计曲线的数值模拟结果和实验结果一致性较好,仅在统计曲线两端略有差距。那么从统计角度讲,数值模拟获得的超高速撞击破片群符合实验破片群的特征尺寸分布规律。

    开展典型93钨弹体超高速撞击Q345钢板问题的数值模拟,弹体直径3.2 mm,长径比5,靶板厚度1.5 mm,撞击速度为3.0 km/s,模型光滑长度取0.1 mm。保存模拟结果中的SPH粒子位置坐标,参照Liang等[26]所使用的广度优先搜索(BFS)算法,识别破片并统计分析相关物理量的空间分布情况。

    图78为破片群的数量和质量在各轴向比速度段、横向比速度段内的累计情况。在每个柱形图内,不同颜色代表不同质量范围破片的累积结果。vx为破片轴向速度,v0为弹体初始速度,vr为破片横向速度,vr,max为破片群的最大横向速度。从图7可以看出,在轴向方向上破片群前端的破片数量高、总质量占比高。而低速段的大质量破片主要为靶板扩孔末期剥落的低速破片。从图8中可以看出,在横向方向上,各比速度段的数量基本相同。大质量破片主要集中在低横向速度段内,即破片群的前端和尾端。

    图  7  破片数量、质量关于破片群轴向比速度的分布
    Figure  7.  Distribution of fragments number and mass in relation to fragment specific axial velocity
    图  8  破片数量、质量关于破片群横向比速度的分布
    Figure  8.  Distribution of fragments number and mass in relation to fragment specific lateral velocity

    图9为破片轴向动量、横向动量(绝对值)在各轴向比速度段内的累积结果。其中,不同颜色代表不同动量范围内的累积情况。可以看出,随轴向比速度的增大,轴向动量呈指数增长形式分布,高轴向动量破片集中分布于破片群前端。另一方面,高横向动量也集中于破片群的前端,结合破片群质量分布特征认为,这是由于该区域的破片群累积质量较高。综上,破片群的前端具有较高的质量、轴向动量以及横向动量(绝对值)。

    图  9  破片轴向动量、横向动量关于破片群轴向比速度的分布
    Figure  9.  Distributions of fragments axial momentum and lateral momentum in relation to fragment specific axial velocity

    为建立适用于超高速撞击问题的修正金属本构模型,引入GRAY三相物态方程描述材料相态变化,采用Johnson-Cook强度模型描述撞击后期材料的力学行为,综合封加波损伤演化模型以及Johnson-Cook失效模型描述层裂、剪切失效机制,并引用曹祥断裂演化模型,描述材料失效后应力归零的过程。通过柱形93钨超高速撞击Q345钢板数值模拟结果与实验结果的对比,验证了本构模型的适用性。分析了典型弹靶撞击条件下破片群的空间分布特征,得出结论如下。

    (1)基于修正金属本构模型,开展了不同撞击条件的柱形93钨超高速撞击Q345钢板数值模拟,获得的靶板穿孔直径、弹体侵蚀长度、破片群扩展速度结果与实验结果一致。

    (2)GRAY三相物态方程能够相对准确地计算93钨合金撞击Q345钢靶过程中弹靶材料所处的相态。数值模拟发现,撞击速度在2~3.5 km/s时,弹靶材料并不会发生大范围熔化,但当剩余弹体和破片群撞击第2层靶板时,破片撞击坑将产生大范围熔化,该现象与实验结果是一致的。

    (3)封加波损伤演化模型能够准确得出超高速撞击过程中靶板是否产生层裂破坏,对于不同靶板厚度与弹体直径比值的撞击问题,数值模拟获得的靶板是否产生层裂破坏的结果与实验结果一致。

    (4)综合封加波损伤演化模型、Johnson-Cook失效模型以及曹祥提出的断裂演化模型后,典型条件下数值模拟获得的破片群对后效靶板穿孔的面积与累积数量统计曲线和实验结果一致。

    (5)基于修正金属本构模型,获得了典型条件下柱形93钨弹体超高速撞击Q345靶板破片群空间分布统计结果,破片群的前端具有较高的质量、轴向动量以及横向动量(绝对值)。

  • 图  1  数值模拟首层靶后破片群形貌及熔化状态

    Figure  1.  The morphology and melting state of the fragment group after the first layer of target obtained by numerical simulation

    图  2  实验第1层和第2层靶板断口的SEM结果

    Figure  2.  SEM results of fracture surface of the first and second layers of the target plates

    图  3  数值模拟弹体穿透第2层靶板时弹靶位置及靶板撞击坑熔化状态

    Figure  3.  Numerical simulation of projectile target position and melting state of impact crater when projectile penetrates the second layer of the target plate

    图  4  实验第2层靶板破片撞击坑SEM形貌图

    Figure  4.  SEM morphologies of impact crater of fragments in the second layer of target

    图  5  实验与数值模拟后效靶穿孔形貌对比

    Figure  5.  Comparison of experiment and numerical simulation on the perforation morphology of aftereffect target

    图  6  实验与数值模拟后效靶穿孔累积数量与穿孔面积统计曲线对比

    Figure  6.  Comparison of experiment and numerical simulation on the statistical curves of cumulative perforation number and perforation area

    图  7  破片数量、质量关于破片群轴向比速度的分布

    Figure  7.  Distribution of fragments number and mass in relation to fragment specific axial velocity

    图  8  破片数量、质量关于破片群横向比速度的分布

    Figure  8.  Distribution of fragments number and mass in relation to fragment specific lateral velocity

    图  9  破片轴向动量、横向动量关于破片群轴向比速度的分布

    Figure  9.  Distributions of fragments axial momentum and lateral momentum in relation to fragment specific axial velocity

    表  1  93钨合金、Q345钢GRAY物态方程及损伤演化模型材料参数

    Table  1.   The material parameters of the GRAY equation of state and the damage evolution model for 93 tungsten alloy and Q345 steel

    材料ρ0/(kg·cm−3)Sγ0c0/(m·s−1)ageTm0VJ/V0Vb/V0
    93钨合金17.61.23[22]1.67[22]4040[22]1.231 0904.521.340.50
    Q345钢7.831.49[23]2.17[23]4569[23]1.49 9602.381.400.51
    材料ayMa/(g·mol−1)ω/(J·m−2)σ0/GPaD0DcβWe/(kJ·m−2)
    93钨合金10.2183.85 3001.41.0×10−40.0110.024
    Q345钢4.4 55.382 0500.71.0×10−40.0210.041
    下载: 导出CSV

    表  2  93钨合金、Q345钢Johnson-Cook强度模型材料参数

    Table  2.   The material parameters of the Johnson-Cook yield criterion for 93 tungsten alloy and Q345 steel

    材料G/GPaQ1/MPaQ2/MPaQ3Q4Q5Tm/K˙ε0/s−1
    93钨合金132.00600.812000.49440.0590.820336830.001
    Q345钢[24]80.47374795.70.45450.015860.885617950.001
    下载: 导出CSV

    表  3  93钨合金、Q345钢Johnson-Cook失效模型材料参数

    Table  3.   The material parameters of the Johnson-Cook failure model for 93 tungsten alloy and Q345 steel

    材料D1D2D3D4D5Tm/K˙ε0/s−1
    93钨合金0.091000.40900−2.4600−0.218001.61336830.001
    Q345钢0.094361.15059−3.05970.016450.24017950.001
    下载: 导出CSV

    表  4  靶板穿孔直径的数值模拟结果与实验结果对比

    Table  4.   Comparison of numerical simulation results and experimental results of perforation diameter of the target plate

    实验弹直径/mm弹长/mm靶厚/mm撞击速度/(m·s−1)实验穿孔直径/mm数值模拟穿孔直径/mm误差/%
    12.3223.21.528475.3575.4932.54
    22.9214.61.531406.4206.7615.31
    32.3223.21.523005.0084.907−2.02
    42.9214.61.521805.5665.7553.40
    53.4510.51.530006.9177.2775.20
    63.4510.51.530106.8087.2796.92
    73.4510.51.531606.8267.3707.97
    下载: 导出CSV

    表  5  破片群扩展速度、弹体侵蚀长度数值模拟结果与实验结果对比

    Table  5.   Comparison of numerical simulation results and experimental results of fragment group expansion speed and erosion length of projectile

    实验弹直径/
    mm
    弹长/
    mm
    靶厚/
    mm
    初速度/
    (m·s−1)
    实验数值模拟误差
    vx,max/(m·s−1)vr,max/(m·s−1)ΔL/mmvx,max/(m·s−1)vr,max/(m·s−1)ΔL/mmε(vx,max) /%ε(vr,max)/%εL)/%
    73.4510.51.5316031438373.4031168753.60−0.864.545.88
    82.9214.61.5297529537383.2129477802.99−0.205.69−6.85
    92.3223.22.0218021234603.1021494582.841.22−0.44−8.39
    下载: 导出CSV

    表  6  不同靶厚与弹体直径比的破片群实验结果与数值模拟结果对比

    Table  6.   Comparison of experimental results and numerical simulation results of fragment group with different target thickness and projectile diameter ratio

    实验tt /dp是否有“尖端”高速摄像结果数值模拟结果
    70.435
    80.514
    90.862
    下载: 导出CSV
  • [1] 邸德宁, 陈小伟, 文肯, 等. 超高速碰撞产生的碎片云研究进展 [J]. 兵工学报, 2018, 39(10): 2016–2047. DOI: 10.3969/j.issn.1000-1093.2018.10.018.

    DI D N, CHEN X W, WEN K, et al. A review on the study of debris cloud produced by normal hypervelocity impact upon a thin plate [J]. Acta Armamentarii, 2018, 39(10): 2016–2047. DOI: 10.3969/j.issn.1000-1093.2018.10.018.
    [2] 宋光明, 武强, 李明, 等. 超高速撞击下空间碎片形状效应研究进展 [J]. 装备环境工程, 2020, 17(3): 45–52. DOI: 10.7643/issn.1672-9242.2020.03.008.

    SONG G M, WU Q, LI M, et al. Research progress in shape effect of space debris under hypervelocity impact [J]. Equipment Environmental Engineering, 2020, 17(3): 45–52. DOI: 10.7643/issn.1672-9242.2020.03.008.
    [3] 焦文俊, 陈小伟. 长杆高速侵彻问题研究进展 [J]. 力学进展, 2019, 49(1): 201904. DOI: 10.6052/1000-0992-17-021.

    JIAO W J, CHEN X W. Review on long-rod penetration at hypervelocity [J]. Advances in Mechanics, 2019, 49(1): 201904. DOI: 10.6052/1000-0992-17-021.
    [4] RANDLES P W, LIBERSKY L D. Smoothed particle hydrodynamics: some recent improvements and applications [J]. Computer Methods in Applied Mechanics and Engineering, 1996, 139(1): 375–408. DOI: 10.1016/S0045-7825(96)01090-0.
    [5] 周旭, 张雄. 物质点法数值仿真(软件)系统及应用 [M]. 北京: 国防工业出版社, 2015: 6–9.

    ZHOU X, ZHANG X. Material point method simulation system [M]. Beijing: National Defense Industry Press, 2015: 6–9.
    [6] LI B, KIDANE A, RAVICHANDRAN G, et al. Verification and validation of the optimal transportation meshfree (OTM) simulation of terminal ballistics [J]. International Journal of Impact Engineering, 2012, 42: 25–36. DOI: 10.1016/j.ijimpeng.2011.11.003.
    [7] 经福谦. 实验物态方程导引 [M]. 北京: 科学出版社, 1999: 25–29.
    [8] 唐蜜. 基于欧拉方法的超高速撞击程序研制及碎片云相分布数值模拟 [D]. 四川绵阳: 中国工程物理研究院, 2015: 45–51.

    TANG M. Development of hypervelocity impact codes based on Euler method and numerical study of the phase distribution in debris cloud [D]. Mianyang, Sichuan: China Academy of Engineering Physics, 2015: 45–51.
    [9] ZHANG X T, JIA G H, HUANG H. A fast numerical approach for Whipple shield ballistic limit analysis [J]. Acta Astronautica, 2014, 93: 112–120. DOI: 10.1016/j.actaastro.2013.06.014.
    [10] CHI R Q, PANG B J, GUAN G S, et al. Analysis of debris clouds produced by impact of aluminum spheres with aluminum sheets [J]. International Journal of Impact Engineering, 2008, 35(12): 1465–1472. DOI: 10.1016/j.ijimpeng.2008.07.009.
    [11] BEISSEL S R, GERLACH C A, JOHNSON G R. A quantitative analysis of computed hypervelocity debris clouds [J]. International Journal of Impact Engineering, 2008, 35(12): 1410–1418. DOI: 10.1016/j.ijimpeng.2008.07.059.
    [12] COLLÉ A, LIMIDO J, VILA J P. An accurate SPH scheme for dynamic fragmentation modelling [J]. The European Physical Journal Conferences, 2018, 183: 01030. DOI: 10.1051/epjconf/201818301030.
    [13] 邸德宁, 陈小伟. 碎片云SPH方法数值模拟中的材料失效模型 [J]. 爆炸与冲击, 2018, 38(5): 948–956. DOI: 10.11883/bzycj-2017-0328.

    DI D N, CHEN X W. Material failure models in SPH simulation of debris cloud [J]. Explosion and Shock Waves, 2018, 38(5): 948–956. DOI: 10.11883/bzycj-2017-0328.
    [14] POVARNITSYN M E, KHISHCHENKO K V, LEVASHOV P R. Simulation of shock-induced fragmentation and vaporization in metals [J]. International Journal of Impact Engineering, 2008, 35(12): 1723–1727. DOI: 10.1016/j.ijimpeng.2008.07.011.
    [15] EFTIS J, CARRASCO C, OSEGUEDA R A. A constitutive-microdamage model to simulate hypervelocity projectile-target impact, material damage and fracture [J]. International Journal of Plasticity, 2003, 19(9): 1321–1354. DOI: 10.1016/S0749-6419(02)00036-0.
    [16] FENG J P, JING F Q, ZHANG G R. Dynamic ductile fragmentation and the damage function model [J]. Journal of Applied Physics, 1997, 81(6): 2575–2578. DOI: 10.1063/1.363921.
    [17] SEISSON G, HÉBERT D, BERTRON I, et al. Dynamic cratering of graphite: experimental results and simulations [J]. International Journal of Impact Engineering, 2014, 63: 18–28. DOI: 10.1016/j.ijimpeng.2013.08.001.
    [18] CORBETT B M. Numerical simulations of target hole diameters for hypervelocity impacts into elevated and room temperature bumpers [J]. International Journal of Impact Engineering, 2006, 33(1): 431–440. DOI: 10.1016/j.ijimpeng.2006.09.086.
    [19] ZHANG X T, LI X G, LIU T, et al. Element fracture technique for hypervelocity impact simulation [J]. Advances in Space Research, 2015, 55(9): 2293–2304. DOI: 10.1016/j.asr.2015.01.040.
    [20] 曹祥, 汤佳妮, 王珠, 等. 损伤演化对韧性金属碎裂过程的影响 [J]. 爆炸与冲击, 2020, 40(1): 013102. DOI: 10.11883/bzycj-2019-0041.

    CAO X, TANG J N, WANG Z, et al. Effect of damage evolution on the fragmentation process of ductile metals [J]. Explosion and Shock Waves, 2020, 40(1): 013102. DOI: 10.11883/bzycj-2019-0041.
    [21] ROYCE E B. GRAY, a three-phase equation of state for metals: UCRL-51121 [R]. Livermore: Lawrence Livermore National Laboratory, California University, 1971.
    [22] 兰彬, 文鹤鸣. 钨合金长杆弹侵彻半无限钢靶的数值模拟及分析 [J]. 高压物理学报, 2008, 22(3): 245–252. DOI: 10.11858/gywlxb.2008.03.004.

    LAN B, WEN H M. Numerical simulation and analysis of the penetration of tungsten-alloy long rod into semi-infinite armor steel targets [J]. Chinese Journal of High Pressure Physics, 2008, 22(3): 245–252. DOI: 10.11858/gywlxb.2008.03.004.
    [23] GMX-6 Group. Selected Hugoniots: LA-4167-MS [R]. Los Alamos, NM: Los Alamos Scientific Laboratory, 1969.
    [24] 于文静. 导管架海洋平台钢结构在爆炸和火灾作用下的力学性能研究 [D]. 上海: 上海交通大学, 2012: 36–37.

    YU W J. Study on mechanical properties of steel jacket offshore platform in blast and fire [D]. Shanghai: Shanghai Jiaotong University, 2012: 36–37.
    [25] 钱伟长. 穿甲力学 [M]. 北京: 国防工业出版社, 1984: 336–340.
    [26] LIANG S C, LI Y, CHEN H, et al. Research on the technique of identifying debris and obtaining characteristic parameters of large-scale 3D point set [J]. International Journal of Impact Engineering, 2013, 56: 27–31. DOI: 10.1016/j.ijimpeng.2012.07.004.
  • 期刊类型引用(4)

    1. 宋俊柏,刘振皓,吴振强,刘武刚,王龙,邢睿思. 隔热瓦组件低速撞击损伤特性及仿真研究. 航天器环境工程. 2025(01): 39-45 . 百度学术
    2. 刘慕皓,张先锋,谈梦婷,包阔,韩国庆,李逸,孙伟境. 考虑微结构特征的陶瓷材料含损伤本构模型. 爆炸与冲击. 2024(01): 62-78 . 本站查看
    3. 周若璞,曾治鑫,张雄. 超高速碰撞下相变效应的交错网格物质点法研究. 计算力学学报. 2024(01): 81-90 . 百度学术
    4. 曹进,陈春林,马坤,高鹏飞,田洪畅,冯娜,钱秉文. 球形含能结构材料弹体超高速撞击多层薄钢靶的毁伤特性. 含能材料. 2023(08): 786-796 . 百度学术

    其他类型引用(2)

  • 加载中
图(9) / 表(6)
计量
  • 文章访问数:  484
  • HTML全文浏览量:  238
  • PDF下载量:  176
  • 被引次数: 6
出版历程
  • 收稿日期:  2021-07-30
  • 修回日期:  2022-01-04
  • 网络出版日期:  2022-04-07
  • 刊出日期:  2022-09-29

目录

/

返回文章
返回