孫月銖，白雪衛(wèi)，李永奎

(沈陽(yáng)農(nóng)業(yè)大學(xué) 工程學(xué)院，沈陽(yáng)　110161)

?

玉米秸稈粉料單?？讐嚎s成型過(guò)程有限元分析

孫月銖，白雪衛(wèi)，李永奎

(沈陽(yáng)農(nóng)業(yè)大學(xué) 工程學(xué)院，沈陽(yáng)110161)

摘要：為了在有限元分析軟件中模擬玉米秸稈粉料，建立了秸稈粉料的D-P材料本構(gòu)模型，并通過(guò)三軸壓縮試驗(yàn)確定了秸稈粉料的粘聚力及內(nèi)摩擦角。在三維建模軟件中建立了秸稈粉料和單模孔成型模具的簡(jiǎn)化分析模型，利用有限元分析軟件中的熱應(yīng)力耦合分析來(lái)模擬秸稈粉料在模具模孔中被擠壓的全過(guò)程。根據(jù)分析得到的玉米秸稈粉料模型被擠壓后的變形及受力的規(guī)律，確定了秸稈粉料的變形特點(diǎn)、受力特點(diǎn)和時(shí)間-應(yīng)變關(guān)系，并進(jìn)行了試驗(yàn)驗(yàn)證，為模具的設(shè)計(jì)提供了依據(jù)。

關(guān)鍵詞：玉米秸稈粉料壓塊；有限元分析；D-P本構(gòu)模型；熱應(yīng)力耦合

0引言

我國(guó)是一個(gè)農(nóng)業(yè)生產(chǎn)大國(guó)，每年農(nóng)田玉米秸稈的產(chǎn)量約為2.2億t[1]。其中，用作還田的秸稈量?jī)H為15%，超過(guò)50%的秸稈被廢棄或就地燃燒[2-3]，既浪費(fèi)了大量生物質(zhì)能資源，又造成了環(huán)境的污染。玉米秸稈松散的形態(tài)及較低的容重使其在流通和貯存過(guò)程中成本增加。將秸稈粉碎后利用壓塊設(shè)備將玉米秸稈粉料壓制成高密度成型塊，其密度提高9～12倍，能大大減少運(yùn)輸成本和儲(chǔ)藏空間，便于秸稈資源商業(yè)化流通[4]。但是，目前秸稈的壓縮成型存在著很多的技術(shù)瓶頸，主要是對(duì)秸稈成型過(guò)程粘結(jié)機(jī)理的研究有待進(jìn)一步深入。

丁寧、孫勇[5]利用有限元分析軟件ANSYS對(duì)玉米秸稈粉料在平模中的成型過(guò)程進(jìn)行了二維模型的靜力分析。孫啟新[6]運(yùn)用非線性大變形分析對(duì)秸稈類生物質(zhì)成型徐變過(guò)程進(jìn)行有限元模擬分析。寧彭輝、于新奇[7]利用ANSYS，使用D-P材料本構(gòu)模型對(duì)秸稈壓縮成型過(guò)程進(jìn)行了分析。目前，對(duì)玉米秸稈擠壓過(guò)程的有限元分析主要集中在二維模型的靜力分析，且沒有考慮到溫度對(duì)成型過(guò)程的影響。鑒于此，本文在研究玉米秸稈粉料的本構(gòu)模型基礎(chǔ)上，利用有限元分析軟件對(duì)模具和玉米秸稈粉料進(jìn)行了三維建模及熱應(yīng)力耦合分析，來(lái)揭示秸稈粉料模型在壓制成型過(guò)程中的受力及應(yīng)變。

1玉米秸稈粉料本構(gòu)模型及參數(shù)測(cè)試

在分析中，玉米秸稈進(jìn)入?？浊?，經(jīng)粉碎機(jī)粉碎，其本構(gòu)模型(即應(yīng)力應(yīng)變關(guān)系)十分復(fù)雜[8]，單純的線彈性材料模型已經(jīng)無(wú)法模擬秸稈粉料，經(jīng)分析后采用了適用于混凝土、巖石、土壤等顆粒狀物料的Drucker-Prager模型。其屈服準(zhǔn)則是D-P準(zhǔn)則，是在米塞斯準(zhǔn)則的基礎(chǔ)上考慮平均主應(yīng)力對(duì)顆粒狀材料抗剪強(qiáng)度的影響而發(fā)展的一種準(zhǔn)則，是對(duì)Mohr-Coulomb準(zhǔn)則的近似。

M-C準(zhǔn)則在三維空間的屈服面為不規(guī)則的六角形截面的角錐體表面。引入應(yīng)力洛德角參數(shù)θσ，M-C屈服準(zhǔn)則可表示為

(1)

(2)

根據(jù)定義可知：θσ能夠反映一點(diǎn)的受力狀態(tài)的形式，即主應(yīng)力分量之間的比例關(guān)系。因而，不同的θσ可以反映材料的不同受力狀態(tài)。θσ= π/6時(shí)的M-C屈服準(zhǔn)則對(duì)應(yīng)的D-P準(zhǔn)則通常稱為外角點(diǎn)外接圓DP1準(zhǔn)則，如圖1 所示。

圖1　π平面內(nèi)的M-C準(zhǔn)則與D-P準(zhǔn)則屈服面

在國(guó)外的ANSYS、MARC、NASTRAN等有限元分析軟件中，D-P材料的屈服準(zhǔn)則通常選用的都是DP1準(zhǔn)則[10-11]。

在主應(yīng)力空間中，D-P屈服準(zhǔn)則表可表示成

(3)

其中，α、k為材料參數(shù)，是c、φ的函數(shù)[12]。

將θσ= π/6代入式(1)，與式(3)對(duì)比得

(4)

(5)

由式(4)和式(5)可知：建立材料D-P本構(gòu)模型時(shí)，必須設(shè)定材料的粘聚力c與內(nèi)摩擦角φ。本研究中，通過(guò)三軸壓縮試驗(yàn)[13]獲得粘聚力c與內(nèi)摩擦角φ。試驗(yàn)時(shí)，以含水率為18%的玉米秸稈粉料為原材料，將粒徑范圍如圖2所示的粉碎玉米秸稈制成3個(gè)φ39.1mm×80mm圓柱狀試樣，在(22±2)℃室溫下以1.5mm/min加載速率對(duì)試樣進(jìn)行軸向壓縮，待其軸向應(yīng)變達(dá)到15%時(shí)，按國(guó)家標(biāo)準(zhǔn)即可認(rèn)為達(dá)到破壞狀態(tài)。

圖2　粉碎玉米秸稈粉料粒徑級(jí)配

試驗(yàn)過(guò)程中，恒定圍壓，即小主應(yīng)力σ3分別設(shè)置為100、200、300kPa。三軸壓縮試驗(yàn)數(shù)據(jù)利用p-q法[14]計(jì)算得到含水率和粒徑范圍為該狀態(tài)下的秸稈粉粒的粘聚力為25.29kPa，內(nèi)摩擦角為27.76°。

2有限元分析

2.1材料參數(shù)設(shè)置

以細(xì)碎化處理后含水率為18%的玉米秸稈粉料為原材料，在(22±2)℃溫度下采用3次壓縮試驗(yàn)取均值的方法，試驗(yàn)數(shù)據(jù)處理后得出玉米秸稈粉料的基本參數(shù)如表1所示。另外，經(jīng)查閱，瞬態(tài)熱分析時(shí)需要設(shè)置的玉米秸稈粉料的比熱為1 620J/(kg·℃)[15]，玉米秸稈粉料的熱傳導(dǎo)率與水稻秸稈粉料近似，在22～100℃之間取熱傳導(dǎo)率數(shù)值為0.05～0.08W/(m·℃)[16]。模具材質(zhì)為結(jié)構(gòu)鋼，在材料庫(kù)中直接選取。

表1　玉米秸稈粉料參數(shù)設(shè)置

2.2分析模型的建立

在有限元分析中，非線性問(wèn)題主要包括狀態(tài)改變的接觸非線性行為，材料具有非線性的應(yīng)力-應(yīng)變關(guān)系，以及大變形導(dǎo)致幾何結(jié)構(gòu)明顯變化引起的非線性行為。在該模型中，3種類型的非線性同時(shí)發(fā)生，會(huì)占用大量的計(jì)算機(jī)資源。為了有效地進(jìn)行分析計(jì)算，得到相對(duì)可靠的結(jié)果，將模型進(jìn)行了適當(dāng)?shù)暮?jiǎn)化，如圖3所示。

圖3　簡(jiǎn)化后的模型

2.3分析前處理

2.3.1網(wǎng)格劃分

網(wǎng)格選擇的是自動(dòng)網(wǎng)格劃分，程序基于幾何體的復(fù)雜程度，自動(dòng)檢測(cè)實(shí)體，對(duì)可以掃掠的實(shí)體劃分六面體網(wǎng)格，對(duì)不能掃掠劃分的實(shí)體劃分四面體。在網(wǎng)格劃分之后，對(duì)物料與模具接觸錐面進(jìn)行了細(xì)化以提高計(jì)算精度。網(wǎng)格劃分結(jié)果如圖4所示。

圖4　模型有限元網(wǎng)格

2.3.2接觸設(shè)置

在玉米秸稈粉料的擠壓過(guò)程中，粉料和模孔會(huì)產(chǎn)生接觸，因?yàn)閮烧咧g存在著摩擦，所以玉米秸稈粉料與?？字g的接觸類型選擇的是面與面之間的摩擦接觸(Frictional)，摩擦因數(shù)取多次試驗(yàn)得到的平均值0.25。

2.3.3施加約束并求解

玉米秸稈主要由纖維素、半纖維和木質(zhì)素構(gòu)成,在高溫下秸稈中的木質(zhì)素和纖維素會(huì)軟化，軟化后的木質(zhì)素和纖維素具有粘結(jié)作用，由此秸稈物料在成型模孔內(nèi)被不斷壓緊、粘結(jié)。在試驗(yàn)過(guò)程中，通過(guò)溫控系統(tǒng)(YL-6SD，汕頭億隆電氣儀表有限公司)保證?？卓妆跍囟缺３衷?00℃。試驗(yàn)期間，實(shí)驗(yàn)室內(nèi)溫度為(22±2)℃，玉米秸稈在未進(jìn)入模孔時(shí)處于常溫狀態(tài)，玉米秸稈粉料進(jìn)入環(huán)模中進(jìn)行壓縮的過(guò)程中環(huán)模與秸稈粉料接觸部位會(huì)發(fā)生熱傳導(dǎo)。因?yàn)榻斩挿哿喜粩嗟剡M(jìn)入?？撞⑶已杆倜撾x?？?，所以最終選擇進(jìn)行瞬態(tài)熱分析。瞬態(tài)熱分析時(shí)，分析設(shè)置打開自動(dòng)時(shí)間步進(jìn)，設(shè)置子步時(shí)長(zhǎng)為0.01s。

在靜力分析中，設(shè)置模具外表面和下表面為固定支撐，在秸稈粉料模型上表面上施加沿Y軸向下16mm的強(qiáng)制位移，對(duì)秸稈粉料模型外表面施加遠(yuǎn)程位移，使物料不會(huì)沿著Y軸產(chǎn)生轉(zhuǎn)動(dòng)。施加重力加速度，因?yàn)樵趯?shí)際生產(chǎn)過(guò)程中物料會(huì)在重力的作用下產(chǎn)生一定的下移趨勢(shì)。

靜力學(xué)分析時(shí)，分析設(shè)置打開自動(dòng)時(shí)間步進(jìn)，由子步數(shù)目來(lái)間接控制時(shí)間步長(zhǎng)。由于在分析過(guò)程中秸稈粉料會(huì)產(chǎn)生位置移動(dòng)，所以必須打開大變形。

3結(jié)果與分析

3.1瞬態(tài)熱結(jié)果分析

?？字饕窍掳氩糠种慰妆诒３?00℃，瞬態(tài)熱分析玉米秸稈粉料模型的溫度變化主要集中模型的邊緣部分，在秸稈模型的縱切面靠近孔壁處向內(nèi)截取1.5mm×2.5mm的一部分模型放大16倍；之后取同一圖例，在該圖例下依次取t= 0.25T、0.5T、0.75T、T這4個(gè)時(shí)刻的溫度梯度，如圖5所示。

該溫度梯度明確顯示出：①由于秸稈粉料與模具的初始溫度不相同，在秸稈粉料與模具的接觸位置發(fā)生了熱傳導(dǎo)，秸稈粉料靠近模具的位置溫度升高，而其他不直接接觸的位置溫度變化較小，這與實(shí)際也是相符合的。②根據(jù)4個(gè)時(shí)刻的秸稈粉料模型的溫度變化可以明顯看出：在t= 0～ 0.25T時(shí)間范圍內(nèi)，秸稈粉料模型吸收的熱量較多，溫度變化速率較快；在t= 0.25T～T時(shí)間范圍內(nèi)，秸稈粉料模型吸收的熱量明顯減少。③在t=T時(shí)刻，玉米秸稈粉料模型溫度達(dá)到80℃的深度為0.45mm，最終擠壓出的物料表層硬化的深度在0～1.0mm之間。

圖5　玉米秸稈粉料的溫度梯度

3.2熱應(yīng)力耦合結(jié)果分析

熱應(yīng)力耦合分析選擇的是進(jìn)行順序耦合，是將瞬態(tài)熱分析求得的玉米秸稈粉料和模具的節(jié)點(diǎn)溫度作為體載荷施加到靜力分析中。

3.2.1秸稈粉料Y方向位移

在同一圖例下，t= 0.25T、0.5T、0.75T、T這4個(gè)時(shí)刻時(shí)秸稈粉料模型的Y方向位移變化情況如圖6所示。根據(jù)圖6可以分析得到：①在玉米秸稈粉料與?？椎慕佑|面位置存在著面與面之間的摩擦，所以在秸稈粉料的擠壓過(guò)程中接觸面會(huì)產(chǎn)生阻礙粉粒向下移動(dòng)的摩擦力；在摩擦力的作用下，接觸面的秸稈粉料的變形要比其內(nèi)部的物料滯后 。②因?yàn)榇怪庇阱F面的壓力大于其他位置，而內(nèi)摩擦因數(shù)各處都取實(shí)驗(yàn)均值0.25，所以在每個(gè)時(shí)刻錐面上受到的摩擦力都是最大的，導(dǎo)致的結(jié)果是靠近錐面的物料只微量移動(dòng)。③秸稈粉料在擠壓過(guò)程中，從模孔進(jìn)口至出口擠壓變形逐漸減小，同時(shí)物料的擠壓是分層擠壓，每層都是弧形的薄片，物料中間部分平鋪，越靠近?？自较蜇Q直方向扭曲，這在很大程度上會(huì)導(dǎo)致物料擠壓后的密度不均勻。

圖6　玉米秸稈粉料Y方向位移圖

3.2.2玉米秸稈粉料模型受力分析

玉米秸稈粉料擠壓過(guò)程中最大摩擦應(yīng)力隨位移變化曲線如圖7所示。從圖7中可以看出：在擠壓到12mm之前，最大摩擦應(yīng)力近似成線性增加；在壓縮秸稈粉料模型12mm之后，摩擦力增加幅度加快。經(jīng)過(guò)分析，物料模型上表面在Y方向向下的反力隨位移變化關(guān)系(見圖8)發(fā)現(xiàn)：在這個(gè)位移范圍反力的增幅也加大。繪制出應(yīng)變與時(shí)間的關(guān)系曲線，如圖9所示。

圖7　摩擦力-位移關(guān)系曲線

圖8　反力-位移關(guān)系曲線

圖9　應(yīng)變-時(shí)間(周期)關(guān)系曲線

應(yīng)變?cè)陂_始擠壓階段，隨時(shí)間的增加而逐漸增大，但是增加速率慢慢減??；擠壓12mm之后，應(yīng)變幾乎不再增大。這一現(xiàn)象產(chǎn)生的原因是秸稈粉料進(jìn)入衰減蠕變階段。在t= 0.75T之后，物料應(yīng)變?cè)黾虞^小，均勻下壓導(dǎo)致Y方向反力增加速率逐漸增大；又考慮到僅僅是少數(shù)區(qū)域進(jìn)入衰減蠕變階段，還有一些區(qū)域應(yīng)變未達(dá)到最大值，所以反力增幅緩慢加大，同時(shí)接觸正壓力增大引起摩擦應(yīng)力緩慢增加。

4驗(yàn)證與討論

圖10為實(shí)際成型試驗(yàn)壓制含水率為18%的秸稈粉料制成的秸稈顆粒。其截面形態(tài)沿壓縮方向，一端凸起、一端凹陷，與有限元分析中的秸稈粉料模型形態(tài)變化趨勢(shì)相一致。壓制成的顆粒表層硬化，表面顏色變深(變成深褐色)，硬化深度低于1.0mm。

圖10　成型試驗(yàn)壓制的顆粒

圖11為實(shí)際玉米秸稈粉料壓縮成型試驗(yàn)時(shí)得到的應(yīng)力-應(yīng)變關(guān)系曲線。該曲線明確顯示：在應(yīng)變低于0.35時(shí)，隨著應(yīng)變?cè)黾?，?yīng)力增加幅度較??；試驗(yàn)進(jìn)行過(guò)程中應(yīng)變達(dá)到0.4之后，需要較大應(yīng)力才能繼續(xù)產(chǎn)生變形，與有限元分析得到的數(shù)據(jù)變化趨勢(shì)接近，應(yīng)變數(shù)量級(jí)一致且數(shù)值接近。

圖11　試驗(yàn)加載過(guò)程應(yīng)力-應(yīng)變曲線

5結(jié)論

1)本文建立的玉米秸稈粉料的D-P本構(gòu)模型，比以往單純給定物料彈性模量、泊松比等參數(shù)建立的線彈性模型更有利于粉粒物料的模擬。

2)建立的玉米秸稈粉料和單?？啄＞叩臒釕?yīng)力耦合有限元模型，相比冷成型分析，更準(zhǔn)確地模擬了玉米秸稈粉料在?？字惺艿綌D壓成型的過(guò)程。

3)通過(guò)實(shí)際成型試驗(yàn)壓制的顆粒和應(yīng)力-應(yīng)變關(guān)系曲線驗(yàn)證，認(rèn)為D-P本構(gòu)模型與熱應(yīng)力耦合分析結(jié)合在一起，可以應(yīng)用于秸稈粉料壓縮成型的分析和研究中。

參考文獻(xiàn)：

[1]杜謀濤，袁曉東，郭和軍. 我國(guó)生物質(zhì)秸稈資源利用現(xiàn)狀及展望[J]. 能源與環(huán)境，2008(2)：76-78.

[2]石元春. 中國(guó)生物質(zhì)原料資源[J]. 中國(guó)工程科學(xué)，2011，13(2)：16-23.

[3]孟軍，張明偉，王紹斌，等. 農(nóng)林廢棄物碳化還田技術(shù)的發(fā)展與前景[J].沈陽(yáng)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2011，42(4)：387-392.

[4]Theerarattananoona K,Xu F,Wilson J,et al.Physical properties of pellets made from sorghum stalk, corn stover, wheat straw, and big bluestem[J]. Industrial Crops and Products，2011，33(2)：325-332.

[5]丁寧，邢碩燕，孫勇，等. 玉米秸稈致密成型燃料擠壓過(guò)程有限元分析[J]. 農(nóng)機(jī)化研究，2014，36(2)：48-51.

[6]孫啟新，張仁儉，董玉平. 基于ANSYS 的秸稈類生物質(zhì)冷成型仿真分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2009， 40(12)：131-134.

[7]寧鵬輝，于新奇.環(huán)模式秸稈壓塊機(jī)致密成型機(jī)理研究[D]. 天津：河北科技大學(xué)，2011.

[8]Faborode M O, O’CALLAGHAN J R. A Rheological Model for the Compaction of Fibrous Agricultural Materials[J]. Journal of Agricultural Engineering Research，1989，42：165-178.

[9]鄧楚鍵，何國(guó)杰，鄭穎人. 基于M-C準(zhǔn)則的D-P系列準(zhǔn)則在巖土工程中的應(yīng)用研究[J]. 巖土工程學(xué)報(bào)，2006，28 (6)：735-736.

[10]陳鵬，徐博侯. 基于φ取值區(qū)間的D-P系列準(zhǔn)則與M-C準(zhǔn)則匹配方法[J].力學(xué)季刊，2012，33(2)：270.

[11]趙尚毅，鄭穎人，劉明維. 基于Drucker-Prager準(zhǔn)則的邊坡安全系數(shù)定義及其轉(zhuǎn)換[J]. 巖石力學(xué)與工程學(xué)報(bào)，2006，25(S1)：2730.

[12]劉金龍，欒茂田，許成順，等.Drucker-Prager準(zhǔn)則參數(shù)特性分析[J].巖土工程學(xué)報(bào)，2006，25(S2)：4009-4011.

[13]歐陽(yáng)明松，徐連民. 某廢渣三軸試驗(yàn)與修正劍橋模型的仿真模擬[J].吉林農(nóng)業(yè)，2010(7)：48.

[14]陳力宏，唐松濤，張洪濤. 常規(guī)三軸試驗(yàn)數(shù)據(jù)處理的電子表格法[J].北京交通大學(xué)學(xué)報(bào)，2010，34(1)：54-57.

[15]李志合，易維明，劉煥衛(wèi). 陶瓷球固體熱載體與粉狀生物質(zhì)的熱平衡[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2009，40(7)：103-107.

[16]陸琳，陸方，羅永浩，等. 探針?lè)y(cè)量水稻秸稈熱導(dǎo)率[J].上海交通大學(xué)學(xué)報(bào)，2009，43(9)：1461-1464.

Finite Element Analysis of Corn Stalk Powder’s Extrusion Process Based on Single Die Hole Compression Mould

Sun Yuezhu， Bai Xuewei， Li Yongkui

(College of Engineering，Shenyang Agricultural University，Shenyang 110161，China)

Abstract：In order to simulate the corn stalk powder in the finite element analysis software , this article established the constitutive model of corn stalk powder which based on the D-P yield criteria , and the corn stalk powder’s cohesive force and it’s friction angle were got by triaxial compression tests. The simplified analysis model of corn stalk powder and single die hole mould was created in three-dimensional (3D) modeling software , and the thermal stress coupling analysis of this mould was done in the finite element analysis software to simulate the whole extrusion process of the corn stalk powder. According to the deformation and the force reaction and the fictional stress of this model , the deformation features and the mechanical characteristics of the corn stalk powder and the time-strain curve were determined and tested which provide foundations to the design of the mould.

Key words：corn stalk powder extrusion；finite element analysis；D-P constitutive model；thermal stress coupling

文章編號(hào)：1003-188X(2016)04-0237-06

中圖分類號(hào)：S216

文獻(xiàn)標(biāo)識(shí)碼：A

作者簡(jiǎn)介：孫月銖(1990-)，女，河北衡水人，碩士研究生，(E-mail)sunyz_1990@163.com。通訊作者：李永奎(1963-)，男，遼寧營(yíng)口人，教授，博士，(E-mail)sauyklee@163.com。

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(51405311)

收稿日期：2015-03-24