摘要：

為解決大面積的淺水藕田地區(qū)的蓮藕采收機(jī)執(zhí)行部件高壓噴頭的個(gè)數(shù)較多，容易出現(xiàn)各噴頭間流速差過大的現(xiàn)象，導(dǎo)致蓮藕采收率的下降和蓮藕品質(zhì)的參差不齊的問題，對(duì)蓮藕采收機(jī)射流系統(tǒng)執(zhí)行部件進(jìn)行優(yōu)化。噴頭作為執(zhí)行機(jī)構(gòu)的重要組件，其射流的穩(wěn)定性對(duì)蓮藕采收的效率以及沖刷效果有著重大的影響。為揭示高壓噴頭的連接工藝和噴頭形狀、口徑對(duì)射流穩(wěn)定性的影響，基于EDEM-Fluent軟件對(duì)內(nèi)部流場進(jìn)行環(huán)境建模以及對(duì)執(zhí)行部件的噴頭進(jìn)行耦合仿真模擬。試驗(yàn)發(fā)現(xiàn)：改變噴頭的形狀增加倒角，調(diào)整各噴頭出水口直徑到恰當(dāng)?shù)牟钪祵?duì)噴頭出水穩(wěn)定性有明顯的改觀，流速波動(dòng)標(biāo)準(zhǔn)差由0.164m/s改變?yōu)?.020m/s；平均挖掘深度由33.50cm改變?yōu)?6.07cm，挖掘深度標(biāo)準(zhǔn)差由1.71cm改為0.69cm。

關(guān)鍵詞：蓮藕采收機(jī)；射流系統(tǒng)；EDEM-Fluent耦合仿真；出水穩(wěn)定性；流體力學(xué)

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

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

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

-06

收稿日期：2022年8月19日" 修回日期：2022年10月18日

基金項(xiàng)目：國家自然科學(xué)基金（51875531）

第一作者：陳強(qiáng)，男，1997年生，浙江義烏人，碩士研究生；研究方向?yàn)檗r(nóng)業(yè)機(jī)械自動(dòng)化。E-mail： 602409669@qq.com

通訊作者：錢孟波，男，1981年生，寧波人，博士，教授，博導(dǎo)；研究方向?yàn)殚g隙機(jī)構(gòu)非線性動(dòng)力學(xué)。E-mail： 46704548@qq.com

Design and optimization of jet coupling machine based on EDEM-Fluent coupling simulation

Chen Qiang， Qian Mengbo， Yu Lang， Sun Fuxing， Pan Jiaxuan

（College of Optics， Mechanical and Electrical Engineering， Zhejiang Agriculture and Forestry University，

Hangzhou， 310000， China）

Abstract：

In order to solve the problem that there are a large number of high-pressure nozzles in the executive parts of the lotus root harvester in a large area of shallow water field， and the flow velocity difference between the nozzles is too large， which leads to the decline of the recovery efficiency of lotus root and the uneven quality of lotus root， this paper optimizes the executive parts of the lotus root harvester jet system. As an important component of the actuator， the stability of the jet has a great influence on the harvesting efficiency and scouring effect of lotus root. In order to reveal the influence of the connection process of the high pressure nozzle and the shape and diameter of the nozzle on the jet stability， based on the Edem-Fluent software， this paper carries out the environment modeling of the internal flow field and the coupling simulation of the nozzle of the executing part. It is found that by changing the shape of the nozzle to increase the chamfer and adjusting the diameter of each nozzle to the appropriate difference value， the stability of the nozzle is obviously improved. The standard deviation of the flow velocity fluctuation is changed from 0.164 m/s to 0.020 m/s. The average excavation depth is changed from 33.50 cm to 36.07cm， and the standard deviation of excavation depth is changed from 1.71 to 0.69.

Keywords：

lotus root harvester; firing system; EDEM-Fluent" coupling simulation; water stability; fluid mechanics

0 引言

蓮藕主要分布于我國南方以及安徽、福建等地區(qū)［1］。蓮藕具有很好的藥用價(jià)值，有除熱清胃、清血化瘀的功效［2］。根據(jù)生長位置的不同可以分為淺水藕和深水藕，國內(nèi)蓮藕種植面積穩(wěn)步增大［3］。

近些年隨著農(nóng)村剩余勞動(dòng)力的快速轉(zhuǎn)移，農(nóng)業(yè)機(jī)械化得到快速發(fā)展［4， 5］，但目前蓮藕采收的方式還是以人工采挖為主［6］，這種方式不僅采收的蓮藕品質(zhì)參差不齊，而且工作效率低、人力成本高，而隨著教育水平的提高，務(wù)農(nóng)人員越來越少，農(nóng)戶因不能及時(shí)采收而腐爛在藕田里造成的經(jīng)濟(jì)損失備受困擾［7］。對(duì)農(nóng)戶擴(kuò)大種植規(guī)模造成極大的阻礙。傳統(tǒng)的人工挖藕方式已經(jīng)不能滿足農(nóng)民的需求。如何才能又好又快的將蓮藕挖出來是當(dāng)務(wù)之急需要解決的問題。實(shí)現(xiàn)農(nóng)業(yè)機(jī)械化，減輕勞動(dòng)強(qiáng)度，以此提高生產(chǎn)效率［8］。

蓮藕在東南亞國家是一種常見的蔬菜，尤其是日本對(duì)蓮藕采收機(jī)已有一定的研究基礎(chǔ)［9］，日本蓮藕的種植面積達(dá)到了600hm2，收獲季節(jié)作業(yè)非常繁重，勞動(dòng)力卻不足，因而迫切要求蓮藕采收作業(yè)的機(jī)械化，早在20世紀(jì)80年代，日本便研制出了噴流式挖藕機(jī)［5］；我國對(duì)蓮藕采收機(jī)也有一定的研究，主要的蓮藕采收機(jī)分為挖掘式和射流式兩種［10， 11］。國內(nèi)高校陸續(xù)研制出多款射流式船體挖藕機(jī)，在2017年對(duì)驅(qū)動(dòng)進(jìn)行優(yōu)化，研制出液力驅(qū)動(dòng)挖藕機(jī)。2018年華中農(nóng)業(yè)大學(xué)工學(xué)院設(shè)計(jì)的自旋射流式挖藕機(jī)，發(fā)明射流自旋式管路結(jié)構(gòu)，其主要特點(diǎn)是簡化挖藕機(jī)的整體結(jié)構(gòu)，減少作業(yè)時(shí)所需要的能源動(dòng)力，并在提高作業(yè)寬幅和效率的基礎(chǔ)上降低制造成本［12］。

為解決大面積淺水藕田地區(qū)蓮藕采收機(jī)構(gòu)執(zhí)行部件高壓噴頭個(gè)數(shù)較多，出現(xiàn)噴頭建流速差過大的現(xiàn)象導(dǎo)致蓮藕采收率下降采收品質(zhì)層次不齊的問題，本文基于EDEM-Fluent軟件對(duì)內(nèi)部流場進(jìn)行環(huán)境建立，并對(duì)執(zhí)行部件進(jìn)行耦合仿真。

1 蓮藕采收機(jī)結(jié)構(gòu)及工作原理

1.1 蓮藕采摘機(jī)整體結(jié)構(gòu)設(shè)計(jì)

目前國內(nèi)挖藕多采用人工挖藕或者半自動(dòng)挖藕，其不僅工作繁重效率低，而且蓮藕破損率也比較高。本文設(shè)計(jì)一種射流式蓮藕采收一體機(jī)，其整體結(jié)構(gòu)如圖1所示。

射流式蓮藕采收一體機(jī)主要應(yīng)用于中等規(guī)模的淺水藕田環(huán)境。結(jié)合藕的物料特性分析，查閱相關(guān)資料，詢問導(dǎo)師以及專家意見，確定了車體式蓮藕采收機(jī)的主要性能參數(shù)如表1所示。

1.2 蓮藕采收機(jī)采收的工作原理

蓮藕采收機(jī)由機(jī)架、除莖裝置、高壓水槍裝置、收集裝置等組成。首先，由工作人員駕駛車體按工作路徑駕駛小車，啟動(dòng)采收按鈕，采收電機(jī)啟動(dòng)，拔禾器將蓮莖、葉卷入，由滾刀切碎從左側(cè)通道排出；水泵啟動(dòng)，高壓水槍噴射水流沖碎泥土，藕從淤泥中浮出水面，如圖2所示；藕通過鏟斗倒入傳送帶中，傳送帶上的水槍與毛刷清洗藕上淤泥，最后通過傳送帶進(jìn)入到收集框中，挖藕過程到此結(jié)束。

2 基于EDEM-Fluent的射流系統(tǒng)出水均勻性仿真分析

噴頭作為挖藕機(jī)執(zhí)行機(jī)構(gòu)的重要組件，其噴頭噴射水流的穩(wěn)定性對(duì)挖藕的效率以及沖刷效果有著重大的影響。郭洋民［6］對(duì)不同的進(jìn)水方式對(duì)水流的流動(dòng)狀態(tài)進(jìn)行了分析，肖科星對(duì)噴頭安裝角度、水流量、流速對(duì)挖藕深度和范圍的影響進(jìn)行了研究分析，但都未對(duì)噴頭的孔徑以及形狀對(duì)各噴頭出水穩(wěn)定性進(jìn)行分析。本文以降低噴頭水流能量損耗為目標(biāo)，仿真分析了不同形狀內(nèi)徑下各噴頭流速平穩(wěn)程度和挖掘深度的情況，為射流車體式挖藕機(jī)的設(shè)計(jì)提供參考。

2.1 有限元模型建立

2.1.1 土壤模型建立

Hertz-Mindlin with JKR模型是一種考慮了接觸區(qū)域內(nèi)范德華力影響的具有內(nèi)聚力的接觸模型［13］，主要運(yùn)用于黏性材料的分析，藕田環(huán)境下的淤泥具有一定的黏性所以采用該模型進(jìn)行建模分析。

在EDEM中建立土壤顆粒模型時(shí)，根據(jù)文獻(xiàn)調(diào)研，建立藕田泥土的堆積角平均值為42°，體積密度為1714kg/m3，土壤彈性模量為1143.2kPa，土壤剪切模量為408.3kPa，泊松比為0.4；選擇JKR模型，設(shè)置其JKR系數(shù)為13.05，恢復(fù)系數(shù)為0.5，靜摩擦因素為1.06，動(dòng)摩擦因素為0.15［14］。

以上述建立的泥土顆粒為基本顆粒，采用顆粒工廠自動(dòng)填充建立一個(gè)大?。ㄩL×寬×高）為0.65m×0.3m×0.40m的基床儲(chǔ)存土壤，建立好的土壤顆粒床如圖3所示。

2.1.2 高壓噴頭模型建立

噴水器主要由橫主管、進(jìn)水口、出水口三部分組成，橫主管外徑為100mm，長度為1220mm，厚度t=2mm，有效作業(yè)幅寬為1000mm，設(shè)置了5個(gè)噴頭以及1個(gè)進(jìn)水口，噴頭均勻的安裝在分水橫主管上。

為了便于網(wǎng)格劃分，提高運(yùn)算速度，對(duì)幾何模型進(jìn)行簡化，刪除了出水口除了管道外的輪廓模型，忽略了橫管與出水口和入水口的連接形狀和方式對(duì)水壓損耗的影響。僅考慮出水口噴頭的形狀以及直徑對(duì)水壓損耗的影響，其劃分好網(wǎng)格后的簡化模型如圖4所示。

2.2 材料屬性和邊界條件設(shè)定

定義材料屬性，打開Fluent軟件，導(dǎo)入已經(jīng)在mesh里劃分好網(wǎng)格的模型，假設(shè)噴水裝置里的water-liquid處于理想狀態(tài)，設(shè)定基于Pressure Based（壓力求解器），求解水的Steady（定常流動(dòng)），設(shè)定重力加速度大小為9.8m/s，方向?yàn)閅軸負(fù)方向。選擇普適性的k-epsilon（2eqn）Realizable 湍流模型，該模型適用于普遍的湍流運(yùn)動(dòng)，在滿足了雷諾應(yīng)力的約束條件上更符湍流的真實(shí)運(yùn)動(dòng)，可以更精確的計(jì)算擴(kuò)散速度，計(jì)算結(jié)果更符合實(shí)際狀態(tài)，得到的結(jié)果更加精確［15］。

定義邊界條件，設(shè)定進(jìn)水口為Velocity inlet（速度入口），設(shè)定湍流的方式為Intensity and Hydraulic Diameter（給定湍流強(qiáng)度和水力直徑）。設(shè)置進(jìn)水口速度為2.4m/s，湍流強(qiáng)度值為4.2%，水力直徑為70mm。對(duì)模型以及設(shè)立的條件設(shè)置9500步長進(jìn)行收斂性檢驗(yàn)，模型殘差圖如圖5所示，圖中共有6條數(shù)據(jù)線分別代表：質(zhì)量殘差、x、y、z三個(gè)方向速度殘差、湍流動(dòng)能殘差k、湍動(dòng)能耗散率殘差ε。通過計(jì)算入水口與出水口的質(zhì)量流量差結(jié)果為0.044%小于0.5%，所以數(shù)據(jù)結(jié)果為收斂。至此仿真模型的前處理和邊界條件設(shè)置完畢，模型收斂程度檢驗(yàn)完成。接下來進(jìn)行EDEM-Fluent耦合仿真求解計(jì)算。

3 仿真結(jié)果與優(yōu)化分析

3.1 噴頭形狀對(duì)射流強(qiáng)度的影響

通過Fluent軟件對(duì)模型進(jìn)行多次仿真計(jì)算得到收斂后的平均數(shù)據(jù)結(jié)果。為了試驗(yàn)的有效性，先對(duì)3種數(shù)據(jù)結(jié)果的進(jìn)出口流量差進(jìn)行計(jì)算，得到方形出水口噴頭的進(jìn)出口流量差為6.8×10-4 kg/s，圓形出水口噴頭的進(jìn)出口流量差為2.9×10-4 kg/s，帶倒角圓形噴頭出水口流量差為1.4×10-4 kg/s，都遠(yuǎn)小于他們的總流量，滿足質(zhì)量受恒定律，驗(yàn)證了參數(shù)設(shè)定及模型建立的正確性。三種形狀噴頭的中心對(duì)稱面處壓力、速度分布云圖如圖6所示。由圖6可以看到，三種噴頭除了進(jìn)水口附近有較強(qiáng)水流外，其余內(nèi)部流場分布均勻且流速低。對(duì)于方形出水口其速度和壓強(qiáng)的損耗相比圓形出水口要大的多。給噴頭從左到右依次編號(hào)1～5，3種不同形狀噴頭出水口的最大流速和平均流速數(shù)據(jù)統(tǒng)計(jì)結(jié)果如表2所示，方形噴頭5個(gè)出水口的最大速度和平均速度的平均值分別是4.77m/s和4.58m/s，標(biāo)準(zhǔn)差分別為0.34m/s與0.32m/s，圓形噴頭五個(gè)出口的最大和平均速度的平均值為7.15m/s和5.82m/s，最大和平均標(biāo)準(zhǔn)差為0.43m/s與0.16m/s，帶倒角圓形出水口5個(gè)出口的最大和平均速度的平均值是6.19m/s和5.66m/s，標(biāo)準(zhǔn)差為0.28m/s與0.20m/s。

分析壓力分布云圖與不同孔形狀出水口最大速度、平均速度對(duì)比表得出水流壓力波動(dòng)較大的時(shí)間為剛進(jìn)水的時(shí)候，之后就趨于平緩，在其他條件都相同的情況下，不同孔的形狀對(duì)水流速度的影響較大，比較方形孔噴頭與圓形孔噴頭，圓形孔噴頭能量損失更小。為了更明顯地看到水流波動(dòng)情況，繪制各噴頭水流速度波動(dòng)圖，如圖7所示，發(fā)現(xiàn)在噴頭與橫水管連接處使用倒角工藝，能使各噴頭水流速度更加穩(wěn)定，波動(dòng)更小。

由平均速度云圖可以看出，各噴頭平均流速波動(dòng)趨勢關(guān)于進(jìn)水口位置呈對(duì)稱分布，各噴頭平均流速的最大值與最小值的差值達(dá)到0.42m/s，波動(dòng)較大。為了更明顯的看到水流的波動(dòng)，將各噴頭的流速標(biāo)準(zhǔn)差做成折線圖（圖7），可以看出，各噴頭流速波動(dòng)的范圍大小，在雷諾數(shù)保持不變的情況下，出水口流速波動(dòng)大的噴頭沖刷的區(qū)域容易因?yàn)闆_刷力過大而導(dǎo)致蓮藕的表皮破損，影響了蓮藕的出售品質(zhì)；而出水口流速小的噴頭沖刷區(qū)容易因?yàn)闆_刷力過小而導(dǎo)致淤泥沒有被沖開，蓮藕漏采的現(xiàn)象，造成采收率低的結(jié)果。

3.2 調(diào)整噴頭孔徑對(duì)射流穩(wěn)定性的影響

查閱文獻(xiàn)得知，在壓強(qiáng)、水管內(nèi)徑不變的情況下，存在最合適的噴頭內(nèi)徑使射流強(qiáng)度最大，除此外射流強(qiáng)度逐漸減小。在一定直徑范圍內(nèi)，噴嘴出口直徑的增大可以使密集段長度和切割能力增加［16］。對(duì)于單入水口的橫水管，保持水泵壓強(qiáng)不變，通過試驗(yàn)發(fā)現(xiàn)調(diào)整噴頭直徑可以使5個(gè)出水口的壓強(qiáng)與速度變得更平穩(wěn)與均勻。其中噴頭直徑增減調(diào)整如表3所示。

由圖8可以看出，通過調(diào)整各噴頭的噴頭直徑可以有效提高噴頭流速的穩(wěn)定性。其中第5組的流速更均勻與穩(wěn)定，符合本次優(yōu)化目的。將標(biāo)準(zhǔn)組和第五組作為試驗(yàn)對(duì)照組進(jìn)行仿真分析，為了方便觀察，改設(shè)入水口流速為10m/s進(jìn)行EDEM-Fluent耦合仿真。因其結(jié)構(gòu)對(duì)稱性，取其中3個(gè)噴頭從外側(cè)到中間依次標(biāo)號(hào)1～3，其噴頭間挖掘深度差如圖9所示。

其中標(biāo)準(zhǔn)組的平均挖掘深度為33.50cm，挖掘深度標(biāo)準(zhǔn)差為1.71cm；第5組的平均挖掘深度為36.07cm，挖掘深度標(biāo)準(zhǔn)差為0.69cm。仿真分析可得，通過調(diào)整各噴頭的噴孔內(nèi)徑，有效的減小了能量的損耗，加大了挖掘深度，并且各噴頭間流速變得更加穩(wěn)定，挖掘深度標(biāo)準(zhǔn)差明顯減小，所以該方案符合優(yōu)化設(shè)計(jì)的目的，取得了一定的優(yōu)化效果。

4 結(jié)論

1） 對(duì)采藕機(jī)的總體方案進(jìn)行設(shè)計(jì)，利用SolidWorks對(duì)蓮藕除莖、挖掘、收集一體機(jī)各主要機(jī)構(gòu)進(jìn)行三維建模和虛擬裝配，對(duì)挖掘裝置中的射流系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì)。

2） 運(yùn)用Fluent軟件對(duì)射流系統(tǒng)的橫水管和高壓噴頭內(nèi)場流體域進(jìn)行仿真分析，對(duì)噴頭的形狀、直徑、連接工藝進(jìn)行流體力學(xué)仿真分析，改變噴頭形狀會(huì)大幅改變流速以及能量的損耗，在橫水管與噴頭連接處使用倒角工藝能有效減小能量的損失，使水流更穩(wěn)定，波動(dòng)更小。

3） 基于EDEM-Fluent對(duì)調(diào)整各噴頭內(nèi)徑后作用于土壤的整個(gè)過程進(jìn)行耦合仿真，仿真發(fā)現(xiàn)通過調(diào)整各噴頭的內(nèi)徑大小，能有效地使各噴頭的流速更均勻與穩(wěn)定，更有力地保障蓮藕的采收效率和質(zhì)量。流速波動(dòng)標(biāo)準(zhǔn)差由0.164m/s改變?yōu)?.020m/s；平均挖掘深度由33.50cm改變?yōu)?6.07cm，挖掘深度標(biāo)準(zhǔn)差由1.71cm改為0.69cm，在一定程度上解決蓮藕采收品質(zhì)參差不齊以及漏采率高的問題。

參 考 文 獻(xiàn)

［1］ 張佳， 李鋒霞， 潘江如， 等. 蓮藕收獲機(jī)械研究現(xiàn)狀與發(fā)展趨勢［J］. 中國農(nóng)機(jī)化學(xué)報(bào)， 2022， 43（10）： 12-17.

Zhang Jia， Li Fengxia， Pan Jiangru， et al. Research status and development trend of lotus root harvesting machinery ［J］. Journal of Chinese Agricultural Mechanization， 2022， 43（10）： 12-17

［2］ 謝晉， 韓迪， 王靖， 等. 中國蓮藕產(chǎn)業(yè)發(fā)展現(xiàn)狀及展望［J］. 農(nóng)業(yè)展望， 2017， 13（12）： 42-45， 51.

Xie Jin， Han Di， Wang Jing， et al. Development status quo and prospect of China’s lotus root industry ［J］. Agricultural Outlook， 2017， 13（12）： 42-45， 51.

［3］ 余有成， 郝應(yīng)昌， 王黨黨. 北方蓮藕種植技術(shù)［J］. 陜西農(nóng)業(yè)科學(xué)， 2014， 60（1）： 120-122.

［4］ 吳雪峰， 盧秉福， 韓衛(wèi)平. 農(nóng)業(yè)機(jī)械化制度變遷的博弈分析［J］. 中國農(nóng)機(jī)化學(xué)報(bào)， 2013， 34（1）： 31-34.

Wu Xuefeng， Lu Bingfu， Han Weiping. Game analysis of institutional innovation of agricultural mechanization in China ［J］. Journal of Chinese Agricultural Mechanization， 2013， 34（1）： 31-34.

［5］ 王計(jì)劃， 王艾倫， 楊誠， 等. 基于等效流體模型的蓮藕超聲波探測技術(shù)研究［J］. 中國農(nóng)機(jī)化學(xué)報(bào)， 2020， 41（5）： 125-133.

Wang Jihua， Wang Ailun， Yang Cheng， et al. Research on ultrasonic detection technology of lotus root based on equivalent fluid moded ［J］. Journal of Chinese Agricultural Mechanization， 2020， 41（5）： 125-133

［6］ 郭洋民. 水田蓮藕采挖機(jī)設(shè)計(jì)與研究［D］. 武漢： 華中農(nóng)業(yè)大學(xué)， 2015.

Guo Yangming. Design and research of lotus root digging machine in paddy field ［D］. Wuhan： Huazhong Agricultural University， 2015.

［7］ 夏文水， 姜啟興， 張家驪. 蓮藕方便食品加工技術(shù)的研究［J］. 食品與機(jī)械， 2007（1）： 139-142.

Xia Wenshui， Jiang Qixing， Zhang Jiali. Development of processing technology of convenient food from btus root ［J］. Food and Machinery， 2007（1）： 139-142.

［8］ 盧秉福. 黑龍江省農(nóng)業(yè)機(jī)械化發(fā)展與農(nóng)村剩余勞動(dòng)力轉(zhuǎn)移互動(dòng)性研究［J］. 中國農(nóng)機(jī)化學(xué)報(bào)， 2014， 35（3）： 268-271.

Lu Bingfu. Study on interaction between agricultural mechanization development and rural surplus labor transfer in Heilongjiang Province ［J］ Journal of Chinese Agricultural Mechanization， 2014， 35（3）： 268-271

［9］ 季大淼， 楊沿喜. 挖藕機(jī)［J］. 糧油加工與食品機(jī)械， 1982： 48-49.

［10］ 劉義滿， 柯衛(wèi)東， 黃新芳. 蓮藕人工采挖和機(jī)械采挖技術(shù)［J］. 長江蔬菜， 2014： 10-15.

［11］ 王維. 4SWO-1.2型船式水力挖藕機(jī)的研制［D］. 南京： 南京農(nóng)業(yè)大學(xué)， 2009.

Wang Wei. Researth of 4SWO-1.2 boat-type water power lotus root digging machine ［D］. Nanjing： Nanjing Agricultural University， 2009.

［12］ 吳昊， 夏俊芳， 張國忠， 等. 基于EDEM-Fluent仿真的自旋射流式挖藕機(jī)的設(shè)計(jì)與試驗(yàn)［J］. 農(nóng)業(yè)工程學(xué)報(bào)， 2018， 34（5）： 9-14.

Wu Hao， Xia Junfang， Zhang Guozhong， et al. Design and experiment of spin-jet flow type lotus root digging machine based on EDEM-Fluent ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2018， 34（5）： 9-14.

［13］ Gilabert F A， Roux J N， Castellanos A. Computer simulation of model cohesive powders： Influence of assembling procedure and contact laws on low consolidation states ［J］. Physical Review E， 2007， 75（1）： 011303.

［14］ 賀一鳴， 向偉， 吳明亮， 等. 基于堆積試驗(yàn)的壤土離散元參數(shù)的標(biāo)定［J］. 湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2018， 44（2）： 216-220.

He Yiming， Xiang Wei， Wu Mingliang， et al. Parameters calibration of loam soil for discrete element simulation based on the repose angle of particle heap ［J］. Journal of Hunan Agricultural University （Natural Science Edition）， 2018， 44（2）： 216-220.

［15］ 熊莉芳， 林源， 李世武. k-e湍流模型及其在FLUENT軟件中的應(yīng)用［J］. 工業(yè)加熱， 2007（4）： 13-15.

Xiong Lifang， Lin Yuan， Li Shiwu. k-e Turbulent model and its application to the FLUENT ［J］. Industrial Heating， 2007（4）： 13-15.

［16］ 梁煒明， 許鵬， 黃齊兵， 等. 高壓水切割混凝土影響因素試驗(yàn)分析［J］. 湖南交通科技， 2021， 47（4）： 99-103， 147.

Liang Weiming， Xu Peng， Huang Qibing， et al. Experimental analysis of influence factors of high pressure water cutting concrete ［J］. Hunan Communications Science and Technology， 2021， 47（4）： 99-103， 147.