紫花苜蓿圓捆機卷壓過程的應(yīng)力松弛特性與模型建立

房佳佳，張 永，楊明韶，王砹敏，王 健，劉德環(huán)，高 津，李海軍

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

0 引 言

牧草具有粘彈性，在加工過程中會產(chǎn)生應(yīng)力松弛現(xiàn)象，研究應(yīng)力松弛現(xiàn)象對降低圓捆機功耗、提高草捆質(zhì)量具有重要現(xiàn)實意義[1-3]。國內(nèi)外學(xué)者曾以方捆機為研究對象，針對牧草、農(nóng)作物秸稈等進(jìn)行了許多開式和閉式壓縮過程中一般工程力學(xué)和流變學(xué)研究，為壓縮工程優(yōu)化設(shè)計提供理論基礎(chǔ)[4-7]。圓捆機和方捆機壓縮方式不同，為滾卷壓縮，纖維物料在纏繞的壓縮過程中不斷滾動。壓縮方式、物料狀態(tài)不同，物料流變特性也不同[8-12]。近幾年，國內(nèi)學(xué)者以圓捆機為研究對象，進(jìn)行了草物料卷壓過程中流變學(xué)試驗研究。2012年，李汝莘等[13]在后倉門拉環(huán)處布置了應(yīng)變片，通過受力分析得到了草捆拉環(huán)方向受到的徑向力，根據(jù)卷壓過程中碎玉米秸稈的流變學(xué)特征，構(gòu)建了伯格斯模型；2015年，雷軍樂等[14]通過在捆繩機構(gòu)所處位置布置壓力傳感器的方法，測取了捆室上半部分鋼輥對草捆的徑向力，進(jìn)行了完整稻稈應(yīng)力松弛試驗研究，以應(yīng)力松弛時間和平衡彈性模量為指標(biāo)，進(jìn)行了試驗參數(shù)優(yōu)化。以上研究多是以小型圓捆機（直徑在1 m以下）為研究對象，針對玉米秸稈及稻稈進(jìn)行流變學(xué)試驗研究，針對牧草卷壓過程中流變行為的研究還未見報道。

另外，市場上的圓捆機一般為大圓捆機，生產(chǎn)干草草捆直徑為1.2～1.8 m，草捆密度為100～120 kg/m3，草捆質(zhì)量150～400 kg，喂入量5～10 t/h。草物料本身的物理特、捆室大小都會影響草物料的流變特性[15-16]。因此，本文針對大圓捆機，選取典型牧草——紫花苜蓿進(jìn)行卷壓試驗，研究牧草卷壓過程的應(yīng)力松弛行為，建立應(yīng)力松弛模型，獲取流變參數(shù)及本構(gòu)方程，為優(yōu)化打捆工藝、降低打捆功耗、提高草捆質(zhì)量提供理論基礎(chǔ)及技術(shù)支持。

1 材料與方法

1.1 試驗材料

紫花苜蓿是世界上栽培面積最廣、最主要的豆科牧草之一，有“牧草之王”之稱。中國苜?？偡N植面積約為133萬hm2，不僅可作為優(yōu)質(zhì)飼料，而且能夠固氮改土、改善生態(tài)環(huán)境建設(shè)[17-18]。

本文試驗采用內(nèi)蒙古呼和浩特市中國農(nóng)業(yè)科學(xué)院草原研究所試驗田八月份下旬收獲的紫花苜蓿為試驗物料。購入紫花苜蓿3 t，經(jīng)測定其含水率為30%，平均株高 80 cm。試驗時采用自然干燥法將苜蓿含水率調(diào)至20%，整株進(jìn)行試驗。含水率計算公式如下

式中Hc為牧草含水率，%；Gsc為牧草濕質(zhì)量，g；Ggc為牧草質(zhì)量，g。

1.2 試驗設(shè)備

試驗用呼倫貝爾市蒙拓農(nóng)機科技股份有限公司生產(chǎn)的 9YG-1.3型鋼輥外卷式圓捆機，所打草捆尺寸為φ1 300 mm×1 505 mm。圓捆機的總體結(jié)構(gòu)尺寸為2 610 mm× 2 330 mm×2 610 mm，主要由撿拾器、卷捆室、捆繩機構(gòu)、液壓系統(tǒng)和傳動系統(tǒng)組成。其中卷捆室尺寸固定，分為前后 2個部分，中間鉸接，后半部分稱為后倉門。14個鋼制輥子（直徑為300 mm）沿卷捆室周向布置，鋼輥的外表面有12條凸起的棱，有利于提高鋼輥對草物料的摩擦力。圓捆機工作時由鏈輪驅(qū)動各個鋼輥轉(zhuǎn)動，在鋼輥對草物料摩擦力和擠壓力作用下完成卷捆過程。圓捆機配套動力為30 kW電機，鋼輥額定轉(zhuǎn)速126 r/min，生產(chǎn)牧草草捆質(zhì)量小于300 kg，生產(chǎn)效率為12～30捆/h。

以蒙拓 9YG-1.3型鋼輥外卷式圓捆機為試驗設(shè)備，配置圓捆機支撐框架、皮帶輸送機、變頻器等組成圓捆機卷壓試驗平臺。圓捆機試驗平臺結(jié)構(gòu)示意圖如圖 1所示。皮帶輸送機的總體結(jié)構(gòu)尺寸為6 000 mm×2 130 mm×500 mm。用2臺電動機（型號為YVP225M-6、YE2112M-4）分別驅(qū)動圓捆機、輸送機。另配 2臺變頻器（型號為ZK1800G-37 kWG/45 KWP、EN600-4T0055P）分別控制鋼輥轉(zhuǎn)速、牧草喂入速度。

圖1 圓捆機卷壓試驗平臺Fig.1 Illustration of rotary compression test platform for round baler

試驗時，帶式輸送機將牧草以一定的速度輸送至圓捆機撿拾器下方，經(jīng)撿拾器彈齒的撿拾、喂入叉的導(dǎo)送進(jìn)入卷捆室。牧草在各鋼輥作用下在卷捆室內(nèi)不斷上升和下落，逐漸形成旋轉(zhuǎn)草芯。起初鋼輥對牧草只有摩擦力，為填充階段。隨著牧草不斷喂入，牧草填滿卷捆室，鋼輥開始對草捆產(chǎn)生徑向擠壓力，進(jìn)入卷壓階段，草捆密度快速增大。當(dāng)草捆密度達(dá)到設(shè)定值時，操作員停止喂入牧草，控制捆繩機構(gòu)動作進(jìn)行捆繩，當(dāng)捆繩圈數(shù)達(dá)到設(shè)定值時，切刀動作將捆繩切斷，液壓系統(tǒng)控制后倉門開啟放出草捆，形成外緊內(nèi)松的圓柱形草捆。

1.3 試驗設(shè)置

試驗用的紫花苜蓿含水率為 20%，草條密度為3 kg/m2[19]，鋪設(shè)紫花苜蓿幅寬為1 m，喂入紫花苜?？傎|(zhì)量為300 kg，喂入速度為0.83 m/s，鋼輥轉(zhuǎn)速為126 r/min。試驗時，將紫花苜蓿均勻喂入卷捆室，記錄鋼輥變形數(shù)據(jù)。當(dāng)草物料充滿捆室，鋼輥開始有變形讀數(shù)出現(xiàn)，草捆開始在鋼輥徑向力作用下被擠壓變形。隨后，每隔0.5 s測試 1次鋼輥變形數(shù)據(jù)。喂料結(jié)束后，草捆繼續(xù)在卷捆室內(nèi)運轉(zhuǎn)170 s（該時間由預(yù)試驗分析結(jié)果及相關(guān)文獻(xiàn)[20]確定），以測定其應(yīng)力松弛狀態(tài)，試驗重復(fù)3次，結(jié)果取平均值。

1.4 數(shù)據(jù)采集及處理

1.4.1 數(shù)據(jù)采集系統(tǒng)

由于鋼輥一直處于運動狀態(tài)且草捆對鋼輥的徑向反作用力非常大，采用接觸測量時有許多困難，如壓力傳感器布置、測量頭容易磨損、接觸狀態(tài)、信號難以采集等問題，容易造成測量誤差。因此，本文選用無線電阻應(yīng)變儀，獲取鋼輥受到的徑向力，即鋼輥對草捆的徑向作用力。

鋼棍結(jié)構(gòu)如圖2所示，內(nèi)部有5個條形支撐板，與12棱滾筒板焊接在一起，通過鏈輪傳動帶動整個鋼輥轉(zhuǎn)動。

圖2 鋼輥結(jié)構(gòu)圖Fig.2 Structure schematic of steel roll

將一個應(yīng)變片粘貼在被測鋼輥的一條凸棱上，距離中間條形支撐板216.5 mm，應(yīng)變片沿軸向粘貼。無線電阻應(yīng)變儀應(yīng)變測量范圍為–15 000～15 000微應(yīng)變，測量精度±0.2%。試驗時，采用四分之一橋方式將應(yīng)變片與無線電阻應(yīng)變儀相連，無線電阻應(yīng)變儀通過磁座吸在鋼輥內(nèi)部，無線網(wǎng)關(guān)與電腦相連，試驗數(shù)據(jù)在BeeData軟件中顯示和保存并在Matlab8.0中進(jìn)行處理，圖3為數(shù)據(jù)采集系統(tǒng)布置圖，圖4為測試系統(tǒng)流程圖。

圖3 數(shù)據(jù)采集系統(tǒng)布置圖Fig.3 Arrangement diagram of data acquisition system

圖4 測試系統(tǒng)流程圖Fig.4 Test system flow chart

1.4.2 草捆所受的徑向力

為通過鋼輥應(yīng)變數(shù)據(jù)得到鋼輥所受的徑向作用力N，利用ANSYS軟件建立橋梁模型并進(jìn)行靜力分析。將被測凸棱及其兩側(cè)各一半滾筒板看作連續(xù)梁，5個條形支撐板看做滑動鉸支座建立模型。12棱滾筒板所用材料為Q235A，Q235A彈性模量為2.1×105 MPa，楊氏模量為0.3，密度為 7 850 kg/m3，連續(xù)梁尺寸 1 486 mm×74.5 mm×6 mm。施加應(yīng)力均布載荷（與實際卷壓過程中鋼輥所受載荷形式一致）求解，得到應(yīng)變結(jié)果，ANSYS模擬結(jié)果如圖5所示。徑向應(yīng)力σ的求解流程為：試驗得到變形數(shù)據(jù)→ANSYS模擬得到變形與試驗變形比較→不斷調(diào)整載荷大小至模擬變形與試驗變形的統(tǒng)計誤差在 0.2%以內(nèi)。其中，草捆受到的徑向應(yīng)力σ為徑向力N與受力面積A之比，計算公式如式（2）。

式中s（t）為任意時刻鋼輥受到的徑向應(yīng)力，kPa；N（t）為任意時刻鋼輥受到的徑向力，N；A為連續(xù)梁面積，mm2。

圖5 ANSYS模擬結(jié)果Fig.5 ANSYS simulation results

1.5 草捆所受徑向力分析

在草捆卷壓過程中，對草捆進(jìn)行受力分析如圖 6所示，圖中鋼輥編號為1～14，G為草捆重力，Ni(t)為1～14號鋼輥對草捆的徑向力，i=（1,2,3,…,14），F(xiàn)i(t)為1～14號鋼輥對草捆的摩擦力，i=（1,2,3,…,14），θi為正下方鋼輥與其他下半部分鋼棍之間的夾角i=（1,2,4,5,6）。草捆在鋼輥對其沿切線方向的摩擦力Fi(t)作用下做旋轉(zhuǎn)運動，在徑向力Ni(t)作用下壓縮變形。由于鋼輥分布位置不同，鋼輥對草捆的徑向力也不同。其中，上半部分7～14號鋼輥對草捆的徑向力為壓縮力，下半部分1～6號鋼輥對草捆既有壓縮力又承受著草捆重力，草捆重力分散到卷捆室下半部分鋼輥上。在卷壓過程中，可認(rèn)為鋼輥對草捆的壓縮力相等，若知道上半部分鋼輥和底部正下方 3號鋼輥對草捆的徑向力，可根據(jù)正下方鋼輥與其他下半部分鋼棍之間的夾角得到各個鋼輥對草捆施加的徑向力。

根據(jù)圖 6幾何關(guān)系可得出，卷捆過程中上半部分鋼輥對草捆施加的徑向力為

卷捆過程中下半部分鋼棍對草捆施加的徑向力

式中Nm為上半部分鋼棍對草捆的徑向力（m=7，8，9，10，11，12，13，14），N；Nn為下半部分鋼棍對草捆的徑向力（n=1，2，4，5，6），N；N1(t)～N14(t)為 1～14號鋼輥對草捆的徑向力，N；Ny(t)為鋼輥對草捆的壓縮力，N；Ng3(t)為正下方3號鋼輥受到的草捆重力，N。

當(dāng)草捆密度逐漸增大時，后倉門克服自身重力產(chǎn)生順時針旋轉(zhuǎn)的趨勢，捆室兩側(cè)液壓活塞拉桿開始對后倉門產(chǎn)生側(cè)拉力，使后倉門保持靜止?fàn)顟B(tài)。此時，可根據(jù)后倉門力矩平衡狀態(tài)方程以及各個鋼輥對草捆的徑向力，求得卷壓過程中使后倉門保持靜止的臨界側(cè)拉力值。

本文研究捆室上半部分鋼輥對草捆的徑向力。今后可通過測取捆室上部分鋼輥和正下方鋼輥的徑向力，再根據(jù)式（1）～（3）得到各個鋼輥對草捆的徑向力。

圖6 卷壓過程草捆受力分析Fig.6 Force analysis of bale in rotary compression

1.6 數(shù)據(jù)處理

用Matlab軟件對變形數(shù)據(jù)處理，剔除異常值，以每5 s 的數(shù)據(jù)平均值作為鋼輥變形的測試結(jié)果。最后，用ANSYS軟件將變形數(shù)據(jù)轉(zhuǎn)換為徑向應(yīng)力，得到草捆所受徑向應(yīng)力的變化過程。

2 結(jié)果與分析

2.1 試驗結(jié)果

根據(jù)測得的數(shù)據(jù)，繪制徑向應(yīng)力隨時間變化曲線，如圖7所示，前40 s為草物料填充階段，鋼輥與草捆之間主要為摩擦力作用，并不產(chǎn)生擠壓力，牧草松散；80 s后，進(jìn)入卷壓階段（AB段），草捆在鋼輥徑向力的作用下擠壓變形，草捆密度不斷增大。到160 s，喂料結(jié)束，草捆繼續(xù)在卷捆室內(nèi)運轉(zhuǎn)，進(jìn)入應(yīng)力松弛階段（BC段），草捆不再變形，應(yīng)力不斷減小。

圖7 徑向應(yīng)力隨時間變化曲線Fig.7 Change curve of radial stress with time

2.2 牧草保持變形條件下的應(yīng)力松弛

2.2.1 應(yīng)力松弛模型建立

農(nóng)業(yè)物料的應(yīng)力松弛現(xiàn)象是流變行為之一。不同物料種類、捆室結(jié)構(gòu)、加工方式的物料應(yīng)力松弛行為不同。探究卷壓過程中物料應(yīng)力松弛行為對圓捆機的優(yōu)化設(shè)計及草捆質(zhì)量的提高具有指導(dǎo)意義[21-25]。山東農(nóng)業(yè)大學(xué)李汝莘通過分析碎玉米秸稈卷壓過程中的變形特征，選用伯格斯模型模擬了應(yīng)力松弛過程，但伯格斯模型一般用來模擬蠕變過程，不能作為應(yīng)力松弛過程的模擬模型；東北農(nóng)業(yè)大學(xué)選用廣義Maxwell模型，通過Matlab對應(yīng)力松弛曲線進(jìn)行回歸分析，確定了方程式，利用數(shù)學(xué)軟件進(jìn)行分析可以方便地求出方程式中的任何參量，但是很難從計算結(jié)果對松弛過程中的流變學(xué)參量進(jìn)行分析[26-27]。

應(yīng)力松弛過程實質(zhì)是保持變形不變，彈性變形力恢復(fù)的過程，彈性變形為非彈性變形取代，彈性變形恢復(fù)力與其彈性變形恢復(fù)對應(yīng)。Maxwell模型是一個彈性元件和一個粘性元件的串聯(lián)，在保持其變形不變的條件下，彈簧的變形恢復(fù)力拉動串聯(lián)的阻尼器隨時間進(jìn)行松弛[28-30]。模擬一般的應(yīng)力松弛過程都以 Maxwell模型為基礎(chǔ)，且從試驗曲線，確定并聯(lián)基礎(chǔ)模型數(shù)及參量，對分析流變學(xué)過程具有直接指導(dǎo)作用。因此本文擬從試驗曲線出發(fā)，以Maxwell模型為基礎(chǔ)建立應(yīng)力松弛模型。

如圖7中BC段，當(dāng)t→∞時，應(yīng)力松弛不為0，存在著平衡應(yīng)力。應(yīng)力松弛過程的應(yīng)力曲線包括可松弛應(yīng)力和不可松弛應(yīng)力2部分，平衡應(yīng)力在模型中用1個并聯(lián)的彈簧表示，在應(yīng)力松弛中不能恢復(fù)。當(dāng)t→∞時，σBD(t)松弛到0。將BD段應(yīng)力松弛曲線在半對數(shù)坐標(biāo)中表示，如圖 8所示，是一條曲線，則該曲線由n個Maxwell模型并聯(lián)組成的n階Maxwell模型。其原理是將Maxwell模型用半對數(shù)坐標(biāo)（lnσ－t）表示，是一條斜率等于應(yīng)力松弛時間倒數(shù)的直線，即 lnσ（t）=lnσ+(–t/T)。因此選用 1 個n階 Maxwell模型和1個彈簧元件并聯(lián)組成的廣義Maxwell模型對應(yīng)力松弛過程進(jìn)行模擬，并進(jìn)一步通過殘數(shù)剩余法確定是由幾個Maxwell模型并聯(lián)組成的n階Maxwell模型以及模型中的參數(shù)。

圖8 第一Maxwell直線Fig.8 First Maxwell straight line

2.2.2 參數(shù)求解

n階Maxwell模型的表達(dá)式為

式中σBD（t）為任意時刻徑向應(yīng)力，kPa；i為 Maxwell模型的單元個數(shù)；Ai為第i個Maxwell模型的起始應(yīng)力，kPa；Ti為第i個Maxwell模型的應(yīng)力松弛時間，s；ih為第i個Maxwell模型中的粘性系數(shù)；Ei為第i個Maxwell模型中的彈性模量；ie為第i個Maxwell模型中的應(yīng)變，%。

在lnσBD尾端取幾個近似成直線關(guān)系的點，做曲線的擬合直線，并沿時間坐標(biāo)向應(yīng)力坐標(biāo)軸延長與其交至0.631（即為lnσ1），得到A1=e0.631=1.88；這條直線相當(dāng)于最大應(yīng)力為1.88的應(yīng)力松弛曲線；其斜率就是這條曲線的應(yīng)力松弛時間T1的倒數(shù)，T1=76.92 s。于是得出第一個Maxwell的應(yīng)力松弛方程式：

用原始應(yīng)力松弛曲線減去第一Maxwell直線得到應(yīng)力松弛剩余曲線，即σ2(t)=BD(t)－σ1(t)。做 lnσ2－t圖，得到第二條斜線，如圖 9所示，根據(jù)直線的斜率和截距得出第二個Maxwell的應(yīng)力松弛方程

圖9 第二Maxwell直線Fig.9 The second Maxwell straight line

至此，應(yīng)力松弛松弛已經(jīng)接近于0。原始應(yīng)力松弛曲線應(yīng)為以上 Maxwell模型和平衡應(yīng)力的疊加，應(yīng)力松弛方程式為

應(yīng)力松弛的模擬模型是由 2個麥克斯韋爾元件和 1個彈簧元件并聯(lián)組成的，如圖10所示。

圖10 應(yīng)力松弛模型Fig.10 Stress relaxation model

2.2.3 方程驗證

采用殘差平方和、相關(guān)指數(shù)對曲線模型的擬合效果進(jìn)行評價。如果殘差平方和占總平方和比例很小，即決定系數(shù)越接近于1，說明曲線擬合效果越好。殘差平方和和決定系數(shù)計算公式分別為

式中s為應(yīng)力實際測量值，kPa；為由方程式算得的應(yīng)力預(yù)測值，kPa；SS總為總平方和。

應(yīng)力測量值和應(yīng)力預(yù)測值如圖11所示，最終得到殘差平方和為7.01，決定系數(shù)為0.92。

圖11 應(yīng)力測量值和預(yù)測值Fig.11 Stress measurements and prediction values

對應(yīng)力松弛各參量進(jìn)行分析，應(yīng)力松弛模量，即 2個并聯(lián)的 Maxwell模型的模量分別為E1=1.88/ε；E2=6.1/ε，松弛過程中物料模量不變。應(yīng)力松弛延滯時間譜為T1=η1/E1=76.92s;T2=η2/E2=23.25s，反應(yīng)了物料松弛快慢程度且各元件應(yīng)力松弛時間相互獨立。應(yīng)力松弛起始應(yīng)力為σ=σ1+σ2+3.65=11.63 kPa，平衡應(yīng)力為 3.65 kPa，說明草捆內(nèi)部還存在著一定的平衡應(yīng)力。

3 結(jié) 論

通過牧草卷壓試驗，在對應(yīng)力松弛試驗曲線和模型原理分析的基礎(chǔ)上，確定了牧草應(yīng)力松弛行為可用兩個Maxwell元件和一個彈簧元件并聯(lián)組成的的廣義麥克斯韋爾模型模擬。獲取了牧草應(yīng)力松弛過程中的流變參數(shù)，各元件應(yīng)力松弛時間分別為76.92 s、23.25 s，最小平衡應(yīng)力為3.65 kPa。該研究對進(jìn)一步分析喂入量、含水率、鋼輥轉(zhuǎn)速、壓縮力、草捆密度等參數(shù)之間的關(guān)系，從而優(yōu)化打捆過程及壓縮設(shè)備，獲取較低功耗下較高質(zhì)量的草捆提具有指導(dǎo)意義。

[1]王春光. 鋼輥外卷式圓捆機結(jié)構(gòu)與原理[J]. 農(nóng)業(yè)機械，2008(34)：43－45.Wang Chunguang. Structure and principle of steel roller winding type round baler[J]. Farm Machinery, 2010(34):43－45. (in Chinese with English abstract)

[2]王德福，張全超，楊星，等. 秸稈圓捆機捆繩機構(gòu)的參數(shù)優(yōu)化與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(14)：55－61.Wang Defu, Zhang Quanchao, Yang Xing, et al. Parameter optimization and experiment of rope-winding mechanism of straw round balers[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(14): 43－45. (in Chinese with English abstract)

[3]李葉龍，王德福，王沫，等. 稻稈圓捆機輥盤式卷捆機構(gòu)的設(shè)計及參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(6)：27－34.Li Yelong, Wang Defu, Wang Mo, et al. Design and parameters optimization of roll-disk baling mechanism for rice straw round baler[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2017, 33(6): 43－45. (in Chinese with English abstract)

[4]楊明韶，張永，李旭英.粗纖維物料壓縮過程的一般流變規(guī)律的探討[J]. 農(nóng)業(yè)工程學(xué)報， 2002，18(1)：135－137.Yang Mingshao, Zhang Yong, Li Xuying. Rheological law of the crop stem fibrous material during compression process[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2002, 18(1): 135－137. (in Chinese with English abstract)

[5]馬彥華，宣傳忠，武佩，等.玉米秸稈振動壓縮過程的應(yīng)力松弛試驗[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(19)：88－94.Ma Yanhua, Xuan Chuanzhong, Wu Pei, et al. Experiment on stress relaxation of corn stover during compression with assisted vibration[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(19): 88－94. (in Chinese with English abstract)

[6]李旭英，楊明韶，魯國成，等.苜蓿壓捆過程中壓縮與恢復(fù)應(yīng)力傳遞規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(16)：61－67.Li Xunying, Yang Mingshao, Lu Guocheng, et al. Transfer rule of compression and springback stress in compression process of alfalfa[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014,30(16): 61－67. (in Chinese with English abstract)

[7]楊明韶，杜建民，李旭英. 新鮮草物料壓縮過程中應(yīng)力松弛的試驗研究[J]. 農(nóng)機化研究， 2006，2：100－103.Yang Mingshao, Du Jianmin, Li Xunying. An experimental study of fresh herbage material stress relaxation in the process of compression[J]. Journal of Agricultural Mechanization Research, 2006, 2: 100－103. (in Chinese with English abstract)

[8]宮澤奇. 青貯玉米秸稈螺旋致密成型工藝及壓縮過程流變學(xué)試驗研究[D]. 北京：中國農(nóng)業(yè)大學(xué)，2017.Gong Zeqi. The Experimental Study on Rheology of Silagecorn Stalk in the Mechanized Spiral Dense Forming Process[D]. Beijing: Chian Agricultural University, 2017. (in Chinese with English abstract)

[9]朱凱，牛智有.粉體棉花秸稈的應(yīng)力松弛特性與模型建立[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報，2014，33(5)：130－134.Zhu Kai, Niu Zhiyou. Compression stress relaxation characteristic of cotton stalk particles and its model building of cot[J]. Journal of Huazhong Agricultural University, 2014,33(5): 130－134. (in Chinese with English abstract)

[10]廖娜，陳龍健，黃光群，等. 玉米秸稈木質(zhì)纖維含量與應(yīng)力松弛特性關(guān)聯(lián)度研究[J]. 農(nóng)業(yè)機械學(xué)報，2011，42(12)：127－132.Liao Na, Chen Longjian, Huang Guangqun, et al. Grey relation analysis of lignocellulose content and compression stress relaxation of corn stalk[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(12): 127－132.(in Chinese with English abstract)

[11]郭文斌，王春光，劉百順. 馬鈴薯應(yīng)力松弛特性[J]. 農(nóng)業(yè)機械學(xué)報，2008(2)：205－207.Guo Wenbin, Wang Chunguang, Liu Baishun. Stress relaxation characteristics of potato[J]. Transactions of the Chinese Society for Agricultural Machinery, 2008(02): 205－207. (in Chinese with English abstract)

[12]許永亮，熊善柏，趙思明. 蒸煮工藝和化學(xué)成分對米飯應(yīng)力松弛特性的影響[J]. 農(nóng)業(yè)工程學(xué)報，2007，23(10)：235－240.Xu Yongliang, Xiong Shanbo, Zhao Siming. Effect of cooking technology and chemical components on stressrelaxation property of cooked rice[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(10): 235－240. (in Chinese with English abstract)

[13]李汝莘，耿愛軍，趙何，等. 碎玉米秸稈卷壓過程的流變行為試驗[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(18)：30－35.Li Ruxin，Geng Aijun，Zhao He，et al. Rheologic behavior of chopped corn stalks during rotary compression[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(18): 30－35. (in Chinese with English abstract)

[14]雷軍樂，王德福，張全超，等. 完整稻稈卷壓過程應(yīng)力松弛試驗[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(8)：76－83.Lei Junle, Wang Defu, Zhan Quanchao, et al. Experiment on stress relaxation characteristics of intact rice straw during rotary compression[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015,31(8): 30－35. (in Chinese with English abstract)

[15]李旭英. 草物料壓縮流變動力學(xué)特性研究及壓縮過程優(yōu)化[D]. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2006.Li Xuying. Rheodynamic Characteristics of Herbage Materials in Compression and Optimization of Compression Process[D]. Hohhot: Inner Mongolia Agricultural University,2006. (in Chinese with English abstract)

[16]Hu J, Lei T. Stress relaxation in straw molding[J].BioResources, 2009, 4(3): 1158－1167.

[17]周慶安，楊培志. 我國苜蓿飼草產(chǎn)業(yè)的發(fā)展現(xiàn)狀與思考[J].中國奶牛，2017(5)：54－56.Zhou Qingan, Yang Peizhi. Development situation and think of alfalfa forage grass industry in china[J]. China Dairy Cattle,2017(5): 54－56. (in Chinese with English abstract)

[18]張玲，李菲菲，孫艷梅，等. 國內(nèi)主要品種苜蓿營養(yǎng)價值分析與育種現(xiàn)狀及發(fā)展歷程[J]. 中國奶牛，2017(11)：43－49.Zhang Ling, Sun Feifei, Sun Yanmei, et al. Evaluation on the nutritional value, breeding and development history of main alfalfa (Medicago Sativa)varieties in China[J]. China Dairy Cattle, 2017(11): 43－49. (in Chinese with English abstract)

[19]中華人民共和國機械行業(yè)標(biāo)準(zhǔn).牧草撿拾器：JB/T5160-2010[S]. 北京：機械電子工業(yè)部，1991.

[20]趙何. 碎玉米秸稈卷壓特性的研究[D]. 泰安：山東農(nóng)業(yè)大學(xué)，2010.Zhao He. Study on Rotary Compression Characteristic of Chopped Corn Stalks[D]. Taian: Shandong Agricultural University, 2010. (in Chinese with English abstract)

[21]聶毓琴，馬洪順，韓志武. 薇菜壓縮應(yīng)力松弛與蠕變力學(xué)特性研究[J]. 農(nóng)業(yè)機械學(xué)報，2005(4)：89－91.Nie Yuqin, Ma Hongshun, Han Zhiwu. Research on stress relaxation and creep mechanics behaviors of osmunda cinnamomea[J]. Transactions of the Chinese Society for Agricultural Machinery, 2005(4): 89－91. (in Chinese with English abstract)

[22]王菊霞，崔清亮，李紅波，等. 基于流變特性的不同品種蘋果果皮質(zhì)地評價[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(21)：305－314.Wang Junxia, Cui Qingliang, Li Hongbo, et al. Evaluation on peels texture of different apple cultivars based on rheological properties[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(21): 305－314. (in Chinese with English abstract)

[23]馬彥華，宣傳忠，武佩，等. 玉米秸稈振動壓縮過程的應(yīng)力松弛試驗[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(19)：88－94.Ma Yanhua, Xuan Chuanzhong, Wu Pei, et al. Experiment on stress relaxation of corn stover during compression with assisted vibration[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(19): 88－94. (in Chinese with English abstract)

[24]楊曉清，王春光. 果品靜載流變特性的研究進(jìn)展[J]. 農(nóng)業(yè)工程學(xué)報，2005，21(9)：178－182.Yang Xiaoqing, Wang Chunguang. Review of research advances in rheological properties of fruits under static loading[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005,21(9): 178－182. (in Chinese with English abstract)

[25]Aline Bauer, Luo Gu, Brian Kwee, et al. Hydrogel substrate stress-relaxation regulates the spreading and proliferation of mouse myoblasts[J]. Acta Biomaterialia, 2017.

[26]楊明韶. 農(nóng)業(yè)物料流變學(xué)[M]. 北京：中國農(nóng)業(yè)出版社. 2010.

[27]楊明韶，韓寶生. 草類物料壓縮過程中流變學(xué)研究的最新進(jìn)展[J]. 農(nóng)機化研究，2005(3)：15－20.Yang Mingshao, Han Baosheng. The new progress of the rheological study in the herbage material compression process[J]. Journal of Agricultural Mechanization Research,2005, 5(3): 15－20. (in Chinese with English abstract)

[28]楊明韶，李旭英，楊紅蕾. 牧草壓縮過程的研究[J]. 農(nóng)業(yè)工程學(xué)報，1996，14(1)：60－64.Yang Mingshao, Li Xuying, Yang Honglei. Research on the hay compression process[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),1996, 14(1): 60－64. (in Chinese with English abstract)

[29]杜建明. 新鮮草物料壓縮過程的流變學(xué)研究[D]. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2005.Du Jianmin. The Experimental Study on Rheology of Fresh Herbage Materials in the Process of the compression[D].Hohhot: Inner Mongolia Agricultural University, 2005. (in Chinese with English abstract)

[30]王春光，楊明韶，高煥文，等.牧草在高密度壓捆時的應(yīng)力松弛研究[J]. 農(nóng)業(yè)工程學(xué)報，1997，13(3)：48－52.Wang Chunguang, Yang Mingshao, Gao Huanwen, et al.Study on stress-relaxation of hay during baling under high density[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 1997, 13(3): 48－52. (in Chinese with English abstract)