袁啟龍 李春霞 許 力 李 言 嚴(yán)和平

（1. 西安理工大學(xué)機(jī)械與精密儀器工程學(xué)院，陜西西安，710048；2. 陜西理工機(jī)電科技有限公司，陜西西安，710018）

造紙工業(yè)作為我國(guó)基礎(chǔ)原材料產(chǎn)業(yè)之一，對(duì)我國(guó)國(guó)民經(jīng)濟(jì)的發(fā)展起著重要作用，與國(guó)民經(jīng)濟(jì)有序發(fā)展和社會(huì)文明建設(shè)息息相關(guān)。造紙技術(shù)的發(fā)展水平在一定程度上體現(xiàn)了國(guó)家經(jīng)濟(jì)的發(fā)展程度，是經(jīng)濟(jì)發(fā)展中具有循環(huán)經(jīng)濟(jì)特征的重要基礎(chǔ)原材料產(chǎn)業(yè)和重要經(jīng)濟(jì)增長(zhǎng)點(diǎn)［1-2］。濕法造紙生產(chǎn)過(guò)程中需要大量水資源，容易造成水污染和空氣污染，且其生產(chǎn)的紙制產(chǎn)品難以滿(mǎn)足一些應(yīng)用場(chǎng)合的特殊需求［3-4］。干法造紙技術(shù)具有能耗低、出紙效率高、無(wú)水污染、紙漿纖維利用率高、投資成本低等優(yōu)點(diǎn)，被廣泛應(yīng)用于各種生活用紙、衛(wèi)生巾、醫(yī)療床單、靠墊、桌布、餐巾紙等生產(chǎn)領(lǐng)域，正逐步成為造紙行業(yè)的主流技術(shù)［5-6］。負(fù)壓箱作為成型布料器的重要組成部件，其內(nèi)部流場(chǎng)分布的均勻性影響干法造紙過(guò)程中紙幅成型的均勻度。本課題對(duì)干法造紙?jiān)O(shè)備（紙幅寬度1.6 m）中的負(fù)壓箱進(jìn)行流體動(dòng)力學(xué)仿真計(jì)算，研究風(fēng)道結(jié)構(gòu)及排風(fēng)量對(duì)箱體內(nèi)流場(chǎng)分布均勻性的影響，對(duì)負(fù)壓箱結(jié)構(gòu)進(jìn)行優(yōu)化，以期為生產(chǎn)高質(zhì)量的無(wú)塵紙?zhí)峁├碚撘罁?jù)。

1 干法造紙技術(shù)基本原理

干法造紙技術(shù)依據(jù)空氣動(dòng)力學(xué)基本原理，利用多風(fēng)道氣流篩粉成型技術(shù)，將粉碎后的紙漿纖維通過(guò)抽吸氣流在成型網(wǎng)上形成所要求的層狀纖維制品［7］。目前，干法造紙成型箱可分為圓筒篩分式結(jié)構(gòu)和攪拌過(guò)濾式結(jié)構(gòu)［8］，其中圓筒篩分式結(jié)構(gòu)具有結(jié)構(gòu)簡(jiǎn)單、成型速度快、工作效率高等特點(diǎn)。基于圓筒篩分式成型布料器的干法造紙?jiān)砣鐖D1 所示，其基本原理是：纖維原料和空氣混合氣流由送風(fēng)風(fēng)機(jī)分別從兩個(gè)旋轉(zhuǎn)的圓筒一端送入，兩個(gè)送料入口分別設(shè)置在箱體兩側(cè)，圓筒壁布滿(mǎn)一定尺寸的小孔，圓筒內(nèi)部裝有打散輥，經(jīng)過(guò)圓筒與打散輥相反的轉(zhuǎn)動(dòng)作用，將纖維原料打散分離，使符合沉降標(biāo)準(zhǔn)的纖維顆粒隨氣流穿過(guò)圓筒的小孔進(jìn)入成型箱體內(nèi)部并呈漂浮狀態(tài)，在上部送風(fēng)與負(fù)壓箱抽吸的共同作用下，纖維顆粒沉降于移動(dòng)的成型網(wǎng)上形成無(wú)強(qiáng)度的纖維層，隨著成型網(wǎng)的移動(dòng)可實(shí)現(xiàn)不間斷布料，經(jīng)過(guò)后續(xù)熱壓緊、噴膠、烘干膨化等工藝的處理，最終形成無(wú)塵原紙［9］。

圖1 圓筒篩分式成型布料器干法造紙?jiān)韴DFig.1 Schematic diagram of air-laid papermaking for cylindrical sieve forming distributor

圓筒篩分式成型布料器結(jié)構(gòu)如圖2所示，其主要由成型箱（包括圓筒篩、打散輥、托輥體、壓輪、機(jī)架等）、成型網(wǎng)和負(fù)壓箱等組成［10-11］。負(fù)壓箱是干法造紙成型布料過(guò)程中的關(guān)鍵部件，通過(guò)負(fù)壓風(fēng)機(jī)的抽吸作用將纖維物料吸附于成型網(wǎng)上，該結(jié)構(gòu)的流場(chǎng)特性決定纖維物料沉降的有序性及紙幅成型的均勻度。

圖2 圓筒篩分式成型布料器結(jié)構(gòu)Fig.2 Structure of cylindrical sieve forming distributor

2 負(fù)壓箱計(jì)算模型的建立

2.1 計(jì)算域的建立與網(wǎng)格劃分

根據(jù)不可壓縮流體連續(xù)性方程可知，當(dāng)排風(fēng)量恒定時(shí)，與單風(fēng)道結(jié)構(gòu)相比，多風(fēng)道結(jié)構(gòu)負(fù)壓箱內(nèi)不易形成渦流，負(fù)壓箱內(nèi)流場(chǎng)更均勻，即風(fēng)道越多時(shí)，負(fù)壓箱內(nèi)流場(chǎng)分布越均勻。負(fù)壓箱排風(fēng)口與負(fù)壓風(fēng)機(jī)抽風(fēng)管道相連，為了使抽風(fēng)管道與排風(fēng)口連接方便，通常將排風(fēng)口形狀設(shè)計(jì)為矩形，封閉式五風(fēng)道負(fù)壓箱結(jié)構(gòu)如圖3所示。

圖3 封閉式五風(fēng)道負(fù)壓箱結(jié)構(gòu)Fig.3 Structure of enclosed five-duct negative pressure box

考慮工程設(shè)計(jì)制造的實(shí)際需求及仿真計(jì)算的局限性，為確保負(fù)壓箱風(fēng)道內(nèi)流量均勻分配并形成均勻穩(wěn)定的流場(chǎng)，設(shè)計(jì)3種負(fù)壓箱風(fēng)道結(jié)構(gòu)。依據(jù)流場(chǎng)均勻性進(jìn)行負(fù)壓箱結(jié)構(gòu)的優(yōu)化選擇，使纖維能在成型網(wǎng)上均勻沉降。圖4 為3 種負(fù)壓箱Z截面示意圖；其中，圖4（a）為單風(fēng)道結(jié)構(gòu)，即在負(fù)壓箱內(nèi)部設(shè)置單一隔板，將箱體內(nèi)部分割為左右單風(fēng)道；圖4（b）為封閉式五風(fēng)道結(jié)構(gòu)，在單風(fēng)道結(jié)構(gòu)的基礎(chǔ)上，左右兩部分設(shè)置等間距折彎隔板（厚度2 mm），將負(fù)壓箱內(nèi)部分割為左右各有5個(gè)風(fēng)道的多風(fēng)道結(jié)構(gòu)；圖4（c）為敞開(kāi)式五風(fēng)道結(jié)構(gòu)，在封閉式五風(fēng)道結(jié)構(gòu)的基礎(chǔ)上不封閉左側(cè)五風(fēng)道隔板末端，使負(fù)壓箱入口氣流快速進(jìn)入下層并從排風(fēng)口排出，避免氣流在負(fù)壓箱上層分流造成流場(chǎng)分布不均勻。

圖4 負(fù)壓箱Z截面示意圖Fig.4 Schematic diagram of Z-section of negative pressure box

為方便表達(dá)仿真結(jié)果，將五風(fēng)道負(fù)壓箱內(nèi)風(fēng)道按照?qǐng)D5所示進(jìn)行命名并標(biāo)明坐標(biāo)系，左右兩部分從外到內(nèi)依次為left1～left5和right1～right5。

圖5 五風(fēng)道負(fù)壓箱Y截面圖Fig.5 Y-section diagram of five-duct negative pressure box

采用ANSYS Workbench 軟件中的ICEM 模塊對(duì)負(fù)壓箱流體域模型進(jìn)行網(wǎng)格劃分，考慮到計(jì)算域模型的體積較大及內(nèi)部隔板布置的復(fù)雜性及不規(guī)律性，為兼顧計(jì)算精度和運(yùn)算效率，采用四面體網(wǎng)格單元對(duì)負(fù)壓箱進(jìn)行離散化處理。在負(fù)壓箱的排風(fēng)口與進(jìn)風(fēng)口處采用較密的網(wǎng)格單元，其他區(qū)域采用相對(duì)稀疏的網(wǎng)格單元，共產(chǎn)生網(wǎng)格單元4290669 個(gè)，節(jié)點(diǎn)總數(shù)802297個(gè)，負(fù)壓箱網(wǎng)格劃分模型如圖6所示。

圖6 負(fù)壓箱網(wǎng)格劃分模型Fig.6 Grid division model of negative pressure box

2.2 邊界條件的設(shè)置與求解設(shè)置

負(fù)壓箱工作時(shí)采用抽出式排風(fēng)方式，負(fù)壓風(fēng)機(jī)在排風(fēng)口抽吸作用下，使氣流由箱體頂部進(jìn)風(fēng)口流入并流經(jīng)內(nèi)部多個(gè)風(fēng)道，最后由排風(fēng)口流出。將劃分好網(wǎng)格的負(fù)壓箱離散模型導(dǎo)入Fluent 中進(jìn)行流體仿真計(jì)算，計(jì)算域進(jìn)口采用自由出入模式，出口采用速度出口，初始出口速度為7.94 m/s（排風(fēng)量20000 m3/h），垂直于出口面；負(fù)壓箱壁面設(shè)置為無(wú)滑移條件；流體介質(zhì)為空氣，其密度為1.225 kg/m3。

利用湍流模型進(jìn)行計(jì)算時(shí)，需要給定入口邊界條件中的湍流參數(shù)，通過(guò)給定湍流強(qiáng)度（I）、水力直徑（DH）、湍流動(dòng)能（k）和湍流耗散率（ε）等來(lái)定義流場(chǎng)邊界上的湍流［12-13］。根據(jù)流體力學(xué)的經(jīng)典理論，湍流強(qiáng)度（I）通過(guò)式（1）計(jì)算。

湍流動(dòng)能（k）和湍流耗散率（ε）可由式（2）和式（3）計(jì)算得到。

式中，Cμ=0.09，ρ為空氣密度。

因此，湍流動(dòng)能（k）和湍流耗散率（ε）取值范圍 分 別 為 ： 7.5×10-4～1.9×10-2m2/s2和 1.8×10-4～2.3×10-2m2/s2。

在初始情況下，求解采用基于壓力的SIMPLE 非耦合隱式算法［14］，梯度項(xiàng)差分方法采用Least Squares Cell Based模式，壓力項(xiàng)采用Second Order格式，其余控制項(xiàng)采用默認(rèn)設(shè)置。設(shè)置求解控制的松弛因子為0.3，殘差收斂設(shè)置為0.0001，迭代3000 步后計(jì)算收斂。

3 負(fù)壓箱流場(chǎng)特性分析

3.1 單風(fēng)道負(fù)壓箱的流場(chǎng)均勻性分析

負(fù)壓箱的壓力和速度流場(chǎng)分布沿著Y、Z方向變化有所不同，Y方向的流場(chǎng)可明確反映各風(fēng)道的流場(chǎng)分布特點(diǎn)。Z方向的流場(chǎng)分布與紙幅的成型有密切關(guān)系。為便于分析，在成型布料器Z方向上分別建立Z=142、590、900 mm 和Y=-800、2、800 mm 等代表性截面作為研究目標(biāo)面。圖7為單風(fēng)道負(fù)壓箱內(nèi)流場(chǎng)仿真得到的不同截面速度分布云圖。從圖7（a1）～圖7（a3）可以看出，Z截面的速度場(chǎng)整體分布較均勻，且速度較小。從圖7（b1）～圖7（b3）可以看出，Y截面速度場(chǎng)分布出現(xiàn)瞬變性，風(fēng)速由進(jìn)風(fēng)口到排風(fēng)口逐漸增大，風(fēng)速最大為13 m/s，排風(fēng)口區(qū)域有效面積逐漸減小，加速集聚于排風(fēng)口，導(dǎo)致風(fēng)速增大。在此截面處的左右風(fēng)道產(chǎn)生明顯的大尺度渦流，這是由于風(fēng)道內(nèi)空氣可流動(dòng)面積較大，且流動(dòng)只受到壁面的阻礙，因而在兩側(cè)壁面之間出現(xiàn)大面積的渦流狀態(tài)，嚴(yán)重影響風(fēng)道速度均勻性，不利于纖維在成型網(wǎng)上的均勻成型。

圖7 單風(fēng)道負(fù)壓箱內(nèi)的速度分布云圖Fig.7 Cloud diagram of velocity distribution in single-air-channel negative pressure box

負(fù)壓箱結(jié)構(gòu)Y=-800 mm、Z=590 mm 截面均便于表達(dá)各風(fēng)道位置，且便于計(jì)算各風(fēng)道流量，因此，本課題選取Y=-800 mm截面計(jì)算各個(gè)風(fēng)道的流量。利用CFD后處理軟件求得單風(fēng)道負(fù)壓箱內(nèi)部Y=-800 mm截面的質(zhì)量流量與體積流量，結(jié)果如表1 所示。從表1可以看出，左右兩個(gè)風(fēng)道的質(zhì)量流量差異較小，體積流量比為1∶1.06，說(shuō)明兩風(fēng)道風(fēng)速差異較小，這是因?yàn)閮蓚€(gè)風(fēng)道的有效面積一致。

表1 單風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道流量Table 1 Air channel flow of single-air-channel negative pressure box at cross section Y=-800 mm

3.2 封閉式五風(fēng)道負(fù)壓箱的流場(chǎng)均勻性分析

圖8為封閉式五風(fēng)道負(fù)壓箱內(nèi)流場(chǎng)仿真得到的不同截面的速度分布云圖。由圖8（a1）～圖8（a3）可知，Z截面各風(fēng)道速度場(chǎng)分布比較均勻，但折彎隔板的拐角處速度較大，靠近隔板處速度減小，這是由空氣在兩側(cè)隔板之間發(fā)生接觸碰撞所致。由圖8（b1）～圖8（b3）可知，Y=-800 mm截面處氣流進(jìn)入負(fù)壓箱后繞過(guò)兩側(cè)的斜隔板沿著左右各風(fēng)道流動(dòng)，右側(cè)第五風(fēng)道出現(xiàn)較明顯的渦旋狀態(tài)，而左側(cè)風(fēng)道內(nèi)的流場(chǎng)分布相對(duì)較均勻。

圖8 封閉式五風(fēng)道負(fù)壓箱內(nèi)的速度分布云圖Fig.8 Cloud diagram of velocity distribution in enclosed five-duct negative pressure box

封閉式五風(fēng)道負(fù)壓箱Y=-800 mm截面的各風(fēng)道流量如表2 所示。由表2 可知，左側(cè)各風(fēng)道質(zhì)量流量與體積流量均相差不大、較均勻；右側(cè)各風(fēng)道質(zhì)量流量和體積流量波動(dòng)范圍較大，最小體積流量與最小質(zhì)量流量均出現(xiàn)在right1 風(fēng)道，與其他各風(fēng)道差異較大。左右風(fēng)道總體積流量比為1.117∶1，表明左右風(fēng)道流量分配不均衡，風(fēng)道速度不均勻。

圖9 為利用表2 中的各風(fēng)道體積流量計(jì)算獲得的Y=-800 mm 截面的風(fēng)道平均速度變化曲線(xiàn)圖。從圖9可以看出，左側(cè)與右側(cè)從外到內(nèi)的風(fēng)道速度都呈現(xiàn)增大趨勢(shì)，左側(cè)風(fēng)道速度變化范圍比右側(cè)的小，除了右側(cè)第一風(fēng)道速度略小于左側(cè)第一風(fēng)道速度外，右側(cè)其他各風(fēng)道速度均大于左側(cè)各風(fēng)道速度。這是由于考慮實(shí)際安裝空間限制和成本，左右風(fēng)道抽風(fēng)口設(shè)置在負(fù)壓箱同側(cè)，使得右側(cè)風(fēng)道與排風(fēng)口距離較近，沿程阻力小，因而右側(cè)風(fēng)道速度大于左側(cè)風(fēng)道速度。

圖9 封閉式五風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道平均速度Fig.9 Air channel average velocity of enclosed five-duct negative pressure box at cross section Y=-800 mm

表2 封閉式五風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道流量Table 2 Air channel flow of enclosed five-duct negative pressure box at cross section Y=-800 mm

圖10 為Y=-800 mm 截面的左右風(fēng)道平均壓力變化曲線(xiàn)圖。從圖10 可以看出，除左側(cè)第一風(fēng)道的壓力小于右側(cè)第一風(fēng)道壓力外，左側(cè)其他風(fēng)道壓力均大于右側(cè)各風(fēng)道壓力，且左側(cè)各風(fēng)道的壓力變化幅度較小，基本都維持在-80 Pa左右，表明左側(cè)風(fēng)道的壓力分布較均勻。right4風(fēng)道的壓力最小，為-297 Pa，而right1風(fēng)道的壓力基本接近正壓，壓力出現(xiàn)較大差異，表明右側(cè)風(fēng)道壓場(chǎng)不均勻，不利于纖維在成型網(wǎng)上的均勻著落。

圖10 封閉式五風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道平均壓力Fig.10 Air channel average pressure of enclosed five-duct negative pressure box at cross section Y=-800 mm

由以上分析可知，沿各風(fēng)道順序，負(fù)壓箱內(nèi)流場(chǎng)壓力變化與速度變化的整體趨勢(shì)相反。但風(fēng)道速度變化較小時(shí)壓力變化也相對(duì)穩(wěn)定，速度變化較大時(shí)的風(fēng)道壓力波動(dòng)也較大。因此可知，封閉式五風(fēng)道負(fù)壓箱流場(chǎng)right1 風(fēng)道速度場(chǎng)和負(fù)壓場(chǎng)的均勻性較差，應(yīng)在此區(qū)域采取措施以適當(dāng)降低壓力，增大速度，從而達(dá)到排風(fēng)均勻性的要求。

3.3 敞開(kāi)式五風(fēng)道負(fù)壓箱的流場(chǎng)均勻性分析

圖11 為敞開(kāi)式五風(fēng)道負(fù)壓箱內(nèi)流場(chǎng)仿真獲得的不同截面的速度分布云圖。從圖11（a1）～圖11（a3）可以看出，沿著+Z方向（Z=142 mm～Z=900 mm），Z截面各風(fēng)道速度場(chǎng)分布越均勻。Z=590 mm 截面各風(fēng)道氣流分布較均勻，在隔板折彎處速度較大，這是由空氣在兩側(cè)隔板之間發(fā)生接觸碰撞所致。對(duì)比封閉式五風(fēng)道負(fù)壓箱內(nèi)的速度分布云圖（見(jiàn)圖8（a1））可知，敞開(kāi)式五風(fēng)道負(fù)壓箱的左側(cè)風(fēng)道流場(chǎng)分布更加均勻。從圖11（b1）可以看出，Y=-800 mm截面處的左側(cè)風(fēng)道內(nèi)的流場(chǎng)分布相對(duì)右側(cè)風(fēng)道更均勻，速度變化梯度較小。結(jié)合圖11（b1）～圖11（b3）可知，負(fù)壓箱上層速度分布整體小于下層，且在負(fù)壓箱入口處速度分布更均勻。

圖11 敞開(kāi)式五風(fēng)道負(fù)壓箱內(nèi)的速度分布云圖Fig.11 Cloud diagram of velocity distribution in open five-duct negative pressure box

敞開(kāi)式五風(fēng)道負(fù)壓箱Y=-800 mm截面的各風(fēng)道流量如表3 所示。由表3 可知，左右風(fēng)道總體積流量比為1.08∶1，right1 風(fēng)道與right5 風(fēng)道的體積流量相差800 m3/h；相較于封閉式五風(fēng)道負(fù)壓箱，敞開(kāi)式五風(fēng)道負(fù)壓箱右側(cè)各風(fēng)道體積流量的差值有所減?。ǔ齬ight1風(fēng)道外），這是因?yàn)槌ㄩ_(kāi)式五風(fēng)道負(fù)壓箱的左右風(fēng)道相通后，平衡了左右風(fēng)道體積流量的分配比。

表3 敞開(kāi)式五風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道流量Table 3 Air channel flow of open five-duct negative pressure box at cross section Y=-800 mm

圖12 為敞開(kāi)式五風(fēng)道負(fù)壓箱Y=-800 mm 截面的左右各風(fēng)道平均速度變化曲線(xiàn)圖。由圖12 可知，左側(cè)與右側(cè)風(fēng)道從外到內(nèi)風(fēng)道的速度都呈增大趨勢(shì)，左側(cè)風(fēng)道速度變化范圍比右側(cè)的小，除右側(cè)第一風(fēng)道速度略小于左側(cè)的外，右側(cè)其他各風(fēng)道速度均大于左側(cè)各風(fēng)道速度。

圖12 敞開(kāi)式五風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道平均速度Fig.12 Air channel average velocity of open five-duct negative pressure box at cross section Y=-800 mm

圖13 為敞開(kāi)式五風(fēng)道負(fù)壓箱Y=-800 mm 截面的左右各風(fēng)道平均壓力變化曲線(xiàn)圖。由圖13 可知，左側(cè)各風(fēng)道的壓力變化較小，均保持在-160 Pa 左右，除right1 風(fēng)道的壓力較大之外，右側(cè)其他風(fēng)道壓力與左側(cè)各風(fēng)道壓力差異較小。

圖13 敞開(kāi)式五風(fēng)道負(fù)壓箱Y=-800 mm截面的風(fēng)道平均壓力Fig.13 Air channel average pressure of open five-duct negative pressure box at cross section Y=-800 mm

由以上分析可知，敞開(kāi)式五風(fēng)道結(jié)構(gòu)不僅有效平衡了負(fù)壓箱內(nèi)左右風(fēng)道流量分配比，且其速度場(chǎng)與壓力場(chǎng)分布更均勻，左右風(fēng)道的速度與壓力變化波動(dòng)也較小，即敞開(kāi)式五風(fēng)道負(fù)壓箱的流場(chǎng)比封閉式的更加均勻，有利于干法造紙過(guò)程中纖維均勻沉降。

3.4 排風(fēng)量對(duì)負(fù)壓箱流場(chǎng)均勻性的影響

在干法造紙過(guò)程中，為了適應(yīng)不同的生產(chǎn)需求，負(fù)壓箱通常要在不同排風(fēng)量條件下工作。分別對(duì)排風(fēng)量Q=30000、40000、50000 和 60000 m3/h 時(shí)的敞開(kāi)式五風(fēng)道負(fù)壓箱的流場(chǎng)進(jìn)行仿真計(jì)算，分析負(fù)壓箱內(nèi)流場(chǎng)速度變化規(guī)律，結(jié)果如圖14 所示（取Y=-800 mm截面作分析）。從圖14可以看出，由外到內(nèi)，左右風(fēng)道的速度都呈遞增趨勢(shì)，且隨著排風(fēng)量的增大，各風(fēng)道速度也相應(yīng)增大。

圖14 不同排風(fēng)量的風(fēng)道速度變化Fig.14 Variation of air channel velocity at different exhausting volume

4 負(fù)壓箱流場(chǎng)分布性能的評(píng)價(jià)

評(píng)價(jià)流場(chǎng)均勻性的標(biāo)準(zhǔn)有多種，本課題以速度不均勻系數(shù)評(píng)價(jià)負(fù)壓箱內(nèi)部流場(chǎng)的均勻性，速度不均勻系數(shù)愈小，速度場(chǎng)愈均勻，流場(chǎng)分布也越均勻［15-16］。

速度不均勻性系數(shù)是在負(fù)壓箱內(nèi)選擇n個(gè)測(cè)點(diǎn)，分別測(cè)得各點(diǎn)風(fēng)速，根據(jù)式（4）和式（5）求得其算術(shù)平均值和均方根偏差。

式中，n為測(cè)點(diǎn)總數(shù)；vi為對(duì)應(yīng)測(cè)點(diǎn)的速度；為算術(shù)平均值；σv為均方根偏差。

速度不均勻系數(shù)（kv）可由式（6）求得。

選取X=-900、-600、-300、300、600、900 mm，Y=-600、 -300、 2、 300、 600 mm 和Z=200、 450、680、900 mm 截面的交點(diǎn)（共120 個(gè)）為計(jì)算點(diǎn)，利用后處理軟件在流場(chǎng)仿真結(jié)果中提取計(jì)算點(diǎn)處的速度值，并計(jì)算速度不均勻系數(shù)。

圖15為不同結(jié)構(gòu)負(fù)壓箱Z截面的速度不均勻系數(shù)變化曲線(xiàn)圖。從圖15可以看出，Z=200 mm截面上的速度不均勻系數(shù)較大，Z=900 mm 截面上速度不均勻系數(shù)較小，表明距離負(fù)壓箱進(jìn)風(fēng)口越近，速度不均勻系數(shù)越小，速度場(chǎng)均勻性越好，纖維在成型網(wǎng)上的沉降越均勻。單風(fēng)道負(fù)壓箱各截面速度不均勻系數(shù)在0.39～0.55 之間波動(dòng)，整體小于另外兩種五風(fēng)道結(jié)構(gòu)的負(fù)壓箱，但單風(fēng)道負(fù)壓箱內(nèi)易出現(xiàn)大尺度渦流，不利于纖維均勻沉降。敞開(kāi)式五風(fēng)道負(fù)壓箱各截面速度不均勻系數(shù)均比封閉式的小，其值在0.50～0.55波動(dòng)，說(shuō)明敞開(kāi)式五風(fēng)道負(fù)壓箱內(nèi)流場(chǎng)分布均勻性更好，有利于形成均勻度較好的紙張。

圖15 不同結(jié)構(gòu)負(fù)壓箱的速度不均勻系數(shù)Fig.15 Velocity inhomogeneity coefficients of different negative pressure boxes

圖16 為敞開(kāi)式五風(fēng)道負(fù)壓箱在不同排風(fēng)量時(shí)各截面處（Z截面）的速度不均勻系數(shù)變化曲線(xiàn)圖。從圖16可以看出，不同排風(fēng)量下，各截面的速度不均勻系數(shù)變化趨勢(shì)相同；在各截面處（不包括Z=200 mm、Z=900 mm截面），各排風(fēng)量對(duì)應(yīng)的速度不均勻系數(shù)相差極小；在Z=200 mm、Z=900 mm截面處，各排風(fēng)量對(duì)應(yīng)的速度不均勻系數(shù)范圍分別為0.535～0.554、0.478～0.505，變化范圍也較??；由此可知，排風(fēng)量的變化對(duì)敞開(kāi)式五風(fēng)道負(fù)壓箱的流場(chǎng)均勻性影響較小。因此，在保證紙張成型質(zhì)量的前提下，可適當(dāng)降低風(fēng)機(jī)功率以節(jié)約能源。

圖16 排風(fēng)量對(duì)速度不均勻系數(shù)的影響Fig.16 Effect of exhausting volume on velocity inhomogeneity coefficient

5 結(jié) 論

本課題運(yùn)用Fluent軟件研究了干法造紙負(fù)壓箱結(jié)構(gòu)及排風(fēng)量對(duì)其流場(chǎng)分布特性的影響，并以速度不均勻系數(shù)作為流場(chǎng)分布均勻性的評(píng)價(jià)指標(biāo)，分析了負(fù)壓箱內(nèi)部流場(chǎng)均勻性變化規(guī)律，主要結(jié)論如下。

（1）單風(fēng)道負(fù)壓箱的左右風(fēng)道體積流量比為1∶1.06，但風(fēng)道內(nèi)出現(xiàn)大尺度渦流，不利于纖維均勻沉降。

（2）封閉式五風(fēng)道負(fù)壓箱內(nèi)流場(chǎng)大尺度渦流現(xiàn)象消失，左右風(fēng)道總體積流量比為1.117∶1，流量分配不均勻，且速度不均勻系數(shù)變化范圍較大，其值在0.470～0.635之間波動(dòng)，表明流場(chǎng)分布不均勻。

（3）敞開(kāi)式五風(fēng)道負(fù)壓箱內(nèi)的左右風(fēng)道總體積流量比為1.08∶1，流量分配較均勻，并且速度不均勻系數(shù)在0.535～0.554 之間波動(dòng)，變化范圍較小，表明其流場(chǎng)分布均勻性較好。

（4）隨著排風(fēng)量增大，敞開(kāi)式五風(fēng)道負(fù)壓箱內(nèi)各風(fēng)道的速度相應(yīng)增大，但流場(chǎng)速度不均勻系數(shù)變化較小，因此，工程實(shí)踐中可適當(dāng)降低風(fēng)機(jī)功率以節(jié)約能源。