劉迪輝 楊三思 高鶴萱 張銀海 郭照宇

通過Fluent流場分析計算出冷卻水道管壁各節(jié)點處的對流傳熱系數(shù)，計算過程中考慮了水道結(jié)構(gòu)所導致的流場分布差異對傳熱能力的影響;利用LS-DYNA完成板料成形及保壓淬火過程的仿真模擬，獲取板料上的溫度數(shù)據(jù);以A柱邊梁為研究對象將仿真分析結(jié)果同熱沖壓生產(chǎn)線上的實際冷卻效果進行對比. 結(jié)果表明，考慮流場影響的模溫分析結(jié)果中高溫區(qū)域的分布和實測一致，12個溫度監(jiān)測點處最大溫度誤差為3.6 ℃，平均溫度誤差為0.96 ℃，而采用平均傳熱系數(shù)的仿真中溫度分布及最大溫差點發(fā)生改變，監(jiān)測點最大溫度誤差為9.1 ℃，平均溫度誤差為4.96 ℃. 仿真與試驗對比表明：流速的改變對傳熱系數(shù)的影響顯著，基于流場分析的模溫分析方法考慮了熱沖壓模具冷卻水道中流速差異性分布的影響，能夠用于復雜零件及多腔模具的模溫分析.

關(guān)鍵詞：熱沖壓;流場;對流傳熱系數(shù);冷卻系統(tǒng);溫度場分布

Abstract：A mold temperature analysis method for the hot stamping dies of complex formed parts considering the influence of different flow field distribution is proposed. First， the convective heat transfer coefficients at each node of the cooling water pipe wall were calculated by Fluent flow field analysis， and the influence of the differences in the flow field distribution caused by the water channel structure on the heat transfer capacity was considered in this process. Then， the processes of sheet forming and holding and quenching were simulated by LS-DYNA to obtain the temperature data of the blank. Finally， the A-pillar side beam was used as the research object to compare the simulation analysis results with the actual cooling effect on the hot stamping production line. The results show that the distribution of high temperature areas in the analysis result of the mold temperature based on the flow field analysis is consistent with the actual measurement. The maximum error at the 12 temperature monitoring points is 3.6 ℃，and the average error is 0.96 ℃. However， the temperature distribution and the maximum temperature difference point change in the simulation using the average heat transfer coefficient. The maximum error of the monitoring point is 9.1 ℃，and the average error is 4.96 ℃.The comparison between the simulation and experiment shows that the change of the flow rate has a significant effect on the heat transfer coefficient， and the mold temperature analysis method based on flow field analysis takes into account the effect of the differential distribution of the flow rate in the cooling channel of the hot stamping mold，which can be used for the analysis of mold temperature of complex parts and multi-cavity molds.

大量使用熱成形件是實現(xiàn)車身輕量化設(shè)計的重要途徑之一. 超高強度鋼板具有較高的比強度，在增加車身強度的同時能有效減輕質(zhì)量. 在熱成形過程中板料加熱到920～950 ℃時完全奧氏體化，具有良好的成形性，成形后在模具中保壓淬火完成向馬氏體組織的轉(zhuǎn)變，強度大大提高. 由于熱成形模具的水道結(jié)構(gòu)復雜，零件局部可能由于水道設(shè)計不合理導致淬火冷卻速率不夠，不能形成馬氏體，這會導致零件各個部位的材質(zhì)不同. 冷卻水道的設(shè)計是熱成形模具開發(fā)中的關(guān)鍵技術(shù)，模溫仿真分析方法可以用來評估冷卻水道的冷卻效果，優(yōu)化水道設(shè)計.

近年來，研究者對熱沖壓成形件起皺、開裂及微觀組織等成形質(zhì)量進行了廣泛的研究，研究表明，通過調(diào)整冷卻水流速、模具導熱參數(shù)，更改局部淬火速率等能夠完成對成形件微觀組織的控制，實現(xiàn)軟硬分區(qū)或差異性強化[1-3].? 文獻[4-6]通過仿真模擬的方法探究了水道結(jié)構(gòu)對保壓結(jié)束溫度的影響，仿真結(jié)果表明，冷卻水道的結(jié)構(gòu)設(shè)計會影響流場分布從而影響傳熱. 文獻[7-10]在研究沖壓成形參數(shù)對溫度分布、起皺、破裂及回彈等影響的過程中，對于冷卻水道和模具的換流傳熱均采用常數(shù)值定義，處理方式趨于理想化，而在實際生產(chǎn)中熱成形模具內(nèi)部冷卻水道結(jié)構(gòu)復雜，局部流速存在差異，仿真時采用平均傳熱系數(shù)不再合適.

仿真分析能否獲得和生產(chǎn)實際一致性較好的結(jié)果主要取決于仿真邊界條件的設(shè)置，模具同冷卻水之間的對流傳熱系數(shù)作為傳熱邊界條件之一，它的設(shè)置至關(guān)重要. 本文首先將傳熱系數(shù)仿真值和理論計算結(jié)果進行對比，證明采用基于Fluent流場分析求解的傳熱系數(shù)來定義冷卻水同模具的對流傳熱的可行性;其次將基于流場分析求解對流傳熱系數(shù)的模溫分析結(jié)果同基于傳熱系數(shù)平均值的分析結(jié)果及熱沖壓實際生產(chǎn)線上的冷卻效果進行對比驗證，對比結(jié)果表明，采用考慮流場影響的模溫分析方法能夠提高仿真精度，對于指導模具冷卻水道設(shè)計、縮短模具開發(fā)周期具有重大意義.

1? ?考慮實際流場的熱成形模溫分析方法

考慮實際流場影響的熱成形模溫分析方法關(guān)鍵在于通過流體仿真分析及數(shù)據(jù)處理獲取如圖1所示的A～D等節(jié)點處的對流傳熱系數(shù)，用于定義模溫分析過程中模具同冷卻水的對流傳熱.

分析流程如圖2所示. 首先，采用Fluent軟件對水道模型進行流場分析，在結(jié)果中輸出水管壁上各節(jié)點處的對流傳熱系數(shù);然后，通過二次開發(fā)軟件進行數(shù)據(jù)處理，將計算的對流傳熱系數(shù)文件轉(zhuǎn)換為LS-DYNA可識別的關(guān)鍵字數(shù)據(jù)格式文件，用于定義模具和冷卻水在管壁處的對流傳熱;最后，在LS-DYNA中完成沖壓成形和保壓淬火的仿真計算，獲得熱沖壓模具及板料的溫度場數(shù)據(jù).

上述分析方法發(fā)揮了Fluent軟件在流體分析方面的優(yōu)勢，同時利用LS-DYNA軟件在成形仿真和模溫分析方面的長處.

2? ?對流傳熱系數(shù)理論與仿真對比

2.1? ?仿真分析

實際沖壓過程中板料在模具中的保壓時間短，冷卻系統(tǒng)出入口溫差較低，一般為6～7 ℃，水的黏度、熱導率等受溫度影響的參數(shù)在低溫差條件下隨溫度變化并不明顯，文中主要考慮復雜水道結(jié)構(gòu)所導致的流速差異對水道表面對流傳熱系數(shù)的影響. 基于Fluent求解的對流傳熱系數(shù)可通過對比簡單水管模型表面?zhèn)鳠嵯禂?shù)的理論計算值進行準確性驗證，水管模型長300 mm，直徑8 mm，入口流速3 m/s，水溫10 ℃，出口壓力0 MPa（Fluent中壓力值均為表壓），10 ℃時水的相關(guān)參數(shù)見表1. 基于Fluent流速分析計算出的測試模型表面各處的對流傳熱系數(shù)分布云圖如圖3所示.

2.2? ?對流傳熱系數(shù)準確性驗證

參考表1中數(shù)值，根據(jù)公式（2）（3）（4）計算努賽爾數(shù)，代入公式（1）計算對流傳熱系數(shù). 對流傳熱系數(shù)理論計算結(jié)果為10 556.4 W/（m2·K），基于Fluent流場分析求得的水道壁面對流傳熱系數(shù)同理論計算值相近，誤差為5.45%，說明將軟件求解的對流傳熱系數(shù)用于模溫分析可靠可信.

3? ?實際零件模溫分析試驗與仿真

3.1? ?熱成形試驗

熱沖壓流程如圖5所示，首先根據(jù)目標零件形狀對成形件進行板料線估計后經(jīng)落料機完成落料. 由于未經(jīng)處理的板料在加熱過程中產(chǎn)生的氧化皮會影響成形性及淬火效果，所以對板料進行防氧化處理后，將板料推進加熱爐加熱300 s至930 ℃，保溫4～6 min以保證板料能夠均勻奧氏體化;經(jīng)機械手完成板料進模過程，本例實際入模溫度為687 ℃，在靠近板料過程中上模運動速度為600 mm/s，成形速度為120 mm/s;冷卻水入口平均溫度為8 ℃，出口平均溫度為15 ℃，出入口壓差為0.3 MPa;單腔理論沖壓力150 T，模具初溫20 ℃，板料厚度1.5 mm，保壓時間8 s.

3.2? ?熱成形仿真

3.2.1? ?有限元模型創(chuàng)建

目前冷卻管道的開設(shè)方式主要分為：鑄造式、鑲拼式及鉆孔式[12]. 鑄造式可以實現(xiàn)模具隨形冷卻，但制造難度大、成本高;鉆孔式加工難度低，但加工方式受限，只能直孔加工. 本文以某車A柱邊梁模具為例，采用鉆孔、鑲拼相結(jié)合的方式進行冷卻管道設(shè)計，可以在降低加工難度的同時，在一定程度上實現(xiàn)“隨形”冷卻. 仿真分析時水管和模面涉及流體計算及界面?zhèn)鳠?，故?yīng)采用較小的單元尺寸以提高求解精度，同時應(yīng)適當增大側(cè)壁和底面的單元尺寸以減少單元總數(shù).單元類型均為三角形，綜合考慮求解精度及計算效率，網(wǎng)格具體尺寸見表2.

流速及模溫分析中所使用的四面體網(wǎng)格均基于封閉的殼網(wǎng)格區(qū)域自動生成，網(wǎng)格模型分別如圖6和圖7所示.

3.2.2? ?冷卻系統(tǒng)流場分析

熱沖壓模具的散熱主要依賴于模具內(nèi)部的冷卻系統(tǒng)：冷流體快速流經(jīng)管道，高熱板料傳遞到模具中的熱量通過同冷流體的對流傳熱被帶走. 如果熱成形模具冷卻系統(tǒng)設(shè)計不合理，連續(xù)沖壓時模具體中將產(chǎn)生嚴重的熱聚集，使得板料對應(yīng)區(qū)域奧氏體向馬氏體的轉(zhuǎn)化受到影響，導致力學性能變差. 同時也更容易出現(xiàn)因冷卻速率不一致而產(chǎn)生翹曲變形. 因此熱成形模具冷卻系統(tǒng)的設(shè)計至關(guān)重要，能否通過模溫分析發(fā)現(xiàn)其設(shè)計缺陷以及驗證冷卻系統(tǒng)的設(shè)計能否滿足冷卻效率的要求主要取決于模具體同冷流體在壁面處對流傳熱系數(shù)的設(shè)置. 由于整個冷卻系統(tǒng)內(nèi)部的最高溫差為7 ℃，隨溫度變化的水的熱物性參數(shù)對傳熱系數(shù)的影響較小，為了進一步弱化水的參數(shù)變化的影響，本文在仿真分析過程中采用水在10 ℃時的參數(shù)進行流場分析，并計算對流傳熱系數(shù).

基于Fluent流場分析求解的對流傳熱系數(shù)以表格文件導出用于模溫分析，局部節(jié)點對流傳熱系數(shù)見表3，在相鄰區(qū)域傳熱系數(shù)值比較接近，不同區(qū)域會因流速不同導致傳熱系數(shù)差值較大.

由表3可知，冷卻水道表面對流傳熱系數(shù)分布存在差異，聯(lián)立公式（1）～（4）得到對流傳熱系數(shù)h和相關(guān)物理參數(shù)之間的關(guān)系如表4所示.

由表4可知，當流體類型選定后，可改變的只有流速和管道直徑，其中傳熱系數(shù)h正比于流速的0.8次方，因此流速的改變對傳熱系數(shù)的影響十分顯著. 表3和表4綜合表明，對熱成形模具進行傳熱分析研究時，應(yīng)考慮由管道結(jié)構(gòu)導致的流速差異的影響.

3.2.3? ?模溫分析

采用2種方法對熱沖壓模具進行模溫分析：基于Fluent流速分析求解得到的對流傳熱系數(shù)直接用于模溫分析;對基于Fluent求解得到的對流傳熱系數(shù)求平均值，再進行仿真分析. 兩者區(qū)別在于，方法1中的流場分析過程考慮了水道結(jié)構(gòu)的影響，對于復雜零件的熱沖壓模具或多腔模具而言模具內(nèi)部管道之間的轉(zhuǎn)接、變徑及流量分配等都會對流場產(chǎn)生影響，導致流場分布存在差異性，突出表現(xiàn)在局部管道流速偏低，流速的變化將導致水管表面上對流傳熱系數(shù)表現(xiàn)出明顯的差異，最終影響保壓結(jié)束后板料上的溫度分布. 圖8為冷卻水道上一節(jié)點處溫度隨時間變化曲線.

由圖8可知，首次沖壓時冷卻水主要用于給冷模具降溫，結(jié)束后節(jié)點溫度維持在22 ℃左右，接近模具初溫，表明首次沖壓過程中，冷卻水并未發(fā)揮太大的作用;第4次沖壓后節(jié)點溫度趨于接近，繼續(xù)多次沖壓后傳熱過程將趨于平衡，本文中選擇10次沖壓結(jié)束后的仿真值和實測值進行對比. 圖9為保壓結(jié)束后采用FLIR熱像儀拍攝的A柱邊梁大頭背側(cè)的溫度分布，高溫區(qū)域特征明顯，呈W形分布. A柱邊梁大頭背側(cè)溫度分布圖中呈W形分布的深色區(qū)域面積大且位于側(cè)壁中間，便于熱像儀測量，故在此區(qū)域選取監(jiān)測點，溫度監(jiān)測點分布如圖10所示.

用于模溫分析的網(wǎng)格質(zhì)量檢查結(jié)果見表5，其中模具四面體網(wǎng)格單元6 168 248個，模具傳熱分析主要研究模具及板料上的溫度分布，忽略模具自身在成形過程中的微形變，將上下模具定義為剛體. 設(shè)置求解時間步為5×10-5，容差為10-4，仿真中能量比率為1.004 5，誤差為0.45%.

4? ?結(jié)果與討論

采用2種模溫分析方法得到的仿真結(jié)果如圖11所示. 由圖11可知，基于Fluent流場分析的仿真結(jié)果呈W形分布和實際冷卻效果吻合較好. 在實際熱沖壓過程中導致板料溫度偏高的因素主要包括：模具鑲塊拼接處的管道間隙、水管到形面的距離以及冷卻水道結(jié)構(gòu)所導致的局部流速偏低. 熱沖壓模具一般采用獨立鑲塊加工再拼裝的方式制造，在拼接處管距較大，冷卻效果不佳;出于對模具強度考慮，凸凹模圓角處水管到模面的距離會適當增大，熱阻增大，熱量集聚在模具圓角處導致局部溫度偏高;模具鑲塊獨立加工，局部結(jié)構(gòu)的不合理設(shè)計以及較多的轉(zhuǎn)接接頭會使得整體水道結(jié)構(gòu)復雜、沿程壓力損失大，導致的局部流速偏低而形成局部高溫區(qū)域.

水管壁面處鄰近區(qū)域?qū)α鱾鳠嵯禂?shù)相近，忽略單元尺寸可認為傳熱系數(shù)是連續(xù)變化的，故鄰近區(qū)域溫度值變化不大，10次沖壓結(jié)束后，在仿真模型中選取和實際成形板料位置相近區(qū)域內(nèi)的多個節(jié)點，提取節(jié)點溫度并求其平均值，作為監(jiān)測點對應(yīng)的溫度仿真值，用于同實際板料上12個監(jiān)測點的溫度實測值進行對比. 采用2種分析方法的仿真值和實際測量結(jié)果對比如圖12所示.

圖12中，仿真1表示基于Fluent流場分析的仿真結(jié)果，仿真2表示基于對流傳熱系數(shù)平均值的仿真結(jié)果，對應(yīng)誤差分別為誤差1和誤差2. 由圖12可知，仿真1中的結(jié)果和實際測量結(jié)果一致性較好，最大誤差為3.6 ℃，平均誤差為0.96 ℃;仿真2中最大誤差為9.1 ℃，平均誤差為4.96 ℃. 相比于仿真1及實測結(jié)果，仿真2中整體溫度下降，表明熱沖壓模具在監(jiān)測區(qū)域水道流速偏低，導致模具和冷卻水的傳熱減弱，這是溫度實測值比較高的原因. 而仿真2由于對傳熱系數(shù)求平均后增大了監(jiān)測區(qū)域的傳熱系數(shù)值，傳熱能力增強，使溫度仿真值下降.

5? ?結(jié)? ?論

1）基于Fluent流速分析求解的對流傳熱系數(shù)和理論計算值相比誤差較小，為定義模具同冷卻水在壁面處的對流傳熱提供了新途徑，特別是針對復雜形狀零件及多腔模具的傳熱仿真分析.

2）采用基于Fluent流場分析求解對流傳熱系數(shù)的模溫分析方法能夠考慮由管道結(jié)構(gòu)所導致的流場差異對冷卻效果的影響，仿真中高溫區(qū)域及12個溫度監(jiān)測點的溫度值和實際結(jié)果吻合較好，12個監(jiān)測點最大誤差為3.6 ℃，平均誤差為0.96 ℃，相比于采用平均傳熱系數(shù)的仿真方法，仿真精度得到提高.

參考文獻

[1]? ? 張德軍. 熱成型高強鋼技術(shù)及成型模具的分析[J]. 內(nèi)燃機與配件，2019（19）：71—73.

ZHANG D J. Research on hot stamping high strength steel technology and forming die[J]. Internal Combustion Engine & Parts，2019（19）：71—73. （In Chinese）

[2]? ? 劉雪飛，賈卓翰，張謙，等. 超高強度鋼板熱沖壓成形同步淬火相變不均勻性影響因素分析[J]. 塑性工程學報，2019，26（2）：245—252.

LIU X F，JIA Z H，ZHANG Q，et al. Analysis on influencing factors of phase transformation inhomogeneity of synchronous quenching in hot stamping for ultra-high-strength steel plates[J]. Journal of Plasticity Engineering，2019，26（2）：245—252. （In Chinese）

[3]? ? 黃超群. 模具溫度對汽車B柱變強度熱成形零件過渡區(qū)的影響[J]. 鍛壓技術(shù)，2019，44（7）：52—57.

HUANG C Q. Influence of die temperature on transition zone of hot stamping parts with variable strength for automotive B-pillar[J]. Forging and Stamping Technology，2019，44（7）：52—57. （In Chinese）

[4]? ? 劉迪輝，萬雨嫻，張文甲. 熱沖壓模具冷卻系統(tǒng)流速均勻優(yōu)化研究[J]. 鍛壓技術(shù)，2017，42（12）：104—108.

LIU D H，WAN Y X，ZHANG W J. Research on flow rate uniform optimization of hot stamping die cooling system[J]. Forging and Stamping Technology，2017，42（12）：104—108. （In Chinese）

[5]? ? LAI C G，MAN C，WEN K P，et al. Numerical study of cooling characteristics for blanking in hot stamping [J]. Hydromechatronics Engineering，2015，43（6）：27—31.

[6]? ? 辛志宇，高乃平，朱彤. 熱沖壓模具冷卻系統(tǒng)數(shù)值模擬研究[J]. 熱加工工藝，2012，41（1）：170—174.

XIN Z Y，GAO N P，ZHU T. Numerical simulation investigation on cooling system of die in hot stamping? ?process[J]. Hot Working Technology，2012，41（1）：170—174. （In Chinese）

[7]? ? 薛克敏，孫大智，李萍，等. 22MnB5超高強鋼熱沖壓成形工藝及試驗[J]. 中國機械工程，2017，28（12）：1498—1503.

XUE K M，SUN D Z，LI P，et al. Hot stamping forming process and experiments of 22MnB5 ultra-high strength steels[J]. China Mechanical Engineering，2017，28（12）：1498—1503. （In Chinese）

[8]? ? 鄒偉，張立強，暫志偉，等. 超高強度鋼熱沖壓成形過程仿真研究[J]. 鍛壓技術(shù)，2018，43（8）：75—79.

ZOU W，ZHANG L Q，ZAN Z W，et al. Simulation study on hot stamping process of ultra high strength steel[J]. Forging and Stamping Technology，2018，43（8）：75—79. （In Chinese）

[9]? ? 張輝，潘愛瓊，李世云. BR1500HS板料熱沖壓成形工藝參數(shù)影響分析與試驗研究[J]. 熱加工工藝，2019，48（19）：97—99.

ZHANG H，PAN A Q，LI S Y. Analysis and experimental study on process parameters of hot stamping of BR1500HS sheet metal[J]. Hot Working Technology，2019，48（19）：97—99. （In Chinese）

[10]? GUO Y H，MA M T，F(xiàn)ANG G，et al. Numerical simulation of hot stamping of side impact beam[J]. Engineering Sciences，2012，10（6）：62—66.

[11]? 張興中，黃文，劉慶國. 傳熱學[M]. 北京：國防工業(yè)出版社，2011：89—95.

ZHANG X Z，HUANG W，LIU Q G. Heat transfer[M]. Beijing：National Defense Industry Press，2011：89—95. （In Chinese）

[12]? 王東生. 高強度鋼板熱沖壓成形模具設(shè)計規(guī)范[J]. 中國金屬通報，2019（4）：105.

WANG D S. Design specification of hot stamping forming die for the high strength steel plate[J]. China Metal Bulletin，2019（4）：105. （In Chinese）

ds：hot stamping;flow fields;convective heat transfer coefficients;cooling systems;temperature distribution