基于ANSYS Workbench的復(fù)合式發(fā)酵塔裝置重要零部件有限元分析

石義棟 馮天玉 陸偉安

摘要：以復(fù)合式發(fā)酵（好氧厭氧結(jié)合）塔裝置為研究對象，通過Solidworks建立復(fù)合式發(fā)酵塔三維實(shí)體模型，并對復(fù)合式發(fā)酵塔進(jìn)行了簡化。利用ANSYS Workbench有限元分析軟件識別結(jié)合面特征，建立復(fù)合式發(fā)酵塔裝置重要零部件的有限元模型。通過理論計算與有限元分析的結(jié)合，分析得出螺旋桿軸受力較小，葉片受力大，葉片兩端受力最大，隔板右上端面厚度太小以至于承受不住物料的壓力。

關(guān)鍵詞：復(fù)合式發(fā)酵塔裝置重要零部件;受力分析;有限元分析;強(qiáng)度;ANSYS Workbench

中圖分類號：S232;TH16? ? ? ? ?文獻(xiàn)標(biāo)識碼：A

文章編號：0439-8114（2019）10-0150-04

Abstract： Taking the composite fermentation （aerobic anaerobic combined） tower device as the research object， the solid three-dimensional solid model of the composite fermentation tower was established by Solidworks， and the composite fermentation tower was simplified. The ANSYS Workbench finite element analysis software was used to identify the characteristics of the joint surface， and the finite element model of the important parts of the composite fermentation tower was established. Through the combination of theoretical calculation and finite element analysis， it is found that the screw shaft has less force， the blade is subjected to large force， and the force at both ends of the blade is the largest. The thickness of the upper right end of the diaphragm is too small to withstand the pressure of the material.

Key words： compound fermentation tower equipment important parts; force analysis; finite element analysis; strength; ANSYS Workbench

隨著畜禽養(yǎng)殖業(yè)的迅速發(fā)展，畜牧養(yǎng)殖規(guī)?；?、集約化和產(chǎn)業(yè)化的提高，中國接連幾年在禽肉、豬肉、蛋類產(chǎn)量均位列世界第一，但帶來的環(huán)境問題也相對不斷增加，越來越多的畜禽排泄物未能得到及時處理，占用了大量土地，導(dǎo)致中國城郊、農(nóng)村、城鎮(zhèn)環(huán)境污染問題日趨嚴(yán)重，已逐漸發(fā)展成不可忽視的環(huán)境和生態(tài)保護(hù)問題。據(jù)《經(jīng)濟(jì)日報》報道，中國畜禽養(yǎng)殖每年產(chǎn)生38億t畜禽糞便，有效處理率不足50%;湖北省每年各類農(nóng)作物秸稈理論資源量3 680萬t，可收集的資源量3 300萬t，但每年仍然有大量秸稈被焚燒或丟棄，得不到利用[1]。目前，國內(nèi)普遍采用床式發(fā)酵（在畜禽居住區(qū)放入墊料，加入菌種攪拌均勻）或堆肥發(fā)酵（用各種植物殘體、作物秸稈、雜草、樹葉、泥炭、垃圾以及其他廢棄物等）為主要原料，混合人畜糞尿經(jīng)堆制腐解），該類技術(shù)雖能實(shí)現(xiàn)畜禽糞便轉(zhuǎn)化有機(jī)肥，但均存在較大的問題：場地占用土地面積大，無組織排放，危害周邊養(yǎng)殖人員、畜禽以及居住人員，污染水資源和大氣環(huán)境，發(fā)酵周期長。

針對上述問題，研發(fā)一種復(fù)合式發(fā)酵裝置，在塔內(nèi)分層管理，進(jìn)行好氧、厭氧的全密閉式發(fā)酵，降解抗生素，脫硫脫硝，除臭，鈍化重金屬殺滅病菌及蟲卵，以期徹底解決上述問題。以復(fù)合式發(fā)酵塔裝置為例，利用Solidworks建立復(fù)合式發(fā)酵塔裝置重要零部件的三維模型，并用有限元分析軟件ANSYS Workbench對螺旋桿進(jìn)行了變形、應(yīng)力計算，獲得了合理的分析結(jié)果，既為復(fù)合式發(fā)酵塔重要零部件的進(jìn)一步優(yōu)化提供了理論基礎(chǔ)，也為復(fù)合式發(fā)酵塔裝置重要零部件的正常運(yùn)行提供了保障。

1? 復(fù)合式發(fā)酵塔重要零部件三維建模

螺旋葉片桿：該螺旋葉片桿全長5 960 mm，螺距160 mm，葉片高度100 mm，材料為45鋼，三維模型如圖1所示。

隔板：中間內(nèi)外圓直徑分別為100、108 mm，垂直方向兩邊長度分別為80、312 mm;垂直和斜邊方向厚度分別為12、4 mm，中心桿長2 800 mm，板長2 776 mm，三維模型如圖2所示。

2? 螺旋葉片的基本參數(shù)

3? 螺旋葉片桿及隔板的理論計算

3.1? 螺旋葉片理論計算

該螺旋葉片桿配套的電機(jī)功率4 kW，轉(zhuǎn)速為30 r/min，每小時處理1.5 t的畜禽糞污量，螺旋桿葉片的受力情況如圖4所示。畜禽糞便對葉片的法向力Fn，將Fn分解為Fx和徑向力Fr，再將Fx分解為畜禽糞便對葉片的軸向力Fa和圓周力Ft[3-5]。

查閱相關(guān)資料[6-8]可查到畜禽糞便與螺旋葉片的動摩擦系數(shù)f=0.6，d為螺旋軸直徑，T為扭矩，T=9 549·P1/n1（N·m），P1為工作時的功率，n1為工作時的轉(zhuǎn)速，n1=30 r/min，正常工作時，轉(zhuǎn)換的功率為額定功率的98%，即P1=P×98%=4 kW×98%=3.92 kW，螺旋螺齒的推物料面傾角α為10°，背面的傾角β為27°。根據(jù)螺旋葉片的面積，可以求得At、Ar、Aa分別為30 558、30 558、251 327 mm2，通過計算可得螺旋的單位圓周力、單位軸向力、單位徑向力分別為0.778、12.657和0.888 MPa。

3.2? 隔板理論計算

該復(fù)合式發(fā)酵塔結(jié)構(gòu)呈長方體，以其中一層進(jìn)行分析，長6 000 mm，寬2 800 mm、高2 200 mm，計算出體積V=6 m×2.8 m×2.2 m=36.96 m3，在現(xiàn)場測得畜禽糞便的密度ρ為800 kg/m3，則該層的畜禽糞便質(zhì)量m=ρ·V=800 kg/m3×36.96 m3=2.956 8×104 kg，一層中由12個隔板連續(xù)排列組成，畜禽糞便在塔內(nèi)堆放，對隔板的壓力為靜壓力。因此，隔板承受的壓力p=F/A，其中，F(xiàn)為畜禽糞便的總質(zhì)量的壓力，F(xiàn)=m·g=2.956 8×104 kg×9.8 N/kg=2.897 6×105 N（g為重力加速度），利用Solidworks的測量功能可以測得A=10.46 m2，求得p 總=2.77×104 Pa，發(fā)酵塔內(nèi)有12個隔板共同承受壓力，所以單個隔板受到的壓力p=p 總/12=2 308 Pa。

4? ANSYS workbench受力分析

4.1? 網(wǎng)格劃分

通過Solidworks建立螺旋桿和隔板的三維模型，定義材料為45號鋼，主要力學(xué)性能[9]：彈性模量E=2.09×e11Pa，泊松比μ=0.269，密度ρ=7 850 kg/m3，抗拉強(qiáng)度為600 MPa，屈服點(diǎn)為355 MPa，選用Solid187完成網(wǎng)格劃分[2]，該模型比較大，所以選擇螺旋桿網(wǎng)格尺寸為30 mm，隔板網(wǎng)格尺寸為20 mm，網(wǎng)格劃分如圖5、圖6所示。

4.2? 定義邊界條件

螺旋桿用于傳輸物料，電機(jī)的旋轉(zhuǎn)帶動螺旋桿的旋轉(zhuǎn)，因此螺旋桿的兩端是固定的，螺旋桿軸做旋轉(zhuǎn)運(yùn)動，對螺旋桿兩端添加約束。如果僅對兩端添加約束，會導(dǎo)致兩端面圓周應(yīng)力過大，需對螺旋桿兩端面添加除Z軸旋轉(zhuǎn)方向外其余位置全約束。

隔板以水平方向靜置時，僅受到畜禽糞便的壓力，分析隔板受到畜禽糞便的壓力，隔板的中間軸固定在兩端，對隔板軸添加全約束。

4.3? 施加載荷及后處理

螺旋桿帶動物料運(yùn)動，物料對螺旋桿有反作用力，分解為物料對螺旋桿的圓周力和軸向力以及底徑圓錐上的徑向力。根據(jù)螺旋桿推動物料運(yùn)動時受到的單位圓周力，單位軸向力以及單位徑向力，分別施加相對應(yīng)的載荷，處理結(jié)果如圖7至圖10所示，其中圖10為總應(yīng)力。

除分析應(yīng)力外，還應(yīng)對螺旋桿進(jìn)行模態(tài)分析，螺旋桿在復(fù)合式發(fā)酵塔的底部隨著電機(jī)的運(yùn)轉(zhuǎn)而運(yùn)轉(zhuǎn)，如果發(fā)生共振，將會縮短螺旋桿軸和發(fā)酵塔的使用年限，存在一定的安全隱患。利用ANSYS Workbench的modal模塊，對螺旋桿施加同樣的約束，進(jìn)行6階模態(tài)分析[10]，處理結(jié)果如圖11所示。

隔板水平靜置，僅受到畜禽糞便對隔板兩上端面的壓力，對隔板兩上端面施加壓力，處理結(jié)果如圖12所示。

由圖10可知，螺旋桿的最小應(yīng)力為1.26 MPa，最大應(yīng)力為591.12 MPa，其中所受最大應(yīng)力存在于螺旋桿葉片端面與螺旋軸接口的位置上，原因是螺旋桿在傳輸物料時，兩端面與物料直接擠壓，受到的反作用力比較大，且圓周力、軸向力以及徑向力都匯聚在葉片與軸的接口處，使其應(yīng)力集中。綜合圖7、圖8、圖9可知，螺旋桿在除了兩端面的其他位置上受力均勻，螺旋軸受力最小，葉片受力比較大，這是因?yàn)槲锪想S著螺旋桿運(yùn)動時，葉片與物料摩擦，圓周力、軸向力、徑向力作用在葉片上，而軸受到的力比較小，這與實(shí)際情況相吻合。從圖11中可以看到，前六階固有頻率分別為2.486 9×e-4、372.91、671.00、809.18、854.28、871.29，該固有頻率值都比較大，因此不會發(fā)生共振現(xiàn)象。

本研究所用螺旋桿材料為45鋼，其材料的屈服強(qiáng)度為355 MPa，此情況下的最大應(yīng)力值超過了屈服應(yīng)力，因此螺旋葉片會產(chǎn)生永久變形，縮短螺旋桿使用壽命。

從圖12可以看到，隔板受到的力不大，但隔板右側(cè)上端面已經(jīng)發(fā)生嚴(yán)重變形，左上端面沒有變形，原因是右側(cè)上端面為斜邊，且非常薄，僅為4 mm，承受的力非常有限，隔板軸最上端受力較大，原因是隔板軸上方與畜禽糞便屬于線接觸，受到的壓力就大;物料畜禽糞便對隔板的壓力遠(yuǎn)遠(yuǎn)超出了右上端面能夠承受的力，在此情況下隔板也會產(chǎn)生永久變形，進(jìn)而縮短隔板的使用壽命，存在安全隱患。

5? 總結(jié)

通過理論計算分析，螺旋桿葉片承受的單位軸向力遠(yuǎn)遠(yuǎn)大于單位圓周力和單位徑向力，結(jié)合ANSYS Workbench分析結(jié)果，確定螺旋桿葉片兩端承受的力最大，螺旋軸承受的力相對較小，因此對螺旋桿進(jìn)行結(jié)構(gòu)優(yōu)化時，總體應(yīng)加強(qiáng)螺旋桿葉片的強(qiáng)度，特別加強(qiáng)葉片兩端的強(qiáng)度;通過對隔板的理論計算和ANSYS Workbench分析，確定隔板右上端面承受的力十分有限，針對隔板的優(yōu)化，應(yīng)進(jìn)一步增加隔板右上端面的厚度，同時應(yīng)加強(qiáng)強(qiáng)度。

參考文獻(xiàn)：

[1] 喬金亮.我國畜禽養(yǎng)殖每年產(chǎn)生38億噸畜禽糞便，有效處理率卻不到50%——治理養(yǎng)殖污染？教你幾招[J].新農(nóng)村（黑龍江），2016（20）：50-52.

[2] 朱文堅.機(jī)械設(shè)計[M].北京：高等教育出版社，2015.

[3] 曹? 冬，李慶忠.基于ANSYS Workbench的螺旋擠渣機(jī)螺旋的有限元分析[J].機(jī)械制造與自動化，2018，47（1）：136-138，151.

[4] 李? 妍.基于ANSYS 軟件的接觸問題分析及在工程中的應(yīng)用[D].長春：吉林大學(xué)，2004.

[5] 鄧全得，盧博友.基于ANSYS的農(nóng)業(yè)物料螺旋壓榨機(jī)螺旋強(qiáng)度分析[J].農(nóng)機(jī)化研究，2013（8）：58-60.

[6] 魏學(xué)鋒，蘇紅鳳，劉艷軍，等.基于ANSYS Workbench的大型擠壓造粒機(jī)雙螺桿芯軸有限元分析[J].橡塑技術(shù)與裝備，2014，40（21）：46-49.

[7] 查太東，楊? 萍.基于Ansys Workbench的固定支架優(yōu)化設(shè)計[J].煤礦機(jī)械，2012，33（2）：28-30.

[8] 趙經(jīng)文，王宏鈺.結(jié)構(gòu)有限元分析[M].北京：科學(xué)出版社，2001.

[9] 王迪迪，宋新旺，陳福洋，等.基于ANSYS Workbench傳動軸模態(tài)分析[J].煤炭技術(shù)，2015，34（9）：287-289.

[10] 唐? 鑫，鄧遠(yuǎn)超.基于ANSYS Workbench和SolidWorks整體硬質(zhì)合金絲錐的模態(tài)分析[J].工具技術(shù)，2011（4）：33-34.