您好,欢迎来到筏尚旅游网。
搜索
您的当前位置:首页基于子模型的电子封装焊点寿命预测方法研究

基于子模型的电子封装焊点寿命预测方法研究

来源:筏尚旅游网
E 如 H哆 环境适应性和可靠性_■ 基于子模型的电子封装焊点寿命预测方法研究 刘功桂 肖慧。李晓延。 (1.中国电器科学研究院有限公司,广州 510300 2.北京工业大学,北京 100022) 摘要:在向无铅化过渡的过程中,封装材料与工艺改变所带来的最突出的问题之一就是无铅焊点的可靠性问题。由于焊点 区非协调变形导致的热疲劳失效是电子封装焊点的主要失效形式。到目前为止,仍无公认的焊点寿命和可靠性的评价方法。 本文采用统一型粘塑性本构模型描述SnAgCu焊点的变形行为,在对焊点可靠性分析的基础上采用基于能量的寿命模型进行 焊点寿命预测;在有限元方面采用子模型的处理技术,并采用节点平均化以及体积平均化方法计算寿命模型相关参数。力 图从理论方面和有限元处理技术方面改进以达到提高焊点寿命预测效率及预测精度的目的。 关键词:电子封装焊点:热疲劳;寿命预测:子模型 中图分类号:TG407 文献标识码:A 文章编号:1004—7204(2012)03—0026—05 Life Prediction Research for Solder Joints in Electronic Package Based on the Submodel Method LIU Gong—gui ,XIAO Hui ,LI Xiao-yan (1.China Electric Apparatus Research Institute Co.,Ltd,Guangzhou 510300 2.Beijing University of Technology,Beijing 100022) Abstract:In the transition to 1ead—free process,one of the most prominent probleros as the changes of packaging materialS and procesS has brought about i S the 1ead—free solder j0int reliabi1ity. Thermomechanical fatigue failure caused by nonconforming deformation of the soldering area iS the main failure mode for solder joints in electronic package.Up to now,there iS no 1 ife prediction and reliability evaluation method for solder joint which iS generally accepted. In thi S paper,a unified viscoplastiC constitutive model was used to describe the deformation behavior of SnAgCu solder joint, and the solder joint 1ife was predicted based on solder joint reliabi1ity analysiS.For finite element analysiS,submodel method was adopted,and the 1ife mode1 parameters were obtaine—d by data average processing,such as node average and element average.Based on the theory optimization and effective data processing for finite element analysiS,the efficiency as wel1 as the preciSion of]ife prediction for solder joint may be improved. Key words:solder joints in electronic package:submodel method;1ife prediction;submodel 前言 无铅化组装已成为电子产品的必然选择,封装材料 与工艺改变所带来的最突出的问题之一就是无铅焊点的 可靠性问题Ⅲ。疲劳破坏是电子封装焊点破坏的最主要 的形式,总的来说,对焊点破坏行为及寿命的研究已有 有关SnPb焊点疲劳模型和寿命预测方法的研究也有不 少文献可以参考。目前,虽然国内外有关无铅焊料开 发、焊接工艺和焊点可靠性的研究已有不少报道,但能 否将用于描述SnPb焊点破坏和寿命的规律直接应用到 SnAgCu系无铅焊料焊点上,仍无理论和实践的保证。无 疑,开展无铅焊点破坏行为和寿命规律的研究对于理解 多年的基础,许多寿命预测方法被不同的研究者采用, 基金项目:国家自然科学基金(50871004),北京市自然科学基金(2112005) Environmental Technology・June 201 2 26 nvironmen talAdaptabi1i哆 I环境适应性和可靠性 无铅焊料焊点的破坏机理、预测其可靠服役寿命有重要 的理论和实际意义。 目前,对焊点破坏和寿命的研究已引起国内外学者 的广泛关注,提出了SnPb焊料和SnAgCu焊料寿命的不 同模型 】,这些模型或者以应变幅为表征参量,或者以 累积蠕变应变为表征参量,或者以累积蠕变应变能密度 为表征参量。在这些模型中,美国Amkor公司Syed ̄ 提 出的基于累积蠕变应变和基于累积蠕变应变能密度的寿 命模型得到较多关注和应用。根据Syed模型,蠕变是控 制SnAgCu焊料损伤的主要机制,蠕变是描述焊料损伤 行为的单一参量,因而寿命预测也应以蠕变变形为基础。 此外还有Darveaux和Baneqi川,Pang John H.L.等嘲提出 蠕变和塑性分离的唯象学模型,将与时间无关的塑性应 变和与时间相关的蠕变应变分开处理。但是,基于非线 性连续介质力学方法,非弹性应变(塑性应变、蠕变变形) 速率与位错运动速率相关,产生于同一机制。此外,在 表征焊点应力应变关系时,忽略与时间无关的塑性变形, 或者将塑性应变与蠕变相分离的本构理论虽然可以在某 些程度上使模型得到简化,但在处理材料循环塑性变形、 蠕变以及它们的相互作用时表现较差。 在电子组装设计中如何运用疲劳模型进行寿命预测 也是企业十分关心的问题。由于电子产品市场变化加 快,不断要求设计出新的组装产品,企业几乎没有时间 进行疲劳寿命试验。对新产品设计工程师来说,就需要 一种快速的方法确定疲劳寿命。应该指出,随着计算机 技术及计算模型的发展,微组装焊点的力学响应特征的 FEM数值模拟工作已经取得了很大发展,如ANSYS、 ABAQUS、MARC、ADINA等大型商用有限元分析软件, 它们为电子产品的可靠性结构分析提供了强有力的技术 支持。但是,当用有限元数值模拟方法分析无铅焊点可 靠性及焊点寿命时还存在许多不足,如无铅材料热物理 性能及无铅焊点可靠性实验数据库的不充分,无铅钎料 的热力本构关系研究不深入等,必将影响寿命预测的准 确性。 本文运用Anand本构模型描述SnAgCu无铅焊点的 27 70】?年06月・环境技术 粘塑性材料特性,对通孔插装焊点以及BGA焊点封装结 构在热循环载荷下的应力应变行为进行了ANSYS有限元 分析,采用子模型处理技术并且应用基于能量方法对无 铅焊点的疲劳寿命进行了预测。力图提供一种能够充分 利用现有无铅组装材料的性能数据,更加精确、有效预 测电子封装无铅焊点的疲劳寿命的计算方法。 1分析模型 1.1结构模型及有限元模型 本文针对如今普遍应用的通孔插装焊点以及BGA封 装焊点,采用数值模拟的方法对其在热循环服役条件下 的应力应变响应进行有限元模拟,并进行寿命评估。考 虑到试件结构模型的复杂性,为了提高计算效率,在有 限元计算时本文采用子模型处理技术,对关键位置焊点 进行网格细化以得到较精确的模拟结果 l。 子模型是得到模型部分区域更加精确解的有限单元 技术。子模型方法又称为切割边界位移法或特定边界位 移法,切割边界就是子模型从整个较粗糙的模型分割开 的边界。整体模型切割边界的计算位移值即为子模型的 边界条件。计算时。全局模型和子模型的热循环加载条 件一致。 本文在研究通孔插装焊点热疲劳寿命预测时,首先 根据实际插装PCB的Protel线路图建立试件全局分析模 型,计算得到关键插孔位置;进而在关键插孔位置建立 焊点子模型,研究焊点应力应变场规律。 如图1所示,PCB尺寸为157.5ram×149.5mm×1.63 mm,图中红色小点即代表不同尺寸的插孔,PCB下方的 两个黑圈代表安装孔。由于在三级电子封装中,PCB通 常是通过互连插座与母板连接起来的,如图2所示 , 因而确定如图3所示试件全局有限元模型,并在靠近 PCB安装孔的下底面施加全约束。 在PCB全局有限元模型中,为节省计算时间,网格 划分采用了较稀疏网格,但不影响整个结构的变形分布 规律,对关键插孔位置的确定没有影响。但在进行寿命 预测时,焊点部分的网格密度,对计算结果影响较大, 一 nvironmen ta1A帅tabili哆I环境适应性和可靠性 1.3热循环加载条件 温度加载采用循环温度载荷。根据IPC9701标准 , 选用热循环温度范围为一40~125 ,升/降温时间均为 5min,高/低保温时间均为lOmin,零应力应变参考温度 设为25℃。由于焊点的应力应变在热循环过程中呈周期 性变化,一般认为在第四周期趋向稳定,所以加载四个 周期的温度循环,进而分析焊点的应力应变场变化规律, 热循环曲线如图7所示。 1.4寿命模型 采用基于塑性应变能密度的Mo ̄ow方程计算焊点寿 命,如(4】式: Ⅳ, Vw =A (4) 对于Sn3.8AgO.7Cu,n--O.897,A=311.7MPa㈣。提取 热循环第4周的焊点平均塑性应变能增量作为方程(4) 中Vw 值。 2计算结果分析 2.1通子L插装焊点结果分析 焊点子模型分析时,热循环加载曲线与PCB全局计 算时的完全相同。图8所示为钎料中的粘塑性应变分布 云图。从图中可以看出钎料与引线、以及钎料与PCB的 p 媸 时间t/s 图7热循环加载曲线 图8通孔插装焊点焊料部分蠕变应变分布图及 寿命计算时所取节点 29 201 2年06月・环境技术 交界面处于高应变区,尤其是引线与钎料相界的拐角处, 这与刘等 的实验研究结果刚好吻合。由于热膨胀系数 ∞ 加 ∞ ∞ ∞ 曲 ∞ 0 加 们 不匹配的原因,使得钎料与引线、钎料与PCB的交界面 处于高应力应变区,导致裂纹沿PCB与引线,以及焊料 与引线接触面扩展。 选取焊点界面的节点进行时间历程后处理。得到如 图9所示为焊点塑性功密度变化图。由图9可知,焊点 中的塑性功密度随热循环不断累积。 焊点寿命计算时,提取钎料与引线、钎料与PCB的 交界面节点的塑性功密度,如图8所示。采用节点平均 化的方法计算热循环每周平均塑.『生功密度增量,如(5) 式所示。图10所示为不同重量元器件对应焊点的寿命比 较图,由图可知,随着焊点所承受器件载荷的增大,焊 点寿命减低。 N △ △ 上 (5) 2.2 BGA封装焊点结果分析 对整体模型分析,得到关键焊点的位置位于器件芯 片边缘的正下方。这是由于芯片材料与焊点材料的热膨 胀系数差异较大造成的。子模型分析结果表明,焊点中 最大塑性变形点位于与IMC层相邻的钎料层(这里简称 d 皇 黼 蚤 掣 时间t/s 图9器件重量为0.02N时对应焊点塑性功密度变化 器件重量/N 图1O不同重量元器件对应焊点的热疲劳寿命 nvironmen.tal Ada 1i哆I环境适应性和可靠性l■ 3结论 本文针对如今普遍应用的通孔插装焊点以及BGA封 为焊点连接层)中,如图11所示。 选取焊点连接层中的节点进行时间历程后处理,图 12所示为IMC厚度为4.5Ixm,经历如图7所示热循环时 焊点的塑性应变功密度曲线。由图可知,焊点中的塑性 功密度随热循环增加而不断累积。 寿命计算时,提取与IMC接触的焊点连接层最外侧 装焊点,采用数值模拟的方法对其在热循环服役条件下 的应力应变响应进行有限元模拟,并进行寿命评估,得 到如下结论: (1)在热循环过程中,焊点的最大塑性应变能密度 单元的塑性功密度,如图11所示。采用单元平均化的方 法计算每周平均塑性功密度增量,如(6)式所示。图 13所示为不同厚度IMC焊点的寿命比较图,由图可知, 随着焊点中界面IMC厚度的增加,焊点的寿命下降。 Ⅳ ∑△ ・ △ = 一 f61 、 图11 焊球及焊点连接层最外侧单元环塑性功密度分布 曼 鲁 蜘 督 时间t/s 图12焊点塑性功密度变化 匪 \ 据 艟 图l3不同厚度IMC焊点的热疲劳寿命 点都是出现在焊点连接界面层中,对于通孔插装焊点和 BGA封装焊点均如此,因而预测焊点界面连接层是焊点 首先出现裂纹的地方。 (2)从模拟得到的焊点塑性应变能密度时间历程图 中可以看出,Anand本构模型能够合理地模拟SnAgCu焊 点的时间温度相关的粘塑性材料特性。 (3)采用子模型的分析方法并结合节点平均化以及 单元平均化的数据处理方法能有效提高有限元分析效率 及焊点寿命预测精度。 参考文献 【.】Andersson C,Tegehall P E,Andersson D R.Thermal cycling aging effect on the shear strength,micr3structure intermetallic compounds(IMC)and crack initiation and propagation of reflow soldered Sn一3.8Ag一0 7Cu and wave solcered Sn一3 5Ag ceramic chip components[J】.IEEE Transactions on Components and Packaging Technology,2008,31【2:331—2】344 【2 LeeWW,Ngnyen L2】 T,SelvadurayGS.Solderjointfatigue models:review cln d applicability to chip scale packages[J】. Micr。electrO n1cs Reliability 2000 44:231—244. 【3】Syed A.Accumulated creep strain aqd energy density based thermal fatigue life prediction models fc・r SnAgCu solder joints[A】. 54th ECTC 2004[C],Las Vegas,Nevada,USA,June 1—4,2004.737—746 【4】Darveaux R Banefii K.Constitutive Relations for Tin—Based— Solder Joints[J】IEEE Transactions on Components and Packaging Technologies,1992,15(61:1013-1024. [5】Pang J H_ Xiong B S,Che F.X Modeling Stress Strain Curves For Lead—Free 95.5Sn一3 8Ag一0 7Cu Solder[A1 Thermomechanical Phenomena in Electronic Systems—Proceedings of the Intersociety Conference[C] 2004 2:131・136 [6】陈玉振,陈亚蜂,瞿亦峰.子模型技术在强度分析中的应用[j】.机械 2008,35(8):28—30 【7】Lau J H.Ball grid a ay technologies New York:McGraw—Hill, 1 995. 【8]况延香,朱颂春三级微电子封装技术[J].电子工艺技术,2004,25(3): I34—137. (下转第57页) Environmental Technology・June 201 2 30 E nTvesirt oEnqmeuipnmetal nt 环境试验设备 表1工件表面温度 时间 12:40 13:25 13:50 1#轮温度 2#轮温度 3#轮温度 4#轮温度 42.5 47 37 42.5 47 38 41 45 36 42 47 37 (上接第30页) 【9]RODGERS B,FLO0D B,PUNCH J,et aI Expe rimenfal determina rion and finite element modeI vaIida rion of the 14:20 47 47.5 46 47.5 anand viscoplasticity model constants for SnAgCuIA]Thermal Mechanical and Multi—Physics Simulation and Experiments in 速进行调节,控制精度高、控制指标均匀性好。 Micro—Electronics and Micro—Systems【C】.Piscataway:IEEE, 2005 490—496 (2)溶液槽内溶液能够自动补给、自动调节液位高度, 试验过程中的试验条件数据可自动记录。 (3)试验箱具有高度的用电安全性,具有故障紧急停 机功能。 f10]CHE F X,PANG J H Ll XU L H.Investigation of IMC layer effect on PBGA solder joint thermal fatigue reliabilify[A].Proceedings of 7th Electronics Packaging Technology Conference[C】 Piscataway:IEEE,2005.427—430 【11】CHE F X PANG J H L XU L H.IMC consideratjon in FEA simulation for Pb—free solder joint reliability[A1 Thermomechanical Phenomena in Electronic Systems—Proceedings of the Intersociety Conference[C】.Piscataway:IEEE,2006.1018—1023. 参考文献 【1]LYON S B,THOMPSON G E,JOHNSON J B New Methods for 【1 2]IPC 970I:Performance Test Methods and Qualification Requirements for Surface Mount Solder Attachments,2002 【13】Che F X,Pang J H L.Thermal Fatigue Reliability Analysis for PBGA with Sn一3.8Ag・0.7Cu Solder Joints【A1 2004 Electronics Corrosion Testing of Aluminum AIIoys[M1 Philadelphia:ASTM STP 1134,1 992:20 Packaging Technology Conference【C】.Piscataway,NJ USA: [2】BALDWIN K G,SMITH C F E Accelerated Corrosion Test for IEEE,2004:787—792 Aerospace Materials:Current Limitation and Future Trends[J】 Aircraft Engineering and Aerospace Technology,1 999,71(3): 【14】刘娜热循环条件下焊点失效行为的研究【D】.北京,北京工业大学, 2009. 239—244 [3】徐明德,刘颖周期浸润腐蚀试验方法【K】.金属腐蚀与防护标准手册, 北京:中国环境科学出版社,I996:52 作者简介 刘功桂(1968年一)男,高级工程师,主要从事电气电子产品、 材料等认证检测技术、环境适应性、可靠性和标准化研究 【4】王振尧金属材料大气腐蚀研究动态[J】全面腐蚀控制|_995 [4 4】:1—3 [5】Feliu S,Morcillo M Chico B Effect of distance from sea on atmosphericcorrosion rafe[J】.Corrosion,1999,55f9):883—891 【6】孙志华 李金桂,李牧铮.金属材料大气腐蚀加速试验研究的发展趋 势[J]材料工程,1995,(12):41—42. 【7】金营,唐其环,彭长灏等大气腐蚀的模拟加速试验方法(摘要川】. 腐蚀科学与防护技术,1995,7(3):214—215 【8】Khobaib M Chang F C,Kepple E E,et al AcceIerated atmosphe ric corrosion testing【A】A Symposium Sponsored by ASTM Committee G一1 on Corrosion of Metals【C】.Denver Colo 1980. 【91 Miranda L R,Safhier LNogueira R,et a1.Atmospheric corrosion tests in Brazilian legal Amazon,field and laboratory tests[J]. Mater Corros.,2000,51(3):182—185. 【10】Johansson L G,Vannerberg N G.The atmospheric corrosion of unprotected carbon stee1.A comparison between field study andlaboratorytest【J].Werkstoffe Und Korrosion,1981,32:265—268. 【I 1】李牧铮等周期浸润腐蚀试验方法的研究,技术资料汇编,第三机 械工业部第621研究所,198I. 作者简介 邱碧涛,工程师,主要从事设备开发与设备改造相关工作。 57 201 2年06月・环境技术 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- efsc.cn 版权所有 赣ICP备2024042792号-1

违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务