NCKU复合材料实验室

1、1,利用有限元素與田口方法探討FCCSP構裝無鉛錫球之最佳化疲勞壽命 Investigation of the Fatigue Life of Lead-Free Solders for Flip Chip Chip Scale Package by Finite Element and Taguchi Method,南台科技大學 副教授,2,報 告 流 程,緒論 理論基礎 模型建立與評估錫球疲勞壽命 以Surface Evolver預測錫球形狀 田口氏品質工程方法 結論 未來研究方向,3,前 言,晶粒尺寸封裝(Chip Scale Package，CSP) ，泛指以各種方式封裝後的IC，若封

2、裝體邊長比內含的晶片邊長，大20%以內，或封裝體的面積內含晶片面積的1.5倍以內，都可稱之為CSP封裝。 由於覆晶封裝具有良好電氣特性、高I/O接點密度，且能縮小IC尺寸增加每片晶圓產出，為未來極具潛力之構裝方式。 錫球在覆晶晶粒尺寸封裝(FCCSP)構裝中之功能主要包含有傳電、導熱及吸收上下元件之膨脹差等，因此錫球的可靠度佔有極為重要之影響。 錫球在室溫下已超過其熔點的一半，易有潛變現象。,4,研究動機與目的,過去的封裝型態，多半以打線作為內部接合，但打線製程花費時間相當長，成為接合技術的最大瓶頸。許多高階、可攜式電子產品，需要使用高腳數、散熱性較佳，或較為輕薄短小的封裝，則需使用覆晶技術。

3、 本文採用氧化鋁陶瓷基板，由於緻密性高，對水分子滲透有優良的阻絕能力、強度良好、散熱性佳、與耐高壓性等優點， 故被廣泛應用在需求高可靠度的IC構裝中。 半導體產業唯有靠不斷地創新及研發，才能在激烈的競爭環境中繼續保有領先之地位，期望透過本研究能對業界有所幫助。,5,研 究 之 方 法,本文以96.5Sn3.5Ag無鉛錫球與氧化鋁陶瓷基板之FCCSP封裝模式，建構在FR-4電路板上。 先利用Surface Evolver預測錫球迴焊後之形狀，再引入ANSYS 7.0有限元素分析軟體建立三維模型。 採用葛拉佛拉-阿瑞尼阿斯潛變模式(Garofalo-Arrhenius Creep) ，透過3-D幾

4、何模型的建立、網格分割及計算求解等步驟，進行分析模擬。 探討不同材質的封裝結構在-20110溫度循環下，錫球之疲勞壽命、應力應變及遲滯曲線等機械行為之變化情形。 以等效潛變應變範圍代入Modified Coffin-Manson計算公式，預估錫球之疲勞壽命。,6,覆晶封裝製程,晶圓進入機台前需確認是否受污染。 以電鍍或印刷植球進行銲錫(或金塊)凸塊製程。 之後需經迴焊(Reflow)製程，使錫球成型。 凸塊完成後進行晶粒切割。 然後由吸嘴吸住晶粒背面沾上助焊膏(Flux Paste)，暫時將晶粒固定在基板上。 置放完成後，進入迴焊爐內形成接合點。助焊膏的殘餘物須以清洗劑清除。 此即所謂的C4

5、(Controlled Collapse Chip Connection) ，此一技術取代了傳統的打線接合。,7,可 控 塌 陷 晶 片 (F4) 連 接 技 術,8,覆 晶 的 主 要 優 點,可降低晶片與基板間的電子訊號傳輸距離，適用在高速元件的封裝。 可縮小晶片封裝後的尺寸，使得晶片封裝前後大小差不多。,9,本 文 討 論 之 變 形 理 論,本文所採用之錫球為錫銀(96.5Sn3.5Ag)之材料。 在高溫且恆溫狀態、穩態負載條件下具有潛變效應。 同時承受反覆熱循環負載，使得材料進入塑性範疇。 本文即討論潛變與塑性變形。,10,文 獻 回 顧,1996年S.M. Heinrich等以63

6、Sn37Pb共晶錫球為分析模型。假設熔融狀態之錫球表面為以一等曲率圓弧，推導出一顯函數解析解，可精確地預測錫球於迴焊過程後之形狀。 2000年Mertol等以有限元素法，模擬晶圓級封裝在溫度循環下應力、應變與變形，並使用田口式方法，探討設計參數對錫球的疲勞壽命之影響。 2001年Pang 以彈性、塑性、潛變之材料，探討覆晶構裝模型在溫度循環下，錫球之疲勞壽命。,11,材料在室溫環境下發生潛變,1.初始(暫態)潛變。 2.穩態潛變，潛變過程時間最長，最主要的部份。 3.加速潛變。 因此通常以穩態潛變方程式描述潛變行為。 本文以ANSYS 7.0對96.5Sn3.5Ag無鉛錫球，以葛拉佛拉-阿瑞尼

7、阿斯(Garofalo-Arrhenius)潛變模式進行分析。以探討構裝體在承受-20至110之熱循環負載時，所發生之錫球破壞行為。,12,葛拉佛拉-阿瑞尼阿斯潛變模式,註1：單位採用MPa；T 單位採用絕對溫度K。 註2： C1 及C2 之參數值為溫度T 之函數，其值會隨溫度而變化，而此處T是以絕對溫度K為單位。,註1,C1 = 18(553-T)/T，單位:1/sec C2 = 1/(43.99-0.079T)，單位:1/MPa C3 = 5.5 C4 = 5802，單位:k,註2,13,Tresca 與 von Mises 降伏準則,目前常用的兩種延性材料之破壞準則為特雷斯卡準則(Tre

8、sca Yield Criterion)及米澤斯降伏準則(von Mises Yield Criterion) 。 由於特雷斯卡準則太過於保守，故本文選用米澤斯降伏準則為依據。,14,塑 性 行 為 模 式,本文選用德國von Mises所提出降伏準則為依據。 依照此準則，當材料中每單位體積扭曲能量之最大值達到材料於拉伸試驗下所發生破壞之扭曲能量時，材料即開始破壞。 米澤斯破壞準則在卡式座標分量下可表示為 當式左項等於或大於右項時，則材料產生破壞，左項為材料之等效應變，右項為材料拉伸試驗下所發生之破壞應變。,15,多 線 性 等 向 性 硬 化 法 則 (Multilinear Isotrop

9、ic Hardening),錫球承受一循環熱應力，而發生降伏進入塑性行為區域時，降伏面會隨著熱循環中塑性應變之增加而改變 。 本文採用多線性等向硬化模式。 理論上在某一溫度下，材料塑性行為應力與應變之關係，應為一光滑曲線，多線性等向硬化模式是以數個線段趨近此曲線 。 等向硬化是指材料在進入塑性變形以後，加載曲面在各方向呈現均勻向外擴張之現象。 此模式中也假設材料受反方向壓縮負載降伏應力值等於拉伸降伏應力值，忽略包新格(Bauschinger)效應。,16,低 循 環 疲 勞 壽 命,結構體承受反覆負載時，將形成循環應變與應力，即使應力小於材料本身的降服強度，亦會導致結構整體的破壞，稱之為疲勞破

10、壞（Fatigue）。 早期的研究著重高循環疲勞(N104)的探討，但低循環疲勞(N104)在工程上的應用也很重要，電子構裝的疲勞的破壞即為低循環疲勞破壞。 其破壞位置通常發生於錫球最大的應變範圍處而非最大應變處。 最常見的疲勞模型為Coffin-Manson公式。,17,Coffin-Manson疲勞壽命估算公式 (共晶錫球),18,Modified Coffin-Manson疲勞壽命估算公式 (無鉛錫球96.5Sn3.5Ag),0.325, -0.41299,24 cycles/day,),(,(-20+110) 45,19,介 紹 Surface Evolver,德國數學教授Ken Br

11、akke在1989年首次寫出Surface Evolver程式軟體 。 Surface Evolver是一套利用能量法(Energy-Based Method)來分析液滴表面形狀的程式軟體。 利用一系列的小三角形面積元素，以能量梯度下降法來計算熔融錫球表面張力的能量、重力位能及錫球固化體積改變所產生的外力能量，進而改變元素位置，使表面總能量達到平衡狀況的最小值而模擬出固化形狀。 文中錫球固化形狀是由熔融錫球的表面張力及上下墊片之濕潤現象所控制，因此Surface Evolver適合模擬錫球迴銲過程後的形狀。,20,Surface Evolver 模 型 建 立,錫球基本假設條件 錫球材料為均質

12、等向性。 固化時錫球為靜態平衡。 錫球固化時，錫球墊為圓形，且完全對準。 融熔狀態的錫球，其表面輪廓為軸對稱。 錫球表面的經線方向為圓弧線。 錫球墊和錫球是完全接觸。 錫球和銲墊為完全黏著。 迴焊前後錫球之體積不變。,21,Surface Evolver 分 析 流 程,先確立研究主題為罩幕定界(Solder Mask Defined, SMD)錫球。並建構錫球之幾何分析模型。 指定輸入參數包括錫球體積、錫球密度、表面張力、外加負載、初始高度及錫球墊片半徑等因子。 當執行Surface Evolver程式完畢之後，需檢查所得到之恢復力F值是否與外加負載F值(重力)一致。 當恢復力與外加負載不一

13、致時，則必須修正所輸入的錫球初始高度，一直至恢復力等於外加負載為止，則此時的高度即是錫球迴焊過程後之高度。,22,分 析 流 程,錫球之中心軸及座標原點 (上),分析錫球形狀前先給定一初始形狀 (右上),罩幕定界錫球於迴焊過程中之受力狀況 (右下),23,分析構裝體模型圖,本文FCCSP分析模型為7.0mm X 7.0mm X 0.41mm之晶片，四週佈植有100個金凸塊。 下接7.45 mmX 7.45mm X 0.25mm Al2O3 陶瓷基板。 在晶片與基板間填膠，基板下面為1010陣列之96.5Sn3.5Ag無鉛錫球，其接墊之直徑為0.3mm，間距為0.65mm 。 錫球下接8.0mm

14、 X 8.0 mmX1.0 mm玻璃樹脂(Epoxy Glass) 所製成之FR-4印刷電路板。,24,FCCP 模 型 圖,25,構裝體分析之基本假設,構裝體所有材料皆為均質等向性，受拉與受壓行為相同。 模型中所有界面皆為完全接觸(Perfectly Adhesive)。 模型內部的溫度與外界環境溫度相同，且假設溫度場不隨空間變化，即在瞬間皆已達平衡狀態，即(T(x,t) = T(t)。 模型在初始狀態25時為應力自由狀態(Stress Free)，忽略製程中所產生之殘留應力(Residual Stress)。 忽略錫球經過迴焊過程後所引起之材質變化。 錫球之材料性質為潛變模式，其餘部分皆為

15、彈性材料。,26,構 裝 體 八 分 之 一 對 稱 模 型,為節省電腦模擬計算時間，由兩對角線的斷面上，截取分析結構右上半段，即八分之一對稱模型進行分析 。 邊界條件設定為Y軸上所有節點為對稱邊界條件。 在模型結構左下端點（原點）設為固定點。 以使分析模型受力變形後能自由翹曲，其餘邊界皆假設為自由端。 兩對角線之兩面設為對稱面。,27,FCCSP構裝體八分之一對稱分析模型圖,對稱面,對稱面,28,溫 度 循 環 負 載,施予構裝體-20低溫至110高溫之溫度循環負載。 初始溫度由25開始，並假設為無應力狀態，在588秒內升至110。 每一溫度循環為60分鐘，從低溫升至高溫（Up-Ramp）需

16、時15分鐘，並在高溫區維持20分鐘之恆溫(Dwell) 。 再以15分鐘降溫至-20(Down-Ramp)，然後再維持10分鐘之恆溫(Dwell)，至此完成一個循環。 共進行五次溫度循環負載，計須18,588秒 。,29,溫度循環負載曲線,15min 20min 15min 10min 5cycles,溫度(),時間（秒）,30,以下分四主題討論,模型之網格收斂分析 模型在熱循環負載之探討 五種模型疲勞壽命之比較 田 口 方 法,31,三 維 條 狀 模 型 網 格 分 割,有限元素模型切割的網格密度，會影響數值收斂趨勢。但網格過多時，電腦運算時間也相對變長，因此必須找到適當的網格元素數目。

17、為了探討有限元素模型之收斂性，分別進行五種不同網格的切割配置。 配置方法主要於模型Y方向，採取不同線段切割數目以改善模型收斂性。,32,於模型Y方向線段切割數目,33,最外側錫球中心點於1788秒高溫與2688秒低溫下之各分佈圖,構裝體之變形分佈圖,最外側錫球等效應力分佈圖,最外側錫球等效應變分佈圖,34,網 格 元 素 數 目,高溫時網格元素數目與構裝體等效應變之關係圖,低溫時網格元素數目與構裝體等效應變之關係圖,35,模 型 之 網 格 收 斂 分 析,五個案例分析結果之構裝體變形，及錫球之等效應力與等效應變之分佈如圖所示 (Case 4)。 最大翹曲位置發生在構裝體最遠處，即印刷電路板對

18、角線上方處。 最大等效應變發生在構裝體最外側錫球處。 由圖可知元素數目為45487個(Case 4)時已趨近於收斂，即能得到相當精確的結果。,36,模型在熱循環負載之探討,構裝體最大翹曲與時間關係圖。 最外側錫球等效應力與時間關係圖。 最外側錫球等效應變與時間關係圖。 而最外側錫球中心點剪應力與時間關係圖。 最外側錫球中心點剪應變與時間關係圖。 由最外側錫球中心點剪應力與剪應變之遲滯曲線，可知當達到第五個溫度循環時已收斂良好 。,37,模型在熱循環負載之關係圖,最外側錫球等效應力與時間關係圖,最外側錫球等效應變與時間關係圖,38,最 外 側 錫 球 中 心 點 各 種 關 係 圖,最外側錫球中

19、心點剪應力與時間關係圖 (上),最外側錫球中心點剪應變與時間關係圖 (右上),最外側錫球中心點剪應力與剪應變之遲滯曲線圖(右下),39,五 種 FCCSP 構 裝 體 模 型 計算外側錫球之疲勞壽命,原設計模型 模型基板為BT-Resin 模型基板為High- Ceramics 模型封膠為Mold Resin (Epoxy Molding Compound) 模型封膠為Potting Resin等五種模型 其材料機械性質如表所示 並利用ANSYS進行模擬，計算在溫度循環負載下，最外側錫球之疲勞壽命。,40,五 種 構 裝 體 分 析 模 型,原設計模型,模型封膠為Mold Resin,模型封膠

20、為Potting Resin,基板為High- Ceramics模型,基板為BT-Resin模型,41,各元件材料機械性質,註1：無鉛錫球96.5Sn3.5Ag之楊氏模數在-40、25、50及125下分別為47,200 MPa、29,525 MPa、24,700 MPa及16,850 MPa。 註2：機械性質為溫度T之函數，其值會隨溫度而變化，而T是以為單位。,42,高溫時構裝體之變形分佈圖 (上),高溫時最外側錫球等效應力分佈圖 (右上),高溫時最外側錫球等效應變分佈圖 (右下),基板為TResin高溫時分佈圖,43,低溫時構裝體之變形分佈圖 (上),低溫時最外側錫球等效應力分佈圖 (右上)

21、,低溫時最外側錫球等效應變分佈圖 (右下),基板為TResin低溫時分佈圖,44,基板為TResin各種關係圖,45,五種模型疲勞壽命之比較,46,五種模型疲勞壽命之比較,原設計模型之基板為Al2O3陶瓷，其楊氏模數比BT-Resin與High- Ceramics大，而熱膨脹係數卻較小，比較不易熱變形，與其他元件之不匹配較大，故其疲勞壽命比模型基板為BT-Resin與High- Ceramics小很多。 BT-Resin的熱膨脹係數比High- Ceramics的熱膨脹係數大，但楊氏模數卻小很多，因此不匹配較小而其疲勞壽命略大。 至於模型封膠為Mold Resin與Potting Resin，

22、其疲勞壽命為最小，由於此兩種模型比前三者多出上面的封膠，故與其他元件不匹配更大，使得疲勞壽命最小。,47,Taguchi Method 找出最佳之水準組合 提升構裝體之可靠度,選擇以L18（21x37）直交表，進行田口實驗配置。 再利用ANSYS 7.0進行模擬，並預估錫球之疲勞壽命。 最後根據模擬結果，求得對疲勞壽命影響最顯著之控制因子，並找出最佳之水準組合，以助提升構裝體之可靠度。,48,田 口 方 法 實 驗 設 計,選定品質特性(Quality Characteristic) :選定最外側錫球的疲勞壽命為品質特性，探討各控制因子對其之影響。 判定品質特性之理想機能(Ideal Func

23、tion):品質特性屬於望大特性，其理想機能是無限大。 列出影響品質特性之因子(Factors) :構裝體外形尺寸，環境因素與構裝體材料性質。 決定主要控制因子(Control Factors)及水準(Levels):本文只考慮基板厚度、錫球墊半徑、楊氏模數與熱膨脹係數對FCCSP構裝體可靠度之影響，選出下列八個設計參數作為控制因子。,49,控制因子及其水準表,50,選定適當之直交表(Orthogonal Arrays),直交表之作用在於以最少之實驗次數獲得最佳之分析，其各項因子水準之自由度為1（2-1）+7（3-1）15 因此選用自由度不小於15且實驗次數最少之L18（21x37）。,51,

24、L18 直 交 表,52,十八組模型等效潛變應變範圍 與最外側錫球之疲勞壽命值,本文之三維模型即利用在五次溫度循環負載下，對最外側錫球所求出之累積等效潛變應變(等效潛變應變範圍值) 。 代入Modified Coffin-Manson計算公式，求得錫球疲勞壽命以預測十八組錫球封裝結構所能承受的壽命次數，藉以評估各封裝體之可靠度趨勢。 十八組 FCCSP封裝結構於最外側錫球之疲勞壽命值如下表所示，其中以第一組有最大值，疲勞壽命為3098次循環。,53,封 裝 結 構 於 最 外 側 錫 球 之 疲 勞 壽 命,54,目 標 函 數 的 S/N 比,對品質特性為構裝體之疲勞壽命而言，理想機能越大越

25、好，即是望大特性。 可透過目標函數的S/N比，判斷各因子對目標函數的影響大小。即各控制因子的水準中，具有較大的S/N比者，表示此水準對構裝體之疲勞壽命值有較好的影響。 望大特性的S/N比公式如下:,55,三 維 模 型 實 驗 結 果,56,三維模型S/N比反應圖,57,三 維 模 型 S/N 比 反 應 表,58,控制因子最佳組合,59,三 維 模 型 變 異 分 析,變異分析(Analysis of Variance, ANOVA) 目的是評估實驗誤差。 若控制因子影響的變異量遠大於實驗誤差所產生的變異量，則可判斷此因子效應是重要的，非實驗誤差所造成的。 相反地，其它無影響力的因子其效應可

26、視為實驗誤差所造成的偶發效應。,60,三維模型各因子變異量與總變異量的平方和,61,第 一 次 誤 差 統 合,進行第一次誤差統合。因其它項的變異量最小，所以將其併入實驗誤差項，並計算其它因子的信心水準。在此過程中，選擇信心水準至少需99以上。 控制因子C底膠(Underfill)之楊氏模數E可視為無影響力的因子。亦即，因控制因子C所造成的變異可視為因實驗誤差所造成的偶發現象，將於下一過程統合至誤差項。 控制因子H(FR4之熱膨脹係數)之信心水準最大為99.99，可視為最重要的因子。,62,三 維 模 型 第 一 次 誤 差 統 合,63,第 二 次 誤 差 統 合,進行第二次誤差統合。將計算

27、其它因子的信心水準。在此過程中，選擇信心水準至少需99以上。 可把控制因子H(FR4之熱膨脹係數)視為最具影響力的因子。 其次為控制因子F(基板之熱膨脹係數)。,64,三 維 模 型 第 二 次 誤 差 統 合,65,原始製程參數設計與最佳製程參數設計之S/N預測值及S/N確認實驗值之比較,依據表中之數據可製作原始製程參數設計與最佳製程參數設計在99的信心水準下。 由S/N預測值及S/N確認實驗值之比較圖，可觀察出，預測值與確認實驗值二者有相互重疊的部份，因此可以認定預測值與確認實驗值是相當接近的，因子效應的估計是可靠的 。 原始製程參數設計之構裝體之疲勞壽命值是49次。 最佳製程參數設計之構

28、裝體之疲勞壽命值是3185次，疲勞壽命值約提昇65倍，對構裝體可靠度有明顯之改善。,66,三維模型原始與最佳設計之S/N預測值及確認實驗值比較,67,三維模型原始製程參數設計預測值與確認實驗值的信心區間,68,三維模型最佳製程參數設計預測值與確認實驗值的信心區間,69,預測值與確認實驗值,由圖顯示，預測值與確認實驗值二者有相互重疊的部份，因此可以認定預測值與確認實驗值相當接近， 即因子效應的估計是可靠的。,70,結 論 (一),構裝體最大翹曲位置發生在構裝體最遠處，即印刷電路板對角線上方處，而最大等效應變發生在構裝體最外側錫球處。 構裝體網格分割元素數目為45487個(Case 4)時已趨近於收斂，因此元素數目大約為45487個時即能得到相當精確的結果，為本文使用之元素數目，而節點數目為52892個。 構裝模型在熱循環負載中，由外側錫球中心點剪應力與剪應變畫成遲滯曲線圖，可知在第五個溫度循環時已收斂良好，故在以後之分析只施與構裝體五個溫度循環，即可得合理之結果。,71,結 論 (二),原設計模型之基板為Al2O3陶瓷，其楊氏模數比BT-Resin與High- Ceramics大，而CTE比BT-Resin與High- Ceramics小，比較不易熱變形，引起與其他元件之不匹配(Mismatch)