一维应变层裂过程中孔洞聚集行为研究

钟政 蒋招绣 王永刚

钟政, 蒋招绣, 王永刚. 一维应变层裂过程中孔洞聚集行为研究[J]. 高压物理学报, 2023, 37(2): 024201. doi: 10.11858/gywlxb.20220655
引用本文: 钟政, 蒋招绣, 王永刚. 一维应变层裂过程中孔洞聚集行为研究[J]. 高压物理学报, 2023, 37(2): 024201. doi: 10.11858/gywlxb.20220655
ZHONG Zheng, JIANG Zhaoxiu, WANG Yonggang. Coalescence Behavior of Voids during One-Dimensional Strain Spallation[J]. Chinese Journal of High Pressure Physics, 2023, 37(2): 024201. doi: 10.11858/gywlxb.20220655
Citation: ZHONG Zheng, JIANG Zhaoxiu, WANG Yonggang. Coalescence Behavior of Voids during One-Dimensional Strain Spallation[J]. Chinese Journal of High Pressure Physics, 2023, 37(2): 024201. doi: 10.11858/gywlxb.20220655

一维应变层裂过程中孔洞聚集行为研究

doi: 10.11858/gywlxb.20220655
基金项目: 国家自然科学基金(11972202)
详细信息
    作者简介:

    钟 政(1997−),男,硕士研究生,主要从事金属层裂损伤演化研究. E-mail:zhongzhengdd@163.com

    通讯作者:

    王永刚(1976−),男,博士,教授,主要从事冲击动力学研究. E-mail:wangyonggang@nbu.edu.cn

  • 中图分类号: O346.1

Coalescence Behavior of Voids during One-Dimensional Strain Spallation

  • 摘要: 延性金属的层裂行为包含孔洞成核、长大和聚集过程,其中孔洞聚集难以通过实验方法直接观测。采用有限元分析方法,研究了延性金属层裂过程中孔洞之间的聚集行为和竞争机制,讨论了孔洞间初始韧带距离、孔洞直径和孔洞位置分布对孔洞聚集的影响。通过实时统计孔洞长大过程中直径的变化,定量分析孔洞聚集的起始时刻。计算结果显示:当初始韧带距离由20 μm增加至50 μm时,相同孔径孔洞间聚集的起始时间不断延长,聚集时的直径增长加速度由1.717 Gm/s2降至0.602 Gm/s2;具有不同孔径比的孔洞之间发生聚集时,小孔优先向大孔靠拢;45°分布方式的孔洞在长大阶段的直径增长加速度最大,约为3.179 Gm/s2,且最早发生聚集。计算结果显示,在三孔聚集中,初始韧带距离和孔径相同的孔洞之间几乎同时发生聚集并贯通,且孔洞间聚集的起始时间随着初始韧带距离的增加而延长,大孔向附近小孔聚集的起始时刻较迟。计算结果揭示了层裂实验中难以观察到的孔洞长大和聚集的细观物理过程,对了解材料层裂的物理过程和本质具有重要的参考价值。

     

  • 冲击载荷下延性金属的层裂破坏是一个复杂的损伤演化过程,目前被广泛认同的延性金属损伤发展模式是孔洞首先在材料内部成核,随后孔洞长大,当孔洞长大到一定程度时,孔洞之间开始聚集,最终在材料内部形成宏观断裂面[1-2],整个过程和物理机制十分复杂。孔洞聚集是层裂过程中的关键环节,因此,研究孔洞聚集行为对揭示材料层裂的整个物理过程和本质具有重要意义。

    近年来,研究人员针对影响层裂行为的多种因素进行了研究,如外载荷条件[3-4]、材料微观结构[5]等。Curran等[6]综合各种实验结果,总结出两种微孔洞的聚集方式:一种是孔洞之间的直接聚集,由塑性应变场相互作用造成,一般发生在延性金属中;另一种是孔洞之间材料的局部塑性失稳,孔洞之间出现一个很窄的局部化区域,孔洞通过该区域连接,从而发生聚集。同时,也有学者为了描述孔洞成核和长大过程,建立了一些微观层裂模型,如NAG模型[1]、VG模型[7]等。目前,孔洞聚集的物理机制研究相对较少,部分学者[8-9]通过逾渗理论描述微孔团簇之间的相互作用,但这也是基于平均场理论建模的。

    随着计算能力的不断提高,多尺度数值模拟方法越来越多地运用于延性金属中的孔洞聚集研究中。Tu等[10]通过分子动力学方法模拟了单晶铝中两个孔洞的聚集过程,结果表明,随着压缩应变的增大,屈服强度呈现先降低后升高的趋势,压缩阶段孔洞表面产生的位错阻碍了其在拉伸阶段的聚集。Yang等[11]对铝中的孔洞聚集进行了分子动力学模拟,结果表明,孔洞聚集的持续时间随温度的升高而延长。Zhao等[12]对高应变率拉伸载荷下单晶铜中孔洞长大和聚集进行了分子动力学模拟,结果表明,应力三轴度的峰值先随初始韧带距离(initial ligament distance,ILD)的增大而增大,直到达到孔洞聚集的临界点,随后减小,应力三轴度可能是表征孔洞是否发生聚集的一个重要指标。王永刚等[13]采用有限元分析方法,研究了延性金属动态拉伸断裂过程中孔洞的贯通行为,发现孔洞之间开始贯通的临界韧带距离约为0.5倍孔洞的实时直径,并指出临界孔洞间韧带距离判据和临界孔洞体积分数判据在物理本质上是一致的。Hure[14]采用有限元方法,利用有效各向同性屈服应力,提出了一个多孔单晶材料的聚集准则。Holte等[15]通过对应变梯度增强材料基体中嵌入离散空隙的单胞模型进行有限元分析,结果表明,在没有应变梯度的影响下,减小孔洞间韧带距离会增强孔洞的聚集倾向。

    由于实验过程中难以实时观测材料的内部信息,人们对孔洞聚集行为的认识仍不够深刻,数值模拟方法便成为研究层裂过程中孔洞聚集行为的有效途径之一。考虑到分子动力学模拟技术所涉及的尺度相对较小,无法真实还原层裂实验中孔洞聚集至靶板断裂的整个损伤演化过程,为此,我们早期通过有限元数值模拟,以厚度比为1∶2建立了一维应变下的平板撞击模型[16],进行了层裂过程中的损伤演化分析。但是,对于随机孔洞的损伤演化而言,其聚集行为仍很难用定量的方式剖析孔洞之间的相互作用。为定量分析孔洞在层裂过程中的聚集行为,本研究拟建立平板撞击条件下高纯铜层裂损伤演化的有限元模型,重点探讨孔洞间的初始韧带距离、孔洞直径以及孔洞分布位置对孔洞聚集行为的影响规律,为建立合理的孔洞聚集模型奠定基础,进一步揭示延性金属层裂的物理机制。

    在延性金属的初始层裂实验中,通过软回收靶板,观察到孔洞之间的聚集现象。图1(a)、图1(b)、图1(c)是早期实验照片[17],可以观察到相同孔洞(孔径相差不大)、不同孔洞以及三孔洞聚集后的现象;从图1 (d)[18]中可以观察到不同位置的孔洞聚集后的现象(图中均用红圈标出)。基于这些实验现象,本研究将通过有限元方法进行对比分析。

    图  1  软回收靶板中孔洞的光学照片:(a) 相同孔洞,(b) 不同孔洞,(c) 3个孔洞,(d) 不同位置孔洞
    Figure  1.  Optical photo of voids in soft recycling target: (a) same voids, (b) different voids, (c) three voids, (d) voids in different position

    在ABAQUS有限元软件中建立高纯铜飞片和靶板模型,为减少计算时间、节约计算成本,这里建立的飞片和靶板模型均为二维模型,所建模型在厚度方向上与实验保持一致,如图2所示,其中飞片的宽度为2.0 mm,厚度为0.4 mm,靶板的尺寸为4.0 mm×0.4 mm。为了模拟一维应变条件,对飞片和靶板施加轴对称周期边界条件。在对称碰撞条件下,由波传播理论分析可知,当飞片以一定的速度撞击靶板时,接触面会向飞片和靶板分别传入冲击压缩波,冲击压缩波分别在飞片和靶板的自由面上反射形成卸载波,相向的卸载波在靶板中相遇,产生拉伸应力区。当飞片和靶板的厚度比为1∶2时,靶板中心及其附近区域拉伸应力的持续时间最长,损伤主要集中在该区域内,因此在靶板中心位置设置边长为0.4 mm的正方形作为损伤区域。在损伤区域的中心布置微米级的填充孔洞作为初始成核源。

    图  2  高纯铜层裂的几何模型
    Figure  2.  Geometric model of high-purity copper spallation

    在压缩过程中,孔洞具有与周围靶板相同的材料性能,孔洞的成核应力pc[19]可以估算如下

    pc=23(σ0+kydG)[1ln32(σ0E+kyEdG)]
    (1)

    式中:E为杨氏模量;dG为晶粒度,这里设为78 μm[16]σ0ky为模型参数,对于高纯铜,σ0=200 MPa,ky=0.14 MPa·m1/2

    计算中,高应变率加载下高纯铜材料的本构关系采用理想流体弹塑性本构关系,即

    σij=p+Sij
    (2)

    式中:σij为应力张量分量,p为静水压力,Sijt时刻下的偏应力张量分量。若von Mises屈服函数满足

    ϕ=12SijSijσ2y20
    (3)

    则偏应力张量为

    St+1ij=2G(εijεv)
    (4)

    反之,若von Mises屈服函数满足

    ϕ=12SijSijσ2y2>0
    (5)

    则偏应力张量为

    St+1ij=σySSij,S=(23SijSij)1/2
    (6)

    式中:St+1ijt+1时刻下的偏应力张量分量,G为剪切模量,εv为体积应变,σy为屈服应力。σy满足Johnson-Cook本构,即

    σy=(A+Bεnp)(1+Cln˙εp)(1Tm)
    (7)

    式中:εp为等效塑性应变;˙εp=(˙εp/˙ε0)˙εp为等效塑性应变率,˙ε0为参考应变率;A为参考应变率和参考温度下的初始屈服应力;B为材料应变硬化模量;C为材料应变率强化参数;n为硬化指数;m为材料热软化指数;T=(TTroom )/(Tmelt Troom )Troom为室温,Tmelt为熔点。采用Johnson-Cook损伤模型模拟材料的损伤和失效,即

    εf=[d1+d2exp(d3σ)](1+d4ln˙εp)(1+d5T)
    (8)

    式中:σ=p/σeq,其中σeq为Mises等效应力; d1d2d3d4d5为材料损伤系数。

    静水压力p由Mie-Grüneisen状态方程给出

    p=Kη(1sη)2(1Γ0η2)+Γ0ρ0emη=ρρ01
    (9)

    式中:K为体积模量,ρ为当前的密度,ρ0为材料初始密度,s为材料常数, Γ0为材料的Grüneisen系数,em为材料的比内能。高纯铜的J-C本构、损伤及状态方程参数见表1

    表  1  高纯铜的Johnson-Cook本构、损伤及状态方程参数[20-21]
    Table  1.  Johnson-Cook constitutive, damage and equation of state parameters of high-purity copper[20-21]
    ρ/(kg·m-3)A/MPaB/MPaCmn˙ε0/s–1Troom/KTmelt/K
    8910902920.0251.090.3113001356
    C0/(m·s-1)sΓ0G/GPad1d2d3d4d5
    39101.51246.60.544.893.030.0141.12
    下载: 导出CSV 
    | 显示表格
    3.1.1   初始韧带距离的影响

    采用上述有限元模型对两个相同孔洞的聚集行为进行数值模拟。初始成核源孔洞直径设置为5 μm,如图2(a)所示,飞片的撞击速度设置为100 m/s,对这两个孔洞在不同孔洞间初始韧带距离条件下的聚集行为进行分析,在一系列模拟中,孔洞间的初始韧带距离dIL0分别设置为20、30、40和50 μm。王永刚等[13]的研究表明:在孔洞独立长大阶段,孔洞为“球对称”长大,即图1(a)中孔洞上AB两点的位移是一致的;AB两点的位移发生明显改变,表明孔洞之间的塑性应变场交叠影响了B点的位移增长速率,判断此时孔洞开始聚集。其中,B点的位移增长速率减小直接影响孔洞直径的变化,因此通过实时跟踪孔洞长大、聚集过程中孔洞直径的变化可以判断孔洞何时发生聚集。

    以孔洞间初始韧带距离dIL0=20 μm的两孔为例,图3给出了孔洞上A点和B点的位移时程曲线以及孔洞直径增长速率时程曲线。从AB两点的位移时程曲线可以发现:t=1.57 μs左右时,AB点的位移确实存在一个偏差点,此时孔洞开始发生聚集。从孔洞直径增长速率曲线可以发现:孔洞直径增长速率起初缓慢增长,随后进入近似线性增长阶段,此阶段的直径增长加速度Kg约为2.995 Gm/s2;当AB两点位移出现偏差时,对应的孔洞直径增长速率出现非线性变化,此阶段的直径增长加速度Kc约为1.717 Gm/s2。孔洞间塑性应变场的交叠会影响孔洞直径的增长,孔洞在长大和聚集过程中的直径增长加速度存在明显变化,以此可判定孔洞的聚集情况;在t=1.63 μs后,孔洞直径出现了减速增长的变化趋势,对应图4中两孔发生聚集的一侧应变局部化严重,受两孔之间的应力联结作用,最终导致两孔贯通。

    图  3  A和点B的位移时程曲线以及孔洞直径增长速率时程曲线
    Figure  3.  Time history curves of displacement of points A and B and void diameter growth rate
    图  4  不同时刻下dIL0=20 μm的孔洞聚集演化云图
    Figure  4.  Evolution distribution of voids coalescence in dIL0=20 μm at different time

    图4给出了不同时刻dIL0=20 μm的孔洞聚集的等效塑性应变演化云图。可见,孔洞的长大和聚集强烈依赖塑性应变。图4中两孔洞各自形成几乎均匀的球对称塑性应变场;当t=1.57 μs时,孔洞开始发生聚集,孔洞周围的塑性应变场相互交叠,孔洞仍为规则的圆形;当t=1.63 μs时,对应图3中孔洞直径的增长速率呈非线性增长,孔洞之间的塑性应变场严重交叠,阻碍了孔洞的增长,孔洞形状也发生明显变化;当t=1.65 μs时,孔洞之间发生贯穿。孔洞周围产生的塑性应变场相互作用形成新的塑性流动路径,这一路径最终发展为宏观的断裂面。图4显示的结果与图1(a)的实验结果相似。

    图5(a)给出了不同孔洞间初始韧带距离下孔洞直径增长速率时程曲线,根据上述判据判断孔洞之间开始聚集的起始时刻(图5(a)中用虚线标出)。当初始韧带距离dIL0为20、30、40、50 μm时,孔洞长大阶段的直径增长加速度相差不大,分别为2.995、2.770、2.315、2.117 Gm/s2,发生聚集后的加速度均出现骤降,分别为1.717、1.483、0.984、0.602 Gm/s2。随着dIL0的增大,孔洞直径增长速率峰值持续时间也越来越长,出现减速增长趋势的时间点也越来越延后,意味着孔洞贯通时间也延后。图5(b)将图5(a)中各组孔洞开始聚集的起始时间tndIL0=20 μm的孔洞开始聚集的起始时间t1进行归一化,结果显示:随着dIL0的增加,孔洞之间开始聚集的起始时间逐渐延后。图6给出了不同dIL0条件下实时孔洞间的韧带距离dIL随时间的演化曲线,根据孔洞之间开始聚集的起始时间,容易确定不同dIL0条件下孔洞开始聚集的临界韧带距离(dILc),可见,dILcdIL0的增大而增大。

    图  5  不同dIL0下孔洞直径增长速率时程曲线以及归一化聚集时间曲线
    Figure  5.  Time history curve of void diameter growth rate with different dIL0 and curve of normalized coalescence time
    图  6  实时韧带距离演化时程曲线
    Figure  6.  Time history curves of evolution of ILD
    3.1.2   孔洞位置分布的影响

    在层裂实验中,成核孔洞随机分布在材料内部,本研究通过保持孔洞直径不变、改变孔洞位置,探究孔洞位置分布对孔洞聚集的影响。设置孔洞分别为水平(0°)分布、垂直(90°)分布和呈45°角分布,如图2(a)所示,孔洞间初始韧带距离dIL0均设置为30 μm,飞片的撞击速度为100 m/s。统计3种位置分布方式的孔洞直径增长速率,结果如图7所示。从图7中可以看到:0°分布和90°分布的两个孔洞直径增长速率无明显差异,而45°分布的两个孔洞在聚集后期存在差异;在孔洞长大阶段,45°分布的孔洞直径增长加速度最大,约为3.179 Gm/s2,90°分布方式次之,约为2.900 Gm/s2, 0°分布方式最慢,约为2.770 Gm/s2;当孔洞直径增长速率出现非线性变化时,45°分布、90°分布、0°分布的孔洞直径增长加速度分别降至2.112、1.500和1.483 Gm/s2。通过上述孔洞发生聚集的判据,可以发现,45°分布、90°分布、0°分布孔洞出现聚集的时间分别为1.57、1.60和1.58 μs,即在这3种孔洞位置分布中,45°分布方式的孔洞更容易发生聚集。

    图  7  不同位置分布的孔洞直径增长速率时程曲线
    Figure  7.  Time history curves of voids diameter growth rate with different arrangement

    图8为3种几何构型孔洞聚集过程的演化云图。0°分布的两个孔洞左右两侧的塑性应变场几乎呈对称分布,孔洞间的拉伸破坏使孔洞间出现贯通的损伤带,最终导致两孔洞贯通;90°分布的两个孔洞的塑性应变场则几乎呈上下对称,最终在竖直方向上出现孔洞贯通;45°分布的两个孔洞的塑性应变场在垂直于45°方向上几乎呈对称分布,最终两孔洞沿45°方向贯通。图8显示的结果与图1(d)中90°和45°分布孔洞的实验结果相似。

    图  8  不同时刻下不同位置分布的孔洞聚集演化云图
    Figure  8.  Evolution distribution of voids coalescence in different arrangement at different time
    3.1.3   孔洞直径的影响

    层裂实验中,成核孔洞的大小是随机的,为了探究随机大、小两个孔洞之间的聚集行为,模拟了孔径不同的两个孔洞之间的聚集行为,保持其中一个初始成核孔洞的直径d1为5 μm不变,另一个孔洞直径d2分别设置为1、3和7 μm,则孔径比分别为1/5、3/5、7/5,孔洞间的dIL0设置为30 μm,飞片的撞击速度设置为100 m/s。

    图9(a)给出了不同孔径比条件下各孔洞直径增长速率时程曲线。可以发现:孔洞直径增长速率均表现为先缓慢增长,随后进入近似线性增长阶段,接着出现非线性增长阶段,最后孔洞直径减速增长。当孔径比为1/5时,小孔和大孔的直径增长速率差异明显,随着孔径比的逐渐增大,这种差异越来越不明显,说明大孔抑制小孔的生长,但是这种抑制效果随孔径比的增大而变弱。根据上述孔洞发生聚集的判据发现:当孔径比为1/5时,小孔出现聚集的时间约为1.57 μs,直径增长加速度约为2.228 Gm/s2,大孔出现聚集的时间约为1.59 μs,直径增长加速度约为0.869 Gm/s2;当孔径比为3/5时,小孔出现聚集的时间约为1.58 μs,直径增长加速度约为1.166 Gm/s2,大孔出现聚集的时间约为1.59 μs,直径增长加速度约为1.141 Gm/s2;当孔径比为7/5时,大孔和小孔出现聚集的时间均约为1.58 μs,直径增长加速度约为1.441 Gm/s2。统计孔径比为1/5的孔洞在不同时刻下水平长度与竖直长度的比值(x/y,简称横纵比),如图9(b)所示,结果显示:小孔的横纵比变化明显大于大孔的横纵比变化,即小孔更早出现破坏,这一结果在图10中表现为靠近小孔一侧的韧带破坏更加严重。

    图  9  不同孔洞的直径增长速率时程曲线及孔洞横纵比曲线
    Figure  9.  Time history curves of different voids diameter growth rate and curves of voids horizontal to vertical ratio
    图  10  不同时刻不同孔径比的孔洞聚集演化云图
    Figure  10.  Evolution distribution of voids coalescence in different aperture ratios at different time

    图10为不同孔径比条件下孔洞聚集的演化云图。从图10中可以看到:当孔径比为1/5时,随着时间的推移,大孔几乎呈球形变化,小孔出现了严重的破坏变形。在初期阶段,孔洞周围形成几乎呈球形的塑性应变场,但是这两个孔洞的塑性应变场不对称,与相同孔洞之间韧带的均匀破坏不同,不同孔径两孔之间的韧带出现小孔被大孔拖拽,发生非对称破坏的现象。图10显示的结果与图1(b)中大小孔洞聚集的实验结果相似。

    基于上述两孔聚集行为的模拟和分析,采用三孔几何模型探讨孔洞聚集时的竞争机制。当两孔发生聚集时,孔洞只有一侧与另一孔洞相互作用;当3个孔洞聚集时,中间孔洞的左右两侧均会受到来自其他两个孔洞的作用,导致中间孔洞与附近孔洞发生聚集时存在一定的先后顺序和竞争关系,以下将通过改变孔洞间初始韧带距离和孔洞直径来研究三孔聚集时的竞争机制。

    3.2.1   初始韧带距离的影响

    三孔的几何模型相比于上述两孔的几何模型做了一些调整,在损伤区域内设置3个直径相同的孔洞,直径大小为5 μm,如图2(b)所示。通过改变孔洞间初始韧带距离来研究孔洞之间聚集的竞争机制,这里选取的3组孔洞中,孔洞1与孔洞2之间的韧带距离统一设置为30 μm,孔洞2与孔洞3之间的韧带距离分别设置为30、50、70 μm,飞片的撞击速度设置为100 m/s。

    图11给出了不同孔洞间初始韧带距离下各孔洞的直径增长速率曲线。从图11可以看出:当3个孔洞之间的初始韧带距离均为30 μm时,孔洞之间同时发生聚集,且孔洞之间聚集的起始时间约为1.57 μs,直径增长加速度约为2.484 Gm/s2;当孔洞2与孔洞3之间的韧带距离增加至50 μm时,孔洞1与孔洞2之间聚集的起始时间约为1.57 μs,直径增长加速度约为2.118 Gm/s2,孔洞3向孔洞2聚集的起始时间稍有延迟,约为1.58 μs,直径增长加速度约为1.168 Gm/s2;当孔洞2与孔洞3之间的韧带距离增加至70 μm时,也可观察到与50 μm时相似的情况。通过对比发现,孔洞之间聚集的起始时间随着初始韧带距离的增加而延长。孔洞2由于受两侧塑性应变场的影响,其直径增长速率与孔洞1和孔洞3存在明显区别,但是随着初始韧带距离的增加,孔洞2受孔洞3的影响变小,其直径增长速率越来越接近其他两个孔洞。

    图  11  不同ILD下各孔洞直径增长速率时程曲线
    Figure  11.  Time history curves of void diameter growth rate with different ILD

    图12给出了不同时刻、不同初始韧带距离下3个孔洞聚集的演化云图。当t=1.54 μs时,损伤区域的拉伸应力达到孔洞成核应力(阈值),孔洞成核。当dIL0=30 μm时,孔洞周围的塑性应变场几乎呈对称变化,3个孔洞之间同时发生贯穿,其中孔洞1和孔洞3的最终形状十分相似,而孔洞2则受两侧塑性应变场的影响,在水平方向上被拉长,且最终的形状也与孔洞1和孔洞3存在区别;随着孔洞2与孔洞3之间初始韧带距离的增加,可以看到孔洞1与孔洞2之间的塑性应变场优先相互作用,即初始韧带距离越短,孔洞越容易发生聚集。相同初始韧带距离下的孔洞1和孔洞2发生贯通的时间均为1.69 μs,而孔洞 2与孔洞3之间贯通的时间随着初始韧带距离的增加而延迟。

    图  12  不同时刻不同初始韧带距离下三孔聚集演化云图
    Figure  12.  Evolution distribution of three voids coalescence in different ILD at different time
    3.2.2   孔洞直径的影响

    通过改变成核孔洞的大小研究不同孔径对孔洞聚集时竞争机制的影响。对损伤区域内成核孔洞的直径进行调整,孔洞1和孔洞2的直径相等,设置为5 μm,孔洞3的直径d3分别设置为5、10、15 μm,孔洞间初始韧带距离设置为50 μm,飞片的撞击速度设置为100 m/s。

    图13给出了三孔聚集时各孔洞直径增长速率时程曲线,结果显示:在d3为5、10、15 μm的条件下,孔洞2与孔洞1几乎同一时刻发生聚集,分别为1.58、1.57、1.57 μs,直径增长加速度分别约为1.916、2.022和2.887 Gm/s2;当d3=5 μm时,孔洞2与孔洞3也几乎同一时刻发生聚集,随着d3的增大,孔洞3的聚集时间为1.58 μs,稍有推迟,直径增长加速度约为0.695 Gm/s2。通过对比发现,孔洞2优先与孔洞1发生聚集,即孔洞优先与相同直径的孔洞聚集,不同直径孔洞之间发生聚集的时刻稍有延迟。

    图  13  不同d3下孔洞直径增长速率时程曲线
    Figure  13.  Time history curves of void diameter growth rate with different d3

    图14给出了不同时刻、不同d3下3个孔洞的长大、聚集演化云图。当t=1.54 μs时,各孔洞成核;随着时间的推移,当d3=5 μm时,孔洞周围的塑性应变场几乎呈对称变化,并且3个孔洞同时贯通;当d3逐渐增大时,孔洞1与孔洞2的塑性应变场优先交叠,并且优先发生贯通,孔洞1与孔洞3之间的塑性应变低于前者。图14显示的结果与图1(c)的实验结果相似。

    图  14  不同时刻不同d3条件下三孔聚集演化云图
    Figure  14.  Evolution distribution of three voids coalescence with different d3 at different time

    采用ABAQUS动力学有限元程序,对二维模型中双孔和三孔的聚集行为进行了数值模拟,讨论了孔洞间初始韧带距离、孔洞直径和孔洞位置分布的影响,得到了如下4点认识。

    (1) 可以通过孔洞直径增长速率判断孔洞之间何时开始聚集,数值模拟结果显示,该方法具有一定的合理性。

    (2) 相同孔径孔洞之间开始聚集的临界韧带距离随孔洞间初始韧带距离的增加而增大,相同孔径孔洞之间聚集时,孔洞各自均匀长大,不同孔径孔洞之间聚集时,小孔洞优先开始向大孔洞聚集。

    (3) 相同孔洞条件下,在孔洞长大阶段,45°分布方式孔洞的增长加速度相比于0°和90°分布时孔洞的增长加速度更大,且45°分布方式的孔洞之间聚集的起始时间最早。

    (4) 三孔聚集时,相同韧带距离下的相同孔洞之间几乎同时聚集且同时贯通,孔洞间聚集的起始时间随初始韧带距离的增加而延后;孔径不同时,大孔向附近小洞聚集的起始时间较迟,相同孔洞之间优先发生聚集并贯通,不同孔洞之间贯通较慢。

  • 图  软回收靶板中孔洞的光学照片:(a) 相同孔洞,(b) 不同孔洞,(c) 3个孔洞,(d) 不同位置孔洞

    Figure  1.  Optical photo of voids in soft recycling target: (a) same voids, (b) different voids, (c) three voids, (d) voids in different position

    图  高纯铜层裂的几何模型

    Figure  2.  Geometric model of high-purity copper spallation

    图 A和点B的位移时程曲线以及孔洞直径增长速率时程曲线

    Figure  3.  Time history curves of displacement of points A and B and void diameter growth rate

    图  不同时刻下dIL0=20 μm的孔洞聚集演化云图

    Figure  4.  Evolution distribution of voids coalescence in dIL0=20 μm at different time

    图  不同dIL0下孔洞直径增长速率时程曲线以及归一化聚集时间曲线

    Figure  5.  Time history curve of void diameter growth rate with different dIL0 and curve of normalized coalescence time

    图  实时韧带距离演化时程曲线

    Figure  6.  Time history curves of evolution of ILD

    图  不同位置分布的孔洞直径增长速率时程曲线

    Figure  7.  Time history curves of voids diameter growth rate with different arrangement

    图  不同时刻下不同位置分布的孔洞聚集演化云图

    Figure  8.  Evolution distribution of voids coalescence in different arrangement at different time

    图  不同孔洞的直径增长速率时程曲线及孔洞横纵比曲线

    Figure  9.  Time history curves of different voids diameter growth rate and curves of voids horizontal to vertical ratio

    图  10  不同时刻不同孔径比的孔洞聚集演化云图

    Figure  10.  Evolution distribution of voids coalescence in different aperture ratios at different time

    图  11  不同ILD下各孔洞直径增长速率时程曲线

    Figure  11.  Time history curves of void diameter growth rate with different ILD

    图  12  不同时刻不同初始韧带距离下三孔聚集演化云图

    Figure  12.  Evolution distribution of three voids coalescence in different ILD at different time

    图  13  不同d3下孔洞直径增长速率时程曲线

    Figure  13.  Time history curves of void diameter growth rate with different d3

    图  14  不同时刻不同d3条件下三孔聚集演化云图

    Figure  14.  Evolution distribution of three voids coalescence with different d3 at different time

    表  1  高纯铜的Johnson-Cook本构、损伤及状态方程参数[20-21]

    Table  1.   Johnson-Cook constitutive, damage and equation of state parameters of high-purity copper[20-21]

    ρ/(kg·m-3)A/MPaB/MPaCmn˙ε0/s–1Troom/KTmelt/K
    8910902920.0251.090.3113001356
    C0/(m·s-1)sΓ0G/GPad1d2d3d4d5
    39101.51246.60.544.893.030.0141.12
    下载: 导出CSV
  • [1] ANTOUN T, CURRAN D R, SEAMAN L, et al. Spall fracture [M]. Berlin: Springer, 2003.
    [2] MALLICK D D, PARKER J, WILKERSON J W, et al. Estimating void nucleation statistics in laser-driven spall [J]. Journal of Dynamic Behavior of Materials, 2020, 6(3): 268–277. doi: 10.1007/s40870-020-00248-6
    [3] JONES D R, FENSIN S J, MARTINEZ D T, et al. Effect of peak stress and tensile strain-rate on spall in tantalum [J]. Journal of Applied Physics, 2018, 124(8): 085901. doi: 10.1063/1.5045045
    [4] XIE P C, WANG Y G, SHI T Y, et al. Damage evolution and spall failure in copper under complex shockwave loading conditions [J]. Journal of Applied Physics, 2020, 128(5): 055111. doi: 10.1063/5.0009521
    [5] 陈伟, 谢普初, 刘东升, 等. 晶粒尺寸对高纯铝板材层裂特性的影响 [J]. 爆炸与冲击, 2021, 41(4): 043102.

    CHEN W, XIE P C, LIU D S, et al. Effect of grain size on the spall behaviors of high-purity aluminum plates [J]. Explosion and Shock Waves, 2021, 41(4): 043102.
    [6] CURRAN D R, SEAMAN L, SHOCKEY D A. Dynamic failure of solids [J]. Physics Reports, 1987, 147(5/6): 253−388.
    [7] JOHNSON J N. Dynamic fracture and spallation in ductile solids [J]. Journal of Applied Physics, 1981, 52(4): 2812–2825. doi: 10.1063/1.329011
    [8] 王永刚. 延性金属动态拉伸断裂及其临界损伤度研究 [D]. 绵阳: 中国工程物理研究院, 2006: 39−47.

    WANG Y G. Study on dynamic tensile fracture and critical damage degree of ductile metals [D]. Mianyang: China Academy of Engineering Physics, 2006: 39−47.
    [9] ZUREY A K, THISSELL W R, JOHNSON J N, et al. Micromechanics of spall and damage in tantalum [J]. Journal of Materials Processing Technology, 1996, 60(1): 261−267.
    [10] TU R C, WEI N, PEI Y M, et al. The effect of compression on the void coalescence under strong dynamic loading [J]. Advances in Materials Science and Engineering, 2022: 1–11.
    [11] YANG X, ZENG X, WANG J, et al. Atomic-scale modeling of the void nucleation, growth, and coalescence in Al at high strain rates [J]. Mechanics of Materials, 2019, 135(1): 98–113.
    [12] ZHAO L Y, LIU Y. Investigation on void growth and coalescence in single crystal copper under high-strain-rate tensile loading by atomistic simulation [J]. Mechanics of Materials, 2020, 151(1): 103615.
    [13] 王永刚, 徐东明. 高应变率下延性金属中微孔洞贯通行为的数值分析 [J]. 兵工学报, 2012, 33(9): 1095–1100.

    WANG Y G, XU D M. Simulation of void coalescence in ductile metals under high strain rate [J]. Acta Armamentarii, 2012, 33(9): 1095–1100.
    [14] HURE J. A coalescence criterion for porous single crystals [J]. Journal of the Mechanics and Physics of Solids, 2019, 124(1): 505–525.
    [15] HOLTE I, SRIVASTAVA A, MARTINEZ-PAEDA E, et al. Interaction of void spacing and material size effect on inter-void flow localization [J]. American Society of Mechanical Engineers Digital Collection, 2021, 88(2): 021010.
    [16] JIANG Z X, ZHONG Z, XIE P C, et al. Characteristics of the damage evolution and the free surface velocity profile with dynamic tensile spallation [J]. Journal of Applied Physics, 2022, 131(12): 125104. doi: 10.1063/5.0082361
    [17] 陈伟. 晶粒尺寸对高纯铝动态力学行为与层裂特性的影响 [D]. 宁波: 宁波大学, 2020: 53−58.

    CHEN W. Effect of grain size on dynamic mechanical behavior and spallation characteristics of high-purity aluminum [D]. Ningbo: Ningbo University, 2020: 53−58.
    [18] 祁美兰. 高纯铝拉伸型动态破坏的临界行为研究 [D]. 武汉: 武汉理工大学, 2007: 44−46.

    QI M L. Critical behavior in dynamic tensile fracture of high purity aluminum [D]. Wuhan: Wuhan University of Technology, 2007: 44−46.
    [19] WILKERSON J W, RAMESH K T. Unraveling the anomalous grain size dependence of cavitation [J]. Physical Review Letters, 2016, 117(21): 215503. doi: 10.1103/PhysRevLett.117.215503
    [20] 徐金中, 汤文辉. SPH方法在层裂损伤模拟中的应用 [J]. 强度与环境, 2009, 36(1): 1–7. doi: 10.3969/j.issn.1006-3919.2009.01.001

    XU J Z, TANG W H. Applications of SPH method to simulate spalling damage [J]. Structure & Environment Engineering, 2009, 36(1): 1–7. doi: 10.3969/j.issn.1006-3919.2009.01.001
    [21] 汤铁钢, 刘仓理. 高应变率拉伸加载下无氧铜的本构模型 [J]. 爆炸与冲击, 2013, 33(6): 581–586. doi: 10.3969/j.issn.1001-1455.2013.06.004

    TANG T G, LIU C L. On the constitutive model for oxygen-free high-conductivity copper under high strain-rate tension [J]. Explosion and Shock Waves, 2013, 33(6): 581–586. doi: 10.3969/j.issn.1001-1455.2013.06.004
  • 加载中
图(14) / 表(1)
计量
  • 文章访问数:  247
  • HTML全文浏览量:  163
  • PDF下载量:  46
出版历程
  • 收稿日期:  2022-09-15
  • 修回日期:  2022-10-20
  • 网络出版日期:  2023-04-17
  • 刊出日期:  2023-04-05

目录

/

返回文章
返回