摘要：土壤沙化問題是目前我國制約農業(yè)發(fā)展的因素之一。設計一款農田沙化土壤改良機，利用SolidWorks建立其整機模型，針對其關鍵部件攪拌筒建立離散元模型，并使用EDEM軟件進行離散元分析，以攪拌軸旋轉速度、旋轉方向和攪拌時長為試驗因素，以土壤和改良液混合均勻程度為優(yōu)化指標，分析改良裝置攪拌筒工作情況。通過正交仿真試驗結果得出：最優(yōu)參數為攪拌架旋轉速度40 r/min、攪拌時間15 s、兩攪拌架轉動方向相反時，其攪拌均勻度為92.61%，變異系數為0.074 9。

關鍵詞：沙化土壤；改良機；攪拌效果；離散元分析；均勻度

中圖分類號：S156" " " 文獻標識碼：A" " " 文章編號：2095?5553 （2024） 09?0285?06

Design of farmland sandy soil improvement machine and numerical simulation

of its mixing effect

Li Bin， Pang Mengyu

（College of Mechanical and Electrical Engineering， Northeast Forestry University， Harbin， 150040， China）

Abstract： The problem of soil desertification is one of the factors restricting the development of agriculture in China. A farmland desertification soil improvement machine was designed. The whole machine model was established by SolidWorks， and the discrete element model was established for the key component mixing cylinder. The EDEM software was used for discrete element analysis. The rotation speed， rotation direction and stirring time of the mixing shaft were used as experimental factors， and the mixing uniformity of soil and improved liquid was used as the optimization index to analyze the working condition of the mixing cylinder of the improved device. Through the orthogonal simulation test results， it is concluded that when the stirring frame rotation speed is 40 r/min， the stirring time is 15 s， and the rotation direction of the two stirring frames is opposite， the stirring uniformity is 92.61%， and the coefficient of variation is 0.074 9， which is the optimal parameter.

Keywords： desertified soil; improved device; stirring effect; discrete element analysis; uniformity

0 引言

土壤沙化泛指良好的土壤或可利用的土地變成含沙很多的土壤甚至變成沙漠的過程。傳統(tǒng)的治理模式工程周期長、工程花費高，且在實施時易受多方面影響。

目前，土壤改良主要通過改善土壤結構，促進土壤蓄水保墑。閆少鋒等[1]通過豎井抽排水試驗，探討天然降水條件下地下水位動態(tài)及土壤脫鹽規(guī)律；劉長江等[2]通過田間示范試驗，研究了不同深松深度對蘇打鹽堿化田土壤理化性狀及作物產量的影響；魯向暉等[3]研究發(fā)現，機械深耕深松能有效解決土壤的板結。截至目前，國內土壤改良方法多種多樣，但小型化機械化沙化土壤改良技術鮮見報道。因此，設計一種智能化農田沙化土壤改良裝置是很有必要的，改良裝置不僅可重新賦予沙化土壤蓄水能力和養(yǎng)分保持能力，且工程周期短、工程花費較低、工程規(guī)模小[4]。

本文以沙化土壤改良機為研究目標，通過SolidWorks2021建立整機三維模型，并將攪拌筒機構細化后將其導入EDEM 2022軟件建立離散元仿真模型，對攪拌過程進行分析。以攪拌軸的旋轉速度、旋轉方向以及攪拌時長為試驗因素，以沙粒和土壤改良液的混合均勻程度為目標指標，設計正交試驗，分析各因素對試驗結果的影響，驗證運動方程的正確性并尋求最佳參數組合。

1 沙土改良機的整機結構與工作原理

1.1 整機結構

該裝置三維建模如圖1所示，結構如圖2所示，主要由履帶輪底盤、PLC電控裝置、攪拌裝置、注液裝置、輸料裝置等部分組成。

1.2 工作原理

該裝置主要工作流程：由螺旋上沙裝置將目標沙子輸送至攪拌筒內部，通過土壤溫濕度傳感器和土壤pH傳感器識別沙化土壤性質后將數據上傳至PLC控制器，實現準確檢測目標土壤的成分及含量，控制器識別處理數據并根據用戶需求發(fā)出信號，可控閥門打開，向攪拌筒注入改良液，攪拌電機啟動，帶動攪拌架將纖維提取液和沙土均勻攪拌，攪拌完成后，電機停機，閥門關閉，達到預定標準后，出料口擋板開啟，將改良后的土壤由出料口排出，均勻覆蓋到地面之上。

2 攪拌筒結構及主要參數

本文重點研究對象為土壤攪拌筒，主要由外筒、攪拌架、套軸、動力裝置和進出沙口構成，為了便于裝置的仿真和計算，對該結構進行簡化，如圖3所示。

攪拌筒內的攪拌裝置由一根實心階梯攪拌軸及其葉片、一根空心攪拌軸及其葉片和齒輪箱、攪拌筒組成，采用反向攪拌搭配線形攪拌葉片的工作模式，將實心攪拌軸部分軸體插入空心攪拌軸。該攪拌結構的混合原理：在攪拌軸上分布按一定角度構成連續(xù)大螺旋葉片，主軸旋轉帶動大螺旋葉片連續(xù)循環(huán)，使得外側土壤向中間流動，內側土壤向兩端運動，葉片和攪拌臂的合理布置可使被混合物料形成連續(xù)的對流、擴散、揉搓等混合運動，產生三維攪拌效果。旋轉攪拌架的結構如圖4所示。

設計攪拌筒為圓柱型，從而有效防止改良后的土壤卡縫、板結，造成工作效率降低，攪拌筒的主要參數如表1所示。

3 離散元模型建立和參數選擇

3.1 容器模型建立

將SolidWorks建立的攪拌筒簡化模型導入EDEM，進行參數設置，并對攪拌螺帶的參數進行計算和優(yōu)化。在確定螺旋帶螺旋升角前做出如下假設：假定飼料在螺旋帶上的單位阻力系數為一常數，因此可以假定飼料跟螺旋帶之間的壓力分布在螺旋攪拌帶的平均半徑上。同時忽略螺旋升角差異，將螺旋帶展開后可得到一個連續(xù)斜面。這樣可將復雜的螺旋帶摩擦問題轉換為斜面摩擦問題來考慮計算[5]。

如圖5所示，對某一在螺帶斜面上的假定沙粒質點進行受力分析，其中P為該質點所受的總驅動力，N；F為總阻力，N；N為螺帶面給予該沙粒的總支撐力，N；α為螺帶的螺旋升角，（°）；δ為沙粒質點與鋼制螺帶之間的摩擦角，（°）；μ為沙粒與鋼之間的摩擦系數。查摩擦系數表可得，一般沙粒與鋼之間的摩擦系數為0.45～0.5。

3.2 容器模型建立

沙化土壤受到風化和水蝕的影響，其沙粒結構復雜，但由于長期的磨損，大多接近球形和橢球形，直徑大約在1 mm，為了減少計算量和提高計算精度，對沙粒模型全部設置為直徑1.1 mm的球體；由于改良液是液體，其顆粒模型建立較為復雜，為了方便計算和仿真，在改良液中加入甲基紫示蹤顆粒，其顆粒定義為球形，直徑大約在0.85 mm，沙粒模型和示蹤顆粒模型如圖8、圖9所示。

3.3 模型力學參數確定

在EDEM仿真過程中，沙粒、示蹤顆粒以及攪拌筒容器之間一直發(fā)生著碰撞、粘結等作用，故應當選擇合適的仿真參數。在本試驗中，由于改良液是液體，附著在沙粒表面之上，使之具有了JKR表面能，因此應當選用Hertz-Mindlin with JKR模型，該模型通過設置顆粒表面的表面能在模擬含水量，進而得到了顆粒之間的粘結狀態(tài)。

接觸參數指顆粒與顆粒之間、顆粒與幾何體之間的恢復系數、靜摩擦系數和滾動摩擦系數[6]。因為摩擦系數和材料的形狀、濕度和材質關聯性較大，而仿真中模型全部采用球型顆粒，實際中物料的形狀是不規(guī)則的，因此接觸系數并不能采用真實狀態(tài)下的值，需要通過虛擬參數標定試驗來獲得[7]。沙粒、示蹤顆粒、攪拌筒的材料參數和力學參數如表2、表3所示。

需要說明的是，由于沙粒等模型的復雜性[8]，表2、表3內的值均為平均值，對提高仿真的效率和可及性有重要的意義[9]。

3.4 顆粒工廠的建立

在攪拌筒的兩端分別建立示蹤物質和沙粒的顆粒工廠，為滿足攪拌筒的充盈系數，分別計算得示蹤物質顆粒和沙子顆粒的生成數量以及生成速度，為了簡化仿真模型、減小誤差，設置兩種顆粒在3 s內加載完成，然后關閉顆粒工廠，兩種顆粒工廠的設置參數如表4所示。

4 正交試驗設計

4.1 試驗因素選取

選取對混合均勻度影響較為顯著的三個因素：攪拌架轉速、攪拌時間、攪拌架的轉動方向，設計三因素三水平的正交試驗。攪拌架的轉速取值是最重要的參數，若轉速太慢會使得物料在攪拌筒空間里的運動不夠充分，攪拌不夠均勻；如果其轉速太快而超過了臨界值，物料在離心力的作用下會附著在攪拌筒內壁上，從而無法達到充分攪拌的目的。本樣機的工作方式為階段式工作，沙化土壤的改良過程由上沙、加改良劑、攪拌及出料四個主要階段構成[10]。攪拌時間指攪拌階段的持續(xù)時間，攪拌時間是影響攪拌質量的因素之一，攪拌時間應在保證攪拌質量的前提下盡量的短。本設備采用兩個獨立的螺帶式攪拌架[11]，安裝方向共有三種情況，即均為正裝、均為反裝和一正一反，對沙化土壤的改良均勻度也有著不可忽視的影響，應當對其進行試驗檢驗[12]。

另外，選取本次試驗的評價指標為攪拌后樣本的混合均勻度M?；旌暇鶆蚨仁怯脕砻枋龌旌衔飫蛸|性的量化指標，即對樣本總體進行隨機抽樣，然后測定每個子樣本的特性指標，運用數理統(tǒng)計方法，對一組數據進行方差分析，計算出特性指標X與標示量理論值p的偏離程度。

具體取樣標準：在混合機上中下部間隔抽取樣本，或等效采用其他方法取樣；取樣數量：每批試驗物料的樣本數量不得少于10個[13]。

樣本標準差

4.2 確定各試驗因素水平

由于沙粒的運動趨勢和受力方向相同，根據牛頓第二定律，可得

將上文所得數據，取α=15°代入式（11），可得n≤67.95 r/min。約定逆時針為攪拌架旋轉正向，對試驗因素確定如表5所示，分別確立3個水平。

為了提高攪拌均勻度，試驗的指標為單指標均勻度。本例是一個三水平的試驗，共有三個因素，在對正交試驗做方差分析時，必須估計隨機誤差，而隨機誤差是通過正交表上空白列得到的。由于空白列中沒有因素作用，因此正好反映隨機因素所引起的誤差[6]。為了減小誤差，最好留有一個空列，所以選用正交表L9（34）來進行試驗。

5 仿真試驗后處理結果分析

在按照正交試驗方案進行試驗后，需要使用EDEM中的后處理（Analyst）功能對試驗所得數據進行統(tǒng)計和分析。

5.1 試驗取樣網格劃分

后處理中Create Grid Bin Groups功能可以對試驗結束的工況進行網格劃分，本變異系數試驗的取樣要求網格數不少于10個，在仿真試驗中，對攪拌筒進行4×4×4共計64個網格的劃分，網格大小確定為0.1 m×0.1 m×0.2 m，網格編號以（x，y，z）命名，最左下角網格為（1，1，1），分別以坐標軸正向遞增。

5.2 試驗結果統(tǒng)計

以試驗號1為例，即攪拌架轉速20 r/min，攪拌時間60 s，兩攪拌架轉向均為正向的試驗，對試驗數據進行分析計算，仿真時間結束時，利用EDEM的Analyst功能對攪拌筒內的顆粒情況進行查看[14]，如圖10所示。

從圖10中可以看出，示蹤物質顆粒與沙粒相對均勻的混合在一起，下面對其試驗數據進行計算，從而分析其攪拌均勻度，去除其中顆粒較少的網格[15]，取攪拌筒中部16個網格樣本，每個網格內的顆粒數量如表6所示。

沙粒顆粒和初始示蹤物質比值為5.00，對試驗號1所取樣本數據進行計算，結果如表7所示。

5.3 正交試驗結果

將9個試驗結果進行正交試驗數據分析，將數據填入表8，可以得出各試驗因素的主次和優(yōu)水平，并得出本試驗的最優(yōu)解。

5.4 正交試驗最優(yōu)解驗證

由于最優(yōu)解不在前述正交試驗的設計中，因此需要對最優(yōu)解進行仿真試驗[17]，即取攪拌架轉速為40 r/min，攪拌時間15 s，攪拌架轉動方向為正+反，利用EDEM對該最優(yōu)解進行驗證[18]，其結果如表9所示。通過對最優(yōu)解的仿真試驗驗證，其均勻度為92.61%，為最優(yōu)結果，符合正交試驗的試驗結果，驗證了選優(yōu)方法和計算的正確性。

6 結論

1） 仿真結果論證土壤改良攪拌倉以及螺旋攪拌架的設計符合土壤改良的原理，并驗證正交試驗設計思路的正確和可行性。

2） 通過EDEM仿真試驗驗證本沙化土壤改良機攪拌筒的設計符合預期要求，對各試驗參數進行選優(yōu)，最終確定攪拌架轉速為40 r/min，攪拌時間15 s，攪拌架轉動方向正+反時，變異系數為0.074 9，混合均勻度為92.61%，為最優(yōu)結果，符合沙化土壤改良的要求。

參 考 文 獻

[ 1 ] 閆少鋒， 吳玉柏， 俞雙恩， 等. 江蘇沿海地區(qū)豎井排鹽試驗研究[J]. 節(jié)水灌溉， 2014（8）： 42-44.

Yan Shaofeng， Wu Yubai， Yu Shuang'en， et al. Experimental of desalination by wells in the coastal area of Jiangsu province [J]. Water?saving Irrigation， 2014（8）： 42-44.

[ 2 ] 劉長江， 李取生， 李秀軍. 深松對蘇打鹽堿化旱田改良與利用的影響[J]. 土壤， 2007， 39（2）： 306-309.

Liu Changjiang， Li Qusheng， Li Xiujun. Effect of deep tillage on amelioration and utilization of soda?alkaline upland fields [J]. Soil， 2007， 39（2）： 306-309.

[ 3 ] 魯向暉， 高鵬， 王飛， 等. 寧夏南部山區(qū)秸稈覆蓋對春玉米水分利用及產量的影響[J]. 土壤通報， 2008， 39（6） ： 1248-1251.

Lu Xianghui， Gao Peng， Wang Fei， et al. Effect of straw mulch on yield and water utilization of spring maize in the southern mountainous area of Ningxia[J]. Soil Bulletin， 2008， 39（6）： 1248-1251.

[ 4 ] 褚雅紅. 小流域治理中水土保持措施及效益分析[J]. 內蒙古水利， 2020（2）： 56-57.

[ 5 ] 劉倬昀. 基于流固耦合分析下分區(qū)土石壩優(yōu)化研究[D]. 楊凌： 西北農林科技大學， 2015.

Liu Zhuoyun. Optimization study for districting dam based on fluid?structure coupling analysis [D]. Yangling： Northwest A amp; F university， 2015.

[ 6 ] 滕海英， 祝國強， 黃平， 等. 正交試驗設計實例分析[J]. 藥學服務與研究， 2008（1）： 75-76.

[ 7 ] 朱澤華， 吳笑宇， 張恒運. 鋰離子電池模組熱擴散結構的優(yōu)化[J]. 機械設計與制造， 2023（1）： 257-263.

Zhu Zehua， Wu Xiaoyu， Zhang Hengyun. Thermal structural optimization of lithium?ion battery for improving thermal performance [J]. Mechanical Design and Manufacture， 2023 （1）： 257-263.

[ 8 ] 易志堅. 沙漠“土壤化”生態(tài)恢復理論與實踐[J]. 重慶交通大學學報（自然科學版）， 2016， 35（S1）： 27-32.

Yi Zhijian. Theory and practice of desert “soilization” for ecological restoration [J]. Journal of Chongqing Jiaotong University （Natural Science Edition）， 2016， 35（S1）： 27-32.

[ 9 ] 邵輝. 渭河流域水土流失變化對梯田措施響應的模擬研究[D]. 楊凌： 西北農林科技大學， 2013.

Shao Hui. Simulation of soil and water loss variation toward terrace practice in the Weihe river basin [D]. Yangling： Northwest A amp; F University， 2013.

[10] 蔡明嫻. 防沙治沙造林技術的實踐探析[J]. 農家參謀， 2021（2）： 174-175.

[11] 盧微坦. 干混砂漿攪拌機攪拌方式的研究[J]. 建設機械技術與管理， 2022， 35（4）： 112-116.

[12] 金偉亮. 螺旋帶式動物飼料攪拌機的結構設計與研究[D]. 合肥： 安徽農業(yè)大學， 2015.

[13] 張翼夫， 李問盈， 胡紅， 等. 鹽堿地改良研究現狀及展望[J]. 江蘇農業(yè)科學， 2017， 45（18）： 7-10.

[14] 童世合， 邵明璽， 曹猛， 等. 基于DEM的玉米秸稈離散元模型參數標定[J]. 中國農機化學報， 2023， 44（2）： 69-75.

Tong Shihe， Shao Mingxi， Cao Meng， et al. Parameter calibration of corn straw discrete element model based on DEM [J]. Journal of Chinese Agricultural Mechanization， 2023， 44（2）： 69-75.

[15] 胡天亮， 胡靖明， 畢陽， 等. 基于EDEM數值模擬的沙棘滾筒篩設計[J]. 中國農機化學報， 2021， 42（9）： 122-129.

Hu Tianliang， Hu Jingming， Bi Yang， et al. Design of sea?buckthorn roller screen based on EDEM numerical [J]. Journal of Chinese Agricultural Mechanization， 2021， 42（9）： 122-129.

[16] 張凡， 李玲玲， 馬學東. 基于EDEM?FLUENT耦合礦物顆粒分離效果的模擬研究[J]. 礦產綜合利用， 2022（6）： 159-166.

[17] 李景恩， 聶少平， 邵燈寅， 等. 江香薷揮發(fā)油提取率與提取時間的動力學關系[J]. 食品與生物技術學報， 2014， 33（3）： 308-314.

Li Jing'en， Nie Shaoping， Shao Dengyin， et al. Dynamic relationship between extraction yield and extraction time of essential oil from mosla chinensis maxim cv. Jiangxiangru [J]. Journal of Food and Biotechnology， 2014， 33（3）： 308-314.

[18] 牛瀟瀟. 油莎豆混流式烘干塔設計與試驗研究[D]. 鄭州： 河南農業(yè)大學， 2022.