基于堆積試驗的玉米包衣種子離散元參數(shù)標(biāo)定

李飛翔， 王鵬， 王云飛， 葛越鋒， 唐凱懌， 李得志

（洛陽智能農(nóng)業(yè)裝備研究院有限公司，河南 洛陽 471000）

氣力式排種器因具有排種精度高、對種子損傷少等優(yōu)點而被廣泛應(yīng)用。在實際排種作業(yè)中，種子處于氣固兩相環(huán)境中，運動情況較為復(fù)雜，離散元（discrete element method，DEM）與計算流體動力學(xué)（computational fluid dynamics，CFD）耦合法能夠有效分析排種器中流場的變化以及種子運動的受力情況，為排種器的結(jié)構(gòu)設(shè)計及參數(shù)優(yōu)化提供依據(jù)[1-3]。為提高玉米的產(chǎn)量及質(zhì)量，常常對玉米種子進(jìn)行包衣處理。種子包衣是將殺蟲劑、微生物肥料、著色劑等非種子材料包裹在種子表面，以提高種子抗逆性、抗病性，對農(nóng)業(yè)機(jī)械化、現(xiàn)代化及農(nóng)業(yè)可持續(xù)發(fā)展有著重要作用。研究人員基于離散元法對包衣種子已做了大量研究：閆聰杰[4]設(shè)計一種大豆包衣裝置，并基于離散元法對其包衣過程進(jìn)行分析；曲芳[5]針對現(xiàn)有大豆種子包衣攪拌裝置包衣不均、破碎率高的問題，基于離散元法設(shè)計一種攪拌葉片并優(yōu)化設(shè)計；胡建平等[6]以番茄磁粉包衣種子為對象，基于離散元法研究了種箱振動頻率、振幅對種群運動規(guī)律及種箱供種性能的影響。

DEM-CFD 常采用 EDEM 與 Fluent 軟件的耦合，要求EDEM 中顆粒體積小于Fluent 中網(wǎng)格體積。因此，EDEM 中常采用黏結(jié)顆粒模型（bonded particle model，BPM），用較小的顆粒黏結(jié)成物料顆粒模型[7]。BPM 需設(shè)置單位面積法向與切向剛度、極限法向與切向應(yīng)力、黏接半徑等參數(shù)[8]，現(xiàn)階段對于玉米種子的仿真參數(shù)設(shè)置研究較少。基于離散元的Hertz-Mindlin 模型以堆積角為響應(yīng)，王云霞等[9]以玉米顆粒種間靜摩擦與滾動摩擦系數(shù)為因素，設(shè)計試驗對其進(jìn)行標(biāo)定；王美美等[10]對玉米籽粒離散元仿真參數(shù)進(jìn)行標(biāo)定，為玉米籽粒離散元仿真提供參考；李秀清[11]基于離散元法對玉米粉碎過程進(jìn)行分析，采用測定的楔形玉米籽粒建模，借鑒玉米籽粒壓縮試驗數(shù)據(jù)，將極限法向應(yīng)力、極限切向應(yīng)力作為固定值，選取法向剛度、切向剛度、黏接半徑3個因素進(jìn)行標(biāo)定，得到玉米籽粒的仿真接觸參數(shù)。

為得到玉米包衣種子BPM 離散元仿真所需的精確的種間接觸參數(shù)，本研究將利用激光掃描儀得到的玉米包衣種子輪廓作為顆粒模板，采用BPM 以小顆粒黏結(jié)方式生成玉米包衣種子元顆粒模型，以堆積角作為響應(yīng)值，基于Isight 軟件的近似模型與DOE 聯(lián)合模塊，通過因素顯著性篩選、最優(yōu)逼近以及響應(yīng)面優(yōu)化對玉米包衣種子的仿真接觸參數(shù)進(jìn)行標(biāo)定，以期為玉米包衣種子氣力式排種器的兩相耦合仿真提供參考。

1 材料與方法

1.1 試驗材料獲取

供試玉米品種為鄭單958，純度＞96%，凈度＞90%，發(fā)芽率＞90%，含水率＞13.5%，由河南金博士種業(yè)股份有限公司提供。由于種子形狀多樣，將鄭單958 玉米種子分為大扁形、長扁形、類球形3類（圖1），比例為2∶3∶5。

圖1 玉米包衣種子形狀Fig.1 Shape of corn coated seed

1.2 接觸模型

為方便與Fluent 軟件耦合，EDEM 顆粒間接觸選用BPM。該模型接觸情況下，黏結(jié)鍵可以抵抗切向和法向運動，直至達(dá)到最大的法向和切向極限應(yīng)力為止。黏結(jié)作用產(chǎn)生后，元顆粒所受的力和力矩設(shè)置為0，并根據(jù)每個時間步長進(jìn)行增量調(diào)整[12-13]。

式中，δFn為法向力增量（N）；δFτ為切向力增量（N）；δMn為法向力矩增量（N·m）；δMτ為切向力矩增量（N·m）；Y為基礎(chǔ)區(qū)域面積（m2）；Rb為黏結(jié)半徑（m）；J為截面極慣性矩；Sn為法向剛度（N·m-3）；Sτ為切向剛度（N·m-3）；vn為法向速度（m·s-1）；vτ為切向速度（m·s-1）；wn為法向角速度（r·s-1）；wτ為切向角速度（r·s-1）；δt為時間步長（s）。

當(dāng)法向和切向剪應(yīng)力超過預(yù)定值σmax和τmax時，黏結(jié)就會斷裂。

黏結(jié)應(yīng)力與剛度值的不同設(shè)置代表了不同的材料屬性，高的剛度值將產(chǎn)生高的結(jié)合力和應(yīng)力。小顆粒在黏結(jié)之前，根據(jù)Hertz-Mindlin 模型相互作用，當(dāng)達(dá)到鍵形成時間時，被定義的小顆粒將會黏結(jié)在一起形成元顆粒。

1.3 堆積角測定方法

依據(jù)已有對堆積角的研究[14-18]，以漏斗法對玉米包衣種子的堆積角進(jìn)行測定，如圖2所示，測定5次取其平均值，測得鄭單958玉米包衣種子的堆積角為29.56o±0.27o。由于標(biāo)定參數(shù)為顆粒間接觸，為消除顆粒與部件的影響，采用亞克力材質(zhì)加工落料底盤，底盤直徑（L）為150 mm，底盤凹面深度為30 mm。待底盤落滿種子后，繼續(xù)落種；待玉米包衣種子溢出底盤后，觀察落料中顆粒堆積高度的變化；待堆積高度無顯著性變化時，用鋼尺測定包衣種子的堆積高度（h），由式（9）計算玉米種子堆積角。

圖2 堆積試驗Fig.2 Stacking test

1.4 仿真試驗

玉米種子為非規(guī)則形狀顆粒，為精確建立種子輪廓模型以提高標(biāo)定精準(zhǔn)性，選擇不同種類形狀的種子，采用海克斯康Tango-S手持式激光掃描儀對種子的外輪廓進(jìn)行掃描。得到種子的點云數(shù)據(jù)后，通過CATIA 軟件對種子點云數(shù)據(jù)進(jìn)行處理，最終得到玉米包衣種子的輪廓模型。

仿真模型采用與實際試驗?zāi)Ｐ拖嗤叽纾凑蘸Y選的大扁形、長扁形、類球形種子2∶3∶5的比例生成顆粒。顆粒堆積仿真試驗如圖3 所示，為節(jié)省仿真時間，先以圓柱底盤為虛擬工廠生成顆粒，待底盤表面落滿顆粒后，導(dǎo)出仿真時間為0時的文件；再以漏斗為虛擬工廠，進(jìn)行堆積角仿真試驗，待底盤溢出一定數(shù)量顆粒后，停止生成顆粒。采用后處理自帶的量角器，測量堆積角。

圖3 玉米種子顆粒堆積的模擬仿真Fig.3 Simulation test model for repose angle

參照玉米種子BPM 參數(shù)的設(shè)置[19-22]，本研究各參數(shù)取值范圍如表1 所示。仿真本征參數(shù)設(shè)定為玉米種子泊松比0.4，密度1 197 kg·m-3，剪切模量1.37×108Pa[23-26]。

表1 離散元仿真參數(shù)Table 1 Parameter in DEM simulation

1.5 標(biāo)定方法

以玉米包衣種子實際堆積角為響應(yīng)值，基于Isight 軟件的近似模型和DOE 聯(lián)合模塊，通過因素顯著性篩選、最優(yōu)逼近以及響應(yīng)面優(yōu)化對玉米包衣種子的仿真接觸參數(shù)進(jìn)行標(biāo)定。本研究的Plackett-Burman 設(shè)計以玉米包衣種子堆積角為響應(yīng)值對BPM 參數(shù)的顯著性進(jìn)行篩選。低水平設(shè)定為最初原始水平，高水平設(shè)為低水平的2倍，為方便對試驗誤差進(jìn)行分析，同時設(shè)定虛擬參數(shù)，試驗參數(shù)如表2所示。利用Isight軟件的近似模型與DOE模塊聯(lián)合對Plackett-Burman結(jié)果進(jìn)行方差分析，運行前，添加ANOVA Table、主效應(yīng)圖、Pareto圖對各參數(shù)的顯著性進(jìn)行分析。Plackett-Burman試驗后，以對玉米包衣種子堆積角影響顯著的參數(shù)進(jìn)行最陡爬坡，以便能快速進(jìn)入到最優(yōu)值的附近區(qū)域。最陡爬坡試驗根據(jù)Plackett-Burman試驗所得的回歸系數(shù)來確定爬坡步長。

2 結(jié)果與分析

2.1 輪廓模型掃描分析

圖4 為玉米包衣種子掃描的輪廓模型，可以看出，通過點云處理后得到的不同種類種子輪廓與實際種子輪廓一致，能夠用于仿真分析。

2.2 種子模型填充分析

通過激光掃描得到種子輪廓后，導(dǎo)入EDEM軟件中作為顆粒模板。結(jié)合玉米氣力式排種器研究進(jìn)展[27-28]，此處以0.5 mm 的小球顆粒對玉米種子顆粒模板進(jìn)行填充，如圖5 所示。對顆粒模板進(jìn)行上半部透明化顯示，填充的顆粒與模板貼合程度良好，通過小顆粒間形成黏結(jié)鍵建立玉米種子元顆粒黏結(jié)模型，建立的顆粒模型與玉米包衣種子外形相似。

圖5 玉米包衣種子仿真模型Fig.5 Simulation model of corn coated seeds

2.3 參數(shù)標(biāo)定結(jié)果分析

2.3.1 Plackett-Burman 試驗結(jié)果分析 Plackett-Burman 設(shè)計及結(jié)果如表3 所示，堆積角范圍在23.03°～30.31°之間，各仿真值沒有特別跳躍性情況。

表3 Plackett-Burman 試驗設(shè)計及結(jié)果Table 3 Design and results of Plackett-Burman test

2.3.2 Plackett-Burman 試驗顯著性分析 主效應(yīng)、Pareto結(jié)果如圖6所示，可清晰觀察各個因素對響應(yīng)的影響趨勢；ANOVA結(jié)果如表4所示，可知對玉米包衣種子堆積角影響極其顯著的參數(shù)為：玉米種子-玉米種子靜摩擦系數(shù)、法向剛度、切向剛度。其余參數(shù)結(jié)合已有研究進(jìn)展[29-30]取值（玉米種子-玉米種子恢復(fù)系數(shù)0.18、玉米種子-玉米種子滾動摩擦系數(shù) 0.05、法向應(yīng)力 1.1×107N·m-3、切向應(yīng)力4.1×106N·m-3、黏結(jié)半徑0.7 mm）進(jìn)行最陡爬坡以及響應(yīng)面試驗設(shè)計。

表4 Plackett-Burman 試驗參數(shù)顯著性分析Table 4 Analysis of significance of parameters in Plackett-Burman test

圖6 玉米種子對顆粒堆積響應(yīng)分析Fig.6 Response analysis of corn seed to grain accumulation

2.3.3 最陡爬坡試驗結(jié)果分析 最陡爬坡試驗設(shè)計及結(jié)果如表5所示，2號水平仿真堆積角與試驗值的誤差最小。由此選取2 號水平為中心點，1號、3號水平為低、高水平進(jìn)行后續(xù)優(yōu)化設(shè)計。

表5 最陡爬坡試驗及結(jié)果Table 5 Design and results of steepest ascent test

2.3.4 Box-Behnken 試驗結(jié)果分析 根據(jù)最陡爬坡試驗結(jié)果，選取顯著性參數(shù)進(jìn)行3個水平試驗設(shè)計。選取3個中心點對誤差進(jìn)行評估，結(jié)果如表 6 所示，基于 Isight 軟件 RSM 優(yōu)化模塊，建立3個參數(shù)與堆積角間的二階回歸方程。

表6 Box-Behnken 試驗及結(jié)果Table 6 Design and results of Box-Behnken test

Box-Behnken 試驗方差分析結(jié)果如表7 所示，玉米包衣種間靜摩擦系數(shù)(H)、法向剛度(J)、切向剛度（K）、法向剛度-切向剛度（JK）、靜摩擦系數(shù)-法向剛度（HJ）、靜摩擦系數(shù)二次項（H2）以及切向剛度二次項（K2）的P值＜0.05，表明這些參數(shù)對堆積角的影響顯著。失擬項P=0.093＞0.05，變異系數(shù)CV=0.801%，表明模型良好，試驗有較高的可靠性。決定系數(shù)R2=0.997 4，校正決定系數(shù)RAdj2=0.992 8，預(yù)測決定系數(shù)RPre2=0.958 9，表明模型能真實地反映實際情況。試驗精密度Adep precision=46.973，表明模型良好的精確度。

表7 Box-Behnken 試驗方差分析Table 7 ANOVA of quadraticl of Box-Behnken test

2.3.5 回歸模型交互效應(yīng)分析 根據(jù)優(yōu)化回歸模型方差分析結(jié)果可知，種間靜摩擦系數(shù)-法向剛度（HJ）以及法向剛度-切向剛度（JK）交互項對玉米包衣種子堆積角影響顯著（P＜0.05）。在玉米種間靜摩擦系數(shù)（H）為0.2 以及切向剛度（K）為0.5 情況下，應(yīng)用Isight 軟件對靜摩擦系數(shù)-法向剛度（HJ）以及法向剛度-切向剛度（JK）交互效應(yīng)的響應(yīng)曲面進(jìn)行繪制，如圖7所示，可以直觀反映交互項對堆積角的影響。由HJ曲面可知，相對玉米包衣種間靜摩擦系數(shù)（H），法向剛度（J）的效應(yīng)面曲線比較陡，表明其對堆積角影響較為顯著。由JK曲面可知，相對于法向剛度（J），切向剛度（K）的效應(yīng)面曲線比較陡，表明其對堆積角影響較為顯著。

圖7 HJ與JK的交互效應(yīng)Fig.7 Interaction effect of HJ and JK

2.3.6 試驗優(yōu)化結(jié)果分析 應(yīng)用Isight軟件RSM優(yōu)化模塊，以玉米包衣種子實際測定的堆積角29.56o為目標(biāo)，對回歸方程尋優(yōu)求解得到玉米包衣種間靜摩擦系數(shù)為0.269，法向剛度為2.54×108N·m-3，切向剛度為5.93×107N·m-3。用所得最佳參數(shù)組合進(jìn)行堆積仿真試驗，結(jié)果如圖8 所示。仿真試驗所得堆積角為29.85o，與實際試驗值29.56o誤差為0.98%，表明最佳參數(shù)組合可行。

圖8 仿真與實際試驗對比Fig.8 Comparison between simulation and actual test

3 討論

現(xiàn)有堆積試驗多采取物料直接堆積于材料表面，實際仿真標(biāo)定中，物料與部件間的參數(shù)設(shè)置對堆積結(jié)果影響較大，為減小試驗誤差，堆積試驗采用凹槽底盤，使得物料落滿底盤后，在物料表面形成堆積。本研究選用鄭單958 玉米包衣種子，在對種子進(jìn)行分類篩選后，通過激光掃描儀對輪廓較好的種子進(jìn)行掃描得到點云數(shù)據(jù)，采用逆向工程技術(shù)，通過CATIA 軟件對點云數(shù)據(jù)處理，最終得到玉米種子仿真模型。

根據(jù)主效應(yīng)和Pareto 圖分析，篩選出對堆積角影響顯著的參數(shù)：玉米種子-玉米種子靜摩擦系數(shù)、法向剛度、切向剛度。基于RSM 模塊，建立3個顯著性參數(shù)與堆積角間的二次回歸模型，由方差分析結(jié)果可知，交互項靜摩擦系數(shù)-法向剛度（HJ）、法向剛度-切向剛度（JK）以及3個顯著性參數(shù)的二次項（H2）對玉米包衣種子堆積角影響極顯著。

以玉米包衣種子實際測定堆積角為目標(biāo)，對回歸方程進(jìn)行尋優(yōu)求解，得到參數(shù)的最佳組合為：種間靜摩擦系數(shù)0.269，法向剛度2.54×108N·m-3，切向剛度5.93×107N·m-3。用所得最佳參數(shù)組合進(jìn)行仿真試驗，與實際值誤差為0.98%，仿真與實際試驗值差異不顯著，表明標(biāo)定所得參數(shù)可用于玉米包衣種子離散元BPM 仿真，為玉米氣力式排種器的耦合模擬提供參考。