BGA热疲劳仿真报告

1、.1.1电子封装发展过程为了便于晶体管在电路中使用和焊接， 要有外壳外接引脚； 为了固定半导体芯片，要有支撑它的外壳底座； 为了防护芯片不受大气环境污染， 也为了使其坚固耐用，就必须有把芯片密封起来的外壳等， 这样产生了微电子封装技术。 概括起来，微电子封装技术经历了以下几个发展阶段：第一阶段：20世纪 60、70年代，IC芯片的制造还处于初始阶段， 集成度很低，对 IC封装没有更多的要求。此阶段采用了以双列直插 (DIP) 为主，以单列直插式(SU）与针栅阵列 (PGA)为辅的封装。特点是封装尺寸大，占印刷电路板 (PCB)面积大，集成度和频率难以提高。第二阶段：进入 80年代，出现了表面组

2、装技术 (SMT)，该技术以回流焊代替波峰焊，进一步提高了 PCB成品率，对 IC的封装提出了新要求，开发出了塑封有引线芯片载体 (PLCC)、四边引出线扁平封装 (QFP)的紧凑型封装。第三阶段： 90年代中前期， 随着个人计算机的普遍使用， 带来了计算机产业质与量上的重大变化， 原有的 PLCC、QFP、封装上的系统 (SOP)已不能满足它的发展要求，从而进一步引入了更小更薄的封装形式： 窄间距小外形封装 (SSOP)、窄间距四边引出线扁平封装 (SQFP)、引线的球栅阵列 (BGA)封装及壳系统封装 (SIP) 。特别是 BGA封装形式，使 IC引出脚大大增加。至此，多年来一直大大滞后于

3、芯片发展的微电子封装，由于 BGA的开发成功而终于能够适应芯片发展的步伐。第四阶段： IT 产业的不断繁荣，促进了电子设备向高性能、高集成、高可靠性方向发展，而支持其发展的关键技术就是 IC组装技术。封装技术进入了高速发展期，先进的封装技术和形式不断涌现，如多芯片组件技术 (MCM)、芯片尺寸封装技术 (CSP)、芯片直接贴装技术 （DCA)、晶圆规模集成技术 (WSI) 等。其中，CSP主要是由 BGA向小型化、薄型化方向发展而形成的一类崭新的封装形式。1.2电子封装分级从一个晶体管到几个集成电路板，电子封装技术可以分为以下五级：(1) 零级封装。包括芯片粘接和芯片互连技术。芯片粘接只需将芯

4、片固定安装在基板上，方法有 Au-Si 合金共熔法、 Pb-Sn合金片焊接法、导电胶粘接法、有机树脂基粘接法；芯片互连技术包括引线键合 (Wire Bonding) 、载带自动焊 (TapeAutomated Bonding) 和倒装焊腰（ Flip Chip Bonding)-种。(2) 一级封装。它是将一个或多个 IC芯片用适宜的材料 ( 金属、陶瓷、塑料或它们的组合 ) 封装起来，同时在芯片的焊区与封装的外引脚间用上述三种芯片互.专业资料 .连方法连接起来， 使之成为有实用功能的电子元器件或组件。 该级封装包括封装外壳制作在的单芯片组件 (SCM)和多芯片组（ MCM)两大类，在各个不同

5、的发展时期都有相应的封装形式(3) 二级封装。这一级封装实际上是组装， 即将上一级各种微电子封装产品、各类元器件及板上芯片 (Chip On Broad)- 同安装到 PWB或其他基板上。该级组装技术包括通孔安装技术 (THT) 、表面安装技术（ SMT）和芯片直接安装技术 (DCA)三部分。(4) 三级封装。这是一级密度更高、功能更全、更加庞大复杂的组装技术，是由二级组装的各个插板或插卡再共同插装在一个更大的母板上构成的， 是一种立体组装技术。(5) 三维 (3D) 封装。以上所述的几种微电子封装技术都是在 xy 平面实现的二维 (2D) 封装。由于电子产品小型化、轻型化、薄型化等高密度封装

6、的要求，电子封装必然要在二维的基础上向 z方向发展，形成三维封装技术。 3D封装主要包括三种方法：埋置型、有源基板型和叠层法。在零级和一级封装中，倒装焊 (FCB)技术得到了广泛应用，该技术是最有前途的一种芯片互连技术，它比 WB和TAB的综合性能高，是正在迅速发展和应用的高新技术。图 1-1图1-1 给出了印刷电路板上的塑料球栅阵列 (PBGA)、芯片尺度封装 (CSP)和直接芯片贴装 (DCA)结构示意图。 PBGA和CSP都是高密度器件封装，而 DCA则省略了器件封装这一环节。2.1 BGA 封装简介球形触点阵列 (BGA)是近几年才迅速发展起来的最新表面安装技术中的一类封装形式。它是在

7、基板的背面按阵列方式制出球形触点作为引脚，在基板正面装配集成电路，而形成的一种表面贴装技术。按基板材料的类型， BGA的封装类型主要有：塑料球栅阵列(Plastic BallGrid Array ，简称 PBGA)、陶瓷球栅阵列 (Ceramic Ball Grid Array ，简称 CBGA)、陶瓷柱栅阵列 (Ceramic Column Grid Array，简称 CCGA)、载带自动键合球珊阵.专业资料 .列 (Tape Automatic Ball Grid Array ，简称 TBGA)。BGA焊球直径有 0.89mm和 0.76mm等多种，间距通常由 1.27mm和1.Omm。这

8、种封装器件具有以下特点：1器件引脚很短， 使信号路径短，减小了引线电感和电容， 增强了电性能；2改善共面问题，再流焊时具有自对正效应；3能与原有的 SMT贴装工艺和设备兼容， 原有丝印机、 贴片机和再流焊设备都能使用；4BGA适合 MCM的封装需要，有利于实现MCM的高密度、高性能。当然， BGA本身尚具有一些不尽人意之处：1焊点隐藏在封装之下，焊点质量检测困难：2BGA的返修需要专用设备；3BGA的可靠性问题。与铜或合金引脚 ( 线) 相比，球形焊点更便宜，也非常适应SMT生产技术的要求，因此各大芯片制造公司( 如Motorola ，IBM，Hitachi ，NEC等) 纷纷将 PGA(Pi

9、nGrid Array)针栅阵列封装中的插针用焊球来代替，从而演变出CBGA陶瓷球栅阵列封装。陶瓷封装因其良好的气密性，而主要应用于军事领域。 CBGA采用了多层陶瓷基底，并使用了两种不同的焊点连接方式。一种是在小一点的陶瓷基片上，采用了球形焊点： 另一种连接方法是更大的陶瓷基片上采用柱形焊点，它现在被称为陶瓷柱形网格阵列封装(Ceramic Column Grid Array， CCGA)。图 2-1 典型 CBGA结构示意图.专业资料 .由于 BT树脂 ( 双马来酞亚胺三氮杂苯树脂 ) 或聚酰亚氨 (PI) 比陶瓷便宜，并且有很好的电性能， 因此被用来代替 CBGA中的陶瓷基底， 并用塑料

10、将基底与晶片密封起来，形成塑料球栅阵列封装 (PBGA)，应用于民用。钎料球 (Solder ball) 组成为共晶钎科 63Sn37Pb或 62Sn36Pb2Ag，目前钎料球直径分为 0.76mm、0.5mm和0.3mm几种，间距一般为 1.5mm、 1.27mm和1mm。早在 1989年， Motorola 公司就由类似 PBGA形式的封装器件成功用于便携式的电子产品中。图 2-2 典型 PBGA结构示意图图2-2 为典型的 PBGA封装结构示意图。封装时首先采用导电的芯片粘接剂将Si 芯片与 BT树脂基板粘接起来，再通过印刷铜导线和过孔将电路引至BT基片底部的焊盘阵列。 为了保护芯片不受

11、损坏， 采用模压化合物对芯片进行密封，最后将共晶钎料球与基板上的 cu焊盘 (Pad) 连接起来形成钎料凸点 (Solderbump)，这些以阵列方式排列的钎料凸点即为PBGA封装器件的 I/0 引脚。本构模型在BGA热疲劳寿命预测中的应用黏性形变是指黏性物体在剪切应力作用下发生不可逆的流动形变， 该形变随时间的增大而增大， 具有率相关的特性； 塑性是指材料在外应力去除后仍能保持部分应变的特性， 材料发生塑性形变而不断裂的能力称为延展性， 塑性变形与时间无关，具有率无关性； 蠕变是指在恒定应力作用下材料的应变随时间增加而逐渐增大的现象， 低温下表现为脆性的材料， 在高温下往往具有不同程度的蠕变

12、性为，蠕变与时间相关，具有率相关性。自然界中实际存在的材料， 其形变一般介于理想弹性固体与理想黏性液体之间，既具有固体的弹性，又具有液体的黏性，即黏弹性。黏弹性材料的力学性质与时间有关，具有力学松弛的特性，常见的力学松弛现象有蠕变、应力松弛、滞后和力损耗等。.专业资料 .焊点的寿命预测一直是焊点可靠性问题的重要容， 已经提出了多种寿命预测模型，如基于应变围的 Coffin-Manson 经验模型、基于断裂力学的裂纹扩展模型和基于损伤累积的能量模型等。 当前普遍认为， SnPb类焊点在热循环条件下的失效机制是蠕变疲劳的交互作用， SnPb焊点的失效断口既有疲劳断裂的疲劳辉纹，又具有蠕变断裂特征的

13、沿晶辉纹。 这是因为焊点的熔点较低， 在常温时就已经发生了不可忽略的蠕变。 所以单纯应用传统的疲劳分析模块是不足以预测焊点的热疲劳寿命的，在进行有限元分析时，必须考虑蠕变效应。在温度循环加载或者功率循环加载实验中， 对于 SnPb焊点钎料通常采用两种不同的理论来描述其非弹性行为。 一种是采用与时间无关的塑性形变和与时间有关的蠕变形变分开处理的模型： 另一种则透过建立统一型粘塑性本构模型来对其进行分析。与前者相比，采用统一型粘塑性本构模型避免了大量材料常数的定义，而且从连续介质热力学角度出发， 塑性形变和蠕变产生于同一机制， 即位错理论，并考虑到 SnPb基钎料具有较强的温度和加载速率的相关性，

14、 所以采用统一型粘塑性本构关系能更准确逑反映材料力学行必。 统一型粘塑性本构方程与传统塑性本构理论不同的是， 它考虑材料部状态变化对形变的影响， 导致在本构方程中引入部变量，并用准确的演化方程描述部变量。 目前，已发展了多种形式的粘塑性本构模型，其中比较有代表性的有： (1)Miller 模型； (2)Bodncr-Partom 模型；(3)Anand 模型。由于 Anand本构方程同时描述塑性和蠕变形变时准确可靠，并且在 ANSYS中可以方便地定义与调用， 因此本文采用 Anand本构模型来分析钎料的力学行为。1982年， Anand等人根据金属的高温成形采用单一变量描述材料与时闻有关及无关

15、的非弹性形变，建立了金属热成形的粘塑性本构模型 ANAND模型。和其他材料模式相比，它具有两个基本特征： (1) 在应力空间没有明确的屈服面，因此在形变过程中，不需要加载卸载准则， 塑性形交在所有非零应力条件下产生；(2) 采用单一部变量描述材料部状态对非弹性流动的阻抗。 部变量 ( 或称形变阻抗 )用 s标记，具有应力量纲。 Anand本构模型可以反映粘塑性材料与应变速率、 温度相关的形变行为，以及应变率历史效应、应变硬化和动态回复等特征。在ANSYSWorkbench中可以方便地定义 ANAND本构模型的材料参数， ANAND本构模型的相关参数一共有九个，常见的 63Sn37Pb的ANAN

16、D参数如下表所示：S0(MPa)Q/RA(s-1)h0(MPa)ms(MPa)na56.33108301.49E7112640.750.30380.4150.02311.34表2-1 63Sn37Pb由于金属的高温成形与钎料的热循环响应产生于相似的位错运动机制，具有相同的粘塑性特征，所以它在钎料力学性能的研究方面的到广泛的应用。Anand.专业资料 .模型在商用软件 ANSYS程序中可以获得，且使用对需要较小的CPU消耗 ( 如对于同一个蠕变问题， Anand模型的计算时间是幂指数型的六分之一) 。在Workbench中可以调用 Transient Structural模块仿真在热循环载荷下的

17、具有粘塑性材料参数的 BGA焊球的力学行为。图 2-3 焊球的材料参数科芬 - 曼森疲劳寿命预测模型通过寿命预测模型， 可利用有限元模拟的结果预测出封装组件的寿命。 针对封装组件来说，依照不同的破坏信息，可将寿命预测模型分成主要四类：(1) 以塑性形变为基础的预测模型。该模型着重于与时间无关的塑性效应，比较有代表性的是 Coffin-Manson 疲劳模型、Soloman疲劳模型和 Engelmaier 疲劳模型。这些模型提供了破坏循环数与每一循环焊点剪切塑性应变大小的经验关系；(2) 以蠕变形变为基础的预测模型。该模型考虑与时间相关的效应，蠕变形变是影响焊点寿命的重要因素， 其机制相当复杂，

18、 至今仍无模型能完全预测其整个过程。比较有代表性的有 Kencht-Fox 模型和 Syed模型，但因忽略塑性形变， 故无法准确预测所有的产品；(3) 以能量为基础的预测模型。 该模型较新，考虑到应力与应变的迟滞能量，通过计算每一次循环的应变能或者应变能密度， 从而得出与焊点寿命之间的联系。比较有代表性的有 Akayi 模型和 Darveaux模型；其中， Darveaux模型在 ABAQUS仿真分析中应用较多；.专业资料 .(4) 以断裂参量为基础的预测模型。该模型以断裂力学为基础，计算裂纹的扩展，累积其过程所造成的破坏效应， 已经在表征工程材料的弹塑性断裂和疲劳中得到了成功的应用。 在焊点

19、的寿命研究方面， 具有代表性的有 Lau和 Pao等人的模型。电子封装 SnPb焊点失效是低周疲劳失效。 目前，焊点热循环失效的寿命模式主要以低周疲劳的 Coffin-Manson 方程 ( 简称 C-M方程 ) 为基础，即材料的低周疲劳寿命 (Nf) 和塑性应变围 ( ) 之间符合如下经验关系：（2-1 ）式中，f 是疲劳韧性系数； c是疲劳韧性指数。 该方程是研究材料低周应变疲劳的重要公式。Solomon采用 C-M方程研究了 Sn60Pb40钎料的等温疲劳性能，得到在-55 +125围，钎料的疲劳寿命和剪切塑性应变围间关系为：（2-2 ）由于 SnPb焊点热循环失效的影响因素很多，在具体

20、应用中， C-M方程被修正成多种形式。例如： Engelmaier 同时考虑热循环温度和频率的影响， C-M方程被修正为：（2-3 ）其中， c为疲劳韧性指数，与频率和温度有关：（2-4 ）Tm为热循环平均温度，单位为；f 为热循环频率，单位为 cycle/day ；3.1基于 Workbench 的 BGA疲劳寿命的仿真预测常见的非线性分析分为三种： 几何大变形， 材料非线性与接触非线性 （状态非线性）。由于考虑到粘塑性问题，所以此类分析属于材料非线性，同时由于计算了 5个温度循环周期， 总时间 10080秒，所以这个问题属于瞬态动力学畴， 综合来看，此类问题属于非线性瞬态动力学问题。.专业

21、资料 .关于非线性需要注意的问题：（1）打开大变形选项，允许程序考虑大变形效应；（2）ANSYS解决非线性问题采用的是牛顿- 拉夫逊迭代算法：ANSYS求解非线性问题的思路是将问题分解成多个载荷步，在每个载荷步将非线性问题转化为线性问题， 并且通过线性求解的结果来调整刚度子矩阵，因此，非线性问题在程序上的表现就是模型的刚度矩阵不再是固定不变的，而是与空间和时间相关的， 而程序的任务就是调整刚度矩阵，使其在每一个载荷步接近该载荷步刚度矩阵的真实值。图3-1 牛顿 - 拉普森迭代示意但是，纯粹的增量近似不可避免地要随着每一个载荷增量积累误差， 导种结果最终失去平衡，产生非常大的误差。 ANSYS程

22、序通过使用牛顿拉普森平衡迭代克服了这种困难， 它迫使在每一个载荷增量的末端解达到平衡收敛 （在某个容限围），在程序上的实现方法是利用切向刚度矩阵不断迭代调整当前刚度矩阵，将每一次迭代计算产生模型力与外力相减， 得到一个残余力的向量， 并对这个向量求数，这个数就是迭代残差，迭代残差接近 0，就代表这力与外力平衡，此时刚度矩阵最接近于真实值。在非线性仿真时经常遇到不收敛的情况， 此时选择输出最近几步的牛顿残差，观察哪些部分牛顿迭代残差较大，就可以知道是哪里出现了问题，以修改模型。（3）关于接触刚度有时不适合的接触刚度会导致迭代不收敛， 如下图所示，在进行非线性求解时发生了不收敛， 检查牛顿残差时发

23、现在接触区域出现了较大的值， 此时可以断定是接触出了问题，所以可以通过降低接触刚度来解决。 接触刚度的默认值为 1，最小可以调到 0.0001 ，而接触形式选择纯惩罚形式也会改善收敛性。.专业资料 .图3-2 接触刚度过高（4）合理控制时间步与网格密度在进行非线性仿真时还容易出现以下提示信息：An Internal SolutionMagnitude LimitWasExceeded.或者 The solverengine was unable to convergeon a solution for the nonlinear problem as constrained.出现这两种错误的原

24、因有以下几种是采用了不适当的约束或者过大的载荷，在BGA分析中出现此类问题很可能是由于在某些时间步中载荷变化过大，比如在温度上升下降的部分如果分的时间步长较大， 就容易出现这种情况， 这时可以在需要细化时间步的位置单独设置一个载荷步，在此载荷步设置较多的时间步数就可以解决这个问题；另外如果在一些关键部位网格划分过于粗糙也会导致结果不收敛，此时需要考虑细化网格以改善收敛性。（5）建模错误导致的特殊情况图3-3建模错误导致的特殊情况在图 3-3 所示图中，给出了当迭代不收敛时的牛顿残差分布，发现在芯片上出现了较大的迭代残余力， 一般来说应该是在芯片表面与塑封体的接触上出了问.专业资料 .塑封体焊球

25、PCB.题，但是经过检查后发现，导致结果不收敛的原因是模型存在错误：焊球与PCB板在建模时由于装配不当存在空隙，导致基板、芯片和塑封体“悬空”而自由热膨胀。所以当结果不收敛检查牛顿迭代残余力时要注意不一定是残余力大的部位导致不收敛，错误原因也可能出现在其他部位。在本模型中，当去除了焊球与 PCB板间的间隙之后，求解就收敛了。3.2 BGA 热疲劳仿真结果内部芯片基板图 3-4 BGA模型示意图焊球的 Anand材料参数如前所述，其氏模量与泊松比定义为与温度相关的材料参数。模型中其他材料参数如表 3-1 所示：表3-1材料参数名称氏模量（ GPa）泊松比热膨胀系数（ -1 ）基板220.281.

26、8E-5焊球上下铜垫片1080.3251.67E-5芯片1310.32.8E-6塑封体260.37E-6PCB板220.281.8E-5仿真的循环温度载荷如图3-5 所示：图3-5模型的循环温度载荷条件仿真的总体位移与 von-Mises 应力结果如图 3-6 所示：.专业资料 .图3-6 总体位移云图与 von-Mises 应力云图从位移云图可以得知越远离芯片中心， 基板、塑封体和焊球的位移越大， 基板翘曲越厉害，由此我们可以得知越远离芯片中心焊球的受力情况越恶劣， 所以，从焊球阵列的角度来看， 焊球阵列中首先出现热失效的焊球应该位于焊球阵列的外围。从von-Mises 应力云图来看，焊球上

27、应力最大的部位位于焊球的顶端和低端，所以，从单个焊球的角度来看， 焊球上首先出现裂纹的地方应该是在焊球与铜垫的结合部。所以综合以上两点，在远离芯片中心的焊球的上下端是芯片热失效的危险部位，这与实验结果和他人仿真结果是一致的。图 3-7 位移历史曲线图3-7 给出了模型的位移历史曲线，从图中可以看出，随着时间的推移、温度循环载荷的施加， 在每一次循环中位移也出现了循环的现象，而且随着循环次数的增加，最大位移的峰值与谷值也在增加，这是由于在一次次的循环过程中累积了粘塑性变形的结果。 位移有两部分组成： 弹性变形与粘塑性变形， 弹性变形在每次循环中是固定的，而粘塑性变形是不断增加的。相应的，我们可以

28、得到每一次循环过程中的塑性应变， 而根据塑性应变我们可以得到 BGA的热疲劳寿命。.专业资料 .图 3-8模型的塑性应变图3-8 给出了模型的塑性应变，图中只有焊球存在塑性应变，而其他部件没有显示，这是因为只有焊球定义了粘塑性材料，而其他部件只定义了弹性材料，所以没有塑性变形。 并且在图中可以观察到， 最大的塑性应出现在焊球阵列外围的焊球的上下端，这就是热疲劳失效的危险部位，这与前面的结论是一致的。图 3-9 Report 中的塑性应变图3-9 给出了 ANSYS自动生成的报告中塑性应变的历史数据，随着时间的推移，每个周期的塑性应变围会趋于稳定，而塑性应变围正是预测热疲劳寿命的关键数据。报告中

29、给出了每个时间步的塑性应变的详细数据值，可以很方便的得到塑性应变围。在此要注意， 科芬曼森疲劳寿命预测模型中要求给出的是剪切塑性应变围，而此处求出的是塑性应变围，需要做一下后处理。最后根据科芬曼森疲劳寿命预测模型求出BGA的剪切塑性应变围为 0.034977 ，热疲劳寿命是 767周，这比较接近中南大学林丹华的仿真结果，他是在 -55 到 125.专业资料 .之间循环，他仿真的结果是剪切塑性应变围为0.035161 ，寿命为 744周。我们计算的误差为 3.09%。（参见中南大学林丹华硕士学位论文PBGA封装热可靠性分析及结构优化一文）3.3 基板材料与 PCB板材料的热配合对热寿命的影响在能够预测 BGA的疲劳寿命之后，所做的进一步工作是研究基板和 PCB板的材料选择对焊球会产生什么样的影响。 限于时间原因， 这一部分没有完全做完， 但是也得到了一些有用的结论。根据之前的仿真经验，可以假设基板与PCB板对焊球的影响主要取决于热配合，当基板与PCB板的热配合性较好时也就是热膨胀系数较为接近时，那么基板与PCB板在热膨胀时能够基本保持一致，使焊球受到的剪切力较小；另外当基板与PCB板的热膨胀系数都比较小时，也应该