半導(dǎo)體制冷箱的宏觀設(shè)計與結(jié)構(gòu)優(yōu)化

于俊達，彭寧濤，張 華，閆獻國，張光新，李 根

（1.太原科技大學(xué) 機械工程學(xué)院，山西 太原 030024；2.中機國際工程設(shè)計研究院有限責任公司，湖南 長沙 410007；3.山西省正誠礦山安全技術(shù)研究所，山西 太原 030012）

0 引言

半導(dǎo)體制冷又叫做溫差制冷或熱電制冷，是一種近年發(fā)展起來的新興制冷技術(shù)。它具有結(jié)構(gòu)簡單、維修方便、無磨損、無噪聲、壽命長、無需制冷液、絕對環(huán)保等優(yōu)點［1］，可應(yīng)用于電子、醫(yī)藥、科研實驗室等場合，有著廣闊的發(fā)展前景。

與傳統(tǒng)的機械壓縮式制冷相比，半導(dǎo)體制冷效率低，且其制冷效果受環(huán)境溫度的影響很大，因此本文針對某自制的半導(dǎo)體制冷箱（容積為10L，保溫層材料為聚苯乙烯泡沫塑料板，制冷片型號為TEC1－12706，數(shù)量為6片，且均勻布置在制冷箱的頂面，環(huán)境溫度為30℃），為了盡可能地增大其制冷效率，對其制冷箱電路、制冷箱內(nèi)膽形狀及保溫層的厚度進行了優(yōu)化。

1 制冷箱電路的設(shè)計

當溫差不變時，制冷量Qc與電流I 成二次函數(shù)，且函數(shù)曲線為開口向下的拋物線。制冷量Qc剛開始隨電流的增大而增大，當電流增大到某一值時制冷量達到最大值，然后隨著電流的增大制冷量反而下降，并逐漸降為零。因此可以用求極大值的方法求出制冷量最大時對應(yīng)的最佳電流Iopt。

制冷箱電路的設(shè)計應(yīng)遵守的條件及需實現(xiàn)的功能為：①盡可能使通過制冷片的最大電流為制冷量最大時的最佳電流；②由于半導(dǎo)體制冷片本身具有一定的電阻，各制冷片并聯(lián)時通過制冷片的電流更大，因此制冷片在電路中采用并聯(lián)連接；③整個電路設(shè)置一個電源開關(guān)，制冷片和風扇分別單獨設(shè)置一個電源開關(guān)；④制冷箱能同時實現(xiàn)制冷／加熱功能，設(shè)置切換開關(guān)實現(xiàn)制冷箱制冷／加熱間的切換，并設(shè)置指示燈，使制冷時為藍色LED 燈亮，加熱時為紅色LED 燈亮；⑤制冷箱工作時功率可以調(diào)節(jié)，并根據(jù)功率的大小設(shè)置3個檔位；⑥為降低制冷箱工作時箱內(nèi)溫度的不均勻性，制冷箱功率的大小通過調(diào)節(jié)流過所有制冷片總電流的大小來控制。根據(jù)上述6條規(guī)則設(shè)計出的制冷箱電路圖見圖1。

圖1 制冷箱電路布置圖

（1）電源開關(guān)：K0為控制整個電路通斷的電源開關(guān)，K1、K3分別為控制風扇和制冷片開啟／停止的開關(guān)。

（2）制冷／加熱切換：通過電流換向開關(guān)K2控制半導(dǎo)體制冷片中電流的流向來切換制冷箱的制冷／加熱。當換向開關(guān)K2與觸點A、B2連接時，若開關(guān)K0、K3、K4閉合，電流的流動路線為：電源正級—觸點A—經(jīng)并聯(lián)分流后分別流向6片半導(dǎo)體制冷片和發(fā)光二級管LED1—觸點B2—電源負級，此時流過半導(dǎo)體制冷片的電流為正向電流，且LED1燈亮（藍色），制冷箱處于制冷狀態(tài)；同理，當換向開關(guān)K2與觸點B1、C連接時，流過半導(dǎo)體制冷片的電流為反向，且LED2燈亮（紅色），制冷箱為加熱態(tài)。

（3）工作功率控制：通過旋轉(zhuǎn)開關(guān)K4控制接入電路中電阻的個數(shù)，達到改變流過每片制冷片電流的大小來控制制冷箱功率的目的。當K4與觸點I連接時，R2和R3兩個電阻都接入電路中，此時電路中的電阻最大，從而流過每片制冷片的電流最小，制冷箱工作功率最??；當K4與觸點J連接時，只有R3一個電阻接入電路中，此時電路中的電阻變小，從而流過每片制冷片的電流變大，功率增大；當K4與觸點K 連接時，R2和R3都沒有接入電路中，此時電路中的電阻最小，從而流過每片制冷片的電流最大（即最佳電流），功率最大。

2 制冷箱結(jié)構(gòu)尺寸的有限元優(yōu)化

有限元優(yōu)化時的基本假設(shè)為：①采用水冷散熱器散熱，制冷箱內(nèi)膽容積為10L，6片制冷片均勻分布在箱體頂面，保溫層的材料為聚苯乙烯泡沫塑料板；②環(huán)境溫度為30℃，制冷片熱端采用水冷散熱且散熱良好時冷端溫度能降到－15 ℃；③制冷箱門封保溫良好，在用ANSYS模擬時不考慮門封的冷量損失；④制冷過程的熱傳遞主要靠箱體內(nèi)膽的熱傳導(dǎo)和內(nèi)膽中空氣的熱對流，忽略熱輻射的影響［2］；⑤邊界條件采用第一類邊界條件［3］，取單元長0.05mm 均勻劃分網(wǎng)格［4］。

2.1 制冷箱內(nèi)膽形狀的優(yōu)化

設(shè)制冷箱內(nèi)膽的長、寬、高的初始值分別為250 mm、200mm、200mm。

將制冷箱內(nèi)膽的長和高作為輸入?yún)?shù)（即設(shè)計變量），箱內(nèi)溫度的最大值作為輸出參數(shù)（即狀態(tài)變量）。由于容積為定值10L，因此只要長和高確定了，寬也會確定下來［5］。

優(yōu)化時取1 000個樣本點進行比較，把最小化箱內(nèi)溫度最大值作為優(yōu)化目標。經(jīng)ANSYS 仿真優(yōu)化后，得到了制冷箱內(nèi)膽的長和高對箱內(nèi)溫度最大值的權(quán)衡圖，如圖2所示（圖2中越左邊的點表示越好，越右邊的點表示越差）［6］。

由圖2可知：箱體內(nèi)膽的高為125mm～135mm時比較好，其中130.48mm 最好；長對制冷箱制冷的影響不大，220 mm～300 mm 之間都比較好，其中244.48mm 最好。

為了方便取放東西，箱體內(nèi)膽的高度應(yīng)不小于150mm，又因為內(nèi)膽的長對制冷效果的影響很小幾乎可以忽略，考慮到箱體的外觀和制造等因素，經(jīng)ANSYS優(yōu)化后，制冷箱內(nèi)膽的長、寬、高分別取320mm、210mm 和150mm。

2.2 制冷箱保溫層厚度的優(yōu)化

令制冷箱保溫層的厚度為設(shè)計變量，令箱內(nèi)目標點（即中心點）C（160，105，75）的溫度為狀態(tài)變量，保溫層厚度的設(shè)計范圍為10mm～30mm，間斷性變化（遞增量為1）。優(yōu)化時取20 個點作為優(yōu)化比較的樣本點，經(jīng)ANSYS仿真優(yōu)化后，將每個樣本C 點的溫度值導(dǎo)入Origin軟件中可以得到保溫層厚度與C 點溫度的響應(yīng)曲線圖，如圖3所示。

圖2 內(nèi)膽箱尺寸與溫度最大值權(quán)衡圖

圖3 保溫層厚度與C 點溫度響應(yīng)圖

由圖3可知，當保溫層厚度小于20 mm 時，厚度的變化對制冷箱內(nèi)C 點溫度影響較大，當厚度大于25 mm 時，厚度對C 點溫度的影響相對較小。因此在后邊的模擬和實驗中取制冷箱保溫層的厚度為25mm。

3 制冷箱制冷過程有限元模擬

經(jīng)ANSYS優(yōu)化分析后，確定了制冷箱的結(jié)構(gòu)尺寸。此時，可利用專業(yè)的CFD 軟件Fluent進行制冷箱制冷過程的精確模擬。用Fluent模擬時將time選項設(shè)置為非定常模型unsteady，采用SIMPLE 方式離散壓力－速度，壓力修正松弛因子設(shè)定為1.0，以達到加快收斂速度的目的。制冷模擬時間取30min，時間步長取0.1s，并選擇Autosave，設(shè)置為每1 200步長保存一次（即每2min自動保存一次）。后處理時，取C 點作為目標點。

由Fluent模擬及實驗測量得到C 點溫度的時間歷程對比圖，如圖4所示。

由Fluent模擬也可以得到制冷溫度穩(wěn)定后的制冷箱溫度分布云圖，如圖5所示。

圖4 箱體C 點溫度模擬與實驗的對比

圖5 制冷箱Fluent分析溫度云圖

由圖4和圖5可知，經(jīng)優(yōu)化后，制冷箱內(nèi)的大部分空間在20min之內(nèi)降到2 ℃～4 ℃左右，當C 點溫度穩(wěn)定后，F(xiàn)luent模擬的數(shù)值為1.2 ℃。

4 結(jié)論

本文用ANSYS有限元軟件對半導(dǎo)體制冷箱的結(jié)構(gòu)進行優(yōu)化分析，其分析結(jié)果與實驗結(jié)果接近，應(yīng)用有限元方法進行制冷箱的研制，大大縮短了研發(fā)周期，節(jié)約了開發(fā)費用，同時其分析結(jié)果對于其他類型制冷箱的研發(fā)設(shè)計也具有一定的參考價值。

［1］ 簡棄非，梁榮光，張勇，等.半導(dǎo)體制冷效率及空間冷量傳遞特性試驗研究［J］.華南理工大學(xué)學(xué)報（自然科學(xué)版），2001，29（5）：72－75.

［2］ 代彥軍，戴維涵，王如竹.半導(dǎo)體冰箱冷熱端散熱條件實驗研究［J］.工程熱物理學(xué)報，2005，26（s1）：221－224.

［3］ 楊世銘，陶文銓.傳熱學(xué)［M］.北京：高等教育出版社，1998.

［4］ 朱紅均，林元華，謝龍漢.FLUENT 流體分析及仿真實用教程［M］.北京：人民郵電出版社，2010.

［5］ 徐昌貴，彭寧濤，閆獻國，等.半導(dǎo)體制冷箱的三維有限元優(yōu)化與實驗研究［J］.低溫與超導(dǎo)，2012，40（11）：70－75.

［6］ 胡國良，任繼文，龍銘.ANSYS 13.0 有限元分析實用基礎(chǔ)教程［M］.北京：國防工業(yè)出版社，2012.