楊慶璐 李洪文 何 進(jìn) 盧彩云 王英博 王慶杰
(1.中國農(nóng)業(yè)大學(xué)工學(xué)院, 北京 100083; 2.中國農(nóng)業(yè)大學(xué)現(xiàn)代農(nóng)業(yè)裝備優(yōu)化設(shè)計北京市重點(diǎn)實驗室, 北京 100083)
東北黑土區(qū)是我國重要的商品糧生產(chǎn)基地,在保障國家糧食安全和促進(jìn)國民經(jīng)濟(jì)發(fā)展中具有重要作用[1]。保護(hù)東北黑土區(qū)耕地資源安全,對于守住“谷物基本自給、口糧絕對安全”戰(zhàn)略底線具有重要意義[2]。但自東北黑土大規(guī)模開墾以來,由于長期過度墾殖、高強(qiáng)度產(chǎn)出和保護(hù)性措施缺失,導(dǎo)致黑土耕地開始明顯退化,土壤肥力逐年下降[3],黑土層明顯變薄、有機(jī)質(zhì)銳減、犁底層增厚,嚴(yán)重影響了東北黑土區(qū)農(nóng)業(yè)可持續(xù)發(fā)展[4-5]。解決東北黑土區(qū)耕地質(zhì)量退化的關(guān)鍵是變革耕作制度,將掠奪式開墾轉(zhuǎn)變?yōu)楸Wo(hù)性耕作。
玉米是我國重要的糧食作物,東北黑土區(qū)播種面積約占我國玉米種植總面積的30%[6]。玉米是喜肥作物,生育期長、產(chǎn)量高、需肥量大,合理施肥對提高玉米的產(chǎn)量和減少環(huán)境污染具有關(guān)鍵性作用[7-9]。不合理的施肥會造成化肥大量揮發(fā)、流失,利用率降低,并對環(huán)境、土壤產(chǎn)生危害[10-11]。
目前在東北地區(qū)玉米施肥仍以傳統(tǒng)分期施肥為主,一般采用基肥和追肥相結(jié)合的方式,不僅費(fèi)時、費(fèi)力,而且大部分農(nóng)民根據(jù)傳統(tǒng)經(jīng)驗追肥,淺施、撒施現(xiàn)象普遍,易造成肥料淋溶和揮發(fā),降低肥料利用率[12-13]。研究表明,合理的深松和肥料深施能改善土壤理化性質(zhì),促進(jìn)作物養(yǎng)分吸收和根系生長,提高肥料利用率,增加作物產(chǎn)量[14-16]。因此相關(guān)專家學(xué)者將深松和深施肥相結(jié)合,提出了機(jī)械化分層施肥技術(shù),即利用施肥機(jī)具將作物所需肥料一次性按比例分層施入土壤,后期不再追肥,能有效提高肥料利用率,同時減少施肥作業(yè)環(huán)節(jié),節(jié)本增效[17-18]。相關(guān)研究表明,分層施肥能夠促進(jìn)作物根系生長,增強(qiáng)養(yǎng)分吸收能力,增加作物產(chǎn)量,同時延長肥效,減少肥料流失[19-20]。黃幸媛[21]通過研究發(fā)現(xiàn)當(dāng)分層施肥的上、中、下3層的施肥量比例為3∶3∶4時,既不會造成“燒苗”,又能保證玉米的不同生長階段養(yǎng)分的供應(yīng),從而獲得更高的產(chǎn)量和經(jīng)濟(jì)效益。分層施肥機(jī)具研究方面,劉進(jìn)寶等[22]設(shè)計了一種對行開溝分層深施肥鏟,上下兩層施肥口前后分置,分別連接不同的排肥盒,實現(xiàn)了兩層肥量獨(dú)立可調(diào)。頓國強(qiáng)等[23]設(shè)計了一種由主開溝器及淺層尾管組成的分層施肥開溝器,通過淺層尾管的位置調(diào)節(jié),實現(xiàn)分層施肥位置的橫向控制。楊然兵等[24]結(jié)合曲面排肥與V形防堵結(jié)構(gòu),設(shè)計了一種深度可調(diào)的分層施肥開溝器,可在兩行種薯塊中間進(jìn)行開溝施肥作業(yè),利用開溝器自動回土實現(xiàn)肥料的覆土掩埋作業(yè)。
作物對肥料的利用不僅與施肥技術(shù)相關(guān),與肥料種類也有很大關(guān)系。傳統(tǒng)普通復(fù)合肥易溶解、肥效快,施用后能立即被作物吸收利用,但肥效持續(xù)時間短,大量施肥后,作物吸收不及時,會造成養(yǎng)分流失,破壞生態(tài)環(huán)境[25]。緩釋肥肥效持久穩(wěn)定,在作物生長期能持續(xù)提供養(yǎng)分,且不易導(dǎo)致土壤鹽分過高而“燒苗”[26],但由于養(yǎng)分釋放緩慢,不能在短時間內(nèi)提供大量養(yǎng)分。目前相關(guān)學(xué)者對分層施肥技術(shù)的研究,大多集中于單一種類肥料的施用,難以達(dá)到最佳的施肥效果。目前東北地區(qū)實行保護(hù)性耕作后,作物秸稈覆蓋還田,傳統(tǒng)的追肥作業(yè)難以實現(xiàn),因此本文根據(jù)東北黑土地玉米免少耕播種與施肥需求,提出一種普通復(fù)合肥與緩釋肥分層施用的方法,并設(shè)計一種分層深施肥裝置,利用計算流體動力學(xué)與離散元耦合仿真的方法,對裝置進(jìn)行參數(shù)優(yōu)化,以期實現(xiàn)普通復(fù)合肥與緩釋肥的分層施用,并形成明顯穩(wěn)定的分層,為氣力式排肥系統(tǒng)的設(shè)計與優(yōu)化提供技術(shù)參考和理論支撐。
氣力式排肥系統(tǒng)主要由分層深施肥裝置、肥料均勻分配裝置、肥料變比例分配裝置、氣-肥分離裝置、供肥裝置、氣-肥混合裝置、高壓風(fēng)機(jī)和輸送管道等組成[27]。其中肥料分配裝置分為肥料均勻分配裝置和肥料變比例分配裝置,包括肥料均勻分配器、肥料變比例分配器、波紋管、排肥管等,供肥裝置分為2套,分別對應(yīng)肥料均勻分配裝置和肥料變比例分配裝置,包括肥箱、排肥電機(jī)、外槽輪排肥器等。氣力式排肥系統(tǒng)結(jié)構(gòu)如圖1所示。
高壓風(fēng)機(jī)通過氣流輸送管道,分別與2個氣-肥混合裝置的進(jìn)風(fēng)口相連;2個外槽輪排肥器分別由2個排肥電機(jī)帶動,出口與氣-肥混合裝置進(jìn)肥口相連;2個氣-肥混合裝置出口經(jīng)氣-肥輸送管道與2個波紋管相連接;肥料均勻分配器入口與波紋管上端相連接,出口經(jīng)排肥管與分層深施肥鏟上層排肥管相連接;肥料變比例分配器入口與波紋管上端相連接,相鄰兩出口經(jīng)氣-肥分離裝置和排肥管,分別與分層深施肥鏟中層和下層排肥管相連接。
氣力式排肥系統(tǒng)利用高速氣流將肥料分配并輸送到分層深施肥鏟各排肥口,分層深施肥鏟將肥料分層施入土壤。工作過程主要分為供肥、混肥、分肥、排肥4個階段,系統(tǒng)工作時,高壓風(fēng)機(jī)產(chǎn)生高速氣流,經(jīng)氣流輸送管道分別進(jìn)入2個氣-肥混合裝置;2個肥箱分別對應(yīng)2個肥料分配裝置,肥箱Ⅰ對應(yīng)肥料均勻分配裝置,肥箱Ⅱ?qū)?yīng)肥料變比例分配裝置,2個肥箱中肥料種類不同,肥箱Ⅰ中放置普通復(fù)合肥,易溶解便于作物吸收利用,肥箱Ⅱ中放置緩釋肥料,肥效長可在作物生長期持續(xù)提供養(yǎng)分;不同種類肥料分別從對應(yīng)的外槽輪排肥器定量排出,依靠自身重力和氣-肥混合裝置產(chǎn)生的壓差進(jìn)入氣-肥混合裝置混合腔;氣-肥混合裝置利用文丘里效應(yīng),將高速氣流與肥料混合,形成氣-肥混合流;氣-肥混合流經(jīng)氣-肥輸送管道進(jìn)入波紋管,在波紋管的作用下,氣流與肥料均勻混合,形成均勻的氣-肥混合流;肥料均勻分配器將氣-肥混合流均勻分配至每一出口[28],經(jīng)排肥管從分層深施肥鏟上層排肥口排出;肥料變比例分配器將相鄰兩出口的肥料按預(yù)定比例分配[27,29],經(jīng)氣-肥分離裝置和排肥管,分別從分層深施肥鏟中層和下層排肥口排出;氣-肥分離裝置將分配完成后的氣-肥混合流中的部分氣流排出系統(tǒng),降低分層深施肥鏟排肥口處氣流速度,減少肥料彈跳,提高分層效果;分層深施肥鏟切開土壤,將肥料分層施入土壤中,且上層肥料為普通復(fù)合肥,中層和下層肥料為緩釋肥,實現(xiàn)普通復(fù)合肥和緩釋肥分層施用。
分層深施肥作業(yè)是將玉米整個生長周期所需的肥料,一次性按需分層施入土壤中,且肥料在土壤中呈半包圍種子的分布狀態(tài),上、中、下3層肥料呈“倒三角形”分布,與玉米根系形狀相吻合,更利于植株對養(yǎng)分側(cè)吸收,提高肥料利用率[21,30]。上層和中層肥料分施兩側(cè),避免“燒種”的同時,更利于玉米根系吸收。上層肥料為普通復(fù)合肥,可在短時間內(nèi)為玉米生長前期提供充足的養(yǎng)分;玉米生長的前期根系不發(fā)達(dá),無法觸及中層和下層肥料,即在玉米生長前期中層和下層肥料還未開始為植株提供養(yǎng)分,但在玉米生長前期這個階段,中層和下層的肥料已經(jīng)在播種同時施入土壤中,為避免肥料養(yǎng)分的流失,同時避免在玉米生長后期“脫肥”,因此中層和下層所施肥料為緩釋肥,能夠在玉米生長的中期和后期為植株提供穩(wěn)定充足的養(yǎng)分,增加玉米產(chǎn)量。玉米播種深度H1為30~50 mm,為避免“燒種”,同時考慮到分層深施肥鏟的結(jié)構(gòu)布置,上層肥料距離種子的橫向距離S1確定為120 mm,距離地表的深度H2為80~100 mm,上層肥料兩側(cè)分施,兩側(cè)肥料以種子所在的豎直垂線左右對稱。中層肥料距離上層肥料的深度H3為100 mm,即距離地表的深度為180~200 mm,與種子的橫向距離S2為80 mm,中層肥料也為兩側(cè)分施,兩側(cè)肥料以種子所在的豎直垂線左右對稱。下層肥料距離中層肥料的深度H4為100 mm,即距離地表的深度為280~300 mm,位于種子正下方。3層肥料呈半包圍種子的分布狀態(tài),如圖2所示。
分層深施肥裝置是實現(xiàn)肥料分層施用的關(guān)鍵部件,主要由深松鏟和肥料分層施肥器組成,其中深松鏟包括鏟尖、鏟柄和安裝座等部件;肥料分層施肥器包括上層、中層和下層排肥管,上層肥料左右兩側(cè)翼鏟和排肥口,中層肥料左右翼鏟和排肥口,下層肥料排肥口,上層和中層肥料Y形分配器,擋土板等部件,其結(jié)構(gòu)如圖3所示。
肥料均勻分配裝置出口與分層深施肥裝置的上層肥料排肥管相連接;肥料變比例分配裝置相鄰兩出口分別與分層深施肥裝置的中層和下層肥料排肥管相連接;上層和中層肥料排肥管在上層和中層肥料翼鏟處均連接有Y形分配器,將上層和中層排肥管分為左右兩側(cè),側(cè)翼安裝在鏟柄兩側(cè),將左右兩側(cè)排肥口包裹,避免堵塞;下層肥料排肥口位于鏟柄正后方,為避免排肥口堵塞,安裝有擋土板。為保證肥料順利排出并降低作業(yè)時的阻力,上層和中層肥料翼鏟向下、向后有一定的傾斜角度。
作業(yè)前,根據(jù)玉米分層施肥各層施肥質(zhì)量比為3∶3∶4,調(diào)整兩供肥裝置的外槽輪排肥器,使普通復(fù)合肥和緩釋肥的比例為3∶7,然后再調(diào)整肥料變比例分配器的調(diào)節(jié)裝置,將相鄰兩出肥口的排肥比例調(diào)整為3∶4,使排肥系統(tǒng)上、中、下3層的施肥比例為3∶3∶4。作業(yè)時,鏟尖入土,分層深施肥裝置隨著機(jī)具前進(jìn)到達(dá)預(yù)定作業(yè)深度,排肥系統(tǒng)開始工作,經(jīng)肥料均勻分配裝置分配后的普通復(fù)合肥進(jìn)入分層深施肥裝置的上層肥料排肥管,在Y形分配器的作用下,肥料被均勻分配至兩側(cè),上層肥料左右兩側(cè)翼鏟隨機(jī)具前進(jìn),將深松鏟鏟柄兩側(cè)的土壤橫向切開,肥料從排肥口排出,排出的肥料落在側(cè)翼劃開的溝槽底部,避免肥料落至下層,造成分層不明顯。經(jīng)肥料變比例分配裝置分配后的緩釋肥進(jìn)入分層深施肥裝置的中層和下層肥料排肥管,中層肥料在Y形分配器的作用下,被均勻分配至兩側(cè),經(jīng)排肥口排出,落入中層肥料側(cè)翼切開的溝槽底部。下層肥料經(jīng)排肥口排出,落入深松鏟剖開的溝槽底部。實現(xiàn)普通復(fù)合肥和緩釋肥分3層施用,且上層和中層肥料兩側(cè)分施,3層肥料形成半包圍種子的分布狀態(tài)。
分層深施肥裝置中深松部件為直鏟式深松鏟,基本參數(shù)如圖4所示,鏟尖寬度S7為45 mm,鏟柄厚度為20 mm,鏟尖入土角θ為23°[31]。肥料分層施肥器排肥管安裝在深松鏟鏟柄后方,3個排肥管的直徑D為20 mm。為保證分層深施肥裝置在土壤中前進(jìn)時的穩(wěn)定性,使上層和中層肥料左右兩側(cè)翼鏟入土角γ與深松鏟鏟尖入土角θ相等,為23°,避免因入土角度不同,分層深施肥裝置受力不均衡,造成施肥深度不穩(wěn)定。上層和中層肥料左右兩側(cè)翼鏟的長度S3、S4相等,為110 mm,上層肥料翼鏟水平長度S5與上層肥料距種子的橫向距離S1相等,為120 mm,中層肥料翼鏟水平長度S6與中層肥料距種子的橫向距離S2相等,為80 mm,上層肥料翼鏟與中層肥料翼鏟的距離H5與上層肥料距離中層肥料的深度H3相等,為100 mm,中層肥料翼鏟與深松鏟鏟尖距離H6與中層肥料距離下層肥料的深度H4相等,為100 mm。上層和中層肥料翼鏟安裝在深松鏟鏟柄的兩側(cè),為保證肥料順利下落,翼鏟需向下傾斜角度β,該角度不僅影響肥料的順利排出,同時影響施肥裝置的回土,進(jìn)而影響分層深施肥作業(yè)效果,是分層深施肥裝置關(guān)鍵參數(shù)之一,該傾角的兩個極限位置是垂直于鏟柄和平行于鏟柄,因此其理論取值范圍為0°~90°,根據(jù)深松鏟柄長度和排肥口距種子的橫向距離S1、S2,結(jié)合肥料落下需有一定傾角,該傾角應(yīng)小于肥料堆積角的余角,因此將該取值范圍調(diào)整為50°~64°。翼鏟在前進(jìn)過程中需切開土壤,為降低翼鏟所受阻力,翼鏟向后傾斜一定角度α,使側(cè)翼傾斜切開土壤,該角度會影響施肥裝置的回土?xí)r間,進(jìn)而影響分層深施肥作業(yè)效果,是分層深施肥裝置關(guān)鍵參數(shù)之一,該傾角的理論取值范圍為0°~90°,根據(jù)鏟柄的寬度和排肥口距種子的橫向距離S1、S2,將該取值范圍調(diào)整為50°~90°。
2.2.1翼鏟切開土壤受力分析
分層施肥裝置入土后,上層和中層肥料出口對應(yīng)的兩側(cè)翼鏟,隨著機(jī)具的前進(jìn)切開土壤,翼鏟上部土壤隨翼鏟前進(jìn)被向上抬升,肥料從排肥口排出后落入切開的土壤中,翼鏟上部土壤自然回落,覆蓋肥料,將肥料固定在相應(yīng)深度。將土壤顆粒視為散粒體,分析翼鏟前端切開土壤過程,認(rèn)為分層深施肥鏟在施肥過程中運(yùn)動為勻速運(yùn)動,深松和施肥深度穩(wěn)定。如圖5所示,翼鏟前端土壤受到翼鏟的支持力FN以及翼鏟對土壤顆粒的摩擦力μFN,翼鏟向后傾斜角α在水平面的投影角度為α′,作業(yè)速度為v,土壤受到的合力為Fc,將土壤所受的力在水平方向和豎直方向分解,得到
(1)
整理得
(2)
式中μ——土壤顆粒與翼鏟的摩擦因數(shù)
其中,α的取值范圍為50°~90°,因此α在水平面的投影角α′小于90°,由式(2)可知,土壤顆粒受到的合力Fc隨α′的增大而增大,因此,側(cè)翼向后傾斜角α?xí)绊懸礴P對土壤的切削,進(jìn)而影響分層施肥作業(yè)效果。
2.2.2翼鏟上方土壤運(yùn)動分析
土壤被翼鏟切開后,翼鏟上方土壤在翼鏟的作用下與翼鏟產(chǎn)生相對運(yùn)動,以正視基準(zhǔn)面將翼鏟剖開,以翼鏟上方的土壤顆粒為研究對象,其運(yùn)動如圖6所示,該土壤顆粒的運(yùn)動微分方程為
(3)
(4)
其中
F′N=mgsinα″
(5)
整理得
(6)
對式(6)求積分,得
vx=-gtsinα″-μgtcosα″+C1
(7)
整理得
(8)
對式(8)求積分,得
x=-gt2cosα″+μgt2sinα″+C1t+C2
(9)
式中m——土壤顆粒質(zhì)量,kg
g——重力加速度,m/s2
α″——翼鏟向后傾斜角在正視基準(zhǔn)面的投影角,(°)
F′N——土壤受到翼鏟的支持力,N
vx——土壤顆粒沿x軸的速度,m/s
C1、C2——在該時刻下土壤顆粒的運(yùn)動常數(shù)
x——土壤沿x軸的位移,m
t——時間,s
由式(9)可知,土壤顆粒在翼鏟上方的運(yùn)動受到側(cè)翼向后傾斜角α在正視基準(zhǔn)面的投影角α″的影響,因此,側(cè)翼向后傾斜角α?xí)绊懲寥李w粒的運(yùn)動,進(jìn)而影響分層施肥作業(yè)效果。
2.2.3排肥管內(nèi)肥料顆粒受力分析
排肥系統(tǒng)中經(jīng)分配器分配后的肥料進(jìn)入分層深施肥裝置排肥管,上層和中層肥料經(jīng)Y形分配器被分配至兩側(cè),肥料顆粒在自身重力、摩擦力、管壁的支持力以及氣流的曳力作用下,運(yùn)動至排肥口。肥料顆粒在兩側(cè)排肥管內(nèi)的受力情況如圖7所示。
肥料顆粒在沿斜面方向和垂直于斜面方向的受力方程為
(10)
整理得
Fr=Fk-μ1m1gsinβ+m1gcosβ
(11)
式中Fr——肥料顆粒受到的合力,N
Fk——?dú)饬髯饔迷诜柿项w粒上的曳力,N
F1N——肥料顆粒受到管壁的支持力,N
m1——肥料顆粒質(zhì)量,kg
μ1——肥料顆粒與管壁的摩擦因數(shù)
其中,β的取值范圍為50°~64°,由式(11)可知,肥料顆粒受到的合力Fr隨著β的增大而減小,因此側(cè)翼向下傾斜角β會影響肥料顆粒的運(yùn)動,進(jìn)而影響分層施肥作業(yè)的效果。
分層深施肥裝置作業(yè)過程中,各層肥料深度處的回土速度是影響分層施肥作業(yè)效果的關(guān)鍵因素之一,分層深施肥裝置鏟尖和翼鏟切開土壤后,土壤依靠自身重力下落,完成回土作業(yè),施肥作業(yè)過程中隨著分層深施肥裝置的前進(jìn),被鏟尖和翼鏟抬升破碎的土壤在慣性的作用下相對于分層施肥裝置繼續(xù)向后運(yùn)動,當(dāng)作業(yè)速度增大時,破碎的土壤相對于分層施肥裝置運(yùn)動的距離越大,回土?xí)r間增大,會造成施肥后回土不及時,影響分層施肥效果,因此施肥作業(yè)速度v也是影響分層施肥作業(yè)效果的重要因素,參照文獻(xiàn)[31-32],深松作業(yè)速度范圍為2~4 km/h。
綜上,影響分層深施肥裝置作業(yè)效果的因素主要有:翼鏟向后傾斜角α、翼鏟向下傾斜角β和機(jī)具作業(yè)速度v。各因素的取值范圍:α為50°~90°,β為50°~64°,v為2~4 km/h。
分層深施肥裝置是氣力式排肥系統(tǒng)的關(guān)鍵部件,其分層深施肥作業(yè)效果是評價排肥系統(tǒng)性能的關(guān)鍵指標(biāo),本研究以各層實際施肥量與理論值的誤差來評價分層深施肥作業(yè)效果。采用計算流體力學(xué)和離散元方法分別建立分層深施肥裝置、土壤和肥料顆粒的仿真模型,進(jìn)行耦合仿真,分析分層深施肥裝置的翼鏟向后傾斜角α、翼鏟向下傾斜角β和機(jī)具作業(yè)速度v對分層深施肥裝置作業(yè)效果的影響。
目前計算流體力學(xué)和離散元法耦合仿真已廣泛應(yīng)用于農(nóng)業(yè)工程領(lǐng)域,該方法可準(zhǔn)確模擬物料顆粒在復(fù)雜受力情況下的運(yùn)動[33]。選用ANSYS Fluent 18.2和EDEM 2019軟件進(jìn)行耦合,模擬分析分層深施肥裝置的作業(yè)過程。
3.1.1仿真模型建立
利用三維建模軟件建立分層深施肥裝置三維模型,導(dǎo)入EDEM軟件中,設(shè)置仿真模型中鏟尖的材料為65Mn鋼。在EDEM軟件中建立土壤模型,選擇Hertz-Mindlin無粘接觸模型作為土壤顆粒間接觸模型,建立尺寸(長×寬×高)為2 000 mm×600 mm×400 mm的土槽,其中土壤模型由半徑5 mm的球形顆粒隨機(jī)排布組成。建立的分層深施肥裝置和土壤模型如圖8所示。
分層深施肥作業(yè)所用肥料為顆粒狀復(fù)合肥,其中普通復(fù)合肥選用中國-阿拉伯化肥有限公司生產(chǎn)的普通型復(fù)合肥(N、P2O5、K2O質(zhì)量比為15∶15∶15),緩釋肥選用山東金正大生態(tài)工程股份有限公司生產(chǎn)的緩釋型復(fù)合肥(N、P2O5、K2O質(zhì)量比為28∶8∶8),從兩種肥料中各隨機(jī)選取5組一定質(zhì)量的肥料顆粒,利用細(xì)鹽填充法測定兩種肥料的密度,取5組數(shù)據(jù)的平均值,測得普通復(fù)合肥的顆粒密度為1 476 kg/m3,緩釋肥顆粒密度為1 564 kg/m3,再從兩種肥料中各隨機(jī)選取50粒肥料顆粒測量三軸尺寸,測得普通復(fù)合肥平均三軸尺寸為3.686 mm×3.369 mm×3.076 mm,緩釋肥的平均三軸尺寸為3.731 mm×3.527 mm×3.384 mm,且兩種肥料的球形率均大于0.9,因此可將普通復(fù)合肥顆粒簡化為直徑為3.38 mm的球體,緩釋肥顆粒簡化為直徑為3.55 mm的球體。
經(jīng)預(yù)試驗和參照文獻(xiàn)[31,34],分層深施肥裝置、土壤和肥料顆粒的模型參數(shù)如表1所示。
表1 模型參數(shù)Tab.1 Model parameters
分層深施肥裝置、土壤和肥料顆粒之間的接觸參數(shù)如表2所示。
表2 接觸參數(shù)Tab.2 Contact parameters
3.1.2仿真參數(shù)設(shè)定
氣力式排肥系統(tǒng)的分配器將肥料分配后,氣流與肥料一同進(jìn)入分層深施肥裝置上、中、下3層肥料的排肥管,肥料在自身重力和氣流的作用下排出分層深施肥裝置,排肥管中氣流的運(yùn)動為不可壓縮流體的湍流運(yùn)動,遵循動量守恒和質(zhì)量守恒,利用軟件中的標(biāo)準(zhǔn)κ-ε模型瞬態(tài)求解的Eulerian-Eulerian耦合算法,計算氣流和肥料之間的動量和能量交換,描述氣流和肥料之間的相互影響關(guān)系。耦合時,兩個軟件中計算Time Step和Data Save的頻率應(yīng)成整數(shù)倍匹配,設(shè)置EDEM中Time Step為1.1×106s,F(xiàn)luent中Time Step為1.1×103s。經(jīng)前期預(yù)試驗,確定排肥管入口氣流速度為5.91 m/s[27],在此氣流速度下,可以滿足輔助肥料從排肥管內(nèi)順利排出,且在排肥口處不會對土壤產(chǎn)生擾動,影響分層作業(yè)效果。EDEM中設(shè)置上層排肥管排出普通復(fù)合肥,中層和下層排肥管排出緩釋肥,上、中、下3層排肥管排出肥料量質(zhì)量比為3∶3∶4,根據(jù)最小作業(yè)速度選擇總仿真時長為4 s。
為研究翼鏟的不同角度和作業(yè)速度對分層深施肥作業(yè)的影響,以分層深施肥裝置各層實際施肥量與理論施肥量誤差為試驗指標(biāo),以翼鏟向后傾斜角α、翼鏟向下傾斜角β和機(jī)具作業(yè)速度v為試驗因素,進(jìn)行二次正交旋轉(zhuǎn)組合仿真試驗。各因素的水平范圍為:翼鏟向后傾斜角50°~90°、翼鏟向下傾斜角50°~64°、作業(yè)速度2~4 km/h。玉米施肥作業(yè)施肥量一般為250~500 kg/hm2[35],為滿足施肥作業(yè)要求,選取最大施肥量500 kg/hm2進(jìn)行仿真試驗,玉米行距為600 mm,根據(jù)作業(yè)速度得到分層深施肥鏟施肥速率(g/s),約為8.33v,將施肥速率按3∶3∶4的比例分配到3根排肥管,上、中、下3層排肥管的施肥速率分別為:2.5v、2.5v、3.33v。
參考GB/T 9478—2005《谷物條播機(jī) 試驗方法》中相關(guān)試驗方法,在仿真試驗完成后的土槽中隨機(jī)選取10個長度為100 mm的區(qū)域,重復(fù)3次,通過EDEM后處理模塊,在選取區(qū)域中,添加質(zhì)量傳感器,分別得到每個區(qū)域內(nèi)土壤深度80~100 mm、180~200 mm、280~300 mm范圍內(nèi)肥料顆粒的質(zhì)量,即上、中、下3層的實際施肥量,如圖9所示,計算每層的實際施肥量與理論施肥量的偏差穩(wěn)定性系數(shù),評價分層深施肥作業(yè)效果。施肥后土壤的截面如圖10所示。
通過耦合仿真試驗,對影響分層深施肥作業(yè)效果的因素進(jìn)行顯著性分析,并根據(jù)實際需求對各參數(shù)組合進(jìn)行優(yōu)化。試驗因素編碼如表3所示。
表3 試驗因素編碼Tab.3 Experimental factors and codes
確定二次正交旋轉(zhuǎn)組合的試驗方案并進(jìn)行仿真試驗,試驗完成后計算上、中、下3層施肥量偏差穩(wěn)定性系數(shù),并進(jìn)行方差分析。試驗方案及結(jié)果如表4所示。
表4 試驗方案與結(jié)果Tab.4 Test plan and experimental data
3.3.1試驗結(jié)果分析與回歸模型建立
將試驗結(jié)果導(dǎo)入Design-Expert中,進(jìn)行二次回歸分析,并進(jìn)行多元回歸擬合,得到上、中、下3層施肥量偏差穩(wěn)定性系數(shù)Y1、Y2、Y3與翼鏟后傾角、翼鏟下傾角、作業(yè)速度之間的數(shù)學(xué)關(guān)系模型,檢驗其顯著性,并對交互作用規(guī)律進(jìn)行分析。
(1)上層施肥量偏差穩(wěn)定性系數(shù)
通過對試驗數(shù)據(jù)的分析和擬合,上層施肥量偏差穩(wěn)定性系數(shù)Y1方差分析如表5所示。由表5可知,試驗數(shù)據(jù)整體模型極顯著(P<0.01),線性主效應(yīng)項中翼鏟后傾角A、翼鏟下傾角B和作業(yè)速度C對試驗指標(biāo)上層施肥量偏差穩(wěn)定性系數(shù)Y1影響極顯著(P<0.01),交互項中翼鏟后傾角和翼鏟下傾角的交互項AB對試驗指標(biāo)影響較顯著(0.05≤P<0.1),二次主效應(yīng)項中翼鏟下傾角B的二次項對試驗指標(biāo)影響較顯著(0.05≤P<0.1),其他各項不顯著,各因素對上層施肥量偏差穩(wěn)定性系數(shù)影響由大到小依次為A、B、C。將方差分析中不顯著的交互項和二次主效應(yīng)項的平方和、自由度并入殘差,剔除不顯著因素,再次進(jìn)行方差分析,結(jié)果如表5所示。得到各因素水平對上層施肥量偏差穩(wěn)定性系數(shù)Y1影響的回歸方程為
Y1=-115.004 31-1.185 56A+4.602 25B+
3.665 03C+0.025 001AB-0.050 123B2
(12)
對上述回歸方程進(jìn)行失擬項檢驗,如表5所示,得到失擬項P為0.574 9,不顯著(P>0.1),表明不存在影響試驗指標(biāo)的其他因素,且試驗指標(biāo)與試驗因素之間存在顯著的二次關(guān)系。
(2)中層施肥量偏差穩(wěn)定性系數(shù)
通過對試驗數(shù)據(jù)的分析和擬合,中層施肥量偏差穩(wěn)定性系數(shù)Y2方差分析如表5所示。由表5可知,試驗數(shù)據(jù)整體模型極顯著(P<0.01),線性主效應(yīng)項中翼鏟后傾角A、翼鏟下傾角B和作業(yè)速度C對試驗指標(biāo)中層施肥量偏差穩(wěn)定性系數(shù)Y2影響極顯著(P<0.01),二次主效應(yīng)項中翼鏟下傾角B的二次項對試驗指標(biāo)影響極顯著(P<0.01),交互項中翼鏟后傾角和翼鏟下傾角的交互項AB對試驗指標(biāo)影響較顯著(0.05≤P<0.1),其他各項不顯著,各因素對中層施肥量偏差穩(wěn)定性系數(shù)影響由大到小依次為B、A、C。將方差分析中不顯著的交互項和二次主效應(yīng)項的平方和、自由度并入殘差,剔除不顯著因素,再次進(jìn)行方差分析,結(jié)果如表5所示。得到各因素水平對中層施肥量偏差穩(wěn)定性系數(shù)Y2影響的回歸方程為
Y2=-182.310 93-0.908 53A+6.673 85B+
3.112 72C+0.019 496AB-0.064 965B2
(13)
對上述回歸方程進(jìn)行失擬項檢驗,如表5所示,得到失擬項P為0.686 7,不顯著(P>0.1),表明不存在影響試驗指標(biāo)的其他因素,且試驗指標(biāo)與試驗因素之間存在顯著的二次關(guān)系。
(3)下層施肥量偏差穩(wěn)定性系數(shù)
通過對試驗數(shù)據(jù)的分析和擬合,下層施肥量偏差穩(wěn)定性系數(shù)Y3方差分析如表5所示。由表5可知,試驗數(shù)據(jù)整體模型極顯著(P<0.01),線性主效應(yīng)項中翼鏟后傾角A、翼鏟下傾角B和作業(yè)速度C對試驗指標(biāo)上層施肥量偏差穩(wěn)定性系數(shù)Y3影響極顯著(P<0.01),交互項中翼鏟后傾角和翼鏟下傾角的交互項AB對試驗指標(biāo)影響較顯著(0.05≤P<0.1),二次主效應(yīng)項中翼鏟下傾角B的二次項對試驗指標(biāo)影響較顯著(0.05≤P<0.1),其他各項不顯著,各因素對下層施肥量偏差穩(wěn)定性系數(shù)影響由大到小依次為A、B、C。將方差分析中不顯著的交互項和二次主效應(yīng)項的平方和、自由度并入殘差,剔除不顯著因素,再次進(jìn)行方差分析,結(jié)果如表5所示。得到各因素水平對上層施肥量偏差穩(wěn)定性系數(shù)Y3影響的回歸方程為
Y3=-118.248 64-1.263 84A+4.752 45B+
3.928 90C+0.026 668AB-0.052 115B2
(14)
對上述回歸方程進(jìn)行失擬項檢驗,如表5所示,得到失擬項P為0.522 6,不顯著(P>0.1),表明不存在影響試驗指標(biāo)的其他因素,且試驗指標(biāo)與試驗因素之間存在顯著的二次關(guān)系。
表5 上、中、下3層施肥量偏差穩(wěn)定性系數(shù)方差分析Tab.5 Variance analysis of fertilization amount deviation stability coefficient of upper, middle and lower layers
3.3.2響應(yīng)曲面分析
對試驗結(jié)果進(jìn)行分析處理,根據(jù)建立的上、中、下3層施肥量偏差穩(wěn)定性系數(shù)回歸模型,各因素之間的交互作用對試驗指標(biāo)的影響效應(yīng)響應(yīng)曲面如圖11所示。
作業(yè)速度為3 km/h時,翼鏟后傾角A和翼鏟下傾角B對上層施肥量偏差穩(wěn)定性系數(shù)Y1的交互影響如圖11a所示。當(dāng)翼鏟后傾角一定時,翼鏟下傾角與上層施肥量偏差穩(wěn)定性系數(shù)呈正相關(guān),較優(yōu)的翼鏟下傾角范圍為52.84°~54.23°。當(dāng)翼鏟下傾角一定時,翼鏟后傾角與上層施肥量偏差穩(wěn)定性系數(shù)呈正相關(guān),較優(yōu)的翼鏟下傾角范圍為58.11°~66.04°。
作業(yè)速度為3 km/h時,翼鏟后傾角A和翼鏟下傾角B對中層施肥量偏差穩(wěn)定性系數(shù)Y2的交互影響如圖11b所示。當(dāng)翼鏟后傾角一定時,翼鏟下傾角與上層施肥量偏差穩(wěn)定性系數(shù)呈正相關(guān),較優(yōu)的翼鏟下傾角范圍為52.84°~57.00°。當(dāng)翼鏟下傾角一定時,翼鏟后傾角與上層施肥量偏差穩(wěn)定性系數(shù)呈正相關(guān),較優(yōu)的翼鏟下傾角范圍為58.11°~64.05°。
作業(yè)速度為3 km/h時,翼鏟后傾角A和翼鏟下傾角B對中層施肥量偏差穩(wěn)定性系數(shù)Y3的交互影響如圖11c所示。當(dāng)翼鏟后傾角一定時,翼鏟下傾角與上層施肥量偏差穩(wěn)定性系數(shù)呈正相關(guān),較優(yōu)的翼鏟下傾角范圍為52.84°~57.59°。當(dāng)翼鏟下傾角一定時,翼鏟后傾角與上層施肥量偏差穩(wěn)定性系數(shù)呈正相關(guān),較優(yōu)的翼鏟下傾角范圍為58.11°~75.95°。
為獲得分層深施肥裝置較優(yōu)的作業(yè)效果,利用Design-Expert軟件的優(yōu)化模塊對3個回歸模型進(jìn)行優(yōu)化求解,根據(jù)分層深施肥裝置的實際作業(yè)條件和工作要求,選擇目標(biāo)函數(shù)的約束條件。
目標(biāo)函數(shù)及約束條件為
(15)
根據(jù)約束條件,對目標(biāo)函數(shù)進(jìn)行優(yōu)化求解,其中作業(yè)速度以最大值為優(yōu)化目標(biāo),即在保證滿足作業(yè)效果的前提下,作業(yè)速度越大越好,求解得到多種優(yōu)化后的參數(shù)組合,結(jié)合實際分層深施肥裝置施肥作業(yè)的農(nóng)藝要求,從中選取較優(yōu)的參數(shù)組合:翼鏟后傾角為58.11°、翼鏟下傾角為52.84°、作業(yè)速度為3.38 km/h時,對應(yīng)的上、中、下3層施肥量偏差穩(wěn)定性系數(shù)分別為8.50%、6.54%和9.10%。
為驗證仿真試驗優(yōu)化得到的參數(shù)組合,根據(jù)優(yōu)化結(jié)果加工了分層深施肥裝置,并進(jìn)行田間試驗驗證。為便于實際加工,將優(yōu)化得到的參數(shù)進(jìn)行了取整,取翼鏟后傾角為58°,翼鏟下傾角為53°,分層深施肥裝置如圖12a所示,試驗樣機(jī)如圖12b所示。
試驗所需普通復(fù)合肥選用中國-阿拉伯化肥有限公司生產(chǎn)的普通型復(fù)合肥(N、P2O5、K2O質(zhì)量比為15∶15∶15),緩釋肥選用山東金正大生態(tài)工程股份有限公司生產(chǎn)的緩釋型復(fù)合肥(N、P2O5、K2O質(zhì)量比為28∶8∶8)。試驗前對試驗地土壤的基本物理特性參數(shù)進(jìn)行了測量,測得含水率為7.93%、容重為1.33 g/cm3。施肥前對排肥系統(tǒng)施肥量進(jìn)行標(biāo)定,使施肥量約為500 kg/hm2。啟動氣力式排肥系統(tǒng)各部件,進(jìn)行分層深施肥作業(yè),控制機(jī)具作業(yè)速度約為3.4 km/h,試驗完成后,根據(jù)GB/T 9478—2005《谷物條播機(jī) 試驗方法》中相關(guān)試驗方法,在施肥的行內(nèi),隨機(jī)選取30個點(diǎn),每個點(diǎn)截取100 mm長度區(qū)域,將深度80~100 mm、180~200 mm、280~300 mm范圍的土壤和肥料取出,將肥料和土壤篩分,稱量肥料的質(zhì)量,計算各層施肥量偏差穩(wěn)定性系數(shù)。圖13為施肥作業(yè)后肥料分層效果。
經(jīng)試驗測得上、中、下3層施肥量偏差穩(wěn)定性系數(shù)分別為9.52%、8.35%和9.87%,與仿真試驗得到的結(jié)果相比有一定差異,但差異不大,在1.81個百分點(diǎn)以內(nèi)。分析差異產(chǎn)生的原因主要是實際田間試驗中,土壤狀況相較于仿真試驗復(fù)雜多變,土壤中較大的沙礫等對試驗結(jié)果產(chǎn)生了影響,但作業(yè)效果仍能滿足設(shè)計要求。
(1)設(shè)計了一種氣力式排肥系統(tǒng)分層深施肥裝置,可將肥料按比例分層施入土壤中,其中上層和中層肥料兩側(cè)分施,3層肥料形成半包圍種子的分布狀態(tài),利于玉米對肥料養(yǎng)分的吸收。
(2)分層深施肥裝置實現(xiàn)了不同種類肥料的分層深施。上層施用普通復(fù)合肥,為玉米生長前期提供快速且充足養(yǎng)分,中層和下層施用緩釋肥,為玉米生長中后期提供長效穩(wěn)定的養(yǎng)分,避免因肥料利用不及時造成的養(yǎng)分流失和環(huán)境污染,同時減少作業(yè)環(huán)節(jié),降低對土壤的擾動,保護(hù)耕地土壤。
(3)通過理論分析確定了影響分層深施肥裝置施肥作業(yè)效果的主要因素,以各層施肥量偏差穩(wěn)定性系數(shù)為試驗指標(biāo),以翼鏟向后和向下傾斜角、作業(yè)速度為試驗因素,利用計算流體力學(xué)和離散元法耦合仿真的方法,進(jìn)行了二次正交旋轉(zhuǎn)組合試驗。通過對試驗結(jié)果的擬合和優(yōu)化分析,得到了翼鏟向后和向下傾斜角分別為58.11°和52.84°,作業(yè)速度為3.38 km/h時,上、中、下3層施肥量偏差穩(wěn)定性系數(shù)分別為8.50%、6.54%和9.10%,以優(yōu)化得到的參數(shù)進(jìn)行了裝置加工和田間試驗,試驗結(jié)果與仿真試驗優(yōu)化所得結(jié)果相吻合,滿足設(shè)計要求。