稀土硼化物LaB<sub>6</sub>晶体材料的弹性性质和力学性能

张飞鹏; 施加利; 张静文; 包黎红; 秦国强; 张光磊; 杨新宇; 张久兴

doi:10.11858/gywlxb.20180668

稀土硼化物LaB₆晶体材料的弹性性质和力学性能

doi: 10.11858/gywlxb.20180668

1.
河南城建学院数理学院建筑光伏一体化河南省工程实验室，河南平顶山　467036
2.
石家庄铁道大学材料科学与工程学院，河北石家庄　050043
3.
合肥工业大学材料科学与工程学院新型功能材料与器件安徽省重点实验室，安徽合肥　230009
4.
内蒙古师范大学物理与电子信息学院，内蒙古呼和浩特　010022

基金项目: 国家自然科学基金（51572066, 51302129, 51662034）；河南省科技计划项目（132300410071）；河南省自然科学基金（162300410007）

详细信息

作者简介:
张飞鹏（1980－），男，博士，副教授，主要从事功能材料研究. E-mail: zhfp163@163.com

中图分类号: O521.2; TN34
计量
- 文章访问数: 9683
- HTML全文浏览量: 4466
- PDF下载量: 38
出版历程
- 收稿日期: 2018-10-19
- 修回日期: 2018-11-14

Elastic and Mechanical Properties of Rare Earth Boride LaB₆ Crystalline Material

1.
Henan Provincial Engineering Laboratory of Building-Photovoltaics, Institute of Sciences, Henan University of Urban Construction, Pingdingshan 467036, China
2.
School of Materials Sciences and Engineering, Shijiazhuang Tiedao University, Shijiazhuang 050043, China
3.
Anhui Provincial Key Laboratory of Advanced Functional Materials and Devices, College of Materials Science and Engineering, Hefei University of Technology, Hefei 230009, China
4.
College of Physics and Electronic Information, Inner Mongolia Normal University, Hohhot 010022, China

摘要

摘要: 利用密度泛函理论和Birch-Murnaghan物态方程，系统分析了LaB₆晶体材料的弹性常数参数、体弹性模量、剪切弹性模量及其他力学性能。结果表明：LaB₆晶体具有较大的弹性常数参数C₁₁，说明在此主轴应力方向上具有较大的弹性常数；同时它还具有较大的体弹性模量，并且体弹性模量具有各向同性，剪切弹性模量具有各向异性；LaB₆晶体的杨氏模量为227.85 GPa，泊松比为0.26，体剪弹性模量比值达到1.44，表明其脆性较强，不易发生弹性形变；LaB₆晶体的硬度达到11.56 GPa，平均弹性波速达4.87 km/s。LaB₆的带隙宽度为0.20 eV，呈金属性，内部电子具有较强的局域性，La和B之间具有较强的共价键成分。
- LaB₆ /
- 弹性性质 /
- 力学性能 /
- 弹性常数
Abstract: The elastic constants parameters, bulk modulus, shear modulus and mechanical properties of the LaB₆ crystalline material have been systematically studied with density functional theory and the Birch-Murnaghan formation. The results show that the LaB₆ has larger C₁₁, indicating its larger elastic constant along this direction. The LaB₆ also has large bulk modulus, and the bulk modulus is isotropic, whereas the shear modulus is moderately anisotropic. The LaB₆ crystalline material has large Young’s modulus of 227.85 GPa, so it is hard for elastic deformation. The Poisson’s ratio is 0.26 and the B/G value is 1.44, these values indicate the moderate brittleness of this material. The hardness of LaB₆ is 11.56 GPa and the mean elastic velocity is 4.87 km/s. The LaB₆ has narrow band gap of 0.20 eV, showing that it is metallic type material. The electrons of LaB₆ have strong localization interaction, and there is covalent bond between La and B.
- LaB₆ /
- elastic properties /
- mechanical properties /
- elastic constant

HTML全文

冲击载荷下延性金属的层裂破坏是一个复杂的损伤演化过程，目前被广泛认同的延性金属损伤发展模式是孔洞首先在材料内部成核，随后孔洞长大，当孔洞长大到一定程度时，孔洞之间开始聚集，最终在材料内部形成宏观断裂面^[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. 实验现象

在延性金属的初始层裂实验中，通过软回收靶板，观察到孔洞之间的聚集现象。图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

下载: 全尺寸图片幻灯片

2. 有限元模型

在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

下载: 全尺寸图片幻灯片

在压缩过程中，孔洞具有与周围靶板相同的材料性能，孔洞的成核应力 ${p_{\rm c}}$ ^[19]可以估算如下

${p\rm_c} = \frac{2}{3}\left( {{\sigma _0}+\frac{{{k\rm_y}}}{{\sqrt {{d\rm_G}} }}} \right)\left[ {1 - \ln \frac{3}{2}\left( {\frac{{{\sigma _0}}}{E}+\frac{{{k\rm_y}}}{{E\sqrt {{d\rm_G}} }}} \right)} \right]$

(1)

式中： $E$ 为杨氏模量； ${d\rm_G}$ 为晶粒度，这里设为78 μm^[16]； ${\sigma _0}$ 和 ${k_{\rm{y}}}$ 为模型参数，对于高纯铜， ${\sigma _0}$ =200 MPa， ${k\rm_y}$ =0.14 MPa·m^1/2。

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

${\sigma _{ij}} = - p+{S_{ij}}$

(2)

式中： ${\sigma _{ij}}$ 为应力张量分量， $p$ 为静水压力， ${S_{ij}}$ 为t时刻下的偏应力张量分量。若von Mises屈服函数满足

$\phi = \frac{1}{2}{S_{ij}}{S_{ij}} - \frac{{\sigma _{\rm{y}}^2}}{2} \leqslant 0$

(3)

则偏应力张量为

$S_{ij}^{t+1} = 2G\left( {{\varepsilon _{ij}} - {\varepsilon _v}} \right)$

(4)

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

$\phi = \frac{1}{2}{S_{ij}}{S_{ij}} - \frac{{\sigma _{\rm{y}}^2}}{2} > 0$

(5)

则偏应力张量为

$S_{ij}^{t+1} = \frac{{{\sigma _{\rm{y}}}}}{{{S^*}}}{S_{ij}},\;\;\;\;\;\;{S^*} = {\left( {\frac{2}{3}{S_{ij}}{S_{ij}}} \right)^{1/2}}$

(6)

式中： $S_{ij}^{t+1}$ 为t+1时刻下的偏应力张量分量， $G$ 为剪切模量， ${\varepsilon _v}$ 为体积应变， ${\sigma \rm_y}$ 为屈服应力。 ${\sigma \rm_y}$ 满足Johnson-Cook本构，即

${\sigma }_{\text{y}}=\left(A+B{\varepsilon }_{\rm p}^{n}\right)\left(1+C\mathrm{ln}\,{\dot{\varepsilon }}_{\rm p}^{*}\right)\left(1-{T}^{*}{}^{m}\right)$

(7)

式中： ${\varepsilon \rm_p}$ 为等效塑性应变； ${\dot{\varepsilon }}_{\rm p}^{*}\text=({\dot{\varepsilon }}_{\rm p}/{\dot{\varepsilon }}_{0})$ ， ${\dot \varepsilon \rm_p}$ 为等效塑性应变率， ${\dot \varepsilon _0}$ 为参考应变率； $A$ 为参考应变率和参考温度下的初始屈服应力； $B$ 为材料应变硬化模量； $C$ 为材料应变率强化参数； $n$ 为硬化指数； $m$ 为材料热软化指数； ${T}^{*}=(T-{T}_{\text{room }})/({T}_{\text{melt }}-{T}_{\text{room }})$ ，T_room为室温，T_melt为熔点。采用Johnson-Cook损伤模型模拟材料的损伤和失效，即

${\varepsilon _{\rm f}} = \left[ {{d_1}+{d_2}\exp \left( {{d_3}{\sigma ^*}} \right)} \right]\left( {1+{d_4}\ln {{\dot \varepsilon }\rm_p}^*} \right)\left( {1+{d_5}{T^*}} \right)$

(8)

式中： ${\sigma ^*} = - p/{\sigma _{\rm eq}}$ ，其中 ${\sigma _{\rm eq}}$ 为Mises等效应力； ${d_1}$ 、 ${d_2}$ 、 ${d_3}$ 、 ${d_4}$ 、 ${d_5}$ 为材料损伤系数。

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

$p = \frac{{K\eta }}{{{{(1 - s\eta )}^2}}}\left( {1 - \frac{{{\varGamma _0}\eta }}{2}} \right)+{\varGamma _0}{\rho _0}{e\rm_m}\eta = \frac{\rho }{{{\rho _0}}} - 1$

(9)

式中： $K$ 为体积模量， $\rho$ 为当前的密度， ${\rho _0}$ 为材料初始密度， $s$ 为材料常数, ${\varGamma _0}$ 为材料的Grüneisen系数， ${e\rm_m}$ 为材料的比内能。高纯铜的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/MPa	B/MPa	C	m	n	${\dot \varepsilon}_0$ /s^–1	T_room/K	T_melt/K
8910	90	292	0.025	1.09	0.31	1	300	1356

C₀/(m·s^－1)	s	Γ₀	G/GPa	d₁	d₂	d₃	d₄	d₅
3910	1.51	2	46.6	0.54	4.89	3.03	0.014	1.12

下载: 导出CSV

| 显示表格

3. 计算结果与分析

3.1 双孔的聚集行为

3.1.1 初始韧带距离的影响

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

以孔洞间初始韧带距离d_IL0=20 μm的两孔为例，图3给出了孔洞上A点和B点的位移时程曲线以及孔洞直径增长速率时程曲线。从A、B两点的位移时程曲线可以发现：t=1.57 μs左右时，A、B点的位移确实存在一个偏差点，此时孔洞开始发生聚集。从孔洞直径增长速率曲线可以发现：孔洞直径增长速率起初缓慢增长，随后进入近似线性增长阶段，此阶段的直径增长加速度K_g约为2.995 Gm/s²；当A、B两点位移出现偏差时，对应的孔洞直径增长速率出现非线性变化，此阶段的直径增长加速度K_c约为1.717 Gm/s²。孔洞间塑性应变场的交叠会影响孔洞直径的增长，孔洞在长大和聚集过程中的直径增长加速度存在明显变化，以此可判定孔洞的聚集情况；在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 不同时刻下d_IL0=20 μm的孔洞聚集演化云图

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

下载: 全尺寸图片幻灯片

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

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

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

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

下载: 全尺寸图片幻灯片

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

Figure 6. Time history curves of evolution of ILD

下载: 全尺寸图片幻灯片

3.1.2 孔洞位置分布的影响

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

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

图9(a)给出了不同孔径比条件下各孔洞直径增长速率时程曲线。可以发现：孔洞直径增长速率均表现为先缓慢增长，随后进入近似线性增长阶段，接着出现非线性增长阶段，最后孔洞直径减速增长。当孔径比为1/5时，小孔和大孔的直径增长速率差异明显，随着孔径比的逐渐增大，这种差异越来越不明显，说明大孔抑制小孔的生长，但是这种抑制效果随孔径比的增大而变弱。根据上述孔洞发生聚集的判据发现：当孔径比为1/5时，小孔出现聚集的时间约为1.57 μs，直径增长加速度约为2.228 Gm/s²，大孔出现聚集的时间约为1.59 μs，直径增长加速度约为0.869 Gm/s²；当孔径比为3/5时，小孔出现聚集的时间约为1.58 μs，直径增长加速度约为1.166 Gm/s²，大孔出现聚集的时间约为1.59 μs，直径增长加速度约为1.141 Gm/s²；当孔径比为7/5时，大孔和小孔出现聚集的时间均约为1.58 μs，直径增长加速度约为1.441 Gm/s²。统计孔径比为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.2 三孔聚集机制

基于上述两孔聚集行为的模拟和分析，采用三孔几何模型探讨孔洞聚集时的竞争机制。当两孔发生聚集时，孔洞只有一侧与另一孔洞相互作用；当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/s²；当孔洞2与孔洞3之间的韧带距离增加至50 μm时，孔洞1与孔洞2之间聚集的起始时间约为1.57 μs，直径增长加速度约为2.118 Gm/s²，孔洞3向孔洞2聚集的起始时间稍有延迟，约为1.58 μs，直径增长加速度约为1.168 Gm/s²；当孔洞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时，损伤区域的拉伸应力达到孔洞成核应力（阈值），孔洞成核。当d_IL0=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的直径d₃分别设置为5、10、15 μm，孔洞间初始韧带距离设置为50 μm，飞片的撞击速度设置为100 m/s。

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

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

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

下载: 全尺寸图片幻灯片

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

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

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

下载: 全尺寸图片幻灯片

4. 结　论

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

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

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

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

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

图 1 LaB₆的立方结构示意图

Figure 1. Schematic cubic structure of LaB₆

下载: 全尺寸图片幻灯片

图 2 LaB₆晶体材料的能带结构和分态密度

Figure 2. Band structure and partial density of states for the LaB₆ crystalline material

下载: 全尺寸图片幻灯片

表 1 LaB₆晶体材料的晶格参数

Table 1. Structural parameters of LaB₆crystalline material

Method a/nm b/nm c/nm α/(°) β/(°) γ/(°)

Experiment 0.415 49 0.415 49 0.415 49 90 90 90
Calculation 0.420 205 0.420 205 0.420 205 90 90 90

下载: 导出CSV

表 2 LaB₆晶体材料的体弹性模量和剪切弹性模量

Table 2. Bulk modulus and shear modulus of LaB₆ crystalline material

B_V/GPa B_R/GPa B_H/GPa G_V/GPa G_R/GPa G_H/GPa

160.55 160.55 160.55 111.49 68.85 90.17

下载: 导出CSV

表 3 LaB₆晶体材料的力学性能参数和弹性波速

Table 3. Mechanical parameters and elastic velocities of LaB₆ crystalline material

E/GPa γ λ A_B A_G H/GPa v_l/(km·s^–1) v_t/(km·s^–1) v_m/(km·s^–1)

227.85 0.26 1.44 0 0.24 11.56 7.72 4.38 4.87

下载: 导出CSV

表 4 LaB₆晶体材料的电荷转移

Table 4. Charge distributions of LaB₆crystalline material

Atom Charge distribution/e
s orbital p orbital d orbital Total charge

B 0.88 2.54 0.00 –0.42
La 1.50 5.46 1.51 2.53

下载: 导出CSV

参考文献(14)

[1]	LAFFERTY J M. Boride cathodes [J]. Journal of Applied Physics, 1951, 22(3): 299–309. doi: 10.1063/1.1699946
[2]	HOSSAIN F M, RILEY D P, MURCH G E. Ab initio calculations of the electronic structure and bonding characteristics of LaB₆ [J]. Physical Review B, 2005, 72(23): 235101. doi: 10.1103/PhysRevB.72.235101
[3]	CHEN C H, AIZAWA T, IYI N, et al. Structural refinement and thermal expansion of hexaborides [J]. Journal of Alloys and Compounds, 2004, 366(1/2): L6–L8.
[4]	张宁, 张玖兴, 包黎红. 悬浮区域熔炼法制备REB₆(LaB₆、CeB₆)单晶体及其表征 [J]. 功能材料, 2012, 43(2): 178–180 doi: 10.3969/j.issn.1007-4252.2012.02.016 ZHANG N, ZHANG J X, BAO L H, et al. Floating zone growth and characterization of single crystal REB₃ (LaB₆, CeB₆) cathode [J]. Chinese Journal of Functional Materials, 2012, 43(2): 178–180 doi: 10.3969/j.issn.1007-4252.2012.02.016
[5]	包黎红, 那仁格日乐, 特古斯, 等. 放电等离子烧结原位合成La_xCe_{1- x}B₆化合物及性能研究 [J]. 物理学报, 2013, 62: 196105 doi: 10.7498/aps.62.196105 BAO L H, NARENGERILE, TEGUS O, et al. Synthesis and properties of La_xCe_1–xB₆ compounds by in-situ spark plasma sintering [J]. Acta Physica Sinica, 2013, 62: 196105 doi: 10.7498/aps.62.196105
[6]	刘洪亮, 张忻, 王杨, 等. 单晶LaB₆阴极材料典型晶面的电子结构和发射性能研究 [J]. 物理学报, 2018, 67: 048101 doi: 10.7498/aps.67.20172187 LIU H L, ZHANG X, WANG Y, et al. Surface electronic structures and emission property of single crystal LaB₆ typical surfaces [J]. Acta Physica Sinica, 2018, 67: 048101 doi: 10.7498/aps.67.20172187
[7]	高瑞兰, 于化顺, 于普涟, 等. LaB₆多晶材料的制备工艺研究 [J]. 山东大学学报(工学版), 2002, 32(6): 593–596 doi: 10.3969/j.issn.1672-3961.2002.06.024 GAO R L, YU H S, YU P L, et al. Preparation of LaB₆ polycrystalline materials [J]. Journal of Shandong University (Engineering Science), 2002, 32(6): 593–596 doi: 10.3969/j.issn.1672-3961.2002.06.024
[8]	NYE J F. Physical properties of crystals: their representation by tensors and matrices [M]. Oxford: Oxford University Press, 1957.
[9]	HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society Section A, 1952, 65(5): 349. doi: 10.1088/0370-1298/65/5/307
[10]	TIAN Y, XU B, ZHAO Z. Microscopic theory of hardness and design of novel superhard crystals [J]. International Journal of Refractory Metals and Hard Materials, 2012, 33: 93–106. doi: 10.1016/j.ijrmhm.2012.02.021
[11]	BORN M, HUANG K. Dynamical theory of crystal lattices [M]. Oxford: Clarendon Press, 1968.
[12]	HAINES J, LEGER J M, BOCQUILLON G. Synthesis and design of superhard materials [J]. Annual Review of Materials Research, 2001, 31(1): 1–23. doi: 10.1146/annurev.matsci.31.1.1
[13]	PANDA K B, CHANDRAN K S R. Determination of elastic constants of titanium diboride (TiB₂) from first principles using FLAPW implementation of the density functional theory [J]. Computational Materials Science, 2006, 35(2): 134–150. doi: 10.1016/j.commatsci.2005.03.012
[14]	ANDERSON O L. A simplified method for calculating the Debye temperature from elastic constants [J]. Journal of Physics and Chemistry of Solids, 1963, 24(7): 909–917. doi: 10.1016/0022-3697(63)90067-2

施引文献

资源附件(0)

访问统计

点击查看大图

图(2) / 表(4)

计量

文章访问数: 9683
HTML全文浏览量: 4466
PDF下载量: 38

1. 实验现象
2. 有限元模型
3. 计算结果与分析
3.1 双孔的聚集行为
3.2 三孔聚集机制
4. 结　论

Method	a/nm	b/nm	c/nm	α/(°)	β/(°)	γ/(°)
Experiment	0.415 49	0.415 49	0.415 49	90	90	90
Calculation	0.420 205	0.420 205	0.420 205	90	90	90

B_V/GPa	B_R/GPa	B_H/GPa	G_V/GPa	G_R/GPa	G_H/GPa
160.55	160.55	160.55	111.49	68.85	90.17

E/GPa	γ	λ	A_B	A_G	H/GPa	v_l/(km·s^–1)	v_t/(km·s^–1)	v_m/(km·s^–1)
227.85	0.26	1.44	0	0.24	11.56	7.72	4.38	4.87

Atom	Charge distribution/e
Atom	s orbital	p orbital	d orbital	Total charge
B	0.88	2.54	0.00	–0.42
La	1.50	5.46	1.51	2.53

稀土硼化物LaB6晶体材料的弹性性质和力学性能

doi: 10.11858/gywlxb.20180668

作者简介: 张飞鹏（1980－），男，博士，副教授，主要从事功能材料研究. E-mail: zhfp163@163.com

计量

出版历程