李　灝，錢新明，陳　威

(爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，北京理工大學(xué)，北京　100081)

?

LED太陽花散熱器正交試驗(yàn)?zāi)M優(yōu)化設(shè)計(jì)

李灝，錢新明，陳威

(爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，北京理工大學(xué)，北京100081)

摘要：利用仿真軟件ICEPAK對一款大功率LED太陽花散熱器進(jìn)行了熱仿真模擬。實(shí)驗(yàn)數(shù)據(jù)曲線和模擬結(jié)果曲線的相似度分析結(jié)果表明利用ICEPAK進(jìn)行太陽花散熱器熱仿真模擬的可行性。為了提高太陽花散熱器的散熱效率以及控制散熱器重量，使用正交試驗(yàn)法對太陽花散熱器進(jìn)行模擬優(yōu)化設(shè)計(jì)，從散熱器圓柱半徑、熱沉高度、翅片厚度、翅片數(shù)量、翅片長度等5個(gè)方面考察了太陽花散熱器的尺寸參數(shù)對散熱器重量與LED最高溫度的影響。得到了一個(gè)散熱效率最優(yōu)化結(jié)果，優(yōu)化后LED的最高溫度較優(yōu)化前下降了11℃，重量基本不變。

關(guān)鍵詞：散熱器；正交試驗(yàn)；優(yōu)化設(shè)計(jì)；模擬

引言

相對于其他照明方式，LED照明具有明顯的優(yōu)勢：體積小、功耗低、綠色環(huán)保、使用壽命長、安全、豎固耐用，其在照明領(lǐng)域得到了很大的應(yīng)用[1]。但是隨著尺寸的減小以及功率的大幅提高，導(dǎo)致LED 結(jié)溫居高不下，引起了光強(qiáng)降低、光譜偏移、色溫升高、熱應(yīng)力增高、元器件加速老化等一系列問題，大大降低了 LED 的使用壽命[2]。對于大功率LED芯片來說，一般超過80%的輸入功率被轉(zhuǎn)化為熱能，其余部分被轉(zhuǎn)化為光能[3]。而半導(dǎo)體元件溫度的上升會(huì)造成其可靠性和安全性的下降[4]。研究顯示，電子器件的故障率隨元件溫度的升高呈指數(shù)關(guān)系增加，器件溫度在70℃～80℃水平上每增加1℃，可靠性下降5%[5]。目前，LED熱管理廣泛采用的技術(shù)是利用外部散熱器散熱。LED主要散熱部件的散熱器在散熱能力與散熱成本上的合理配置，是散熱設(shè)計(jì)的首要目標(biāo)。散熱設(shè)計(jì)是保證LED照明燈具散熱能力的首要保障措施。常見的大功率LED散熱器分為直肋式和太陽花式。楊桂婷[6]對直肋式散熱器的翅片長度、寬度和數(shù)量對散熱器溫度的影響進(jìn)行了模擬仿真。梁融[7]在對LED 肋片散熱器結(jié)構(gòu)數(shù)值分析基礎(chǔ)上，提出了開縫對散熱器的影響。劉靜[8]利用ICEPAK軟件對大功率LED路燈進(jìn)行建模仿真，研究了肋間距、肋厚度、換熱面積等結(jié)構(gòu)的優(yōu)化。應(yīng)用ICEPAK對太陽花結(jié)構(gòu)進(jìn)行熱模擬研究較少。本文將正交試驗(yàn)設(shè)計(jì)方法用于太陽花散熱器的優(yōu)化設(shè)計(jì)，通過ICEPAK仿真模擬，考察了太陽花散熱器各項(xiàng)尺寸參數(shù)對給定功率的芯片達(dá)到穩(wěn)態(tài)時(shí)散熱器溫度的影響，獲得了優(yōu)化尺寸參數(shù)組合。提高了太陽花散熱器的散熱效率。

1太陽花散熱器實(shí)驗(yàn)測量及仿真建模

選用一款60W的LED太陽花散熱器作為研究對象。首先利用對太陽花散熱器的實(shí)物測量和模擬對比，驗(yàn)證模擬的可靠性。太陽花散熱器的實(shí)物與Solidworks模型如圖1所示。圖1(b)中由上至下有3部分，分別為LED基板、導(dǎo)熱硅脂和熱沉，熱沉內(nèi)部稱為圓柱，外部稱為翅片。熱沉的材料為鋁。芯片采用COB封裝大功率LED，其總功率為60W，散熱器模型的結(jié)構(gòu)和材料參數(shù)如表1所示。利用紅外熱像儀進(jìn)行溫度測量，實(shí)驗(yàn)時(shí)室溫為25℃。

表1　太陽花散熱器的結(jié)構(gòu)和材料參數(shù)

圖1　太陽花散熱器實(shí)物與仿真模型Fig.1　Physical and simulation model of sunflower radiator

COB芯片不斷產(chǎn)生熱量， 導(dǎo)熱硅脂和散熱器進(jìn)行熱傳導(dǎo)，太陽花翅片與周圍空氣發(fā)生自然對流換熱。因此，近似地把問題看作是三維、穩(wěn)態(tài)、常物性、有內(nèi)熱源的導(dǎo)熱和對流換熱的耦合問題。在計(jì)算過程中由溫差引起的輻射換熱忽略不計(jì)，由于溫差而引起的浮生力作用， ICEPAK在計(jì)算中引入了Boussinesq假設(shè)[9]：①流體中的粘性耗散項(xiàng)忽略不計(jì)；②除密度外其他物性皆為常數(shù)；③密度僅考慮動(dòng)量方程中與體積力有關(guān)的項(xiàng)，其余各項(xiàng)中的密度作常數(shù)處理。計(jì)算區(qū)域采用有限容積法在同位網(wǎng)格上進(jìn)行控制方程的離散，k-ε雙方程模型求解。壓力、動(dòng)量和能量方程中的對流項(xiàng)均采用二階迎風(fēng)格式?？刂品匠虨椋?/p>

(1)

(2)

(3)

(4)

(5)

在Solidworks中進(jìn)行散熱器建模，為了建模及劃分網(wǎng)格的可靠性，刪除了定位孔及對波浪形翅片進(jìn)行了簡化。將模型導(dǎo)入ICEPAK中進(jìn)行模擬仿真，根據(jù)所采用的LED規(guī)格設(shè)定LED的光電轉(zhuǎn)換效率η為20%，即散熱器功率為48W，采用面光源。電源為外置電源，因此在散熱建模仿真時(shí)電源發(fā)熱量不計(jì)入散熱器功率中。環(huán)境溫度為溫度測量時(shí)的25℃，基座與實(shí)驗(yàn)臺接觸的面為絕熱且氣流不能穿過的平面，另外5個(gè)面為opening。實(shí)驗(yàn)與模擬結(jié)果的對比如圖2所示。利用Matlab的corrcoef函數(shù)進(jìn)行相關(guān)性分析，相關(guān)系數(shù)在0到1 范圍內(nèi)，越接近1，說明兩條曲線的相似度越高，對穩(wěn)態(tài)下實(shí)驗(yàn)數(shù)據(jù)曲線和模擬結(jié)果曲線分析得到相似度結(jié)果為0.9746。這表明了利用ICEPAK進(jìn)行太陽花散熱器熱仿真模擬的可行性。

圖2　60W太陽花實(shí)驗(yàn)與模擬溫度對比Fig.2　Comparison of experimental result with simulation result for 60W sunflower radiator

2太陽花散熱器正交試驗(yàn)優(yōu)化設(shè)計(jì)

2.1正交試驗(yàn)法

正交試驗(yàn)設(shè)計(jì)法具有完成試驗(yàn)要求所需的實(shí)驗(yàn)次數(shù)少、數(shù)據(jù)點(diǎn)分布均勻、可用相應(yīng)的極差分析方法等優(yōu)點(diǎn)[10]。

將正交試驗(yàn)優(yōu)化法引用到太陽花散熱器設(shè)計(jì)中，可優(yōu)化LED太陽花散熱器結(jié)構(gòu)，提高散熱能力，控制散熱器的質(zhì)量。其主要方法是對散熱器進(jìn)行結(jié)構(gòu)優(yōu)化設(shè)計(jì)，確定各相關(guān)尺寸因素的最佳參數(shù)組合，綜合考慮散熱效果和散熱器重量。

2.2　太陽花散熱器正交試驗(yàn)優(yōu)化設(shè)計(jì)

選定太陽花散熱器5個(gè)關(guān)鍵影響參數(shù)，分別為A:圓柱半徑、B:熱沉高度、C:翅片長度、D:翅片厚度、E:翅片數(shù)度。每個(gè)影響參數(shù)有五個(gè)水平值。根據(jù)影響參數(shù)和水平值選用L25(56)的正交試驗(yàn)表來安排模擬仿真，正交試驗(yàn)表如表2所示。按照表2中的實(shí)驗(yàn)數(shù)據(jù)分別進(jìn)行建模仿真，得出參數(shù)在不同水平組合下的LED最高溫度(Tmax)及散熱器重量(G)。

表2　太陽花散熱器正交試驗(yàn)設(shè)計(jì)安排

整理后得出各個(gè)關(guān)鍵參數(shù)對應(yīng)的溫度和質(zhì)量極差及其變化曲線圖，如圖2~圖6。通過對圖2~圖6的分析得出，翅片長度和圓柱直徑對溫度的影響最高，而翅片厚度對溫度的影響最小，各參數(shù)影響順序依次是A>C>B>D>E；而圓柱半徑對重量的影響最大，翅片厚度對重量的影響最小，各參數(shù)影響順序依次是A>B>D>E>C。由極差分析結(jié)果可知，圓柱半徑的大小對散熱器最高溫度和散熱器重量來說都是影響最大的因素。

圖2　圓柱半徑與LED最高溫度和散熱器重量的關(guān)系Fig.2　Effect of cylindrical radius on the maximum temperature of LED and the weight of radiator

圖3　熱沉高度與LED最高溫度和散熱器重量的關(guān)系Fig.3　Effect of radiator height on the maximum temperature of LED and the weight of radiator

圖4　翅片長度與LED最高溫度和散熱器重量的關(guān)系Fig.4　Effect of fin thickness on the maximum temperature of LED and the weight of radiator

圖5　翅片厚度與LED最高溫度和散熱器重量的關(guān)系Fig.5　Effect of fin thickness on the maximum temperature of LED and the weight of radiator

圖6　翅片數(shù)量與LED最高溫度和散熱器重量的關(guān)系Fig.6　Effect of the number of fins on the maximum temperature of LED and the weight of radiator

對流換熱能力與對流換熱系數(shù)及換熱面積有關(guān)。隨著圓柱半徑的增加，散熱器溫度呈下降趨勢，但是下降趨勢放緩。圓柱半徑過小，導(dǎo)致翅片根部間距小，影響空氣流動(dòng)，隨著半徑增加，空氣流動(dòng)速度在增大，但是繼續(xù)增加半徑會(huì)造成熱阻增大。增加翅片長度是有效增加對流換熱面積的方式，過長的翅片會(huì)造成翅片根部的空氣流速減小，同時(shí)造成熱阻增大。增加熱沉高度也是增加對流換熱面積，但高度增加到一定程度散熱器會(huì)出現(xiàn)明顯的溫度差，下端的散熱器溫度較低。增加翅片數(shù)量和增加翅片厚度對溫度影響較小，一味增加翅片數(shù)量或者厚度反而會(huì)削弱翅片間的空氣流速，進(jìn)而減小傳熱系數(shù)。

圖7　優(yōu)化前仿真結(jié)果Fig.7　Simulation results before the optimization

圖8　優(yōu)化后仿真結(jié)果Fig.8　Simulation results after the optimization

根據(jù)優(yōu)化試驗(yàn)結(jié)果綜合考慮，選定散熱器結(jié)構(gòu)參數(shù)組合為A3C4B4D2E2。圖7和圖8是優(yōu)化前后，散熱器模擬溫度的對比圖。優(yōu)化前的模型尺寸為：圓柱半徑40mm、翅片長度40mm、熱沉高度65mm、翅片厚度2mm、翅片數(shù)量36；優(yōu)化后的模型尺寸為圓柱半徑40mm、翅片長度45mm、熱沉高度75mm、翅片厚度1.6mm、翅片數(shù)量32。通過仿真模擬得出其優(yōu)化后散熱器最高溫度為64℃，散熱性能大幅提升。優(yōu)化后的LED最高溫度下降11℃。優(yōu)化后的散熱器重量為1390g，與優(yōu)化前的1397g相比略微下降。

3結(jié)論

本文通過對實(shí)驗(yàn)數(shù)據(jù)曲線和模擬結(jié)果曲線的相似度分析，表明了ICEPAK應(yīng)用于太陽花散熱器熱模擬的可行性。應(yīng)用正交試驗(yàn)結(jié)合ICEPAK對一款常用的60W太陽花散熱器尺寸(圓柱半徑40mm、翅片長度40mm、熱沉高度65mm、翅片厚度2mm、翅片數(shù)量36)進(jìn)行了模擬優(yōu)化設(shè)計(jì)，在維持散熱器重量基本不變的基礎(chǔ)上，優(yōu)化后的太陽花散熱器(圓

柱半徑40mm、翅片長度45mm、熱沉高度75mm、翅片厚度1.6mm、翅片數(shù)量32)能將LED芯片穩(wěn)態(tài)溫度較優(yōu)化前降低11℃，顯著提高了太陽花散熱器的散熱效率。

參考文獻(xiàn)

[1] ZUKAUSKAS A, SHUR M S, GASka R. Introduction to solid-state lighting[M]. New York: John Wiley & Sons Inc, 2002：21-29.

[2] 付賢政，胡良兵.LED 燈的散熱問題研究[J].照明工程學(xué)報(bào)，2011，22(3):73-77．

[3] CHENG Ting,LUO Xiaobing. Thermal analysis and optimization of multiple LED packaging based on a general analytical solution[J]. International Joournal of Thermal Sciences,2010,49(1):196-201.

[4] HUANG B,HSU P.,WU M.. Study of system dynamics model and control of a high-power LED lighting luminaire[J]. Energy,2007,32(11):2187-2198.

[5] 過增元.國際傳熱學(xué)研究前沿—微尺度傳熱[J].力學(xué)進(jìn)展，2000，30(1)：1-6.

[6] 楊桂婷，劉一兵.大功率LED燈具散熱裝置的設(shè)計(jì)[J].照明工程學(xué)報(bào)，2014，25(6)：91-94.

[7] 梁融，聶宇宏，聶德云，等.LED散熱器散熱特性分析及結(jié)構(gòu)優(yōu)化[J].照明工程學(xué)報(bào)，2014，25(1)：94-97.

[8] 劉靜，劉生春.大功率LED照明器的熱設(shè)計(jì)[J].光學(xué)與光電技術(shù).2008,6(5)：91-93.

[9] 陶文銓. 數(shù)值傳熱學(xué)[M].第2版．西安：西安交通大學(xué)出版社，2003.

[10] 張琦，陳旭.LED路燈熱分析及散熱結(jié)構(gòu)設(shè)計(jì)[J].電子與封裝，2009，9(5):44-48.

本刊聲明

為適應(yīng)我國信息化建設(shè)，擴(kuò)大本刊及作者知識信息交流渠道，本刊已被俄羅斯《文摘雜志》(AJ, VINITI)、美國《烏利希期刊指南》(Ulrich’s Periodicals Directory，UPD)、波蘭《哥白尼索引》(IC，Index of Copernicus)，《萬方數(shù)據(jù)庫》、《中國學(xué)術(shù)期刊網(wǎng)絡(luò)出版總庫》及CNKI系列數(shù)據(jù)庫、《中文科技期刊數(shù)據(jù)庫》、《中國學(xué)術(shù)期刊文摘數(shù)據(jù)庫》、中國郵政中郵閱讀網(wǎng)、超星數(shù)字圖書館等多家機(jī)構(gòu)收錄，其作者著作權(quán)使用費(fèi)與本刊稿酬一次性給付。免費(fèi)提供作者文章引用統(tǒng)計(jì)分析資料。若作者不同意文章被收錄，請?jiān)趤砀鍟r(shí)向本刊聲明，本刊將做適當(dāng)處理。

Optimum Design of LED Sunflower Radiator Based

on Orthogmal Experiment

Li Hao，Qian Xinming，Chen Wei

(StateKeyLaboratoryofExplosionScienceandTechnology,BeijingInstituteofTechnology,Beijing100081,China)

Abstract：Thermal simulation of a large power LED sunflower radiator was carried out by using simulation software ICEPAK. Comparing the experimental data curves with simulation results, it is feasible to simulate the thermal behavior of sunflower radiator by using ICEPAK. In order to improve the heat dissipation efficiency of sunflower radiator and control the weight of the radiator, using the orthogonal test method to optimize the design of the radiator. The influence of 5 factors, such as cylindrical radius, heat sink height, fin thickness, fin number and fin length, is analyzed. An optimized result is obtained, which is the maximum temperature of LED decreased by 11 degree centigrade, with the weight unchanging.

Key words:radiator; orthogonal experiment; optimum design; simulation

通訊作者：錢新明，E-mail:qsemon@bit.edu.cn

中圖分類號：TM923

文獻(xiàn)標(biāo)識碼：A

DOI:10.3969j.issn.1004-440X.2016.01.022