李海英，趙榮煊，滕軍華，張　滔

(華北理工大學 冶金與能源學院，河北 唐山　063000)

?

秸稈環(huán)模式壓塊機模具優(yōu)化研究

李海英，趙榮煊，滕軍華，張?zhí)?/p>

(華北理工大學 冶金與能源學院，河北 唐山063000)

摘要：針對秸稈環(huán)模式壓塊機的主要構(gòu)件環(huán)模模具在使用中出現(xiàn)的磨損問題，對其進行優(yōu)化處理。根據(jù)模具的工作原理，應用UG和ANSYS軟件對模具進行建模與靜力學分析，得出秸稈與模具應力與應變模擬分布情況。依據(jù)模擬結(jié)果，對模具進行4種組合的倒角分析，得到了不同倒角時模具與秸稈的應力與應變分布云圖，分析了模具與秸稈接觸間應力和應變趨勢及與倒角變化之間的關(guān)系，并對優(yōu)化前與優(yōu)化后的模具進行疲勞分析，模擬出安全系數(shù)云圖。結(jié)果表明：模具在使用過程中容易出現(xiàn)應力、應變集中現(xiàn)象，集中區(qū)域容易造成模具損壞；模具倒角處理后，發(fā)現(xiàn)應力、應變集中區(qū)域隨著倒角半徑的增大而分散；當?shù)菇墙咏趫A形時模具優(yōu)化效果最好，模具整體安全系數(shù)提高，具有較好的實用性能，為環(huán)模結(jié)構(gòu)進一步優(yōu)化提供了依據(jù)。

關(guān)鍵詞：環(huán)模模具；倒角；疲勞分析；有限元模擬；秸稈

0引言

秸稈作為可再生的生物質(zhì)資源，取之不盡用之不竭；但作為燃料有密度小、質(zhì)地疏松、不便運輸?shù)热毕?，嚴重制約了生物質(zhì)能的健康發(fā)展。利用生物質(zhì)固化成型技術(shù)，可將廢棄的秸稈回收、加工制成具有一定密度、便于儲存和運輸?shù)纳镔|(zhì)固體燃料[1-2]。環(huán)模式秸稈壓塊機作為生物質(zhì)固化成型的主要設備，具有產(chǎn)量高、節(jié)能環(huán)保等優(yōu)點[3]。作為壓塊型設備關(guān)鍵部件的環(huán)模，其壽命是影響環(huán)模壓塊機推廣的關(guān)鍵因素。

秸稈原料與環(huán)模模具之間的相互作用力是影響環(huán)模使用壽命的重要因素，因此對秸稈與環(huán)模進行數(shù)值分析具有重要的意義[4]。為了提高壓塊機環(huán)模模具的使用壽命，分析模具在擠壓秸稈過程中的受力情況是非常必要。目前，主要的研究集中在對環(huán)模壓塊機的模具結(jié)構(gòu)分析，現(xiàn)有文獻對秸稈壓塊機模孔與倒角之間的受力變化關(guān)系、疲勞分析較少。為此，本文在研究壓塊機環(huán)模結(jié)構(gòu)分析的基礎上，研究環(huán)模模具在壓塊過程中與秸稈之間相互的受力，明確環(huán)模模具受力部位與受力分布情況，并對其優(yōu)化；運用UG和ANSYS軟件，研究模具在工作過程中的應力及應變分布，同時模擬模具不同倒角的應力應變情況，為提高環(huán)模使用壽命，優(yōu)化環(huán)模結(jié)構(gòu)提供一定的依據(jù)。

1壓塊機環(huán)模結(jié)構(gòu)及工作原理

唐山新順達農(nóng)機裝備有限公司生產(chǎn)的9Y-70型環(huán)模秸稈壓塊機由固定的單列方?？篆h(huán)模和2個嚙合的偏心壓輪組構(gòu)成。環(huán)模是由單個“工”字形模具拼接成，?？淄ǔ榉叫?，截面積的邊長一般約為40mm，環(huán)模一周有32個?？住簤K機環(huán)模如圖1、圖2所示，采用42CrMo合金鋼鍛造而成。

圖1　環(huán)模模具

圖2　環(huán)模示意圖

工作時，粉碎后的秸稈經(jīng)過喂料口進入環(huán)模與兩個偏心輪壓輪組成的成型腔中，隨著主軸旋轉(zhuǎn)帶動兩個壓輪的轉(zhuǎn)動；壓輪的轉(zhuǎn)動將秸稈推壓進環(huán)模腔和表面的溝槽中，由于兩個偏心壓輪沿著環(huán)模溝槽內(nèi)切公轉(zhuǎn)和摩擦自轉(zhuǎn)，將物料擠壓進環(huán)?？字校粌蓚€偏心壓輪每完成1個周期的公轉(zhuǎn)，就將環(huán)模溝槽內(nèi)的物料擠壓進?？變?nèi)，形成秸稈壓塊的一個簡單壓層；隨著偏心壓輪的公轉(zhuǎn)、自轉(zhuǎn)及秸稈的不間斷擠壓，接連地形成秸稈壓層，擠入模孔中，秸稈在模孔中不斷擠壓最終成型[5]。

2成型受力分析

秸稈在環(huán)模模孔內(nèi)的擠壓成型過程是環(huán)模壓塊機工作的主要過程，也是環(huán)模模具磨損的主要過程。其受力主要為秸稈在模具內(nèi)與?？椎哪Σ亮烷_口角形成的擠壓力，共同作用實現(xiàn)秸稈的壓縮變形[6]。根據(jù)作用力與反作用力間的關(guān)系，秸稈受到的推力主要是作用在偏心輪上，所以推力的大小直接影響著秸稈與模具的擠壓力、摩擦力的大??；而秸稈與?？妆诘哪Σ亮τ绊懼斩挼某尚图澳？妆诘哪p，對壓塊機環(huán)模的使用壽命影響很大。成型過程力學分析如圖3所示[7]。

P.擠壓力　FN.內(nèi)壁對原料正壓力　Ff.原料與內(nèi)壁摩擦力

由環(huán)模的截面圖分析可知：在水平方向上的受力為

P=Ffcosθ+FNsinθ

PnSmax=Ff(D2-d2)cosθ+2FN(D+d)Lsinθ

式中Pn—對原料的單位擠壓力；

Smax—?？走M口端面積；

ε—側(cè)壓系數(shù)；

μ—摩擦因數(shù)。

由式(4)可看出：秸稈性能參數(shù)為ε、μ，成型壓力主要克服摩擦力，當其它參數(shù)給定時，開口角度θ的大小是影響環(huán)模?？啄p的關(guān)鍵因素[8]。

3壓縮成型過程的有限元分析

3.1幾何建模

選用秸稈壓塊環(huán)模成型機作為研究對象，如圖4所示。

圖4　模具幾何模型

ANSYS要求導入的模型參數(shù)應與建模軟件中的單位統(tǒng)一，本文采用的單位都是mm。采用UG三維軟件進行建模，將其模型導入ANSYS有限元軟件中進行分析。由于環(huán)模屬于環(huán)狀，因此可以選擇一個環(huán)?？鬃鳛檠芯繉ο?；同時模孔又具有軸對稱的特征，關(guān)于兩個坐標面對稱，為了減少計算的時間，因此只需要取出1/4進行分析即可，如圖5所示。

圖5　簡化幾何模型

環(huán)模模具采用42CrMo，彈性模量為2.06GPa，泊松比為0.3，抗拉強度為6.2E+08Pa；秸稈的彈性模量1.6GPa，泊松比為0.3。根據(jù)實際情況，對模型適當簡化，模具簡化為剛形體，秸稈為塑性體，忽略其內(nèi)部的物理、化學反應。

網(wǎng)格劃分是有限元分析的必不可少的一部分，由于秸稈成型過程中產(chǎn)生大變形的塑性變形，因此對網(wǎng)格劃分的疏密程度要求較高。網(wǎng)格質(zhì)量的好壞直接影響到計算結(jié)果的精度和正確性，不合理的網(wǎng)格不僅導致計算時間過長、結(jié)果精度差，甚至可能導致無法求解。本文為了對節(jié)點的分布進行有效控制，采用對這兩個物體按照掃描方式劃分網(wǎng)格，對模具與秸稈的接觸表面進行網(wǎng)格細化，同時結(jié)合手動控制相關(guān)邊界處的單元長度。劃分網(wǎng)格總數(shù)為21 168，節(jié)點數(shù)為93 742，網(wǎng)格劃分效果如圖6所示。

(a)　結(jié)構(gòu)網(wǎng)格

(b)　網(wǎng)格質(zhì)量

3.2邊界條件及求解

對模型施加邊界條件，秸稈和?？捉佑|的各個面的接觸定義為摩擦接觸，秸稈物料和環(huán)模模具相互接觸、相互擠壓，在此過程中秸稈發(fā)生較大的變形，很明顯這是一對剛?cè)峤佑|[9-10]。

本文主要是模擬秸稈成型過程中環(huán)模模具應力、應變分布，環(huán)模模具為剛性體作為目標面，秸稈模型為柔性體作為接觸面，摩擦因數(shù)取0.2，建立接觸對。秸稈物料被擠壓時發(fā)生很大的幾何變形和塑性變形，環(huán)模模具發(fā)生忽略不計的彈性變形如圖6所示。因此，根據(jù)ANSYS中的規(guī)定，秸稈與模具不能發(fā)生互相滲透，可傳遞法向壓縮力和切向摩擦力，但不能傳遞法向拉伸力。

如圖7所示，接觸的細節(jié)面板設置如下：

1)說明接觸類型是帶摩擦的接觸，摩擦因數(shù)是0.2，是非對稱接觸；

2)指明法向接觸面的剛度因子是0.1。

圖7接觸設置

Fig.7Contact establishment

由于本例是非線性問題，在大變形問題的摩擦接觸中，使用增強拉格朗日法(Augmented Lagrange)，這是因為增強拉格朗日公式增加了額外的控制自動減少滲透功能。

如圖8所示，進行求解設置如下：

1)意味著只有一個載荷步，該載荷步也只有一個載荷子步，關(guān)閉了自動時間步長，該載荷步結(jié)束的時間是100步；

2)是打開大變形開關(guān)[11]。

圖8　求解設置

3.3結(jié)果與分析

3.3.1應力應變分析

由模擬結(jié)果可以看出：在環(huán)模的應力場分布云圖中最大應力出現(xiàn)在環(huán)模模具的拐角處，?？鬃畲髴?04.66MPa。根據(jù)所選環(huán)模材料的物理參數(shù)可知：環(huán)模所受最大應力遠小于其材料的彈性極限，故壓縮成型過程中環(huán)模是比較安全的。同時，出現(xiàn)了應力集中區(qū)域，也就是這個區(qū)域最容易發(fā)生磨損、破壞。從秸稈的運動軌跡來看，秸稈邊緣受到模具接觸而產(chǎn)生摩擦力影響速度較慢。因此，在成型塊的截面上存在剪切應力，由圖9可以看出：模具最大應力出現(xiàn)在?？捉孛婀战翘?有限模擬與理論分析相吻合。應力沿軸向不均勻分布且越與進料口接近所受到的應力越大，導致環(huán)?？椎墓战翘幰壮霈F(xiàn)磨損，應力集中對構(gòu)件的疲勞壽命影響很大，使模具產(chǎn)生疲勞裂紋。

圖9　應力場模擬

通過對環(huán)模模具部分進行ANSYS分析可知：在壓塊過程中，環(huán)模受到秸稈的擠壓和摩擦力的雙重作用， 會產(chǎn)生微量的應變。應變量的大小影響模具的磨損程度，模擬結(jié)果如圖10所示。模具應變是最大的為0.000 566 03mm/mm。應變云圖和應力云圖相近，應力、應變集中區(qū)域都是在模具的拐角處，進一步證明了這個區(qū)域是環(huán)模模具磨損最大的區(qū)域。

圖10　應變場模擬

3.3.2模型優(yōu)化

根據(jù)對模具的實際運行狀況的模擬，秸稈在成型過程中產(chǎn)生較大的擠壓力，由于作用力與反作用力，模具同時也受到較大的擠壓力。綜上分析可得：為減輕秸稈成型過程中環(huán)模?？桌饨翘帒鞋F(xiàn)象，對?？捉Y(jié)構(gòu)進行改進。通過閱讀與環(huán)模數(shù)值模擬相關(guān)的國內(nèi)外文獻，提出對環(huán)模的改進模擬。對?？桌饨沁M行倒圓處理，對模具應力和應變集中的區(qū)域進行優(yōu)化，對參數(shù)進行設置并求解。該過程與未進行倒角處理的模型完全相同，圖11～圖14所示為應力場模擬。

圖11　倒角半徑為5mm應力場模擬

圖12　倒角半徑為10mm應力場模擬

圖13　倒角半徑為15mm應力場模擬

圖14　倒角半徑為20 mm應力場模擬

通過應力云圖分析表明：應力沿軸向不均勻分布且越與進料口接近所受到的應力越大；模具的最大應力一般集中在環(huán)模孔的拐角處，將會造成該處出現(xiàn)磨損。應力集中對構(gòu)件的疲勞壽命影響很大，使模具產(chǎn)生疲勞裂紋。未經(jīng)過優(yōu)化過的?？?，倒角半徑為0mm應比較集中，最大應力為104.66MPa。隨著倒角半徑的逐漸增大，從5、10、15、20mm應力69集中現(xiàn)象分散，最大應力也在同時減少。

圖15～圖18為應變場模擬圖。通過對模塊的環(huán)模部分進行ANSYS分析，在不同加載倒角環(huán)模的受力都有所減小，可知在壓塊過程中環(huán)模受到物料的擠壓和摩擦力的作用下，會有應變產(chǎn)生。

圖15　倒角半徑為5mm應變場模擬

圖16　倒角半徑為10mm應變場模擬

圖17　倒角半徑為15mm應變場模擬

圖18　倒角半徑為20mm應變場模擬

應變量的大小影響模具的磨損程度。模擬結(jié)果表明：未經(jīng)過優(yōu)化處理過的模具，應變是最大的為0.000 566 03mm/mm，隨著倒角半徑的增加，應變先是逐漸減小，最小的應變范圍出現(xiàn)在倒角為20mm，接近于圓形。

對環(huán)模成型過程中應力有限元模擬，結(jié)果表明：秸稈致密成型過程對環(huán)模最大的應力集中產(chǎn)生于棱角處，而對?？讓嵤├饨堑箞A，應力集中情況得到明顯改善；圓角的半徑越大，最大與最小的應力值指差最小，最小值為22.786MPa。另外，對環(huán)模模具進行了應變分析，結(jié)果表明：由于秸稈成型過程中與模孔壁面擠壓，環(huán)模產(chǎn)生應力形變，應力最大同時又是應力集中的部位形變量是最大的；就局部和整體而言，倒角的環(huán)模應變場分布較為平緩，最大應變量為0.000 118 63mm/mm，比未倒角的環(huán)模形變量減少0.000 447 4mm/mm，因此更能延長環(huán)模的使用壽命。

根據(jù)上述模擬得出表1，并繪制出擬合曲線如圖19所示。其應力變化隨著倒角半徑的增加而減少，而應變與倒角之間不成線性關(guān)系。當?shù)菇前霃椒秶鸀?～13.5時，物料接觸面上的應力應變分布更為均勻合理，可以有效緩解?？桌饨翘幍膽鞋F(xiàn)象。比較優(yōu)化后與優(yōu)化前的模具可以發(fā)現(xiàn)，優(yōu)化后的模具在與秸稈物料接觸的整個接觸面上，應變分散同時應變量也在減小，應力分布更為均勻合理。因此，倒圓角是改善環(huán)模模孔應力集中的一種行而有效的方法。

表1　應力、應變關(guān)系

圖19　應力與應變曲線

3.4疲勞分析

當模具受到多次重復變化的載荷作用后，在應力值始終沒有超過材料的強度極限，甚至比彈性極限還低的情況下就有可能發(fā)生破壞。秸稈對模具應力波動的每個周期都會或多或少地損壞物體，在循環(huán)一定數(shù)量的周期之后，模具會變得越來越衰弱。

根據(jù)模擬得出來的結(jié)果，對其環(huán)模進行疲勞分析，如圖20和圖21所示。

圖20　未優(yōu)化模具安全系數(shù)圖

圖21　優(yōu)化后模具安全系數(shù)圖

由此可知：優(yōu)化后的明顯比優(yōu)化前的模具安全系數(shù)高；優(yōu)化前模具出現(xiàn)了非常低的安全系數(shù)SF值，該值出現(xiàn)在模具的拐角處區(qū)域，說明該區(qū)域在工作中極易出現(xiàn)磨損，從而導致模具失效；優(yōu)化后模具安全系數(shù)SF值在1～4之間，很明顯在拐角處無較低安全系數(shù)SF值，相對優(yōu)化前的模具安全系數(shù)有了較大的提高。

4結(jié)論

分析了秸稈壓塊環(huán)模機的成型過程，研究了環(huán)模的受力，并應用有限元分析軟件ANSYS數(shù)值模擬了成型過程中環(huán)模的應力應變場。結(jié)果表明：環(huán)?？讬M截面倒角對秸稈壓塊機?？资褂脡勖绊戄^大，且拐角處應力最大。為了減少應力集中現(xiàn)象，對?？讓嵤├饨堑箞A處理。結(jié)果表明：倒角接近于圓形時，經(jīng)過優(yōu)化后的模具，棱角處的應力不僅減小而且分散，集中情況得到明顯改善，可以有效緩解?？桌饨翘幍膽鞋F(xiàn)象，為以后改善環(huán)模的磨損、延長環(huán)模的使用壽命提供理論依據(jù)。

參考文獻：

[1]韓芳．我國可再生能源發(fā)展現(xiàn)狀和前景展望[J].可再生能源，2010，13(4)：23-25.

[2]楊俊成. 秸稈壓餅技術(shù)與應用[J].農(nóng)牧產(chǎn)品開發(fā)，1998(4)：22-23.

[3]歐陽雙平，侯書林，趙立欣，等.生物質(zhì)固體成型燃料環(huán)模成型技術(shù)研究進展[J].可再生能源，2011,29(1):14-18.

[4]霍麗麗，侯書林，田宜水，等.生物質(zhì)固體燃料成型機壓輥磨損失效分析[J]．農(nóng)業(yè)工程學報，2010，26(7):102-106．

[5]歐陽雙平，侯書林，趙立欣，等．生物質(zhì)固體成型燃料環(huán)模成型技術(shù)研究進展[J]．可再生能源，2011,29(1) :14-18，22．

[6]王誠輝，嚴永林．木質(zhì)顆粒成型機環(huán)模受力分析[J]．中南林業(yè)科技大學學報，2011，31(12):187-191．

[7]王慧.基于生物質(zhì)碾壓成型機理的成型能耗影響因素研究[D].濟南:山東大學，2011．

[8]焦安勇.基于有限元模擬分析的生物質(zhì)壓縮成型機的研發(fā)[D].長春：吉林大學，2009.

[9]武凱，施水娟，彭斌彬，等.環(huán)模制粒擠壓過程力學建模及影響因素分析[J]．農(nóng)業(yè)工程學報，2010，26(12):142-147．

[10]施水娟，武凱，蔣愛軍．制粒過程中環(huán)模力學模型的建立及有限元分析[J]．機械設計與制造，2011(1):38-40．

[11]袁越錦，徐英英，張艷華.ANSYS Workbench 14.5建模與仿真從入門到精通[M].北京：電子工業(yè)出版社.2014:384-390.

Design of the Divided No-till Wheat Planter Stalks Ring Pattern Briquetting Machine Mold Bevel Edge Optimization Research

Li Haiying, Zhao Rongxuan, Teng Junhua, Zhang Tao

(College of Metallurgy and Energy, North China University of Science and Technology, Tangshan 063000,China)

Abstract：In view of the straw stalk ring pattern briquetting machine primary member ring mold, the attrition question which appears in the use, carries on optimized processing to it. According to the mold principle of work, carries on the modelling and the statics analysis using UG and the ANSYS software to the mold, obtains the straw stalk and the mold stress and the strain simulation distributed situation. Based on the analogue result, carries on 4 kind of combinations to the mold the bevel edge analyses, obtained when the different bevel edge mold and straw stalk stress and strain distributed cloud chart. Has analyzed the mold and the straw stalk contact the stress and the strain tendency as well as with the bevel edge change between relations before, and to optimizes carries on the fatigue analysis after the optimized mold, simulates the safety coefficient cloud chart. The result indicated that, the mold easy to appear the stress, the strain centralism phenomenon in the use process, the centralized region easy to create the mold damage.After mold bevel edge processing, discovered the stress, the strain concentrate the region to increase along with the bevel edge radius disperse, when the bevel edge approaches in the circular the mold optimizes the effect to be best. The mold whole safety coefficient enhances, has the good practical performance, carried on the structure for the link mold further to optimize has provided the basis.

Key words：ring mold; chamfer angle; fatigue analysis; finite element modeling; stalks

文章編號：1003-188X(2016)05-0065-07

中圖分類號：S817.11+5；TK6

文獻標識碼：A

作者簡介：李海英(1971-),女，河北唐山人，教授，碩士生導師，(E-mail)zhguijie@vip.sina.com。通訊作者：趙榮煊(1988-)，男，河北秦皇島人，碩士研究生，(E-mail)291462226@qq.com。

基金項目：河北省科技廳計劃項目(14227309D)

收稿日期：2015-06-29