仇 義，陳 智，侯占峰，弭龍凱，邵志威，陳利杰

(內(nèi)蒙古農(nóng)業(yè)大學(xué) 機(jī)電工程學(xué)院，呼和浩特 010018)

0 引言

草原是我國重要的植被類型之一，也是重要的可更新資源和草地畜牧業(yè)基地。近年來，隨著牧區(qū)人口和社會(huì)經(jīng)濟(jì)的發(fā)展，對(duì)草原的利用強(qiáng)度日益加大，加之不合理的利用方式，使大面積的天然草原不斷退化，生產(chǎn)力下降[1-3]。內(nèi)蒙古自治區(qū)退化草原面積已占全區(qū)草原面積的31.77%，正趨向于成為一個(gè)不能自我維持且不可持續(xù)發(fā)展的系統(tǒng)[4-5]。因此，加快退化草地的治理、恢復(fù)草原生態(tài)系統(tǒng)是當(dāng)前十分緊迫的任務(wù)。

目前，退化草地補(bǔ)播是國內(nèi)外對(duì)退化草地植被恢復(fù)與重建的主要措施之一，但依靠人工補(bǔ)播牧草種子已經(jīng)不能滿足大面積草地植被恢復(fù)與重建的需要，只有利用飛機(jī)散播和噴播機(jī)噴播，才能最大限度地提高作業(yè)效率。不管是飛播還是噴播，都要求對(duì)種子進(jìn)行丸粒化包衣，以確保作業(yè)后種子的發(fā)芽率和成活率[6-9]。因此，研發(fā)牧草種子特別是小粒牧草種子丸?；略O(shè)備，對(duì)于采取工程措施恢復(fù)與重建退化草地植被，進(jìn)一步改善草原生態(tài)環(huán)境，實(shí)現(xiàn)畜牧業(yè)可持續(xù)發(fā)展具有十分重要的意義。

1 材料與方法

1.1 試驗(yàn)儀器設(shè)備

利用本課題自主設(shè)計(jì)研發(fā)冰草種子丸?；略囼?yàn)臺(tái)開展試驗(yàn)，如圖1所示。試驗(yàn)所需主要其他儀器設(shè)備主要包括電子天平(測(cè)量精度 0.1g)、測(cè)速儀、變頻器，以及用于種衣劑稀釋的燒杯、玻璃棒等。

1.粉料進(jìn)料口 2.進(jìn)料口閥門 3.稱重系統(tǒng) 4.步進(jìn)電機(jī) 5.粉料輸送管道 6.藥液管道 7.藥液儲(chǔ)存罐 8.高壓泵 9.電機(jī) 10.旋轉(zhuǎn)步進(jìn)電機(jī) 11.包衣鍋電機(jī) 12.出料口 13.電動(dòng)激振器 14.包衣鍋 15.噴頭 16.分流板 17.控制器 18.種子進(jìn)料口

1.2 試驗(yàn)材料

選取適宜在干旱半干旱地區(qū)自然氣候條件下生長的天然牧草種子(即冰草)為試驗(yàn)對(duì)象，種衣劑和藥粉作為試驗(yàn)材料。為了提高冰草種子包衣合格率、單籽率及單仔抗壓強(qiáng)度，選擇合適的包衣配方及包衣材料就顯得尤為重要。本課題組選用的包衣材料及試驗(yàn)組分[10]如下：①粉料 1.0份；②粘合劑 0.1～0.5份；③種衣劑水溶液 0.05～0.5份。其中，粉料組分及質(zhì)量百分比含量為：①硅藻土 25%～50%；②膨潤土 25%～50%；③滑石粉 10%～40%。粘合劑組分及百分比含量如下：羧甲基纖維素0.1%～10%，聚乙烯醇0.1%～5%，水適量。

1.3 試驗(yàn)響應(yīng)參數(shù)

以冰草種子包衣合格率J為冰草包衣作業(yè)質(zhì)量考核指標(biāo)，按中華人民共和國行業(yè)標(biāo)準(zhǔn)種子包衣機(jī)試驗(yàn)方法(JB/T 7730-2011)，以包衣劑包敷在冰草種子上的面積大于80%的冰草認(rèn)定為包衣合格，按照下式計(jì)算包衣合格率，即

(1)

其中，J為包衣合格率(%)；Zh為種衣劑包敷冰草面積大于80%的冰草粒數(shù)(粒)；Zb為包衣不合格冰草數(shù)(粒)。

2 數(shù)值模擬模型

2.1 接觸模型

為了更好地研究顆粒的運(yùn)動(dòng)軌跡，將包衣鍋內(nèi)的顆粒分成若干個(gè)相互獨(dú)立的單元，根據(jù)牛頓第二定律，循環(huán)迭代計(jì)算，確定每個(gè)時(shí)間步長內(nèi)每個(gè)單元受力及位移，進(jìn)而得到每個(gè)單元在每個(gè)時(shí)步下的平動(dòng)及轉(zhuǎn)動(dòng)關(guān)系；并實(shí)時(shí)更新單元位置，通過對(duì)每個(gè)單元進(jìn)行力、位移跟蹤計(jì)算，最終得到包衣鍋內(nèi)粒子群的宏觀運(yùn)動(dòng)規(guī)律[11-12]。在包衣鍋內(nèi)某時(shí)步下單個(gè)橢球顆粒單元i的平動(dòng)運(yùn)動(dòng)方程為

(2)

(3)

(4)

Ft=-Stδ

(5)

(6)

2.2 數(shù)值模擬

基于以上的顆粒運(yùn)動(dòng)學(xué)及接觸力學(xué)模型，本文采用離散單元法(Discrete Element Method, DEM)進(jìn)行數(shù)值模擬建模分析，從而可以有效地解決散粒體運(yùn)動(dòng)學(xué)及動(dòng)力學(xué)問題[13-15]。為了更好地模擬冰草種子與包衣藥粉彼此間的運(yùn)動(dòng)、碰撞及摩擦規(guī)律。利用Pro-E建立包衣鍋三維實(shí)體模型，并將該模型導(dǎo)入EDEM軟件中(見圖2)進(jìn)行仿真，在EDEM的后處理模塊中可得到每個(gè)顆粒的速度、運(yùn)動(dòng)軌跡及種粉在包衣過程中的顆粒的離散程度。在離散元EDEM軟件仿真Dynamics模塊下，將振動(dòng)引入到包衣鍋上。設(shè)定在振動(dòng)頻率為25Hz、振幅為3mm、包衣鍋轉(zhuǎn)速為40r/min，在仿真時(shí)長6s的條件下進(jìn)行仿真。采用Vector矢量型的表示方法，可以有效地得到冰草種子和包衣藥粉的下一時(shí)刻的運(yùn)動(dòng)及相應(yīng)的運(yùn)動(dòng)軌跡。同時(shí)，由矢量圖可以清楚地看到冰草種子和藥粉在振動(dòng)與旋轉(zhuǎn)的復(fù)合作用下能更充分地混合，進(jìn)而可達(dá)到較理想的丸化包衣效果。

從圖2可以看出：在振動(dòng)、旋轉(zhuǎn)及離心力的復(fù)合作用下，冰草種子和藥粉的分布主要集中在包衣鍋的右側(cè)。當(dāng)包衣鍋振動(dòng)頻率為25Hz、振幅為3mm、包衣鍋轉(zhuǎn)速為40 r/min時(shí)，種粉能達(dá)到理想的離散狀態(tài)且能夠充分接觸混合。在碰撞與摩擦力的作用下，種粉開始降速丸化包衣，達(dá)到理想的丸化包衣效果。因此，由仿真矢量圖可以清楚地看到冰草種子和藥粉每一時(shí)刻的時(shí)空分布、速度分布及冰草種子和藥粉的宏觀混合規(guī)律，進(jìn)而為后續(xù)分析處理提供可靠依據(jù)。

圖2 顆粒運(yùn)動(dòng)矢量圖Fig.2 Particle motion vector diagram

3 數(shù)值模擬過程及分析

3.1 單顆粒軌跡分析

為了明確顆粒的隨機(jī)運(yùn)動(dòng)，同時(shí)研究單粒冰草種子在包衣過程中與藥粉間是否具有跟隨性，以及是否冰草種子和藥粉具有一致的運(yùn)動(dòng)軌跡及趨勢(shì)。在散粒體仿真軟件EDEM后處理中的Selection模塊下的Grid Bin Group下對(duì)整個(gè)包衣鍋內(nèi)的粒子進(jìn)行網(wǎng)格劃分，分成若干個(gè)長方體區(qū)域，每個(gè)區(qū)域內(nèi)均有一定量的冰草種子和包衣藥粉，進(jìn)而分析每個(gè)小長方體區(qū)域內(nèi)種粉的運(yùn)動(dòng)及混合情況。同時(shí)，在離散元EDEM軟件仿真Dynamics模塊下，設(shè)定振動(dòng)頻率為25Hz，振幅為3mm，包衣鍋轉(zhuǎn)速為40r/min，在仿真總時(shí)長6s的條件下進(jìn)行仿真。選取其中一個(gè)長方體內(nèi)僅有1粒冰草種子和若干藥粉顆粒的小區(qū)域，在EDEM軟件中Manual Selection模塊下對(duì)冰草種子在整個(gè)仿真時(shí)長(6s)下進(jìn)行標(biāo)記。本文隨機(jī)選取某一時(shí)刻粒子群內(nèi)單顆粒在整個(gè)運(yùn)動(dòng)過程的軌跡圖，如圖3所示。

3.2 仿真結(jié)果與討論

從圖3可以看出：標(biāo)記的冰草種子與藥粉在包衣鍋中的某個(gè)限定的小區(qū)域內(nèi)，從產(chǎn)生的那一時(shí)刻開始，在包衣鍋振動(dòng)加旋轉(zhuǎn)的復(fù)合作用下與包衣鍋一起做圓周運(yùn)動(dòng)。此時(shí)，藥粉可在包衣鍋的某個(gè)區(qū)域內(nèi)一直跟隨著冰草種子運(yùn)動(dòng)；當(dāng)冰草種子和藥粉的喂入量逐漸增加的時(shí)候，冰草種子和藥粉的碰撞摩擦次數(shù)也在逐漸增多，促進(jìn)冰草種子和藥粉丸?；?。因此，在整個(gè)運(yùn)動(dòng)過程中包衣藥粉與冰草種子有足夠的時(shí)長進(jìn)行充分接觸，進(jìn)而提高冰草種子丸?；滦Ч鞍潞细衤?。

圖3 單顆粒運(yùn)動(dòng)軌跡Fig.3 Trajectory of single particle

4 正交驗(yàn)證

以包衣合格率為試驗(yàn)性能指標(biāo)，以包衣鍋轉(zhuǎn)速、包衣時(shí)間、振動(dòng)頻率、包衣鍋傾角為試驗(yàn)因素，每個(gè)因素分別取3個(gè)水平進(jìn)行正交試驗(yàn)[16-17]，因素水平編碼表如表1所示。

表1 因素水平編碼表Table 1 Factors and levels encode table

每組性能試驗(yàn)進(jìn)行3次，每次間隔時(shí)間分別為10、15、20min，試驗(yàn)結(jié)果取平均值，試驗(yàn)物料每次數(shù)量大于50 000粒。根據(jù)預(yù)試驗(yàn)得到的種衣劑配比要求進(jìn)行包衣，從每份試驗(yàn)樣品中分出200粒試樣用5倍放大鏡觀察每粒試樣，再分出包衣劑包敷的冰草面積大于80%的冰草認(rèn)定為包衣合格，將其包衣合格率記錄在表2中。

計(jì)算極差并確定因素的主次順序，即

(7)

其中，Rj稱為第j列的極差或其所在因素的極差。極差的計(jì)算結(jié)果見表2。

根據(jù)極差大小排出3個(gè)因素分布對(duì)3個(gè)指標(biāo)影響的重要性的主次順序?yàn)镃>A>B。

由正交試驗(yàn)可知，影響包衣合格率的主次因素分別為包衣鍋傾角>轉(zhuǎn)速>振頻。因此，通過正交試驗(yàn)可知，在振動(dòng)和旋轉(zhuǎn)的復(fù)合作用下，冰草種子具有較好的丸粒包衣效果，其包衣合格率可達(dá)98.2%。

表2 正交試驗(yàn)方案及試驗(yàn)結(jié)果Table 2 Orthogonal test scheme and test results

5 結(jié)論

1)利用Pro-E對(duì)包衣鍋進(jìn)行三維實(shí)體建模，同時(shí)采用三維仿真軟件EDEM模擬了冰草種子與藥粉的包衣過程和運(yùn)動(dòng)規(guī)律。從單粒種子運(yùn)動(dòng)軌跡曲線圖可知：包衣藥粉與牧草種子在丸?；逻^程具有跟隨性，種子與藥粉會(huì)隨包衣鍋?zhàn)鰣A周運(yùn)動(dòng)。

2)將振動(dòng)引入到冰草種子丸粒化包衣機(jī)內(nèi)，種粉能在不同的時(shí)間點(diǎn)、相同的空間區(qū)域分布得較為離散，達(dá)到理想的種粉混合效果，并可提高冰草種子的包衣合格率。

3)對(duì)影響丸?；碌膮?shù)進(jìn)行正交試驗(yàn)，確定了影響包衣合格率的主次因素為：包衣鍋傾角>包衣鍋轉(zhuǎn)速>包衣鍋振動(dòng)頻率。