田間機(jī)收元胡離散元模型參數(shù)標(biāo)定與試驗(yàn)

摘要：

為提高機(jī)收元胡所用離散元仿真模型參數(shù)準(zhǔn)確性，通過物理試驗(yàn)與仿真試驗(yàn)結(jié)合的方法對(duì)元胡離散元參數(shù)進(jìn)行標(biāo)定。首先，采用物理試驗(yàn)確定元胡基本物性參數(shù)范圍；其次，依據(jù)物理試驗(yàn)測(cè)定的參數(shù)范圍進(jìn)行休止角的顯著性試驗(yàn)，得到對(duì)元胡休止角影響由大到小的順序?yàn)樵g滾動(dòng)摩擦系數(shù)、元胡與鋼板間滾動(dòng)摩擦系數(shù)、元胡間碰撞恢復(fù)系數(shù)；再通過最陡爬坡試驗(yàn)縮小顯著性參數(shù)取值范圍，并通過響應(yīng)面試驗(yàn)建立休止角與顯著因素的最優(yōu)二階回歸方程；以物理試驗(yàn)休止角為目標(biāo)獲取最佳參數(shù)組合：元胡間滾動(dòng)摩擦系數(shù)0.08、元胡與鋼板間滾動(dòng)摩擦系數(shù)0.32、元胡間碰撞恢復(fù)系數(shù)0.49。最后，在對(duì)比標(biāo)定模型參數(shù)下仿真堆積體與實(shí)際堆積體，測(cè)得二者的休止角和堆體外形輪廓無(wú)明顯差異，表明該最佳參數(shù)組合可作為元胡機(jī)收過程中離散元仿真模型參數(shù)。

關(guān)鍵詞：元胡；離散元模型；參數(shù)標(biāo)定；休止角；收獲機(jī)械

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

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

文章編號(hào)：2095-5553 （2024） 04-0032-07

收稿日期：2022年12月6日" 修回日期：2023年5月16日

基金項(xiàng)目：陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃項(xiàng)目（2023—JC—YB—018）

第一作者：劉珂，男，1994年生，西安人，碩士研究生；研究方向?yàn)楝F(xiàn)代設(shè)計(jì)理論及方法。E-mail： liukeqwq@163.com

通訊作者：魏偉鋒，男，1978年生，西安人，博士，講師；研究方向?yàn)榫軝C(jī)械CAD/CAM。E-mail： sxhzwwf@163.com

Parameter calibration and test of rhizoma corydalis discrete element model of

field machine harvester

Liu Ke1， Wei Weifeng1， 2， Zhu Jiyue1， 2， Li Gaowei1， Gao Zhennan1

（1. School of Mechanical Engineering， Shaanxi University of Technology， Hanzhong， 723000， China；

2. Shaanxi Key Laboratory of Industrial Automation， Hanzhong， 723000， China）

Abstract：

In order to improve the accuracy of the discrete element simulation parameters used in the machine-received element， the discrete element parameters of the Rhizoma corydalis were calibrated by combining the physical test with the simulation test. Firstly， the range of basic physical property parameters of Rhizoma corydalis was determined by physical tests. Secondly， the significance test of the angle of repose was carried out according to the parameter range determined by the physical test. The order of the influence on the angle of repose of the Rhizoma corydalis was as follows： the coefficient of rolling friction between the Rhizoma corydalis， the coefficient of rolling friction between the Rhizoma corydalis and the steel plate， the coefficient of collision recovery between the Rhizoma corydalis. And according to the parameter range of physical test， the significance test of angle of repose was carried out， and the significant parameters that had an effect on the characteristic parameters were determined as follows： the rolling friction coefficient between Rhizoma corydalis and Rhizoma corydalis， the rolling friction coefficient between the Rhizoma corydalis and the steel plate， and the collision recovery coefficient between Rhizoma corydalis and Rhizoma corydalis. Again through the steepest uphill test narrowing the scope of significant parameter selection， and the angle of repose and significant factors based on response surface tests the optimal quadratic regression equation， then aiming at physical test rest angle to obtain the best parameter combination： Rhizoma corydalis rolling friction coefficient between 0.08， Rhizoma corydalis， and steel rolling friction coefficient between 0.32， Rhizoma corydalis collision between recovery coefficient of 0.49. At last， the simulation test was carried out under the calibrated model parameters and compared with the stack test， and no significant difference was found between the angle of repose and the contour of the pile， which showed that the optimal parameter combination" could be used as the discrete element simulation model parameters in the process of Rhizoma corydalis machine harvesting.

Keywords：

rhizoma corydalis; DEM; parameter calibration; angle of repose; harvesting machinery

0 引言

元胡主要種植地為陜西省漢中市，是一味大宗中藥材［1］，目前其采收方式主要為人工采收。人工用釘耙刨挖，再將土壤和元胡混合物利用篩網(wǎng)篩分。該挖掘過程耗時(shí)長(zhǎng)、效率低［2］。因此，元胡的機(jī)械化收獲勢(shì)在必行。元胡收獲機(jī)提升機(jī)構(gòu)是類似于鏤空的傳送帶，若傳送速度過高則易使元胡被拋出不利于收集，若傳送速度過低易造成堵塞，提升機(jī)構(gòu)作為元胡收獲機(jī)的關(guān)鍵技術(shù)之一，其提升效率直接影響整機(jī)工作效率。

在收獲過程中，元胡與提升機(jī)構(gòu)間的運(yùn)動(dòng)機(jī)理十分復(fù)雜。利用離散元法對(duì)元胡提升機(jī)構(gòu)系統(tǒng)的仿真研究，有利于分析元胡在提升機(jī)構(gòu)中的運(yùn)動(dòng)機(jī)理，進(jìn)而實(shí)現(xiàn)對(duì)元胡收獲機(jī)提升機(jī)構(gòu)機(jī)械參數(shù)優(yōu)化。本文以元胡為研究對(duì)象，對(duì)其離散元仿真參數(shù)進(jìn)行標(biāo)定，以推進(jìn)離散元法在元胡收獲裝備研發(fā)中的應(yīng)用。

離散元仿真時(shí)，需定義仿真模型的物性參數(shù)，主要包括密度、剪切模量、泊松比、顆粒之間的一系列接觸參數(shù)等［3， 4］。其中密度、剪切模量、泊松比為可直接測(cè)量的顆粒基礎(chǔ)物理參數(shù)。由于仿真顆粒與真實(shí)顆粒之間幾何形態(tài)存在差異，使得直接測(cè)量的接觸參數(shù)不能直接用于離散元仿真中，需要對(duì)其重新標(biāo)定。目前接觸參數(shù)主要是通過直接測(cè)量和虛擬標(biāo)定兩種方法［5， 6］。然而，如顆粒間動(dòng)摩擦、顆粒間碰撞恢復(fù)系數(shù)等難以實(shí)際測(cè)量，仿真顆粒與真實(shí)顆粒幾何形態(tài)存在差異［7， 8］，因此，直接測(cè)量與虛擬標(biāo)定兩種手段結(jié)合已經(jīng)成為主流。

目前，國(guó)內(nèi)外農(nóng)作物顆粒離散元參數(shù)標(biāo)定主要包括玉米、水稻、豌豆籽粒、小麥等大宗植物。Wang等［9］通過試驗(yàn)方法標(biāo)定玉米離散元模型的摩擦因數(shù)和剛度，并通過筒倉(cāng)卸料和斗裝試驗(yàn)驗(yàn)證；張榮芳等［10］通過休止角標(biāo)定水稻離散元模型參數(shù)，并通過排種器排種試驗(yàn)驗(yàn)證參數(shù)；張克平等［11］通過試驗(yàn)的方法確定豌豆籽粒參數(shù)范圍，再通過堆積試驗(yàn)標(biāo)定豌豆籽粒離散元模型參數(shù)；劉凡一等［12］通過堆積試驗(yàn)標(biāo)定小麥離散元模型參數(shù)，以實(shí)際休止角為目標(biāo)值，通過離散元仿真，建立回歸模型，求解最佳參數(shù)組合；而對(duì)元胡顆粒離散元仿真參數(shù)標(biāo)定的研究較少。元胡與玉米、豌豆籽粒相比在質(zhì)量、幾何形狀上差異較大，如果采用直接測(cè)量的方法標(biāo)定其離散元仿真參數(shù)，會(huì)引起較大誤差。

本文采用物理試驗(yàn)與仿真試驗(yàn)相結(jié)合的方法，以物理試驗(yàn)測(cè)量的參數(shù)范圍作為離散元仿真試驗(yàn)參數(shù)范圍選取依據(jù)，以實(shí)際自然休止角為目標(biāo)，依次采用顯著性試驗(yàn)、響應(yīng)面試驗(yàn)對(duì)元胡的離散元參數(shù)進(jìn)行標(biāo)定和校準(zhǔn)，獲取元胡離散元參數(shù)，為元胡收獲過程中，提升機(jī)構(gòu)、篩分機(jī)構(gòu)等作業(yè)仿真提供參數(shù)依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料與物理參數(shù)測(cè)定

本研究以秦巴山區(qū)人工種植，自然氣候條件下生長(zhǎng)的元胡為試驗(yàn)對(duì)象，收獲后的元胡為散體顆粒，其基本物理參數(shù)包括元胡外形尺寸（長(zhǎng)度、寬度、厚度）、密度、含水率、泊松比、彈性模量和剪切模量。

1.1.1 基本物理參數(shù)測(cè)定

元胡物理參數(shù)中，易獲取的參數(shù)包括外形尺寸、密度、含水率，因此首先對(duì)上述參數(shù)進(jìn)行測(cè)量。對(duì)元胡表面土壤清理后，隨機(jī)選4顆元胡，采用精度為0.1 g的電子天平測(cè)量質(zhì)量，采用排水法測(cè)量其密度。采用精度為0.02 mm的游標(biāo)卡尺測(cè)量元胡的外形尺寸，將元胡高不足10 mm稱為扁平形，高大于10 mm稱為類球形；采用型號(hào)為DSH-50A-1的水分測(cè)定儀測(cè)量元胡含水率。

上述試驗(yàn)均重復(fù)5次取平均值，元胡物理參數(shù)平均值如表1所示。

1.1.2 泊松比和彈性模量

元胡外形尺寸差異較大，其泊松比難以采用常規(guī)方法測(cè)量。本研究將元胡切成長(zhǎng)方體，利用Brookfield CT3質(zhì)構(gòu)分析儀對(duì)元胡進(jìn)行壓縮試驗(yàn)，測(cè)量元胡彈性模量，如圖1所示。

Brookfield CT3質(zhì)構(gòu)分析儀所用壓力測(cè)量范圍為0～25 kN，元胡壓縮試驗(yàn)時(shí)，將長(zhǎng)方體元胡放置在平板上，將長(zhǎng)方體側(cè)面面積較大一側(cè)放置在平板上，以30 mm/min進(jìn)行加載，探頭半徑為5 mm，垂直向下加載3 s后停機(jī)。上述試驗(yàn)重復(fù)5次取平均值，將軟件后處理模塊得到元的胡壓縮試驗(yàn)過程中的位移—載荷數(shù)據(jù)，通過式（1）計(jì)算可得元胡彈性模量為1.06×107 Pa。

E=σε

（1）

式中：

E——元胡彈性模量，Pa；

σ——最大壓應(yīng)力，Pa；

ε——線應(yīng)變。

1.探針 2.食品物性分析儀 3.試驗(yàn)平面

泊松比是指顆粒在單向受拉或受壓時(shí)，橫向正應(yīng)變與軸向正應(yīng)變的比值，可以反映顆粒變形時(shí)的彈性系數(shù)，其計(jì)算如式（2）所示。

τ=εcεt=ΔL/LΔH/H

（2）

式中：

τ——泊松比；

εc——長(zhǎng)方塊元胡橫向應(yīng)變；

εt——長(zhǎng)方塊元胡厚度方向應(yīng)變；

ΔL——長(zhǎng)方塊元胡橫向變形量，mm；

L——長(zhǎng)方塊元胡橫向原長(zhǎng)，mm；

ΔH——

長(zhǎng)方塊元胡厚度方向變形，mm；

H——長(zhǎng)方塊元胡厚度方向原長(zhǎng)，mm。

泊松比測(cè)定試驗(yàn)采用Brookfield CT3質(zhì)構(gòu)分析儀進(jìn)行加載，試驗(yàn)前將元胡處理為長(zhǎng)方塊。元胡壓縮試驗(yàn)時(shí)，以30 mm/min進(jìn)行加載，加載3 s，探頭半徑為5 mm，加載方向垂直向下，加載3 s后停機(jī)。測(cè)量元胡外形尺寸，重復(fù)測(cè)量5次，取平均值，通過式（2）得元胡泊松比為0.138。由式（3）可得元胡剪切模量為4.66×106 Pa。

η=E2（1+τ）

（3）

式中：

η——元胡剪切模量，Pa。

1.2 物理試驗(yàn)

離散元仿真過程中元胡顆粒之間，元胡顆粒與設(shè)備之間接觸參數(shù)有：元胡與元胡間碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)、滾動(dòng)摩擦系數(shù)、元胡與鋼材間碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)、滾動(dòng)摩擦系數(shù)［13］。為提高離散元仿真接觸參數(shù)的準(zhǔn)確性，首先通過物理試驗(yàn)測(cè)定上述參數(shù)范圍，再以物理試驗(yàn)結(jié)果為參考進(jìn)行離散元仿真試驗(yàn)。

1.2.1 靜摩擦系數(shù)測(cè)定

斜面滑動(dòng)法是測(cè)定靜摩擦系數(shù)的常用方法，物體靜止在斜板上如圖2所示，隨著斜面角度增大，當(dāng)元胡出現(xiàn)滑動(dòng)趨勢(shì)時(shí)，臨界角與元胡和鋼板間靜摩擦系數(shù)之間的關(guān)系如式（4）、式（5）所示。

Gsinθ=γGcosθ

（4）

γ=tanθ

（5）

式中：

G——元胡重力，N；

θ——臨界角，（°）；

γ——元胡與鋼板間靜摩擦系數(shù)。

物理試驗(yàn)中，采用不銹鋼作為試驗(yàn)斜面，為防止元胡在斜面上滾動(dòng)，將12顆元胡相互黏結(jié)放置在平板上，試驗(yàn)過程中不銹鋼板自由端繞固定端緩慢旋轉(zhuǎn)，當(dāng)元胡在不銹鋼板上開始滑動(dòng)時(shí)，利用角度尺測(cè)量鋼板與水平面的夾角θ，重復(fù)試驗(yàn)20次，夾角平均值為26.2°，由式（5）可得元胡與鋼板間靜摩擦系數(shù)為0.488。

如圖3所示，將上述12顆相互黏結(jié)的元胡顆粒固定在平板，隨機(jī)選取一顆類球形元胡與扁平形元胡，分別將其隨機(jī)放置在黏結(jié)在一起的元胡顆粒上，將平板緩慢抬起，當(dāng)元胡開始滑動(dòng)時(shí)，讀取平板與水平面的夾角，重復(fù)試驗(yàn)20次，放置類球形元胡平板與水平面夾角平均值為29.76°，放置扁平形元胡平板與水平面夾角平均值為37.12°。元胡間元胡靜摩擦系數(shù)為0.572～0.757。

1.2.2 滾動(dòng)摩擦系數(shù)測(cè)定

滾動(dòng)摩擦是物質(zhì)基本屬性之一，其大小與物質(zhì)本身性質(zhì)有關(guān)，如：質(zhì)量、表面形狀、材料等。試驗(yàn)過程中，假設(shè)元胡顆粒運(yùn)動(dòng)為純滾動(dòng)，且不考慮運(yùn)動(dòng)過程中所受靜摩擦力，根據(jù)能量守恒定律

mgSsinβ=μrmg（Scosβ+l）

（6）

式中：

m——元胡質(zhì)量，g；

g——重力加速度，m/s2；

S——

元胡在斜置鋼板上滾動(dòng)距離，mm；

β——斜置鋼板與水平面夾角，（°）；

μr——滾動(dòng)摩擦系數(shù)；

l——元胡在水平面上滾動(dòng)距離，mm。

通過預(yù)試驗(yàn)，確定斜置鋼板與水平面夾角β=25°，斜置鋼板上滾動(dòng)距離S為50 mm。預(yù)試驗(yàn)發(fā)現(xiàn)類球形元和和扁平形元胡滾動(dòng)距離相差較大，因此，本文對(duì)類球形元胡與扁平形元胡分別進(jìn)行滾動(dòng)試驗(yàn)，如圖4所示。根據(jù)上述試驗(yàn)要求進(jìn)行20次重復(fù)滾動(dòng)試驗(yàn)，最終類球形元胡滾動(dòng)試驗(yàn)平均值L1=454.2 mm，扁平形元胡滾動(dòng)試驗(yàn)平均值L2=143.45 mm；將上述測(cè)定值分別代入式（6）中，得到物理試驗(yàn)下元胡與鋼板之間的滾動(dòng)摩擦系數(shù)范圍為0.110 8～0.359 7。

1.3 物理休止角測(cè)定

采用提升無(wú)底圓筒方法進(jìn)行物理休止角試驗(yàn)［14］，測(cè)量裝置由無(wú)底圓筒（內(nèi)徑為100 mm、高度為300 mm）和底板（長(zhǎng)為600 mm、寬為600 mm）組成，材質(zhì)均為不銹鋼，以5 mm/s速度向上提升圓筒，待元胡顆粒穩(wěn)定后，測(cè)量其底部圓半徑和堆積體高度，求得堆積體高度與底部半徑比值，再利用反正切函數(shù)求休止角，取10次試驗(yàn)結(jié)果，得休止角平均值為36.82°。

2 元胡離散元模型建立及表征參數(shù)標(biāo)定

2.1 元胡離散元模型建立

元胡分類及離散元模型如圖5所示。

為提高仿真試驗(yàn)準(zhǔn)確性，建立元胡離散元仿真模型時(shí)，應(yīng)盡可能接近元胡物理形態(tài)［15］。由于元胡顆粒間差異性大，外形輪廓復(fù)雜，按照元胡高度是否高于10 mm分為類球形和扁平形兩類，類球形元胡與扁平形元胡數(shù)量比為3∶1。本研究采用多球形顆粒填充方法建立元胡仿真離散元模型［16］，以物理測(cè)量數(shù)據(jù)為基礎(chǔ)，采用多顆粒填充方法建立元胡離散元仿真模型，離散元仿真中類球形元胡與扁平狀元胡生成比例與實(shí)際相同（圖5）。

2.2 離散元表征參數(shù)標(biāo)定

2.2.1 Plackett-Burman試驗(yàn)

采用Design-Expert軟件進(jìn)行Plackett-Burman試驗(yàn)設(shè)計(jì)，離散元仿真中鋼材密度為7 850 kg/m3，泊松比為0.3，剪切模量7.9×104 MPa［14］；元胡間靜摩擦系數(shù)采用物理試驗(yàn)測(cè)量得平均值0.66，元胡與鋼板碰撞恢復(fù)系數(shù)參考根莖類作物（三七、微型馬鈴薯）取值范圍平均值0.6［17， 18］。其余試驗(yàn)參數(shù)范圍依據(jù)物理試驗(yàn)結(jié)果，如表2所示。表2以元胡堆積后的休止角為目標(biāo)值，通過Plackett-Burman試驗(yàn)篩選出對(duì)目標(biāo)值影響顯著性的參數(shù)。Plackett-Burman試驗(yàn)方案及結(jié)果如表3所示。

利用Design-Expert軟件對(duì)試驗(yàn)結(jié)果進(jìn)行顯著性分析［19］，得到各離散元模型參數(shù)顯著性結(jié)果，如表4所示。由表4可知，元胡間滾動(dòng)摩擦系數(shù)的Plt;0.01，對(duì)仿真試驗(yàn)休止角的影響極顯著；元胡與鋼板間滾動(dòng)摩擦系數(shù)、元胡間碰撞恢復(fù)系數(shù)的Plt;0.05，對(duì)仿真試驗(yàn)休止角影響顯著；其他仿真試驗(yàn)參數(shù)的Pgt;0.05，對(duì)仿真試驗(yàn)休止角無(wú)明顯影響。

2.2.2 最陡爬坡試驗(yàn)設(shè)計(jì)

以Plackett-Burman試驗(yàn)中3個(gè)有顯著影響的離散元參數(shù)進(jìn)行最陡爬坡試驗(yàn)，其目的是縮小離散元仿真時(shí)顯著性參數(shù)取值范圍，獲取最優(yōu)離散元模型參數(shù)組合。最陡爬坡試驗(yàn)結(jié)果如表5所示，第4組仿真試驗(yàn)與自然休止角相對(duì)誤差最小，因此在后續(xù)響應(yīng)面試驗(yàn)中，分別以最陡爬坡中第3組和第5組作為響應(yīng)面試驗(yàn)參數(shù)取值范圍。

2.2.3 Box-Behnken試驗(yàn)

在Design-Expert軟件中，以最陡爬坡仿真試驗(yàn)序號(hào)4作為中心點(diǎn)，序號(hào)3、序號(hào)5分別作為低、高水平進(jìn)行響應(yīng)面設(shè)計(jì)，響應(yīng)面試驗(yàn)設(shè)計(jì)方案及結(jié)果如表6所示，仿真試驗(yàn)中，其余非顯著性參數(shù)取物理測(cè)量值平均值。

響應(yīng)面試驗(yàn)方差分析結(jié)果如表7所示，其中元胡間滾動(dòng)摩擦系數(shù)、元胡-鋼板間滾動(dòng)摩擦系數(shù)對(duì)休止角影響極其顯著，元胡間碰撞恢復(fù)系數(shù)的二次項(xiàng)對(duì)元胡休止角影響顯著，其余項(xiàng)均對(duì)元胡休止角無(wú)明顯影響。該優(yōu)化后二階回歸模型Plt;0.01，失擬項(xiàng)P=0.109 3gt;0.05，表明該二階回歸模型擬合較好，無(wú)失擬現(xiàn)象發(fā)生。綜上，該回歸模型顯著，能夠準(zhǔn)確可靠的反映真實(shí)情況，可用于離散元仿真對(duì)元胡自然休止角預(yù)測(cè)分析。

通過Design-Expert軟件對(duì)響應(yīng)面試驗(yàn)結(jié)果建立回歸模型，并對(duì)其模型進(jìn)行優(yōu)化，得到仿真試驗(yàn)休止角與影響因素的二階回歸方程

θ=

37.53+0.99A+0.70B+0.18C-0.29AC-0.28A2-0.58C2

（7）

2.2.4 對(duì)比試驗(yàn)

在Design-Expert 8.05軟件中，以物理試驗(yàn)休止角（36.82°）為目標(biāo)，對(duì)優(yōu)化后的二階回歸模型求解，解得最佳仿真參數(shù)組合：元胡間滾動(dòng)摩擦系數(shù)0.08、元胡與鋼板間滾動(dòng)摩擦系數(shù)0.32、元胡間碰撞恢復(fù)系數(shù)0.49。

為驗(yàn)證標(biāo)定后元胡離散元模型參數(shù)的可靠性和準(zhǔn)確性，以上述所確定接觸參數(shù)作為元胡離散元仿真模型接觸參數(shù)，進(jìn)行3次仿真試驗(yàn)，得到元胡休止角分別為36.65°、35.84°、37.25°，平均值為36.58°，與物理試驗(yàn)休止角（36.82°）相對(duì)誤差僅為0.65%，表明所標(biāo)定離散元參數(shù)可用于元胡離散元仿真中，休止角對(duì)比如圖6所示。

3 結(jié)論

1）" 通過物理試驗(yàn)測(cè)定得到元胡的基本物性參數(shù)；元胡間靜摩擦系數(shù)平均值為0.66；元胡與鋼板間靜摩擦系數(shù)0.488、滾動(dòng)摩擦系數(shù)平均值為0.23。

2）" 通過Plackett-Burman試驗(yàn)，得出對(duì)離散元仿真中休止角影響顯著參數(shù)為：元胡間滾動(dòng)摩擦系數(shù)和元胡間碰撞恢復(fù)系數(shù)，元胡與鋼板間滾動(dòng)摩擦系數(shù)。

3）" 通過響應(yīng)面試驗(yàn)，建立休止角與顯著參數(shù)的最優(yōu)二階回歸方程，并以物理試驗(yàn)休止角36.82°為目標(biāo)值進(jìn)行求解，得到最佳離散元模型參數(shù)組合：元胡間滾動(dòng)摩擦系數(shù)為0.08、和元胡間碰撞恢復(fù)系數(shù)0.49、元胡與鋼板間滾動(dòng)摩擦系數(shù)為0.32。

4）" 以試驗(yàn)設(shè)計(jì)所得最佳元胡離散元參數(shù)進(jìn)行休止角仿真，測(cè)得仿真結(jié)果平均值為36.58°，物理試驗(yàn)結(jié)果平均值為36.82°，兩者相對(duì)誤差為0.65%，仿真值與物理試驗(yàn)值相對(duì)誤差在5%以內(nèi)，表明本文所標(biāo)定的元胡離散元參數(shù)可用于后續(xù)元胡收獲機(jī)提升機(jī)構(gòu)離散元仿真。

參 考 文 獻(xiàn)

［1］ 呂美巧. 小型元胡收獲機(jī)的研制［D］. 南京： 南京農(nóng)業(yè)大學(xué)， 2013.

Lü Meiqiao. Development of small-sized corydalis tuber harvester ［D］. Nanjing： Nanjing Agricultural University， 2013.

［2］ 魏萬(wàn)行. 小型元胡收獲機(jī)關(guān)鍵部件的設(shè)計(jì)與研究［D］. 漢中： 陜西理工大學(xué)， 2022.

Wei Wanhang. Design and investigation of key components of small rhizoma corydalis harvesting machine ［D］. Hanzhong： Shaanxi University of Technology， 2022.

［3］ 胡國(guó)明. 顆粒系統(tǒng)的離散元素法分析仿真［M］. 武漢： 武漢理工大學(xué)出版社， 2010.

［4］ Horabik J， Molenda M. Parameters and contact models for DEM simulations of agricultural granular materials： A review ［J］. Biosystems Engineering， 2016， 147： 206-225.

［5］ González-Montellano C， Llana D F， Fuentes J M， et al. Determination of the mechanical properties of corn grains and olive fruits required in DEM simulations ［C］. 2011 Louisville， Kentucky， August 7-10， 2011. American Society of Agricultural and Biological Engineers， 2011： 1.

［6］ Martin C L， Bouvard D， Shima S. Study of particle rearrangement during powder compaction by the discrete element method ［J］. Journal of the Mechanics and Physics of Solids， 2003， 51（4）： 667-693.

［7］ 王憲良， 胡紅， 王慶杰， 等. 基于離散元的土壤模型參數(shù)標(biāo)定方法［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2017， 48（12）： 78-85.

Wang Xianliang， Hu Hong， Wang Qingjie， et al. Calibration method of soil contact characteristic parameters based on DEM theory ［J］. Transactions of the Chinese Society for Agricultural Machinery， 2017， 48（12）： 78-85.

［8］ Santos K G， Campos A V P， Oliveira O S， et al. DEM simulations of dynamic angle of repose of acerola residue： A parametric study using a response surface technique ［J］. Blucher Chemical Engineering Proceedings， 2015， 1（2）： 11326-11333.

［9］ Wang L， Li R， Wu B， et al. Determination of the coefficient of rolling friction of an irregularly shaped maize particle group using physical experiment and simulations ［J］. Particuology， 2018， 38： 185-195.

［10］ 張榮芳， 焦偉， 周紀(jì)磊， 等. 不同填充顆粒半徑水稻種子離散元模型參數(shù)標(biāo)定［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2020， 51（S1）： 227-235.

Zhang Rongfang， Jiao Wei， Zhou Jilei， et al. Parameter calibration and experiment of rice seeds discrete element model with different filling particle radius ［J］. Transactions of the Chinese Society for Agricultural Machinery， 2020， 51（S1）： 227-235.

［11］ 張克平， 侯傳凱， 孫步功， 等. 收獲期豌豆籽粒離散元仿真參數(shù)標(biāo)定［J］. 干旱地區(qū)農(nóng)業(yè)研究， 2022， 40（6）：286-294.

Zhang Keping， Hou Chuankai， Sun Bugong， et al. Discrete element simulation parameter calibration of pea grains at harvest time ［J］. Agricultural Research in the Arid Areas， 2022， 40（6）： 286-294.

［12］ 劉凡一， 張艦， 李博， 等. 基于堆積試驗(yàn)的小麥離散元參數(shù)分析及標(biāo)定［J］. 農(nóng)業(yè)工程學(xué)報(bào)， 2016， 32（12）： 247-253.

Liu Fanyi， Zhang Jian， Li Bo， et al. Calibration of parameters of wheat required in discrete element method simulation based on repose angle of particle heap ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2016， 32（12）： 247-253.

［13］ 侯占峰， 戴念祖， 陳智， 等. 冰草種子物性參數(shù)測(cè)定與離散元仿真參數(shù)標(biāo)定［J］. 農(nóng)業(yè)工程學(xué)報(bào)， 2020， 36（24）： 46-54.

Hou Zhanfeng， Dai Nianzu， Chen Zhi， et al. Measurement and calibration of physical property parameters for agropyron seeds in a discrete element simulation ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2020， 36（24）： 46-54.

［14］ 賈富國(guó)， 韓燕龍， 劉揚(yáng)， 等. 稻谷顆粒物料堆積角模擬預(yù)測(cè)方法［J］. 農(nóng)業(yè)工程學(xué)報(bào)， 2014， 30（11）： 254-260.

Jia Fuguo， Han Yanlong， Liu Yang， et al. Simulation prediction method of repose angle for rice particle materials ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2014， 30（11）： 254-260.

［15］ 張德海， 李艷芹， 謝貴重， 等. 三維光學(xué)掃描技術(shù)逆向工程應(yīng)用研究［J］. 應(yīng)用光學(xué)， 2015， 36（4）： 519-525.

Zhang Dehai， Li Yanqin， Xie Guiguang， et al. Reverse engineering application of 3D optical scanning technology ［J］. Journal of Applied Optics， 2015， 36（4）： 519-525.

［16］ 劉瑞， 李衍軍， 劉忠軍， 等. 包衣玉米種子離散元參數(shù)分析與標(biāo)定［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2021， 52（S1）： 1-8， 18.

Liu Rui， Li Yanjun， Liu Zhongjun， et al. Analysis and Calibration of discrete element parameters of coated maize seeds ［J］. Transactions of the Chinese Society for Agricultural Machinery， 2021， 52（S1）： 1-8， 18.

［17］ 劉文政， 何進(jìn)， 李洪文， 等. 基于離散元的微型馬鈴薯仿真參數(shù)標(biāo)定［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2018， 49（5）： 125-135， 42.

Liu Wenzheng， He Jin， Li Hongwen， et al. Calibration of simulation parameters for potato minituber based on EDEM ［J］. Transactions of Society for Agricultural Machinery， 2018， 49（5）： 125-135， 142.

［18］ 于慶旭， 劉燕， 陳小兵， 等. 基于離散元的三七種子仿真參數(shù)標(biāo)定與試驗(yàn)［J］. 農(nóng)業(yè)機(jī)械學(xué)報(bào)， 2020， 51（2）： 123-132.

Yu Qingxu， Liu Yan， Chen Xiaobing， et al. Calibration and experiment of simulation parameters for panax notoginseng seeds based on DEM ［J］. Transactions of the Chinese Society for Agricultural Machinery， 2020， 51（2）： 123-132.

［19］ 葛宜元. 試驗(yàn)設(shè)計(jì)方法與Design-Expert軟件應(yīng)用［M］. 哈爾濱： 哈爾濱工業(yè)大學(xué)出版社， 2015.