周嘉誠，劉 芳

?

基于有限元的車載電路板多場熱應力分析

周嘉誠，劉 芳*

（武漢紡織大學 機械工程與自動化學院，湖北 武漢 430073）

主要研究汽車在行駛中發(fā)動機發(fā)熱和芯片自身發(fā)熱對車載電路板組件產(chǎn)生的影響。運用有限元軟件ANSYS對發(fā)動機模塊電路板建模，分別施加芯片自身發(fā)熱、發(fā)動機工作溫度以及兩者共同作用三種情形下的熱載荷，進行溫度場熱應力分析，比較三種情況下的溫度場分布和熱變形情況。結(jié)果顯示：在實際車載工作環(huán)境中，發(fā)動機發(fā)熱是車載電路板產(chǎn)生熱變形的主要因素。電路板組件中最大離面位移出現(xiàn)在電路板中間一帶芯片處，上述芯片位置是電路板組件在芯片和發(fā)動機共同發(fā)熱條件下的最危險區(qū)域。

車載電路板；熱載荷；芯片；熱應力

0 引言

車載電子與機載電子在服役過程中，不可避免地會遭受高溫、振動、潮濕等惡劣環(huán)境的影響。根據(jù)美國Air Force Avionics Integrity Program的統(tǒng)計數(shù)據(jù)所列舉出的影響電子設(shè)備可靠性的主要環(huán)境因素中，溫度和受熱所占的比例高達55%[1]。隨著微電子技術(shù)的進步，也使電路板的芯片集成度不斷提高，芯片的功耗和發(fā)熱量也隨之增大。同時，在車載機載環(huán)境中，電子設(shè)備還會受到由發(fā)動機等動力部件工作發(fā)熱的影響[2]。由于材料各部分之間膨脹不均勻而彼此相互約束，熱效應不可避免的產(chǎn)生熱應力[3]。因此，為降低電路板的熱應力，提高其抗變形能力，改善產(chǎn)品的熱性能以及提高可靠性，分析電路板各溫度場特性及其熱變形必不可少。

目前，已經(jīng)有不少的學者利用有限元軟件對電路板進行了受熱和溫度分析[4-8]，尤其是在航空航天的電子設(shè)備的熱分析上，已經(jīng)做出了較為深入的研究，這些仿真模擬和實驗得出的結(jié)論，為電路板后續(xù)研究提供了寶貴的數(shù)據(jù)經(jīng)驗。但是，對于我們?nèi)粘Ｊ褂闷囍械碾娮宇惍a(chǎn)品，目前該方面的研究還比較少。因此，模擬出符合車載環(huán)境下的電路板的發(fā)熱和受熱條件，對在條件下的電路板進行熱應力分析，變得尤為重要。

本文研究對象為汽車車載電子設(shè)備中的印刷電路板。利用有限元軟件ANSYS建立模型，結(jié)合電路板上芯片發(fā)熱與汽車正常工作下的發(fā)熱情況對電路板組件施加熱載荷，分別模擬（1）僅受電路板上芯片發(fā)熱影響；（2）僅受汽車發(fā)動機工作溫度影響；（3）受電路板芯片發(fā)熱與發(fā)動機工作溫度共同影響等三種情況。首先對上述三種情況的電路板組件進行溫度場分析，其后根據(jù)實際情況施加約束條件并展開應力分析。最后觀察和比較的三種情況下電路板組件的溫度場分布與熱變形，找到電路板組件在該溫度場條件下的最危險位置。對于車載環(huán)境下的印制電路板受熱條件和結(jié)構(gòu)變形的仿真模擬，對今后類似情況下的設(shè)計提供參考。

1 研究對象與有限元模型

本文中有限元模型實體來源于某品牌轎車發(fā)動機模塊中的一部分，位于發(fā)動機和駕駛室之間。如圖1所示，印制電路板長173mm，寬112mm，厚度2mm，電路板上主要有兩個邊長分別為15mm和22mm的正方形芯片，和若干長方形小芯片。具體電路板組件模型和芯片的編號如圖2所示。該電路板由印制電路板基板、硅芯片、電容和電線能元件構(gòu)成。其中主要的發(fā)熱和受熱影響的是芯片，電容和電線因其功耗較小，故忽略不計。為了提高分析計算的效率，在有限元建模時將電路板中功率較大，發(fā)熱較高的芯片保留下來。同時，由于受熱時間不長且非熱循環(huán)載荷，故此處不考慮蠕變效應對電路板造成的影響。

圖1 印刷電路板

圖2 有限元模型

根據(jù)上訴條件將模型簡化，建立有限元模型。模型劃分的網(wǎng)格質(zhì)量不僅關(guān)系到有限元計算的效率，而且關(guān)系到分析結(jié)果的準確性。網(wǎng)格劃分后，模型共生成113836個單元，144748個節(jié)點。進行熱應力分析時，可通過熱單元向結(jié)構(gòu)單元轉(zhuǎn)換[9]。電路板及組件的材料性能參數(shù)如表1所示。

表1 電路板及組件材料參數(shù)

2 PCB溫度場下的熱應力分析

利用ANSYS有限元軟件分析時候，采用時序間接耦合的方式，按照相關(guān)場的先后循序進行分析，通過把第一次場分析的結(jié)果作為第二次場分析的載荷實現(xiàn)兩種場的耦合。在仿真分析前，定義分析的單元類型、材料屬性、網(wǎng)格劃分、載荷和約束等參數(shù)。上述條件設(shè)定后，先進行溫度場的熱分析，然后把熱分析結(jié)果數(shù)據(jù)當作結(jié)構(gòu)力學分析的載荷，完成整個溫度場熱應力的分析過程[10]。

2.1 三種不同受熱情況下的PCB溫度場有限元分析

根據(jù)實際的汽車內(nèi)部實際工作環(huán)境，環(huán)境初始溫度設(shè)為25℃，空氣自然對流系數(shù)為20W/(m2·℃)。對車載電路板進行分別模擬（1）僅受電路板上芯片發(fā)熱影響；（2）僅受汽車發(fā)動機工作溫度影響；（3）受電路板芯片發(fā)熱與發(fā)動機工作溫度共同影響等三種情況下的溫度場分析，其分析結(jié)果如圖3所示。

從圖3中可以看出三種受熱情況下PCB溫度分布情況。圖3(a)中顯示了僅受芯片發(fā)熱的電路板組件溫度分布圖，由于芯片自身發(fā)熱，電路板組件中，芯片位置溫度高于電路板其他位置。其中3號芯片功率最大，發(fā)熱最明顯。圖3(b)中顯示了僅受發(fā)動機發(fā)熱的電路板組件溫度分布圖，由于電路板前端靠近發(fā)動機，因此電路板前端區(qū)域溫度和熱量分布較高，后端區(qū)域溫度和熱量相對較低，溫度變化是由前往后逐漸降低的。圖3(c)中顯示了受芯片和發(fā)動機共同發(fā)熱下的電路板組件溫度分布圖，圖中溫度整體分布于圖3(b)趨勢相同，但是芯片處受熱有明顯不同。

通過對車載環(huán)境下PCB模型的溫度場分析，在實際應用中可以根據(jù)溫度的分布情況，選用合適的耐溫材料，在發(fā)動機等熱源與電路板接觸的路徑上增加隔熱保護層，采取必要的散熱設(shè)計[11]。大部分車載電子產(chǎn)品都對溫度很敏感，超過極限溫度時電子元件的性能將會收到嚴重影響，溫度大大超過工作溫度范圍，元件就可能會損壞[12]。溫度過高也會引起焊料的脫焊，大大降低焊點的可靠性。

圖3 三種情況下PCB溫度分布云圖

2.2 PCB變形的有限元分析

在溫度場下的熱分析后進行PCB的結(jié)構(gòu)分析，首先進行單元類型轉(zhuǎn)換，并施加位移條件。在實際發(fā)動機模塊電路板中，在上端兩個頂角和下邊中間部位有四個螺栓孔，用來將電路板固定。因此，模擬為了跟實際情況相符，對上述螺栓孔部位施加合理約束。最后把熱分析獲得的溫度數(shù)據(jù)作為熱載荷導入，進行靜力場求解，得到各方向位移分布云圖如圖4所示[13]。

圖4 PCB位移分布云圖

圖4顯示了三種受熱情況下電路板的整體位移情況。從圖4(a)中可以看出，在僅受芯片發(fā)熱條件下，功率較大，發(fā)熱較高的3號和4號芯片出現(xiàn)明顯的離面位移。說明該處芯片存在較大的熱變形。從圖4(b)和(c)中得知，僅受發(fā)動機發(fā)熱與受芯片共同發(fā)熱作用下的電路板熱變形整體趨勢相同，由此可知，在實際車載環(huán)境中，發(fā)動機發(fā)熱是影響電路板熱變形的主要因素。此外，觀察b和c圖中的芯片情況，發(fā)現(xiàn)芯片整體都呈現(xiàn)出中間凹陷，四周翹起的狀態(tài)。由于該電路板芯片都屬于PQFP封裝形式，芯片引腳都在芯片四周邊緣[14]。因此，出現(xiàn)上圖熱變形，很容易造成芯片四周處引腳焊點過度拉伸甚至斷裂，從而導致電路板芯片失效，使得電路板可靠性降低[15]。

表2 三種情況下Z坐標軸方向的最大最小位移

圖5 PCB Z軸位移分布云圖

從靜力場分析結(jié)果中可以得到面內(nèi)位移X、Y和離面位移Z的變形情況。在電路板組件中最容易引起失效部位是芯片和與芯片連接的焊點。對于芯片和焊點來說，里面位移Z方向的變形情況對其影響最大。如表2所示，下面是電路板在三種受熱情況下Z軸方向的最大最小位移。從表2中可以看出，三種情況下Z軸正向與負向最大值最小值變化差距并不大，但從圖5中觀察可得在芯片與發(fā)動機共同發(fā)熱下，電路板組件中的3號，4號和6號芯片都處于負向最大值區(qū)域，呈現(xiàn)明顯的凹陷狀態(tài)。因此，可以推斷該電路板組件中，由于受熱影響，3號，4號和6號芯片的引腳處焊點很容易失效，導致電路板整體可靠性降低。所以，上述芯片位置是該電路板組件的最危險處，在后續(xù)的設(shè)計與優(yōu)化中需要重點改進。

3 結(jié)論

通過對車載電路板在三種不同受熱情況下的的溫度熱應力仿真分析，觀察電路板組件的溫度分布情況，可以確定不同發(fā)熱條件對電路板熱變形產(chǎn)生的影響。通過觀察電路板組件的熱變形位移情況，找到了工作環(huán)境下電路板組件最危險處。仿真結(jié)果表明：（1）在實際車載工作環(huán)境中，發(fā)動機發(fā)熱是車載電路板產(chǎn)生熱變形的主要因素；（2）電路板組件中最大離面位移出現(xiàn)在3號，4號和6號芯片處，即該電路板中間區(qū)域位置，上述位置是電路板組件該狀態(tài)下的最危險區(qū)域。對車載電路板的熱應力分析，為后續(xù)車載電路板的設(shè)計優(yōu)化和散熱方式的選擇及安裝位置提供了參考。采用有限元模擬的方法能模擬出真實情況下電路板處于的環(huán)境溫度條件，能快速找到問題，對后續(xù)提高印制電路板的熱可靠性有很大意義。

[1] 郭強. 振動沖擊條件下SMT焊點疲勞壽命與可靠性的理論研究 [D]. 上海：上海交通大學，2005.

[2] 姜青龍，胡剛. 電印刷電路板溫度-應力耦合場有限元 [J]. 計算機科學與技術(shù)，2009, 27(2)：37-42.

[3] 于巖，王守緒，等. 基于COSMOSWORKS有限元分析的HDI板熱應力仿真 [J]. 印制電路信息，2012, 494（07）：494-500.

[4] Rodgers, B., Punch, J. and Jarvis, J. Finite Element Modelling of a BGA package subjected to power cycling[C]. International Society Conference on Thermal Phenomena. Piscatway, IEEE, 2002. 993-1000.

[5] Kenny C.Otiaba, R.S. Bhatti and N.N.Ekere. Finite Element Analysis of the Effect of Silver Content for Sn-Ag-Cu Alloy Compositions on Thermal Cycling reliability of Solder Die Attach[J]. Engineering Failure Analysis, 2012, 28(11):192-207.

[6] 馬巖. 印制電路板詳細模型的熱仿真分析[J].機械設(shè)計與制造工程，2016, 45(1): 52-55.

[7] Ubachs, R. and Schreurs, P. Elasto-viscoplastic nonlocal damage modeling of thermal fatigue in anisotropic lead-free solder[J]. Met Mater, 2007, 39(7): 685-701.

[8] 馬鵬，張學玲. 電路板的有限元熱分析及熱變形分析方法[J].科教導刊（電子版），2016，30(3):167-169.

[9] 張劍. 基于ANSYS的液氮生物容器熱應力分析[J]. 醫(yī)療裝備，2014, 27(2):14-16.

[10]張才洪. ANSYS14.0理論解析與工程應用實例[M]. 北京：機械工業(yè)出版社，2013.

[11]黃志亮，周云山. 預變形印刷電路板熱-力耦合結(jié)構(gòu)優(yōu)化設(shè)計[J]. 計算機輔助工程，2015, 24(4)：77-81.

[12]柴樹峰，張學玲. IRF4115S型MOSFET電路板有限元熱分析[J]. 軍事交通學院學報，2015，17(7)：91-95.

[13]Suhir, E. Thermal Stress Failures in Electronics and Photonics: Physics, Modeling, Prevention[J]. Journal of Thermal Stresses, 2013，36(6)：536-563.

[14]王興久. 多芯片雙面PCB的熱應力分析[J]. 南京理工大學學報（自然科學版），2010, 34（2）：170-175.

[15]陳苑明，何為，王守緒. 印制電路基板的溫變熱性能研究[J]. 印制電路信息，2012, 5（1）：489-493.

Thermal Stress Analysis of the PCB Assembly in Vehicle under Different Thermal Filed Based on Finite Element Method

ZHOU Jia-cheng, LIU Fang

(School of Mechanical Engineering and Automation, Wuhan Textile University, Wuhan Hubei 430073, China)

In order to study the influence of thermal produced by engine and chips to the PCB assembly in the vehicle under the driving car, the PCB in vehicle was taken as a study object by ANSYS. The the temperature field analysis and thermal stress analysis was carried out under the –thermal produced by chips, by engine and the both– three different thermal loadings. Comparing the distribution of temperature and thermal deformation under three different thermal loadings, the details of the PCB assembly can be obtained. The results showed that the thermal produced by engine played the major role in the thermal deformation of the PCB assembly under the working condition of the car. The maximum displacement of out-of-plane appeared the areas located on the PCB center. These areas of the chips were the most dangerous on the PCB assembly.

PCB Assembly in vehicle; thermal loading; chips; thermal stress

TN306

A

2095－414X(2017)06－0076－05

通訊作者：劉芳（1976-），女，副教授，博士后，研究方向：振動沖擊分析與控制.

湖北省教育廳科學研究計劃青年人才項目（No. Q20141608)；國家自然科學基金青年項目（No.11102141）.