《一种计算集成电路芯片焊点热疲劳失效概率的方法.pdf》由会员分享，可在线阅读，更多相关《一种计算集成电路芯片焊点热疲劳失效概率的方法.pdf（15页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102955881 A(43)申请公布日 2013.03.06CN102955881A*CN102955881A*(21)申请号 201210422623.9(22)申请日 2012.10.30G06F 17/50(2006.01)(71)申请人温州大学地址 325000 浙江省温州市茶山高教园(72)发明人万毅(74)专利代理机构北京中北知识产权代理有限公司 11253代理人李雪芳(54) 发明名称一种计算集成电路芯片焊点热疲劳失效概率的方法(57) 摘要本发明提供了一种计算集成电路芯片焊点热疲劳失效概率的方法，包括以下步骤：（1）确定集成电路芯片焊点热疲劳可靠性的。
2、影响参数；（2）以正态分布随机地产生n组影响参数，并通过有限单元数值计算与每一组影响参数对应的最大等效的热应力；（3）对n组影响参数进行归一化处理；（4）得到影响参数-热应力之间的强非线性关系；（5）建立芯片动态时变焊料热疲劳强度模型；（6）通过CRS方法计算出芯片焊点的热疲劳失效概率。本发明解决了热疲劳失效及热疲劳可靠度计算难于建模型的问题。该方法速度快、简便易行，为复杂的集成电路芯片焊点热疲劳失效和热疲劳可靠性设计提供了一条新的途径。(51)Int.Cl.权利要求书1页 说明书8页 附图5页(19)中华人民共和国国家知识产权局(12)发明专利申请权利要求书 1 页 说明书 8 页 附图 5。
3、 页1/1页21.一种计算集成电路芯片焊点热疲劳失效概率的方法，其特征在于，包括以下步骤：（1）确定集成电路芯片焊点热疲劳可靠性的影响参数；所述影响参数为芯片本体的弹性模量(Ec)、芯片本体的热传导系数(TCc)、芯片本体的热膨胀系数(Cc)、焊点的弹性模量(Ew)、焊点热传导系数(TCw)、焊点热膨胀系数(Cw)、基板的弹性模量(Eb)、基板热传导系数(TCb)、基板热膨胀系数(Cb)、芯片工作热循环的平均温度(Tc)、芯片的长度(Lc)、芯片的宽度(Wc)、基板的长度(Lb)、基板的宽度(Wb)、焊球的高度(Hw)和焊球半径(Rw)；（2）以正态分布随机地产生n组影响参数，并通过有限单元数。
4、值计算与每一组影响参数对应的最大等效的热应力；其中n为大于1的正整数；（3）对步骤（2）中产生的n组影响参数进行归一化处理；（4）将步骤（3）中经过归一化处理后的数值作为输入，封装芯片焊点最大的热应力作为输出，将步骤（2）中的n组影响参数经过归一化处理后数据和对应的热应力作为支持向量机的训练样本，对最小二乘支持向量机进行训练和优化，得到影响参数-热应力之间的强非线性关系；（5）基于最小二乘支持向量机建立芯片动态时变焊料热疲劳强度模型；所述芯片动态时变焊料热疲劳强度模型以焊料的初始强度、工作的平均温度、温度循环的幅度和热载荷作用的次数作为最小二乘支持向量机的输入，焊料的热疲劳强度作为最小二乘支持。
5、向量机的输出；（6）根据步骤（4）中得到的影响参数-热应力之间的强非线性关系，以及步骤（5）中得到的芯片动态时变的疲劳强度的智能模型，通过CRS方法计算出芯片焊点的热疲劳失效概率。2.根据权利要求1所述的计算集成电路芯片焊点热疲劳失效概率的方法，其特征在于，步骤（6）包括以下子步骤：（6.1）根据影响参数-热应力之间的强非线性关系，将16个影响因素的一组随机数经过归一化处理后输入，得到焊点的最大的热应力Ts；（6.2）根据芯片动态时变焊料热疲劳强度模型，以正态分布产生随机疲劳强度值Tr；（6.3）计算g(Tr,Ts)=Tr-Ts；（6.4）根据CRS算法，不断地重复（6.1）-（6.3）步骤的。
6、计算，累积记录出现g0时的热疲劳失效次数L和计算的总次数N；当计算的次数大于预先设置的数值时，结束计算；（6.5）计算出芯片焊点的热疲劳失效概率Pf=L/N。权 利 要 求 书CN 102955881 A1/8页3一种计算集成电路芯片焊点热疲劳失效概率的方法技术领域0001 本发明属于电子系统集成电路领域，具体涉及一种计算集成电路芯片焊点热疲劳失效概率的方法。背景技术0002 在电子系统集成电路中，芯片是直接焊接到PCB或陶瓷基板表面，通过焊点实现器件和电路板之间的电气与机械连接。随着芯片集成度的提高和工作环境的复杂化，其发热密度也越来越高，集成电路的过热问题已经成为电子系统失效的主要原因。在。
7、随机和非对称循环热载荷作用下，集成电路芯片焊点的疲劳损伤逐渐累加，当累积损伤达到某一数值时，结构就会发生疲劳破坏而产生热疲劳失效，定量地预测和计算集成电路芯片焊点热疲劳失效概率和热疲劳可靠度是预防电路发生故障和保证电子设备稳定可靠工作最根本的途经，也是芯片进行热可靠性设计的前提和基础，它为集成电路芯片改进热设计提供依据，从而有针对性地进行热可靠性优化设计。0003 但是芯片疲劳破坏一般要经历材料局部出现高塑性区、萌生短裂纹、短裂纹扩展、长裂纹扩展到失效的复杂动态物理过程，焊点内部产生的热应力与外部环境、热载荷、芯片及基板几何尺寸及材料特性存在复杂的关系，很难用固定的数学模型进行描述，并且焊料的。
8、强度也是一种动态的时变过程，它的强度会随着使用时间的增加而在应力的作用下逐渐退化，它与焊料的初始强度，工作的平均温度，温度循环的幅度和热载荷作用的次数之间存在复杂的强非线性关系，采用传统的疲劳试验根本很难获得，应用在机械领域的估算模型如Gerber,Goodman,Soderberg,CepeHceH,Morrow和DING氏模型很难适合动态时变强非线性的集成电路芯片焊点的热疲劳可靠性计算和设计。以上原因使得集成电路芯片焊点定量的热疲劳失效概率和热疲劳可靠度的计算存在很大的困难，目前主要是凭经验进行一些定性的预测，一些学者提出一些理论和方法如Schaff提出的剩余强度模型，巴西国立学会的空间研。
9、究院的LuisAntBnio、Waak Bambace等人提出的边界积分方法，波特兰州立大学机械工程学院的马一东提出的分界面热对流数值方法等，它们都采用单一的和确定性的模型，不能从根本上满足集成电路芯片焊点热疲劳可靠性设计和计算的要求。发明内容0004 本发明针对上述现有技术的不足，提供了一种能够满足集成电路芯片焊点热疲劳可靠性设计要求的计算集成电路芯片焊点热疲劳失效概率的方法。0005 本发明是通过如下技术方案实现的：0006 一种计算集成电路芯片焊点热疲劳失效概率的方法，包括以下步骤：0007 （1）确定集成电路芯片焊点热疲劳可靠性的影响参数；所述影响参数为芯片本体的弹性模量(Ec)、芯片。
10、本体的热传导系数(TCc)、芯片本体的热膨胀系数(Cc)、焊点的弹性模量(Ew)、焊点热传导系数(TCw)、焊点热膨胀系数(Cw)、基板的弹性模量(Eb)、基板热传导系数(TCb)、基板热膨胀系数(Cb)、芯片工作热循环的平均温度(Tc)、芯片的长度(Lc)、芯片说 明 书CN 102955881 A2/8页4的宽度(Wc)、基板的长度(Lb)、基板的宽度(Wb)、焊球的高度(Hw)和焊球半径(Rw)；0008 （2）以正态分布随机地产生n组影响参数，并通过有限单元数值计算与每一组影响参数对应的最大等效的热应力；其中n为大于1的正整数；0009 （3）对步骤（2）中产生的n组影响参数进行归一化。
11、处理；0010 （4）将步骤（3）中经过归一化处理后的数值作为输入，封装芯片焊点最大的热应力作为输出，将步骤（2）中的n组影响参数经过归一化处理后数据和对应的热应力作为支持向量机的训练样本，对最小二乘支持向量机进行训练和优化，得到影响参数-热应力之间的强非线性关系；0011 （5）基于最小二乘支持向量机建立芯片动态时变焊料热疲劳强度模型；所述芯片动态时变焊料热疲劳强度模型以焊料的初始强度、工作的平均温度、温度循环的幅度和热载荷作用的次数作为最小二乘支持向量机的输入，焊料的热疲劳强度作为最小二乘支持向量机的输出；0012 （6）根据步骤（4）中得到的影响参数-热应力之间的强非线性关系，以及步骤（。
12、5）中得到的芯片动态时变的疲劳强度的智能模型，通过CRS方法计算出芯片焊点的热疲劳失效概率。0013 本发明专利提出了一种基于优化支持向量机（OSVM）、有限单元(FET)和计算机随机模拟算法(CRS)相融合的集成电路芯片焊点热疲劳失效概率和热疲劳可靠度定量计算的新方法（简称OFC法），利用均方差交叉验证网格搜索参数的优化支持向量机强大的学习、容错和泛化能力，建立了影响集成电路芯片焊点热疲劳可靠性参数与热应力之间的强非线性关系和焊料动态时变的热疲劳强度的智能模型，在些基础上，采用有限单元(FET)和计算机随机模拟算法(CRS)计算出集成电路芯片焊点热疲劳失效概率和热疲劳可靠度，本发明解决了集成。
13、电路芯片焊点内部产生的热应力与外部环境、热载荷、芯片及基板几何尺寸及材料特性之间复杂的强非线性关系、焊料强度复杂的动态过程以及热失效联合概率密度函数无法确定而使得热疲劳失效及热疲劳可靠度计算难于建模型的问题。该方法速度快、简便易行，为复杂的集成电路芯片焊点热疲劳失效和热疲劳可靠性设计提供了一条新的途径。附图说明0014 图1为单元几何形状图；0015 图2为芯片的焊料动态时变的非线性疲劳强度支持向量机智能结构图；0016 图3为OFC算法的流程图；0017 图4为OFC算法的总体结构图；0018 图5为焊球有限单元的热应力分布图；0019 图6为优化支持向量机对焊球热应力的计算误差图；0020。
14、 图7为焊点热疲劳可靠度随着使用时间的变化图；0021 图8为焊点热疲劳失效概率随着使用时间的变化图。具体实施方式0022 下面结合附图和具体实施方式对本发明做进一步详细的说明。说 明 书CN 102955881 A3/8页50023 本发明提供了一种计算集成电路芯片焊点热疲劳失效概率的方法，包括以下步骤：0024 （1）确定影响集成电路芯片焊点热疲劳可靠性的因素；0025 影响集成电路芯片焊点热疲劳可靠性和失效概率的因素有16个，它们包括芯片封装以及焊点的结构几何形状尺寸、材料特征参数、热载荷参数等。具体的参数如下：0026 芯片本体的弹性模量(Ec)，芯片本体的热传导系数(TCc)，芯片本。
15、体的热膨胀系数(Cc)，焊点的弹性模量(Ew)，焊点热传导系数(TCw)，焊点热膨胀系数(Cw)，基板的弹性模量(Eb)，基板热传导系数(TCb)，基板热膨胀系数(Cb)，芯片工作热循环的平均温度(Tc)，芯片的长度(Lc)，芯片的宽度(Wc)，基板的长度(Lb)，基板的宽度(Wb)，焊球的高度(Hw)，焊球半径(Rw)。0027 （2）进行芯片焊点热应力的有限单元数值计算；对各影响参数以正态分布N(,)随机地产生n组，通过有限单元计算芯片焊点每一组参数对应的最大等效的热应力。0028 为了更合理地选取影响因素随机参数，本专利采用(0，1)均匀试验设计法对参数进行取值，它忽略了正交设计中的整齐。
16、可比的特性，只考虑试验点在系统中的均匀分散性，这样既可以扩大因素的试验范围，又使试验次数增加较少，尤其对因素和水平数较多的复杂试验系统，均匀设计被认为是一种理想的试验设计方法，因此本专利采用均匀试验设计法来确定有限单元分析的样本，然后将(0，1)区间均匀分布的随机数转换成指定正态分布的随机抽样N(,)。转化的方法是：0029 0030 式中r1、r2二个相互独立的(0，1)均匀分布的随机变量，是芯片及焊点的热载荷、结构参数、材料参数的均值，2是芯片及焊点的热载荷、结构参数、材料参数的方差。0031 采用以上方法随机选取16个影响因素的n组热载荷、几何尺寸及材料参数，将它们作为有限单元分析程序的。
17、n组输入值，通过有限单元分析计算，可得到相应的n组最大的热应力响应值。有限单元热应力计算的方法如下：0032 （2.1）结构的离散化0033 建立基板、芯片和焊点的几何模型，给各实体赋材料属性，将芯片、基板和焊点的整体结构物分割成有限个单元体，使相邻单元体仅在结点处连接，而以这样的单元结合体去代替原来的结构。本发明采用八节点六面体的SOLID70单元，在结构分析中，对塑性材料本发明采用八节点六面体的SOLID45单元，对具有粘塑性的焊点材料，采用八节点六面体的VISC0107结构分析单元，每个节点都有三个自由度，分别为X,Y,Z三个方向上的节点位移，其单元示意图如图1所示。0034 （2.2）。
18、选择热位移模式0035 选定芯片、基板和焊点的整体结构所用单元后，导出用结点位移表示单元体内任一点热位移的关系式：0036 w=Nu(e)(1)0037 式中，w是单元内任一点的热位移列向量；u(e)是单元的结点的热位移列向量；而N称为形函数矩阵，其元素是位置坐标的函数。说 明 书CN 102955881 A4/8页60038 （2.3）分析单元的热力学特性0039 利用几何方程，由位移表达式(1)导出用结点热位移表示单元热应变的关系式:0040 =Bu(e)(2)0041 式中，是单元任一点的热应变列向量；B是单元热应变矩阵。0042 利用本构方程，由热应变的表达式(2)导出用结点表示单元热。
19、应力的关系式：0043 =DBu(e)(3)0044 式中，是单元内任一点的热应力列向量；D是与单元材料的热弹性矩阵。0045 利用变分原理，建立作用于单元上的结点力和结点位移之间的关系式，0046 即单元的平衡方程：0047 P(e)=k(e)u(e)(4)0048 式中，k(e)=BTDBdxdydz称为单元刚度矩阵；P(e)是单元等效结点力向量。0049 （2.4）建立热平衡方程，求解结点位移，计算单元中的热应力0050 集合所有单元的热平衡方程，利用最小势能原理建立结构的结点热荷载和结点热位移之间的关系式，即热的平衡方程:0051 0052 式中，是整体刚度矩阵；是热荷载列阵；是热位移。
20、列阵；符号上的“”表示在整体坐标系下，此前无“”各量表示在单元坐标系下。线性代数方程组(5)在代入边界条件后，经解算可求得所有未知的结点热位移。依据求得的结点热位移，由式(2)、式(3)可求出单元中任意单元的热应力。0053 （3）归一化处理；考虑到各种参数数值的差异性及分散性，为了降低它们之间由于量值差异过大造成的影响，对以上n组影响芯片及焊点热可靠性的参数通过线性变换映射到0，1范围内的相对含量，以降低各特征参数数据之间的互斥性。归一化处理方法如下：0054 0055 式中Ri表示芯片及焊点热疲劳可靠性的参数值，i=1,2,16。0056 （4）采用优化支持向量机建立芯片本体的弹性模量(E。
21、c)，芯片本体的热传导系数(TCc)，芯片本体的热膨胀系数(Cc)，焊点的弹性模量(Ew)，焊点热传导系数(TCw)，焊点热膨胀系数(Cw)，基板的弹性模量(Eb)，基板热传导系数(TCb)，基板热膨胀系数(Cb)，芯片工作热循环的平均温度(Tc)，芯片的长度(Lc)，芯片的宽度(Wc)，基板的长度(Lb)，基板的宽度(Wb)，焊球的高度(Hw)和焊球半径(Rw)的强非线性智能关系。输入的个数为16，分别为影响芯片及焊点热疲劳可靠性的16个因素经过归一化处理后的数值，输出为封装芯片焊点最大的热应力，将有限单元分析得到的n组影响参数-热应力数据作为支持向量机的训练样本，对最小二乘支持向量机进行训。
22、练和优化，以便得到支持向量机的最佳参数值，通过支持向量机训练确定影响参数-热应力之间的强非线性关系。0057 训练支持向量机的过程如下：设给定N个训练样本(xk,yk)|k=1,2,N,xkRd,ykR，xk是输入数据，yk是输出数据。函数估计问题可以描述求解下面问题：0058 说 明 书CN 102955881 A5/8页70059 0060 其中，是核空间映射函数，权矢量wRd，误差变量ekR，b是偏差量。损失函数J是SSE误差和规则化量之和，是控制对错分样本惩罚的程度的可调参数，它控制对超出误差样本的惩罚的程度，实现在训练误差和模型复杂度之间的折衷。与标准SVM不同，这里只有等式约束，且。
23、优化目标的损失函数是误差ek的二范数。核空间映射函数的目的是从原始空间中抽取特征，将原始空间中的样本映射为高维特征空间中的一个向量，以解决原始空间中线性不可分的问题。引入Lagrange函数，把有约束优化问题转化成无约束优化问题。0061 0062 即：0063 0064 其中，拉格朗日乘子kR，k0。对上式进行优化求解，根据KKT最优条件，把上式分别对w，b，ek和k求偏微分并令它们等于0，即：0065 0066 0067 得到：0068 0069 对于k=1,n消去w和ek，得到如下方程：0070 0071 其中0072 I=1;1T，Y=y1;ynT,=1;nT0073 即得到如下方程组。
24、：0074 说 明 书CN 102955881 A6/8页80075 优化支持向量机的过程如下：采用均方差交叉验证网格搜索法确定最小二乘支持向量机的最佳参数对(,)。具体的讲，包括如下子步骤：0076 (a1)确定合适的正则化参数集和核参数集。实验发现，按照指数增长方式生成两种参数集是一种有效的方法，例如，2-2,2-0,.,210，.，2-6,2-4,.,20,.，网格搜索简单直接，因为每一个参数对(,)是独立的，可以并行地进行网格搜索。0077 (a2)应用网格搜索法在正则化参数集和核参数集中选择一个参数对(,)，用该参数对进行交叉验证。首先用一个步长为22的(,)组合,得到学习精度最高的。
25、和的值。然后在这两个值旁的一定范围内进行一次更细致的网格搜索。在两个值的一定范围内把样本集D分为S组G1,G2,GS，把任意的S-1组作为训练集，剩余的一组作为验证集。通过选择不同的验证集，可重复S次。泛化性能可通过下式评价：0078 0079 式中：Gi是第i组验证集；yv是验证集的样本；是用D-Gi作为训练样本时得到的参数向量；是LS-SVM的输出。0080 (a3)循环选择参数对进行交叉验证，计算每个参数对的MSEcv，直到网格搜索停止。使得MSEcv最小的参数对(,)是最佳的，应用不完全交叉验证方法选择参数能够避免过拟合问题。0081 本发明通过样本集的学习，完成训练过程寻找到支持向量。
26、集，确定拉格朗日乘子k和偏移因子b，得到优化后的支持向量机。0082 （5）由于芯片的焊料在热载荷的反复作用下，它的疲劳强度会随着使用时间的增加逐渐退化，是一个复杂时变的动态过程，本发明提出了一种基于最小二乘支持向量机的动态时变焊料热疲劳强度模型。此模型把芯片的焊料所承受热载荷的作用时间（次数）和加载的热载荷大小与幅度进行了综合，得到了芯片的焊料动态时变的非线性疲劳强度智能模型。建模方法如下：0083 模型以焊料的初始强度Tr0，工作的平均温度T，温度循环的幅度T和热载荷作用的次数t，作为最小二乘支持向量机的输入，输入数为4，焊料的热疲劳强度Tr(t)，作为最小二乘支持向量机的输出，输出数为1。
27、，模型的结构如图2所示，支持向量机采用高斯核，即：0084 K(x，xi)=exp(-|x-xi|2/2)0085 核参数，正则化参数采用不完全交叉验证网格搜索优化策略进行优化。0086 （6）进行芯片焊点热疲劳失效概率的CRS计算和分析。在建立了芯片焊点的热载荷、结构参数、材料参数与热应力之间的支持向量机模型和焊料动态时变的疲劳强度的智能模型之后，就可通过CRS方法测试芯片焊点的热疲劳可靠度和热疲劳失效概率。0087 进行测试时，首先必须对16个影响因素及芯片焊料热疲劳强度影响因素进行随机抽样，以获得各随机因素的抽样值。随机抽样时，通常先在(0，1)区间上产生均匀分布随说 明 书CN 102。
28、955881 A7/8页9机数，然后将(0，1)区间均匀分布的随机数转换成指定正态分布的随机抽样。0088 然后将16个影响因素的一组随机数经过归一化处理后输入到已经训练好的优化支持向量机中，可得到焊点的最大的热应力Ts，并按芯片焊料热疲劳强度的智能模型以及正态分布产生的随机疲劳强度值Tr。把所得的随机值代入功能函数g(Tr,Ts)=Tr-Ts，计算g值。0089 根据CRS算法，不断地重复以上步骤的计算，累积记录出现g0热疲劳失效次数L和计算的总次数N。当计算的次数足够多时，如N100/Pf（Pf为芯片焊点的热疲劳失效概率）时，已经满足了精度，则计算结束。最后由Pf=L/N计算出芯片焊点的热。
29、疲劳失效概率，如果大于一定的临界值，则此芯片焊料热疲劳失效，本专利中，L取值500。图3和图4分别为OFC算法的流程图和OFC算法的总体结构图。0090 采用本发明专利中的融合OFC算法对微电子集成电路球栅阵列尺寸封装PBGA焊点(Sn-3.5Ag焊料)进行了实施例。表1为部分焊点热应力有限单元法的计算结果。图5为焊球有限单元的热应力分布，图6为支持向量机对焊球热应力的计算误差。图7为焊点热疲劳可靠度随着使用时间的变化，图8为焊点热疲劳失效概率随着使用时间的变化。通过融合OFC方法计算出当焊料在热载荷循环1500次后的热疲劳可靠度为0.7，热疲劳失效概率为30%。0091 表1焊点热应力有限单元的计算结果（部分数据）说 明 书CN 102955881 A8/8页100092 说 明 书CN 102955881 A10。