陳作炳，項(xiàng)勤，呂銀雷，鄒遠(yuǎn)志，印世杰

?

立式有機(jī)廢物發(fā)酵系統(tǒng)攪拌器的阻力分析

陳作炳，項(xiàng)勤，呂銀雷，鄒遠(yuǎn)志，印世杰

(武漢理工大學(xué)機(jī)電工程學(xué)院，湖北武漢430070)

研究開發(fā)了一種立式有機(jī)廢棄物發(fā)酵系統(tǒng)，系統(tǒng)由主軸攪拌器、上料裝置、通風(fēng)系統(tǒng)、傳動(dòng)機(jī)構(gòu)、卸料裝置組成。系統(tǒng)通過上料裝置進(jìn)料并接入適量菌種，由主軸攪拌器攪拌混合，通風(fēng)系統(tǒng)提供足量新鮮空氣，傳動(dòng)機(jī)構(gòu)為系統(tǒng)提供動(dòng)力，待發(fā)酵完成由卸料裝置完成卸料。針對(duì)系統(tǒng)在固定工況條件下主軸攪拌器發(fā)酵攪拌過程中的阻力展開研究，運(yùn)用現(xiàn)代化離散元模型(EDEM)模擬和分析物料在筒體內(nèi)的運(yùn)動(dòng)和攪拌器受力情況，并建立攪拌器數(shù)學(xué)模型，對(duì)攪拌阻力進(jìn)行理論計(jì)算。結(jié)果表明：EDEM數(shù)值模擬計(jì)算得到攪拌器轉(zhuǎn)矩平均值為84 927.6 N?m，理論計(jì)算得到攪拌器理論轉(zhuǎn)矩計(jì)算值為86 549 N?m，比模擬值略大，與數(shù)值模擬計(jì)算結(jié)果具有一致性。

有機(jī)廢物發(fā)酵系統(tǒng)；攪拌器；阻力分析

目前，有機(jī)廢棄物處理方式主要有衛(wèi)生填埋法和焚燒法和傳統(tǒng)堆肥法[1]。衛(wèi)生填埋法容易引起環(huán)境污染，打破環(huán)境平衡，惡化土質(zhì)；焚燒法產(chǎn)生煙氣劇毒物質(zhì)，焚燒廠選址困難[2–4]；傳統(tǒng)堆肥法存在堆肥反應(yīng)時(shí)間長(zhǎng)、成本高、堆肥效果不理想等問題[5]。筆者根據(jù)有機(jī)廢棄物處理需遵循的資源化、無害化、減量化原則[6–7]，開發(fā)設(shè)計(jì)出一種立式有機(jī)廢棄物發(fā)酵系統(tǒng)，系統(tǒng)5~7 d可快速完成發(fā)酵，成本較低。筆者針對(duì)系統(tǒng)關(guān)鍵結(jié)構(gòu)主軸攪拌器攪拌過程中的阻力進(jìn)行分析，運(yùn)用離散元模型(EDEM)模擬和分析物料的運(yùn)動(dòng)情況，現(xiàn)將結(jié)果報(bào)道如下。

1　發(fā)酵系統(tǒng)總體結(jié)構(gòu)

系統(tǒng)由主軸攪拌器、上料裝置、通風(fēng)系統(tǒng)、傳動(dòng)機(jī)構(gòu)以及卸料裝置5部分組成，如圖1所示。上料裝置進(jìn)料并接入適量菌種后，由主軸攪拌器攪拌混合，通風(fēng)系統(tǒng)通入足量新鮮空氣，傳動(dòng)機(jī)構(gòu)為系統(tǒng)提供動(dòng)力，待發(fā)酵完成由后卸料裝置完成卸料。

圖1　立式有機(jī)廢物發(fā)酵系統(tǒng)的結(jié)構(gòu)

發(fā)酵系統(tǒng)筒體內(nèi)徑2.7 m，筒體高4.6 m，總高8.2 m。主軸攪拌器為中空不銹鋼管，可通過主軸通入新鮮空氣，主軸攪拌器轉(zhuǎn)速設(shè)定為4 r/h。轉(zhuǎn)速可以根據(jù)有機(jī)廢棄物材料特性和外部環(huán)境的需要進(jìn)行調(diào)節(jié)，調(diào)節(jié)范圍2~8 r/h。

2　攪拌器的力學(xué)分析

主軸攪拌器是系統(tǒng)的關(guān)鍵部件，針對(duì)主軸攪拌器發(fā)酵攪拌過程中的阻力展開研究。將有機(jī)廢棄物料視為一系列獨(dú)立運(yùn)動(dòng)的離散單元[8]，通過離散元模型模擬物料運(yùn)動(dòng)情況和受力分析。

2.1　幾何模型的建立

通過CREO建立發(fā)酵系統(tǒng)筒體簡(jiǎn)化模型，如圖2所示。對(duì)簡(jiǎn)化模型進(jìn)行網(wǎng)格劃分，如圖3所示。

圖2　系統(tǒng)簡(jiǎn)化模型

圖3　網(wǎng)格劃分模型

2.2　計(jì)算模型的選取

有機(jī)廢棄物料發(fā)酵過程中，用Hertz– Mindlin (no–slip)接觸模型[9–10]模擬其內(nèi)部的運(yùn)動(dòng)過程。

2.3　材料特性參數(shù)確定

在進(jìn)行有機(jī)廢棄物料EDEM仿真計(jì)算時(shí)，經(jīng)過相似理論轉(zhuǎn)換后的模型需要滿足計(jì)算機(jī)的解算條件限制。實(shí)際物理模型和轉(zhuǎn)換后的EDEM模型參數(shù)基本保持一致，同時(shí)需要對(duì)有機(jī)廢棄物料參數(shù)進(jìn)行確定。通過查閱相關(guān)文獻(xiàn)資料、材料手冊(cè)和參考類似物料性質(zhì)[11–14]，將有機(jī)廢棄物料和發(fā)酵菌種的特征參數(shù)確定如表1所示。

表1　有機(jī)廢棄物料和發(fā)酵菌種的特性參數(shù)

實(shí)際生產(chǎn)中，發(fā)酵菌種需要量很小，約為0.1%，但在進(jìn)行EDEM數(shù)值模擬的時(shí)候，考慮到模型的代表性和典型性，物料含量過小會(huì)導(dǎo)致模擬難以正常進(jìn)行，分析結(jié)果也會(huì)不準(zhǔn)確，因此在EDEM數(shù)值模擬過程中將發(fā)酵菌種與有機(jī)廢棄物料顆粒數(shù)調(diào)整為1∶8，以保證數(shù)值模擬過程的順利進(jìn)行。

2.4　結(jié)果與分析

經(jīng)過EDEM數(shù)值模擬計(jì)算，發(fā)酵系統(tǒng)攪拌效果如圖4所示。開始攪拌前，物料分批次進(jìn)料，發(fā)酵菌種和有機(jī)廢棄物料均勻散落在每一層表面。圖4–a中，綠色代表有機(jī)廢棄物料，紅色代表發(fā)酵菌種。筒體內(nèi)2種物料總共分為4層，界限明顯；攪拌完成后，發(fā)酵菌種均勻分布在整個(gè)筒體，達(dá)到了有機(jī)廢物與發(fā)酵菌種攪拌均勻的效果，如圖4–b所示。

a　攪拌前 b　攪拌后

圖5為攪拌主軸各部位受力圖，是由Hertz– Mindlin(no–slip)接觸模型中的法向力和切向力經(jīng)過數(shù)值模擬計(jì)算得到的結(jié)果。圖5–a為攪拌軸受力云圖，由于物料主要集中在下部區(qū)域，因此底層攪拌葉片既需要推動(dòng)物料作旋轉(zhuǎn)運(yùn)動(dòng)，又需要承受來自上部物料的壓力作剪切運(yùn)動(dòng)，受力較大；隨著有機(jī)廢物高度的上升，攪拌葉片受力逐漸減小，因此需要對(duì)底部葉片設(shè)置加強(qiáng)筋板來保證攪拌葉片強(qiáng)度。圖5–b為攪拌軸受力處理云圖，由圖可知，主要受力集中部位在攪拌葉片頂端，這是由于攪拌葉片離主軸越遠(yuǎn)，旋轉(zhuǎn)半徑越大，需要帶動(dòng)的物料也越多，因此受力越大；受力面主要為與物料直接接觸的剪切面，因而需要考慮攪拌葉片主切割板的材質(zhì)，減輕剪切面磨損。

a攪拌軸受力圖

b攪拌軸受力處理圖

圖5攪拌軸受力

Fig.5Force diagram of agitator’s shaft

3　攪拌器阻力的理論計(jì)算

發(fā)酵系統(tǒng)容積為27 m3，可裝入21 t左右有機(jī)廢棄物原料，攪拌裝置對(duì)物料進(jìn)行攪拌混合過程，主軸所受扭矩較大。通過對(duì)攪拌裝置建立工作的數(shù)學(xué)模型，可計(jì)算出整機(jī)所需要的扭矩，以便對(duì)傳動(dòng)系統(tǒng)中關(guān)鍵部件進(jìn)行力學(xué)分析和優(yōu)化設(shè)計(jì)。

3.1　攪拌機(jī)構(gòu)數(shù)學(xué)模型

攪拌主軸葉片對(duì)有機(jī)廢棄物料和發(fā)酵菌種進(jìn)行攪拌、混合[15]，是主要的受力構(gòu)件，葉片在運(yùn)動(dòng)時(shí)受力情況如圖6所示。葉片在沿中心線回轉(zhuǎn)作攪拌運(yùn)動(dòng)時(shí)，主要承受上部物料施加給葉片的重力、葉片帶動(dòng)運(yùn)動(dòng)方向物料的擠壓力和葉片附近物料對(duì)葉片的摩擦力F，攪拌主軸上呈螺旋90°均勻布置8層葉片，為使攪拌效果更好，葉片軸向分布為下部緊密，上部稀疏，因此每層葉片受力不同，同一葉片沿徑向方向各點(diǎn)對(duì)主軸產(chǎn)生的轉(zhuǎn)矩也不相同。

圖6　運(yùn)動(dòng)時(shí)攪拌主軸葉片的受力

攪拌軸由底部傳動(dòng)系統(tǒng)驅(qū)動(dòng)，主軸轉(zhuǎn)動(dòng)時(shí)，物料對(duì)葉片產(chǎn)生作用力，葉片對(duì)主軸產(chǎn)生轉(zhuǎn)矩，轉(zhuǎn)矩的大小對(duì)于傳動(dòng)系統(tǒng)選型、動(dòng)力設(shè)計(jì)有重要指導(dǎo)意義。以筒體內(nèi)所有物料為研究對(duì)象，攪拌裝置推動(dòng)物料在筒體內(nèi)運(yùn)動(dòng)，攪拌裝置簡(jiǎn)化后的機(jī)構(gòu)運(yùn)動(dòng)簡(jiǎn)圖，如圖7所示。

圖7　攪拌器

3.2　攪拌主軸扭矩的理論計(jì)算

筒體內(nèi)部物料對(duì)主軸產(chǎn)生的扭矩可分為2部分：一是葉片向前推動(dòng)物料進(jìn)行攪拌混合時(shí)對(duì)主軸產(chǎn)生的扭矩；二是葉片運(yùn)動(dòng)過程中與物料之間摩擦阻力所產(chǎn)生的扭矩。

3.2.1推動(dòng)物料運(yùn)動(dòng)產(chǎn)生的扭矩

分析攪拌裝置在筒體內(nèi)的運(yùn)動(dòng)時(shí)，為了分析方便，假設(shè)物料是不可壓縮的，采用瞬時(shí)穩(wěn)態(tài)假定方法，將筒體內(nèi)物料均勻劃分為個(gè)微單元，取筒體內(nèi)任意一微元物料作為研究對(duì)象，如圖8所示。微元距離筒體中心線為，微元長(zhǎng)為d，高為d，寬為dα。由于筒體內(nèi)葉片沿軸線方向按層分布，葉片在筒體內(nèi)運(yùn)動(dòng)時(shí)，只能帶動(dòng)部分物料運(yùn)動(dòng)，主要作用在葉片上下附近的物料，2層葉片之間部分物料沒有運(yùn)動(dòng)，記1為物料的運(yùn)動(dòng)系數(shù)，即筒體內(nèi)葉片推動(dòng)的物料體積與總物料體積之比。

圖8　積分單元

式中：1為葉片可撥動(dòng)筒體內(nèi)物料的高度(m)；1為筒體內(nèi)物料的高度(m)；1為攪拌葉片所能帶動(dòng)物料的高度(m)；2為刮料葉片所能帶動(dòng)物料的高度(m)；為攪拌葉片層數(shù)。

由于1=3.9 m，1=0.1 m，2=0.047 m，=8，因而1=0.847 m，1=0.217。

則物料對(duì)軸的扭矩為：

式中：d1為微元對(duì)于軸的扭矩(N?m)；1為筒體內(nèi)物料對(duì)于軸的扭矩(N?m)；為微元距離軸線的距離(m)；d為微元長(zhǎng)度，m；d為微元高度(m)；為微元夾角(°)。

對(duì)(2)式積分得：

3.2.2葉片與物料之間摩擦阻力產(chǎn)生的扭矩

當(dāng)攪拌裝置在筒體內(nèi)推動(dòng)物料運(yùn)動(dòng)時(shí)，葉片與上下層物料發(fā)生相對(duì)運(yùn)動(dòng)，物料與葉片之間產(chǎn)生摩擦力，摩擦力對(duì)主軸產(chǎn)生扭矩。把攪拌裝置作為一個(gè)整體來分析，則葉片承受筒體內(nèi)所有物料重力，簡(jiǎn)化筒體內(nèi)葉片受力如圖9，物料對(duì)葉片的摩擦力作用于半徑中點(diǎn)。

圖9　摩擦阻力產(chǎn)生的扭矩

葉片受到的正壓力等于筒體內(nèi)物料重力。

葉片所受摩擦力為：

則摩擦力產(chǎn)生的扭矩為：

即：

所以，攪拌裝置的主軸承受扭矩為：

(9)

式中：1為材料與物料之間的摩擦系數(shù)，1=0.4。

3.3　攪拌阻力

3.3.1數(shù)值計(jì)算

圖10為攪拌軸在攪拌過程中所承受的扭矩變化曲線。0~5 s時(shí)，為下料階段，物料在下落過程中對(duì)攪拌葉片有沖擊，攪拌葉片受到很小的轉(zhuǎn)矩；下料完成后，攪拌軸開始轉(zhuǎn)動(dòng)，葉片帶動(dòng)物料在筒體內(nèi)作回轉(zhuǎn)運(yùn)動(dòng)，受到的阻力增加，導(dǎo)致轉(zhuǎn)矩急劇上升達(dá)到76 800.1 N·m，當(dāng)攪拌軸攪拌速度穩(wěn)定后，攪拌軸承受總轉(zhuǎn)矩在85 000 N·m上下浮動(dòng)，這是由于物料在攪拌過程中，物料在筒體內(nèi)分布不均勻造成的，因此攪拌葉片承受的轉(zhuǎn)矩?cái)?shù)值以一穩(wěn)定值為基礎(chǔ)上下浮動(dòng)。由EDEM導(dǎo)出的數(shù)據(jù)顯示，攪拌軸在攪拌過程中承受的總轉(zhuǎn)矩平均值為84 927.6 N·m。

圖10　攪拌過程中攪拌軸轉(zhuǎn)矩的變化

3.3.2理論計(jì)算結(jié)果分析

發(fā)酵系統(tǒng)正常工作時(shí)，筒體上方需預(yù)留出40~50 cm高的空間，用于氣體流動(dòng)、熱交換、水分蒸發(fā)等。當(dāng)加入筒體內(nèi)物料高度為3.9 m時(shí)，攪拌裝置主軸所受扭矩為：

=86 549 N·m (10)

理論計(jì)算值略大于模擬計(jì)算值，具有一致性。在理論計(jì)算時(shí)，對(duì)葉片受力進(jìn)行了假設(shè)，把攪拌葉片作為一個(gè)整體來分析，葉片承受筒體內(nèi)所有物料重力，而實(shí)際情況物料之間也存在相互作用力，因此葉片只承受筒體內(nèi)物料大部分重力，因此數(shù)值模擬計(jì)算值比理論計(jì)算值略小，符合生產(chǎn)實(shí)際。

[1] 趙秀玲，朱新萍，羅艷麗，等．溫度與秸稈比例對(duì)牛糞好氧堆肥的影響[J]．環(huán)境工程學(xué)報(bào)，2014，8(1)：334–340．

[2] ZHU B，GIKAS P，ZHANG R，et al．Evalation of the rotary drum reactor as municipal solid waste pretreament technology for biomass utilization[J]. International Conference on Biomass Energy Technologies，2008 (10)：160–168．

[3] 張海霞，李愛民，鞠茂偉．小型好氧堆肥設(shè)備處理有機(jī)垃圾[J]．環(huán)境工程學(xué)報(bào)，2013，7(5)：1939–1944．

[4] 紀(jì)丹鳳，夏訓(xùn)峰，劉駿，等．北京市生活垃圾處理的環(huán)境影響評(píng)價(jià)[J]．環(huán)境工程學(xué)報(bào)，2011，5(9)：2101–2107．

[5] 王玉新，任勇翔，王曉昌，等．塔式自然通風(fēng)好氧堆肥反應(yīng)器的開發(fā)與應(yīng)用[J]．環(huán)境工程學(xué)報(bào)，2014，8(4)：1631–1636．

[6] 黃華，黃懿梅，劉學(xué)玲，等．覆蓋處理對(duì)豬糞秸稈堆肥中氮素轉(zhuǎn)化和堆肥質(zhì)量的影響[J]．環(huán)境科學(xué)學(xué)報(bào)，2013，33(3)：780–786．

[7] BHARAT Bhushan，MADHUSUDAN Singh．Adaptive control of DC motor using bacterial foraging algorithm[J]. Applied Soft Computing，2011，11：4913–4920．

[8] 胡國(guó)明．顆粒系統(tǒng)的離散單元法分析仿真[M]．武漢：武漢理工大學(xué)出版社，2010．

[9] 王國(guó)強(qiáng)，郝萬軍，王繼新．離散單元法及其在EDEM上的實(shí)踐[M]．西安：西北工業(yè)大學(xué)出版社，2010．

[10] Paul W Cleary．Industrial particle flow modeling using discrete element method[J]．Engineering Conputationgs，2009，26(6)：698–743．

[11] 李季，彭生平．堆肥工程實(shí)用手冊(cè)[K]．2版．北京：化學(xué)工業(yè)出版社，2011．

[12] HAUG R T．The Practical Handbook of Compost Engineering[K]．Boca Raton Fla：Lewis Publishers，1993．

[13] 高云航，勾長(zhǎng)龍，王雨瓊，等．低溫復(fù)合菌劑對(duì)牛糞堆肥發(fā)酵影響的研究[J]．環(huán)境科學(xué)學(xué)報(bào)，2014，34(12)：3166–3170．

[14] BAO D J，SHU Y W，LI Q X，et al．The effect of salinity on waste activated sludge alkaline fermentation and kinetic analysis[J]．Journal of Environmental Sciences，2016(5)：80–90．

[15] 陳志平．?dāng)嚢枧c混合設(shè)備設(shè)計(jì)選用手冊(cè)[K]．北京：化學(xué)工業(yè)出版社，2004．

責(zé)任編輯：羅慧敏

英文編輯：吳志立

Resistance analysis of the agitator in the vertical organic waste fermentation treatment system

CHEN Zuobing, XIANG Qin, Lü Yinlei, ZOU Yuanzhi, YIN Shijie

(School of Mechanical and Electronic Engineering, Wuhan University of Technology, Wuhan 430070, China)

An vertical organic waste fermentation treatment system has been developed, which composed of main agitator, feeding device, ventilation system, transmission, and discharge part. The organic waste was fed through the feeding device as well as the appropriate bacteria, and then mixed by the main agitator. The ventilation system provided adequate fresh air to ensure the organic waste fermenting completely. The transmission transmitted power to the whole system. When the fermentation is finished, the organic waste was discharged by the discharge part. Resistance of the main agitator was studied during the fermentation process under normal condition. The mathematical model of the agitator was established to simulate and analyze the movement of the material in the cylinder and the force on the agitator by discrete element soft of EDEM. The results show that the average torque on the main agitator was 84 927.6 N?m obtained by simulation, which slightly larger than the value of 86 549 N?m from theoretical calculation. It is found that the simulation by EDEM is consistency with theoretical calculation.

organic waste fermentation system; agitator; resistance analysis

X705

A

1007-1032(2017)05-0570-05

2017–02–24

2017–07–05

中國(guó)農(nóng)業(yè)科學(xué)院水體污染控制與治理科技重大專項(xiàng)(2017ZX07603–002)

陳作炳(1962— )，男，湖北孝感人，教授，主要從事節(jié)能環(huán)保技術(shù)及其應(yīng)用研究，zbchen01@126.com

投稿網(wǎng)址：http://xb.hunau.edu.cn