刘 明 侯冬杨 高诚辉
(福州大学机械工程及自动化学院,福州 350116)
断裂韧性是指材料抵抗裂纹扩展的能力[1],通常用临界应力强度因子KIC作为衡量指标,单位为字母K为应力场强度因子,反映的是裂纹尖端区域应力场强弱;字母C 指的是裂纹扩展的临界情况;下标罗马数字I 是指裂纹扩展形式为张开型,脆性材料的裂纹扩展类型即为I 型.研究材料的KIC可以用来评估典型工程结构(例如石化容器和储罐、石油和天然气管道以及汽车、船舶和飞机结构[2])的材料性能以及用来确定脆性材料(如石英玻璃)的微观参数[3].
目前,测量材料KIC的方法主要有:山形切口梁法(CNB)[4-5]、单边预裂梁法(SEPB)[4-5]、表面弯曲裂纹法(SCF)[4]、单边切口梁法(SENB)[4-5]、单边V 形切口梁法(SEVNB)[5]、短V 形切口杆法(SR)[5]、双扭法(DT)[4-5]、双悬臂梁法(DCB)[4-5]、微米划痕法[6-8]、压痕法(如纳米压痕法[9-10]和维氏压痕法[11-12])等.SR,DCB 和SEPB 法的测试试样难生产、成本高,故无法得到广泛使用;SENB,SEVNB 和CNB 法加工试样缺口较为困难;DT 法试件的几何尺寸会对测量值产生影响;SCF 法必须要去除足够深度的表面层来消除残余应力场,才能保证KIC不被高估[4-5];微米划痕法需要考虑压头的磨损以确保测试结果的准确性[13];而压痕法具有制备试样简单、测试效率高、以及综合成本低等优点,已被广泛地应用于陶瓷材料[14-15]、硬质合金[16-17]和玻璃材料[18-19]的硬度、弹性模量和韧性测量.
最初,Lawn 等[11]在1980 年基于弹−塑性断裂力学提出了计算KIC的公式(称为LEM 模型).随后,Anstis 等[12]使用了许多较大KIC范围(0.74~12 MPa·)的参考材料,得出了公式中的无量纲常数χ=0.016 ± 0.004.基于Lawn 等[11]和Anstis 等[12]的研究,研究者们提出了许多确定KIC的公式[20-30],见表1.这些公式主要是考虑两种不同的裂纹类型得到的,一种是基于半椭圆型的巴氏裂纹(Palmqvist crack);另一种是基于半月状的中位裂纹(median crack),而第三种是基于曲线拟合方法得到的,同时适用于两种裂纹系统.
近年来,使用压痕法确定KIC的有效性引起了激烈的争论.文献[31-32]认为,计算KIC的公式在实际应用中有许多的限制或问题,例如裂纹系统(径向、中位、横向等)的多样性引起的复杂性,压头下方应力场的不确定性以及有时会得到不可靠的KIC计算结果.因此,对于陶瓷材料的断裂韧性测试,压痕法不建议被使用.然而,Marshall 等[33]认为压痕法是基于Griffith-Irwin 平衡断裂力学得到的,并且压痕是反映材料机械性能独一无二的特征,因此,压痕法仍然是有效的方法.但是,使用压痕法评估材料的KIC仍然存在许多不足(如对低泊松比的材料并不适用[34]),如何根据不同的材料、不同的压头选择适合的公式,成为当前压痕法在实际应用中需要解决的困难之一.
单晶硅作为一种性能优越的半导体材料,现已被广泛应用于芯片、光学元件和太阳电池[35]等领域的加工制造.碳化硅作为第三代半导体材料,由于其具有更大的禁带宽度,非常适合用于制造5G 射频器件和高电压功率器件[36],裂纹扩展进而引起的脆性断裂是其主要失效形式之一[37-39].在一般工程构件中,这些材料发生脆断前并没有明显的塑性变形,而断裂韧性可以很好地衡量这种特性.本工作基于单晶硅(111)和碳化硅(4H-SiC)维氏微米硬度和玻氏纳米压痕实验KIC测试结果,与划痕实验测试结果进行比对分析后,讨论了使用表1 总结的13 个压痕公式评估半导体材料KIC的适用性,并希望获得有效的压痕实验条件和特定压头下材料的适用公式,以期为压痕法表征半导体材料的断裂韧性提供参考.
表1 利用维氏硬度法计算材料断裂韧性的方程Table 1 Equations for calculating fracture toughness of materials by Vickers hardness method
实验压痕试样涉及两种半导体材料,单晶硅(111)和碳化硅(4H-SiC,0001 面),试样尺寸为块体,满足平面应变条件.实验之前对试样表面进行多级水磨,并通过机械抛光去除表面划痕直到达到镜面效果,然后在无水乙醇中对试样进行超声波清洗,以清除试样表面残留的研磨颗粒等杂质,材料的弹性模量和泊松比见表2.
表2 样品参数和维氏压痕尺寸Table 2 Parameters of samples and Vickers indentation
1.2.1 压痕实验
维氏压痕实验使用数显硬度计(MHVKD-100,上海钜晶仪器制造有限公司)将维氏压头按确定的载荷垂直压入样品表面,保载15 s.为提高测量结果的可靠性和重复性,在不同载荷下,对同一种材料进行测试的数量为30 次.玻氏压痕实验使用纳米压痕仪(型号NHT2,Anton Paar 公司)在力控制模式下将压头垂直压入样品表面[43].实验的压痕形貌图如图1所示.
1.2.2 划痕实验
使用微米划痕仪(型号MST2,Anton Paar 公司)在渐进加载模式下利用半径为R=100µm 的Rockwell 压头进行划痕实验[13].单晶硅和碳化硅最大加载载荷分别为6 N 和10 N,划痕速率和长度均为4 mm/min 和2 mm.
图1 单晶硅和碳化硅维氏和玻氏压痕的光学图片Fig.1 Optical micrographs of the Vickers and Berkovich indentation for Si(111)and 4H-SiC
图2 维氏压头下Si(111)和4H-SiC 材料裂纹长度c 的概率密度函数图和正态性检验图Fig.2 Probability density function diagram and normal inspection diagram of crack length c of Si(111)and 4H-SiC under Vickers indenter
图2 给出了维氏压头下单晶硅在载荷0.98 N 和1.96 N、碳化硅在载荷2.94 N 和4.90 N 时c的经验概率密度函数图,其中虚线是使用高斯(Gauss)拟合得到的.单晶硅在0.98 N 和1.96 N 载荷下c的范围分别分布在15~21µm 和26~35µm 之间,最大值均约为最小值的1.4 倍,对于碳化硅而言是1.2 倍且数值分布都较为离散,这说明对c进行统计分析是十分有意义的,说明c的测量数据服从高斯分布.对于单晶硅和碳化硅材料的其他载荷下也可以得到相同的结论.表2 列出了维氏压头下材料在不同载荷下c的平均值cA和正态统计分析值cG,结果表明,两者几乎相同,说明在较大样本量(压痕测试次数为30 次)的情况下,对c取平均值的处理方法是合理的,这与Gong 等[44]的表述相吻合.
图3 是材料的纳米压痕硬度(HB)与维氏硬度(HV)随载荷的变化图,需要注意的是,这里硬度的定义均是载荷与残余压痕面积的比值,而不是载荷与投影面积的比值.结果表明,在玻氏压头下,随着载荷的增大,材料的压痕硬度整体表现出下降的趋势,表现出压痕尺寸效应(ISE)[45].碳化硅材料在较小区间内会出现硬度略微上升的现象(见图中箭头指向),维氏压头下也有相同现象产生,微小上升的原因归因于表面形貌的影响或测量带来的误差.研究发现,单晶硅的HV在0.49 N,0.98 N 和1.96 N 时均为9 GPa,与Lawn 等[11]的测试结果(9 GPa)一致;碳化硅的HV在1.96 N,2.94 N 和4.90 N 时分别为25.3 GPa,27.9 GPa 和28.1 GPa,与文献[46]的测试结果相接近.
图3 样品的维氏和纳米硬度随载荷的变化图Fig.3 Variation of hardness with the load under Vickers and Berkovich indenter
由于计算材料KIC所使用的压痕公式建立在不同的裂纹系统的基础上,有必要对其裂纹类型进行识别.图4 给出了两种实验材料的c/a值随载荷的变化图(图中还包含了两种裂纹系统模型图).研究发现,材料的c/a值随着压痕载荷P的增大而增大.在玻氏压头下,单晶硅和碳化硅材料的c/a值均小于2.5,根据经验判据(若c/a<2.5,材料表现为巴氏裂纹系统;若c/a≥2.5,则表现为中位裂纹系统),均表现为巴氏裂纹系统.在维氏压头下,碳化硅在测试载荷下均表现为中位裂纹系统,而单晶硅的裂纹系统与载荷有关,在0.49 N 时表现为巴氏裂纹系统,在0.98 N 和1.96 N 时则表现为中位裂纹系统.结果表明,如果不考虑压头形状对裂纹类型的影响,材料的裂纹类型则取决于实验施加的载荷,在低载荷下表现为巴氏裂纹,在高载荷下表现为中位裂纹,这与文献[47-48]的结论相吻合.
图4 样品的c/a 值及巴氏和中位裂纹模型Fig.4 c/a of samples and model of Palmqvist and Median crack
由于测试试样尺寸较小,无法使用常规宏观测试方法测量其KIC,这里是通过与划痕实验KIC测试结果进行比较分析,并依此来评价各种压痕经验公式的适用性.根据线弹性断裂力学(LEFM)理论,微米划痕法计算材料KIC的公式为[49-50]
式中Ft为切向力,单位为mN;p和A分别为压头与试样接触部分的投影周长和投影面积,单位分别为µm 和µm2.
图5(a)和图5(b)分别为单晶硅和碳化硅材料基于LEFM 模型求解的KIC随压入深度dp的变化图.KIC随dp的增大而增大,随后趋于稳定,之后发生剧烈波动.KIC在dp较大时发生剧烈波动的原因是由于在大载荷下材料发生严重地脆性破坏,并在压头下方产生碎屑导致Ft的测量值不稳定,如图5 插图所示.另外,由于LEFM 模型是在材料达到裂纹失稳的状况下建立的,即其在较大载荷下才能反映出材料的真实断裂韧性值,因此,取KIC-dp曲线稳定段的平均值作为单晶硅和碳化硅材料基于LEFM 模型的计算结果.单晶硅在KIC-dp曲线稳定段的平均值为0.96 MPa·,与四点弯曲法[51](0.91 MPa·)的测试结果接近.碳化硅的平均值为2.89 MPa·,与DT 法[26](3.1 MPa·)的测试结果基本一致,说明LEFM 模型能够很好地表征半导体材料的KIC.
图5 基于LEFM 模型的断裂韧性分析Fig.5 Fracture toughness analysis based on LEFM model
2.4.1 维氏压头下材料断裂韧性的讨论
图6 给出了维氏和玻氏压头下载荷P与c3/2之间关系的检查结果,表明单晶硅和碳化硅材料的载荷P与c3/2之间呈线性关系.两种材料断裂韧性的不同导致了对应直线斜率的差异,且韧性大的材料斜率小,韧性小的材料斜率大.
图6 维氏和玻氏压头下载荷P 与c3/2 的关系图Fig.6 Plots of c3/2 versus load under Vickers and Berkovich indenters
表3 给出了材料在维氏压头下使用表1 公式计算得到的KIC数据.结果表明,使用不同公式计算出来的KIC相差较大,因为这些公式都是由实验得到的经验方程,适用于不同的材料,所以,如果在计算KIC时选择了不适用的公式,则可能产生较大的误差.在不同载荷下,使用同一公式计算出来的KIC也有明显差异.因此,选择出适合材料的压痕公式十分重要.研究发现,对于单晶硅而言,在1.96 N 载荷下,使用P-3 公式计算得到的KIC与划痕实验测试结果一致;P-4 公式的计算结果随载荷的增大而增大,而在大载荷下单晶硅材料会发生严重地脆性破坏,KIC会被低估,故P-3 和P-4 公式并不适合被用于评估单晶硅材料的KIC.在所有测试载荷下,M-2,M-3,M-4 和M-6 公式的计算结果与划痕实验测试结果相接近,但在0.49 N 和0.98 N 载荷下,M-4 公式的计算结果与划痕实验结果几乎一致.因此,在维氏压头下使用M-4 公式评估单晶硅材料的KIC时建议使用此载荷条件.另外,M-2,M-3 和M-4 公式的计算结果相一致,这是因为M-3 和M-4 公式是在经典的Evans-Charles 方程——M-2 公式上进行修正和演化而来的,它们的表达式形式相似,只是系数略有不相同.对于碳化硅,尽管P-2 公式的计算结果与划痕实验结果一致,但其仅在1.96 N 载荷下适用,在测试载荷内的稳定性较差,故认为其不适用于评估碳化硅材料的KIC.在所有测试载荷下,M-5 和C-2 公式的计算结果与划痕结果接近,但在2.94 N 的载荷下,结果更加准确.因此,使用M-5 和C-2 公式评估碳化硅材料的KIC时建议使用2.94 N 的压痕载荷.
表3 维氏压头下单晶硅和碳化硅的KIC计算值Table 3 Fracture toughness of Si(111)and 4H-SiC determined by Vickers indenter
2.4.2 玻氏压头下材料断裂韧性的讨论
表1 中的13 个压痕公式是基于维氏压头发展出的半经验公式,不能直接应用于具有非对称的玻氏压头[43].因此,Ouchterlony[52]发展出了一个修正因子f(n)以说明从承受中心膨胀力的载荷点发出的裂纹数量对KIC的影响.
式中n为裂纹数量,且n≤9 时,维氏压头与玻氏压头的断裂韧性满足一次函数的关系[53].因此,在玻氏压头下使用表1 公式计算材料KIC时,需在公式前乘以修正系数kB=f(4)/f(3)=1.073.
图7 给出了在玻氏压头下计算材料KIC的示意图.通过拟合出在不同载荷下表1 压痕公式中各参数之间的线性关系,得到两种材料拟合直线的斜率kh(Si)和ks(4H-SiC),根据kh和ks与KIC的等式关系(对13 个压痕公式进行形式变换得到的)即可求出材料的KIC.这里需要注意的是,图7(a)和图7(h)中圈出的部分分别是材料在0.2 N 和0.125 N 以下得到的数据点,均出现了偏离线性关系的现象,表明在小载荷下不适合使用P-1,C-1 和C-3 公式评估材料的KIC,因而在拟合直线时这些点将被丢弃,这说明压痕公式都有其局限及适用范围.使用适用超薄玻璃的GB/T 37900—2019[54]中的压痕公式计算了两种实验材料的KIC.结果表明,GB/T 37900—2019[51]可以用来评估单晶硅材料的KIC,但对于碳化硅不适用,验证了上面的结论(见表3 中GB 列).
图7 玻氏压头下计算材料KIC的示意图Fig.7 Schematic diagram of calculating fracture toughness of materials under Berkovich indenter
图7 玻氏压头下计算材料KIC的示意图(续)Fig.7 Schematic diagram of calculating fracture toughness of materials under Berkovich indenter(continued)
表4 给出了玻氏压头下材料的KIC测量结果.表明,对于单晶硅而言,使用M-2、M-3 和M-4 公式计算出的KIC与划痕实验测试结果相接近,但M-4 公式的计算结果更加准确.对于碳化硅而言,适用的公式是P-3 和M-5.研究发现,在玻氏压头下,除P-3 公式外,其他适用公式的KIC计算结果明显小于划痕实验测试结果,说明利用玻氏压头评估半导体材料KIC时,在表1 中的维氏压痕公式的基础上乘以一个理论修正系数1.073 不足以获得准确的KIC结果,需要对其进行校正,校正后的修正系数kB见表5.结果表明,利用玻氏压头评估半导体材料的KIC时的压痕公式的修正系数应为1.13±0.01.
表4 玻氏压头下单晶硅和碳化硅材料的KIC计算值Table 4 The fracture toughness of Si(111)and 4H-SiC determined by Berkovich indenter
表5 利用玻氏压头评估单晶硅和碳化硅材料KIC时适用公式的修正系数kBTable 5 The correction coefficient kBof applicable equations when the fracture toughness of Si(111)and 4H-SiC is determined by Berkovich indenter
对比表3 和表4 的KIC测量结果,发现适用于中位裂纹系统的压痕公式更适合评估半导体材料的KIC,且没有一个公式可以同时适用于单晶硅和碳化硅材料,但对于同一材料,可以获得同时适用于两种压头的KIC计算公式.对于单晶硅而言,两种压头的共同适用公式是M-4,对于碳化硅而言,是M-5 公式.
本工作通过对Si(111)和4H-SiC(0001)开展维氏微米硬度和玻氏纳米压痕实验,讨论了材料在维氏压头下产生的裂纹长度的特点,使用13 个压痕公式计算出了材料的断裂韧性并与划痕实验测量结果进行了比对分析,验证了压痕法评估半导体材料断裂韧性的适用性.根据实验结果和上面的讨论,可以得出以下结论.
(1)由于在维氏压痕实验中产生的裂纹长度分布较为离散,有必要对其进行统计分析.结果表明,在大样本量(压痕次数为30 次)下,可以对裂纹长度进行取平均值的处理.
(2)利用维氏压头评估半导体材料断裂韧性时,单晶硅的建议载荷是0.49 N 和0.98 N,碳化硅是2.94 N.在玻氏压头下,可以通过拟合不同力下压痕公式各参数之间的线性关系来求解半导体材料的断裂韧性.
(3)压痕公式都有其局限及适用范围.对于单晶硅而言,两种压头的共同适用公式是M-4,对于碳化硅而言,是M-5 公式.