基于DEM-FEM 耦合法的斷齒式對輥制粒機仿真研究

王昌德， 高愛民， 任 偉， 王紅松， 劉 龍

（甘肅農(nóng)業(yè)大學(xué)機電工程學(xué)院，甘肅 蘭州 730070）

0 引言

秸稈是我國農(nóng)業(yè)生產(chǎn)中極具利用價值的可再生能源[1]。我國作為一個農(nóng)業(yè)大國，秸稈資源非常豐富，但由于其松散且密度小等特性，導(dǎo)致難以運輸和保存，作物秸稈飼料化技術(shù)可以解決這一問題[2-4]。國內(nèi)主要的顆粒化加工設(shè)備有環(huán)模制粒機、平模制粒機和對輥制粒機等，對輥制粒機相較其他類型制粒機具有結(jié)構(gòu)簡單、成型性能好和制粒效率高等優(yōu)點。目前，國內(nèi)外許多研究人員針對對輥制粒機展開了3 個方面的研究。①利用有限元分析方法對壓輥在擠壓過程中的強度及對輥制粒機結(jié)構(gòu)進行分析。張孟輝等[5]利用ANSYS 軟件建立了壓輥與壓輥軸裝配的有限元模型，根據(jù)壓輥和壓輥軸過盈配合裝配下壓輥和壓輥軸的強度、變形量和接觸應(yīng)力對壓輥軸進行了優(yōu)化。寧延州等[6]利用SolidWorks 軟件仿真分析了對輥制粒機關(guān)鍵零部件的相對運動軌跡，驗證了對輥結(jié)構(gòu)的可行性。②通過設(shè)計研制對輥制粒機實物進行試驗驗證，并分析影響制粒產(chǎn)量的因素。朱德文等[7]針對傳統(tǒng)制粒機在成型過程中出現(xiàn)的問題，設(shè)計了空心對輥式顆粒成型機，并分析確定了對輥式顆粒成型機核心部分結(jié)構(gòu)參數(shù)，試驗表明，其結(jié)構(gòu)和參數(shù)設(shè)計合理，能夠滿足生物質(zhì)顆粒成型的條件。閆昌國等[8]設(shè)計研制了雙環(huán)模顆粒成型設(shè)備，通過設(shè)置環(huán)模的不同轉(zhuǎn)速來研究該因素對顆粒產(chǎn)量的影響，結(jié)果表明，制粒產(chǎn)量并不是隨著轉(zhuǎn)速的增加而一直增大。③對制粒機核心部件的耐磨性進行分析。蔣清海等[9]研究發(fā)現(xiàn)，壓輥的主要磨損機理為拋光磨損、磨粒磨損和疲勞磨損，不同工作部位的磨損機理有所不同，并根據(jù)其磨損機理分析提出改進建議，對提高制粒機穩(wěn)定性和延長其壽命有一定指導(dǎo)意義。

目前對輥制粒機仍然存在許多不足，如壓輥齒形設(shè)計不合理導(dǎo)致物料擠壓不充分，成型顆粒質(zhì)量差、易松散等問題。傳統(tǒng)的對輥制粒機均是連續(xù)直齒形結(jié)構(gòu)，本研究將連續(xù)直齒形結(jié)構(gòu)改為斷齒形結(jié)構(gòu)，為驗證斷齒形壓輥結(jié)構(gòu)的可行性，利用離散元方法，對斷齒與連齒形結(jié)構(gòu)壓輥的擠壓過程進行仿真分析，以期獲得一種成型效果好的對輥制粒設(shè)備，并對該方案壓輥的強度和應(yīng)力疲勞進行預(yù)測分析。

傳統(tǒng)分析壓輥所受應(yīng)力和疲勞壽命的方法是通過對壓輥施加均勻載荷，但壓輥在實際擠壓過程中的受力是不斷變化的，通過傳統(tǒng)分析方法得到的結(jié)果與實際過程存在一定的差異[10]。本研究采用離散元與有限元耦合的方法對兩種齒形壓輥的制粒效果及壓輥的強度與壽命進行分析，先用離散元軟件EDEM 對兩種齒形的壓輥進行擠壓過程模擬仿真，通過研究兩種壓輥在相同仿真時間點的擠壓狀態(tài)及物料顆粒的受力特性，再結(jié)合顆粒重疊量的變化來分析兩種壓輥的制粒效果；并將EDEM 中仿真得到的壓輥受力數(shù)據(jù)耦合到ANSYS Workbench 中研究擠壓過程中壓輥的應(yīng)力及變形特性，再將耦合分析得到的壓輥應(yīng)力數(shù)據(jù)導(dǎo)入nCode Designlife 軟件對斷齒形壓輥進行了應(yīng)力疲勞預(yù)測；驗證了斷齒形壓輥結(jié)構(gòu)的合理性，為對輥制粒機的設(shè)計制造提供參考。

1 對輥制粒機工作過程

1.1 結(jié)構(gòu)及飼料顆粒成型原理

對輥制粒機由給料系統(tǒng)、擠壓制粒系統(tǒng)等組成，其結(jié)構(gòu)如圖1 所示，當(dāng)物料向下墜落至擠壓區(qū)域時，通過兩壓輥之間的擠壓力來制成飼料顆粒[11]。為減少仿真時間，提高計算效率，簡化為3 排模孔排列的壓輥進行計算與分析。圖2 為連齒形對輥結(jié)構(gòu)和斷齒形對輥結(jié)構(gòu)示意圖，為便于分析，使兩種齒形壓輥結(jié)構(gòu)參數(shù)相同，其中壓輥內(nèi)徑、壓輥外徑、?？族F角、模孔有效長度和?？字睆椒謩e設(shè)置為200 mm、170 mm、45°、12 mm 和10 mm。

工作時壓輥在電機的帶動下順時針轉(zhuǎn)動，物料經(jīng)過調(diào)制后不斷地進入兩壓輥之間的擠壓區(qū)，在物料與壓輥之間摩擦力的作用下帶動定軸壓輥相對轉(zhuǎn)動。物料在進入兩壓輥輪處受到物料與壓輥之間產(chǎn)生的摩擦力和壓輥運轉(zhuǎn)時產(chǎn)生的相互擠壓力共同作用下，使物料不斷被強制擠壓入?？字行纬深w粒，并隨壓輥的旋轉(zhuǎn)而被不斷地從?？字袛D出，從壓輥內(nèi)部自然斷開后甩出，形成顆粒狀飼料，如圖3 所示。

圖3 對輥制粒機工作原理Fig.3 Working principle of roller granulator

1.2 物料與壓輥受力分析

1.2.1 物料受力分析

取制粒機在擠壓成型區(qū)域的一段物料進行受力分析（忽略物料重力的影響），如圖4 所示，A、B兩點分別為變形壓緊區(qū)與擠壓成型區(qū)的交界點，引A、B兩點的切線相交與C點。

圖4 物料受力分析Fig.4 Force analysis of material force

1.2.2 壓輥受力分析

在擠壓時，兩壓輥通過物料來傳遞壓力，根據(jù)作用與反作用力公理分析可知，兩壓輥表面的受力情況是對稱的，因此只需對其中一個壓輥進行分析即可。

當(dāng)?？兹肟谔幍臄D壓力為Fi時物料可被擠入?？?，則擠壓成型區(qū)內(nèi)壓輥表面各部位所受到的壓力都近似為Fi，因此在擠壓成型區(qū)內(nèi)壓輥表面各位置處的壓力基本相等，而在變形壓緊區(qū)，物料所受到的壓力逐步增大，其中在接近擠壓成型區(qū)的位置處，壓輥表面與擠壓成型區(qū)的擠壓力基本相等。此外，對輥制粒機在運轉(zhuǎn)過程中，物料堆積量也不能超出壓輥最高處，以堆積量最大為例進行分析，如圖5 所示。

圖5 壓輥受力分析Fig.5 Force analysis of pressure roller

根據(jù)數(shù)學(xué)幾何分析，擠壓成型區(qū)對應(yīng)的圓心角β=a/2，假設(shè)壓輥表面任意一點的圓心角為φ，則擠壓力F表達式為

2 對輥制粒機成型工作過程仿真分析

2.1 離散元接觸模型及EDEM 參數(shù)設(shè)定

2.1.1 離散元接觸模型

離散元模型選取較為常見的制粒材料玉米秸稈為制粒物料，其材料性質(zhì)是各向異性和不連續(xù)的，物料顆粒之間會產(chǎn)生較強的粘結(jié)力，常規(guī)的顆粒接觸模型不能表征這種粘結(jié)力，不利于準(zhǔn)確地模擬秸稈物料的流變性。因此，選用Hertz-Mindlin with JKR 接觸模型模擬物料顆粒之間的相互作用，該模型基于JKR 理論基礎(chǔ)建立，考慮了顆粒間的作用力對其塑性變形的影響。JKR 模型通過添加表面能來表征顆粒間接觸的作用力，較適用于秸稈這類顆粒間具有黏結(jié)力作用的顆粒模型模擬分析[12]。JKR 接觸模型的響應(yīng)方程為

式中R——接觸顆粒半徑，cm

R*——當(dāng)量顆粒半徑，cm

a——承受外載荷時的接觸半徑，cm

a0——外載荷為0 時的接觸半徑，cm

E——彈性模量，MPa

E*——當(dāng)量彈性模量，cm

P——施加的反向作用力，N

ω——表面能，N

δ——兩物體相互接觸的形變

F——兩物體相互接觸間的黏著力，N

K——有效彈性模量，N/m2

2.1.2 EDEM 基本參數(shù)設(shè)定

在EDEM 中需設(shè)置的材料參數(shù)主要有壓輥與秸稈物料顆粒2 部分，壓輥材料選用45#鋼，通過測量設(shè)置相關(guān)參數(shù)，如表1 和表2 所示。

表1 材料屬性參數(shù)Tab.1 Material property parameter

表2 接觸屬性參數(shù)Tab.2 Contact attribute parameter

將SolidWorks 中建立的三維模型導(dǎo)入EDEM 并設(shè)定模型零件的材料參數(shù)，設(shè)置沿Y軸負方向的重力加速度g，大小設(shè)為9.81 m/s2，采用動態(tài)生成的方式，仿真時間設(shè)置為2 s，仿真時間步長設(shè)置為30%，數(shù)據(jù)保存時間設(shè)為0.005 s，顆粒直徑設(shè)為1.8 mm，生成顆粒數(shù)為50 000，如圖6 所示，等顆粒落在擠壓區(qū)并完全穩(wěn)定后進入待攫取狀態(tài)。

圖6 對輥制粒機待攫取狀態(tài)Fig.6 Counterroll granulating machine ready to grab

2.2 仿真結(jié)果與分析

2.2.1 斷齒式與連齒式壓輥擠壓效果分析

在兩組仿真試驗中設(shè)定相同的仿真時間2 s，生成相同數(shù)量的顆粒，參考對輥制粒機實際工作狀態(tài)將壓輥轉(zhuǎn)速設(shè)置為5 rad/s，在相同的仿真條件下進行。連齒式與斷齒式壓輥在1.5 和2.0 s 時的擠壓狀態(tài)如圖7和圖8 所示，物料擠壓過程中，連齒式壓輥較斷齒式壓輥出現(xiàn)了嚴重的漏料現(xiàn)象，顆粒大多以松散狀態(tài)存在，不利于顆粒成型。

圖7 連齒式壓輥不同時間點擠壓狀態(tài)Fig.7 Extrusion state of toothed press roller at different time points

圖8 斷齒式壓輥不同時間點擠壓狀態(tài)Fig.8 Extrusion state of broken gear press roller at different time points

斷齒式壓輥對物料的擠壓成型效果明顯優(yōu)于連齒式壓輥，采用斷齒式壓輥設(shè)計可顯著提升顆粒的成型率，這是由于連齒式壓輥因軸向孔被連齒隔開導(dǎo)致擠壓區(qū)不連續(xù)且孔間隙處物料無法進入孔中而漏出，采用斷齒式壓輥設(shè)計后，使周向孔之間由周向排布的斷齒隔開，如圖2b 所示，當(dāng)兩個斷齒式壓輥嚙合后，兩壓輥的齒相互交錯擠壓以形成連續(xù)擠壓區(qū)，近似將變形壓緊區(qū)分為多個單獨的擠壓區(qū)域，每個區(qū)域由一個孔對應(yīng)一個齒以單獨擠壓，更有利于引導(dǎo)物料，使變形壓緊區(qū)的物料依次擠入模孔形成顆粒；有效避免了大量的物料因落在孔與孔之間的空隙上無法形成有效地擠壓而從兩壓輥間穿過。

圖9 和圖10 分別為斷齒式與連齒式壓輥在擠壓過程中顆粒重疊量和顆粒受力的變化曲線，由于齒形的不同，斷齒式壓輥中每個齒對應(yīng)一個?？祝梢猿浞值貙⑽锪蠑D壓進入?？字?，使顆粒之間更好地發(fā)生黏結(jié)，由圖9 和圖10 可知，斷齒式壓輥相較連齒式壓輥中物料顆粒所受的擠壓力和重疊量的峰值變化更大，這也說明斷齒式相對連齒式更有利于擠壓物料。

圖9 顆粒重疊量變化Fig.9 Variation of particle overlap

圖10 顆粒擠壓力變化Fig.10 Change in particle extrusion pressure

2.2.2 斷齒式對輥制粒機轉(zhuǎn)速分析

通過對斷齒式壓輥設(shè)置多組不同的轉(zhuǎn)速來計算各轉(zhuǎn)速下物料顆粒所受的擠壓力，結(jié)果表明，隨著壓輥轉(zhuǎn)速的增大，物料所受擠壓力也在隨之增大，當(dāng)轉(zhuǎn)速增大至20 rad/s 后，擠壓力增幅逐漸減??；為更直觀反映擠壓力增幅隨壓輥轉(zhuǎn)速的變化過程，從多組仿真試驗中選取轉(zhuǎn)速為5、10、15、20、25 和30 rad/s 6 組來表征這一結(jié)果，圖11 為6 種轉(zhuǎn)速下物料顆粒所受擠壓力隨時間的變化曲線，從曲線中得到各轉(zhuǎn)速下顆粒受到擠壓力的最大峰值及對應(yīng)的時間點如表3 所示，壓輥轉(zhuǎn)速在5～20 rad/s 范圍內(nèi)擠壓力增幅較大，壓輥轉(zhuǎn)速在20～30 rad/s 范圍內(nèi)擠壓力增幅逐漸減??；結(jié)合轉(zhuǎn)速在5～30 rad/s 范圍內(nèi)各轉(zhuǎn)速下擠壓狀態(tài)對比，轉(zhuǎn)速低于20 rad/s 時，由于轉(zhuǎn)速過慢而制粒效率較低，轉(zhuǎn)速超過25 rad/s 時，由于轉(zhuǎn)速過快，部分物料得不到有效擠壓導(dǎo)致擠壓區(qū)域漏料較嚴重。壓輥轉(zhuǎn)速范圍在20～25 rad/s 內(nèi)更有利于制粒，最大轉(zhuǎn)速為25 rad/s。

表3 壓輥在不同轉(zhuǎn)速下顆粒所受最大擠壓力及對應(yīng)的時間點Tab.3 The maximum extrusion pressure and corresponding time point of roller under different speed

圖11 不同轉(zhuǎn)速下顆粒擠壓力隨時間變化Fig.11 Particle extrusion pressure changes with time at different rotational speeds

2.2.3 斷齒式對輥制粒機在最大轉(zhuǎn)速下壓輥的應(yīng)力疲勞分析

通過離散元方法可以較準(zhǔn)確地得到物料顆粒對斷齒式壓輥的作用力，利用離散元和有限元耦合的方法，將離散元軟件EDEM 中計算得到的壓輥受力數(shù)據(jù)導(dǎo)入有限元軟件ANSYS Workbench 中，分析對輥制粒機的壓輥在工作狀態(tài)中由顆粒擠壓與摩擦所產(chǎn)生的應(yīng)力和應(yīng)變特性；將EDEM 與ANSYS 耦合分析得到的壓輥在擠壓過程中所受應(yīng)力的數(shù)據(jù)導(dǎo)入nCode Designlife 軟件，進而對壓輥的應(yīng)力疲勞情況進行預(yù)測分析。

在上述分析中得到對輥制粒機最佳轉(zhuǎn)速范圍為20～25 rad/s，考慮到制粒機對物料顆粒擠壓力的影響，在25 rad/s 時物料顆粒所受擠壓力最大，因而選擇壓輥轉(zhuǎn)速為25 rad/s 時的壓輥受力數(shù)據(jù)進行耦合分析，由圖11e 可知，在0～0.2 s 的時間段，物料的受力產(chǎn)生明顯的波峰，急劇增大，在0.12 s 時物料顆粒所受的擠壓力最大，因此，選取0.12 s 的數(shù)據(jù)分析壓輥的等效應(yīng)力和變形，從EDEM 結(jié)果后處理中將擠壓力數(shù)據(jù)導(dǎo)入ANSYS Workbench 中。導(dǎo)入壓輥的簡化模型，進行網(wǎng)格劃分，此時的模型網(wǎng)格尺寸及大小要與在離散元中分析模型的保持一致，以確保數(shù)據(jù)傳遞的準(zhǔn)確性，最后求解計算，得到壓輥的變形云圖和等效應(yīng)力云圖。

擠壓區(qū)域的物料經(jīng)過兩壓輥最小間隙處，一部分被擠壓進入?？?，一部分繼續(xù)繞著壓輥旋轉(zhuǎn)，由圖12所示，物料經(jīng)過兩壓輥最小間隙時壓輥的最大應(yīng)力為72.279 MPa，而壓輥材料45#鋼的屈服極限為355 MPa，遠大于壓輥所受的最大應(yīng)力；為防止壓輥在轉(zhuǎn)動時發(fā)生干澀，壓輥模型裝配時，兩壓輥間隙設(shè)為5 mm，由圖13 可知，其最大變形量為0.835 mm，遠小于許用變形量，由此說明，該壓輥方案的強度與剛度性能均滿足要求。

圖13 壓輥變形云圖Fig.13 Nephogram of press roll deformation

對輥制粒機在工作過程中，兩壓輥不斷受到相互擠壓力的作用，使壓輥周期性地受到疲勞彎曲應(yīng)力，導(dǎo)致壓輥發(fā)生疲勞破壞。根據(jù)對輥制粒機的實際工作情況，對輥制粒機屬于高周疲勞，故采用應(yīng)力疲勞分析方法。

建立壓輥材料的正應(yīng)力幅值與應(yīng)力循環(huán)次數(shù)曲線（S-N 曲線）；由于壓輥受到的載荷具有交變性，需要施加恒定幅值載荷譜以模擬壓輥在工作過程中受到的交變極限載荷；在nCode Designlife 軟件中根據(jù)壓輥材料參數(shù)來建立材料理論的S-N 曲線，經(jīng)抗拉極限修正，得到修正后的S-N 曲線如圖14 所示[13]。

圖14 壓輥材料的S-N 曲線Fig.14 S-N curve of press roller material

在nCode Designlife 軟件中對壓輥進行應(yīng)力疲勞分析，將材料、應(yīng)力和載荷參數(shù)設(shè)置完成后，分析得到壓輥的疲勞壽命云圖如圖15 所示。最大損傷位置發(fā)生在兩壓輥的最小間隙處，這是由于在擠壓區(qū)域物料不斷被擠壓進入?？祝谀？字車菀自斐蓱?yīng)力集中的現(xiàn)象，導(dǎo)致?？滋幐菀装l(fā)生疲勞破壞，?？滋幤趬勖蚕噍^其他部位更小，最小壽命為7.199e+12 次；在壓輥的兩壁面周圍由于物料無法被擠壓進入?？?，物料堆積導(dǎo)致此處的應(yīng)力過大，對壓輥壁面同樣造成疲勞破壞，相較?？撞课唬淦谄茐牧Χ雀。灰虼?，在對斷齒式對輥制粒機設(shè)計制造時需對壓輥表面模孔部分進行強化處理。

圖15 壓輥疲勞壽命云圖Fig.15 Fatigue life cloud of press roller

3 結(jié)束語

（1）斷齒式壓輥可以更充分地將物料擠壓入模孔，成型率高，并且斷齒式較連齒式壓輥在擠壓過程中顆粒受到的擠壓力和顆粒重疊量更大，由此，斷齒式壓輥的擠壓效果與制粒效率均優(yōu)于連齒式壓輥。

（2）物料受到來自壓輥的擠壓力隨著壓輥轉(zhuǎn)速的增加而增大，當(dāng)轉(zhuǎn)速由5 增至20 rad/s 時擠壓力的增幅變化較大，超過20 rad/s 時，擠壓力的增幅逐漸減?。唤Y(jié)合各轉(zhuǎn)速下模擬擠壓的效果對比，綜合得出最佳轉(zhuǎn)速范圍為20～25 rad/s。

（3）壓輥在工作狀態(tài)下持續(xù)受到物料顆粒的沖擊作用，產(chǎn)生的最大應(yīng)力主要集中在兩個壓輥的最小間隙處，最大應(yīng)力72.279 MPa 遠小于壓輥材料的屈服強度；壓輥在工作狀態(tài)下的最大變形量為0.835 mm，滿足剛度要求。nCode Designlife 軟件對壓輥的疲勞壽命分析表明，制粒機工作過程中在兩壓輥間隙最小時?？撞课坏膲勖钚?，也更容易發(fā)生疲勞損傷。