1. 引言
航空航天材料的结构性能越来越受到关注,航空发动机中的零部件需要在发动机运转时处在高温、高转速、高压的极端环境下安全工作 [1] [2] [3] 。为了研究出更符合此极端环境下的结构材料,需要材料在高温探针台中进行测试 [4] [5] [6] [7] 。而高温探台中的金属电学探针在高温下的物理性能对测试结果有影响,其中探针的杨氏模量影响尤其明显 [8] [9] 。
杨氏模量又称弹性模量。是宏观上描述物体受外作用力作用时,其抗形变的能力的物理量。测量材料的杨氏模量对于研究材料的力学特性有着重要的意义:无论是静态下固体的受力分析,还是研究力学的动态特性,都需要测量杨氏模量。杨氏模量不仅与材料本身参数相关,例如熔点和线性膨胀系数;还与材料的制作工艺和工作环境等因素有关 [10] [11] 。
目前,得益于社会飞速发展,材料能在更高的温度工作,但是高温测量杨氏模量非常困难。因此本文就杨氏模量在仿真模拟实验平台进行温度与杨氏模量的研究,即高温下材料不同位置的杨氏模量的变化规律,并且为高温探针台的探针受温度的影响和力学特性的变化提供一定的参考 [12] [13] 。
2. 理论模型
2.1. 信号频率的影响
图1为金属探针的结构图,探针长70 mm,针杆半径2.5 mm,针尖顶半径为0.1 mm。材料为钨钢。考虑到探针主要受温度和作用力影响,其具体参数见表1。
![](//html.hanspub.org/file/2-1270792x8_hanspub.png?20231212170632287)
Figure 1. Structure diagram of metal probe
图1. 金属探针结构图
![](Images/Table_Tmp.jpg)
Table 1. Specific parameters of metal probes
表1. 金属探针的具体参数
2.2. 网格划分及其精准性
采用COMSOL自带的网格划分功能,序列类型选择物理场控制网格,单元大小为超细化;共生成27,163个单元,并用该网格进行独立性检验(质量越接近1计算结果越准确)。平均单元质量为0.67,能在有限元分析计算中得到较为精准的结果 [14] 。
2.3. 计算假设
为了模拟金属探针在实验中准确的环境条件,采用了固体力学和固体传热两个物理场接口。固体力学主要分析金属探针在高温环境下受到热应力的影响,计算金属探针的应力分布;固体传热则是模拟金属探针在高温探针台里测试过程中温度分布情况。
设置好物理场后,便可在模拟平台进行计算,求解方法通常为有限元法。有限元法是目前最常见、计算能力强、涉及领域最多的一种工程分析算法。当需要模拟仿真和计算特定的工程问题时,建立基本方程和边界条件相对比较简单。但因为模型复杂的几何形状、材料本身的性能参数和对环境因素的考虑,导致精确解难以求出。于是,求近似解成为了一种可靠的解决途径。
本文主要探究金属探针的杨氏模量特性,测量杨氏模量需要观察金属探针的小形变。在这类问题中,运用线弹性有限元法更适合解决问题。金属探针的应力和应变都服从胡可定律的线性关系,微小应变和小位移具有线性关系。线弹性有限元法预算步骤如下:1) 结构离散化;2) 选择位移插值函数;3) 分析单元的力学特性;4) 集成所有单元平衡方程,得到整体结构平整方程;5) 由平衡方程求解点位移;6) 单元应变和应力的计算。
2.4. 边界条件的确定
1) 高温探针台的工作温度往往在1000摄氏度,金属探针的针尖接触高温物体,尾部则是被固定在夹具 [15] [16] 。为了简化模型和减少计算量,设置了整个针尖及倾斜壁为1000摄氏度,而不是针尖的顶端为1000摄氏度。因为针尖的顶端的单元数仅为4个,温度进行传热需要更长的步长,且热量传递效率低。故选择将整个针尖设置温度,有较好的温度传热模型,也减少了计算量。尾部则添加了固定的约束模拟金属探针被固定在夹具上。
控制方程:
(1)
(2)
式(1)中:p为金属探针密度,Cp为热容量,T为温度,为u流速场,
为不同方向的温度梯度,
为热传导速率,Q为热源,Qred为增加热源,式(2)中,q为热通量矢量场,k为物质的导热系数,单位均为国际单位。
2) 为了更一步简化模型,将不考虑环境温度的变化,环境温度固定在20摄氏度。
2.5. 多物理场的耦合
COMSOL Multiphysics除了能建立仿真模型,还能将多个物理场耦合在一起。常见的有固体力学和固体传热的耦合,简称热固耦合。还有流–固耦合(流体流动–固体力学)。对于本文设计的模拟实验则需要热固耦合。固体传热方程形式为瞬态,为了模拟金属探针的热传导的过程,设置了在180秒内每隔6秒就会输出一个结果。通过时间进行一个耦合,结果作为固体力学场的因变量,再和固体传热所设置相同的间隔计算出结果。固体力学的方程形式为稳态。
3. 结果与讨论
3.1. 应力和位移分布和结果分析
图2、图3表明,当金属探针加热到1000摄氏度时,经过180 s后金属探针在温度作用下由于不同位置存在一定的温差产生了热应力,尾部有向外膨胀的趋势。而其他位置整体受到方向向下的应力。在位移场上,因为尾部固定约束的影响,金属探针有向下的位移。这对应着应力向下的方向,越接近尖端位移越大。
见图4,每个点初始位置为参照点,绘制了不同位置的移动位移图。可以看出最大的位移为0.23 mm,而且位移变化趋于稳定。若加热时间增至360 s后,定义金属探针的伸长率为某时间点下金属探针的长度和初始状态的长度,得到图5。图5表明加热到200 s时,金属探针伸长率已经稳定为100.445%,说明探针的长度已稳定。认为金属探针承受温度剧烈的变化,会导致在高温晶圆测试中,探针的伸长会给晶圆带来异常针痕。所以在测试前要对金属探针进行预加热,来解决升温产生的问题。
![](//html.hanspub.org/file/2-1270792x13_hanspub.png?20231212170632287)
Figure 2. Stress distribution at 1000 degrees celsius
图2. 1000摄氏度下应力分布
![](//html.hanspub.org/file/2-1270792x14_hanspub.png?20231212170632287)
Figure 3. Displacement at 1000 degrees celsius
图3. 1000摄氏度下位移大小
![](//html.hanspub.org/file/2-1270792x15_hanspub.png?20231212170632287)
Figure 4. Plot of displacement at different positions of the probe over time
图4. 探针不同位置的位移随时间的变化图
![](//html.hanspub.org/file/2-1270792x16_hanspub.png?20231212170632287)
Figure 5. Graph of elongation of metal probes over time
图5. 金属探针的伸长率随时间的变化图
3.2. 温度分布及分析
设置针尖温度为1000摄氏度,模拟金属探针在高温探针台工作180 s后温度分布情况。且对距离针尖11 mm、20 mm、30 mm、40 mm、50 mm、60 mm、70 mm处的温度变化规律进行测量,得到图6。可以观察到金属探针的温度梯度变化情况。
![](//html.hanspub.org/file/2-1270792x17_hanspub.png?20231212170632287)
Figure 6. The variation of temperature at different positions of the probe over time
图6. 探针不同位置的温度随时间的变化规律
3.3. 杨氏模量计算和变化规律
在温度的作用下,金属探针受到微小应力和形变。通过模拟出不同位置对应的微小应力和形变,即
![](//html.hanspub.org/file/2-1270792x18_hanspub.png?20231212170632287)
Figure 7. Changes in Young’s modulus values at different positions of the probe
图7. 探针不同位置的杨氏模量值变化曲线
可计算出不同位置杨氏模量的数据。杨氏模量计算方程为:
(3)
其中是E金属探针的平均应力,单位为N/m2;ϵ是金属探针平均应变。
根据杨氏模量的计算方程得出下图。图7为金属探针不同位置的杨氏模量在加热过程中的变化曲线图。
对比图6和图7可以得出杨氏模量和温度呈线性关系。随着温度的升高,杨氏模量值逐渐下降。分析认为杨氏模量随温度升高而下降的原因主要是由于温度对分子之间的间距的影响,从而造成分子间引力作用减弱引起的变化。
图8为金属探针的杨氏模量与温度的关系。
![](//html.hanspub.org/file/2-1270792x20_hanspub.png?20231212170632287)
Figure 8. The relationship between the Young’s modulus of metal probes and temperature
图8. 金属探针杨氏模量与温度之间的关系
3.4. 模拟结果分析
通过模拟实验所得出的数据做分析。图3表明,越远离针尖一段温度跨度越大,温度变化的速率逐渐降低。温度分布的不均匀使得金属探针自身产生了热应力和细小的位移,通过分析这两者的数值即可按照公式计算出不同时间不同位置的杨氏模量值。图5表明,随着金属探针温度从常温升高接近1000摄氏度,杨氏模量从常温的706 GPa降到了656 GPa,下降了约7%。杨氏模量反映的是材料刚性,杨氏模量越小代表材料刚性差,易发生形变。考虑认为金属探针温度升高,材料内部原子热运动变得更加活跃,导致分子间的引力作用力下降,从而引起材料刚性变差。因此易形变即杨氏模量逐渐减小。
4. 总结
本文基于通过有限元法计算了金属探针的杨氏模量在高温情况下的变化特性。研究结果表明杨氏模量与加热时间和离针尖距离密切相关。首先,探针温度随时间变化呈非线性上升关系并趋于饱和,杨氏模量与温度几乎线性相关。其次,距离针尖越远,杨氏模量的下降幅度越大并趋于饱和。
因此,在高温电学测试前期的准备中,需要对探针进行预加热来减少探针伸长的影响,最终需寻找一个探针位置的平衡点以降低对实验结果的影响。
基金项目
本论文感谢科技部国家重点研发计划项目(2021YFB3801701)的支持。
参考文献