基于FLUENT的污泥脫水螺旋壓榨機流場分析

張念猛 曲慶文 陳云華

(山東理工大學(xué)機械工程學(xué)院，山東 淄博，255049)

造紙污泥脫水是造紙污水處理過程中的關(guān)鍵環(huán)節(jié)，低效的污泥脫水將帶來很多不良后果，包括環(huán)境污染、水資源浪費等，采用效率更高的設(shè)備能夠有效地提高污泥脫水效率[1］。目前已經(jīng)有學(xué)者在造紙污水處理方面做了深度研究[2］，而在機械脫水領(lǐng)域，將螺旋壓榨機應(yīng)用于造紙污泥脫水已在國外悄然興起，實際生產(chǎn)表明螺旋壓榨機的污泥脫水率要優(yōu)于其他設(shè)備[3］。由于該設(shè)備剛剛進入污泥脫水領(lǐng)域，雖然使用前景廣闊，但由于生產(chǎn)和使用經(jīng)驗不足，導(dǎo)致該設(shè)備的結(jié)構(gòu)設(shè)計方面還存在較多不足[4］。有學(xué)者利用ANASYS軟件對螺旋壓榨機螺旋軸進行了有限元分析[5］，但在其流場分析領(lǐng)域仍顯薄弱。在與螺旋壓榨機結(jié)構(gòu)相似的螺旋輸送機螺旋軸設(shè)計中，有學(xué)者利用了流場流量與流速的計算理論[6］，為簡化復(fù)雜計算，有學(xué)者也提出了合理的條件假設(shè)[7］。本文依據(jù)脫水參數(shù)要求，建立螺旋軸等外徑變螺距螺旋線參數(shù)的計算方程，根據(jù)該方程建立螺旋壓榨機螺旋軸的三維模型，并導(dǎo)入FLUENT軟件對螺旋壓榨機的內(nèi)部流場進行仿真模擬，以期為螺旋壓榨機的結(jié)構(gòu)優(yōu)化提供參考。

1 榨料流量運算方程的建立

在進行榨料流量運算方程推導(dǎo)之前，為方便過程分析，簡化復(fù)雜計算，提出如下假設(shè):①榨料填充連續(xù)、規(guī)則且盡量緊密;②不考慮壓榨的過渡狀態(tài)，即造紙污泥一進入螺旋壓榨機就開始壓榨脫水過程;③榨料沿螺旋軸軸向的運動速度等于螺旋葉片軸向的旋進速度。

1.1 榨料軸向運動速度分析

假設(shè)螺旋軸以均勻的角速度ω旋轉(zhuǎn)，在其截面任一半徑為r的P點處有一榨料質(zhì)點，該榨料質(zhì)點有沿葉片表面的相對滑動和沿螺旋軸軸向的移動，其運動速度分析如圖1所示。

圖1 螺旋葉片上某物料質(zhì)點速度分析圖

由速度三角形可得，矢量PA即表示葉片上P點旋轉(zhuǎn)的圓周速度vo，其中:

vo方向沿P點回轉(zhuǎn)的切線方向;PB為物料質(zhì)點的絕對運動速度va，方向是沿P點螺旋線的法線方向，且與螺旋軸軸向的夾角為α，易知:

將va分解為沿螺旋軸軸向的軸向速度vp和垂直螺旋軸軸向的速度vt，其中:

又由幾何關(guān)系可得半徑r的計算公式為:

聯(lián)立式 (1)、式 (2)、式 (3)、式 (4)可得:

vp即為螺旋軸截面上任一半徑r處的榨料質(zhì)點的軸向運動速度，其表達式為螺旋轉(zhuǎn)角α的函數(shù)，P為本段螺旋線的螺距。根據(jù)前面假設(shè)，如果相鄰兩段螺距變化不大，則可認為在整個螺旋截面區(qū)域內(nèi)榨料的軸向運動速度都相同，其表達式為:

1.2 榨料流量分析

由榨料流量 (Q)的定義可知:

式中，dσ為面積元，且dσ=2πrdr。解式 (7)可求得:

式中，Q為榨料流量 (m3/h);R為螺旋葉片外半徑 (mm);r為螺旋葉片內(nèi)半徑 (mm)。

2 變螺距運算方程的推導(dǎo)

2.1 基于流量的螺距方程

螺旋壓榨機兩端為進料口和出料口，設(shè)進料口流量為Qi，流入污泥的含水率為Ci，出料口流量為Qo，流出污泥的含水率為Co。由于壓榨過程是將多余水分榨出，而固態(tài)污泥總量基本不變，故從進料口到出料口的過程中流量是不斷減小的，且Qi和Qo之間滿足方程:

根據(jù)上述理論以進料端為坐標(biāo)原點，建立如圖2所示坐標(biāo)系。

圖2 螺旋軸坐標(biāo)系

設(shè)距原點L1處的螺旋截面流量為Q1，在L1+ΔL處的流量為Q2，則其流量差為ΔQ=Q1－Q2，當(dāng)ΔL趨于無窮小時，其單位長度的榨料量應(yīng)視為常數(shù)，即:

2.2 基于螺旋轉(zhuǎn)角的螺距方程

對于以變螺距螺旋線為線型的螺旋軸，當(dāng)它以ω的角速度勻速轉(zhuǎn)動時，在Δt的時間間隔內(nèi)，螺旋軸相應(yīng)地轉(zhuǎn)過一個角度Δα，由角速度的定義可知:

由于式 (22)為前期理論推導(dǎo)所得，在實際應(yīng)用時需引入修正系數(shù)λ，則最終得到螺旋線螺距方程為:

其中λ根據(jù)螺旋壓榨機實際工作情況確定。

3 螺旋壓榨機內(nèi)部流場仿真

3.1 螺旋軸三維模型的建立

根據(jù)實際生產(chǎn)要求設(shè)定已知參數(shù)如下:進料口流量Qi=30 m3/h，流入污泥含水率Ci=95%，流出污泥含水率Co=65%，螺旋葉片外半徑R=400 mm，中間軸大端半徑Rb=300 mm，小端半徑Rs=200 mm，螺旋軸長度L=2000 mm，螺旋軸轉(zhuǎn)速n=2 r/min。

由式 (9)可計算出出料口流量Qo=4.3 m3/h，螺旋軸旋轉(zhuǎn)角速度ω=2πn=4π r/min，則由進出口流量及相關(guān)數(shù)據(jù)可求得P0=816 mm，Pn+1=152 mm，這里的Pn+1也是參數(shù)螺距，不代表實際結(jié)構(gòu)尺寸。由式 (14)可求得ε=0.35。由式 (23)得第一段螺距為P1=547 mm，以此類推P2=485 mm，P3=342 mm，P4=241 mm，P5=170 mm。

根據(jù)以上螺旋軸的尺寸參數(shù)，建立螺旋軸的三維模型，如圖3所示。

圖3 螺旋軸的三維模型圖

3.2 內(nèi)部流場的仿真分析

此螺旋壓榨機處理的污水主要含水、泥兩相介質(zhì)，其中污泥體積系數(shù)為5%，密度為1051 kg/m3，黏度為0.02 Pa·s。定義邊界條件為速度入口，且入口處流速為0.087 m/s，出口處為自由流出，且出口處為標(biāo)準(zhǔn)大氣壓。根據(jù)螺旋壓榨機的整機及螺旋軸尺寸，建立內(nèi)部流場的三維模型，并將其導(dǎo)入GAMBIT軟件劃分網(wǎng)格，如圖4所示。

圖4 流場模型網(wǎng)格劃分圖

經(jīng)過FLUENT軟件的仿真計算，得到流場中內(nèi)、外壁面以及中軸面的壓強分布云圖，如圖5所示。

由圖5(a)、圖5(b)可以看出，隨著壓榨過程的進行，壓強自入口向出口逐漸增大;由圖5(c)可知，在螺旋葉片附近的流體壓強要高于其他區(qū)域的壓強，而螺旋葉片根部的壓強又要高于葉片頂部的壓強。螺旋軸及葉片的材料為40Cr，其許用應(yīng)力為σn=400 MPa，取安全系數(shù)為2，則其基本許用應(yīng)力為σn/2=200 MPa，由圖5中可以看出，最大壓強出現(xiàn)在出口處最后一圈螺旋葉片的根部，且該壓強值遠低于材料的許用應(yīng)力。對于分裝式的螺旋葉片和中間軸，對焊接質(zhì)量有較高要求，該機器螺旋葉片的強度總體滿足工作要求，安全性較高，為下一步的優(yōu)化設(shè)計提供了理論依據(jù)。

圖5 流場壓強分布云圖

圖6 流場速度分布云圖

流場的速度分布云圖如圖6所示。從圖6可以看出，機器中的流體流速整體較慢，并且自入口向出口逐漸減小，其壓榨過程主要是依靠流體自身的重力過濾以及不斷減小的螺旋腔對流體進行逐層壓榨，較小的壓榨速度可保證壓榨過程充分進行。

4 結(jié)論

4.1 根據(jù)進出口榨料相等理論，得到了等外徑變螺距螺旋壓榨螺旋軸的數(shù)學(xué)模型，根據(jù)脫水參數(shù)，可求出所需螺旋線參數(shù)方程，為螺旋壓榨機螺旋軸的結(jié)構(gòu)設(shè)計提供了依據(jù)。

4.2 上述模型推導(dǎo)過程表明，螺旋線的螺旋轉(zhuǎn)角α隨軸向L按指數(shù)規(guī)律變化，螺旋線螺距P沿軸向L線性減小。

4.3 通過FLUENT軟件分析，得到流場壓強的整體分布規(guī)律，壓強自入口向出口逐漸增大，在螺旋葉片附近的流體壓強要高于其他區(qū)域的壓強，螺旋葉片根部的壓強高于葉片頂端的壓強，

4.4 通過分析得到流場速度整體較小，且自入口向出口逐漸減小，保證了壓榨過程的充分進行，實用性良好。

[1]ZHANG Yong，CAO Chun-yun，F(xiàn)ENG Wen-ying，et al.Advances in Pollution Controls Science and Technology of China'Paper Industry[J］.China Pulp ＆ Paper，2012，31(2):57.張 勇，曹春昱，馮文英，等.我國制漿污染治理科學(xué)技術(shù)的現(xiàn)狀與發(fā)展[J］.中國造紙，2012，31(2):57.

[2]ZHAI Yang.Current Status and Development Trend of Advanced Treatment Technologies for Waste Water of Paper Industry in China[J］.China Pulp ＆ Paper，2011，30(10):56.翟 陽.我國造紙工業(yè)廢水深度處理的技術(shù)現(xiàn)狀及其發(fā)展趨勢[J］.中國造紙，2011，30(10):56.

[3]REN Tie，SUN Xiao，ZHANG Xiaodong.The application of screw press dewaterer in wastewater treatment[J］.Environmental Engineering，2008，26(2):84.任 鐵，孫 曉，張曉東.螺壓脫水機在污水處理中的應(yīng)用[J］.環(huán)境工程，2008，26(2):84.

[4]ZHAO Chuan-shan，LI Rong-gang，WU Chao-jun.Research on some problems of structural design for screw press[J］.Journal of Shandong Institute of Light Industry，2010，24(4):1.趙傳山，李榮剛，吳朝軍.對螺旋壓榨脫水機結(jié)構(gòu)設(shè)計中幾個問題的探討[J］.山東輕工業(yè)學(xué)院學(xué)報，2010，24(4):1.

[5]DENG Quan-de，LU Bo-you.The strength analysis of screw press applied in agricultural materials based on ANASYS[J］.Journal of Agricultural Mechanization Research，2013(8):58.鄧全得，盧博友.基于ANASYS的農(nóng)業(yè)物料螺旋壓榨機螺旋強度分析[J］.農(nóng)機化研究，2013(8):58.

[6]SONG Ya-ying.Theory study on structure parameters of variable-pith augur[J］.Journal of Suzhou University(Natural Science)，1996，12(4):96.宋亞英.變螺距螺旋結(jié)構(gòu)參數(shù)的理論研究[J］.蘇州大學(xué)學(xué)報，1996，12(4):96.

[7]WANG Lian-dong，Parametric modeling and structural analysis of the press screw with variable diameters and variable pitches[D］.Xi'an:Northwest A＆F University，2010.王連棟.變徑變螺距榨螺的參數(shù)化建模與結(jié)構(gòu)分析[D］.西安:西北農(nóng)林科技大學(xué)，2010.