程傳峰，金 明，王項(xiàng)如，朱英霞，程一峰，盤朝奉，王 園

(1. 安徽新富新能源科技股份有限公司，安徽 安慶 246001)

(2. 安徽環(huán)新集團(tuán)股份有限公司，安徽 安慶 246001)

(3. 江蘇大學(xué)機(jī)械工程學(xué)院，江蘇 鎮(zhèn)江 212013)

(4. 江蘇大學(xué) 汽車工程研究院，江蘇 鎮(zhèn)江 212013)

1 前 言

目前，在航空航天、汽車、通訊、軍工等領(lǐng)域，大量設(shè)備正向著微型化、高功率和結(jié)構(gòu)高密度的方向發(fā)展[1]，因此工作時(shí)的熱流密度遠(yuǎn)大于常規(guī)尺度的設(shè)備，如果散熱不夠理想，其工作性能和壽命會(huì)受到嚴(yán)重影響[2]。因此，傳熱性能優(yōu)異且同時(shí)具備節(jié)能、降本、環(huán)保優(yōu)勢(shì)的薄壁多孔微小通道扁管被越來越多的研究者和制造廠商應(yīng)用于微型化、高功率和結(jié)構(gòu)高密度集成化設(shè)備的散熱系統(tǒng)中。

薄壁多孔微小通道扁管是一種采用精煉鋁棒，通過熱擠壓-熔焊工藝和表面噴鋅防腐處理制造成形的薄壁多孔扁形管[3]。由于薄壁微小通道扁管具有十分突出的環(huán)保、增效、節(jié)能、降本4大優(yōu)勢(shì)，因此早在1996年就在汽車空調(diào)系統(tǒng)中獲得應(yīng)用。與傳統(tǒng)的管翅式換熱器相比，薄壁多孔微小通道扁管換熱器的傳熱效率可以提高14%～33%。在獲得相同制冷效果的前提下，所需制冷劑的量可減少35%[4]。

薄壁多孔微小通道扁管的結(jié)構(gòu)精細(xì)復(fù)雜，加工制造難度較高，常用的材料主要有1050、1060、1100、1A97、3003、3F03、3102和3103鋁合金[5，6]。

多孔微小通道的導(dǎo)熱性能與流道的尺寸和橫截面形狀密切相關(guān)[7]，因此很多學(xué)者針對(duì)多孔微小通道扁管的空間結(jié)構(gòu)展開了研究。夏國棟團(tuán)隊(duì)[8]和謝忱創(chuàng)[9]均研究發(fā)現(xiàn)，波形結(jié)構(gòu)的薄壁多孔微小通道扁管能大大提高新能源汽車電池組的散熱能效。最新的研究成果也表明，較之平行微小通道，波形微小通道能使流通的冷卻液再循環(huán)和回流，進(jìn)而強(qiáng)化扁管的能效傳遞，且相同工況下，波形微小通道扁管的溫度分布更均勻[10-12]。

現(xiàn)有研究說明，多孔微小通道波形扁管將是解決新能源汽車鋰電池組和機(jī)車電控元件這類發(fā)熱密集、溫度敏感型零部件恒溫問題的絕佳材料[10，11]。然而由于缺乏薄壁多孔微小通道扁管的波形成形工藝與技術(shù)，上述研究成果依舊停留在理論研究階段，迫切需要對(duì)薄壁多孔微小通道扁管的波形成形工藝展開研究。

李大永團(tuán)隊(duì)研究了多孔微小圓通道扁管[13]、多層平行流微小通道換熱器[4]和折疊式薄壁多孔扁管[14]的繞彎成形能力。繞彎成形與本文研究的波形沖壓同屬管類的二次塑性成形工藝。李大永團(tuán)隊(duì)的研究結(jié)果表明：① 薄壁多孔微小通道扁管的彎曲半徑和截面尺寸對(duì)繞彎成形結(jié)束后微小通道的截面變形有重要影響；② 截面變形是薄壁多孔微小通道扁管最為顯著的彎曲成形缺陷，也是對(duì)其導(dǎo)熱性能影響最嚴(yán)重的缺陷，且與常規(guī)管材不同的是，其通道的縱截面變形比橫截面變形更顯著。上述研究結(jié)果與工藝方法對(duì)本研究具有較好的借鑒意義，但是其研究對(duì)象接近板材，獲得的數(shù)據(jù)與經(jīng)驗(yàn)規(guī)律不能直接用于管材。

本文以適用于高密度集成化設(shè)備散熱系統(tǒng)的3003Al-H14薄壁多孔微小通道扁管為研究對(duì)象，對(duì)其波形沖壓工藝過程中成形尺寸對(duì)截面變形的影響規(guī)律展開了研究。本研究對(duì)波形微小通道的精確塑性成形發(fā)展具有很好的科研價(jià)值和技術(shù)指導(dǎo)意義。

2 建立3003Al-H14薄壁微小通道扁管的波形沖壓成形有限元模型

2.1 3003Al-H14薄壁微小通道扁管的成形尺寸與材料參數(shù)

3003Al-H14薄壁微小通道扁管的橫截面存在“筋”與孔，孔的個(gè)數(shù)m通常很多。橫截面寬度l遠(yuǎn)大于高度h，因此呈現(xiàn)扁狀。管材壁厚t極薄，通常情況下t≤0.8 mm。如圖1所示，薄壁微小通道扁管的成形尺寸參數(shù)還包括：橫截面孔的內(nèi)徑寬度w，縱截面總長L，內(nèi)彎曲面波形半徑r(至彎曲中性層)，與外彎曲面波形半徑R(至彎曲中性層)，參數(shù)值如表1所示。

表1 3003Al-H14薄壁微小通道扁管的成形尺寸參數(shù)值

圖1 3003Al-H14薄壁微小通道扁管的成形尺寸示意圖Fig.1 The geometric dimensions of the 3003Al-H14 thin-walled micro-channel flat tube

2.2 3003Al-H14微小通道薄壁扁管的材料本構(gòu)模型

為了獲得材料的應(yīng)力-應(yīng)變曲線，將3003Al-H14薄壁微小通道扁管沿橫截面寬度方向切割為完全相同的3部分，分別作為拉伸試樣進(jìn)行了拉伸實(shí)驗(yàn)。實(shí)驗(yàn)按照國標(biāo)GB/T 228.1—2010《金屬材料室溫拉伸試驗(yàn)方法》進(jìn)行，拉伸速度設(shè)置為1 mm/min。拉伸后的試樣斷裂狀態(tài)如圖2所示。

圖2 3003Al-H14薄壁微小通道扁管拉伸斷裂試樣照片F(xiàn)ig.2 Fracture state of the tensile specimen of the 3003Al-H14 thin-walled micro-channel flat tube

拉伸實(shí)驗(yàn)獲得的真應(yīng)力-應(yīng)變曲線如圖3所示，材料力學(xué)性能參數(shù)如表2所示。

表2 3003Al-H14薄壁微小通道扁管的材料力學(xué)性能參數(shù)

圖3 3003Al-H14微小通道扁管拉伸真應(yīng)力-應(yīng)變曲線Fig.3 The ture stress-strain curve of the tensile test of the 3003Al-H14 micro-channel flat tube

采用各向同性硬化模型描述3003Al-H14材料的應(yīng)力應(yīng)變關(guān)系。其屈服條件f如公式(1)所示：

(1)

其中，s為偏應(yīng)力，εp為等效塑性應(yīng)變，硬化參數(shù)A、K和n的取值來自圖3中應(yīng)力應(yīng)變曲線擬合，結(jié)果如表2所示。

2.3 3003Al-H14薄壁微小通道扁管的波形沖壓-回彈有限元模型

所有的截面變形數(shù)據(jù)都取自回彈發(fā)生之后，所以建立了3003Al-H14微小通道薄壁扁管的波形沖壓-回彈有限元模型。該模型包括沖壓成形和回彈2步，分別采用動(dòng)態(tài)顯示算法和靜態(tài)隱式算法，如圖4所示。沖壓過程的結(jié)構(gòu)包括上模、下模和扁管這3部分，如圖4a所示。扁管與上模接觸的成形面為內(nèi)彎曲面，與下模接觸的成形面為外彎曲面。扁管與上、下模之間的摩擦條件，均為機(jī)械油潤滑。扁管的網(wǎng)格采用S4R殼單元。

圖4 薄壁微小通道扁管波形沖壓-回彈有限元模型：(a)沖壓過程，(b)回彈過程Fig.4 The finite element model of corrugated stamping-springback processes of thin-walled micro-channel flat tube：(a) stamping process，(b) springback process

采用3003Al-H14薄壁微小通道扁管的波形沖壓實(shí)驗(yàn)驗(yàn)證所建有限元模型的可靠性。模擬邊界條件設(shè)置和實(shí)驗(yàn)條件對(duì)比如表3所示。

表3 模擬條件設(shè)置和實(shí)驗(yàn)條件對(duì)比

圖5對(duì)比了模擬和實(shí)驗(yàn)獲得的薄壁微小通道扁管波形沖壓件。由圖5a和5b可知，所建立的有限元模型可以準(zhǔn)確模擬管坯橫、縱截面的起皺狀態(tài)、截面變形狀態(tài)。圖5c對(duì)比了某一特定縱截面上的截面變形率預(yù)測情況與實(shí)驗(yàn)數(shù)據(jù)，發(fā)現(xiàn)有限元模擬的平均預(yù)測誤差為16.82%，最大預(yù)測誤差為24.32%，均在合理誤差范圍內(nèi)。綜上，認(rèn)為所建3003Al-H14薄壁微小通道扁管的波形沖壓-回彈有限元模型能夠較為可靠地預(yù)測扁管沖壓過程中的截面變形。

圖5 薄壁微小通道扁管波形沖壓的實(shí)驗(yàn)和模擬結(jié)果對(duì)比：(a) 實(shí)驗(yàn)截面變形狀態(tài)，(b) 模擬截面變形狀態(tài)，(c) 代表縱截面PC上截面變形率的模擬預(yù)測值和實(shí)驗(yàn)值對(duì)比Fig.5 Comparisons of experimental and simulative results of corrugated stamping of thin-walled micro-channel flat tubes：(a) experimental wrinkle state，(b) simulative wrinkle state，(c) comparison of the simulative and experimental values of the cross-section deformation rate on the representative longitudinal section PC

3 3003Al-H14薄壁微小通道扁管波形沖壓的截面變形描述及模擬研究條件

本研究中截面變形量(率)的取值來自兩橫、兩縱四條代表截面，其位置和截面上的節(jié)點(diǎn)編號(hào)如圖6所示?？梢钥吹酱頇M截面PP1是波峰橫截面，CC1是波谷橫截面；代表縱截面PC是位于邊緣孔對(duì)稱線上的縱截面，P1C1是位于中央孔對(duì)稱線上的縱截面。由于肋板的截面變形量非常小，因此橫、縱截面上代表節(jié)點(diǎn)i的選取皆對(duì)應(yīng)于各孔的對(duì)稱線，如圖6所示。

圖6 代表橫截面PP1、CC1和代表縱截面PC、P1C1的截面位置，及其截面上的節(jié)點(diǎn)編號(hào)順序Fig.6 The representative cross sections PP1，CC1 and the representative longitudinal sections PC，P1C1，and the node number sequence on the sections

(2)

(3)

其中，Δh為截面變形量，h1為微小通道薄壁波形扁管沖壓成形后的截面高度，i為代表截面上的節(jié)點(diǎn)編號(hào)。

為了揭示3003Al-H14薄壁微小通道扁管波形沖壓的成形尺寸-截面變形作用規(guī)律，結(jié)合實(shí)際生產(chǎn)中所需扁管的成形尺寸，設(shè)置模擬研究范圍如表4所示。表4中的基礎(chǔ)值是1.3節(jié)中的實(shí)驗(yàn)設(shè)置值，也是模擬研究的基礎(chǔ)數(shù)據(jù)。成形尺寸包括截面高度h、管坯壁厚t、內(nèi)彎曲面波形半徑r，以及內(nèi)、外面相對(duì)彎曲半徑r/R的同步縮放系數(shù)β。其中β的含義如式(4)：

(4)

表4 模擬研究成形尺寸值

4 結(jié)果與討論

4.1 3003Al-H14薄壁微小通道扁管波形沖壓的截面變形分布特征

圖7給出了使用表4基礎(chǔ)模擬研究數(shù)據(jù)的情況下，薄壁微小通道扁管的截面變形率分布情況。從圖7可以看出：① 代表橫截面PP1(波峰)和CC1(波谷)的截面變形率δh分布趨勢(shì)完全一致，邊緣孔的變形最為嚴(yán)重，其它孔的變形率相對(duì)較小且基本一致。② 代表縱截面PC(邊緣孔縱截面)和P1C1(中央孔縱截面)的截面變形率δh分布趨勢(shì)完全一致，呈現(xiàn)兩端點(diǎn)高、中央段低的波動(dòng)分布趨勢(shì)。以上說明，薄壁微小通道扁管的截面變形情況復(fù)雜，橫、縱截面變形分布規(guī)律差異較大。由于橫截面PP1的平均截面變形率要高于CC1，縱截面PC的平均截面變形率要高于P1C1，因此將PP1和PC截面作為后續(xù)研究對(duì)象。

圖7 代表橫截面PP1、CC1和代表縱截面PC、P1C1的截面變形率Fig.7 Distributions of the section deformation rate of the representative cross sections PP1，CC1 and the representative longitudinal sections PC，P1C1

4.2 截面高度對(duì)截面變形率的作用規(guī)律

圖8a給出了不同截面高度h下的波峰橫截面變形情況?？梢钥闯霎?dāng)橫截面的寬度不變時(shí)，h值越大，橫截面的中空度越高，截面剛度越差，因而變形程度越大。當(dāng)h達(dá)到2.5 mm時(shí)，由于筋的折彎并不顯著，整體橫截面尚能保持一定的形狀精度，但邊緣孔率先出現(xiàn)了較為顯著的截面塌陷，這是由于邊緣孔的材料流動(dòng)較之其它孔要自由。而當(dāng)h=4 mm時(shí)，筋發(fā)生顯著的彎折，導(dǎo)致整體橫截面塌陷嚴(yán)重。當(dāng)h=7 mm時(shí)，孔的截面剛度進(jìn)一步變差，筋的彎折加劇，整體的橫截面塌陷也愈加嚴(yán)重，如圖8a和8b所示。

圖8 不同截面高度(h)下扁管截面變形情況：(a) 不同h下的波峰橫截面對(duì)比圖，(b) h=7 mm時(shí)的扁管截面變形圖Fig.8 Sectional deformation of flat tubewith different section heights(h)：(a) comparison of cross-sections of wave crests with different h，(b) cross-section deformation of flat tube with h=7 mm

圖9 不同截面高度(h)下的扁管截面變形率(δh)分布和平均截面變形率變化趨勢(shì)：(a) PP1截面的δh分布，(b) PC截面的δh分布，(c) PP1和PC截面的隨h的變化趨勢(shì)Fig.9 The section deformation rate (δh) distributions and average section deformation rate varation trends of flat tubes with different heights (h)：(a) δh distribution of PP1 section，(b) δh distribution of PC section，(c) variation trends of PP1 and PC sections with h

4.3 壁厚尺寸對(duì)截面變形率的作用規(guī)律

圖10 不同管壁厚度(t)下的扁管截面變形率(δh)分布和平均截面變形率變化趨勢(shì)：(a) PP1截面的δh分布，(b) PC截面的δh分布，(c) PP1和PC截面的隨t的變化趨勢(shì)Fig.10 The section deformation rate (δh) distributions and average section deformation rate variation trends of flat tubes with different wall thicknesses (t)：(a) δh distribution of PP1 section，(b) δh distribution of PC section，(c) variation trends of PP1 and PC sections with t

4.4 內(nèi)外面相對(duì)彎曲半徑對(duì)截面變形率的作用規(guī)律

圖11 不同內(nèi)彎曲面波形半徑r下的平均截面變形率變化趨勢(shì)Fig.11 Variation trends of average section deformation rate with wave radius r of inner curved surface

圖12 r/R=10/6時(shí)，平均截面變形率隨同步縮放系數(shù)(β)的變化趨勢(shì)：(a) PP1和PC截面的隨β的變化趨勢(shì)，(b)不同β下的波峰橫截面對(duì)比Fig.12 Variation trends of average section deformation rate with synchronous zoom factor (β)，when r/R=10/6：(a) variation trends of PP1 and PC sections with β，(b) comparison of cross-sections of wave crests with different β

5 結(jié) 論

(1)薄壁微小通道扁管的截面變形情況復(fù)雜，橫、縱截面呈現(xiàn)完全不同的截面變形分布趨勢(shì)。波峰橫截面和波谷橫截面的截面變形情況基本一致，即：邊緣孔的截面變形量通常較大，其余孔的截面變形相對(duì)較小且基本一致。邊緣孔縱截面和中央孔縱截面的截面變形率分布趨勢(shì)基本一致，呈現(xiàn)波峰、波谷高，中間段低的波動(dòng)特點(diǎn)。