α-Fe 中氦泡极限压强的分子动力学模拟*

李翔 尹益辉 张元章

1) (中国工程物理研究院总体工程研究所, 绵阳 621999)

2) (西南科技大学应用技术学院, 绵阳 621010)

为了深入认识α-Fe 中氦泡冲出位错环的微观机制, 有必要研究α-Fe 中氦泡冲出位错环时的极限压强特性.本文建立金属-氦泡的立方体型代表性体积单元模型, 针对8 种不同初始半径的球形氦泡, 以初始氦空位比为变量, 开展分子动力学模拟, 得到了各模型中位错环开始形成时的氦泡极限压强和临界氦空位比.研究结果表明: 对于无量纲半径介于2—10 的氦泡, 冲出位错环时的氦泡极限压强和临界氦空位比均随着氦泡初始半径的增大而非线性减小; 基体中氦泡冲出位错环时的临界氦空位比具有明显的尺寸效应; 初始时刻(0 ps),在经过立方体型模型中心的横截面上, 氦泡周围Fe 原子阵列的剪应力集中和最大剪应力出现在对角线与氦泡边界交点(即45°)处, 并且关于横截面上平行于边的两条对折线对称分布, 剪应力集中区的范围和最大剪应力均随着初始氦空位比的增大而增大; 位错环冲出方向对应最大剪应力方向.本文的研究加深了对金属中氦泡物理特性的认识, 为后续分析氦泡对材料宏观物理和力学性质的影响奠定了有益的基础.

1 引 言

核反应堆中第一壁和储氚压力容器是核工业中典型的涉氚涉氦结构.被引入到金属材料中的氦, 几乎不溶于金属基体材料, 会在金属材料中聚集, 形成纳米级的氦泡, 从而引起氦效应, 如使金属材料的韧性降低、延伸率下降等[1−5], 并最终改变结构的承载能力, 影响结构的安全性.因此, 氚在金属基体材料中的渗透与扩散及其引起的氚和氦浓度的时空变化、氦泡在基体材料中的形核长大及其对材料力学性能的影响成为核材料和核技术领域的研究焦点[1−10].

在氦泡形核与长大过程中, 氦泡压强是氦泡的主要物理特性指标, 也是氦泡影响金属宏观物理与力学性能的重要参数.因此, 氦泡的压强特性受到了越来越多的研究者关注[11−15].由于金属中氦泡尺寸大多为纳米量级, 因此氦泡压强很难通过实验手段直接测量.现行做法是先采用电子能量损失谱方法(electron energy-loss spectroscopy, EELS)测量氦泡中的He 原子密度, 再利用氦泡的状态方程, 计算得到氦泡压强[11,12].因基体材料存在机械稳定极限[2,13], 就内部含有高压氦泡的金属基体而言, 当氦泡施加在金属基体上的作用强度超过基体的机械稳定极限时, 氦泡将会通过冲出位错环机制而长大, 此时的氦泡压强为其极限压强[2].由于氦泡极限压强可达数十GPa 量级, 因此它不仅是金属基体中氦泡长大机制的特征参数, 而且也是研究氦泡引起的金属材料微观损伤机理的重要物理参数.氦泡冲出位错环现象是快速发生的[14], 且无明显的临界特征, 采用现有技术条件还难以精确掌握冲出位错环时的临界点, 因此对于氦泡冲出位错环时的临界He 原子密度, 也难以采用EELS 方法进行准确测量, 相应地, 也就难以准确计算氦泡极限压强.

Trinkaus[13]和Wolfer[15]分别通过解析理论建模和求解, 给出了单个氦泡冲出位错环时的极限压强公式.与Trinkaus 的公式相比, Wolfer 的公式更为简洁, 即氦泡冲出位错环时的极限压强为基体材料剪切模量的0.2 倍.随着计算软硬件技术的发展, 分子动力学(molecular dynamics, MD)方法成为有力的工具, 已用于诸多的模拟研究中.Zhang 等[16]利用微观-细观-宏观3 个尺度的实验信息, 校准MD 模拟的模型系统, 研究了氢对HR-2 钢弹性性能的影响, 深化了对氢脆机理的认识.Deng 等[17]采用MD 方法模拟了α-Fe 中氦团簇在块体和晶界的扩散特性.Gao 等[18]模拟了α-Fe 中氦浓度和辐照温度对He 原子与移位级联相互作用的影响.关于金属中的氦泡压强, Zhang 等[19]采用MD 方法, 对钛中初始氦空位比为2 (He 原子数与初始孔洞空位数的比值nHe/nV= 2)的氦泡进行模拟, 得到了在不同温度时, 4 种不同初始半径氦泡的压强, 并通过数据拟合, 给出了氦泡的状态方程; Hetherly 等[20]通过在MD 模型中先用1 个He 原子替换掉原位置的Fe 原子, 然后采用每2 ps 向体系中加入1 个He 原子的方式, 不断添加He 原子进行模拟, 结果表明: 当He 原子数为99时, 氦泡冲出位错环现象发生, 而冲出位错环后,氦泡的压强下降到局部最小值27.20 GPa; Caro等[21]通过MD 模拟得到了Fe 和FeCr 中氦泡压强随He 原子数增加的变化过程, 并观察到了氦泡冲出位错环现象; Xie 等[22]采用MD 方法, 模拟了钨中初始半径为1.00 nm 的氦泡, 结果表明: 当位错演化成棱形位错环并离开氦泡时, 氦泡压强从40.50 GPa 下降到35.60 GPa, 氦泡半径从1.16 nm增大到1.18 nm.

如前所述, 通过MD 模拟研究, 文献[20−22]观察到了冲出位错环现象和相应氦泡压强的变化,但是并没有关注和讨论氦泡极限压强及其与初始氦空位比和氦泡尺寸的关联特性等问题.同时, 根据文献[23−25], 非热形成时, 氦泡的平均半径为2.50 nm 或更小.目前, 对该尺寸范围内氦泡的压强特性, 包括氦泡极限压强, 还有待进一步研究.本文采用MD 方法构建含氦泡的模拟模型, 运用LAMMPS[26]软件进行模拟, 研究不同初始氦空位比下氦泡的压强特性.由于涉核结构和压力容器以铁基合金为主, 故本文以α-Fe 为基体材料, 构建基体材料内含球形氦泡的立方体型代表性体积单元 模 型(representative volume element, RVE).为了分析球形氦泡的尺寸效应, 分别考虑了初始半径为0.50, 0.81, 1.00, 1.31, 1.50, 2.00, 2.12, 2.50 nm的8 种情况, 这些尺寸跨越了氦泡的形核阶段(氦泡半径 < 1.00 nm)和非理想气体阶段(1.00 nm ≤氦泡半径 < 100.00 nm)[13].MD 模拟结果表明: 对于初始半径介于0.50—2.50 nm 的氦泡, 冲出位错环时的氦泡极限压强和临界氦空位比都随着氦泡初始半径的增大而非线性减小.基于MD 模拟结果, 拟合得到了冲出位错环时的氦泡极限压强和临界氦空位比随氦泡初始半径变化的关系式, 与Wolfer的解析分析结果吻合.本文的研究深化了对金属中氦泡物理特性的认识, 为后续分析氦泡对材料宏观物理和力学性质的影响奠定了有益的基础.

2 模型与计算方法

2.1 分子动力学模型

MD 建模和模拟均采用软件LAMMPS 实施.本文对含氦泡的α-Fe 采用如下方法建模: 首先建立尺寸为45a0× 45a0× 45a0的α-Fe 原子阵列模型, Fe 原子总数为182250 个, 这里a0为其晶格常数,a0= 2.86 Å (1 Å = 0.1 nm); 然后以模型中心为球心, 删除一定半径范围内的Fe 原子, 形成球形孔洞; 最后在孔洞中填入不同数量的He 原子,即可形成氦泡.该建模方法在文献[19,27]中也有运用.建立的金属-氦泡的立方体型RVE 模型如图1 所示.图1 所示模型中, Fe-Fe, He-He 和Fe-He 之间的相互作用势, 分别采用Ackland 等[28]、Aziz 等[29]和Gao 等[30]开发的势函数.RVE 模型的3 个晶向[100]、[010]和[001]分别对应于图1中直角坐标系的x,y和z轴, 且3 个方向均施加周期性边界条件, 模拟温度保持在300 K, 选用NPT系综, 模拟时间步长为1 fs.

图1 含氦泡的α-Fe 计算模型(图中He → He 泡, Fe →Fe 基体)Fig.1.Model of α-Fe with helium bubble (In Fig.He → helium bubble, Fe → Fe matrix).

2.2 氦泡压强计算方法及其有效性

氦泡压强的计算表达式为[31]

式中:n是氦泡中的He 原子个数;σ11(i),σ22(i) 和σ33(i)为第i个He 原子的原子应力分量, 由LAMMPS 中的原子应力[32]模块计算;V是氦泡体积, 使用Voro++软件包[33], 采用Voronoi 算法计算得到.

Mills 等[34]通过实验手段测量了液氦的压强,拟合得到了温度区间为 7 5 K

图2 MD 模拟结果与Mills 等[34]的实验结果对比(300 K)Fig.2.Comparison of MD simulation results with experimental results of Mills et al.[34] (300 K).

2.3 计算与数据处理过程

对含有不同初始半径氦泡的各RVE 模型, 以初始氦空位比为变量, 开展MD 模拟.首先对各RVE 模型弛豫100 ps; 然后继续模拟30 ps, 并在此过程中每间隔0.30 ps, 以0.002 ps 的频率提取一次氦泡压强, 连续提取100 次, 总计提取6000 次;最后将6000 次的算术平均值作为最终的氦泡压强.

模拟结束后利用可视化软件OVITO[35]中的Dislocation Extraction Algorithm (DXA)[36]模块,分析各模型中位错形成情况.

3 结果与分析

3.1 位错环冲出与氦泡极限压强

MD 模拟所得氦泡压强-初始氦空位比曲线如图3 所示.图中氦泡压强峰值点用实心标识符表示, 其他数据点用空心标识符, 红色实线是基于各模型的模拟数据点的拟合曲线.图3 表明: 对于金属基体中8 种不同初始半径的氦泡, 各自的氦泡压强随初始氦空位比的变化规律总体上是相同的, 都具有类抛物线特征.在压强峰值点之后, 存在小幅度的波动现象, 该现象的产生归因于金属基体材料存在机械稳定极限, 当氦泡压强即孔洞内压超过基体材料的机械稳定极限时, 氦泡周围的基体材料就产生局部塑性变形和破坏, 从而造成了这些波动.这在Zhang 等[19]、Hetherly 等[20]和Caro 等[21]的研究中也有所揭示.

图3 氦泡压强随初始氦空位比的变化Fig.3.Changes of pressure of helium bubble with initial helium-to-vacancy ratio.

文献[2]的研究表明, 当氦泡压强超过基体材料的机械稳定极限时, 基体材料将发生局部塑性变形和破坏, 最有可能的情况是在基体中冲出位错环.而一旦冲出位错环, 氦泡体积就会明显长大而压强下降.因此将氦泡冲出位错环时的压强定义为氦泡极限压强.但值得注意的是, 图3 中的氦泡压强是在RVE 模型弛豫后提取的, 而弛豫过程中可能会发生冲出自间隙原子或位错环现象.因此图3中各条曲线上的氦泡压强峰值点是基体材料内未发生冲出位错环现象时的最大氦泡压强, 而并非冲出位错环时的极限压强.根据图3 中各条曲线的峰值点, 可以得到对应的初始氦空位比, 按氦泡初始半径由0.50 nm 依次增大到2.50 nm 的顺序, 它们分 别 是2.60, 2.40, 2.25, 2.08, 2.05, 1.93, 1.92 和

1.8 4.通过对MD 模拟过程中的各RVE 模型进行DXA 分析, 可以得到氦泡在基体中冲出位错环时的初始氦空位比—临界氦空位比, 同样按氦泡初始半径由0.50 nm 依次增大到2.50 nm 的顺序, 所得结果依次是2.95, 2.45, 2.30, 2.10, 2.08, 1.94, 1.93和1.85.对比可见, 氦泡压强峰值点对应的初始氦空位比与冲出位错环时的临界氦空位比相比有一定偏差, 后者比前者大.已有研究表明, 高压氦泡在基体中冲出的位错环由基体材料的自间隙原子组成[21,22,37].对于形核初期或中期的氦泡(如氦泡初始半径为0.50 nm 的氦泡), 由于其外表面积较小, 因此在氦泡压强峰值点所对应的初始氦空位比时, 氦压作用下冲出的自间隙Fe 原子数目较少,还不足以形成位错环; 而要形成位错环则需更大的初始氦空位比, 以便形成更高的初始氦压, 冲出更多的自间隙Fe 原子; 对于初始半径为0.50 nm 的氦泡, 其冲出位错环所需的临界氦空比比其压强峰值点对应的初始氦空位比约大13.46%.Caro 等[21]的研究也得到了类似结果.但对于形核后期的氦泡(如氦泡初始半径为0.81 nm 的氦泡)和非理想气体阶段的氦泡(如氦泡初始半径为1.00—2.50 nm的氦泡), 各氦泡压强峰值点对应的初始氦空位比与冲出位错环时的临界氦空位比应当相同, 前述偏差产生的原因主要在于: 本模拟中建立的氦泡, 其初始氦空位比采用间隔取值, 因而受到取值间隔精度的影响, 且在图3 中给出的氦泡压强是在弛豫过后再采样提取的, 未能反应出冲出位错环时的极限压强和临界氦空位比.初始半径为0.50 nm 的氦泡与其他7 种尺寸的氦泡, 在冲出位错环时所需临界氦空位比的不同, 正是氦泡尺寸效应的具体表现.

对于特定初始半径的氦泡, 初始氦空位比越大, 初始氦压也就越高.因此可知氦泡冲出位错环时的极限压强肯定高于图3 中峰值点的氦泡压强.氦泡冲出位错环时, 氦泡体积长大而压强下降, 因而可以通过MD 模拟提取得到模拟过程中不同时间点的氦泡压强, 特别是氦泡冲出位错环前后的压强值.在本文的模拟条件下, 为了获得各模型冲出位错环时的氦泡极限压强, 可由DXA 分析得到冲出位错环时的临界氦空位比, 提取模型在该临界氦空位比下氦泡压强的最大值, 作为冲出位错环时的极限压强.

3.2 氦泡极限压强、临界氦空位比随氦泡初始半径的变化关系

根据DXA 分析结果得到临界氦空位比, 提取氦泡极限压强, 绘制它们随氦泡初始半径的变化曲线, 如图4 所示.通过拟合数据点, 得到氦泡极限压强Pultimate和临界氦空位比(nHe/nV)cr与的关系式分别为:

式中:是无量纲量,,R为氦泡初始半径,b为Burgers 矢量a0/2〈111〉的模,b= 0.25 nm;G是α-Fe 的剪切模量,G= 80.40 GPa[38].

图4 氦泡极限压强、临界氦空位比随氦泡初始半径的变化Fig.4.Changes of ultimate pressure and critical helium-tovacancy ratio of helium bubble with initial radius of helium bubble.

由图4、(2)式和(3)式可知, 氦泡极限压强和临界氦空位比均随氦泡初始半径的增大而非线性减小, 其中前者与Trinkaus[13]和Wolfer[15]的研究结果一致, 后者与Wolfer[15,39]的研究结果吻合.

3.3 氦泡周围Fe 基体的原子应力分析

在氦泡压强驱动下, 基体中产生了冲出自间隙原子及位错环的现象[21,22,37].因此, 这里对氦泡周围Fe 原子阵列的原子应力进行研究, 以便从应力角度理解冲出位错环的机制.图5 给出了经过氦泡初始半径为1.31 nm 的MD 模型中心横截面上的原子剪应力分布, 图中紫色代表He 原子, 其余为基体Fe 原子.图5(a)—(c)3 幅分图依次对应初始氦空位比小于、等于和大于临界氦空位比2.10 的情况, 每幅分图包括模拟时间为0 ps 和100 ps 时基体Fe 原子的应力分布图.由图5 可知: 在图1所示MD 模型坐标系中, 氦泡压强使其周围Fe 原子阵列的剪应力集中和最大剪应力出现在45°处,并且关于横截面上平行于边的两条对折线对称分布, 这与Xie 等[22]研究钨时的结果一致.当初始氦空位比分别小于、等于和大于临界氦空位比时, 与0 ps 时的氦泡形状相比, 100 ps 时的氦泡形状依次是基本保持球形、近似椭球体和多面体, 并且氦泡的体积也依次长大.这种氦泡形状的变化原因在于当初始氦空位比大于等于临界氦空位比时, 氦泡周围的Fe 原子被冲出点阵位置而成为自间隙原子,而冲出方向是沿特定的晶面和晶向进行, α-Fe 为bcc晶体结构, 其密排面为{110}、{112}和{123},密排方向为〈111〉, 因此密排面上沿密排方向的Fe 原子被冲出而成为自间隙原子, 而氦泡中的He 原子则占据相应空位, 使得100 ps 时氦泡的形状呈多面体或近似椭球体.初始氦空位比小于临界氦空位比时, 氦泡压强不足以将点阵位置的Fe 原子冲出, 或者冲出的数量少, 因而100 ps 时氦泡的形状基本与0 ps 时相同, 仍保持为球形.结合图5(a)—(c)可知, 与0 ps 时的原子剪应力分布相比, 100 ps 时, 氦泡周围的应力集中情况变得缓和, 原因是高压氦泡将应力集中区域的部分Fe 原子冲出, 成为自间隙原子, 氦泡内的空位数增多,氦泡长大, 氦泡压强得以释放, 进而氦泡周围的Fe 原子应力集中得到消散.同时, 由图5(a)—(d)可知, 在0 ps 时, 基体Fe 原子的最大剪应力和剪应力集中区的范围均随着初始氦空位比的增大而增大.

图5 氦泡周围Fe 原子阵列的剪应力分布 (a) nHe/nV = 1.50; (b) nHe/nV = 2.10; (c) nHe/nV = 2.50; (d) 上下限相同时0 ps 处的原子应力分布Fig.5.Shear stress distribution of Fe atom array around helium bubble: (a) nHe/nV = 1.50; (b) nHe/nV = 2.10; (c) nHe/nV = 2.50;(d) atomic stress distribution at 0 ps with the same upper and lower limit.

图6 为氦泡在基体中冲出的位错环在xy平面上的投影, 其中红色箭头表示Burgers 矢量方向,红绿蓝色表示离氦泡表面距离不同的Fe 原子, 显示时隐去了完美晶格状态的Fe 原子.包绕在氦泡周围的Fe 原子由两部分组成, 一部分是经氦泡压强冲出点阵的自间隙Fe 原子, 另一部分则是紧邻氦泡的Fe 原子.由图6 可见, 位错环起、止都在氦泡表面, 位错环的冲出方向与图5 中Fe 原子阵列的最大剪应力出现位置吻合.由图7 进一步分析可知, 剪应力分布沿〈111〉方向最大, 这种剪应力分布与图6(d)中冲出的Burgers 矢量为1/2〈111〉型的位错相吻合.这表明当氦泡压强足够大时, 就会使其周围Fe 原子阵列的应力达到极限, 从而相应的Fe 原子就被冲出点阵位置而成为自间隙原子, 而随着自间隙原子的不断累计增多, 就形成间隙型位错环, 如图6(a)—(c)所示.这与文献[40]针对钨中氦泡冲出位错环现象的分析结果一致.

图6 初始半径为1.31 nm, nHe/nV = 2.10 时氦泡冲出的位错环在xy 平面投影 (a) 88 ps; (b) 92 ps; (c) 100 ps; (d)删除图6(c)中Fe 原子后带缺陷网格的位错图Fig.6.Projection of dislocation loops produced by helium bubbles in xy plane, when initial radius of helium bubble is 1.31 nm, nHe/nV = 2.10: (a) 88 ps; (b) 92 ps; (c) 100 ps; (d)a dislocation picture with a defect mesh after removing the Fe atoms in Fig.6(c).

3.4 高压氦泡冲出位错环时的极限压强对比分析

图7 (110)面原子剪应力分布(氦泡初始半径为1.31 nm,nHe/nV = 2.10)Fig.7.(110) Surface atomic shear stress distribution (initial radius of helium bubble is 1.31 nm, nHe/nV = 2.10).

Trinkaus[13]根据基体材料的机械稳定极限, 给出了氦泡冲出位错环时极限压强的理论估算公式:式中,γ是表面自由能, 单位为Pa·m.对于α-Fe,取表面自由能γ= 1.65 Pa·m[38].Donnelly[41]指出(4)式的适用范围为

Wolfer[15]分别基于Helmholtz 自由能和基体材料弹性理论解, 得到氦泡冲出位错环时的极限压强公式为

Trinkaus、Wolfer 和本文MD 的结果对比, 如图8 所示.由Wolfer[15]的研究可知, (5)式只适用于>10 的情况; 但Wolfer 的研究中给出了310范围内的极限压强数据, 如图8 中蓝色实线所示.Wolfer[15,39]指出对于小氦泡其极限压强甚至可达到0.50G, 接近于发射自间隙原子所需压力,这与Trinkaus 在文献[42]中指出的极限压强应该为0.20G—0.50G的结论相吻合.值得注意的是虽然图8 中Wolfer[15]的数据是以Ni 为例, 基于Helmholtz 自由能得出的, 但其在随后的研究中指出氦泡压强(初始氦空位比)才是影响位错环冲出与否的主要因素, 金属材料间的差异为次要因素[39].

由图8 可见, 对于 2<<10 的氦泡, 氦泡极限压强随着氦泡初始半径的增大而非线性减小.对于 2<3 的氦泡(如氦泡初始半径为0.50 nm的氦泡), MD 模拟值相对Trinkaus 和Wolfer 的结果偏低, 但是按照Trinkaus[13]的划分, 初始半径为0.50 nm的氦泡, 属于形核阶段的氦泡, 从冲出位错环所需临界氦空位比(nHe/nV=2.95)角度来分析, MD 模拟值与Wolfer[39]的研究结果吻合.在3<10的范围内, Trinkaus 的计算值相对偏低, 而本文MD的模拟值与Wolfer 的结果比较接近.并且结合文献[2,15]可知, 在10 时, 氦泡极限压强Pultimate≤ 0.20G, 因此本文通过MD 模拟得到的初始半径为2.50 nm 氦泡的极限压强在合理范围内, 并且所需临界氦空位比也与Wolfer[15,39]的研究结果吻合.

图8 氦泡极限压强与氦泡初始半径关系Fig.8.Ultimate pressure vs.initial radius of helium bubble.

本文采用的氦泡压强计算方法的准确性已经在2.2 节加以检验, 可认为在临界氦空位比下计算得到的氦泡极限压强是合理的.不可否认的是, 本文中8 种不同初始半径的氦泡, 冲出位错环时的临界氦空位比, 受到初始氦空位比取值间隔精度影响, 进而影响氦泡极限压强的确定.由3.1 节的分析可知, 各尺寸的氦泡均是在图3 氦泡压强峰值点之后开始出现冲出位错环现象, 如果进一步细化峰值点之后的初始氦空位比, 可以得到更为精确的临界氦空位比, 进而确定不受取值间隔精度影响的氦泡极限压强.

4 结 论

采用分子动力学方法, 建立金属-氦泡的立方体型代表性体积单元模型, 针对8 种不同初始半径的球形氦泡, 以初始氦空位比为变量, 开展MD 模拟, 获得了各RVE 模型中位错环开始形成时的氦泡极限压强和临界氦空位比.通过对MD 模拟结果的分析, 得到如下结论:

1)基于MD 模拟给出了冲出位错环时的氦泡极限压强和临界氦空位比随氦泡初始半径非线性变化的公式, 可用于估算α-Fe 中 2<10 的氦泡极限压强和临界氦空位比.

2)对于 2<10 的氦泡, 在金属基体中冲出位错环时的氦泡极限压强和临界氦空位比随着氦泡初始半径的增大而非线性减小.

3)基体中氦泡冲出位错环时的临界氦空位比具有明显的尺寸效应, 对于形核后期的氦泡(如氦泡初始半径为0.81 nm 的氦泡)和非理想气体阶段的氦泡(如氦泡初始半径为1.00—2.50 nm 的氦泡), 其冲出位错环时的临界氦空位比与各自氦泡压强峰值点对应的初始氦空位比相同, 而对于形核初期或中期的氦泡, 如初始半径为0.50 nm 的氦泡, 其冲出位错环所需的临界氦空比比其压强峰值点对应的初始氦空位比约大13.46%.

4)初始时刻(0 ps), 在经过RVE 中心的横截面上, 氦泡周围Fe 原子阵列的剪应力集中和最大剪应力出现在对角线与氦泡边界交点(即45°)处,并且关于横截面上平行于边的两条对折线对称分布, 剪应力集中区的范围和最大剪应力均随着初始氦空位比的增大而增大; 位错环冲出方向对应于最大剪应力方向.

本文的研究深化了对金属中氦泡物理特性的认识, 为后续分析氦泡对材料宏观物理和力学性质的影响奠定了有益的基础.