基于離散元仿真參數(shù)的微型薯物料測定

余參參，段宏兵,，蔡興奎，徐濤 ，姚飛虎， 陳志惠，嚴(yán)福勇

1.華中農(nóng)業(yè)大學(xué)工學(xué)院，武漢 430070; 2.農(nóng)業(yè)農(nóng)村部馬鈴薯生物學(xué)與生物技術(shù)重點(diǎn)實(shí)驗(yàn)室，武漢 430070

我國馬鈴薯種植面積與產(chǎn)量均居世界首位，但單產(chǎn)與世界平均水平仍有差距。實(shí)踐證明，采用脫毒微型馬鈴薯(簡稱“微型薯”)代替切塊薯增產(chǎn)效果顯著，單位面積產(chǎn)量一般可增加30%～50%，高的可增產(chǎn)3～4倍[1-2]。通過對(duì)脫毒種薯的擴(kuò)繁而形成優(yōu)質(zhì)商品種薯，是降低種薯成本的有效手段?，F(xiàn)階段適用于脫毒種薯的機(jī)械化播種機(jī)具較少，因而造成脫毒種薯的種植過程機(jī)械化水平較低。同時(shí)由于脫毒種薯塊莖的形狀差異較大，農(nóng)業(yè)生產(chǎn)中急需微型薯播種與分級(jí)技術(shù)與裝備[3]。隨著數(shù)字仿真技術(shù)的發(fā)展，離散元應(yīng)用在作物與作物之間、作物與機(jī)械設(shè)備之間的仿真研究漸多。研究表明，機(jī)械設(shè)備的離散元仿真分析結(jié)果與種子物理特性存在較大聯(lián)系[4-5]。目前部分學(xué)者針對(duì)種子的物理特征參數(shù)和部分接觸力學(xué)參數(shù)進(jìn)行了標(biāo)定，結(jié)果均表明離散元仿真的平均誤差值較小，基本與試驗(yàn)一致，可以縮短研究周期[6-7]。由于微型薯顆粒個(gè)體差異較大，品種、含水率等對(duì)物料參數(shù)的測定結(jié)果影響巨大，因此，本研究以3種不同品種的微型薯為研究對(duì)象，以離散元軟件EDEM中微型薯機(jī)械設(shè)備研究所需的相關(guān)仿真參數(shù)為目標(biāo)，對(duì)微型薯的物理力學(xué)特征參數(shù)和接觸力學(xué)參數(shù)進(jìn)行測定分析，旨在為微型薯相關(guān)機(jī)械設(shè)備在EDEM中的仿真參數(shù)的設(shè)置提供參考。

1 材料與方法

1.1 試驗(yàn)材料

由于不同品種的微型薯差異較大，因此，試驗(yàn)用薯選取種植相對(duì)廣泛且具有一定代表性的“大西洋”“中薯5號(hào)”和“華薯1號(hào)”，3種種薯均取自華中農(nóng)業(yè)大學(xué)農(nóng)業(yè)農(nóng)村部馬鈴薯生物學(xué)與生物技術(shù)重點(diǎn)實(shí)驗(yàn)室種植基地。采樣方法為：選取種植基地某處大棚內(nèi)長勢較為穩(wěn)定的區(qū)域，隨機(jī)選取半徑2.0 m范圍采挖3個(gè)品種的種薯，直至各收集500個(gè)種薯為止，現(xiàn)場收集的種薯經(jīng)標(biāo)記后進(jìn)行測量并記錄數(shù)據(jù)，測量后置于冷庫儲(chǔ)存?zhèn)溆谩y試結(jié)果顯示，3種種薯含水率分別為71.75%～83.65%、75.45%～82.58%和74.49%～80.48%。

1.2 微型薯特征參數(shù)測定

1)形態(tài)參數(shù)與質(zhì)量的測量。對(duì)采集的“大西洋”“中薯5號(hào)”“華薯1號(hào)”樣品各500個(gè)進(jìn)行三軸尺寸、曲率半徑和質(zhì)量測量。質(zhì)量測定采用“苦竹”電子秤，精度0.01 g；采用“艾瑞澤”數(shù)顯游標(biāo)卡尺測量種薯三軸尺寸，可測量三軸極限尺寸，精度0.01 mm；曲率半徑選用“天目”多功能數(shù)顯半徑測試儀，其可針對(duì)不同的圓弧形面和測量精度選用不同跨度的量爪，精確讀出種薯的外圓弧半徑，精度0.005 mm。三軸尺寸和曲率半徑的測量均分別選取長、寬、高3個(gè)方向進(jìn)行(圖1)。

X為測量時(shí)長度方向，Y為寬度方向，Z為高度方向，測量三軸尺寸時(shí)分別選取每個(gè)方向上的2個(gè)頂點(diǎn)進(jìn)行測量，頂點(diǎn)選取如1、2所示；測量曲率半徑時(shí)，將量爪的中心分別與所測種薯3個(gè)方向的中心對(duì)齊后讀取數(shù)據(jù)。X is the length direction, Y is the width direction, and Z is the height direction. When measuring triaxial dimensions, select two vertices in each direction for measurement. The vertices are selected as sign 1 and 2; when measuring the radius of curvature, measure data is recorded after the center of the paw is aligned with the center of the three directions.

微型薯的三軸尺寸分布統(tǒng)計(jì)見圖2，曲率半徑與質(zhì)量概率分布統(tǒng)計(jì)情況見圖3。圖2顯示，3個(gè)品種的微型薯尺寸多分布在寬度10～28 mm、厚度10～26 mm、長度13～40 mm，均占比90%以上，其中寬度主要集中在12～18 mm，占比60%以上，厚度集中在12～18 mm，占比70%以上，長度集中在16～25 mm，占比50%左右。

圖3顯示，“大西洋”3個(gè)方向的曲率半徑相差較大，符合程度較低，即3個(gè)方向的外圓弧半徑差異較大，長度方向的尺寸大于寬度和高度，且高度方向和寬度方向尺寸較為接近，整體形狀呈長橢圓形?！叭A薯1號(hào)”和“中薯5號(hào)”3個(gè)方向的曲率半徑符合程度較高，即各個(gè)方向上的彎曲程度相近，三軸尺寸差別相對(duì)較小，相對(duì)來說外形偏圓，但“中薯5號(hào)”曲率半徑較“華薯1號(hào)”符合程度更高，圓度更高，截面規(guī)則?！叭A薯1號(hào)”雖外形近圓，但其表面有小凹坑，截面形狀有邊角，近似橢圓形。不同品種的微型薯質(zhì)量分布呈現(xiàn)出相似的變化規(guī)律，3個(gè)品種的種薯峰值均處于3～5 g，質(zhì)量在10 g以下的種薯均占比90%以上；“中薯5號(hào)”質(zhì)量分布范圍最廣，最小的不足1 g，最大的超過40 g。綜上，同一品種的微型薯個(gè)體差異大，但基本符合正態(tài)分布；品種間形狀差異較大，總體變化規(guī)律相似。

圖2 三軸尺寸概率分布

A:“大西洋” “Daxiyang”; B:“中薯5號(hào)” “Zhongshu No.5”; C：“華薯1號(hào)” “Huashu No.1”; D：微型薯質(zhì)量 Weight.

2)密度測量。由于微型薯表面較為光滑，體積較小，無疏松孔隙，液體難以浸潤，可采用量筒法測量相關(guān)參數(shù)。隨機(jī)選取3種微型薯各100個(gè)，選用量程1 500 mL的量筒和精度為0.01 g的電子秤，對(duì)3個(gè)品種的微型薯進(jìn)行體積和質(zhì)量的測量，試驗(yàn)重復(fù)5次，測量結(jié)果取平均值。 “大西洋”“中薯5號(hào)”“華薯1號(hào)”的密度分別為1 068、1 044、1 056 kg/m3。

3)剪切模量、泊松比的測定。本研究討論在不超過屈服載荷下微型薯的彈性變形能力，采用萬能試驗(yàn)機(jī)測量種薯破裂瞬間的相關(guān)數(shù)據(jù)，進(jìn)而計(jì)算剪切模量和泊松比。試驗(yàn)儀器為瑞格爾生產(chǎn)制造的RGB-3005型萬能試驗(yàn)機(jī)，可施加0～5 kN載荷。根據(jù)需要選用90 mm平板壓頭，通過該壓縮儀器獲得計(jì)算所需的橫向變形量、縱向變形量、試驗(yàn)載荷力等數(shù)據(jù)參數(shù)，加載速度固定為10 mm/min。試驗(yàn)用種薯為測量物理力學(xué)特征參數(shù)時(shí)所用種薯，放置1個(gè)月后進(jìn)行壓縮試驗(yàn)測量(圖4)。由于微型薯形狀大小差異較大，考慮該差異可能對(duì)試驗(yàn)所測結(jié)果產(chǎn)生影響，試驗(yàn)時(shí)根據(jù)微型薯種植要求按質(zhì)量小于3 g、3～5 g和6～10 g標(biāo)準(zhǔn)分組進(jìn)行試驗(yàn)。

圖4 壓縮試驗(yàn)

剪切模量是剪切變形時(shí)剪切應(yīng)力與剪切應(yīng)變的比值[8]，其公式為：

(1)

式(1)中：E為彈性模量，MPa；μ為泊松比；G為剪切模量。由式(1)知，剪切模量可通過對(duì)彈性模量、泊松比的數(shù)值計(jì)算得到，泊松比、彈性模量參考文獻(xiàn)[8]方法測定。

將3種微型薯按照質(zhì)量差異分9組進(jìn)行試驗(yàn)，每組測量50個(gè)，分別得9組剪切模量、泊松比和彈性模量。對(duì)各個(gè)品種的微型薯按質(zhì)量分級(jí)段進(jìn)行試驗(yàn)，結(jié)果發(fā)現(xiàn)彈性模量、泊松比與質(zhì)量無關(guān)。如表1所示，僅列出了質(zhì)量不一的“大西洋”的試驗(yàn)結(jié)果，試驗(yàn)結(jié)果統(tǒng)計(jì)僅取平均值。試驗(yàn)結(jié)果顯示，同一品種間不同質(zhì)量的微型薯，其彈性模量、泊松比及剪切模量無明顯變化，但不同品種間差異較明顯。

表1 微型薯的剪切模量和泊松比 Table 1 Shear modulus and Poisson’s ratio for potato minitubers

1.3 微型薯接觸力學(xué)參數(shù)的測定

測定3種微型薯(“大西洋”“中薯5號(hào)”和“華薯1號(hào)”)分別與3種材料(鋼板、ABS和各自品種微型薯薯群)的碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動(dòng)摩擦系數(shù)。為消除微型薯的形狀和質(zhì)量影響，挑選大小相近、質(zhì)量相當(dāng)?shù)姆N薯進(jìn)行試驗(yàn)。當(dāng)微型薯薯群作為試驗(yàn)材料時(shí)需要挑選大小相近，形狀相同的若干較大微型薯粘結(jié)在一起。

1)碰撞恢復(fù)系數(shù)測定。如果物體與固定平面碰撞,則碰撞前后固定平面的速度均為0，此時(shí)碰撞恢復(fù)系數(shù)可表示為：

Cr=v1/v0

(2)

式(2)中：v0為碰撞前微型薯的速率，m/s；v1為碰撞后微型薯的速率，m/s。碰撞試驗(yàn)裝置如圖5所示，以“華薯1號(hào)”與3種材料的碰撞為例。

試驗(yàn)全程采用高速攝影拍攝，頻率為1 000 fps/s。微型薯從同一高度(400 mm)由靜止開始做自由落體運(yùn)動(dòng)，與碰撞材料接觸后碰撞產(chǎn)生回彈力使之向上運(yùn)動(dòng)，整個(gè)運(yùn)動(dòng)過程能量守恒，碰撞恢復(fù)系數(shù)可由公式(3)得：

(3)

式(3)中：g為重力加速度，9.8 m/s2；h2為碰撞后回彈高度，m；h1為碰撞前下落高度，m。當(dāng)碰撞材料為微型薯薯群時(shí)，下落的種薯碰撞到微型薯薯群上任一種薯為有效試驗(yàn)，種薯下落到薯與薯之間的空隙時(shí)無效，不計(jì)入試驗(yàn)。每個(gè)品種的微型薯與3種碰撞材料均取20次有效試驗(yàn)，每次試驗(yàn)重復(fù)5次。試驗(yàn)結(jié)果顯示：“大西洋”“中薯5號(hào)”“華薯1號(hào)”與鋼板、ABS和各自品種微型薯薯群間的碰撞恢復(fù)系數(shù)分別為：0.577、0.482、0.449；0.564、0.480、0.460；0.562、0.459、0.440。

A:“華薯1號(hào)”與ABS碰撞 Collision between “Huashu No.1” and ABS; B:“華薯1號(hào)”與“華薯1號(hào)”碰撞 Collision between “Huashu No.1” and “Huashu No.1”; C:“華薯1號(hào)”與鋼板碰撞 Collision between “Huashu No.1” and steel plate;1.碰撞材料 Collision material；2.刻度尺 Scale；3.碰撞的微型薯 Collision of minitubers; 4.試驗(yàn)臺(tái)架 Experiment shelf.

2)靜摩擦系數(shù)測定。微型薯置于傾斜材料上，重力分解為與斜面平行的力F1和與斜面垂直的力F2，最大靜摩擦力f與F1大小相等，方向相反，即F1=f=mgsinα，此時(shí)斜面的傾角與最大靜摩擦存在以下關(guān)系：

μ=f/F2=mgsinα/(mgcosα)=tanα

(4)

因此，可采用靜摩擦測量裝置進(jìn)行測量。測量時(shí)將3種材料依次置于斜板上固定，微型薯(為消除滾動(dòng)摩擦的影響，挑選各自品種代表且大小相近的4個(gè)種薯粘連在一起為1組)放置在固定好的材料上，放置位置相同，均為斜板300 mm處。緩慢調(diào)節(jié)斜板角度，直至微型薯和材料存在滑動(dòng)趨勢時(shí)停止角度調(diào)節(jié)，記錄此時(shí)底板與材料之間高度h，由α=arctan(h/300)計(jì)算得到傾斜角度，即為靜摩擦系數(shù)。每個(gè)品種的種薯選取10組進(jìn)行試驗(yàn)，每組試驗(yàn)重復(fù)20次取平均值。結(jié)果顯示：“大西洋”“中薯5號(hào)”“華薯1號(hào)”與鋼板、ABS和各自品種微型薯薯群間的靜摩擦系數(shù)分別為：0.469、0.493、0.434;0.386、0.390、0.356；0.426、0.411、0.38。

3)滾動(dòng)摩擦系數(shù)測定。當(dāng)點(diǎn)或線接觸的兩物體在法向和轉(zhuǎn)矩作用下產(chǎn)生相對(duì)滾動(dòng)時(shí)，其表面會(huì)產(chǎn)生阻礙相對(duì)運(yùn)動(dòng)趨勢的形變，可用滾動(dòng)摩擦力衡量[9]。滾動(dòng)摩擦系數(shù)的測定可采用斜面滾動(dòng)法。微型薯的滾動(dòng)距離和滾動(dòng)時(shí)間經(jīng)高速攝像記錄后，由式(5)計(jì)算相應(yīng)的滾動(dòng)摩擦系數(shù)。試驗(yàn)裝置同靜摩擦試驗(yàn)裝置，材料固定在斜板上，微型薯放置在材料上同一位置由靜止釋放，以0 m/s開始在材料上運(yùn)動(dòng)。高速攝影記錄滾動(dòng)過程，利用幀數(shù)變化來計(jì)算滾動(dòng)時(shí)間。設(shè)置相機(jī)幀數(shù)200 fps/s，從微型薯開始滾動(dòng)的幀數(shù)算起，直至在斜面的滾動(dòng)距離為200 mm時(shí)的幀數(shù)為止，幀數(shù)差即為滾動(dòng)時(shí)間，滾動(dòng)距離200 mm。微型薯滾動(dòng)的過程中滿足：

(5)

式(5)中：μ為滾動(dòng)摩擦系數(shù)；θ為斜面夾角，(°)；S為滾動(dòng)距離，m；g為重力加速度，9.8 m/s2；t為滾動(dòng)時(shí)間，s。

微型薯的彈性模量較鋼板、ABS小，微型薯運(yùn)動(dòng)時(shí)可忽略微觀變形和彈性滯后，只有可能產(chǎn)生些許滑動(dòng)影響試驗(yàn)結(jié)果。當(dāng)傾斜角度較大時(shí)，滑動(dòng)影響越大。為了降低這種影響，對(duì)3種微型薯進(jìn)行預(yù)滾動(dòng)，選取3種微型薯都能滾動(dòng)的最小角度(16°)進(jìn)行滾動(dòng)試驗(yàn)；每個(gè)品種進(jìn)行30次重復(fù)滾動(dòng)試驗(yàn)取平均值。結(jié)果顯示：“大西洋”“中薯5號(hào)”和“華薯1號(hào)”與鋼板、ABS間的滾動(dòng)摩擦系數(shù)分別為： 0.119、0.120；0.067、0.089；0.076、0.082。

1.4 仿真試驗(yàn)

1)仿真顆粒模型設(shè)置。根據(jù)各微型薯的物理力學(xué)特征參數(shù)建立微型薯顆粒模型，各品種按比例產(chǎn)生不同尺寸的種薯顆粒(圖6)。

2)堆積角仿真參數(shù)。由于微型薯之間的滾動(dòng)摩擦系數(shù)影響因素太多，難以直接測量，而堆積角與物料的摩擦特性有很大關(guān)系[10-11]，因此，目前通用的方法是采用離散元軟件EDEM建立測量堆積角模型，并輸入微型薯間的靜摩擦系數(shù)、微型薯間的滾動(dòng)摩擦系數(shù)以及微型薯和ABS間滾動(dòng)摩擦系數(shù)進(jìn)行仿真試驗(yàn)，改變EDEM軟件中的各摩擦系數(shù)，得到不同的物料堆積模型，測得對(duì)應(yīng)物料堆積模型的堆積角。選取當(dāng)仿真結(jié)果的堆積角與實(shí)際試驗(yàn)中的堆積角近似相等時(shí)的滾動(dòng)摩擦系數(shù)作為物料間的實(shí)際滾動(dòng)摩擦系數(shù)。

由于微型薯間滾動(dòng)摩擦系數(shù)較小，故前期仿真試驗(yàn)優(yōu)先確定其范圍為0.01～0.11，仿真時(shí)按照優(yōu)先確定的范圍將其等分設(shè)置為3個(gè)水平；微型薯間靜摩擦系數(shù)、微型薯與ABS間滾動(dòng)摩擦系數(shù)的設(shè)置前述試驗(yàn)測得的參數(shù)(表2)為基礎(chǔ)值上下浮動(dòng)各0.05。以堆積角θ為目標(biāo)值，在Design Expert軟件中分別對(duì)3種微型薯開展三因素三水平響應(yīng)面設(shè)計(jì)，每組試驗(yàn)17組，3組試驗(yàn)共計(jì)51組。每組試驗(yàn)重復(fù)3次取平均值。其余各仿真參數(shù)設(shè)置如表2所示。

A:“大西洋” “Daxiyang”; B:“華薯1號(hào)” “Huashu No. 1”; C:“中薯5號(hào)” “Zhongshu No.5”.

表2 堆積角仿真參數(shù)設(shè)置 Table 2 Simulation parameter setting for stacking

3)堆積角試驗(yàn)。堆積角試驗(yàn)采用抽板法(圖7)。仿真接觸模型選擇無滑移Hertz-Mindlin，材料為ABS，在軟件Geometry頁面中創(chuàng)建1個(gè)200 mm×150 mm×150 mm無蓋長方體，一個(gè)虛擬的四邊形蓋板作為種薯下落到長方體內(nèi)的基準(zhǔn)平面。顆粒工廠設(shè)為動(dòng)態(tài)，種薯的總質(zhì)量3.0 kg，根據(jù)種薯的實(shí)際幾何尺寸設(shè)置，種薯按照一定的數(shù)量比例以500個(gè)/s的速度填充長方體內(nèi)部，填充時(shí)間為5 s。待種薯填滿盒內(nèi)后，創(chuàng)建1個(gè)四邊形置于長方體底面充當(dāng)承接掉落種薯的地板，設(shè)置長方體的一側(cè)面以0.1 m/s的速度向上運(yùn)動(dòng)，使種薯能以較緩的速度從盒中自由散落在地板上，仿真時(shí)間設(shè)置為10 s，待種薯靜止后測量堆積角度。真實(shí)試驗(yàn)與軟件仿真模型相同，對(duì)3種微型薯進(jìn)行5次抽板法的重復(fù)試驗(yàn)，得平均值θ，見式(6)～(8)。

A:軟件仿真 Simulation; B:真實(shí)試驗(yàn) Actual test.

(6)

(7)

(8)

2 結(jié)果與分析

2.1 響應(yīng)面試驗(yàn)

響應(yīng)面試驗(yàn)設(shè)計(jì)結(jié)果見表3(以“大西洋”為例)，對(duì)試驗(yàn)結(jié)果進(jìn)行方差分析，發(fā)現(xiàn)自變量中“大西洋”間滾動(dòng)摩擦系數(shù)(B1)極顯著(P<0.000 1)，失擬項(xiàng)為0.038，該方程擬合較好，比較可靠，結(jié)果表明因素B1對(duì)堆積角的結(jié)果影響較大，“大西洋”間靜摩擦系數(shù)(A1)、“大西洋”與ABS間滾動(dòng)摩擦系數(shù)(C1)影響較小。為了得到合理的堆積角，仿真與試驗(yàn)結(jié)果均采用MATLAB對(duì)圖片進(jìn)行處理[12]，提取邊界后進(jìn)行直線擬合，獲取堆積角的真實(shí)角度(圖8)。

表3 “大西洋”微型薯仿真試驗(yàn)設(shè)計(jì)及結(jié)果 Table 3 Simulation test design and results of “Daxiyang”

A:堆積角 Stacking angle; B:提取邊界 Extraction boundary; C:直線擬合 Straight line fitting.

使用軟件的優(yōu)化功能分析確定符合的參數(shù)，以自變量A1、B1、C1的變化范圍為基礎(chǔ)，將實(shí)際堆積角作為響應(yīng)值(目標(biāo))對(duì)回歸模型求解(式6)。式(6)中，A、B、C分別為微型薯間靜摩擦系數(shù)、微型薯間滾動(dòng)摩擦系數(shù)、微型薯與ABS間滾動(dòng)摩擦系數(shù)；1、2、3分別為“大西洋”“中薯5號(hào)”和“華薯1號(hào)”，如A1為“大西洋”間靜摩擦系數(shù)。確定 “大西洋”“中薯5號(hào)”和“華薯1號(hào)”與各自品種間的滾動(dòng)摩擦系數(shù)分別為0.031、0.020和0.024。

2.2 試驗(yàn)結(jié)果驗(yàn)證

在幾何體模塊新建一個(gè)仿真模型(圖9A)，創(chuàng)建一個(gè)500 mm×500 mm的四邊形底面置于圓筒底部，一個(gè)無底無蓋圓筒200 mm(半徑)×350 mm(高)，將前述試驗(yàn)測得的物理力學(xué)特征參數(shù)、接觸力學(xué)參數(shù)和仿真得到的微型薯間的滾動(dòng)摩擦系數(shù)輸入EDEM軟件參數(shù)設(shè)置中，其余相關(guān)設(shè)置參見表2。設(shè)置微型薯顆粒的總質(zhì)量為2.0 kg，仿真時(shí)間10 s。待種薯填滿圓筒且穩(wěn)定后，設(shè)置圓筒以0.1 m/s的速度向上做線性運(yùn)動(dòng)，直至仿真結(jié)束，測得“大西洋”“中薯5號(hào)”“華薯1號(hào)”的仿真堆積角分別為23.377°、19.824°和20.913°。

實(shí)際試驗(yàn)(圖9B)與仿真試驗(yàn)方法一致，進(jìn)行5次重復(fù)試驗(yàn)后得“大西洋”“中薯5號(hào)”“華薯1號(hào)”的堆積角分別為23.778°、19.423°和21.314°，與仿真的堆積角試驗(yàn)結(jié)果相比，相對(duì)誤差分別為1.69%、2.06%、1.88%，表明采用此方法可得到較為相似的結(jié)果，與實(shí)際情況相符，即微型薯離散元顆粒參數(shù)可為EDEM模擬仿真微型薯顆粒與機(jī)械設(shè)備優(yōu)化提供參考。

A:仿真試驗(yàn) Simulation test; B 真實(shí)試驗(yàn) Actual test.

3 討 論

本研究測量了湖北地區(qū)常見3種微型薯(“大西洋”“中薯5號(hào)”“華薯1號(hào)”)的物理力學(xué)特征參數(shù)，包括其基本物理特征參數(shù)(三軸尺寸、質(zhì)量、曲率半徑)、密度、泊松比和剪切模量。測試結(jié)果表明，密度、剪切模量和泊松比與微型薯的基本物理特征參數(shù)相關(guān)性較小，但基本物理特征參數(shù)差異較大，品種差異較大。在設(shè)計(jì)、仿真計(jì)算時(shí)可以根據(jù)用戶需求選擇所需參數(shù)。3種微型薯和各自品種微型薯薯群之間的接觸力學(xué)參數(shù)較鋼板、ABS均低，3種微型薯和ABS的靜、動(dòng)摩擦系數(shù)均高于鋼板，但碰撞恢復(fù)系數(shù)則呈相反趨勢，ABS均低于鋼板。“大西洋”的接觸力學(xué)參數(shù)均高于“中薯5號(hào)”和“華薯1號(hào)”，3種材料中“大西洋”與鋼板的接觸力學(xué)參數(shù)最大，與ABS和微型薯薯群之間的接觸力學(xué)參數(shù)依次降低，表明其彈性恢復(fù)性能最好，摩擦力較大，相對(duì)不容易損傷，但流動(dòng)性較差，機(jī)械化播種性能差。“中薯5號(hào)”和“華薯1號(hào)”的接觸力學(xué)參數(shù)差異較小，“中薯5號(hào)”略大于“華薯1號(hào)”。不同品種對(duì)接觸力學(xué)參數(shù)影響大，在設(shè)計(jì)、仿真時(shí)要充分考慮品種差異，不能忽略。采用響應(yīng)曲面法設(shè)計(jì)試驗(yàn)，以靜摩擦系數(shù)、滾動(dòng)摩擦系數(shù)為試驗(yàn)變量，堆積角為試驗(yàn)指標(biāo)，根據(jù)試驗(yàn)因素與變量的二階回歸方程進(jìn)行方差分析，結(jié)果表明，各品種微型薯之間的滾動(dòng)摩擦系數(shù)差異極顯著。利用Design Expert軟件中的優(yōu)化模塊將試驗(yàn)測得的實(shí)際堆積角作為響應(yīng)值對(duì)得到的回歸模型求解，確定了“大西洋”“中薯5號(hào)”“華薯1號(hào)”與各自品種間的滾動(dòng)摩擦系數(shù)分別為0.031、0.020和0.024。通過圓筒提升試驗(yàn)對(duì)仿真得到的滾動(dòng)摩擦系數(shù)進(jìn)行驗(yàn)證，結(jié)果顯示仿真結(jié)果與試驗(yàn)結(jié)果的差異均小于3%，說明本研究所測定的微型薯顆粒實(shí)際參數(shù)可作為EDEM軟件中仿真參數(shù)使用，仿真結(jié)果誤差較小。