Stress analysis and optimization of POP stacked solder joints under thermal cyclic load
-
摘要: 建立叠层封装(packaging on packaging，POP)堆叠焊点有限元模型，基于ANAND本构方程，分析了热循环载荷下焊点应力分布状态及热疲劳寿命；基于灵敏度法分析了POP封装结构参数对焊点热应力的影响显著性；基于响应面法建立POP堆叠焊点热应力与结构参数的回归方程，并结合粒子群算法对结构参数进行了优化. 结果表明，焊点与铜焊盘接触处应力最大，该处会率先产生裂纹，上层焊点高度和下层焊点高度对POP堆叠焊点热应力影响较为显著；最优结构参数水平组合为上层焊点高度0.35 mm、下层焊点高度0.28 mm、中层印刷电路板厚度0.26 mm，优化后上、下两层焊点的最大热应力分别下降了0.816和1.271 MPa，延长了POP堆叠焊点热疲劳寿命.Abstract: The finite element model of packaging on packaging(POP) stacked solder joints was established, and the stress distribution state and thermal fatigue life of the solder joints under thermal cyclic load were analyzed based on the ANAND constitutive equation. The significance of the influence of the structural parameters of the POP package on the thermal stress of the solder joints was analyzed based on the sensitivity method. The regression equation of the thermal stress of the POP stacked solder joints and the structural parameters was established based on the response surface method, and the structural parameters were optimized by combining the particle swarm algorithm. The results show that the stress is highest at the contact between the solder joint and the copper solder disc, which is the location of the crack initiation. The upper and lower solder joint heights have a significant effect on the thermal stresses in POP stacked joints. The optimal combination of structural parameter levels is 0.35 mm for the upper solder joint height, 0.28 mm for the lower solder joint height and 0.26 mm for the middle PCB thickness, which reduces the maximum thermal stress of the upper and lower solder joints by 0.816 and 1.271 MPa, respectively, and prolongs the thermal fatigue life of POP stacked joints.
-
0. 序言
焊接是轨道车辆关键结构连接的主要工艺，作为一种新型熔化焊方法，激光填丝焊因具有变形小、焊缝成形好、缺陷少以及接头力学性能高等优点，近年来已备受车体结构焊接的青睐[1-2]. Q310NQL2-Q345NQR2耐候钢激光焊对接接头作为轨道车辆部分关键结构连接常用的接头形式，其在服役运行过程中经常承受复杂的动态载荷，进而引发突然性的疲劳断裂，因此其服役过程的疲劳失效一直是领域内关注的重点问题.
目前焊接接头的抗疲劳设计重点关注接头的疲劳极限和疲劳寿命评估[3-4]，传统的疲劳极限评估方法主要为升降法，而疲劳寿命评估则较多采用成组法完成，这些方法不仅历时长、成本高[5-11]，而且结果往往依赖于经验性的数据统计，因而在很大程度上限制了产品性能测试的速度以及研发效率的提升.疲劳实际上是材料内微观结构在循环载荷状态下朝着损伤累积和失效方向逐渐演化的过程，这一过程表现为耗散能的持续释放，同时伴随着热力学的温升响应变化[12]. 随着热像技术发展的日趋成熟，使得借助红外热像仪精确获取疲劳过程的三阶段温升响应成为可能，进而为构建物理意义更为明确的能量耗散模型奠定基础.
国内外学者以能量耗散为基础做了大量的疲劳研究工作，Fan等人[13-14]以Q235钢为研究对象，建立了一种可实现快速预测宏微观疲劳行为的能量耗散模型，基于此模型测定了Q235钢高周疲劳过程的能量耗散，并与其损伤变量建立了函数关系，进而对其疲劳极限进行了快速评估；Teng等人[15]结合RVE(representataive volume element)模型，进一步优化了Fan等人建立的能量耗散模型，并对SAE 1045钢的宏微观疲劳行为进行了快速预测评估; Yang等人[16-17]进一步研发了与微结构运动相关的能量耗散模型，该模型将滞弹性和非弹性行为的分界点作为疲劳极限预测的特征点，实现了FV520B钢疲劳极限的快速预测；Guo等人[18-20]基于一维热传导方程，推导了一种新的能量耗散计算方法，该方法以平均温升为基础参量，可最大程度上减小计算误差，为发展精确的能量耗散计算模型奠定了基础. 然而上述相关模型依赖的基础参量为第二阶段的温升增量，该参量往往需要大量的时间测得；此外上述研究主要针对非焊接金属材料展开，对于非均质焊接接头疲劳研究虽有涉及[5-11, 21-23]，但其能量耗散模型尚未与疲劳过程的微结构行为建立紧密联系.
为了实现快速、精确的焊接接头疲劳性能评估，以Q310NQL2-Q345NQR2激光焊对接接头为研究对象，以第一阶段的初始温升斜率为基础，结合ISV模型，以期建立一种与焊接接头高周疲劳过程微结构运动联系紧密的能量耗散模型，进而为实现从能量耗散角度出发快速评估疲劳性能奠定基础.
1. 基于ISV理论的本构模型
从热动力学的角度来讲，Hemholtz自由能是状态函数，其值可以通过获取当前的状态变量计算得到，考虑在热弹性范围内，Hemholtz自由能有如下定义[24]，即
$$ \psi = \psi ({\boldsymbol{\varepsilon }},T) $$ (1) 式中：
$\psi $ 为自由能；ε 和T分别为应变张量和绝对温度，是决定状态函数的观测状态变量. 从式(1)可以看出，在弹性尺度内，当ε 和T两个观测状态变量(observable state variable, OSV)确定时，其相关状态函数的值也随之确定. 当材料的变形行为更为复杂时，即其变形存在不完全弹性的情况，例如滞弹性和微塑性，往往必须扩展内部状态变量(ISV)实现独立定义现有的热力学状态函数[24].疲劳是一个伴随着温升演化的不可逆热力学过程, 为精确描述疲劳过程中的热力学状态变化，实现对包含非弹性的力学行为表征，这里引入α和β 两个ISV变量, 因此式(1)可以表示为[24]
$$ \psi (T,{\boldsymbol{\varepsilon }},{\boldsymbol{\alpha }},{\boldsymbol{\beta }}) = u - Ts $$ (2) 式中: T, ε, u和s分别为绝对温度、应变张量、内能以及熵流，此外ISV变量α 代表着与可恢复微结构运动相关的滞弹性行为，例如应力驱动的原子振荡和晶界处的内摩擦，而变量β 则表示不可恢复微结构运动相关的微塑性行为，进而引起材料内部的损伤(图1)，这里α 和β 可以定义为
$$ \left\{\begin{gathered} {\boldsymbol{\alpha }} = ({\alpha _1},{\alpha _2},\cdots,{\alpha _i},\cdots,{\alpha _m}) \\ {\boldsymbol{\beta }} = ({\beta _1},{\beta _2},\cdots,{\beta _j},\cdots,{\beta _n}) \\ \end{gathered} \right.$$ (3) 2. 能量耗散模型
根据热力学第一定律，疲劳过程的热平衡方程可以表示为[19]
$$ \rho C\dot T + {\text{div}}{{{J}}_{\text{q}}} = {d_1} + {\theta _{{\text{the}}}} + {\theta _{{\text{inc}}}} + {r_{\text{e}}} $$ (4) 式中: ρ和C分别为焊缝材料密度和比热；
${{J}}_{\text{q}} $ 为热流；d1，θthe，θinc和re分别表示为$$ {d_1}{\text{ = }}\left(\sigma {{ - }}\rho \frac{{\partial \psi }}{{\partial \varepsilon }}\right) \cdot {{\boldsymbol {\boldsymbol{\varepsilon}} }} - \rho \frac{{\partial \psi }}{{\partial \alpha }} \cdot {{ {\boldsymbol{\alpha}} }} - \rho \frac{{\partial \psi }}{{\partial \beta }} \cdot {{ {\boldsymbol{\beta}} }} $$ (5) $$ {\theta _{{\text{the}}}} = \rho T\frac{{{\partial ^2}\psi }}{{\partial T\partial \varepsilon }} \cdot {{\boldsymbol {\boldsymbol{\varepsilon}} }} $$ (6) $$ {\theta _{{\text{ic}}}} = \rho T\frac{{{\partial ^2}\psi }}{{\partial T\partial \alpha }} \cdot {{ {\boldsymbol{\alpha}} }} $$ (7) 根据文献[24]，可以推导得到
$$ \sigma = \rho \frac{{\partial \psi }}{{\partial \varepsilon }} $$ (8) 将式(8)代入式(5)，d1可以表示为
$$ {d_1}{\text{ = }} - \rho \frac{{\partial \psi }}{{\partial \alpha }} \cdot {{\boldsymbol \alpha }} - \rho \frac{{\partial \psi }}{{\partial \beta }} \cdot {{\boldsymbol \beta }} $$ (9) 式中：
$- \rho \dfrac{{\partial \psi }}{{\partial \alpha }} $ 和$- \rho \dfrac{{\partial \psi }}{{\partial \beta }} $ 为与内部变量α 和β 相关的热力学共轭矢量，其被认为是驱动材料内部达到热力学平衡状态的非平衡力，这里将其分别定义为A和B，得到$$ {d_1}{\text{ = }}{\boldsymbol{A}} \cdot {{\boldsymbol \alpha }} + {\boldsymbol{{ B}}} \cdot {{\boldsymbol \beta }} $$ (10) 式中：α 和β 为与微结构运动相关的内变量变化率，通常被视为广义热力学通量，由非平衡力A和B决定. 这里出于简化目的，假设广义热力学通量α 和β 与非平衡力A和B之间的关系服从线性或比例关系，这个假设对于正在经历高周疲劳过程的材料来说是非常合理的，因为其仅是在平衡状态附近波动，即接近一个准静态平衡过程[25]，即
$$ \left\{\begin{gathered} {{\boldsymbol \alpha }} = {m_{\rm{an}}}{\boldsymbol{A}} \\ {{\boldsymbol \beta }} = {m_{\rm{in}}}{\boldsymbol{B}} \\ \end{gathered} \right.$$ (11) 式中：
${m_{\rm{an}}}$ 和${m_{\rm{in}}}$ 是独立于热力学通量和非平衡驱动力的材料常数. 将式(11)带入式(10)，得到$$ {d_1} = {m_{\rm{an}}}{\left\| {\boldsymbol{A}} \right\|^2} + {m_{\rm{in}}}{\left\| {\boldsymbol{ B}} \right\|^2} $$ (12) 式中：||•||表示矢量的欧几里得范数. 从式(12)可以看出，d1由两部分组成，即一部分由滞弹性引导的能量耗散dan(不构成损伤)，而另外一部分则是与非弹性行为相关的能量耗散(引起疲劳损伤增加)din，式(12)变形为
$$ {d_1}{\text{ = }}{d_{{\text{an}}}}{\text{ + }}{d_{{\text{in}}}}{\text{ = }}{{m_{{\rm{an}}}}}{\left\| {\boldsymbol{A}} \right\|^2} + {{m_{{\rm{in}}}}}{\left\| {\boldsymbol{ B}} \right\|^2} $$ (13) 对于疲劳过程来说，与材料微结构运动相关的非平衡力A和B与外加载荷以及频率有关，因此这里假设||A||2和||B||2与应力幅和频率有如下关系，即
$$ \left\{ \begin{gathered} {\left\| {\boldsymbol{A}} \right\|^2} = \left[ {{u_{{\text{an}}}}({\sigma _{\text{a}}} - {\sigma _{\text{0}}})} \right]\left[ {H({\sigma _{\text{a}}} - {\sigma _{{\text{c0}}}})} \right] \\ {\left\| {\boldsymbol{B}} \right\|^2} = ({u_{{\text{in}}}}{\sigma _{\text{a}}}^kf)\left[ {H({\sigma _{\text{a}}} - {\sigma _{{\text{c1}}}})} \right] \\ \end{gathered}\right. $$ (14) 式中：H(x)为海维赛德阶跃函数(heaviside step function)，即当x≥0时，H(x)=1，当x<0时，H(x)=0；k为幂指数；uan和uin为比例系数；σc0和σc1为与材料微结构运动状态相关的两个特征应力幅值；
${\sigma _{\text{a}}} $ 为应力幅. 将式(14)代入式(13)，得到$$ {d_1} = \left\{ \begin{array}{l} {d_{{\rm{an}}}} = \left[ {{F_{{\rm{an}}}}({\sigma _{\rm{a}}} - {\sigma _{{\rm{c0}}}})} \right],{\sigma _{{\rm{c0}}}} < {\sigma _{\rm{a}}} < {\sigma _{{\rm{c1}}}}\\ {d_{{\rm{an}}}} + {d_{{\rm{in}}}} = \left[ {{F_{{\rm{an}}}}({\sigma _{\rm{a}}} - {\sigma _{{\rm{c0}}}})} \right] + {F_{{\rm{in}}}}{\sigma _{\rm{a}}}^k,{\sigma _{{\rm{c1}}}} < {\sigma _{\rm{a}}} \end{array} \right. $$ (15) 式中：Fan(=man uan f )，Fin(=min uin f )，k，σc0以及σc1可以通过试验测定.
3. 疲劳寿命预测模型
3.1 能量耗散响应
疲劳过程的温升进程由3个阶段组成[5]，分别为初始加载的温度迅速增加阶段(Ⅰ)、热产和热耗散实现平衡的稳定温升阶段(Ⅱ)以及因裂纹扩展引起的温升陡然变化阶段(Ⅲ)，如图2所示.
对于高周疲劳过程，其某一应力水平下的能量耗散计算式为[26]
$$ {d_1}{\text{ = }}{\left. {\rho C\frac{{{\text{d}}T}}{{{\text{d}}t}}} \right|_{t = {t_0}}} $$ (16) 式中：
${\dfrac{{\rm{d}}T}{{\rm{d}}t} }$ 为温度变化率；t0为初始加载时间.如第2节所述，高周疲劳过程的机械能耗散主要由滞弹性和非弹性行为引起[19]，其中当应力幅小于焊接接头的疲劳极限时，其微结构运动主要为可恢复的滞弹性行为，例如轻微的原子振荡或内摩擦，不构成损伤累积；而当应力幅高于焊接接头的疲劳极限时，其微结构运动主要为不可恢复的非弹性行为，例如微塑性，构成接头内部的损伤[5]. 基于此，将疲劳过程中对接接头在不同应力幅水平下微结构运动的分界点定义为疲劳极限，具有明确的物理意义.
为了从宏观热力学角度量化微观尺度内的接头微结构运动行为，在第1节里引入了ISV内变量理论，并分别定义了与滞弹性和非弹性相关的内变量α 和β，进而推导得到了高周疲劳过程的能量耗散计算表达式. 结合第2节建立的能量耗散模型，疲劳过程中不同应力幅水平下的能量耗散响应示意图如图3所示.
![]()
图 3 不同应力幅水平下的能量耗散响应Figure 3. Energy dissipation response under different stress amplitude levels3.2 基于损伤累积的疲劳寿命预测模型
如第1节所述，ISV里的内变量β与引导疲劳损伤累积的不可恢复微结构运动相关，而该疲劳损伤的演化正由内部非平衡力B驱动，并在宏观尺度上表现为非弹性耗散din，因此非弹性耗散din可被视为有效的疲劳损伤指标. 基于此可建立相应的累积损伤评估模型，即
$$ D = \frac{{N{F_{{\rm{in}}}}\sigma _{\text{a}}^k}}{{{{\bar E}_{\text{c}}}}}\left[ {H({\sigma _{\text{a}}} - {\sigma _{{\text{c}}1}})} \right] $$ (17) 式中：D为损伤参数；N为与损伤参量D对应的疲劳寿命；
$\bar {E}_{\rm{c}} $ 为能量耗散临界值，其值可以通过试验测定. 对于一个原始试样来说，其损伤值D可认为是0，当载荷逐渐循环直至疲劳断裂出时，D=1，因此式(17)变形为$$ {N_{\text{f}}} = \frac{{{{\bar E}_{\text{c}}}}}{{{F_{{\rm{in}}}}\sigma _{\text{a}}^k}}\left[ {{H^{ - 1}}({\sigma _{\text{a}}} - {\sigma _{{\text{c}}1}})} \right] $$ (18) 式中：H −1(x)为海维赛德阶跃函数的倒数，即当x≥0时，H −1(x)=1，而当x<0时，H −1(x)=∞. 这里需要注意的是，当σa<σc1时，认为对接接头试样拥有无限寿命，因此依据提出的模型，关键应力幅σc1实际上可以被认为是材料内部可恢复和不可恢复微结构运动的分界应力幅，即σc1等于疲劳极限σf.
通过监测试样焊缝中心的温升增量，建立能量耗散模型，进而结合式(18)即可实现Q310NQL2-Q345NQR2激光焊对接接头的疲劳寿命快速预测.
4. 试验方法
4.1 试验材料
采用激光填丝焊对2.5 mm板厚的Q310NQL2耐候钢和4 mm板厚的Q345NQR2耐候钢进行焊接，激光器功率为6 kW的光纤激光器，焊丝为ER50-G，送丝设备选用数字化送丝结构. 焊接过程中，填充焊丝位于激光之前，激光与试板垂线的夹角为10°，保护气体为氩气，坡口形式为I形(全焊透)，激光焊接过程示意图如图4所示. Q310NQL2耐候钢、Q345NQR2耐候钢、ER50-G焊丝的化学成分以及力学性能见表1和表2所示. Q310NQL2-Q345NQR2激光焊对接接头尺寸如图5所示.
表 1 Q310NQL2和Q345NQR2耐候钢以及ER50-G焊丝主要化学成分Table 1. Main chemical composition of Q310NQL2, Q345NQR2 weathering steel, and ER50-G filler wire材料 C Si Mn P S Cu Cr Ni Ti Fe Q310NQL2 ≤0.12 0.25 ~ 0.75 0.20 ~ 0.50 0.06 ~ 0.12 ≤0.02 0.25 ~ 0.50 0.30 ~ 1.25 0.12 ~ 0.65 — 余量 Q345NQR2 ≤0.12 0.25 ~ 0.75 ≤1.0 0.06 ~ 0.15 ≤0.02 0.25 ~ 0.50 0.30 ~ 1.25 0.12 ~ 0.65 — 余量 ER50-G ≤0.10 ≤0.60 0.90 ~ 1.30 ≤0.025 ≤0.02 0.20 ~ 0.50 0.30 ~ 0.90 0.20 ~ 0.60 — 余量 表 2 Q310NQL2和Q345NQR2耐候钢以及ER50-G焊丝力学性能Table 2. Mechanical propertie parameters of Q310NQL2, Q345NQR2 weathering steel, and ER50-G filler wire材料 屈服强度
ReL/MPa抗拉强度
Rm/MPa断后伸长率
A(%)Q310NQL2 ≥310 480 ~ 670 ≥22 Q345NQR2 ≥345 490 ~ 675 ≥22 ER-50G ≥400 ≥500 ≥22 4.2 疲劳试验与热像测试
根据国家标准GB/T 3075—2008《金属材料 疲劳试验 轴向力控制方法》在PLG-100型疲劳试验机上对Q310NQL2-Q345NQR2对接接头进行疲劳测试(图6)，试验频率为106 Hz，载荷形式采用单轴、正弦变化的动态载荷，应力比(R=σmin/σmax)取0.1. 选用Fluke Ti450红外热像仪对焊缝中心疲劳过程的温升响应进行全程监测，其采集图像的频率为9 Hz. 疲劳试验热像测试具体试验步骤如下.
(1) 借助热像仪监测对接接头焊缝中心在不同应力幅水平下(σa=99.5~171 MPa)的初始阶段温升斜率，结合式(16)计算不同应力幅水平的能量耗散.
(2) 结合建立的能量耗散模型，计算得到不同应力幅下的能量耗散响应，进而测定Fan , Fin，k，σc0以及σc1(疲劳极限)等参数.
(3) 进行恒幅疲劳测试获取能量耗散临界值
${\bar E}_{\rm{c}} $ .5. 试验结果
5.1 不同应力幅下的能量耗散响应分析
根据第2节建立的能量耗散模型，结合式(16)计算了不同应力幅下的能量耗散值，如图7所示. 结合式(15)，采用最小二乘法对试验数据实现了较好的拟合，并对σc0，σc1，Fan，Fin以及k等参数进行了确定(表3).
![]()
图 7 Q310NQL2-Q345NQR2对接接头在不同应力幅下的能量耗散响应Figure 7. Energy dissipation response of Q310NQL2-Q345NQR2 butt joints under different stress amplitudes表 3 能量耗散模型参数Table 3. Parameters of energy disspation modelσc0 σc1 Fan Fin k 65 126 4.7 × 10−2 8.20 × 10−23 10.18 从图7中可以看出，随着应力水平的增加，能量耗散值随之增加，且在第一类关键应力(σc0=65 MPa)和第二类关键应力(σc1=126 MPa，即疲劳极限)之间，其能量耗散与应力幅的依赖程度几乎是线性的；而当应力幅值超过第二类关键应力σc1时，其不同应力幅下的能量耗散响应关系呈现为非线性趋势，这与前述模型所述一致(图3). 结合建立的能量耗散模型(式(15))对试验数据拟合，拟合的决策系数R2为0.97，这表明能量耗散模型与试验数据的一致性较高，这种较高的一致性源于两种不同的微结构演化机制，如3.1节所述，当应力幅等级介于σc0和σc1之间，其微结构运动主要为可恢复的滞弹性行为，例如轻微的原子振荡或内摩擦，因此在该阶段其能量耗散的响应几乎是线性的；而当应力幅等级高于σc1(即疲劳极限)时，由于不可逆微塑性行为的出现，且占该阶段微结构运动的绝大部分，使得其能量耗散响应急剧增加，且呈现为非线性的特点.
$${d_1}{\rm{ = }}\left\{ \begin{array}{l} {d_{{\rm{an}}}} = \left[ {0.047 ({\sigma _{\rm{a}}} - 65)} \right],\quad 65 < {\sigma _{\rm{a}}} < 126\\ {d_{{\rm{an}}}}{\rm{ + }}{d_{{\rm{in}}}} = \left[ {0.047 ({\sigma _{\rm{a}}} - 65)} \right]+\\ \qquad\qquad 8.20 \times {10^{ - 23}} {\sigma _{\rm{a}}}^{10.18},\quad {\sigma _{\rm{a}}} > 126 \end{array} \right.$$ (19) 5.2 基于能量耗散的疲劳寿命预测
如3.3节所述，疲劳过程中的能量耗散存在临界值Ec, 因此结合图3和式(19)，对应力幅为148.5 MPa下与不可逆微结构运动相关的能量耗散临界值进行了计算，可确定
${\bar E}_{\rm{c}}$ 为1.35 × 105 J/m3. 从图8中可以看出，试样的不可逆能量耗散值基本分布在平均值加减一个标准差的范围内，这表明该临界值几乎是一个常数，与文献报道的结果一致[17].![]()
图 8 应力幅为148.5 MPa下的对接接头高周疲劳能量耗散及其平均值Figure 8. Energy dissipation and its mean value of butt joints with the stress amplitude equals 148.5 MPa基于试验测定的能量耗散临界值
${\bar E}_{\rm{c}} $ ，结合3.3节建立能量耗散模型(式(18))对Q310NQL2-Q345NQR2对接接头的疲劳寿命进行了快速预测，如图9所示，拟合方法为最小二乘法. 从图9中可以看出，两条拟合S-N曲线非常接近，其中经预测数据拟合的S-N曲线斜率和截距为−10.18和27.22，而经试验数据拟合的S-N曲线的斜率和截距为−10.78和28.52，其相对误差值分别为5.54%和4.77%，均小于10%，从而验证了提出模型具有较高的精度. 因此在试样有限的情况下，结合建立的能量耗散模型，获取其不同应力等级下的能量耗散响应情况，继而可实现Q310NQL2-Q345NQR2对接接头的疲劳寿命快速预测评估. 值得指出的是，尽管该模型对于Q310NQL2-Q345NQR2对接接头的疲劳性能具有较好的预测能力，但是其对于其它类型接头的适用性还需要在接下来的工作中进一步验证.![]()
图 9 经预测数据和试验数据拟合的中值S-N曲线对比Figure 9. Comparison of the median S-N curve fitted by the predicted and tested data6. 结论
(1)结合ISV理论的本构模型，提出了一种新的基于能量耗散理论的对接接头高周疲劳寿命快速预测模型.
(2)模型引入了两个关键应力幅，并将其与可恢复微结构运动(滞弹性)和不可恢复微结构运动(非弹性)开始形成的关键应力幅σc0和σc1(疲劳极限)建立了紧密关系，阐明了疲劳过程中的微结构运动的物理意义.
(3)基于建立的疲劳寿命预测模型，对Q310NQL2-Q345NQR2对接接头的疲劳进行了快速预测，并借助预测数据对其S-N曲线进行了拟合，结果表明，经预测数据和试验数据拟合的中值S-N曲线非常接近，其斜率和截距误差仅为5.54%和4.77%，从而验证了提出模型的有效性.
-
![]()
图 4 上层和下层焊点热应力分布
Figure 4. Thermal stress distribution of upper and lower solder joints. (a) upper solder joints; (b) lower solder joints
![]()
图 5 上层焊点的散点图
Figure 5. Scatter diagram of upper solder joints. (a) thickness of intermediate PCB; (b)thickness of top PCB; (c) height of lower solder joint; (d) height of upper solder joint
![]()
图 7 下层焊点的散点图
Figure 7. Scatter diagram of lower solder joints. (a) thickness of intermediate PCB; (b) thickness of top PCB; (c) height of lower solder joint; (d) height of upper solder joint
![]()
图 10 最优参数组合上层和下层焊点热应力分布
Figure 10. Thermal stress distribution diagram of upper and lower solder joints with optimal parameter combination. (a) upper solder joints; (b) lower solder joints
表 1 POP堆叠焊点尺寸
Table 1 POP stacked solder joint dimensions mm
上层焊点 下层焊点 高度h1 直径d1 间距δ1 高度h2 直径d2 间距δ2 0.30 0.40 0.65 0.23 0.30 0.50 
下载: 导出CSV
表 2 材料参数
Table 2 Material parameters
材料 弹性模量E/GPa 泊松比μ 热膨胀系数α/10−6K−1 SAC305 38.7 − 16.9T 0.35 1.25 芯片 163 0.28 2.5 PCB板 300 0.20 3.5 铜焊盘 117 0.30 14.3
下载: 导出CSV
表 3 SAC305钎料ANAND模型参数
Table 3 SAC305 filler metal ANAND model parameters
初始形变阻抗值 $ {{{S}}_{ \rm{o}}} $ /MPa应力乘子ξ 常数 $A$ /104s−1变形阻力饱和值系数 $ \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{S} $ /MPa激活能常数(Q/R)/K−1 硬化/软化常数 $ {{{h}}_{{{\rm{o}}}}} $ /MPa应变率敏感性指数m 饱和值 应力 软硬化 45.9 2 587 58.3 7460 9350 0.0942 1.5 0.015
下载: 导出CSV
表 4 焊点粘塑性应变能密度增量
Table 4 Viscoplastic strain energy density increment of solder joints J/m3
焊点位置 应变能密度增量ΔWave 上层 82 822.671，83 134.183，83 208.574 下层 26 745.566，26 810.497，26 863.966
下载: 导出CSV
表 5 封装体结构参数取值
Table 5 Value of package structure parameters mm
上层焊点高度HS2 下层焊点高度HS1 中层PCB厚度δP2 顶层PCB厚度δP3 0.25 ~ 0.35 0.18 ~ 0.28 0.26 ~ 0.34 0.10 ~ 0.16
下载: 导出CSV
表 6 因素水平表
Table 6 Factor level table
水平 上层焊点高度HS2 /mm 下层焊点高度HS1 /mm 中层PCB厚度δP2 /mm −1 0.25 0.18 0.26 0 0.30 0.23 0.30 1 0.35 0.28 0.34
下载: 导出CSV
表 7 响应面设计组合与热应力分析结果
Table 7 Response Surface design combination and thermal stress analysis results
试验编号 上层焊点高度HS2 /mm 下层焊点高度HS1 /mm 中层PCB厚度δP2 /mm 下层焊点应力y1 /MPa 上层焊点应力y2 /MPa 1 0.25 0.18 0.30 45.726 41.402 2 0.30 0.23 0.30 43.908 40.422 3 0.35 0.18 0.30 45.731 39.720 4 0.35 0.23 0.26 43.876 39.629 5 0.25 0.23 0.26 43.887 41.488 6 0.35 0.28 0.30 42.661 39.681 7 0.30 0.18 0.26 45.453 40.380 8 0.30 0.28 0.26 42.642 40.363 9 0.30 0.23 0.30 43.908 40.422 10 0.30 0.23 0.30 43.908 40.422 11 0.30 0.18 0.34 45.516 40.416 12 0.30 0.28 0.34 42.682 40.451 13 0.30 0.23 0.30 43.908 40.422 14 0.25 0.23 0.34 43.931 41.311 15 0.25 0.28 0.30 42.669 41.356 16 0.35 0.23 0.34 43.923 39.755 17 0.30 0.23 0.30 43.908 40.422
下载: 导出CSV
表 8 响应面分析结果
Table 8 Response surface analysis results
焊点 P值 回归方程系数R2 回归方程调整系数Ra2 回归方程预测系数Rp2 上层 < 0.0001 0.9990 0.9978 0.9847 下层 < 0.0001 0.9983 0.9962 0.9736
下载: 导出CSV
-
[1] 弗林·卡森. 封装体叠层(PoP)技术及其应用[J]. 集成电路应用, 2015(2): 38 − 41. doi: 10.3969/j.issn.1674-2583.2015.02.012 Flynn Carson. Package-on-Package (PoP) technology and its applications[J]. Application of IC, 2015(2): 38 − 41. doi: 10.3969/j.issn.1674-2583.2015.02.012
[2] Honarvar S, Nourani A, Karimi M. Effect of thermal treatment on fracture behavior of solder joints at various strain rates: Comparison of cyclic and constant temperature[J]. Engineering Failure Analysis, 2021, 128(8): 105636.
[3] 韦何耕, 黄春跃, 梁颖, 等. 热循环加载条件下PBGA叠层无铅焊点可靠性分析[J]. 焊接学报, 2013, 34(10): 91 − 94. Wei Hegeng, Huang Chunyue, Liang Ying, et al. Reliability analysis of plastic ball grid array double-bump lead-free solder joint under thermal cycle[J]. Transactions of the China Welding Institution, 2013, 34(10): 91 − 94.
[4] 王洋. 不同温度环境下POP堆叠封装的可靠性研究[D]. 镇江: 江苏大学, 2016. Wang Yang. Research on the reliability of package-on-package in different temperature environment[D]. Zhenjiang: Jiangsu University, 2016.
[5] 孙磊, 张屹, 陈明和, 等. 3D封装芯片焊点可靠性有限元分析[J]. 焊接学报, 2021, 42(1): 49 − 53. doi: 10.12073/j.hjxb.20201021002 Sun Lei, Zhang Yi, Chen Minghe, et al. Finite element analysis of solder joint reliability of 3D packaging chip[J]. Transactions of the China Welding Institution, 2021, 42(1): 49 − 53. doi: 10.12073/j.hjxb.20201021002
[6] 王红霞, 王蓓, 冉红锋, 等. TMV POP热翘曲变形及可靠性[J]. 半导体技术, 2018, 43(6): 462 − 467. Wang Hongxia, Wang Bei, Ran Hongfeng, et al. Warpage behaviors and reliability in TMV POP[J]. Semiconductor Technology, 2018, 43(6): 462 − 467.
[7] 刘昭云, 杨雪霞. 热循环载荷下POP堆叠电子封装的可靠性研究[J]. 太原科技大学学报, 2018, 39(6): 473 − 478. doi: 10.3969/j.issn.1673-2057.2018.06.011 Liu Zhaoyun, Yang Xuexia. Research on reliability of POP stacked electronic package under thermal cycling load[J]. Journal of Taiyuan University of Science and Technology, 2018, 39(6): 473 − 478. doi: 10.3969/j.issn.1673-2057.2018.06.011
[8] Shen Y A, Zhou S, Li J, et al. Sn-3.0Ag-0.5Cu/Sn-58Bi composite solder joint assembled using a low-temperature reflow process for PoP technology[J]. Materials & Design, 2019, 183: 108144.
[9] Sitek J, Koscielski M, Borecki J, et al. The dependence of reliability and mechanical strength of the solder joints in 3D PoP structures from sizes of solder powders applied in soldering materials[J]. Soldering & Surface Mount Technology, 2017, 29(1): 23 − 27.
[10] Chen L , Lai X , Xiao Y , et al. Analysis of POP solder ball thermal cycling fatigue life based on stress strain[C]// 2012 13th International Conference on Electronic Packaging Technology and High Density Packaging. IEEE, Guilin, China, 2012.
[11] Li D B, Wang J P, Yang B, et al. Thermal performance and fatigue life prediction of POP stacked chip assembly under thermal cycling load[J]. Microelectronics International, 2020, 34(7): 165 − 171.
[12] 田茹玉, 王晨曦, 田艳红, 等. 极限温度下CBGA焊点热冲击疲劳寿命预测[J]. 焊接学报, 2017, 38(10): 93 − 97. doi: 10.12073/j.hjxb.20151208005 Tian Ruyu, Wang Chenxi, Tian Yanhong, et al. Life prediction of CBGA soldered joints under extreme temperature thermal shock[J]. Transactions of the China Welding Institution, 2017, 38(10): 93 − 97. doi: 10.12073/j.hjxb.20151208005
[13] Zhang L, Han J G, Guo Y H, et al. Anand model and FEM analysis of SnAgCuZn lead-free solder joints in wafer level chip scale packaging devices[J]. Microelectronics Reliability, 2014, 54(1): 281 − 286. doi: 10.1016/j.microrel.2013.07.100
[14] Darveaux R. Effect of simulation methodology on solder joint crack growth correlation and fatigue life prediction[J]. Journal of Electronic Packaging, 2002, 124(3): 147 − 154. doi: 10.1115/1.1413764
[15] 赵胜军, 黄春跃, 梁颖, 等. BGA无铅焊点再流焊焊后残余应力分析与优化[J]. 机械工程学报, 2020, 56(10): 86 − 94. doi: 10.3901/JME.2020.10.086 Zhao Shengjun, Huang Chunyue, Liang Ying, et al. Analysis and optimization of residual stress after reflow soldering of BGA lead free solder joints[J]. Journal of Mechnical Engineering, 2020, 56(10): 86 − 94. doi: 10.3901/JME.2020.10.086
[16] 关挺, 钟绍华, 周才, 等. 基于响应面和粒子群算法的悬置多目标优化[J]. 物流技术, 2017(10): 107 − 114. doi: 10.3969/j.issn.1005-152X.2017.10.021 Guan Ting, Zhong Shaohua, Zhou Cai, et al. Multi-objective optimization of Cab suspension based on response surface and PSO[J]. Logistics Technology, 2017(10): 107 − 114. doi: 10.3969/j.issn.1005-152X.2017.10.021
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 292
- HTML全文浏览量: 54
- PDF下载量: 73
- 被引次数: 1
