焊点焊后残余应力预测与优化的变可信度模型确定方法

1.本发明涉及电子元器件封装技术领域，具体涉及焊点焊后残余应力预测与优化的变可信度模型确定方法。


背景技术：

2.随着电子产品逐步向高性能化、高集成化、微小型化的趋势发展，电子封装得到更为广泛的应用，并成为电子信息产业的核心技术之一。其中球栅阵列封装(ball grid array package，bga)因集成度高且性能优异，成为目前电子封装中使用最多的技术之一。bga焊点在球栅阵列封装中起着电气连接、机械支撑和散热等作用，同样也是高密度电子封装中失效率最高的部位。然而，由于bga封装中各组件热膨胀系数不同，在回流焊冷却过程中，焊点内会产生热应力，即焊后残余应力。残余应力过大极易导致出现焊点断裂等缺陷，使得焊点连接可靠性下降，增加服役过程中失效风险。因此，bga焊点残余应力的预测分析与优化，对于降低生产成本和产品缺陷率来说尤为重要。
3.目前bga焊点焊后残余应力存在实测困难、成本较高的问题。上述问题亟待解决，为此，提出焊点焊后残余应力预测与优化的变可信度模型确定方法。


技术实现要素：

4.本发明所要解决的技术问题在于：如何解决bga焊点焊后残余应力实测困难、成本较高的问题，提供了焊点焊后残余应力预测与优化的变可信度模型确定方法，通过对回流焊工艺过程进行仿真，获得焊后的最大残余应力，据此建立焊点高度和焊盘半径等尺寸参数与焊后最大残余应力间的变可信度近似模型，采取嵌套和非嵌套的采样方式组成四种变可信度近似模型，并进行预测效果对比，同时使用粒子群优化算法进行焊点尺寸参数的优化，结果验证了变可信度近似模型的预测效果更优并获得其最佳构建方案。
5.本发明是通过以下技术方案解决上述技术问题的，本发明包括以下步骤：
6.s1：将回流焊过程分为焊点熔化成型和冷却凝固两个阶段，选择对焊后残余应力影响显著的焊点尺寸参数，即焊点高度h和焊盘半径r作为输入，在焊点融化成型阶段，基于交互式有限元分析软件surface evolver模拟单个bga焊点三维形态；
7.s2：将模拟的bga焊点形态转化为相关文件，基于三维数字化软件solidworks构建bga焊点有限元仿真模型；
8.s3：基于数值模拟综合分析软件ansys workbench进行焊后残余应力分析，采用间接耦合方法，以典型回流曲线为边界条件边界条件，对回流焊过程进行仿真计算，获得焊后最大残余应力；
9.s4：对两个阶段分别使用增加细化次数的方式，并使用edge sizing模块细化网格，进行有限元分析的收敛性检验，确定高/低精度样本点的划分方式；
10.s5：使用拉丁超立方采样方法(latin hypercube sampling,lhs)，获取嵌套与非嵌套的高/低精度样本点；
11.s6：根据高/低精度样本点、划分方式以及样本点嵌套/非嵌套关系，构建四种co-kriging变可信度近似模型，同时构建rbf单可信度近似模型；
12.s7：将co-kriging变可信度近似模型和rbf单可信度近似模型进行焊后残余应力预测效果的验证与对比；
13.s8：使用粒子群优化算法，优化焊点尺寸参数方案，以获取bga焊点焊后残余应力的最小值；
14.s9：分别基于co-kriging变可信度近似模型和rbf单可信度近似模型进行全局预测，进行寻优，分别获得最优焊点尺寸方案，并对结果仿真验证，结果表明co-kriging变可信度近似模型的预测性能更优。
15.更进一步地，在所述步骤s1中，基于最小能量原理和有限元数值分析方法，定义初始几何条件、能量形式、约束条件、材料性质参数，采用surface evolver软件模拟bga焊点成型三维形态。
16.更进一步地，在所述步骤s2中，构建bga焊点有限元仿真模型的具体过程如下：
17.s21：利用surface evolver软件使用有限元方法模拟并输出焊点模型的点、面参数文件，将焊点表面分为大量有限的小单元和节点，输出单个焊点的节点坐标参数文本文件，并将其转化为ascii码格式的面片文件，采用solidworks软件实体化焊点模型并装配成为bga组件，转为“x_t”格式文件；
18.s22：根据“x_t”格式文件，采用局部细化edge sizing方法划分网格，基于ansys软件建立bga焊点有限元仿真模型。
19.更进一步地，在所述步骤s3中，采用ansys workbench有限元综合分析软件，使用间接耦合方法，使用瞬态热学模块对bga组件施加温度场，将热分析的温度场结果作为载荷，并载入结构分析中进行热结构耦合分析，热分析过程的边界条件为典型回流曲线，结构分析过程的边界条件和仿真结束时刻模拟真实环境设置。
20.更进一步地，在所述步骤s4中，在两个阶段分别使用增加份数的方式细化网格，随着网格密度的增加，最大残余应力随着网格的细化逐渐趋近于一确定值；其中，对焊点成型模拟过程，利用surface evolver增加细化次数以提高单元密度，使得单元数目呈指数增长，并以细化次数作为收敛性检验的输入参数，最大残余应力为输出参数，得到最终结果；对焊点残余应力分析，仅对焊点阵列四角的焊点进行细化，在ansys workbench的mesh模块内选择edge sizing的细化方法并按份数划分网格，以份数作为收敛性检验的输入参数，最大残余应力为输出参数，得到最终结果。
21.更进一步地，在所述步骤s5中，根据co-kriging变可信度近似模型对高/低精度样本点关系的限定，分为高/低精度样本点非嵌套、嵌套两种，其中对于非嵌套样本，直接用lhs方法分别生成高/低精度样本点；对于嵌套样本，先用传统试验设计方法独立生成高/低精度样本点，然后对高精度样本点平移，以满足嵌套关系。
22.更进一步地，在所述步骤s6中，分别基于所获得的低精度样本数据和高/低精度数据差值构建相应的kriging近似模型，并通过对样本值的最大似然估计求解co-kriging变可信度近似模型的模型参数值，分别为低精度样本点构建的kriging近似模型的重要度参数θc、高/低精度数据差值构建相应的kriging近似模型的重要度参数θd、低精度样本点构建的kriging近似模型的基函数指数pc、高/低精度数据差值构建相应的kriging近似模型的
基函数指数pd、缩放参数ρ，从而确定co-kriging变可信度近似模型的预测估计表达式，其中，构建的co-kriging变可信度近似模型表示如下：
23.ze(x)＝ρzc(x)+zd(x)
24.其中，x为样本点，zc表示低精度样本数据构建的kriging近似模型，zd表示高/低精度数据差值构建的kriging近似模型；
25.对于独立的高/低精度两组数据，高精度数据有ne个样本xe，低精度数据有nc个样本xc，用rc(xc，xe)表示样本数据xc与xe间的相关性矩阵，σc和σd分别表示zc和zd的方差，l表示值全为1的列向量；
26.co-kriging变可信度近似模型的协方差矩阵表示为：
[0027][0028]
co-kriging变可信度近似模型的预测估计表达式为：
[0029][0030]
更进一步地，在所述步骤s7中，选取最大相对误差mre和均方相对误差mrse作为近似模型精度的评价标准，对四种co-kriging变可信度近似模型和rbf单可信度近似模型进行焊后残余应力预测效果的验证与对比；其中，mrse和mre两个指标表达式分别如下：
[0031][0032]
其中，yi为各点真实值，为近似模型预测值。
[0033]
更进一步地，在所述步骤s8中，粒子群优化算法参数设定后，用co-kriging变可信度近似模型进行全局预测，粒子群优化算法产生随机初始粒子和每次迭代产生的新粒子，循环直到满足迭代次数要求，最终可以获得最优解，即最小残余应力。
[0034]
更进一步地，在所述步骤s9中，对焊点设计空间进行网格划分，设置具体参数，分别利用co-kriging变可信度近似模型和rbf单可信度近似模型进行全局预测，并利用粒子群优化算法进行寻优，分别得到优化后的最优焊点尺寸方案；对两组优化结果进行仿真验证，得到co-kriging变可信度近似模型的结果与仿真值误差更小，即验证co-kriging变可信度近似模型预测性能更优。
[0035]
本发明相比现有技术具有以下优点：
[0036]
(1)发现网格粗细划分在焊点形态模拟阶段的非嵌套高/低精度样本点构建的变可信度近似模型预测性能最佳，且基于变可信度近似模型的全局预测更准确，为焊点形态的设计以及焊后残余应力地预测分析与优化提供了一定理论依据。
[0037]
(2)避免了传统设计方法中凭经验设计，大量试验验证工作，能缩短设计周期，减少研发成本。
[0038]
(3)实现了在节约成本的条件下比单精度模型更好的目标效果，兼顾成本与精度、
同时具有较好的优化收敛性，极大地方便了bga焊点焊后的残余应力预测分析与优化。
附图说明
[0039]
图1是本发明实施例二中焊点焊后残余应力预测与优化的变可信度模型确定方法的流程示意图；
[0040]
图2是本发明实施例二中回流焊温度曲线及焊球状态示意图；
[0041]
图3a是本发明实施例二中bga焊点有限元仿真模型中四角焊点网格划分示意图；
[0042]
图3b是本发明实施例二中bga焊点有限元仿真模型中内部焊点网格划分示意图；
[0043]
图3c是本发明实施例二中bga焊点有限元仿真模型中bga组件网格划分示意图；
[0044]
图4是本发明实施例二中回流焊冷却区曲线图；
[0045]
图5是本发明实施例二中焊后温度分布图；
[0046]
图6是本发明实施例二中bga焊点残余应力分布图；
[0047]
图7是本发明实施例二中最大残余应力与网格份数的关系曲线图；
[0048]
图8a是本发明实施例二中非嵌套样本点示意图；
[0049]
图8b是本发明实施例二中嵌套样本点示意图；
[0050]
图9是本发明实施例二中样本点与四种co-kriging变可信度近似模型的关系图；
[0051]
图10a是本发明实施例二中co-kriging变可信度近似模型model1在设计空间的映射图；
[0052]
图10b是本发明实施例二中co-kriging变可信度近似模型model1在设计空间的三维图；
[0053]
图10c是本发明实施例二中co-kriging变可信度近似模型model2在设计空间的映射图；
[0054]
图10d是本发明实施例二中co-kriging变可信度近似模型model2在设计空间的三维图；
[0055]
图10e是本发明实施例二中co-kriging变可信度近似模型model3在设计空间的映射图；
[0056]
图10f是本发明实施例二中co-kriging变可信度近似模型model3在设计空间的三维图；
[0057]
图10g是本发明实施例二中co-kriging变可信度近似模型model4在设计空间的映射图；
[0058]
图10h是本发明实施例二中co-kriging变可信度近似模型model4在设计空间的三维图；
[0059]
图10i是本发明实施例二中rbf单可信度近似模型在设计空间的映射图；
[0060]
图10j是本发明实施例二中rbf单可信度近似模型在设计空间的三维图；
[0061]
图11a是本发明实施例二中基于co-kriging变可信度近似模型的pso算法工作过程中种群初始状态分布图；
[0062]
图11b是本发明实施例二中基于co-kriging变可信度近似模型的pso算法工作过程中种群中间状态分布图；
[0063]
图11c是本发明实施例二中基于co-kriging变可信度近似模型的pso算法工作过
程中种群最终状态分布图；
[0064]
图11d是本发明实施例二中基于co-kriging变可信度近似模型的pso算法工作过程中迭代过程适应度值趋势图；
[0065]
图12a是本发明实施例二中基于rbf单可信度近似模型的pso算法工作过程中种群初始状态分布图；
[0066]
图12b是本发明实施例二中基于rbf单可信度近似模型的pso算法工作过程中种群中间状态分布图；
[0067]
图12c是本发明实施例二中基于rbf单可信度近似模型的pso算法工作过程中种群最终状态分布图；
[0068]
图12d是本发明实施例二中基于rbf单可信度近似模型的pso算法工作过程中迭代过程适应度值趋势图。
具体实施方式
[0069]
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
[0070]
本实施例提供一种技术方案：焊点焊后残余应力预测与优化的变可信度模型确定方法，包括以下步骤：
[0071]
步骤1：将回流焊过程分为焊点熔化成型和冷却凝固两个阶段，选择对焊后残余应力影响显著的焊点尺寸参数，即焊点高度h和焊盘半径r作为输入，在焊点融化成型阶段，基于交互式有限元分析软件surface evolver模拟单个bga焊点三维形态；
[0072]
步骤2：将模拟的bga焊点形态转化为点、面参数文件，并基于三维数字化软件solidworks转为实体文件，与其他组件进行装配，构建bga焊点三维模型；
[0073]
步骤3：基于数值模拟综合分析软件ansys workbench进行焊后残余应力分析，采用间接耦合方法，以典型回流曲线为边界条件边界条件，对回流焊过程进行仿真计算，获得焊后最大残余应力；
[0074]
步骤4：对两个阶段分别使用增加细化次数的方式和使用edge sizing模块细化网格，进行有限元分析的收敛性检验，确定高/低精度样本点的划分方式；
[0075]
步骤5：使用拉丁超立方采样方法(latin hypercube sampling,lhs)，获取嵌套与非嵌套的高/低精度样本点；
[0076]
步骤6：根据高/低精度样本点和划分方式和样本点嵌套/非嵌套关系，构建四种co-kriging变可信度近似模型，构建rbf单可信度近似模型；
[0077]
步骤7：选取最大相对误差mre和均方相对误差mrse作为近似模型精度的评价标准，对变可信度近似模型和单可信度近似模型进行焊后残余应力预测效果的验证与对比；
[0078]
步骤8：使用粒子群优化算法(pso)，优化焊点尺寸参数方案，以获取bga焊点焊后残余应力的最小值；
[0079]
步骤9：分别基于co-kriging变可信度近似模型和rbf单可信度近似模型进行全局预测，进行寻优，分别获得最优焊点尺寸方案，并对结果仿真验证，结果表明co-kriging模型的预测性能更优。
[0080]
进一步地，步骤1中，基于最小能量原理和有限元数值分析方法，定义初始几何条件、能量形式和约束条件，材料性质参数等，采用surface evolver软件模拟bga焊点成型三维形态；
[0081]
优选地，步骤2中，建立bga焊点有限元仿真模型的步骤，包括：
[0082]
利用surface evolver软件使用有限元方法模拟并输出焊点模型的点、面参数文件，将焊点表面分为大量有限的小单元和节点，输出单个焊点的节点坐标参数文本文件，并将其转化为ascii码格式的面片文件，采用solidworks软件实体化焊点模型并装配成为bga组件，转为“x_t”文件。依据所述文件，采用局部细化edge sizing方法划分网格，基于ansys软件建立bga焊点有限元仿真模型。
[0083]
进一步地，步骤3中，采用ansys workbench有限元模拟综合分析软件，使用瞬态热学模块对bga组件施加温度场，将热分析的温度场结果作为载荷，并载入结构分析中进行热结构耦合分析，热分析过程的边界条件为典型回流曲线。结构分析过程的边界条件和仿真结束时刻模拟真实环境设置，对回流焊过程进行仿真计算，获得焊后最大残余应力。
[0084]
进一步地，步骤4中，在两个阶段分别使用增加份数的方式细化网格，随着网格密度的增加，最大残余应力随着网格的细化逐渐趋近于一确定值。对焊点成型模拟过程，利用surface evolver增加细化次数以提高单元密度，使得单元数目呈指数增长，并以细化次数为收敛性检验的输入参数。最大残余应力为输出参数，得到最终结果；对焊点残余应力分析，为提高计算效率并获得更细的网格划分，仅对焊点阵列四角的焊点进行细化，在ansys workbench的mesh模块内选择edge sizing的细化方法并按份数划分网格，以份数作为收敛性检验的输入参数，最大残余应力为输出参数，得到最终结果。
[0085]
进一步地，步骤5中，根据变可信度近似模型对高/低精度样本点关系的限定，可分为高/低精度样本点非嵌套和嵌套两种。对于非嵌套样本，直接用lhs方法分别生成高/低精度样本点；对于嵌套样本，先用传统试验设计方法独立生成高/低精度样本点，然后对高精度样本点平移，以满足嵌套关系。
[0086]
进一步地，步骤6中，分别基于所获得的低精度样本数据和高/低精度数据差值构建相应的kriging近似模型，并通过对样本值的最大似然估计求解co-kriging模型的五个主要模型参数值。分别为重要度参数θc、θd，指数pc、pd，缩放参数ρ，从而确定co-kriging模型的预测估计表达式；构建rbf单可信度近似模型。
[0087]
优选地，步骤6中，构造的co-kriging近似模型和相应协方差矩阵可表示为：
[0088]
ze(x)＝ρzc(x)+zd(x)
[0089][0090]
其中，x为样本点，ze表示co-kriging近似模型，zc表示低精度样本数据构建的kriging近似模型，zd表示高/低精度数据差值构建的kriging近似模型。对于独立的高/低精度两组数据，高精度数据有ne个样本xe，低精度数据有nc个样本xc，用rc(xc，xe)表示样本数据xc与xe间的相关性矩阵。
[0091]
进一步地，步骤7中，将得到的四种co-kriging模型结果与rbf单可信度近似模型结果互相对比，选取mre和mrse作为近似模型精度的评价标准，结果表明变可信度近似模型
的预测性能更好。
[0092]
优选地，步骤7中，验证对比的结果具体包括：几种变可信度近似模型的形状相似，说明了变可信度近似模型的稳定性，焊点尺寸方案易产生较大残余应力，近似模型的变化梯度较大，可见尺寸参数的改变会对残余应力产生显著影响。精度验证结果表明，变可信度近似模型的最大相对误差和均方相对误差更小，预测性能更好，非嵌套高/低精度样本点构建的近似模型误差较低，rbf模型预测效果最差。
[0093]
进一步地，步骤8中，pso算法参数设定后，用近似模型进行全局预测，pso算法产生随机初始粒子和每次迭代产生的新粒子，如此循环直到满足迭代次数要求，最终可以获得最优解，即最小残余应力。
[0094]
进一步地，步骤9中，对焊点设计空间进行网格划分，设置具体参数，分别利用co-kriging模型和rbf模型进行全局预测，并利用粒子群优化算法进行寻优，分别得到优化后的最优焊点尺寸方案。对两组优化结果进行仿真验证，得到co-kriging模型的结果与仿真值误差更小，即验证co-kriging模型预测性能更优。
[0095]
优选地，步骤9中，pso算法具体参数为：初始种群数n＝500，最大迭代次数g＝100，粒子速度v∈[-1.5,1.5]，惯性权重ω＝0.6，自我学习因子c1＝0.6，群体学习因子c2＝0.4。
[0096]
实施例二
[0097]
下面结合附图以及具体实施对本发明做进一步的说明，但本发明的保护范围并不限于此。
[0098]
焊点焊后残余应力预测与优化的变可信度模型确定方法，以bga焊点为例，实施流程图如图1所示，具体实施步骤包括：
[0099]
步骤一：基于最小能量原理和有限元数值分析方法，本实施例选择焊点材料为sn-3.8-0.7cu(sac387)，并针对单个bga焊点进行形态模拟，设置单个焊点焊料量为0.0041mm3，材料密度为7500kg/m3，表面张力取520mn/m。将回流焊过程分为焊点熔化成型和冷却凝固两个阶段，选择对焊后残余应力影响显著的焊点尺寸参数，即焊点高度和焊盘半径作为输入，通过surface evolver软件模拟得到对应的焊点形态。
[0100]
步骤二：利用surface evolver软件有限元方法，将焊点表面分为大量有限的小单元和节点，从而实现对焊点表面形态的模拟，并输出单个焊点的节点坐标参数。将节点参数“txt”文件重写为“stl”文件的ascii码格式，生成焊点模型的面片文件，但其不能直接被用于仿真分析。另将“stl”文件导入三维建模软件solidworks中转为实体文件，并构建芯片与pcb等其他组件进行装配，并转为“x_t”文件用于后续回流焊后焊点残余应力的有限元模拟分析。
[0101]
作为可选地实施方式，步骤二具体为：依据所述文件，基于ansys软件建立bga焊点有限元仿真模型，包括芯片、bga焊点和pcb等组件。其中，芯片尺寸为5mm
×
5mm
×
0.5mm；焊点数量为36个，按6
×
6阵列排列，焊点间距为0.65mm；pcb尺寸为10mm
×
12mm
×
1mm。焊点材料选用sac387，焊点与其他组件材料属性如表1所示。且根据初次仿真得到的结果，最大残余应力位置出现在离阵列中心最远的焊点，即四角焊点，故采用局部细化edge sizing方法划分网格，同时为减小计算量，本实施例仅在离阵列中心最远处的焊点处选择细网格划分，在芯片和pcb上选择粗网格划分，参考图3a-3c，建立bga焊点有限元仿真模型。
[0102]
表1各部分材料属性
[0103][0104]
步骤三，bga焊点焊后残余应力有限元分析采用间接耦合方法，确立边界条件和仿真结束时刻，对回流焊过程进行仿真计算，获得焊后最大残余应力；
[0105]
作为可选地实施方式，步骤三具体为：使用瞬态热学模块对bga组件(bga焊点、pcb、芯片)施加温度场，由于不同组件对流换热系数不同，施加回流焊温度曲线后组件上出现温度梯度，将热分析的温度场结果(参考图5)作为载荷，并载入结构分析中进行热结构耦合分析。热分析过程的边界条件为无铅焊料的典型回流曲线，参考图4，回流焊过程包括预热、保温、再流和冷却四个主要阶段。图4中a点对应凝固温度217℃，冷却结束时及保温阶段温度对应常温25℃。结构分析过程的边界条件为模型中pcb底面z方向位移约束，bga焊点焊后残余应力结果参考图6。
[0106]
步骤四：使用增加份数的方式细化网格，以进行有限元分析的收敛性检验，确定精度；
[0107]
步骤四具体如下：
[0108]
有限元收敛性是指当仿真模型的网格逐渐加密时，有限元求解的结果逐渐稳定、越来越接近真实解的过程。本实施例使用增加份数的方式细化网格，以进行有限元分析的收敛性检验。参考图7，为回流焊仿真模拟中网格份数与最大残余应力的关系，随着网格密度的增加，最大残余应力的值逐渐趋近于一固定值，得到了一条收敛的网格尺寸与残余应力的关系曲线。本实施例选取最大误差为10％，参考图7中虚线所示，仅选取虚线范围内的值。
[0109]
对焊点成型模拟过程，采取多次细化的方式。有限元分析软件surface evolver求解bga焊点三维形态问题的计算过程中，通过增加细化次数增加单元密度，使得单元数目呈指数增长。以细化次数为收敛性检验的输入参数，最大残余应力为输出参数，最终结果为高精度模型细化次数为3次，低精度模型细化次数为2次。
[0110]
对焊点残余应力分析，在ansys workbench的mesh模块内施加局部网格细化命令。根据已有研究结果，csp焊点焊后残余应力的分布不均匀：即距离焊点阵列中心位置越远的焊点，最大残余应力越大，且出现在焊点与焊盘交界处。为提高计算效率并获得更细的网格划分，仅对焊点阵列四角的焊点进行细化。选择edge sizing的细化方法并按份数划分网格，按比例划分偏差增长率为1.5。以份数作为收敛性检验的输入参数，最大残余应力为输出参数，最终结果为高精度模型单边网格份数为6，低精度模型单边网格份数为3。细化次数为3次，低精度模型细化次数为2次。
[0111]
步骤五：采取基于lhs采样方法的试验设计；
[0112]
作为可选地实施方式，步骤五具体为：根据变可信度近似模型对高/低精度样本点关系的限定，可分为高/低精度样本点非嵌套和嵌套两种，即高精度样本点是否为低精度样本点的子集。对于非嵌套样本，直接用lhs方法分别生成高/低精度样本点；对于嵌套样本，
先用传统试验设计方法独立生成高/低精度样本点，然后对高精度样本点平移，以满足嵌套关系。参考图8是两个设计变量的拉丁超立方抽样的结果，可以看出样本点在设计空间内分布均匀，其中图8-a为高/低精度样本点非嵌套的采样结果，图8-b为高/低精度样本点嵌套的采样结果。焊点高度的变化范围是[0.27,0.33]，焊盘半径的变化范围是[0.17,0.20]，低精度样本点数取为40，高精度样本点数取为20。
[0113]
步骤六：构建基于co-kriging变可信度近似模型，并验证近似模型的精度；
[0114]
步骤六具体如下：
[0115]
kriging模型可以表示为：y(x)＝β+z(x)，式中：β为常数；z(x)为均值为零但方差不为零的高斯随机过程。
[0116]
点xi与xj之间z(x)的协方差可表示为：cov(z(xi)，z(xj))＝σ2r(xi，xj)，
[0117]
式中：σ为随机过程z(x)的方差；r(xi，xj)为相关函数。
[0118]
常用的高斯相关函数可以表示为：式中：θk为维度k上未知的超参数；d为设计变量的维度。
[0119]
式中：θk为维度k上未知的超参数；d为设计变量的维度。
[0120]
通过最大化kriging模型在观测点的极大似然函数，求得kriging模型的超参数β，σ以及θk，可得kriging模型在未观测点处的预测响应及预测均方差：
[0121][0122][0123]
式中：r为未观测点x与样本点之间的相关性向量；r为样本点间的相关性矩阵；y为样本点处的响应的向量；l为元素均为1的列向量。
[0124]
co-kriging方法是kriging方法的扩展，构造的近似模型可表示为ze(x)＝ρzc(x)+zd(x)，式中：zc表示低精度样本数据构建的kriging近似模型，zd表示高/低精度数据差值构建的kriging近似模型。对于独立的高/低精度两组数据，高精度数据有ne个样本xe，低精度数据有nc个样本xc，用rc(xc，xe)表示样本数据xc与xe间的相关性矩阵，co-kriging模型的协方差矩阵可表示为：
[0125][0126]
此模型有两种相关函数，需要估计更多的相关参数如θc、θd、pc、pd和放大因子ρ等。co-kriging模型的相关参数θc、pc的估计求解与kriging方法相同。为求解超参数θd、pd和ρ，定义：d＝y
e-ρyc(xe)，式中：yc(xe)为yc在与xe相同位置的值。类似kriging模型的超参数求解方法，可以求得对θd、pd和ρ的估计。co-kriging模型的预测估计表达式为：
[0127][0128][0129]
[0130]
co-kriging模型的均方误差估计表达式为：
[0131][0132]
均方相对误差(mean relative square error，mrse)可以衡量预测值与真实值之间的整体偏差，本实施例的建模目标是获得一个全局精确的模型，选用均方相对误差可以实现这一目标。最大相对误差(max relative error，mre)可以显示近似模型预测的最大偏差，以确定模型预测性能最差的位置。因此，选择mrse和mre两个指标验证近似模型的精度：
[0133][0134]
式中，yi为各点真实值，为近似模型预测值。
[0135]
径向函数是以待测点与样本点之间的距离为自变量的一类函数，以径向函数为基函数，线性组合得到的函数称为径向基函数模型(rbf)。径向函数可表示为：
[0136][0137]
式中，ωi权函数，ri＝||x-xi||，待测点x与样本点xi间的欧氏距离。φ(r)是径向函数，常用径向函数包括高斯函数、多元二次函数、逆多元二次函数、薄板样条函数等。考虑薄板样条函数φ(r)＝r2ln(r)可以避免多项式差值的摆动现象，本发明选择薄板样条函数作为基函数进行rbf单可信度近似模型的构建。
[0138]
步骤七：选取mre和mrse作为近似模型精度的评价标准，对所构建变可信度近似模型进行焊后残余应力预测效果的验证与对比；
[0139]
作为可选地实施方式，步骤七具体为：由收敛性检验得到两组高/低精度样本点，高/低精度样本点满足嵌套或非嵌套的关系，组合可得到四组高/低精度样本点，分别建立co-kriging变可信度近似模型model1、model2、model3和model4，参考图9。径向基函数模型rbf和四种co-kriging模型结果参考图10。取最大相对误差mre和均方相对误差mrse作为近似模型精度的评价标准，co-kriging模型与rbf模型的精度验证结果如表2所示。
[0140]
表2 co-kriging模型与rbf模型的精度结果对比
[0141]
模型mremrsemodel10.16780.082model20.10990.0594model30.16170.0726model40.11980.0600rbf0.27930.1138
[0142]
参考图10a-10j，由近似模型构建结果可以看出，几种变可信度近似模型的形状相似，说明了变可信度近似模型的稳定性，预测的残余应力最大值均在(0.3，0.19)附近，此区域的焊点尺寸方案易产生较大残余应力，近似模型的变化梯度较大，可见尺寸参数的改变会对残余应力产生显著影响。另外，除峰值附近区域，其他部分目标函数值相差不大，可知设计时需要避免此区域并能够大致确定最优尺寸参数的位置。
[0143]
如表3所示，由近似模型精度验证结果可知，变可信度近似模型的最大相对误差和
均方相对误差更小，预测性能更好，能够实现对bga焊点回流焊后最大残余应力的准确预测。mrse可以体现模型的整体预测精度，四种变可信度近似模型的mrse均在10％以内，而单可信度近似模型rbf的mrse为11.38％，可见变可信度近似模型的预测精度更高。最大相对误差mre体现模型的最大偏离程度，非嵌套高/低精度样本点构建的近似模型model2和model4误差较低，model1和model3预测效果稍差，rbf模型预测效果最差。
[0144]
步骤八：将粒子群优化算法与co-kriging变可信度近似模型结合，迭代寻找最优解，即bga焊点焊后残余应力的最小值；
[0145]
作为可选地实施方式，步骤八具体为：pso算法参数设定后，将pso算法与co-kriging模型结合，pso算法产生随机初始粒子和每次迭代产生的新粒子，通过近似模型预测相应的残余应力值，如此循环直到满足迭代次数要求，最终可以获得最优解，即最小残余应力。
[0146]
步骤九：分别利用co-kriging模型和单可信度模型进行全局预测，并利用粒子群优化算法进行寻优，分别获得最优焊点尺寸方案，并对结果仿真验证，结果表明co-kriging模型的预测性能更优。
[0147]
作为可选地实施方式，步骤九具体为：初始设置中对焊点高度和焊盘半径的设计空间(即[0.27,0.33]，[0.17,0.20])分别进行网格划分(粒度为0.001和0.0005)。具体参数设置为：初始种群数n＝500，最大迭代次数g＝100，粒子速度v∈[-1.5,1.5]，惯性权重ω＝0.6，自我学习因子c1＝0.6，群体学习因子c2＝0.4。分别利用co-kriging模型和rbf模型进行全局预测，并利用粒子群优化算法进行寻优，优化过程中种群最初分布状态、中间状态、最终状态及收敛过程分别参考图11a-11d、12a-12d。
[0148]
参考图12a-12b可知，不论是基于变可信度近似模型还是单可信度近似模型的全局预测，pso算法的收敛过程差异较小，适应度值均在20次左右迭代之后收敛，优化效果较为接近。两个近似模型的寻优结果有一定的差异：基于co-kriging模型的优化后最优焊点尺寸方案为h＝0.3172mm，r＝0.1732mm，此时的bga组件焊后最大残余应力为49.36mpa；基于rbf模型的优化后最优焊点尺寸方案为h＝0.2933mm，r＝0.1832mm，此时的bga组件焊后最大残余应力为48.87mpa。对两组优化结果进行仿真验证，得到其与仿真值的误差分别为1.59％与5.30％，可知co-kriging模型的预测性能更优，表3为两种方案的对比。
[0149]
表3焊点尺寸参数最优方案对比
[0150][0151][0152]
综上所述，上述实施例的焊点焊后残余应力预测与优化的变可信度模型确定方法，通过对回流焊工艺过程进行仿真，获得焊后的最大残余应力，据此建立焊点高度和焊盘半径等尺寸参数与焊后最大残余应力间的变可信度近似模型，采取嵌套和非嵌套的采样方式组成四种变可信度近似模型，并进行预测效果对比，同时使用粒子群优化算法进行焊点尺寸参数的优化，结果验证了变可信度近似模型的预测效果更优并获得其最佳构建方案。
[0153]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

技术特征：
1.焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于，包括以下步骤：s1：将回流焊过程分为焊点熔化成型和冷却凝固两个阶段，选择对焊后残余应力影响显著的焊点尺寸参数，即焊点高度h和焊盘半径r作为输入，在焊点融化成型阶段，基于交互式有限元分析软件surface evolver模拟单个bga焊点三维形态；s2：将模拟的bga焊点形态转化为相关文件，基于三维数字化软件solidworks构建bga焊点有限元仿真模型；s3：基于数值模拟综合分析软件ansys workbench进行焊后残余应力分析，采用间接耦合方法，以典型回流曲线为边界条件边界条件，对回流焊过程进行仿真计算，获得焊后最大残余应力；s4：对两个阶段分别使用增加细化次数的方式，并使用edge sizing模块细化网格，进行有限元分析的收敛性检验，确定高/低精度样本点的划分方式；s5：使用拉丁超立方采样方法，获取嵌套与非嵌套的高/低精度样本点；s6：根据高/低精度样本点、划分方式以及样本点嵌套/非嵌套关系，构建四种co-kriging变可信度近似模型，同时构建rbf单可信度近似模型；s7：将co-kriging变可信度近似模型和rbf单可信度近似模型进行焊后残余应力预测效果的验证与对比；s8：使用粒子群优化算法，优化焊点尺寸参数方案，以获取bga焊点焊后残余应力的最小值；s9：分别基于co-kriging变可信度近似模型和rbf单可信度近似模型进行全局预测，进行寻优，分别获得最优焊点尺寸方案，并对结果仿真验证，结果表明co-kriging变可信度近似模型的预测性能更优。2.根据权利要求1所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s1中，基于最小能量原理和有限元数值分析方法，定义初始几何条件、能量形式、约束条件、材料性质参数，采用surface evolver软件模拟bga焊点成型三维形态。3.根据权利要求2所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s2中，构建bga焊点有限元仿真模型的具体过程如下：s21：利用surface evolver软件使用有限元方法模拟并输出焊点模型的点、面参数文件，将焊点表面分为大量有限的小单元和节点，输出单个焊点的节点坐标参数文本文件，并将其转化为ascii码格式的面片文件，采用solidworks软件实体化焊点模型并装配成为bga组件，转为“x_t”格式文件；s22：根据“x_t”格式文件，采用局部细化edge sizing方法划分网格，基于ansys软件建立bga焊点有限元仿真模型。4.根据权利要求3所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s3中，采用ansys workbench有限元综合分析软件，使用间接耦合方法，使用瞬态热学模块对bga组件施加温度场，将热分析的温度场结果作为载荷，并载入结构分析中进行热结构耦合分析，热分析过程的边界条件为典型回流曲线，结构分析过程的边界条件和仿真结束时刻模拟真实环境设置。
5.根据权利要求4所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s4中，在两个阶段分别使用增加份数的方式细化网格，随着网格密度的增加，最大残余应力随着网格的细化逐渐趋近于一确定值；其中，对焊点成型模拟过程，利用surface evolver增加细化次数以提高单元密度，使得单元数目呈指数增长，并以细化次数作为收敛性检验的输入参数，最大残余应力为输出参数，得到最终结果；对焊点残余应力分析，仅对焊点阵列四角的焊点进行细化，在ansys workbench的mesh模块内选择edge sizing的细化方法并按份数划分网格，以份数作为收敛性检验的输入参数，最大残余应力为输出参数，得到最终结果。6.根据权利要求5所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s5中，根据co-kriging变可信度近似模型对高/低精度样本点关系的限定，分为高/低精度样本点非嵌套、嵌套两种，其中对于非嵌套样本，直接用lhs方法分别生成高/低精度样本点；对于嵌套样本，先用传统试验设计方法独立生成高/低精度样本点，然后对高精度样本点平移，以满足嵌套关系。7.根据权利要求6所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s6中，分别基于所获得的低精度样本数据和高/低精度数据差值构建相应的kriging近似模型，并通过对样本值的最大似然估计求解co-kriging变可信度近似模型的模型参数值，分别为低精度样本点构建的kriging近似模型的重要度参数θ
c
、高/低精度数据差值构建相应的kriging近似模型的重要度参数θ
d
、低精度样本点构建的kriging近似模型的基函数指数p
c
、高/低精度数据差值构建相应的kriging近似模型的基函数指数p
d
、缩放参数ρ，从而确定co-kriging变可信度近似模型的预测估计表达式，其中，构建的co-kriging变可信度近似模型表示如下：z
e
(x)＝ρz
c
(x)+z
d
(x)其中，x为样本点，z
c
表示低精度样本数据构建的kriging近似模型，z
d
表示高/低精度数据差值构建的kriging近似模型；对于独立的高/低精度两组数据，高精度数据有n
e
个样本x
e
，低精度数据有n
c
个样本x
c
，用r
c
(x
c
，x
e
)表示样本数据x
c
与x
e
间的相关性矩阵，σ
c
和σ
d
分别表示z
c
和z
d
的方差，l表示值全为1的列向量；co-kriging变可信度近似模型的协方差矩阵表示为：co-kriging变可信度近似模型的预测估计表达式为：8.根据权利要求7所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s7中，选取最大相对误差mre和均方相对误差mrse作为近似模型精度的评价标准，对四种co-kriging变可信度近似模型和rbf单可信度近似模型进行焊后残余
应力预测效果的验证与对比；其中，mrse和mre两个指标表达式分别如下：其中，y
i
为各点真实值，为近似模型预测值。9.根据权利要求8所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s8中，粒子群优化算法参数设定后，用co-kriging变可信度近似模型进行全局预测，粒子群优化算法产生随机初始粒子和每次迭代产生的新粒子，循环直到满足迭代次数要求，最终可以获得最优解，即最小残余应力。10.根据权利要求9所述的焊点焊后残余应力预测与优化的变可信度模型确定方法，其特征在于：在所述步骤s9中，对焊点设计空间进行网格划分，设置具体参数，分别利用co-kriging变可信度近似模型和rbf单可信度近似模型进行全局预测，并利用粒子群优化算法进行寻优，分别得到优化后的最优焊点尺寸方案；对两组优化结果进行仿真验证，得到co-kriging变可信度近似模型的结果与仿真值误差更小，即验证co-kriging变可信度近似模型预测性能更优。

技术总结
本发明公开了焊点焊后残余应力预测与优化的变可信度模型确定方法，属于电子元器件封装技术领域。本发明通过对回流焊工艺过程进行仿真，获得焊后的最大残余应力，据此建立焊点高度和焊盘半径等尺寸参数与焊后最大残余应力间的变可信度近似模型，采取嵌套和非嵌套的采样方式组成四种变可信度近似模型，并进行预测效果对比，同时使用粒子群优化算法进行焊点尺寸参数的优化，结果验证了变可信度近似模型的预测效果更优并获得其最佳构建方案。的预测效果更优并获得其最佳构建方案。的预测效果更优并获得其最佳构建方案。


技术研发人员：王梅 王儒 任好 张楚雷 任海林 黄珊 张云斐 王荧飞 王国新 胡子翔
受保护的技术使用者：北京理工大学
技术研发日：2023.05.10
技术公布日：2023/8/14