基于Ansys Icepak 的負載器熱設計與熱分析

蔡惠華 賈豐鍇 張莉莉 孫德沖 孟翔宇

(北京航天計量測試技術研究所，北京 100076)

1 引 言

隨著人類對太空探索需求的增加，大運載火箭技術已經成為火箭發展的重要方向之一。由于大運載火箭的箭上設備越來越多，而火箭在飛行過程中，所有設備全部需要靠自身攜帶的電源進行供電。為了驗證箭上電源的帶載能力，設計供配電負載器用來模擬箭上負載，與其他配電設備配合使用，驗證配電電源設備等的帶載能力。

根據需要，模擬的負載器額定功率可以達到幾千瓦[1]，在工作過程中，會產生大量的熱量，這些熱量大部分使設備自身溫度升高，因此，負載器設計的重點不僅僅是實現需要的額定功率，而且更要考慮的是工作過程中的散熱問題。判斷負載器的散熱情況主要通過負載器的溫升指標來判斷。工作過程中溫度的測量，一般只能通過散熱孔，將相應的溫度測試端伸到發熱電阻上，對于一個額定功率為幾千瓦的負載器，機內往往有上百個大功率電阻，逐個測試，并找到溫升較高的電阻進行測試比較困難。本文基于負載器設計的理論分析，采用Ansys Icepak對負載器的散熱情況進行模擬仿真，找出負載器機箱內溫升最高的大功率電阻，指導實際測試。另外通過該軟件仿真，對供配電負載器進行了優化設計。

2 熱設計基礎理論及負載器工作原理

2.1 熱設計基礎理論[2]

熱力學第一定律認為，能量不會消失和創造，它只能在形式上進行轉換或者在物質間進行傳遞。其中熱量的傳遞是從高溫向低溫區進行傳遞，熱量傳遞的基本表達式如式(1)所示

Q=KAΔt

(1)

式中：Q——熱流量，W；K——換熱系數，W/(m2·℃)；A——換熱面積，m2；Δt——高低溫物質之間的溫差，℃。

熱傳遞的方式分為傳導、對流、輻射三種方式。其中傳導換熱、對流換熱都需要有導熱介質，而熱輻射無需任何介質，直接在空中進行能量傳遞。

2.1.1 熱傳導是由兩個物體間或同一個物體由于溫度梯度引起的介質間能量交換。熱傳導依據傅利葉定律表征。熱傳導的基本表達式如式(2)所示

Q=hcA(tw-tf)

(3)

式中：hc——對流換熱系數，W/(m2·°C)；A——有效的對流換熱面積，m2；tw——物體表面的溫度，℃；tf——周圍介質的溫度，℃。

2.1.3 輻射換熱是指物體間或者一個物體內熱量的交換方式是通過電磁波的形式向外傳遞能量的過程，也稱為熱輻射。熱輻射計算如式(4)所示

式中：δ0——為玻耳茲曼常數；A——物體輻射熱換的表面積；εxt——系統的發射率；F12——角系數；T1，T2——絕對溫度，K。

2.2 負載器工作原理

負載器的基本原理可以簡化成如圖1所示，根據大功率電阻器的功耗，將電阻器分解為如圖2所示。

圖1 負載構成示意圖 圖2 負載設計原理框圖 Fig.1 Load structuresketch map Fig.2 Load design principleframe map

將4個負載電阻固定在一片散熱器上，假設每個電阻功率是28W，則4個電阻器總功耗P=112W，考慮環境溫度為Ta=35℃，控制負載電阻的最高工作溫度為T=80℃，則需要散熱器的熱阻滿足下式

R=(T-Ta)/P=0.4℃/W

(5)

根據某型號散熱器手冊，在自然對流的情況下，熱阻為0.4℃/W時，需要此種散熱器的長度為220mm，才能滿足設計要求。另外，在機箱側部安裝軸流風機，進行強迫風冷，進一步增強散熱效果。

2.3 負載器的熱分析

2.3.1 發熱熱源分布[3]

如何進行熱源計算是熱分析的基礎，負載器的熱源主要是供電的導線及大功率電阻組成。對于負載器的熱源計算，如果考慮所有的單元細節，會導致計算比較復雜，由于大功率電阻的溫度變化遠遠大于其它熱源，因此在進行熱分析時，以大功率電阻達到熱穩態后的發熱作為主要分析對象。

2.3.2 熱量傳遞分布[4]

當系統達到熱穩態后，負載器的熱量傳遞可以考慮以下幾個方面：(1)負載器內部各器件與散熱器之間；(2)殼體內空氣與外殼之間的導熱；(3)外殼的強制對流散熱和輻射散熱。由于負載器熱源電阻通過散熱器和軸流風機進行強制風冷散熱，因此在散熱分析中主要考慮對流換熱與熱傳導換熱，不需要考慮輻射換熱。

3 負載器熱仿真

Ansys Icepak 軟件仿真的過程和大部分軟件仿真過程類似，主要分為以下四個步驟：首先是建模，即創建軟件分析所需的基本實體模型；其次施加載荷，即設置相關參數；再次是求解，即進行軟件的分析運算；最后查看和分析軟件的仿真結果。

3.1 創建實體模型

創建實體模型是用于創建真實的熱傳遞路徑模型，對于供配電負載器，實體模型包括機架結構、散熱器、導熱硅脂、大功率電阻、螺母、螺絲、把手、接插件等。ANSYS Icepak軟件提供軟件自建模方式和CAD模型導入兩種方式。自建模方式中，軟件提供了常用的熱源、散熱器、開口、邊界、軸流風機等多種方式，滿足大部分常用的電子設備需求。對于供配電負載器產品，已建立CAD模型，因此，本文直接采用CAD模型導入的方法。

3.1.1 CAD模型導入

(1) 模型修復：實體模型中，有一些特征或者器件對散熱不影響或者影響不大的需要進行刪除，這個過程稱為模型修復，模型修復有利于軟件仿真中網格的合理劃分，所以CAD模型導入之前，一般需要進行模型修復。模型修復有兩種使用方法，一種方法是可以使用DM軟件中直接進行修復，另一種方法是先通過ANSYS SCDM軟件提供的平臺進行修復，本文采用后一種方法。

(2) 建立工作圖在Ansys Workbench平臺下建立仿真的工作圖，如圖3所示。雙擊圖3 A中的Geomrtry，打開DM軟件。DM中Tools→Electronics命令，是ANSYS公司專門給Ansys Icepak提供的數據導入轉換接口。

圖3 Ansys Workbench平臺下仿真的工作圖Fig.3 Ansys Workbench emluation work map

(3) 模型導入及轉換：將修復后的負載器CAD模型通過該指令的相關選項，轉換成Ansys Icepak認可的模型，導入Ansys Icepak中，雙擊圖3 B中的Setup，就可行性模型相關設置了。

3.1.2 CAD模型簡化與設置

(1) 模型簡化：為降低模型對計算機內存的要求，需要減少網格數量，本文對負載器仿真模型進行了兩方面簡化，一是由于負載器的結構采用上下層的設計，上下層散熱器、電阻功率完全一致，只對模型上半部分進行仿真；另一方面由于負載器是采用軸流風機進行強迫風冷散熱，因此，忽略了機箱外殼與空氣的自然對流及輻射換熱，使用Icepak軟件自建模提供的Cabinet特征作為機箱的外殼，從而達到了減少網格的數量。

(2) 參數設置：負載器的參數設置主要包括散熱孔開孔率設置、軸流風機風量設置、熱源功率設置、導熱硅脂導熱系數設置等。機箱的散熱孔采用自建模型中的Grille代替，通過計算，計算出負載器的開孔率為0.201，在Grille中的Free area ratio中輸入0.201。負載器選用的軸流風機固定風量為0.023 33m3/s。，根據散熱器數據手冊，軸流風機的在平均風速3m/s情況下，負載器選用的散熱器熱阻長度為220mm，熱阻是0.15℃/W。負載器大功率電阻用相應功率大小的熱源表示，固定功率為28W。導熱硅脂采用自建模型中的plate模擬，導熱系數為0.913，厚度1mm。

3.2 網格劃分及運算

網格劃分關系到計算結果的準確性，網格劃分一般要經歷先生成粗糙網格，再進行評估，再細化，再檢查生成的網格是否滿足要求，再細化，再求解等過程。ANSYS Icpeak軟件提供了三種主要類型的網格：六面體結構化網格、六面體非結構化網格及六面體核心結構化網格。這三種主要的網格類型都可以按實際要求進行局部細化。但是網格的劃分不是越細越好，網格劃分越細，網格數量越多，節點數也越多，雖然計算精度越高，但對計算機的內存也要求越高，仿真時間越長。另外，在相同節點數量情況下，計算精度會受網格疏密承擔的影響。綜合考慮三種網格類型及計算機處理能力，本文選擇六面體核心結構化網格的網格類型進行網格劃分，負載器模型進行網格劃分完成后可直接進行運算。

4 結果分析

4.1 仿真結果

為了分析找到溫升最高的負載電阻，本文通過兩個步驟進行：第一步，對每片散熱器上的負載電阻進行分析、仿真，找出該散熱片上溫升最高的負載電阻；第二步，將每片散熱器上溫升最高的負載電阻再進一步比較。

每片散熱器上有四個負載電阻，這四個負載電阻的結構分布如圖4所示。將這四個負載電阻設置為變量監控點。溫度變化仿真結果如圖5所示。從圖5仿真結果可以看出，source104負載電阻器遠離軸流風機，且介于其他負載電阻器之間所以溫度遠高于其他負載電阻器。

圖4 散熱片上負載電阻分布示意圖Fig.4 Load resistance distribution in radiator sketch map

圖5 四個負載電阻溫度變化Fig.5 Change temperature of four load resistance

通過第一步仿真，明確找出每片散熱器上溫升最高的負載電阻器為source104，接下去就是比較各個散熱片上相應位置的負載電阻，如圖6所示。其溫度變化仿真結果如圖7、圖8所示，各個負載電阻溫度相差在2℃范圍內，在實際測試過程中，由于每片散熱器上溫升最高的負載電阻溫度測量值相差不大，另外，考慮到仿真存在誤差，如果對負載器溫升值有精確要求，仍需要對每片散熱器上溫升最高的負載電阻進行溫度測試。

圖6 各個散熱片上的相應位置的負載電阻Fig.6 All load resistance distribution in every radiator sketch map

圖7 各個散熱片上溫升最高負載電阻的溫度變化Fig.7 Change temperature of temperature rise highest load resistance in every radiator

圖8 各個散熱片上溫升最高負載電阻的溫度變化局部放大圖Fig.8 Zoom out part’s map of change temperature of temperature rise highest load resistance in every radiator

4.2 優化應用

在實際的優化設計中，由于產品的功率是固定的，同時由于產品應用環境的限制以及為了保證產品體積小、重量輕的優勢，供配電負載器的尺寸和結構也無法改變，因此負載器熱特性的優化空間是十分有限的。

目前現有供配電負載器散熱孔采用圓形設計，開孔率僅為0.201。針對散熱孔形狀和大小，本文進行了優化設計。

如果將原先圓形散熱孔改成正六邊形(圖9)，半徑由2mm增加到2.5mm，開孔率能達到0.437時，通過仿真表明，由于熱源功率固定，負載電阻器一直處于發熱狀態，所以負載電阻器溫度變化不大(圖10)，但負載器內部溫升能降低7℃左右，如圖11所示。

圖10 正六邊形開孔時負載電阻的溫度變化Fig.10 Change temperature of load resistance with hexagon holes

5 結束語

本文基于Ansys Icepak軟件對負載器進行熱仿真，簡化仿真模型，降低對計算機內存的需求；分析比較各個負載熱阻溫度變化情況，得到溫升最高的負載熱阻，為實際過程中測試提供仿真依據，避免盲目測試；另外，本文優化散熱孔設計，采用正六邊形開孔，能進一步提高開孔率，以達到更好散熱效果。

(a) 圓形開孔時負載器內部溫度圖(a) Inner temperature in load equipments with round holes

(b) 正六邊形開孔時負載器內部溫度(b) Inner temperature in load equipments with hexagon holes圖11 不同開孔形狀的負載器內部溫度分布圖Fig.11 Inner temperature in load equipments with different holes

[1] 姜學東等.大功率電子負載在通信電源系統中的應用[J].通信電源技術, 2000(2):8～11.

[2] 王永康.ANSYS Icepak 電子散熱基礎教程[M].北京：國防工業出版社, 2014.

[3] 黃世澤,郭其一等.基于Ansys的電磁系統溫度場仿真技術研究[J]. 電器與能效管理技術, 2015(2):32～37.

[4] 吳燕，劉幗巾.基于Ansys的小型汽車繼電器的溫度場仿真[J].電器與能效管理技術, 2015(2):27～31.