CN104820781A - 考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法 - Google Patents
考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法 Download PDFInfo
- Publication number
- CN104820781A CN104820781A CN201510226672.9A CN201510226672A CN104820781A CN 104820781 A CN104820781 A CN 104820781A CN 201510226672 A CN201510226672 A CN 201510226672A CN 104820781 A CN104820781 A CN 104820781A
- Authority
- CN
- China
- Prior art keywords
- load
- fatigue life
- solder joint
- bga
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,包括1.分析温循载荷顺序加载对BGA焊点热疲劳寿命的影响及原因;2.根据温循载荷顺序加载对BGA焊点热疲劳寿命的影响并参考Paris公式对Darveaux模型进行改进,构建考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型;3.确定考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型中各参数的拟合确定方法;4.进行考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型的分析验证,分别进行单温度循环载荷下与Darveaux模型疲劳寿命预测结果的对照分析,以及温循载荷顺序加载下与Miner理论疲劳寿命预测结果的对照分析;5.完成所述预测方法的构建。
Description
技术领域
本发明涉及考虑不同温度循环载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,属于失效物理技术领域。
背景技术
随着电子封装技术的不断发展,应用电子产品的领域日益广阔,电子产品使用过程中所遭受的环境条件愈加复杂、恶劣,可能出现环境突变、工作模式变换等情况,使电子产品处于载荷顺序加载的载荷条件下。BGA焊点(球栅阵列封装焊点)作为保障电子产品可靠性的关键环节,其疲劳寿命的预测评估对产品而言具有重要意义。然而,对于焊点疲劳中最主要的热疲劳,现有温度循环条件下的寿命预测模型不能考虑加载条件变化时载荷顺序的影响,而应用Miner理论进行线性累积结果会有较大误差,相应的研究工作又十分有限,没有便于应用的成熟方法。
鉴于此,需要研究不同温度循环载荷顺序加载对焊点的疲劳寿命的影响,提出考虑该影响的改进寿命预测模型,进而形成一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法。在总结已有焊点疲劳寿命预测模型及累计损失理论的基础上,结合相关试验数据,分析了温循载荷顺序加载对焊点疲劳寿命的影响,并构建了考虑该影响的疲劳寿命预测模型,进而应用试验数据与Miner理论预测结果进行了分析验证,最终提出了一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法。
发明内容
本发明的目的是:提供一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,提高温循载荷顺序加载下BGA焊点热疲劳寿命预测结果的准确性。
首先引入七个定义。
定义1:温循载荷顺序加载:指不同温循载荷条件按照某种顺序和循环时间依次加载的载荷条件。
定义2:Darveaux模型:Darveaux针对BGA封装,将试验测得的低周疲劳裂纹萌芽及生长率和非弹性功联系起来,提出的基于损伤累积的能量模型。
定义3:Miner理论:Palmgren-Miner线性累积损伤准则。
本发明的技术方案:分析温循载荷顺序加载对BGA焊点热疲劳寿命的影响;构建考虑温循载荷顺加载影响的BGA焊点热疲劳寿命预测模型;提出考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型参数拟合方法;应用试验数据并与Miner理论结果对照,进行上述模型的分析验证;结合分析结果,完成考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法的构建。
本发明是一种考虑温循载荷顺加载影响的BGA焊点热疲劳寿命预测方法,其步骤如下:
步骤1,根据试验数据,分析温循载荷顺序加载对BGA焊点热疲劳寿命的影响及原因;
步骤2,根据温循载荷顺序加载对BGA焊点热疲劳寿命的影响并参考Paris公式对Darveaux模型进行改进,构建考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型;
步骤3,结合模型特点,确定考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型中各参数的拟合确定方法;
步骤4,利用试验数据进行考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型的分析验证,分别进行单温度循环载荷下与Darveaux模型疲劳寿命预测结果的对照分析,以及温循载荷顺序加载下与Miner理论疲劳寿命预测结果的对照分析;
步骤5,在上述工作的基础上,完成考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法的构建。
其中,在步骤1中所述的“温循载荷顺序加载对BGA焊点热疲劳寿命的影响及原因”,温循载荷顺序加载对BGA焊点热疲劳寿命的影响为:当温循载荷由高到低顺序加载时,BGA焊点热疲劳试验中测得的疲劳寿命比Miner理论预测的疲劳寿命短,而当温循载荷由低到高顺序加载时,试验测得的疲劳寿命比Miner理论预测的疲劳寿命长;产生上述影响的原因为:焊点热疲劳失效的早期,焊点内部受滑移带形成、塑形变形的影响会产生循环软化现象,从而降低热不匹配产生的应力,延缓焊点热疲劳发展,并且在温循载荷条件较低的情况下该现象更为显著,另外较低的载荷条件能够缓解细微缺陷,而较高的载荷条件会加剧细微缺陷的演化,促进焊点裂纹的形成,综上不同温循载荷顺序加载会对焊点热疲劳寿命产生不同影响。
其中,在步骤2中所述的“根据温循载荷顺序加载对BGA焊点热疲劳寿命的影响并参考Paris公式对Darveaux模型进行改进”,改进内容为:将Darveaux模型裂纹扩展阶段中描述的平均裂纹扩展速率改为与循环数相关的瞬时裂纹扩展速率,引入与承受载荷条件相关的载荷因子以及与裂纹萌生阶段承受载荷历史相关的状态因子。
其中,在步骤3中所述的“各参数的拟合确定方法”,拟合确定方法为:裂纹萌生阶段相关参数直接采用Darveaux模型裂纹萌生阶段相关参数,裂纹扩展阶段相应的载荷因子和状态因子采用利用试验数据拟合的方法确定。
其中,在步骤4中所述的“以及温循载荷顺序加载下与Miner理论疲劳寿命预测结果的对照分析”,其中温循载荷顺序加载为:焊点在低温循载荷条件下循环一定周期数后加载高温循载荷条件直至焊点失效以及高温循载荷条件下循环一定周期数后加载低温循载荷条件直至焊点失效。
其中,在步骤5中所述的“完成考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法的构建”,该寿命预测方法包括:假设适用条件、考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型以及模型中相关参数的确定方法。
本发明与现有技术相比有如下优点:
第一,本发明考虑温循载荷顺序加载对BGA焊点疲劳寿命的影响,相应载荷条件下应用本方法得出的BGA焊点热疲劳寿命与传统方法相比更符合实际结果、更为准确。
第二,本发明给出了考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法实施应用细节,具有工程应用价值。
附图说明
图1为本发明所述的考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法流程图;
图2为温循载荷条件1到温循载荷条件2顺序加载下裂纹扩展循环数与裂纹长度关系示意图;
图3为焊点裂纹长度与裂纹扩展循环数关系示意图;
图4为条件下,考虑温循载荷顺序加载影响的热疲劳模型和Darveaux模型裂纹扩展循环数与裂纹长度拟合结果对照图;
图5为考虑温循载荷顺序加载影响的热疲劳模型和Darveaux模型裂纹扩展循环数与裂纹长度拟合结果对照图;
图6为温循载荷由高到低顺序加载下考虑温循载荷顺序加载影响的热疲劳模型与Miner准则预测结果对比;
图7为温循载荷由低到高顺序加载下考虑温循载荷顺序加载影响的热疲劳模型与Miner准则预测结果对比;
图2中符号说明如下:
n1:当载荷条件变化时,焊点已经历的裂纹扩展阶段对应的循环数;
n2:当载荷条件变化时已完成循环数n1所对应同等裂纹长度在载荷条件2下需完成的循环数;
a1:载荷条件1下循环数n1对应的裂纹长度;
a2:载荷条件2下循环数n2对应的裂纹长度。
具体实施方式
本发明所述的考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,见图1所示,该方法的具体实施方式步骤如下:
步骤1,根据试验数据,分析温循载荷顺序加载对BGA焊点热疲劳寿命的影响及原因。
温循载荷顺序加载对BGA焊点热疲劳寿命的影响为:当温循载荷由高到低顺序加载时,BGA焊点热疲劳试验中测得的疲劳寿命比Miner理论预测的疲劳寿命短,而当温循载荷由低到高顺序加载时,试验测得的疲劳寿命比Miner理论预测的疲劳寿命长。
产生上述影响的原因为:焊点热疲劳失效的早期,焊点内部受滑移带形成、塑形变形的影响会产生循环软化现象,从而降低热不匹配产生的应力,延缓焊点热疲劳发展,并且在温循载荷条件较低的情况下该现象更为显著,另外较低的载荷条件能够缓解细微缺陷,而较高的载荷条件会加剧细微缺陷的演化,促进焊点裂纹的形成,综上不同温循载荷顺序加载会对焊点热疲劳寿命产生不同影响。
步骤2,根据温循载荷顺序加载对BGA焊点热疲劳寿命的影响并参考Paris公式对Darveaux模型进行改进,构建考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型;
Paris公式:
式中,C,n为试验测得的材料常数;△K为应力场强度因子幅;N为裂纹扩展阶段的循环数;a为第N次循环时对应的裂纹长度;Y为几何因子,其值与裂纹形态、试样形状及加载方式有关,一般在1~2之间。
Darveaux模型:
式中:No为裂纹萌生阶段的循环数;△Wave为体积均匀化的粘塑性应变能密度增量;a为断裂特征长度;为裂纹扩展率;Nf为特征寿命(失效概率为63.2%的循环数);k1、k2、k3、k4为相关系数。
步骤201,结合温循载荷顺序加载对BGA焊点热疲劳寿命的影响提出考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型的前提假设条件。
裂纹到达一定长度时焊点发生失效:结合BGA封装器件特性,将疲劳裂纹演化到使器件产生结构功能失效时对应的裂纹长度作为失效判别准则。
焊点疲劳失效分为裂纹萌生和裂纹扩展两个阶段:认为焊点疲劳失效包括裂纹萌生和裂纹扩展两个阶段,即:Nf=No+N,其中No为裂纹萌生寿命,N为裂纹扩展寿命,Nf为裂纹总疲劳寿命。
裂纹萌生阶段存在循环软化现象:认为温度循环载荷条件下,裂纹萌生阶段焊点的弹性模量会随循环数的增加而逐渐降低,即存在循环软化现象,而裂纹扩展阶段无这种循环软化现象。
裂纹扩展行为受裂纹萌生阶段载荷条件影响:认为焊点裂纹萌生阶段的载荷条件会影响其内部缺陷,即低载荷能够缓解内部缺陷,根据断裂力学应力和强度原则,能够延长焊点疲劳寿命;而高载荷不能缓解内部缺陷,甚至促进缺陷的扩大,从而缩短焊点疲劳寿命。
焊点疲劳失效中无过载停滞现象:过载是指材料遭受恒幅疲劳载荷作用中突然受到更高载荷作用,随后又恢复到初始载荷条件的加载过程。裂纹闭合模型对过载这种载荷变化瞬间带来的次序效应进行了分析,认为过载会使裂纹扩展产生短暂停滞,即产生过载停滞现象,从而影响材料疲劳寿命。但相较于材料所受载荷史对其结构参数产生的影响,过载停滞现象影响更为短暂、产生次数也有所限制,因此可忽略该现象产生的影响,模型构建中不对其进行考虑。
步骤202,在上述前提假设条件下对Darveaux模型进行改进,构建考虑温循载荷顺序加载影响的裂纹萌生及裂纹扩展阶段的模型公式。
裂纹萌生阶段:
式中,No为裂纹萌生阶段的循环数;△Wave为体积均匀化的粘塑性应变能密度增量;k1、k2为相关系数。
裂纹扩展阶段:
式中,da/dN为裂纹扩展阶段对应的裂纹瞬时扩展速率;N为裂纹扩展阶段循环数;α,β为参数常量,α为载荷因子,与BGA焊点所承受的温循载荷条件相关,β为状态因子,与BGA焊点裂纹萌生阶段所受的温循载荷相关,载荷越低其值越小。特别地,当β=0时da/dN=α为常量,此时模型与Darveaux模型一致。
单温循载荷条件下的疲劳寿命预测公式:
式中,a为特征裂纹长度;Nf为裂纹扩展到特征裂纹长度a对应的循环数,即焊点疲劳寿命;No为裂纹萌生阶段的循环数;△Wave为体积均匀化的粘塑性应变能密度增量;k1、k2为相关系数;da/dN为裂纹扩展阶段对应的裂纹瞬时扩展速率;N为裂纹扩展阶段循环数;α为载荷因子常量,与BGA焊点所承受的温循载荷条件相关;β为状态因子,与BGA焊点裂纹萌生阶段所受的温循载荷相关。
步骤203,在上述公式模型的基础上,推到考虑温循载荷顺序加载的BGA焊点疲劳寿命预测模型。
以载荷条件1和载荷条件2两种不同温循载荷顺序加载为例,假设BGA焊点在温度循环载荷条件1下进入裂纹扩展阶段并在该阶段经过n1循环后,温度循环载荷条件变化为条件2加载直至焊点失效,其过程见图2,此过程中引入了条件1到条件2载荷顺序(由高载荷到低载荷或低载荷到高载荷)施加的影响,按照之前提出的公式模型,该情况下焊点热疲劳的寿命预测公式推导如下:
BGA焊点在温度循环载荷条件1下进入裂纹扩展阶段后,当温度循环载荷条件发生变化时,若不考虑载荷顺序产生的影响,同等裂纹长度下(a2=a1),载荷条件2对应的裂纹瞬时扩展速率为:
式中:n2为当载荷条件变化时已完成循环数n1所对应同等裂纹长度在载荷条件2下需完成的循环数;da2/dN2代表a2=a1时载荷条件2下的裂纹瞬时扩展速率;α2,β2为载荷条件2下裂纹扩展公式中的常量参数。
但考虑到建模时循环软化及裂纹扩展行为假设,此时受载荷顺序影响状态因子β取决于裂纹萌生阶段原载荷条件,即β=β1(β1为载荷条件1下裂纹扩展公式中的状态因子常量),因此实际裂纹瞬时扩展速率应为:
此时与循环数n1和循环数n2对应的裂纹长度分别为:
又由于循环数n1和循环数n2对应的裂纹长度相等,即a1=a2,因此由式(10)和式(11)联立可求得循环数n2与循环数n1的关系式,将其带入式(9)可得裂纹扩展速率与循环数n1的关系式:
令
假设以后裂纹扩展速率满足上述形式,并且该曲线(图2中载荷1→载荷2曲线横坐标n1之后部分)过图2中点(n1,a1),则载荷条件改变后,将式(12)进行积分可得循环数n1之后对应裂纹增长量△a为:
将上式与式(10)联立,可得总裂纹长度a与裂纹扩展循环数n之间的关系为:
由上式可求得裂纹扩展循环数n关于总裂纹长度a的表达式,用符号Np代替符号n,及裂纹扩展寿命Np为:
所以由式(5)和式(15)可得焊点热疲劳失效总寿命为:
特别地,当β1=β2=0时
综上所述,BGA焊点在温循载荷条件1到温循载荷条件2顺序施加下的改进疲劳寿命预测模型为:
式中: k1,k2为裂纹萌生阶段相关系数,△Wave1为载荷条件1下平均粘塑性应变能密度,n1为载荷条件1下裂纹扩展阶段循环数,α1,α2,β1,β2分别为载荷条件1和载荷条件2下裂纹扩展阶段相关系数。
步骤3,结合模型特点确定考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型中参数k1,k2选用Darveaux模型中给定值,而各载荷条件下对应的参数α,β值由以下方法确定。
基于循环数与裂纹长度观测试验的参数确定方法:
试验中可直接观测循环数与裂纹长度关系时,用一种载荷条件的试验数据即可求得对应载荷条件下的α,β,图3为试验记录的q个相同BGA器件在某种温度循环载荷条件下角焊点(位于器件四角的焊点,通常为热疲劳最先失效之处)观测到的裂纹长度与循环数关系示意图。
设共进行了p次记录,各记录点对应循环数为n1,n2,…,np,第i个记录点记录的各样本裂纹长度为ai1,ai2,…,aiq,则各记录点裂纹长度均值为:
n次循环加载后的裂纹长度a(n)与循环数n间的关系为:
两边取对数,可在对数坐标系中表示为直线关系:
由最小二乘法知参数α,β的最小二乘估计为:
此时,用ni次循环下裂纹长度均值ai替代式中ai,并带入ni即可求得某种温循载荷下改进模型中参数α,β。
基于循环数与电阻值观测试验的参数确定方法:
IPC等标准及文献常将样本相对阻值漂移10%或20%等定为焊点失效的判别标准。由于单个焊点阻值很小,试验中常采用检测多个焊点组成的菊花链结构的阻值变化来判别焊点是否失效。为获得两种载荷条件下对应参数的值,需开展温循载荷条件1、温循载荷条件2、温循载荷由条件1到条件2加载以及温循载荷由条件2到条件1加载四组试验,获取相关的失效数据,其具体确定方法如下。
开展新试验前,由单种温循载荷条件下焊点的失效数据带入Darveaux模型可求得对应条件下焊点的平均寿命,由此可建立下列等式:
另外,开展温循载荷由条件1到条件2加载以及由条件2到条件1加载的两组试验,将相应数据带入式(18),可得如下关系等式:
式中,
将式(25)至式(27)联立可得求解参数β1,β2的超越方程组如下:
其中 可由式(24)、式(25)求出。
由于该方程组为超越方程组,无法求得β1,β2的解析式,此时可应用Matlab软件进行计算求解,其方法如下。
首先编写方程组函数:
function F=myfun(x)
A=1;%表示的值。
B=1;%表示的值。
C=3;%表示的值。
D=3;%表示的值。
E=3;%表示的值。
G=3;%表示的值。
a=x(1);%代表参数β1
b=x(2);%代表参数β2
F=[
(((1+b+a.*(a+1))./((b+1).^2)).*(1-E.^(a+1)).*(D.^b)
.*(E.^(a+1)./(b+1))).^((1+b)./(1+b+a.*(a+1)))-A;
(((1+a+b.*(b+1))./((a+1).^2)).*(1-G.^(b+1)).*(C.^a)
.*(G.^(b+1)./(a+1))).^((1+a)./(1+a+b.*(b+1)))-B;
];
end
其中A,B,C,D,E,G为方程组中常数,可见函数%之后说明,此处随机赋予了量值以便说明。而a,b为待解的参数β1,β2。
编写函数后,应用以下命令可求得β1,β2的值。
[x,fval]=fsolve(myfun,[00]);
在得到β1,β2的值后,可由下式求得参数α1,α2的值。
综上,可求得模型中温循载荷条件1和载荷条件2下α1,β1,α2,β2所有参数的值。
步骤4,利用试验数据进行考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型的分析验证。
此处采用的试验数据为Darveaux温度循环试验所得循环数与裂纹长度试验数据。试验分为和两种循环条件,各条件下又可按照裂纹演化理论分为萌生阶段和扩展阶段,其中各载荷条件下萌生阶段的循环数为根据数据关系图测算所得,扩展阶段循环数不包含萌生阶段的循环数。具体试验数据如下。
表Darveaux温度循环试验所得循环数与裂纹长度试验数据
步骤401,进行单温度循环载荷下虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型与Darveaux模型疲劳寿命模型预测结果的对照分析。
按照步骤2和步骤3中推得模型公式,以及Darveaux模型可得两种载荷条件下相应的表达式如下。
温循载荷条件下,考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型各参数及表达式为:
拟合相关系数r=0.984264。
温循载荷条件下,Darveaux模型各参数及表达式为:
da1/dN1=6.641×10-4,a1=0.000664np-0.03294,Np1=732;
拟合相关系数r=0.986766。
利用以上求得的关系式,可得考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型与Darveaux模型的拟合曲线结果,见图4。可以发现两种模型的拟合结果相差不大,相关系数方面也相差不多,而考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型在前半部分拟合较好,在后半部分在相同裂纹长度下,该模型预测循环数结果要比Darveaux模型略高。
温循载荷条件下,考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型各参数及表达式为:
拟合相关系数r=0.992467。
温循载荷条件下,Darveaux模型各参数及表达式为:
da2/dN2=1.285×10-4,a2=0.0001285np-0.01529,Np2=3644;
拟合相关系数r=0.986792。
利用以上求得的关系式,可得考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型与Darveaux模型的拟合曲线结果,见图5。可以发现两种模型的拟合结果相差不大,相关系数方面也相差不多,考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型在后半部分在相同裂纹长度下,预测循环数结果要比Darveaux模型略高,但相对误差不会超过5%。
以上结果可以证明,在单一温度循环载荷条件下,考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型足以完成响应的寿命预测工作,其结果与Darveaux模型与测结果相近,略微偏高。
步骤402,进行温循载荷顺序加载下考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型与Miner理论疲劳寿命预测结果的对照分析。
将上一步求得的参数值带入式(18)中,可得温循载荷条件有高到低顺序加载时(载荷条件下达到裂纹扩展阶段并在裂纹扩展阶段经历n1次循环后加载载荷条件)焊点的疲劳寿命为:
同理,当温循条件由低到高加载时,焊点疲劳寿命为:
以上两式中,a为裂纹最大长度,约为0.453mm;n1为载荷变化时对应的裂纹扩展循环数。
根据步骤402分析结果可知,载荷条件下预测的裂纹扩展寿命为763次循环,裂纹萌生寿命为35次循环,即疲劳寿命为798次循环;而载荷条件下裂纹扩展寿命为3744次循环,裂纹萌生寿命为300次循环,即疲劳寿命为4044次循环。因此,按照Miner线性累积损伤理论推得不同加载顺序下焊点疲劳寿命分别为:
式中,n1为载荷变化时对应的裂纹扩展循环数。
由此分别绘制由高到低加载和由低到高顺序加载条件下焊点疲劳寿命曲线,见图6和图7。可以看出,当由高到低顺序加载时改进的疲劳寿命预测模型计算得出的疲劳寿命要小于Miner理论结果,若以Miner理论结果为基准,其比值在0.388~1范围内浮动;而由低到高顺序加载时预测寿命要大于Miner理论结果,若以Miner理论结果为基准,其比值在1~1.769范围内浮动。而对于不同载荷顺序加载,各学者普遍认可当载荷由高载荷到低载荷顺序施加时,Miner理论预测的疲劳寿命与试验结果相比偏大,而当载荷由低载荷到高载荷顺序施加时,Miner理论预测的疲劳寿命与试验结果相比偏小,而且各载荷条件下材料损坏时按Miner理论计算的实际损伤在0.25~4之间。由此说明了提出的考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型得出的焊点在不同温循载荷顺序施加下的热疲劳寿命基本符合各研究中的试验结果,证实了该疲劳寿命预测模型的合理性和正确性。
此外可以发现,考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型中当裂纹扩展阶段循环数越小时改进模型预测结果与Miner理论相差越大,这种差异可能由于假设材料软化仅发生在裂纹萌生阶段而被放大了,实际差异可能会相对小些,因此该模型也有一定的适用条件,一般认为该疲劳寿命预测模型在第一种温循载荷条件加载达到相对寿命的30%以上时较为适用。
步骤5,在上述工作的基础上,完成考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法的构建;
考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,主要由步骤2和步骤所述内容组成,即考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型参数的确定以及模型公式的应用计算,该方法主要适用于两种或两种以上不同温循载荷条件顺序加载的情况,当第一种载荷条件加载时间较长,超过焊点相对寿命的30%时,应用该方法预测的BGA焊点疲劳寿命更为准确。
Claims (6)
1.一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,其特征在于:其步骤如下:
步骤1,根据试验数据,分析温循载荷顺序加载对BGA焊点热疲劳寿命的影响及原因;
步骤2,根据温循载荷顺序加载对BGA焊点热疲劳寿命的影响并参考Paris公式对Darveaux模型进行改进,构建考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型;
步骤3,结合模型特点,确定考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型中各参数的拟合确定方法;
步骤4,利用试验数据进行考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型的分析验证,分别进行单温度循环载荷下与Darveaux模型疲劳寿命预测结果的对照分析,以及温循载荷顺序加载下与Miner理论疲劳寿命预测结果的对照分析;
步骤5,在步骤1至步骤4的基础上,完成考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法的构建。
2.根据权利要求1所述的一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,其特征在于:在步骤1中所述的“温循载荷顺序加载对BGA焊点热疲劳寿命的影响及原因”,温循载荷顺序加载对BGA焊点热疲劳寿命的影响为:当温循载荷由高到低顺序加载时,BGA焊点热疲劳试验中测得的疲劳寿命比Miner理论预测的疲劳寿命短,而当温循载荷由低到高顺序加载时,试验测得的疲劳寿命比Miner理论预测的疲劳寿命长;产生上述影响的原因为:焊点热疲劳失效的早期,焊点内部受滑移带形成、塑形变形的影响会产生循环软化现象,从而降低热不匹配产生的应力,延缓焊点热疲劳发展,并且在温循载荷条件较低的情况下该现象更为显著,另外较低的载荷条件能够缓解细微缺陷,而较高的载荷条件会加剧细微缺陷的演化,促进焊点裂纹的形成,综上不同温循载荷顺序加载会对焊点热疲劳寿命产生不同影响。
3.根据权利要求1所述的一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,其特征在于:在步骤2中所述的“根据温循载荷顺序加载对BGA焊点热疲劳寿命的影响并参考Paris公式对Darveaux模型进行改进”,改进内容为:将Darveaux模型裂纹扩展阶段中描述的平均裂纹扩展速率改为与循环数相关的瞬时裂纹扩展速率,引入与承受载荷条件相关的载荷因子以及与裂纹萌生阶段承受载荷历史相关的状态因子。
4.根据权利要求1所述的一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,其特征在于:在步骤3中所述的“各参数的拟合确定方法”,拟合确定方法为:裂纹萌生阶段相关参数直接采用Darveaux模型裂纹萌生阶段相关参数,裂纹扩展阶段相应的载荷因子和状态因子采用利用试验数据拟合的方法确定。
5.根据权利要求1所述的一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,其特征在于:在步骤4中所述的“以及温循载荷顺序加载下与Miner理论疲劳寿命预测结果的对照分析”,其中温循载荷顺序加载为:焊点在低温循载荷条件下循环一定周期数后加载高温循载荷条件直至焊点失效以及高温循载荷条件下循环一定周期数后加载低温循载荷条件直至焊点失效。
6.根据权利要求1所述的一种考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法,其特征在于:在步骤5中所述的“完成考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测方法的构建”,该寿命预测方法包括:假设适用条件、考虑温循载荷顺序加载影响的BGA焊点热疲劳寿命预测模型以及模型中相关参数的确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510226672.9A CN104820781B (zh) | 2015-05-06 | 2015-05-06 | 考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510226672.9A CN104820781B (zh) | 2015-05-06 | 2015-05-06 | 考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104820781A true CN104820781A (zh) | 2015-08-05 |
CN104820781B CN104820781B (zh) | 2017-09-29 |
Family
ID=53731076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510226672.9A Active CN104820781B (zh) | 2015-05-06 | 2015-05-06 | 考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104820781B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156421A (zh) * | 2016-07-01 | 2016-11-23 | 电子科技大学 | 基于脉冲涡流热成像的电子封装焊点热疲劳寿命预测方法 |
CN106568660A (zh) * | 2016-10-14 | 2017-04-19 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
CN107203666A (zh) * | 2017-05-22 | 2017-09-26 | 北京航空航天大学 | 一种bga焊点热疲劳寿命的预测方法及系统 |
CN107273632A (zh) * | 2017-06-28 | 2017-10-20 | 北京航空航天大学 | 一种非累积冲击下多部件系统寿命预测方法及装置 |
CN108287978A (zh) * | 2018-03-26 | 2018-07-17 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 焊点设计参数的优化方法 |
CN108664669A (zh) * | 2017-03-27 | 2018-10-16 | 无锡市五十五度科技有限公司 | 一种bga焊点热寿命预测方法、测试平台及测试机箱 |
CN109116140A (zh) * | 2018-07-16 | 2019-01-01 | 中国航空综合技术研究所 | 一种用于pbga封装器件的测试方法 |
CN109387533A (zh) * | 2017-08-08 | 2019-02-26 | 北京航空航天大学 | 焊点在不规则温循剖面下损伤累积法则实验方法 |
CN110059412A (zh) * | 2019-04-19 | 2019-07-26 | 上海工程技术大学 | 金属材料寿命预测方法 |
CN111414709A (zh) * | 2020-04-14 | 2020-07-14 | 上海交通大学 | 一种基于神经网络的超级电容器循环寿命和容量衰减的早期预测方法 |
CN111521507A (zh) * | 2020-04-30 | 2020-08-11 | 江苏师范大学 | 一种环保不锈钢焊缝热疲劳表面裂纹扩展速率的试验方法 |
CN113189136A (zh) * | 2021-03-29 | 2021-07-30 | 北京航空航天大学 | 一种针对宇航用球栅阵列器件装联工艺的温度循环试验分级评价方法 |
CN113495189A (zh) * | 2021-09-08 | 2021-10-12 | 深圳荣耀智能机器有限公司 | 评估电子装联材料可靠性的测试方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000046905A (ja) * | 1998-07-29 | 2000-02-18 | Sharp Corp | 電子機器の信頼性評価方法、その信頼性評価装置およびその信頼性評価プログラムを記録した媒体 |
WO2001030532A1 (en) * | 1999-10-27 | 2001-05-03 | Caterpillar Inc. | Method and apparatus for providing a simulation of a welding process using integrated models |
CN101122560A (zh) * | 2007-08-23 | 2008-02-13 | 南京航空航天大学 | 机械结构的裂纹扩展率和裂纹扩展寿命预测方法 |
JP4532689B2 (ja) * | 2000-07-21 | 2010-08-25 | 大成建設株式会社 | 凍結止水性の判定方法 |
CN102103658A (zh) * | 2010-12-03 | 2011-06-22 | 北京航空航天大学 | 一种实用的电子产品寿命评估模型参数高精度萃取方法 |
CN102169531A (zh) * | 2011-05-05 | 2011-08-31 | 北京航空航天大学 | 一种带热障涂层的圆管的热疲劳寿命预测方法 |
CN102967512A (zh) * | 2012-10-25 | 2013-03-13 | 中南大学 | 基于非对称循环应力控制加载的低周疲劳寿命预测方法 |
CN102980806A (zh) * | 2012-11-21 | 2013-03-20 | 中南大学 | 一种预测多步加载条件下金属材料低周疲劳寿命的方法 |
CN103778292A (zh) * | 2014-01-23 | 2014-05-07 | 北京航空航天大学 | 一种热振联合载荷下bga焊点疲劳寿命预测方法 |
CN103838975A (zh) * | 2014-03-18 | 2014-06-04 | 唐山轨道客车有限责任公司 | 车体焊缝疲劳寿命计算方法及装置 |
CN104268335A (zh) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | 微组装组件振动疲劳寿命预测方法和系统 |
CN104408325A (zh) * | 2014-12-12 | 2015-03-11 | 温州大学 | 一种微电子封装热膨胀系数可靠性匹配优化方法 |
-
2015
- 2015-05-06 CN CN201510226672.9A patent/CN104820781B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000046905A (ja) * | 1998-07-29 | 2000-02-18 | Sharp Corp | 電子機器の信頼性評価方法、その信頼性評価装置およびその信頼性評価プログラムを記録した媒体 |
WO2001030532A1 (en) * | 1999-10-27 | 2001-05-03 | Caterpillar Inc. | Method and apparatus for providing a simulation of a welding process using integrated models |
JP4532689B2 (ja) * | 2000-07-21 | 2010-08-25 | 大成建設株式会社 | 凍結止水性の判定方法 |
CN101122560A (zh) * | 2007-08-23 | 2008-02-13 | 南京航空航天大学 | 机械结构的裂纹扩展率和裂纹扩展寿命预测方法 |
CN102103658A (zh) * | 2010-12-03 | 2011-06-22 | 北京航空航天大学 | 一种实用的电子产品寿命评估模型参数高精度萃取方法 |
CN102169531A (zh) * | 2011-05-05 | 2011-08-31 | 北京航空航天大学 | 一种带热障涂层的圆管的热疲劳寿命预测方法 |
CN102967512A (zh) * | 2012-10-25 | 2013-03-13 | 中南大学 | 基于非对称循环应力控制加载的低周疲劳寿命预测方法 |
CN102980806A (zh) * | 2012-11-21 | 2013-03-20 | 中南大学 | 一种预测多步加载条件下金属材料低周疲劳寿命的方法 |
CN103778292A (zh) * | 2014-01-23 | 2014-05-07 | 北京航空航天大学 | 一种热振联合载荷下bga焊点疲劳寿命预测方法 |
CN103838975A (zh) * | 2014-03-18 | 2014-06-04 | 唐山轨道客车有限责任公司 | 车体焊缝疲劳寿命计算方法及装置 |
CN104268335A (zh) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | 微组装组件振动疲劳寿命预测方法和系统 |
CN104408325A (zh) * | 2014-12-12 | 2015-03-11 | 温州大学 | 一种微电子封装热膨胀系数可靠性匹配优化方法 |
Non-Patent Citations (9)
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156421B (zh) * | 2016-07-01 | 2019-04-05 | 电子科技大学 | 基于脉冲涡流热成像的电子封装焊点热疲劳寿命预测方法 |
CN106156421A (zh) * | 2016-07-01 | 2016-11-23 | 电子科技大学 | 基于脉冲涡流热成像的电子封装焊点热疲劳寿命预测方法 |
CN106568660A (zh) * | 2016-10-14 | 2017-04-19 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
CN106568660B (zh) * | 2016-10-14 | 2019-06-25 | 中国人民解放军海军航空工程学院青岛校区 | 一种复合材料胶接修补结构的剩余疲劳寿命的预测方法 |
CN108664669A (zh) * | 2017-03-27 | 2018-10-16 | 无锡市五十五度科技有限公司 | 一种bga焊点热寿命预测方法、测试平台及测试机箱 |
CN107203666A (zh) * | 2017-05-22 | 2017-09-26 | 北京航空航天大学 | 一种bga焊点热疲劳寿命的预测方法及系统 |
CN107203666B (zh) * | 2017-05-22 | 2019-06-21 | 北京航空航天大学 | 一种bga焊点热疲劳寿命的预测方法及系统 |
CN107273632B (zh) * | 2017-06-28 | 2019-09-27 | 北京航空航天大学 | 一种非累积冲击下多部件系统寿命预测方法及装置 |
CN107273632A (zh) * | 2017-06-28 | 2017-10-20 | 北京航空航天大学 | 一种非累积冲击下多部件系统寿命预测方法及装置 |
CN109387533A (zh) * | 2017-08-08 | 2019-02-26 | 北京航空航天大学 | 焊点在不规则温循剖面下损伤累积法则实验方法 |
CN108287978A (zh) * | 2018-03-26 | 2018-07-17 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 焊点设计参数的优化方法 |
CN109116140A (zh) * | 2018-07-16 | 2019-01-01 | 中国航空综合技术研究所 | 一种用于pbga封装器件的测试方法 |
CN110059412A (zh) * | 2019-04-19 | 2019-07-26 | 上海工程技术大学 | 金属材料寿命预测方法 |
CN110059412B (zh) * | 2019-04-19 | 2023-08-11 | 上海工程技术大学 | 金属材料寿命预测方法 |
CN111414709A (zh) * | 2020-04-14 | 2020-07-14 | 上海交通大学 | 一种基于神经网络的超级电容器循环寿命和容量衰减的早期预测方法 |
CN111521507A (zh) * | 2020-04-30 | 2020-08-11 | 江苏师范大学 | 一种环保不锈钢焊缝热疲劳表面裂纹扩展速率的试验方法 |
CN113189136A (zh) * | 2021-03-29 | 2021-07-30 | 北京航空航天大学 | 一种针对宇航用球栅阵列器件装联工艺的温度循环试验分级评价方法 |
CN113495189A (zh) * | 2021-09-08 | 2021-10-12 | 深圳荣耀智能机器有限公司 | 评估电子装联材料可靠性的测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104820781B (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104820781A (zh) | 考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法 | |
Kamal et al. | Advances in fatigue life modeling: A review | |
CN104242325B (zh) | 一种电力系统低频振荡模式参数辨识方法 | |
CN100533119C (zh) | 机械结构的裂纹扩展率和裂纹扩展寿命预测方法 | |
CN107991103A (zh) | 一种基于真实路谱的电动汽车电池包结构疲劳寿命预测方法 | |
CN104809311B (zh) | 一种基于多因素融合修正的结构件剩余寿命预测方法 | |
CN110414117A (zh) | 一种软包锂离子电池密封可靠度预测方法 | |
CN103955568A (zh) | 一种基于失效物理的mos器件可靠性仿真评价方法 | |
WO2020125682A1 (zh) | 一种应用随钻数据计算岩石强度的方法和系统 | |
CN105910902A (zh) | 混凝土构件裂缝扩展路径的分形分析方法 | |
RU2758181C1 (ru) | Способ оценки динамической надежности в отношении неисправностей соединения срединного желоба скребкового конвейера | |
CN110987676A (zh) | 随机多轴载荷下考虑裂纹闭合效应的全寿命预测方法 | |
CN113607580B (zh) | 一种金属构件疲劳试验方法及剩余寿命预测方法 | |
CN103364703B (zh) | Led多应力下可靠性快速评测方法 | |
CN106067032A (zh) | 用于数据异常监测的传感器选择方法的评价方法 | |
CN103353295B (zh) | 一种精确预测大坝坝体垂直变形量的方法 | |
CN117594164A (zh) | 基于数字孪生的金属结构剩余疲劳寿命计算评估方法、系统 | |
CN111460702A (zh) | 一种基于正逆向损伤特征融合的结构件损伤识别方法 | |
CN101520746A (zh) | 一种用于多种软件形态的质量评估方法及系统 | |
CN107609223B (zh) | 一种抗拉强度1200MPa的冷轧双相钢动态变形本构模型建立方法 | |
CN103399974B (zh) | 量化比较随机振动仿真数据与实验数据的方法 | |
CN110147643B (zh) | 车钩钩体剩余寿命确定方法和装置 | |
CN115408755B (zh) | 一种考虑时变效应的组合梁桥动力疲劳可靠性评估方法 | |
Kuczynska et al. | Effect of shear and tensile-dominant cyclic loading on failure in SnAgCu solder | |
CN105069214A (zh) | 一种基于非线性相关分析的工艺可靠性评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |