基于流體動(dòng)力學(xué)的熱壓罐框架式模具溫度場優(yōu)化

郭靖宇， 岳光， 呂佳鎂， 任琳， 潘玉田

(1.莫斯科鮑曼國立技術(shù)大學(xué) 特種制造系, 俄羅斯 莫斯科 105005;2.中北大學(xué) 機(jī)電工程學(xué)院, 山西 太原 030051; 3.太原工業(yè)學(xué)院, 山西 太原 030008)

0 引言

近年來，隨著航空業(yè)的快速發(fā)展，相較于傳統(tǒng)材料，具有更大比強(qiáng)度和比模量，疲勞壽命更長的復(fù)合材料在航空航天工業(yè)中有著更加廣闊的應(yīng)用范圍和使用場景，尤其在一些飛行器關(guān)鍵部位，其逐漸取代傳統(tǒng)材料的趨勢愈發(fā)凸顯。

復(fù)合材料具有材料與產(chǎn)品形狀同步生成的特點(diǎn)，并且在后期幾乎不可以進(jìn)行機(jī)械加工。因此，在生產(chǎn)復(fù)合材料的過程中，必須保證材料可以達(dá)到設(shè)計(jì)要求。作為廣泛使用的整體成型復(fù)合材料構(gòu)建技術(shù)，熱壓罐成型技術(shù)還存在著由于溫度分布不均勻?qū)е虏牧显诠袒A段出現(xiàn)固化不均勻的現(xiàn)象，即固化梯度。這將在材料內(nèi)部引起殘余應(yīng)力、變形，甚至是出現(xiàn)氣泡等致命問題，因此研究優(yōu)化熱壓罐成型框架式模具的溫度場具有重要意義。

熱壓罐利用電熱絲或者電熱管等加熱方式將工作氣體加熱至工藝溫度，并經(jīng)過位于工作部位外的風(fēng)扇驅(qū)動(dòng)，以一定的速度流至工作部位。在工作部位內(nèi)通過工作氣體的溫度和流動(dòng)時(shí)的壓力對(duì)復(fù)合材料進(jìn)行作用，在理想條件下，復(fù)合材料基體將以均勻的速度固化，繼而獲得均勻的力學(xué)特性。

在國內(nèi)外的相關(guān)研究中，研究者集中研究了框架式結(jié)構(gòu)對(duì)溫度分布的影響：Gniatczyk等指出，使用框架式模具可以有效提高模具的導(dǎo)熱性；花蕾蕾等研究了框架式結(jié)構(gòu)的散熱孔形狀、支撐板結(jié)構(gòu)的有無、厚度等對(duì)溫度場分布均勻性的影響；Han等提出通過改變模具在熱壓罐中的放置角度來改善溫度分布的方法；張鋮等指出框架結(jié)構(gòu)的溫度分布均勻性隨著熱壓罐內(nèi)熱流的速度的增加而有所改善；張旭升等利用流體力學(xué)分析軟件Fluent對(duì)多種形狀的風(fēng)道進(jìn)行模擬計(jì)算，得出T字形風(fēng)道可以改善模具型面溫度場分布均勻性的結(jié)論；趙一鳴等從熱阻與材料厚度的關(guān)系出發(fā)，研究了模具型板厚度對(duì)型面溫度分布的影響；羅登峰研究樹脂基復(fù)合材料在固化過程中自身因化學(xué)反應(yīng)而產(chǎn)生的熱量及溫度對(duì)自身溫度場的影響；Hu等提出，通過在框架式結(jié)構(gòu)內(nèi)部關(guān)鍵部位使用不同傳熱系數(shù)的材料，來調(diào)節(jié)傳遞到模具型面的溫度。

上述研究者利用改變模具型面厚度分布、使用不同傳熱系數(shù)的材料構(gòu)建模具以及更換模具內(nèi)部框架結(jié)構(gòu)等方式，來改善熱壓罐固化技術(shù)在復(fù)合材料成型過程中型面溫度分布的不均勻性。

在空間允許的范圍內(nèi)，氣體由壓力大的區(qū)域流向壓力小的區(qū)域，并且流速隨著氣壓差的增大而增大。因此可以設(shè)計(jì)一種通過逐漸縮小工作氣體流通面積的方式來減緩工作氣體流速下降，繼而改善模具型面溫度場分布均勻性的新模具。

綜上所述，本文基于流體流速對(duì)對(duì)流換熱的影響，提出一種具有經(jīng)濟(jì)性的結(jié)構(gòu)優(yōu)化方法，即在框架式模具內(nèi)添加導(dǎo)流板，來改善熱壓罐內(nèi)框架式結(jié)構(gòu)溫度分布的不均勻性。并使用多物理場仿真軟件COMSOL來對(duì)上凸型、直板型和下凹型3種變截面構(gòu)型模具進(jìn)行對(duì)比分析，為模具結(jié)構(gòu)設(shè)計(jì)提供參考。

1 數(shù)學(xué)模型的創(chuàng)建

1.1 熱傳遞方程

熱壓罐內(nèi)的工作原理是利用加熱至工藝溫度的工作氣體在工作壓力下以一定速度流經(jīng)模具，在模具表面與其進(jìn)行熱交換，從而達(dá)到固化復(fù)合材料基體的目標(biāo)。同時(shí)由于熱壓罐的結(jié)構(gòu)限制了熱源與模具之間的熱輻射，故而可以假設(shè)在熱壓罐內(nèi)僅有對(duì)流換熱和熱傳導(dǎo)兩種換熱方式。此外假設(shè)熱壓罐內(nèi)氣體為理想氣體，從而方便描述其密度隨溫度、壓力的變換。因此可以建立方程來對(duì)此過程進(jìn)行數(shù)學(xué)描述：

1)流體域?qū)α鲹Q熱。

質(zhì)量守恒方程：

(1)

式中：為熱壓罐內(nèi)工作氣體的密度；為工作氣體流動(dòng)速度。

動(dòng)量守恒方程：

(2)

(3)

(4)

式中：、、分別表示熱壓罐內(nèi)平行于型面且垂直于圓罐軸線的方向、圓罐軸線的方向以及垂直于型面向上的方向；為工作氣體的動(dòng)力黏度；、、分別為動(dòng)量方程在軸、軸、軸方向上的廣義源項(xiàng)；為工作氣體的壓力。

能量守恒方程：

(5)

式中：為工作氣體的對(duì)流傳熱系數(shù)；為工作氣體的傳熱系數(shù)；為工作氣體的溫度；為模具內(nèi)部熱源，=0；為工作氣體的機(jī)械能由于黏性作用轉(zhuǎn)換成其他能量形式的部分。

理想氣體狀態(tài)方程：

=，

(6)

式中：為工作氣體體積；為物質(zhì)的量；為理想氣體狀態(tài)常數(shù)。

牛頓散熱方程：

(7)

式中：為工作氣體與模具間的傳熱功率；為工作氣體與模具間的傳熱系數(shù)；為模具溫度；為模具的熱容；為模具體積。

根據(jù)(7)式可知，流體的速度越快，溫度的分布越均勻，故可以通過逐步減小模具內(nèi)部氣體流經(jīng)的面積，工作氣體將加速流動(dòng)來滿足連續(xù)性原理(質(zhì)量守恒)，繼而補(bǔ)償因橫向筋板阻擋而損失的氣體速度，進(jìn)一步改善溫度的分布。

2) 模具內(nèi)部熱傳導(dǎo)：

(8)

式中：為模具的密度；為模具的傳熱系數(shù)。

1.2 邊界條件與初始值

在研究對(duì)流換熱的過程中，若取整個(gè)熱壓罐為研究對(duì)象，則模型的邊界條件為第1類邊界條件，即已給定相關(guān)邊界條件內(nèi)容，求解模型中的溫度分布，如圖1所示。

圖1 熱壓罐示意圖Fig.1 Schematic diagram of autoclave

在工作氣體入口處，氣體按照工藝流速和工藝溫度流入熱壓罐工作區(qū)，工藝流速為2.5 m/s，工藝溫度如圖2所示；在工作氣體出口處，存在工作壓力0.6 MPa. 熱壓罐側(cè)壁為絕熱壁。

圖2 工藝溫度曲線Fig.2 Process temperature curve

在研究固體換熱的過程中，取模具為研究對(duì)象，其邊界條件為模具與工作氣體所有的接觸面均發(fā)生熱交換。

圖3 熱壓罐框架式模具Fig.3 Frame mold for autoclave

2 建模與仿真研究

建立如圖3所示的熱壓罐框架式模具，其規(guī)格為1 700 mm×1 500 mm×400 mm，其中包含半徑25 mm的半圓形散熱孔和尺寸130 mm×220 mm的矩形散熱孔，散熱孔均勻分布在筋板上，筋板的厚度為4 mm，型面的厚度為7 mm.

熱壓罐工作部分的直徑為2.5 m，長度為7 m. 模具型面下部幾何中心位于熱壓罐幾何中心。熱壓罐內(nèi)的工作氣體為空氣，模具的材料是結(jié)構(gòu)鋼，材料的熱力學(xué)性能分別如表1和表2所示。

表1 空氣的熱力學(xué)性能Tab.1 Thermodynamic properties of air

表2 結(jié)構(gòu)鋼的熱力學(xué)性能Tab.2 Thermodynamic properties of structural steel

根據(jù)熱壓罐工作部分的尺寸和工藝流速，按照(9)式可計(jì)算出熱壓罐內(nèi)工作氣體的雷諾數(shù)：

(9)

式中：為熱壓罐垂直于工作氣體流動(dòng)方向的橫截面面積；為工作氣體的黏性系數(shù)。

由(9)式計(jì)算可得，=259 148≥12 000，繼而可判斷熱壓罐內(nèi)的工作氣體運(yùn)動(dòng)為湍流。

由于研究對(duì)象存在兩種換熱方式，并且與工作氣體的流動(dòng)有關(guān)，故在COMSOL軟件中設(shè)置兩個(gè)物理場，即湍流(-模型，為湍流動(dòng)能，為耗散率)物理場以及固體與流體傳熱物理場，并設(shè)置其耦合物理場為非等溫流動(dòng)。將框架式模具導(dǎo)入COMSOL軟件后，為簡化研究，將模型沿平面剖開，取其中一半為研究對(duì)象。按照12節(jié)中所述的邊界條件和初始值在COMSOL軟件中編輯相關(guān)內(nèi)容。

使用非結(jié)構(gòu)網(wǎng)格對(duì)模型進(jìn)行劃分，得到約50萬個(gè)單元，結(jié)果如圖4所示。在模具型面上表面按照文獻(xiàn)[20]實(shí)驗(yàn)中布置測溫電偶的位置在COMSOL軟件中添加域點(diǎn)探針用以收集采樣點(diǎn)溫度。表3和表4分別為實(shí)驗(yàn)和仿真中5個(gè)具有代表性測量點(diǎn)位的數(shù)據(jù)。

在相同工藝條件下，將本文仿真數(shù)據(jù)與文獻(xiàn)[20]中的實(shí)驗(yàn)數(shù)據(jù)比較。取仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)差別最大和最小點(diǎn)位具體分析，如圖5和圖6所示。由此可以證明本文仿真方法可以較為準(zhǔn)確地模擬實(shí)驗(yàn)過程。

圖4 非結(jié)構(gòu)化網(wǎng)絡(luò)的建立Fig.4 Unstructured mesh establishment

表3 實(shí)驗(yàn)數(shù)據(jù)

表4 仿真數(shù)據(jù)Tab.4 Simulated data K

圖5 誤差最大處數(shù)據(jù)Fig.5 Data at maximum error point

圖6 誤差最小處數(shù)據(jù)Fig.6 Data at minimum error point

在第1節(jié)所述的工藝條件下，對(duì)完整的工藝過程進(jìn)行仿真，通過提取仿真結(jié)果中型面約10 000個(gè)節(jié)點(diǎn)的溫度，可以得到型面溫度的標(biāo)準(zhǔn)差隨時(shí)間的變化規(guī)律：隨著工作氣體溫度的升高，型面表面溫度分布的不均勻性逐漸增加，在溫度升高的結(jié)束時(shí)達(dá)到最大值；后隨著保溫時(shí)間的增加，溫度分布的不均勻性迅速降低，其標(biāo)準(zhǔn)差值逐漸趨于0 K；之后隨著工作氣體的降溫，溫度標(biāo)準(zhǔn)差值開始增加，即型面溫度分布逐漸不均勻，如圖7所示。

圖7 溫度標(biāo)準(zhǔn)差隨時(shí)間的變化規(guī)律Fig.7 Standard deviation of temperature over time

根據(jù)1.1節(jié)所提出的理論方案，建立如圖8的直板型模具，即在原有的框架中添加一個(gè)沿著工作氣體流向方向逐漸上升的平面。其逐漸縮小了工作氣體的流動(dòng)面積，在理論上達(dá)到了增加模具尾部工作氣體流速的目的。同時(shí)在該原理的指導(dǎo)下，可以得到上凸型結(jié)構(gòu)(見圖9)和下凹型結(jié)構(gòu)(見圖10)。將上述結(jié)構(gòu)導(dǎo)入COMSOL軟件中，在垂直于工作氣體流向每一個(gè)筋板的半圓形散熱孔內(nèi)添加域點(diǎn)探針，以測量流經(jīng)型面下表面的工作氣體流速，具體位置如圖11所示。

圖8 直板型框架式模具Fig.8 “Straight” frame mold

圖9 上凸型框架式模具Fig.9 “Upward convex” frame mold

圖10 下凹型框架式模具Fig.10 “Downward concave” frame mold

圖11 域點(diǎn)探針的位置Fig.11 Position of domain point probe

進(jìn)行與未加導(dǎo)流板的框架式模具(以下簡稱標(biāo)準(zhǔn)型)相同內(nèi)容的仿真。

取其型面為研究對(duì)象，在型面溫度分布最不均勻時(shí)刻，即升溫結(jié)束時(shí)刻，3種類型的模具型面溫度分布分別如圖12、圖13、圖14所示。

圖12 直板型模具型面溫度分布Fig.12 Temperature distribution of “straight” mold surface

圖13 上凸型模具型面溫度分布Fig.13 Temperature distribution of “upward convex” mold surface

圖14 下凹型模具型面溫度分布Fig.14 Temperature distribution of “downward concave” mold surface

對(duì)型面上約10 000個(gè)節(jié)點(diǎn)的溫度變化進(jìn)行統(tǒng)計(jì)，計(jì)算同一時(shí)刻溫度分布的標(biāo)準(zhǔn)差，得到如圖15所示曲線。由圖15可以發(fā)現(xiàn)，在標(biāo)準(zhǔn)差最大處，即升溫結(jié)束時(shí)刻，直板型、下凹型和上凸型模具相較于標(biāo)準(zhǔn)型模具，均有效降低了型面溫度分布的不均勻性。

圖15 型面溫度分布標(biāo)準(zhǔn)差Fig.15 Standard deviation of surface temperature distribution

利用提前設(shè)置好的域點(diǎn)探針可以測量得到型面下表面工作氣體的流速，其在不同結(jié)構(gòu)中表現(xiàn)如圖16所示。

圖16 型面下表面工作氣體流速Fig.16 Working gas velocity on the lower surface of the profile

由圖15可知，增加了導(dǎo)流板后，模具型面溫度分布的不均勻性較標(biāo)準(zhǔn)型模具均有所改善，其中以直板型和下凹型尤為顯著。直板型模具平均標(biāo)準(zhǔn)差從原來的0.78 K降低至0.58 K，降幅約為25%；下凹型則降低至0.59 K，降幅約為24%。

由圖16可知，與標(biāo)準(zhǔn)型模具相比，工作氣體在3種添加導(dǎo)流板模具尾部的流速均有所上升。其中，上凸型模具最為明顯，但是在該型模具尾部，工作氣體流速有明顯下降的趨勢，這與要提高模具末端的氣體流速以補(bǔ)償氣體與模具前端熱交換而損失熱量的優(yōu)化目的不一致，故而無法顯著提高型面溫度分布的均勻性。與之相反，直板型和下凹型模具在中段部位存在工作氣體流速較慢的情況，但在模具尾部，流速有明顯的提高，因此可以達(dá)到優(yōu)化目的。

故在工程實(shí)踐中，只需在原有框架式模具中加裝沿工作氣體流動(dòng)方向逐漸向上的導(dǎo)流板，即可提高25%型面溫度分布的均勻性，同時(shí)該方案具有較高的經(jīng)濟(jì)性，也易于在生產(chǎn)中實(shí)現(xiàn)。

在實(shí)際生產(chǎn)過程中，需要進(jìn)行固化的復(fù)合材料大多不是單一的平面，而是存在一些曲面。以火箭燃料罐為例，在使用框架式成型模具的條件下，與模具不接觸的位置，幾乎不會(huì)出現(xiàn)受到因模具結(jié)構(gòu)阻擋而導(dǎo)致溫度分布不均勻的情況。因此研究的關(guān)鍵部位為罐體與模具的接觸部位，遂建立如圖17所示的曲面結(jié)構(gòu)以及相關(guān)的框架式模具，其添加導(dǎo)流板后的結(jié)構(gòu)如圖18所示。

圖17 曲面型面框架式模具Fig.17 Curved surface frame mold

圖18 添加導(dǎo)流板的曲面型面框架式模具Fig.18 Curved surface frame mold with baffle plate

經(jīng)過與第2節(jié)中所述條件相同的仿真過程，并對(duì)相關(guān)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，可得如圖19所示型面溫度分布標(biāo)準(zhǔn)差隨時(shí)間的變化趨勢。從圖19中可以看出，添加導(dǎo)流板后，曲面型面溫度分布標(biāo)準(zhǔn)差有明顯的降低。

圖19 曲面型面溫度分布標(biāo)準(zhǔn)差Fig.19 Standard deviation of temperature distribution on curved surface

綜上所述，對(duì)于以曲面為主的型面，添加導(dǎo)流板可有效改善溫度分布的均勻性，從而擴(kuò)展了添加導(dǎo)流板框架式模具的適用性。

3 結(jié)論

基于牛頓散熱方程所展現(xiàn)的固體與流體換熱過程中換熱效率與流體流速的關(guān)系，本文開展在熱壓罐框架式模具內(nèi)添加沿工作氣體流向逐漸升高的導(dǎo)流板優(yōu)化研究。通過對(duì)比已有實(shí)驗(yàn)數(shù)據(jù)，驗(yàn)證了使用COMSOL軟件對(duì)熱壓罐框架式模具型面溫度分布模擬的準(zhǔn)確性。經(jīng)COMSOL軟件建模仿真，在控制工作氣體流出模具面積的條件下，比較了3種添加了導(dǎo)流板的模具型面溫度分布的不均勻特性——型面溫度分布標(biāo)準(zhǔn)差，得出直板型模具可以有效地降低型面溫度分布的不均勻性，繼而達(dá)到改善復(fù)合材料基體固化時(shí)溫度梯度過大的目的，提高了復(fù)合材料內(nèi)部的均勻特性。同時(shí)擴(kuò)展研究了添加導(dǎo)流板對(duì)曲面型面的影響，證明了導(dǎo)流板對(duì)于曲面型面的模具亦可以有效改善表面的溫度分布。

本文只探討了添加3種不同類型的導(dǎo)流板對(duì)工作氣體流速和型面溫度分布的影響，對(duì)于導(dǎo)流板具體的安裝位置尚未進(jìn)行深入研究?；诂F(xiàn)有理論，可通過更改導(dǎo)流板的位置以及斜率來針對(duì)模具尾部氣體流速下降較快的現(xiàn)狀進(jìn)行進(jìn)一步優(yōu)化和分析。