基于Fluent的紙漿模塑熱壓定型加熱板溫度均勻性分析

吳上生 胡錦榕 周運(yùn)岐

基于Fluent的紙漿模塑熱壓定型加熱板溫度均勻性分析

吳上生1胡錦榕1周運(yùn)岐2

（1. 華南理工大學(xué) 機(jī)械與汽車工程學(xué)院，廣東 廣州 510640；2. 韶關(guān)市宏乾智能裝備科技有限公司，廣東 韶關(guān) 512029）

模內(nèi)熱壓定型是紙漿模塑餐具生產(chǎn)過(guò)程中的一種常用干燥方式，它是將成型后獲得的濕紙模胚在被模具擠壓與抽真空的條件下進(jìn)行加熱。加熱板作為熱壓定型機(jī)的熱源，其工作表面的溫度均勻性影響著制品的干燥質(zhì)量。針對(duì)紙漿模塑熱壓定型過(guò)程中加熱板的溫度不均勻性問(wèn)題，文中提出了一種聯(lián)合仿真與正交試驗(yàn)的優(yōu)化方法。首先，對(duì)加熱板的工作過(guò)程進(jìn)行分析，建立了加熱板的傳熱模型；然后，基于Fluent對(duì)加熱板進(jìn)行溫度場(chǎng)的數(shù)值模擬，根據(jù)溫度場(chǎng)分布結(jié)果將油路結(jié)構(gòu)中的高溫區(qū)域與低溫區(qū)域盡可能交錯(cuò)，設(shè)計(jì)了4種新的迷宮式油路結(jié)構(gòu)；最后，以油路結(jié)構(gòu)、油路平面高度、加熱板厚度和油路截面直徑設(shè)計(jì)了4因素4水平的正交試驗(yàn)，并進(jìn)行了極差分析與方差分析。結(jié)果表明：在實(shí)際干燥過(guò)程中，工作表面的最高溫度為224.47 ℃，最低溫度為209.92 ℃，溫度極差高達(dá)14.55 ℃，溫度標(biāo)準(zhǔn)差為3.01 ℃；加熱板厚度和油路直徑的大小對(duì)溫度極差的影響顯著，油路結(jié)構(gòu)對(duì)溫度標(biāo)準(zhǔn)差的影響顯著。基于以上分析改進(jìn)了加熱板的結(jié)構(gòu)，與原設(shè)計(jì)方案相比，加熱板工作表面的溫度極差降至7.27 ℃，溫度標(biāo)準(zhǔn)差降至1.09 ℃，保證了加熱板溫度的均勻性，提升了紙漿模塑產(chǎn)品的質(zhì)量。

紙漿模塑；加熱板；溫度分布；正交試驗(yàn)

紙漿模塑制品是一種以天然纖維或二次纖維為原料，經(jīng)過(guò)碎解制漿、成型、干燥、整型與切邊工序后制成的產(chǎn)品，廣泛應(yīng)用于包裝領(lǐng)域中［1］。在制品的熱壓定型過(guò)程中，模具通過(guò)與加熱板緊密貼合吸收熱量。因此，加熱板工作表面的溫度分布將直接影響制品的定型質(zhì)量。溫差過(guò)大，會(huì)使同一模具不同型腔內(nèi)的干燥速率不同，生產(chǎn)出的制品質(zhì)量不均。除此之外，溫度較高區(qū)域的制品收縮率增加、尺寸精度下降，溫度較低區(qū)域的制品由于干燥不完全而機(jī)械強(qiáng)度低［2］。因此，對(duì)現(xiàn)有加熱板結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化，提高加熱板工作面的溫度均勻性具有重要意義。

在紙漿模塑熱壓定型過(guò)程中，加熱板的加熱方式主要有油加熱、電加熱及蒸汽加熱3種類型。油加熱方式是通過(guò)導(dǎo)熱油在其內(nèi)部孔道回路中的流動(dòng)來(lái)傳遞熱量，相比于電加熱板，它具有升溫速率快、熱均勻性好的優(yōu)點(diǎn)［3］。近年來(lái)，國(guó)內(nèi)外學(xué)者對(duì)提高加熱板的均勻性進(jìn)行了大量研究?；瑥V軍等［2］利用Ansys對(duì)導(dǎo)熱油加熱板內(nèi)部孔道結(jié)構(gòu)的不同布置形式進(jìn)行溫度場(chǎng)分析并優(yōu)化，減小了工作面81%的溫度極差。徐京豫等［4］利用Fluent模擬了硫化機(jī)熱板溫度場(chǎng)分布及導(dǎo)熱油流速對(duì)整體溫差的影響規(guī)律，對(duì)熱板結(jié)構(gòu)及入口流速進(jìn)行了優(yōu)化，實(shí)驗(yàn)結(jié)果表明，優(yōu)化后的方案相比于原方案降低了33%的溫差。劉鳴宇等［5］針對(duì)膠帶硫化機(jī)熱板的溫度不均勻性問(wèn)題，利用正交試驗(yàn)設(shè)計(jì)的方法并使用Comsol模擬了各方案溫度場(chǎng)，得到了最佳設(shè)計(jì)方案并通過(guò)實(shí)驗(yàn)驗(yàn)證，最佳方案的表面溫差相較于原方案降低了8.2 ℃。Malygin等［6］建立了液壓機(jī)不同類型加熱板優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型，給出了感應(yīng)加熱板和電阻加熱板模型的求解結(jié)果，并與工業(yè)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較，平均偏差在1.5%內(nèi)。Wang等［7］利用Comsol對(duì)熱板的感應(yīng)加熱過(guò)程進(jìn)行了數(shù)值模擬，并采用主成分分析和正交試驗(yàn)相結(jié)合的方法對(duì)熱板的形狀參數(shù)進(jìn)行分析，得到了熱板目標(biāo)表面的最佳溫度分布。目前，國(guó)內(nèi)對(duì)電加熱板的溫度控制方面的研究居多，針對(duì)電熱棒的功率及尺寸排布的優(yōu)化已很成熟［8-9］，而對(duì)導(dǎo)熱油加熱板的結(jié)構(gòu)尺寸的優(yōu)化還很少見。

文中以某公司生產(chǎn)的加熱板為研究對(duì)象，通過(guò)Fluent對(duì)實(shí)際熱壓定型過(guò)程中的加熱板進(jìn)行了溫度場(chǎng)模擬，并以加熱板工作表面的溫度極差及溫度標(biāo)準(zhǔn)差最小為優(yōu)化目標(biāo)，設(shè)計(jì)正交試驗(yàn)優(yōu)化加熱板的油路結(jié)構(gòu)、油路平面高度、加熱板厚度及油路截面直徑等參數(shù)，以改善加熱板工作表面的溫度均勻性，提升產(chǎn)品質(zhì)量。

1　加熱板的溫度場(chǎng)模擬

1.1　加熱板物理模型的建立

如圖1及圖2所示，加熱板的三維（3D）尺寸（長(zhǎng)×寬×高）為950 mm×950 mm×80 mm、厚度=80 mm，油路平面高度=50 mm，加熱板上表面使用石棉板進(jìn)行隔熱。加熱板邊緣設(shè)置有安裝槽，便于與模具鎖合。側(cè)面油路進(jìn)出口設(shè)有與管路法蘭相連接的螺紋孔，機(jī)械臂通過(guò)加熱板中部的吸氣孔抽真空。油路結(jié)構(gòu)為單進(jìn)單出型，進(jìn)出口之間的油路對(duì)稱布置，各導(dǎo)熱通道之間的尺寸為1=122.5 mm，2=125.0 mm，3=105.0 mm，距邊緣距離=65.0 mm，油路截面直徑為25.0 mm。

圖1　加熱板的三維模型及尺寸

圖2　加熱板的油路結(jié)構(gòu)

1.2　數(shù)學(xué)模型

在持續(xù)工作過(guò)程中，導(dǎo)熱油在加熱板內(nèi)的換熱過(guò)程可視為三維無(wú)內(nèi)熱源常物性流體的穩(wěn)態(tài)傳熱，在直角坐標(biāo)系下的控制方程［10-11］為

質(zhì)量守恒方程

動(dòng)量守恒方程

能量守恒方程

式中：為流體速度矢量，為重力加速度矢量，為流體溫度矢量，均為三維矢量；為流體壓力；為流體動(dòng)力黏度；為固體材料導(dǎo)熱系數(shù)；為流體比熱容；為流體密度。

1.3　傳熱模型的簡(jiǎn)化

對(duì)加熱板傳熱過(guò)程進(jìn)行以下基本假設(shè)：

（1）導(dǎo)熱油流動(dòng)過(guò)程中的雷諾數(shù)較大，可視為完全發(fā)展的湍流流動(dòng)；

（2）工作過(guò)程中導(dǎo)熱油進(jìn)、出口溫差較小，視為不可壓縮常物性流體；

（3）忽略材料熱物性參數(shù)及邊界條件隨溫度的變化；

（4）加熱板溫度較低，忽略傳熱過(guò)程中的輻射影響及導(dǎo)熱油的黏性耗散生熱；

（5）周圍環(huán)境溫度恒為25 ℃。

1.4　網(wǎng)格劃分

為了便于網(wǎng)格劃分，對(duì)加熱板螺紋孔、安裝孔及排氣孔、小圓角進(jìn)行了簡(jiǎn)化處理。對(duì)流體區(qū)域采用結(jié)構(gòu)化網(wǎng)格劃分，對(duì)固體區(qū)域采用混合網(wǎng)格劃分。在流道拐角處進(jìn)行局部加密，流固交界面處劃分邊界層，局部網(wǎng)格劃分如圖3所示，總網(wǎng)格數(shù)量為440萬(wàn)。

圖3　加熱板局部網(wǎng)格劃分

1.5　邊界條件設(shè)置

基于上述建立的傳熱模型，選用烷基苯合成導(dǎo)熱油，導(dǎo)熱油初始溫度為240 ℃。加熱板材料為Al 6061，其余熱物性參數(shù)如表1所示［12］。流體與固體的交界面設(shè)置為耦合邊界條件。入口邊界條件設(shè)置為速度入口，油泵體積流量為16.798 m3/h，入口速率由式（4）求得，為4.753 m/s。出口邊界條件設(shè)為自由出流邊界條件。加熱板側(cè)面、石棉頂部表面為第3類邊界條件，與空氣的對(duì)流換熱系數(shù)取為8.45 W/（m2·℃）［3，13］；環(huán)境溫度為25 ℃，為第3類邊界條件。加熱板下表面為工作表面，對(duì)外輸出熱流。熱流密度根據(jù)實(shí)際工況設(shè)為-49 529.09 W/m2。

式中，V為油泵的體積流量，為加熱板個(gè)數(shù)，為加熱板進(jìn)口面積，為加熱板導(dǎo)熱通道管徑。

表1　模擬所用的材料熱物性參數(shù)

Table 1　Simulated thermophysical parameters of materials

名稱材料密度/（kg·m-3）比熱容/（J·（kg·K）-1）熱導(dǎo)率/（W·（m·K）-1）運(yùn)動(dòng)黏度/（10-4 Pa·s） 導(dǎo)熱油烷基苯DR 305715.02 7200.141.465 75 加熱板Al 60612 675.8996206.20

1.6　加熱板溫度場(chǎng)的求解

為了更好地評(píng)估溫度均勻性，除了工作表面的溫度極差外，引入加熱板工作表面的溫度標(biāo)準(zhǔn)差作為第二指標(biāo)，即

加熱板工作表面溫度場(chǎng)結(jié)果如圖4所示，加熱板工作表面的溫度呈現(xiàn)左高右低的趨勢(shì)，且溫度沿著油路的循行方向總體呈下降趨勢(shì)，但在油路的拐角處會(huì)有局部溫升。如圖5所示，導(dǎo)熱油的流動(dòng)在進(jìn)入拐角之前直線段部分為層流，經(jīng)過(guò)拐角處導(dǎo)熱油受壁面的黏性作用出現(xiàn)了二次流漩渦，破壞了流動(dòng)邊界層和熱邊界層的發(fā)展，從而加強(qiáng)了導(dǎo)熱油的對(duì)流換熱作用。由于邊緣處受壁面自然對(duì)流作用的影響較強(qiáng)，加熱板邊緣的溫度梯度高，溫度較低。云圖結(jié)果表明，在實(shí)際干燥過(guò)程中，工作表面的最高溫度為224.47 ℃，最低溫度為209.92 ℃，溫度極差高達(dá)14.55 ℃，溫度標(biāo)準(zhǔn)差為3.01 ℃。工作表面溫度分布不均勻，需要對(duì)加熱板進(jìn)行優(yōu)化。

圖4　加熱板工作表面溫度分布云圖

圖5　導(dǎo)熱油路彎曲部分流場(chǎng)圖

2　加熱板結(jié)構(gòu)參數(shù)單因素影響分析

2.1　加熱板油路結(jié)構(gòu)的優(yōu)化方案

原加熱板的拐角較少，導(dǎo)熱油在加熱板中的換熱不夠充分，兩平行相鄰豎直流道之間溫差較大，需要縮短回路間距。根據(jù)實(shí)際的工程需要，文中提出4種油路結(jié)構(gòu)的改進(jìn)方案，如圖6所示，4種油路結(jié)構(gòu)均為迷宮式結(jié)構(gòu)。流道距邊緣距離為70 mm，部分流道間距0為90 mm，其余流道間距均為70 mm。各油路結(jié)構(gòu)的溫度性能比較及溫度云圖分別如表2和圖7所示。

圖6　4種油路結(jié)構(gòu)優(yōu)化方案

表2　各油路結(jié)構(gòu)的溫度性能比較

Table 2　Comparison of temperature performance affected by oil circuit structure

油路結(jié)構(gòu)溫度/℃ 極差標(biāo)準(zhǔn)差平均值 原油路14.553.01218.21 油路111.971.71221.59 油路212.261.77221.35 油路312.221.35220.86 油路411.631.20221.34

圖7　4種油路結(jié)構(gòu)的溫度云圖

從表2及圖7可以看出，新設(shè)計(jì)的4種油路結(jié)構(gòu)工作表面的溫度極差和溫度標(biāo)準(zhǔn)差均低于原加熱板，平均溫度均高于原加熱板。油路4的溫度極差與溫度標(biāo)準(zhǔn)差均最低，分別為11.63 ℃和1.20 ℃。新油路結(jié)構(gòu)溫度極差相差不大，最高與最低僅相差0.63 ℃，但溫度標(biāo)準(zhǔn)差下降明顯，說(shuō)明在流量一定的情況下，僅改變油路結(jié)構(gòu)就能使加熱板的換熱更加充分，溫度分布更加均勻。迷宮式油路結(jié)構(gòu)的布置方式使油路中靠近入口的高溫導(dǎo)熱油與靠近出口的低溫導(dǎo)熱油相互交錯(cuò)，盡可能地平衡了高溫區(qū)域?qū)Φ蜏貐^(qū)域的影響，使整個(gè)板面的溫度分布更均勻［13］。

2.2　其他尺寸因素的影響

為了探究其他尺寸因素對(duì)溫度均勻性的影響，選擇均勻性最好的油路4作為分析的油路結(jié)構(gòu)［14］。不同因素對(duì)溫度均勻性的影響如圖8所示，可以看出，隨著加熱板厚度及油路直徑的增加，工作表面的溫度極差及溫度標(biāo)準(zhǔn)差都相應(yīng)地降低，但加熱板厚度增加到超過(guò)105 mm后，工作表面的溫度均勻性沒(méi)有得到明顯改善，油路直徑超過(guò)34 mm，工作表面的溫度均勻性開始變差，這是受加熱板材料的傳熱能力所限，隨著工作表面平均溫度的增加，在達(dá)到一定溫度后，低溫區(qū)域的升溫幅度比高溫區(qū)域要慢。油路平面在加熱板中部附近均勻性好，除了油路平面高度為50 mm的點(diǎn)外，溫度極差與溫度標(biāo)準(zhǔn)差隨著各因素的變化趨勢(shì)相同，說(shuō)明兩者是從不同角度衡量溫度均勻性。因此，油路平面高度應(yīng)靠近加熱板中部，厚度與直徑應(yīng)在一定范圍內(nèi)選較大值，如何選定這些參數(shù)的最優(yōu)組合，需要進(jìn)一步進(jìn)行正交試驗(yàn)。

3　加熱板結(jié)構(gòu)參數(shù)的協(xié)同優(yōu)化

3.1　正交試驗(yàn)設(shè)計(jì)

正交試驗(yàn)可以揭示參數(shù)的影響程度，以較少試驗(yàn)次數(shù)確定最佳參數(shù)組合［15］。根據(jù)單因素試驗(yàn)結(jié)果選取合適的水平，以溫度極差、溫度標(biāo)準(zhǔn)差為考察指標(biāo)，不考慮交互作用［16］，根據(jù)標(biāo)準(zhǔn)正交表16（45），安排4因素4水平正交試驗(yàn)，因素水平安排如表3所示，試驗(yàn)方案及結(jié)果如表4所示。

表3　加熱板參數(shù)因素水平表

Table 3　Parameter factor level table of heating plate

水平A油路結(jié)構(gòu)B油路平面高度/mmC加熱板厚度/mmD油路截面直徑/mmE空白列 1油路12980251 2油路23690282 3油路343100313 4油路450110344

表4　正交試驗(yàn)方案及結(jié)果

Table 4　Scheme and results of orthogonal test

序號(hào)ABCDE溫度標(biāo)準(zhǔn)差/℃溫度極差/℃ 1111111.6211.16 2122221.319.38 3133331.178.47 4144441.087.52 5212341.408.69 6221431.4410.22 7234121.438.46 8243211.439.42 9313421.138.17 10324311.047.45 11331241.2410.73 12342131.1910.18 13414231.238.41 14423141.007.91 15432411.158.72 16441321.1710.78

3.2　極差分析與方差分析

分別對(duì)表4中不同參數(shù)水平組合條件下工作表面的最高溫差和溫度標(biāo)準(zhǔn)差數(shù)據(jù)進(jìn)行極差與方差分析，結(jié)果如表5及表6所示，其中k為因素的各水平均值。

表5　極差分析結(jié)果

Table 5　Results of range analysis

性能指標(biāo)水平ABCD 溫度極差k19.1329.10710.7229.428 k29.1988.7409.2409.483 k39.1339.0978.4968.847 k48.9549.4737.9608.659 極差0.2440.7332.7620.825 極差排序4312 優(yōu)組合A4B2C4D4 溫度標(biāo)準(zhǔn)差k11.2971.3481.3691.310 k21.4261.1981.2621.304 k31.1501.2481.1841.196 k41.1371.2161.1941.201 極差0.2890.1500.1850.114 極差排序1324 優(yōu)組合A4B2C3D3

表6　方差分析結(jié)果1）

Table 6　Results of variance analysis

性能指標(biāo)來(lái)源自由度f(wàn)離均差平方和均方F值P值顯著性 溫度極差A(yù)30.1310.0440.150.922不顯著 B31.0750.3581.250.431不顯著 C317.2575.75219.990.017極顯著 D32.0550.6852.380.247不顯著 殘差誤差30.8630.288 合計(jì)1521.381 溫度標(biāo)準(zhǔn)差A(yù)30.2230.0746.750.076顯著 B30.0540.0181.620.351不顯著 C30.0870.0292.630.224不顯著 D30.0470.0161.430.388不顯著 殘差誤差30.0330.011 合計(jì)150.444

1）=3，01（3，3）=5.36，0.05（3，3）=9.28。

極差分析結(jié)果表明：各因素對(duì)工作表面溫度極差的影響順序?yàn)镃>D>B>A，優(yōu)組合為A4B2C4D4；各因素對(duì)溫度標(biāo)準(zhǔn)差的影響順序?yàn)锳>C>B>D，優(yōu)組合為A4B2C3D3，正交試驗(yàn)的16組數(shù)據(jù)中溫度極差的最小值為7.45 ℃，溫度標(biāo)準(zhǔn)差的最小值為1.00 ℃。兩種性能指標(biāo)的優(yōu)組合的模擬值在正交試驗(yàn)表中均無(wú)出現(xiàn)，需要進(jìn)行補(bǔ)充試驗(yàn)。

方差分析結(jié)果表明：在置信度為95%的情況下，因素C對(duì)溫度極差的影響顯著；在置信度為90%的情況下，因素A對(duì)溫度極差的影響顯著。根據(jù)值大小，得出各因素對(duì)溫度極差的顯著性排序?yàn)镃>D>B>A，各因素對(duì)溫度標(biāo)準(zhǔn)差的顯著性排序?yàn)锳>C>B>D。溫度極差的空白列均方值略大于因素A，為次小值，而溫度標(biāo)準(zhǔn)差的空白列均方值為最小值，可認(rèn)為試驗(yàn)中各因素之間的交互作用不顯著［15］。

3.3　補(bǔ)充實(shí)驗(yàn)

對(duì)兩種性能指標(biāo)的參數(shù)優(yōu)組合A4B2C4D4與A4B2C3D3進(jìn)行補(bǔ)充實(shí)驗(yàn)，得到溫度極差的模擬值分別為7.27 ℃和7.39 ℃，溫度標(biāo)準(zhǔn)差的模擬值為1.09 ℃和1.06 ℃，由于因素C和D對(duì)溫度標(biāo)準(zhǔn)差的影響不顯著，故最優(yōu)組合可取A4B2C4D4，溫度標(biāo)準(zhǔn)差模擬值僅高了0.03 ℃。與優(yōu)化前的加熱板相比，最優(yōu)結(jié)構(gòu)的溫度極差降低了50.0%，溫度標(biāo)準(zhǔn)差降低了63.8%。最優(yōu)結(jié)構(gòu)的工作面溫度云圖如圖9所示，可以看出，最優(yōu)結(jié)構(gòu)的高溫部分均勻分布于加熱板中部，達(dá)到較優(yōu)水平，模具應(yīng)在此部分排布較密，在邊緣部分排布較疏。

圖9　最優(yōu)結(jié)構(gòu)的工作面溫度云圖

4　結(jié)論

文中通過(guò)Fluent對(duì)加熱板溫度場(chǎng)進(jìn)行了仿真，采用正交試驗(yàn)設(shè)計(jì)法對(duì)影響加熱板的參數(shù)組合進(jìn)行了極差與方差分析，并根據(jù)分析結(jié)果進(jìn)一步選出最優(yōu)的參數(shù)組合，得出以下結(jié)論：

（1）影響加熱板表面溫度均勻性最顯著的因素是加熱板厚度，其次是油路結(jié)構(gòu)，因此在對(duì)加熱板尺寸進(jìn)行結(jié)構(gòu)設(shè)計(jì)時(shí)，應(yīng)保證加熱板的厚度足夠，高溫區(qū)域與低溫區(qū)域盡可能交錯(cuò)；

（2）優(yōu)化后加熱板表面的溫度極差為7.27 ℃，溫度標(biāo)準(zhǔn)差為1.09 ℃，較優(yōu)化前溫度極差降低了50.0%，溫度標(biāo)準(zhǔn)差降低了63.8%，加熱板工作表面的均勻性得到顯著提高；

（3）可以根據(jù)加熱板工作表面的溫度場(chǎng)的仿真結(jié)果，合理排布模具中不同制品型腔的位置分布，高溫區(qū)域制品排列應(yīng)較密，低溫區(qū)域制品排列應(yīng)較疏。

［1］ 黃俊彥．紙漿模塑生產(chǎn)實(shí)用技術(shù)［M］．2版．北京：文化發(fā)展出版社，2008：22-34.

［2］ 滑廣軍，廖澤順，謝勇，等．基于Ansys的紙漿模塑模具加熱板溫度場(chǎng)分析及優(yōu)化［J］．包裝工程，2016，37（1）：78-82.

HUA Guang-jun，LIAO Ze-shun，XIE Yong，et al．Analysis and optimization of temperature field in pulp molding mold heating plate based on Ansys［J］．Packaging Engineering，2016，37（1）：78-82.

［3］ 楊彬，孫寶壽，張雅靈，等．基于ANSYS Workbench熱板電加熱與油加熱模擬及對(duì)比分析［J］．寧波大學(xué)學(xué)報(bào)（理工版），2017，30（5）：13-18.

YANG Bin，SUN Baoshou，ZHANG Yaling，et al．Comparative analysis：simulation of electrical heating and oil heating of hot plate using ANSYS workbench ［J］．Journal of Ningbo University （Natural Science & Engineering Edition），2017，30（5）：13-18.

［4］ 徐京豫，孫寶壽．油壓機(jī)熱板溫度場(chǎng)均勻性研究［J］．機(jī)械設(shè)計(jì)與研究，2019，35（1）：200-204.

XU Jingyu，SUN Baoshou．Analysis of surface temperature uniformity of heat conducting oil heating plate ［J］．Machine Design and Research，2019，35（1）：200-204.

［5］ 劉鳴宇，李軍霞，張弘玉，等．基于正交試驗(yàn)的膠帶硫化機(jī)熱板溫度場(chǎng)研究［J］．煤炭工程，2021，53（11）：152-157.

LIU Ming-yu，LI Jun-xia，ZHANG Hong-yu，et al．Study on heating plate temperature field of conveyor belt vulcanizer based on orthogonal ［J］．Coal Engineering，2021，53（11）：152-157.

［6］ MALYGIN E，KRASNYANSKIY M，MOKROZUB V，et al．Design optimization of heating plates for hydraulic presses ［J］．MATEC Web of Conferences，2017，129：01008/1-5.

［7］ WANG J，CHEN X，YU L，et al．A novel method improving the temperature uniformity of hot-plate under induction heating ［J］．Proceedings of the Institution of Mechanical Engineers Part C：Journal of Mechanical Engineering Science，2021，235（1）：190-201.

［8］ 張晨曦，王晉強(qiáng)，丁永生．基于Kriging模型的熱板表面溫度均勻性控制［J］．電子工藝技術(shù)，2019，40（4）：227-229，240.

ZHANG Chenxi，WANG Jinqiang，DING Yongsheng．Temperature uniformity control of hot plate surface based on Kriging model ［J］．Electronics Process Technology，2019，40（4）：227-229，240.

［9］ 陳惠賢，馮國(guó)慶，王保增，等．基于溫度場(chǎng)的層壓機(jī)電加熱板優(yōu)化設(shè)計(jì)［J］．新技術(shù)新工藝，2013（9）：63-65.

CHEN Huixian，F(xiàn)ENG Guoqing，WANG Baozeng，et al．Optimal design of electric heating plate of laminating machine based on temperature field ［J］．New Technology & New Process，2013（9）：63-65.

［10］ 楊世銘，陶文銓．傳熱學(xué)［M］．4版．北京：高等教育出版社，2006：205-206.

［11］ 殷術(shù)貴，郭偉科，黃棟，等．流延薄膜傳熱特性及冷卻水量設(shè)計(jì)的仿真研究［J］．華南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2021，49（12）：23-34.

YIN Shugui，GUO Weike，HUANG Dong，et al．Simulation study on heat transfer characteristics and cooling water design for the casting film ［J］．Journal of South China University of Technology（Natural Science Edition），2021，49（12）：23-34.

［12］ MILLS K C．Recommended values of thermophysical properties for selected commercial alloys［M］．Sawston：Woodhead Publishing，2002：64-67.

［13］ 王野平，李翠娟，周韓．基于FLUENT的熱壓板通道布置設(shè)計(jì)優(yōu)化［J］．機(jī)械設(shè)計(jì)，2012，29（3）：50-52.

WANG Ye-ping，LI Cui-juan，ZHOU Han．Optimization of hot platen channel layout design based on FLUENT ［J］．Journal of Machine Design，2012，29（3）：50-52.

［14］ 何偉，黃春躍，梁穎．基于正交設(shè)計(jì)的LTCC基板微流道散熱性能研究［J］．機(jī)械強(qiáng)度，2019，41（5）：1223-1227.

HE Wei，HUANG Chunyue，LIANG Ying．Research on heat dissipation of LTCC substrate microchannel based on orthogonal design ［J］．Journal of Mechanical Strength，2019，41（5）：1223-1227.

［15］ 任露泉.試驗(yàn)優(yōu)化設(shè)計(jì)與分析［M］．2版．北京：高等教育出版社，2003：10-25.

［16］ 謝方偉，徐純潔，萬(wàn)快弟，等．CDC減振器用先導(dǎo)閥參數(shù)協(xié)同優(yōu)化設(shè)計(jì)［J］．華南理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2021，49（2）：131-139.

XIE Fangwei，XU Chunjie，WAN Kuaidi，et al．Cooperative optimization design of pilot valve parameters for CDC shock absorber ［J］．Journal of South China University of Technology （Natural Science Edition），2021，49（2）：131-139.

Analysis of Temperature Uniformity of Pulp Molding Hot Pressing Heating Plate Based on Fluent Simulation

112

（1. School of Mechanical and Automotive Engineering，South China University of Technology，Guangzhou 510640，Guangdong，China；2. Shaoguan Hongqian Intelligent Equipment Technology Co.，Ltd.，Shaoguan 512029，Guangdong，China）

In-mold hot-pressing is a common drying method in the production of pulp moulded tableware. The wet paper mold embryo obtained after moulding is heated under the condition of extrusion and vacuum. Heating plate is the heat source of hot-pressing machine, and the temperature uniformity of its working surface affects the drying quality of products. Aiming at the temperature non-uniformity of the heating plate in the process of pulp molding hot pressing, this paper proposed an optimization method which combined simulation and orthogonal experiment. Firstly, the working process of heating plate was analyzed and the heat transfer model of heating plate was established. Then, based on Fluent, the temperature field was simulated numerically. According to the results of temperature field distribution, the high-temperature area and low-temperature area in the oil circuit structure were staggered as far as possible, and four new labyrinth oil circuit structures were designed. Finally, an orthogonal test with 4 factors and 4 levels was designed based on the structure of oil circuit, the plane height of oil circuit, the thickness of heating plate and the section diameter of oil circuit, and the range analysis and variance analysis were carried out. The results show that during the actual drying process, the maximum temperature and minimum temperature of the working surface are 224.47 ℃ and 209.92 ℃, respectively, and the temperature range is as high as 14.55 ℃ and the temperature standard deviation is 3.01 ℃. The thickness of heating plate and the diameter of oil passage have a significant influence on the temperature difference, and the structure of oil passage has a significant influence on the temperature standard deviation. Based on the above analysis, the structure of heating plate was improved. Compared with the original design, the temperature range of the working surface of heating plate has been reduced to 7.27 ℃ and the temperature standard deviation has been reduced to 1.09 ℃, which ensures the uniformity of temperature of heating plate and improves the quality of pulp molded products.

molded pulp；heating plate；temperature distribution；orthogonal test

Supported by the Science and Technology Planning Project of Guangdong Province （2021440002002290）

10.12141/j.issn.1000-565X.220363

2022?06?10

廣東省科技計(jì)劃項(xiàng)目（2021440002002290）

吳上生（1963-），男，博士，教授，主要從事精密機(jī)械設(shè)計(jì)理論、自動(dòng)化生產(chǎn)線設(shè)計(jì)、機(jī)電一體化技術(shù)及其應(yīng)用研究。E-mail：shshwu@scut.edu.cn

TH122；TK124

1000-565X（2023）05-0122-08