仿旗魚頭部曲線型開溝器設(shè)計與性能試驗

趙淑紅，劉宏俊，譚賀文，曹秀振，張先民，楊悅乾

仿旗魚頭部曲線型開溝器設(shè)計與性能試驗

趙淑紅，劉宏俊，譚賀文，曹秀振，張先民，楊悅乾※

（東北農(nóng)業(yè)大學工程學院，哈爾濱150030）

為了減小播種開溝器的工作阻力，降低其對土壤的擾動，根據(jù)旗魚頭部的流線型曲線，建立了擬合曲線的數(shù)學模型，設(shè)計了仿旗魚頭部曲線型開溝器。為獲得該開溝器工作時工作阻力及其對土壤擾動情況，進行離散元仿真和室內(nèi)土槽試驗。在前進速度為0.5 m/s時，分別以3種土壤含水率、3種開溝器開溝深度為因素進行仿真及室內(nèi)試驗。仿真結(jié)果：在含水率為12%±1%時，工作阻力、對土壤擾動寬度和回土深度隨開溝深度的增加呈增加趨勢；在開溝深度為60 mm時，工作阻力隨含水率的增加而增加，但對土壤擾動寬度、回土深度影響不明顯。工作阻力仿真與試驗結(jié)果對比，兩者變化趨勢大致相同，相對誤差在5.04%～27.08%之間。土壤擾動情況仿真與試驗結(jié)果對比：回土深度的仿真與試驗結(jié)果接近，而土壤擾動寬度的仿真與試驗結(jié)果差別較大，可能是土壤顆粒粒徑較大引起的，但變化趨勢基本一致。由此說明了采用離散元法模擬開溝器的工作阻力及其對土壤擾動情況的可行性。該文可為開溝器的設(shè)計與優(yōu)化提供參考。

離散單元法；仿生；設(shè)計；曲線型開溝器；仿真；試驗

0 引 言

開溝器是播種機上用于開溝的關(guān)鍵部件[1-2]。一個具有較好結(jié)構(gòu)的播種開溝器，不僅能有效減少對土壤結(jié)構(gòu)的破壞，避免造成較大的失墑，還能夠減少其受到的工作阻力[3-4]。中國的播種機上，開溝器多選用尖角式和圓盤式[5-7]。尖角式開溝器較多用于中小型播種機上，但有時其對土壤擾動較大；而圓盤開溝器的土壤擾動性較小，但需要的機具質(zhì)量較大，不適合中小型播種機[8-9]。因此，優(yōu)化設(shè)計開溝器的結(jié)構(gòu)參數(shù)使之更加符合農(nóng)藝要求，不僅能夠降低開溝器的工作阻力，還能夠提高其入土能力，減輕對土壤的擾動。

目前對于播種開溝器的研究多集中于通過試驗的方法來改進設(shè)計開溝器[10-13]；為研究開溝器的性能，首先必須探討開溝器工作時土壤運動規(guī)律和施加于土壤的作用力。離散元法(discrete element method，DEM)可用來模擬顆粒材料和研究材料間的微觀宏觀變形，允許顆粒材料間存在接觸的形成和破壞，也適合仿真土壤和剛性體間的相互作用。學者基于DEM對土壤與耕作機具間的作用過程進行了廣泛研究[14-18]，這些研究都證實離散元仿真能夠模擬耕作過程。在離散元中可以追蹤顆粒的運動，進而對土壤顆粒在耕作過程中的受力和擾動情況進行分析。在對開溝器與土壤相互作用的研究中，于建群等[19]通過對芯鏵式開溝器工作阻力的離散元仿真，但其建立的開溝器和土壤顆粒模型為二維模型過于簡單；Chen等[20]通過建立針對3種不同土壤的三維離散元模型并模擬開溝器與土壤間的相互作用，驗證了模型的準確性；Ucgul等[21-22]分別針對無粘性土壤、粘性土壤采用不同的離散元接觸模型并對開溝過程進行離散元仿真，驗證了離散元模型的可行性；姜珊珊[23]和潘世強[24]運用三維離散元軟件EDEM（enhanced discrete element method）對芯鏵式開溝器的工作阻力進行了仿真與分析，并與試驗結(jié)果進行了對比，通過改進修正離散元模型，優(yōu)化設(shè)計了滑刀式開溝器。這些研究都證實離散元仿真能夠模擬開溝器與土壤相互作用，研究多集中于開溝器的工作阻力仿真上，而開溝器對土壤的擾動情況的研究較少。因此，運用離散元法對開溝器工作時受到的工作阻力及其對土壤擾動情況進行仿真分析，為播種開溝器設(shè)計與改進、提高效率、節(jié)省成本具有不可忽視的現(xiàn)實意義。

大自然已經(jīng)經(jīng)歷了數(shù)十億年的進化，他的鬼斧神工常常蘊藏著精妙的設(shè)計思想，旗魚是短距離游泳速度最快的動物，其頭部具有良好的流線型，在其高速前進時，降低阻力。本文依據(jù)旗魚頭部的流線型曲線，建立擬合曲線的數(shù)學模型，設(shè)計仿旗魚頭部曲線型開溝器，確定開溝器的結(jié)構(gòu)參數(shù)，運用離散元法對開溝器工作時受到的工作阻力及土壤擾動情況進行仿真分析，為播種開溝器的設(shè)計提供可行的研究方法，并進行室內(nèi)土槽試驗，進一步確定開溝器設(shè)計的合理性，為開溝器的設(shè)計與優(yōu)化提供參考。

1 仿旗魚頭部曲線型開溝器的設(shè)計

1.1 擬合曲線數(shù)學模型的建立

根據(jù)文獻[25]獲取了旗魚的三維圖形，如圖1a所示，然后利用CATIA二維制圖獲取旗魚頭部輪廓曲線，并記錄了曲線相對x、y軸參考線位置坐標，如圖1所示。根據(jù)相對位置坐標點采用分段多項式曲線擬合原理[26]進行數(shù)學分析，對每相鄰的4個點依次擬合，分段后的輪廓曲線，如圖1b所示。

圖1 旗魚及其頭部擬合曲線Fig.1 Sailfish and fitting curve of its head

在選取分段點時，為保證分段后相鄰曲線邊界點連續(xù)，采用3或4個數(shù)據(jù)點擬合一條二次曲線。由圖1b可知，首先對輪廓線分成了8段，一共24個數(shù)據(jù)點，然后再對各分段中的數(shù)據(jù)點進行曲線擬合，得到擬合曲線數(shù)學模型，如表1所示。

從分段曲線擬合方程中可以看出，決定系數(shù)R2都大于0.99，擬合程度較高，為開溝器的設(shè)計提供了可靠的數(shù)學模型。

1.2 開溝器的結(jié)構(gòu)參數(shù)的確定

根據(jù)東北地區(qū)播種農(nóng)藝開溝要求，結(jié)合上述擬合曲線數(shù)學模型相關(guān)參數(shù)，確定開溝器結(jié)構(gòu)尺寸。開溝器由開溝器入土部分、鏟柄組成，如圖2a所示。開溝器入土部分結(jié)構(gòu)示意圖，如圖2b、2c所示。根據(jù)擬合曲線數(shù)學方程將開溝器入土部分的入土曲線設(shè)計為擬合曲線型。根據(jù)農(nóng)藝要求及相關(guān)文獻[10-13]可知：1）入土角過大易造成干濕混合；過小，則強度不夠，壽命降低，一般取值范圍為20°～40°；2）入土隙角過小，入土性能變差；過大，土壤提前回落，使種子和肥料不能達到開溝的最大深度，入土隙角一般為5°～14°；3）當橫向?qū)挾刃∮阽P柄的橫向?qū)挾葧r，能夠明顯的降低對土壤擾動量；4）長度不能過長，否則不能保證整個開溝部分沉沒于土中，且對土壤擾動量較大；另一方面開溝器太長，易降低播種機通過性能。根據(jù)上述原則，確定入土角為30°，入土隙角為6°，入土部分長度為175 mm，橫向?qū)挾葹?0 mm，上下寬度相同，材料為65 Mn。鏟柄采用30 mm×50 mm的45號鋼管，使開溝器鏟柄與導(dǎo)種管集于一身，其長度為500 mm。

表1 擬合曲線數(shù)學模型Table1 Mathematical model of fitting curve

圖2 仿旗魚頭部曲線型開溝器結(jié)構(gòu)示意圖Fig.2 Structure diagram of opener based on bionic sailfish head curve

2 開溝器的離散元仿真

開溝器與土壤的相互作用是一個復(fù)雜的運動過程，進行離散元仿真可以從宏觀和細觀角度來研究開溝器與土壤的相互作用。

2.1 土壤離散元模型建立

本文選取含水率和開溝深度為試驗因素，建立不同土壤狀態(tài)及部件工況條件下的離散元顆粒模型。其中土壤含水率因素水平選取與春季播種時期土壤含水率相近，同時在后續(xù)試驗驗證階段將以此含水率的土壤為作業(yè)載體，控制土壤含水率分別為12%±1%、16%±1%、20%± 1%，測得土壤體積密度分別為1 280、1 300、1 330 kg/m3；通過土壤直接剪切試驗測得土壤內(nèi)摩擦系數(shù)為0.536、0.414、0.329，如表2所示；內(nèi)聚力為60～80 kPa；選用Hertz-Mindlin with bonding模型作為土壤顆粒間接觸力學模型[27-28]，用于模擬土壤顆粒間的相互作用。

表2 標定后土壤顆粒密度與顆粒內(nèi)摩擦系數(shù)Table2 Calibrated particle density and internal friction coefficient

2.1.1 參數(shù)校準

真正土壤顆粒太小，因而不可能直接測量其參數(shù)。由于仿真軟件EDEM和仿真時間不可能模擬出真實的不同尺寸和形狀的土壤顆粒。因而需要對仿真顆粒的細觀參數(shù)進行校準，以便給出更接近實際的顆粒模型作為土壤模型。

1）顆粒結(jié)構(gòu)參數(shù)

通過篩分法得出75%的土壤顆粒尺寸介于0.25～5 mm，其它25%小于等于0.25 mm。

通過查閱資料[29]可知，肥沃土壤的結(jié)構(gòu)形態(tài)由土壤粘結(jié)成粒狀和小團塊狀，大體呈球形；文獻[20-22]將土壤顆粒都設(shè)置成10 mm以上的球形，本文為了將仿真時間控制在一個合理的范圍內(nèi)，將土壤顆粒模型設(shè)置成粒徑2～5 mm之間的球形顆粒，在顆粒接觸模型中可以通過設(shè)置粘結(jié)接觸使所有顆粒粘 結(jié)在一起，從而生成其他形狀的顆粒集合體。

2）土壤顆粒密度

土壤顆粒密度也稱為土壤固相密度，與土壤體積密度不同。在離散元仿真分析時，需要調(diào)整土壤顆粒密度，保證仿真土壤體積密度與實際土壤體積密度相一致。

校準方法：土壤顆粒生成后使其在重力作用下達到穩(wěn)定，然后在模型上方施加一定壓力使其壓實后，慢慢釋放壓力，使顆粒再次穩(wěn)定。測得此時土壤顆粒的總質(zhì)量，得到仿真土壤體積密度，并與實際土壤體積密度相比較，如果不一致，再次調(diào)整土壤顆粒密度，直到生成的仿真土壤體積密度與實際土壤體積密度相同。校準后土壤顆粒密度見表2。

3）土壤顆粒間動靜摩擦系數(shù)

在離散元仿真分析時，需要調(diào)整土壤顆粒間動靜摩擦系數(shù)，保證仿真土壤內(nèi)摩擦系數(shù)與實際內(nèi)摩擦系數(shù)相一致。

內(nèi)摩擦系數(shù)的大小是內(nèi)摩擦角的正切值，在對生成的土壤顆粒間添加粘結(jié)模型前，土壤顆粒的休止角近似等于其內(nèi)摩擦角。通過模擬未添加粘結(jié)鍵的土壤顆粒的堆積試驗，對堆積后形成的傾角與測得的內(nèi)摩擦角進行對比（即正切比例關(guān)系）。仿真時，顆粒接觸模型為Hertz-Mindlin模型，顆粒粒徑大小為2～5 mm，顆粒個數(shù)為2 000個，時間步長為10-5s，顆粒生成在一個圓柱體內(nèi)，生成完畢后撤掉圓柱體，讓顆粒自然堆積，模擬結(jié)果，如圖3所示。

圖3 離散元模擬土壤顆粒不同含水率下的內(nèi)摩擦角Fig.3 Internal friction angle of soil particles with different moisture contents by discrete element simulation

采用校準土壤顆粒密度相同的方法，含水率為12%±1%、16%±1%、20%±1%時，土壤顆粒間動摩擦系數(shù)分別為0.25、0.22、0.20；土壤顆粒間靜摩擦系數(shù)分別為0.5、0.4、0.3；校準后與實際內(nèi)摩擦系數(shù)相差無幾，其說明參數(shù)校準的合理性，見表2。

2.1.2 土壤模型生成

在仿真時，用不同含水率的土壤參數(shù)分別對土槽建模；并采用EDEM的顆粒工廠生成土壤顆粒集合體。土壤顆粒及開溝器材料參數(shù)，見表3所示，其中需要靠經(jīng)驗取值的土壤參數(shù)參考文獻[24]。選取的顆粒粒徑大小為2～5 mm，顆粒數(shù)量為100 000個。在重力的作用下，顆粒生成完成后，再在顆粒上方加載校準土壤密度時施加的垂直載荷，壓實土壤模型，最后得到尺寸為700 mm× 300 mm×130 mm土壤模型，如圖4所示。

表3 離散元法仿真的基本參數(shù)表Table3 Basic parameters of discrete element method simulation

圖4 土壤模型Fig.4 Soil model

2.2 開溝器的仿真分析及結(jié)果

仿真分析前，在EDEM軟件中的前處理器模塊依次進行接觸力學模型、仿真參數(shù)、土壤的顆粒模型、開溝器邊界模型和顆粒工廠等的設(shè)置[30]。

結(jié)合東北地區(qū)實際田間播種土壤農(nóng)藝狀態(tài)，進行了預(yù)試驗，發(fā)現(xiàn)各因素之間的交互作用影響不大，故本文主要進行單因素試驗，旨在找到試驗因素（含水率與開溝深度）對開溝器工作性能的影響規(guī)律。根據(jù)文獻[23-24]確定試驗因素：土壤含水率為12%±1%、16%±1%、20%±1%，播種深度為30、60、90 mm。而在多數(shù)年份中，由于春季風大造成表層土壤含水率低，一般為12%±1%，且為了保證作物正常生長，播種時開溝深度多為60 mm，故重點在12%±1%的土壤含水率、60mm的開溝深度條件下進行試驗。仿真試驗方案：1）按照土壤含水率為12%±1%、工作速度為0.5 m/s、開溝深度為30、60、90 mm進行單因素仿真試驗；2）按照開溝深度為60 mm、工作速度為0.5 m/s、土壤含水率為12%±1%、16%±1%、20%±1%進行單因素仿真試驗。按照仿真試驗方案進行了模擬，并在EDEM后處理工具模塊進行仿真模擬的結(jié)果分析。測量溝形參數(shù)時，由于仿真的土壤模型顆粒比實際土壤顆粒要大，故選取溝形的關(guān)鍵參數(shù)：擾動寬度l和回土深度h進行測量，如圖5所示。自然回土后溝底基線與開夠深度基準線的距離為回土深度h。

圖5 模擬的開溝溝形與形狀參數(shù)Fig.5 Ditching shape and shape parameters in simulation

1）開溝深度對工作阻力的影響

采用土壤含水率為12%±1%、工作速度為0.5 m/s、開溝深度為30、60、90 mm進行仿真模擬，擬合曲線開溝器的工作阻力，如表4所示；仿旗魚頭部曲線型開溝器與土壤顆粒的接觸狀態(tài)，如圖6所示。由表4可知，開溝器的水平阻力在19.7～74.6 N，垂直阻力在9.37～26.9 N，工作阻力隨開溝深度增加而增加的趨勢；由圖6可知，土壤顆粒顏色由藍色至綠色至紅色依次加深，表示所受到的力依次增大，隨著開溝深度的增大，仿旗魚頭部曲線型開溝器入土部分上堆積的土壤顆粒逐漸變多，顆粒顏色變化加深，這說明了仿旗魚頭部曲線型開溝器對土壤的擾動在增大，其受到的工作阻力也在相應(yīng)變大。

表4 含水率和開溝深度對開溝器工作阻力的影響Table4 Effects of moisture content and ditching depth on working resistance of opener

圖6 含水率12%±1%時開溝器與土壤顆粒的接觸狀態(tài)Fig.6 Contact state between opener and soil particles of moisture content of 12%±1%

2）土壤含水率對工作阻力的影響

采用開溝深度為60 mm、工作速度為0.5 m/s、土壤含水率為12%±1%、16%±1%、20%±1%、進行仿真模擬，擬合曲線開溝器的工作阻力，如表4所示；仿旗魚頭部曲線型開溝器與土壤顆粒的接觸狀態(tài)，如圖7所示。由表4可知，開溝器的水平阻力在41.9～65.43 N，垂直反力在17.9～26.4 N，開溝器工作阻力隨著含水率的增加而增加；由圖7可知，土壤顆粒顏色由藍色至綠色至紅色依次加深，表示所受到的力依次增大，隨著含水率的增加，土壤顆粒顏色在加深，隨著土壤模型中的粘結(jié)力增大，這也說明了隨著含水率的增加，顆粒變化區(qū)域在加深，開溝器受到的工作阻力也在增加。

3）開溝深度對土壤擾動情況的影響

采用土壤含水率為12%±1%、工作速度為0.5 m/s、開溝深度為30、60、90 mm進行仿真模擬，土壤擾動情況進行分析。在仿旗魚頭部曲線型開溝器在3種開溝深度時，溝形的關(guān)鍵參數(shù)，如表5所示。

由表5可知，開溝器對土壤顆粒的擾動寬度在140～200 mm之間，回土深度在19～80 mm之間，隨著開溝深度的增加，開溝器對土壤的擾動寬度在增加，回土深度也在增加，說明開溝器對土壤的擾動在增大。

4）含水率對土壤擾動情況的影響

取開溝深度為60 mm、工作速度為0.5 m/s，仿旗魚頭部曲線型開溝器在不同的土壤含水率（12%±1%、16%±1%、20%±1%）時對土壤擾動情況進行分析。在3種土壤含水率下，溝形的關(guān)鍵參數(shù)，如表5所示。

圖7 開溝深度為60 mm時開溝器與土壤顆粒的接觸狀態(tài)Fig.7 Contact state between opener and soil particles of ditching depth of 60 mm

表5 含水率和開溝深度對溝形參數(shù)的影響Table5 Effects of moisture content and ditching depth on ditching shape parameters

隨著土壤含水率的增加，開溝器對土壤的擾動寬度以及回土深度變化不明顯，雖然在仿真參數(shù)設(shè)置時，隨著含水率的增加，顆粒間粘結(jié)力在變大，但當開溝器對土壤顆粒擾動后，顆粒間的粘結(jié)力消失，顆粒間的相互作用受到自身摩擦系數(shù)的影響，由于隨著含水率的增加，顆粒間的內(nèi)摩擦系數(shù)減小，這也說明了從溝形參數(shù)上不能準確的分析含水率對土壤顆粒擾動情況的影響。

3 仿真結(jié)果與試驗結(jié)果的對比

為了檢驗離散元法仿真結(jié)果的可信性及開溝器設(shè)計的合理性，在東北農(nóng)業(yè)大學農(nóng)具土槽實驗室內(nèi)進行試驗。

試驗方案：在土壤含水率為12%±1%，開溝深度為30、60、90 mm時，開溝器工作阻力及溝形參數(shù)的試驗值與仿真值進行對比；在開溝深度為60 mm、土壤含水率為12%± 1%、16%±1%、20%±1%時，開溝器工作阻力的試驗與仿真結(jié)果進行對比；并與圓弧型開溝器進行了對比試驗。

試驗條件：土槽土壤為東北地區(qū)典型作業(yè)黑壤土；試驗裝置[31]由試驗臺車、開溝器連接測試系統(tǒng)、仿旗魚頭部曲線型開溝器和驅(qū)動電機等組成，如圖8所示；試驗臺車動力由變頻柜控制電動機提供，主要儀器設(shè)備：變頻器（型號F1000-G055T3C，歐瑞傳動電氣有限公司生產(chǎn)）、三相異步電動機（型號Y2-10L2-4，永策機械設(shè)備有限公司生產(chǎn)）。試驗時將土槽分為啟動區(qū)、測試區(qū)和停止區(qū)3部分，其中穩(wěn)定測試工作區(qū)3 m，前后啟動和停止區(qū)各0.5 mm。

圖8 對比試驗臺架Fig.8 Contrast experiment bench

3.1 工作阻力的對比

3種含水率、3種開溝深度下工作阻力的仿真與試驗結(jié)果對比結(jié)果，如表6所示。

表6 不同含水率和開溝深度下工作阻力仿真與試驗結(jié)果對比Table6 Comparison of opener simulation and experimental results with different moisture content and ditching depth

由表6可知，工作阻力仿真與試驗結(jié)果都隨著含水率及開溝深度的增加而增加，且變化趨勢大致相同。隨著含水率和開溝深度增加，相對誤差有減小的趨勢，相對誤差在5.04%～27.08%之間，由此證明離散元法對工作阻力仿真的可信性。

3.2 土壤擾動情況的對比

對仿真與試驗測得的溝形形狀參數(shù)進行對比，來對比分析仿旗魚頭部曲線型開溝器對土壤的擾動情況。在土壤含水率為12%±1%時，仿旗魚頭部曲線型開溝器在開溝深度為30、60、90 mm時的溝形參數(shù)進行分析，如表7所示。

表7 土壤含水率為12%±1%時的溝形參數(shù)Table7 Ditching shape parameters of moisture content of 12%±1%

由表7可知，隨著開溝深度增加，開溝器對土壤的擾動寬度及回土深度也增加?；赝辽疃确抡娼Y(jié)果與試驗結(jié)果幾乎相同，之間差值在-2.5～0.7 mm，但擾動寬度仿真結(jié)果要比試驗結(jié)果小106.7～140.0 mm，這可能與仿真顆粒粒徑較大有關(guān)，但變化趨勢基本一致，由此確定離散元法對土壤擾動情況仿真的可行性。

3.3 對比試驗

為了驗證擬合曲線性開溝器工作性能，同時保證試驗的可比性，本研究選擇圓弧型開溝器進行了對比試驗。試驗選擇土壤含水率為12%±1%，開溝深度為60 mm，工作速度為0.5 m/s。

由表8可知，仿旗魚頭部曲線型開溝器工作阻力比圓弧型開溝器??；仿旗魚頭部曲線型開溝器的擾動寬度小且回土深度深，這說明仿旗魚頭部曲線型開溝器工作性能較比圓弧型開溝器要好。由上述可知，仿旗魚頭部曲線型開溝器有一定的優(yōu)越性，可以降低工作阻力與土壤擾動，具有一定的應(yīng)用價值。

表8 開溝器類型工作性能對比Table8 Working performance comparison of openers

4 討 論

1）本文運用離散元數(shù)值仿真分析仿旗魚頭部曲線型開溝器工作過程，以工作阻力和土壤擾動為目標，通過開溝器工作阻力數(shù)據(jù)、開溝器與土壤模型接觸狀態(tài)、溝形形狀參數(shù)的分析，定性并定量的分析開溝器在不同含水率土壤模型和不同開溝深度工況下的作業(yè)狀態(tài)，并對比研究此類型開溝器特點[19,23-24]，為開溝器優(yōu)化改進及效率提高提供參考。

2）由仿真與試驗結(jié)果對比可知，仿真結(jié)果均小于試驗結(jié)果。文獻[19,24]在采用離散元研究開溝器作業(yè)性能時中也發(fā)現(xiàn)仿真結(jié)果比試驗結(jié)果小，可能原因是在離散元仿真時，仿真模型的土壤顆粒形狀和數(shù)目與實際情況存在一定差距，但仿真與試驗結(jié)果變化趨勢基本一致，這表明離散元仿真分析驗證開溝器工作性能的可行性。

3）在仿真分析含水率對土壤擾動情況時，含水率對土壤擾動情況影響不顯著，這與實際情況不吻合，該結(jié)果表明本研究所選用土壤-機械間接觸力學模型在模擬不同含水率對土壤擾動情況時存在一定的局限性。在后續(xù)研究中將運用C語言對彈性函數(shù)進行編譯，通過EDEM軟件的應(yīng)用編程接口（API，application programming interface）進行軟件二次開發(fā)，使其更精準地進行離散元仿真模型，以期使能保證在含水率不同對土壤的擾動能與實際情況吻合，同時還將研究開溝器在有無秸稈覆土壤環(huán)境下的離散元模型，研究秸稈對土壤擾動的影響。

5 結(jié) 論

1）本研究依據(jù)旗魚頭部的流線型曲線，采用最小二乘法分段曲線擬合原理，建立了擬合曲線數(shù)學模型，設(shè)計了仿旗魚頭部曲線型開溝器，其能提高機具通過性，減小工作阻力，降低土壤擾動，同時具有良好的回土性能。

2）由仿真與試驗結(jié)果對比可知，開溝器工作阻力仿真與試驗結(jié)果的變化趨勢大致相同，工作阻力都隨開溝深度的增加而增加，隨土壤含水率的增大而增大，其中仿真均比試驗值小，相對誤差在5.04%～27.08%之間。

3）由對比試驗可知，仿旗魚頭部曲線型開溝器可降低工作阻力，并能夠有效的減小土壤擾動，進一步說明開溝器設(shè)計結(jié)構(gòu)的合理性。

[1] 付乾坤，薦世春，賈洪雷，等. 玉米滅茬起壟施肥播種機的設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報，2016，32(4)：9－16.

Fu Qiankun, Jian Shichun, Jia Honglei, et al. Design and experiment on maize stubble cleaning fertilization ridging seeder[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(4): 9－16. (in Chinese with English abstract)

[2] 王慶杰，趙宏波，何進，等. 組合刀式固定道小麥免少耕播種機設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報，2016，32(17)：12－17.

Wang Qingjie, Zhao Hongbo, He Jin, et al. Design and experiment of blades-combined no and minimum-till wheat planter under controlled traffic farming system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(17): 12－17. (in English with Chinese abstract)

[3] 蔣金琳，馬廷璽. 免耕播種機土壤工作部件評述[J]. 農(nóng)業(yè)機械學報，1996，10(27)：68－71.

Jiang Jinlin, Ma Tingxi. Review on the soil interacting tools of no-till planter[J]. Transactions of the Chinese Society for Agricultural Machinery, 1996, 10(27): 68－71. (in Chinese with English abstract)

[4] 王徐建，宋建農(nóng)，劉彩玲，等. 甘草傾斜移栽開溝器的設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報，2016，32(13)：16－23.

Wang Xujian, Song Jiannong, Liu Cailing, et al. Design and experiment on licorice tilt transplanting furrow opener[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(13): 16－23. (in Chinese with English abstract)

[5] 馬云海，馬圣勝，賈洪雷，等. 仿生波紋形開溝器減黏降阻性能測試與分析[J]. 農(nóng)業(yè)工程學報，2014，30(5)：36－41.

Ma Yunhai, Ma Shengsheng, Jia Honglei, et al. Measurement and analysis on reducing adhesion and resistance of bionic ripple opener[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(5): 36－41. (in Chinese with English abstract)

[6] 李寶筏. 農(nóng)業(yè)機械學[M]. 北京：中國農(nóng)業(yè)出版社，2003.

[7] 龔振平，楊悅乾. 作物秸稈還田技術(shù)與機具[M]. 北京：中國農(nóng)業(yè)出版社，2012.

[8] 高煥文，李洪文，姚宗路. 輕型高防堵性能免耕播種機研究[J]. 中國工程科學，2007，9(9)：35－39.

Gao Huanwen, Li Hongwen, Yao Zonglu. Study on the light no-till seeders with high anti-blockage performance[J]. Engineering Sciences, 2007, 9(9): 35－39. (in Chinese with English abstract)

[9] 王繼山. 免耕播種機開溝器的分析與研究[J]. 山西農(nóng)業(yè)大學學報，2007，6(增刊1)：7－10.

Wang Jishan. Analysis and research on the drill opener of no tillage planter[J]. Journal of Shanxi Agricultural University, 2007, 6(Supp.1): 7－10. (in Chinese with English abstract)

[10] 趙淑紅，楊悅乾，閆以勛，等. 鴨嘴式開溝器的設(shè)計試驗[J]. 東北農(nóng)業(yè)大學學報，2013，44(11)：113－117.

Zhao Shuhong, Yang Yueqian, Yan Yixun, et al. Design of duck mouth type opener[J]. Journal of Northeast Agricultural University, 2013, 44(11): 113－117. (in Chinese with English abstract)

[11] 茍文，馬榮朝，楊文鈺，等. 小麥免耕播種機開溝器的設(shè)計[J]. 農(nóng)業(yè)工程學報，2012，28(增刊1)：21－25.

Gou Wen, Ma Rongchao, Yang Wenyu, et al. Design of opener on no-till wheat seeder[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(Supp.l): 21－25 (in Chinese with English abstract)

[12] 顧耀權(quán)，賈洪雷，郭慧，等. 滑刀式開溝器設(shè)計與試驗[J]．農(nóng)業(yè)機械學報，2013，44(2)：38－42．

Gu Yaoquan, Jia Honglei, Guo Hui, et al. Design and experiment of sliding knife furrow opener[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(2): 38－42. (in Chinese with English abstract)

[13] 莊健，賈洪雷，馬云海，等.具有滑刀式缺口的圓盤開溝器設(shè)計與試驗[J]. 農(nóng)業(yè)機械學報，2013，44(增刊1)：83－88.

Zhuang Jian, Jia Honglei, Ma Yunhai, et al. Design and experiment of sliding-knife-type disc opener[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(Supp.l): 83－88. (in Chinese with English abstract)

[14] 張銳，李建橋，李因武. 離散單元法在土壤機械特性動態(tài)仿真中的應(yīng)用進展[J]. 農(nóng)業(yè)工程學報，2003，19(1)：16－19.

Zhang Rui, Li Jianqiao, Li Yinwu. Development of simulation on mechanical dynamic behavior of soil by distinct element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2003, 19(1): 16－19. (in Chinese with English abstract)

[15] 張銳，李建橋，李因武，等. 部件復(fù)雜表面影響土壤擾動行為的離散元宏細觀分析[J]. 吉林大學學報:工學版，2009，39(5): 1218－1223.

Zhang Rui, Li Jianqiao, Li Yinwu, et al. DEM macroscopic and mesoscopic analysis in disturbed behavior of soil acted by part with complex surface [J]. Journal of Jilin University :Engineering and Technology Edition, 2009, 39(5): 1218－1223. (in Chinese with English abstract)

[16] 徐泳，孫其誠，張凌，等. 顆粒離散元法研究進展[J]. 力學進展，2003，33(2)：251－260.

Xu Yong, Sun Qicheng, Zhang Ling, et al. Advances in discrete element methods for particulate materials[J]. Advances in Mechanics, 2003, 33(2): 251－260. (in Chinese with English abstract)

[17] Jiang M J，Yu H S，Harris D. Discrete element modeling of deep penetration in granular soils[J]. Numerical and Analytical Methods Geomechanics, 2006, 30(4): 271－305.

[18] Shmulevich I, Asaf Z, Rubinstein D. Interaction between soil and a wide cutting blade using the discrete element method[J]. Soil & Tillage Research, 2007, 97(1): 37－50.

[19] 于建群，錢立彬，于文靜，等. 開溝器工作阻力的離散元法仿真分析[J]. 農(nóng)業(yè)機械學報，2009，40(6)：53－57.

Yu Jianqun, Qian Libin, Yu Wenjing, et al. DEM analysis of the resistances applied on furrow openers[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(6): 53－57. (in Chinese with English abstract)

[20] Chen Y, Lars J M, Tavs N. A discrete element model for soil-sweep interaction in three different soils[J]. Soil & Tillage Research, 2013, 126: 34－41.

[21] Ucgul M, John M F, Chris S. 3D DEM tillage simulation: validation of a hysteretic spring (plastic) contact model for a sweep tool operation in a cohesionless soil[J]. Soil &Tillage Research, 2014, 144: 220－227.

[22] Ucgul M, John M F, Chris S. Three dimensional discrete element modeling (DEM) of tillage: accounting for soil cohesion and adhesion[J]. Biosystem Engineering, 2015, 129: 298－306.

[23] 姜珊珊. 基于離散元法的開溝器試驗研究[D]. 長春：吉林農(nóng)業(yè)大學，2014.

Jiang Shanshan. Experimental Research of Furrow Openers Based on DEM[D]. Changchun; Jilin Agricultural University, 2014. (in Chinese with English abstract)

[24] 潘世強. 基于離散元法的芯鏵式開溝器優(yōu)化設(shè)計與試驗研究[D]. 長春：吉林大學，2015.

Pan Shiqiang. Research on the Optimization design and the Experiment of the Core Ploughshare Furrow Opener Based on the Discrete Element Method[D]. Changchun: Jilin University, 2015. (in Chinese with English abstract)

[25] 齊文溥. 仿生設(shè)計在高速列車頭型、車燈設(shè)計中的應(yīng)用[D].成都：西南交通大學，2012.

Qi Wenbo. Bionic Design in the Head Type and Lights of High-speed Trains[M]. Chengdu: Southwest Jiaotong University, 2012. (in Chinese with English abstract)

[26] 劉霞，王運鋒. 基于最小二乘法的自動分段多項式曲線擬合方法研究[J]. 科學技術(shù)與工程，2014，14(3)：55－58.

Liu Xia, Wang Yunfeng. Research of automatically piecewise polynomial curve-fitting method based on leastsquare principle[J]. Science Technology and Engineering, 2014, 14(3): 55－58. (in Chinese with English abstract)

[27] 翟力欣，姬長英. 基于離散單元法的土壤力學接觸模型的建立[J]. 江西農(nóng)業(yè)學報，2008，20(9)：108－111.

Zhai Lixin, Ji Changying. Foundation of dynamic soil model based on distinct element method[J]. Acta Agriculturae Jiangxi, 2008, 20(9): 108－111. (in Chinese with English abstract)

[28] Potyondy D O，Cundall P A. A bonded-particle model for rock[J]. Journal of Rock Mechanics and Minimg Sciences.2004, 41(8): 1329－1364.

[29] 邵明安，王全九，黃明斌. 土壤物理學[M]. 北京：高等教育出版社，2006.

[30] 胡國明. 顆粒系統(tǒng)的離散元素法分析仿真[M]. 武漢：武漢理工大學出版社，2010.

[31] 王金武，張春鳳，周文琪，等. 基于MATLAB 仿真的非規(guī)則齒輪行星系扎穴機構(gòu)的優(yōu)化設(shè)計[J]. 農(nóng)業(yè)工程學報，2016，32(3)：22－28.

Wang Jinwu, Zhang Chunfeng, Zhou Wenqi, et al. Optimization design of non-circular planetary gear trains pricking hole mechanism based on MATLAB[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(3): 22－28. (in Chinese with English abstract)

Design and performance experiment of opener based on bionic sailfish head curve

Zhao Shuhong, Liu Hongjun, Tan Hewen, Cao Xiuzhen, Zhang Xianmin, Yang Yueqian※
(College of Engineering, Northeast Agricultural University, Harbin 150030, China)

In order to reduce working resistance of seeding opener and soil disturbance, the opener based on bionic sailfish head curve was designed. The two-dimensional (2D) drafting software of CATIA (computer aided three-dimensional interactive application) is used for getting the sailfish head profile curve. The head profile curve segments of sailfish is obtained based on the relative positional coordinates of points. Using the agronomic requirements and the parameters of the head profile curve segments of sailfish determine the structural parameters of opener. The main structure parameters of the opener are penetration angle of 30° and penetration clearance angle of 6°. Due to the interaction of the soil and the opener during movement is a complex process, the discrete element simulation from microscopic point of view is performed to study the interaction between the opener and the soil. The discrete element model of soil is established. In order to ensure the accuracy of test data, the related parameters are calibrated. Soil resistance is influenced by working depth and soil moisture content. At the same time, the relationship between working depth and soil disturbance is obtained. The software of EDEM (enhanced discrete element method) is used for the simulation analysis of working resistance and soil disturbance of the opener. Soil moisture content is 12%±1%, 16%±1%, and 20%±1%, working depth is 30, 60, and 90 mm, respectively, and forward speed of the opener is 0.5 m/s. Under these working conditions, the working resistance of the opener and its disturbance to soil are studied. Simulation results show that the working resistance of the opener increases with the increasing of soil moisture content, and also increases with the increasing of ditching depth, the horizontal resistance is between 19.70 and 74.60 N, the vertical resistance is between 9.37 and 26.90 N, the width of the disturbance in the opener is between 140 and 200 mm, and the soil returning depth is between 19 and 80 mm. In order to prove the credibility of the simulation results of the EDEM and the rationality of the opener, the experiments are carried out in the agricultural soil test laboratory of the northeast agricultural university. Simulation of work resistance and test results showed that the trend of change in work resistance was roughly the same, the relative error were between 5.04% and 27.08%. Compared with the experimental results, the simulation results were almost the same as the experimental results, and the difference of the ditching depth was between -2.5 and 0.67 mm, but the simulation results of the disturbance width were between 106.7 and 140 mm, smaller than the experimental results, which may be related to the larger particle size of the simulated particles. The trend was basically the same. The opener based on bionic sailfish head curve and the arc-type opener were selected for contrast experiment. The results showed that the performance of the opener based on bionic sailfish head curve was better than the arc-type opener. The research shows the feasibility of using the discrete element method to simulate opener’s working resistance and soil disturbance, and provides a convenient method for the opener design.

discrete element method; bionic; design; curve opener; simulation; experiment

10.11975/j.issn.1002-6819.2017.05.005

S225+2

A

1002-6819(2017)-05-0032-08

趙淑紅，劉宏俊，譚賀文，曹秀振，張先民，楊悅乾. 仿旗魚頭部曲線型開溝器設(shè)計與性能試驗[J]. 農(nóng)業(yè)工程學報，2017，33(5)：32－39.

10.11975/j.issn.1002-6819.2017.05.005 http://www.tcsae.org

Zhao Shuhong, Liu Hongjun, Tan Hewen, Cao Xiuzhen, Zhang Xianmin, Yang Yueqian. Design and performance experiment of opener based on bionic sailfish head curve[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(5): 32－39. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.05.005 http://www.tcsae.org

2016-05-29

2016-12-05

“十二五”國家科技支撐計劃項目（2014BAD11B01）；黑龍江省重大科技招標項目（GA14B101）；東北農(nóng)業(yè)大學博士啟動基金（2010RCB45）

趙淑紅，教授，博士生導(dǎo)師，主要從事田間農(nóng)業(yè)機械及力學特性的研究。哈爾濱 東北農(nóng)業(yè)大學工程學院，150030。Email：shhzh091@sina.com※通信作者：楊悅乾，研究員級工程師，主要從事保護性農(nóng)業(yè)裝備研究。哈爾濱 東北農(nóng)業(yè)大學工程學院，150030。Email：yangyueqian@126.com