摘要：為保證離散元仿真模型的精確性，以玉米種子為研究對象，采用實際物理試驗與仿真模擬試驗相結(jié)合的方法，對玉米種子離散元模型進行參數(shù)標定。首先，利用實際物理試驗分別測量玉米種子外形尺寸及接觸參數(shù)，根據(jù)基本物性參數(shù)測定實際堆積角為24.95°。根據(jù)實際物理試驗獲取仿真模擬參數(shù)取值范圍，進行Plackett-Burman試驗設(shè)計，挑選顯著性參數(shù)：玉米-玉米靜摩擦系數(shù)、玉米-不銹鋼靜摩擦系數(shù)和玉米-不銹鋼滾動摩擦系數(shù)，通過最陡滑坡試驗，確定顯著性因素取值范圍。利用Box-Behenken試驗建立顯著性參數(shù)與堆積角之間的二次回歸模型并剔除不顯著參數(shù)對其進一步優(yōu)化，最終得到最優(yōu)數(shù)值模擬參數(shù)：玉米-玉米靜摩擦系數(shù)為0.273、玉米-不銹鋼靜摩擦系數(shù)為0.28、玉米-玉米滾動摩擦系數(shù)為0.06、玉米-玉米恢復(fù)系數(shù)為0.23、玉米-不銹鋼滾動摩擦系數(shù)為0.06、玉米-不銹鋼恢復(fù)系數(shù)為0.675，仿真試驗堆積角平均值為24.48°，相對誤差為1.88%，驗證了仿真模擬的準確性。

關(guān)鍵詞：玉米種子；離散元；參數(shù)標定；接觸參數(shù)；堆積角

中圖分類號：S229 文獻標志碼：A

文章編號：1002-1302（2024）18-0209-06

收稿日期：2023-09-11

基金項目：國家自然科學基金（編號：51872254）。

作者簡介：繆雨龍（2000—），男，江蘇如東人，碩士研究生，主要從事農(nóng)業(yè)機械方面的研究。E-mail：2511832234@qq.com。

通信作者：趙毅紅，博士，副教授，碩士生導師，研究方向為機電控制及智能機器人。E-mail： yihongzhao99@sina.com。

散狀物料的傳輸與轉(zhuǎn)載在礦業(yè)、建材、化工和農(nóng)業(yè)等領(lǐng)域都有廣泛應(yīng)用^［1］。隨農(nóng)業(yè)生產(chǎn)機械化效率的提高，散狀物料玉米的傳輸?shù)玫街匾?sup>［^2］。通過離散元軟件EDEM模擬玉米在轉(zhuǎn)載過程中的運動特性及物料間的相互作用，減少研發(fā)周期，降低研發(fā)成本。物料的物性參數(shù)對EDEM仿真模擬結(jié)果的可信度至關(guān)重要，包括物料與材料本身基本參數(shù)（泊松比、密度、剪切模量等）以及物料與材料之間的接觸參數(shù)（靜摩擦系數(shù)、滾動摩擦及恢復(fù)系數(shù)等）^［3-5］。近年來，許多學者對散狀物料離散元參數(shù)進行了研究。劉大為等采用物理試驗與仿真試驗相結(jié)合的方法，以堆積角為響應(yīng)值對玉針香父母本種子顆粒模型間的摩擦參數(shù)進行標定^［6］。鄭效帥等對油莎豆種子離散元參數(shù)進行仿真試驗標定，通過排種仿真及物理試驗獲得最佳參數(shù)組合^［7］。周濱等對顆粒鉀肥測量三軸尺寸和質(zhì)量，通過無底圓筒堆積物理試驗，獲得顆粒鉀肥的休止角，通過仿真驗證試驗效果與標定后仿真結(jié)果的一致性^［8］。

在玉米種子模型建立及離散元參數(shù)標定模擬方面也有少量文獻描述^［9-14］。本試驗以商用玉米為研究對象，結(jié)合物理試驗與仿真模擬，確定玉米種子的基本物理參數(shù)及接觸參數(shù)，建立離散元仿真模型。運用EDEM軟件仿真玉米種子堆積過程，利用Plackett-Burman試驗篩選、最陡滑坡試驗及響應(yīng)面試驗確定玉米種子最優(yōu)參數(shù)，為后期玉米運輸及轉(zhuǎn)載仿真提供參數(shù)依據(jù)。

1 材料與方法

1.1 試驗材料與儀器

選用2023年7月江蘇谷瑞斯機電工程有限公司提供的商用玉米為試驗物料，在揚州大學機械工程學院農(nóng)機實驗室完成試驗。本試驗主要儀器：數(shù)顯游標卡尺（精度0.01 mm）、自制摩擦系數(shù)測量儀、JA5003電子天平（量程500 g，精度0.01 g）、SX2-4-10型烘干箱（電壓220 V，功率4 KW）、無底不銹鋼圓筒和304不銹鋼平面。

1.2 物理參數(shù)標定

玉米種子按形狀可分為馬齒形、半球形和球錐形，如圖1所示。玉米種子屬于散狀物料，其基本物理參數(shù)主要包括玉米種子外形尺寸（長×寬×厚）、密度、含水率、泊松比和剪切模量等^［15］。材料成型后其物理特性幾乎不發(fā)生變化，固有參數(shù)（表1）可通過參考文獻獲取^［16-18］。通過數(shù)顯游標卡尺測量玉米種子平均尺寸（表2）。

1.3 接觸參數(shù)標定

離散元仿真試驗主要考慮參數(shù)：玉米種子之間和玉米種子與不銹鋼板之間的碰撞恢復(fù)系數(shù)、靜摩擦系數(shù)和滾動摩擦系數(shù)，利用實際物理試驗確定所需參數(shù)范圍，為離散元試驗參數(shù)確定提供參考。

1.3.1 靜摩擦系數(shù)測定 利用靜力學原理對摩擦系數(shù)進行測定，玉米種子摩擦系數(shù)測量原理如圖2所示。玉米種子受力分析滿足式（1）至式（3），當斜面傾角逐漸增大到使玉米種子有下滑趨勢時，此時的斜面傾角即為玉米種子的靜摩擦角，根據(jù)式（4）可得玉米種子靜摩擦系數(shù)^［19］。

f=F=mg·sinα；（1）

F₂=mg·cosα；（2）

f=μ₁.F₂；（3）

μ₁=f/F₂=tanα。（4）

式中：α為玉米種子靜摩擦角，°；μ₁為玉米種子靜摩擦系數(shù)；m為玉米種子重量，kg；g為重力加速度，9.8 N/s²。

本試驗利用摩擦系數(shù)測量儀（圖3）測量玉米種子間與接觸面間的靜摩擦系數(shù)。測量玉米種子與不銹鋼板的靜摩擦系數(shù)時，為了防止玉米種子在試驗中發(fā)生滾動影響試驗結(jié)果，將2粒玉米種子黏在一起作為一組^［20］，共設(shè)置10組（圖4）。試驗過程：將不銹鋼板置于斜面測量儀平面上，使其水平，保證與鋼材充分接觸，再將黏好的種子放置上面，緩慢上升斜板，當觀察到種子出現(xiàn)滑動現(xiàn)象時將斜板固定，利用測量儀自帶的角度尺測量角度，最后根據(jù)摩擦系數(shù)公式計算得到種子與鋼板之間的靜摩擦系數(shù)，每組試驗重復(fù)5次。測定玉米種子間靜摩擦系數(shù)時，將玉米種子粘在平面上，制成玉米板，然后重復(fù)上述試驗，最終得出玉米種子與不銹鋼板之間的靜摩擦系數(shù)在0.26～0.40之間，玉米種子之間的靜摩擦系數(shù)在0.23～0.34之間。

1.3.2 滾動摩擦系數(shù)測定 玉米種子在斜面上發(fā)生滾動時所產(chǎn)生的斜面傾角為玉米種子滾動摩擦角。滾動摩擦系數(shù)由摩擦系數(shù)測量儀確定，試驗分為5組，每組試驗重復(fù)10次，玉米種子之間的滾動摩擦系數(shù)在0.05～0.07之間，玉米種子與不銹鋼之間的滾動摩擦系數(shù)在0.03～0.09之間。

1.3.3 碰撞恢復(fù)系數(shù)測定 碰撞恢復(fù)系數(shù)是形容物體變形后恢復(fù)原本形狀的能力，定義恢復(fù)系數(shù)為碰撞結(jié)束時接觸點的瞬時法向分離速度與碰撞前法向接近速度比值^［21］。本研究根據(jù)參考文獻［22-24］確定，玉米種子間碰撞恢復(fù)系數(shù)范圍為0.18～0.28，玉米種子與不銹鋼間恢復(fù)系數(shù)為0.62～0.73。

1.3.4 堆積角測定 選用直徑60 mm、高180 mm的304不銹鋼圓筒。試驗時圓筒垂直放置于不銹鋼板上，玉米種子確保無損填充300 g置圓筒，隨后將圓筒以0.05 m/s的速度勻速上升，底板上生成錐形玉米堆，試驗重復(fù)5次。為了得到較準確的試驗數(shù)值，采用高清相機拍攝玉米堆積圖像，利用Matlab圖像處理技術(shù)，分別對圖像進行灰度處理和二值化處理，堆積角邊界線數(shù)據(jù)提取，最后利用Origin對邊界線數(shù)據(jù)進行擬合，擬合斜率即為所測角度正切值，最終獲得玉米種子堆積角平均值為24.95°（圖5）。玉米種子堆積角邊界擬合如圖6所示。

2 玉米種子離散元模型建立及仿真參數(shù)標定

2.1 玉米種子離散元模型建立

兼顧仿真時間和效率，采用9球、3球和4球在EDEM軟件中建立馬齒形、半球形和球錐形的離散元模型^［25］，如圖7所示。

EDEM仿真模擬堆積角試驗?zāi)Ｐ托璋凑諏嶋H堆積角測定裝置繪制。為保證仿真試驗與實際試驗條件相同，將圖7所示玉米種子模型按實際比例進行設(shè)置，馬齒形、半球形、球錐形比例分別為5∶2∶3。顆粒接觸模型選擇Hertz Mindlin（no slip）模型，在不銹鋼筒上方口徑生成顆粒工廠，用于玉米顆粒生產(chǎn)，顆粒共生成300 g，顆粒尺寸采用固定形式，仿真試驗從開始到結(jié)束總時長設(shè)置為5 s。瑞利時間步長取20%，網(wǎng)格尺寸大小取最小顆粒半徑的3倍^［26］。仿真開始時，顆粒填充滿圓筒需要3.5 s，3.5 s后設(shè)置圓筒以0.05 m/s速度勻速向上提升，經(jīng)過5 s仿真時長種子顆粒在底板上堆積完成，通過EDEM軟件量角器進行測量。

2.2 離散元仿真參數(shù)標定

2.2.1 確定顯著性參數(shù) 利用Design-Expert 13.0軟件進行Plackett-Burman設(shè)計，根據(jù)實際物理試驗范圍，以實際玉米種子堆積角為目標值，將6種因素分為高低2種水平，設(shè)置高水平為+1，低水平為-1。通過Plackett-Burman設(shè)計對仿真因數(shù)進行試驗篩選，試驗水平及因素見表3，最后對以下6種因素進行篩選，得到仿真參數(shù)最顯著因素，試驗結(jié)果如表4所示，參數(shù)顯著性分析如表5所示。由表5可知，玉米-玉米靜摩擦系數(shù)（X₁）、玉米-不銹鋼靜摩擦系數(shù)（X₄）、玉米-不銹鋼滾動摩擦系數(shù)（X₅）對玉米堆積角具有顯著影響（P＜0.05），其他參數(shù)遠大于0.05，表明其他因素對堆積角影響暫不顯著。

2.2.2 最陡爬坡試驗 由表5 所得結(jié)論，對篩選出的3個顯著參數(shù)X₁、X₄及X₅進行最陡爬坡試驗，對其他不顯著參數(shù)取平均值，通過離散元仿真分析實際堆積角與仿真堆積角相對誤差，確定仿真參數(shù)范圍，設(shè)計方案及結(jié)果如表6所示。

2.2.3 顯著參數(shù)響應(yīng)面優(yōu)化 由表6仿真試驗結(jié)果可知，2號試驗相對誤差最小，隨著參數(shù)X₁、X₄及X₅的增大，相對誤差呈先減小后增大的趨勢，分別取1號為低水平（-1），取2號為中心點（0），3號為高水平（+1），進行顯著參數(shù)Box-Behnken試驗，其余參數(shù)取值與最陡爬坡試驗設(shè)計相同，設(shè)計方案如表7所示。

采用Design-Expert軟件進行二次回歸建模，二次多項式模型方差分析，得出二次多項式方程

θ=186.615-340.178X₁-829.114X₄-891.020X₅+2 255.556X₁X₄+2 304.167X₁X₅-143.750X₄X₅-632.083X²₁+566.528X²₄+3 255.937 5X²₅。

模型方程分析結(jié)果如表8所示，其擬合模型P=0.005＜0.01，說明其自變量與因變量相關(guān)性極顯著；變異系數(shù)（CV）=5.68%＜10%，說明試驗依據(jù)具有可靠性；失擬項P=0.907 6＞0.05，說明方程擬合度較好，無失擬現(xiàn)象發(fā)生；決定系數(shù)R²=0.916 2接近于1；R²_adj=0.808 5，R²_pred=0.728 1且R²_adj-R²_pred＜0.2，說明該模型具有合理性；精密度10.364>4，說明模型精密度良好。

為進一步優(yōu)化模型數(shù)據(jù)，在保證模型顯著性良好及其決定系數(shù)高的前提下，去除不顯著因素，優(yōu)化回歸模型方差分析如表9所示，得到新方程：

θ=145.872-553.653X₁-507.715X₄+2 255.556X₁X₄。

由表9優(yōu)化模型分析可得，交互項X₁和X₄的P<0.01，說明對堆積角影響極顯著，X₁X₄的P<0.05，說明對堆積角影響顯著，X₁和X₄交互項影響效果如圖8所示。R²_adj-R²_pred交互項差值更小，模型合理性更高，精密度由10.364提高至14.310，說明可以用來預(yù)測玉米種子堆積角。

2.2.4 試驗證明 利用Design-Expert軟件將實際物理堆積角度帶入，求出最優(yōu)參數(shù)玉米-玉米靜摩擦系數(shù)X₁=0.273、玉米-不銹鋼靜摩擦系數(shù)X₄=0.28，其余不顯著參數(shù)取平均值（玉米-玉米滾動摩擦系數(shù)X₂=0.06、玉米-玉米恢復(fù)系數(shù)X₃=0.23、玉米-不銹鋼滾動摩擦系數(shù)X₅=0.06、玉米-不銹鋼恢復(fù)系數(shù)X₆=0.675）。為驗證仿真模擬試驗準確性，采用上述參數(shù)進行5次仿真，仿真堆積角分別為24.28°、25.40°、24.53°、25.40°、22.79°，得到平均仿真堆積角為24.48°，與實際物理堆積角24.95°相對誤差為1.88%，試驗對比如圖9所示，驗證了仿真模擬的準確性。

3 結(jié)論

通過實際物理試驗分別測量玉米種子馬齒形尺寸12.11 mm×8.76 mm×4.77 mm，球錐形尺寸11.54 mm×6.50 mm×5.47 mm，半球形尺寸8.30 mm×8.70 mm×8.20 mm，玉米種子之間靜摩擦系數(shù)、滾動摩擦系數(shù)分別為0.23～0.34、0.05～0.07，玉米種子與鋼板間靜摩擦系數(shù)、滾動摩擦系數(shù)分別為0.26～0.40、0.03～0.09，測得玉米種子實際物理堆積角為24.95°。

以實際試驗所得的參數(shù)為依據(jù)，進行Plackett-Burman試驗參數(shù)顯著性分析，挑選對堆積角存在顯著性影響的參數(shù)：玉米-玉米靜摩擦系數(shù)、玉米-不銹鋼靜摩擦系數(shù)和玉米-不銹鋼滾動摩擦系數(shù)，采用最陡滑坡試驗，確定顯著因素的取值范圍。

利用Box-Behenken試驗建立顯著性參數(shù)與堆積角之間的二次回歸模型并剔除不顯著參數(shù)對其進一步優(yōu)化，試驗表明玉米-玉米靜摩擦系數(shù)、玉米-不銹鋼靜摩擦系數(shù)兩者交互項對堆積角具有顯著影響。

利用響應(yīng)面試驗得到最優(yōu)參數(shù)：玉米-玉米靜摩擦系數(shù)0.273，玉米-不銹鋼靜摩擦系數(shù)0.28，玉米-玉米滾動摩擦系數(shù)0.06、玉米-玉米恢復(fù)系數(shù)0.23，玉米-不銹鋼滾動摩擦系數(shù)0.06，玉米-不銹鋼恢復(fù)系數(shù)0.675，利用以上參數(shù)進eqOp/G7xqVZfkEI5CmHSAKBltXKqiQV17erOUzoWItY=行仿真試驗，得到平均仿真堆積角為24.48°，與實際物理堆積角24.95°，相對誤差為1.88%，驗證了仿真模擬的準確性。

參考文獻：

［1］宋偉剛，陳洪亮，李勤良，等. 散狀物料轉(zhuǎn)載沖擊載荷的DEM仿真［J］. 東北大學學報（自然科學版），2011，32（11）：1631-1634.

［2］張國成. 玉米種植現(xiàn)狀與新技術(shù)應(yīng)用的效率探究［J］. 農(nóng)家參謀，2021（7）：35-36.

［3］曾智偉，馬 旭，曹秀龍，等. 離散元法在農(nóng)業(yè)工程研究中的應(yīng)用現(xiàn)狀和展望［J］. 農(nóng)業(yè)機械學報，2021，52（4）：1-20.

［4］馬令然，曹樹坤，鐘文鎮(zhèn)，等. 基于EDEM的清選裝置作業(yè)參數(shù)仿真研究［J］. 農(nóng)業(yè)技術(shù)與裝備，2017（7）：80-83.

［5］蔣恩臣，孫占峰，潘志洋，等. 基于CFD-DEM的收獲機分離室內(nèi)谷物運動模擬與試驗［J］. 農(nóng)業(yè)機械學報，2014，45（4）：117-122.

［6］劉大為，段佳鵬，陳星宇，等. 基于EDEM的稻種離散元摩擦參數(shù)標定［J］. 沈陽農(nóng)業(yè)大學學報，2023，54（2）：189-195.

［7］鄭效帥，何曉寧，尚書旗，等. 油莎豆種子籽粒離散元仿真參數(shù)標定與試驗［J］. 農(nóng)機化研究，2024，46（2）：172-178.

［8］周 濱，周 嶺，萬 暢，等. 基于休止角試驗的顆粒鉀肥離散元模型參數(shù)標定［J］. 農(nóng)機化研究，2024，46（1）：46-53.

［9］張榮芳，周紀磊，劉 虎，等. 玉米顆粒粘結(jié)模型離散元仿真參數(shù)標定方法研究［J］. 農(nóng)業(yè)機械學報，2022，53（增刊1）：69-77.

［10］李飛翔，王 鵬，王云飛，等. 基于堆積試驗的玉米包衣種子離散元參數(shù)標定［J］. 中國農(nóng)業(yè)科技導報，2022，24（7）：97-107.

［11］王美美，王萬章，楊立權(quán)，等. 基于響應(yīng)面法的玉米籽粒離散元參數(shù)標定［J］. 華南農(nóng)業(yè)大學學報，2018，39（3）：111-117.

［12］陳國強，閆俊霞，張連仁. 玉米顆粒堆積特性的離散元法模擬研究［J］. 輕工機械，2020，38（6）：11-14.

［13］戴繼光. 基于PFC^3D的玉米籽粒離散元參數(shù)研究［J］. 農(nóng)業(yè)與技術(shù)，2021，41（16）：22-26.

［14］王美美，王萬章，楊立權(quán)，等. 基于EDEM的玉米子粒建模方法的研究［J］. 河南農(nóng)業(yè)大學學報，2018，52（1）：80-84，103.

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

［16］金秀南，張軍昌，薛晶峰，等. 玉米種子與橡膠帶離散元接觸參數(shù)標定［J］. 農(nóng)機化研究，2022，44（7）：39-43.

［17］張文龍，武文斌，呂少杰，等. 制粉中的小麥力學參數(shù)和接觸參數(shù)試驗標定［J］. 河南工業(yè)大學學報（自然科學版），2023，44（2）：105-112.

［18］吳文江，郭 斌，高占鳳，等. 基于DEM的沙土顆粒建模及參數(shù)標定研究［J］. 中國農(nóng)機化學報，2019，40（8）：182-187.

［19］徐效偉，魏 海，顏建春，等. 花生莢果離散元仿真參數(shù)標定［J］. 中國農(nóng)機化學報，2022，43（11）：81-89.

［20］丁文波，朱繼平，陳 偉，等. 基于EDEM的青稞接觸參數(shù)仿真標定［J］. 中國農(nóng)機化學報，2021，42（9）：114-121.

［21］彭才望，周 婷，宋世圣，等. 基于Hertz接觸理論的黑水虻幼蟲碰撞恢復(fù)系數(shù)測定［J］. 農(nóng)業(yè)機械學報，2021，52（11）：125-134.

［22］周文秀. 玉米籽粒的物理力學特性研究［D］. 哈爾濱：東北農(nóng)業(yè)大學，2016.

［23］Dai F，Zhao W Y，Han Z S，et al. Experiment on Poisson’s ratio determination about corn kernel［J］. Advanced Materials Research，2013，781-784：799-802.

［24］Co?kunM B，Yal??n I·，?zarslan C. Physical properties of sweet corn seed （Zea mays saccharata Sturt.）［J］. Journal of Food Engineering，2006，74（4）：523-528.

［25］Shi L R，Zhao W Y，Yang X P. Effects of typical corn kernel shapes on the forming of repose angle by DEM simulation［J］. International Journal of Agricultural and Biological Engineering，2022，15（2）：248-255.

［26］胡國明. 顆粒系統(tǒng)的離散元素法分析仿真：離散元法的工業(yè)應(yīng)用與EDEM軟件簡介［M］. 武漢：武漢理工大學出版社，2010.