• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于氣力分配的不同肥料分層深施裝置設(shè)計與試驗

    2021-11-09 08:35:50楊慶璐李洪文盧彩云王英博王慶杰
    關(guān)鍵詞:作業(yè)

    楊慶璐 李洪文 何 進(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)

    0 引言

    東北黑土區(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ù)參考和理論支撐。

    1 氣力式排肥系統(tǒng)結(jié)構(gòu)與工作原理

    1.1 總體結(jié)構(gòu)

    氣力式排肥系統(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ù)合肥和緩釋肥分層施用。

    1.2 分層深施肥肥料分布

    分層深施肥作業(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所示。

    1.3 分層深施肥裝置結(jié)構(gòu)及工作原理

    分層深施肥裝置是實現(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)。

    2 分層深施肥裝置設(shè)計

    2.1 基本參數(shù)

    分層深施肥裝置中深松部件為直鏟式深松鏟,基本參數(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 作業(yè)過程分析

    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。

    3 分層深施肥裝置仿真試驗

    分層深施肥裝置是氣力式排肥系統(tǒng)的關(guān)鍵部件,其分層深施肥作業(yè)效果是評價排肥系統(tǒng)性能的關(guān)鍵指標(biāo),本研究以各層實際施肥量與理論值的誤差來評價分層深施肥作業(yè)效果。采用計算流體力學(xué)和離散元方法分別建立分層深施肥裝置、土壤和肥料顆粒的仿真模型,進(jìn)行耦合仿真,分析分層深施肥裝置的翼鏟向后傾斜角α、翼鏟向下傾斜角β和機(jī)具作業(yè)速度v對分層深施肥裝置作業(yè)效果的影響。

    3.1 DEM-CFD耦合仿真模型與參數(shù)設(shè)定

    目前計算流體力學(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。

    3.2 仿真試驗設(shè)計

    為研究翼鏟的不同角度和作業(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 仿真試驗結(jié)果與分析

    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%。

    4 試驗驗證

    4.1 試驗裝置

    為驗證仿真試驗優(yōu)化得到的參數(shù)組合,根據(jù)優(yōu)化結(jié)果加工了分層深施肥裝置,并進(jìn)行田間試驗驗證。為便于實際加工,將優(yōu)化得到的參數(shù)進(jìn)行了取整,取翼鏟后傾角為58°,翼鏟下傾角為53°,分層深施肥裝置如圖12a所示,試驗樣機(jī)如圖12b所示。

    4.2 試驗方法與結(jié)果

    試驗所需普通復(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è)計要求。

    5 結(jié)論

    (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è)計要求。

    猜你喜歡
    作業(yè)
    作業(yè),我終于打敗你了!
    小主人報(2022年1期)2022-08-10 08:28:44
    讓人羨慕嫉妒恨的“作業(yè)人”
    作業(yè)聯(lián)盟
    我愿作業(yè)少一點(diǎn)
    快來寫作業(yè)
    一次特殊的作業(yè)
    誰沒交作業(yè)
    修改“作業(yè)”
    跟一群抄作業(yè)的講垂直進(jìn)步?
    能源(2016年2期)2016-12-01 05:10:46
    作業(yè)
    故事大王(2016年7期)2016-09-22 17:30:08
    久久久a久久爽久久v久久| 成人国产av品久久久| 国产欧美日韩一区二区三区在线 | 尤物成人国产欧美一区二区三区| 亚洲欧美一区二区三区国产| 精品人妻熟女av久视频| 国产人妻一区二区三区在| 精品久久久噜噜| 成人美女网站在线观看视频| 在线观看人妻少妇| 老师上课跳d突然被开到最大视频| 九九爱精品视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产亚洲av天美| 欧美成人精品欧美一级黄| 亚洲天堂国产精品一区在线| 亚洲人成网站高清观看| 欧美潮喷喷水| 国产精品99久久久久久久久| 三级经典国产精品| 夫妻性生交免费视频一级片| av网站免费在线观看视频| 久久国产乱子免费精品| 亚洲av.av天堂| 97超视频在线观看视频| 国产成年人精品一区二区| 欧美精品国产亚洲| 国产男人的电影天堂91| 国国产精品蜜臀av免费| 久久人人爽av亚洲精品天堂 | 国产精品久久久久久av不卡| 嫩草影院新地址| 欧美成人一区二区免费高清观看| 国内精品宾馆在线| 国产美女午夜福利| 免费黄网站久久成人精品| 午夜精品一区二区三区免费看| 亚洲性久久影院| 神马国产精品三级电影在线观看| 国产视频首页在线观看| 免费看av在线观看网站| 亚洲最大成人中文| 少妇人妻 视频| 69av精品久久久久久| 晚上一个人看的免费电影| 亚洲,欧美,日韩| av在线观看视频网站免费| 丰满乱子伦码专区| 99热全是精品| 亚洲在久久综合| 久久国产乱子免费精品| 91精品伊人久久大香线蕉| 亚洲精品国产av蜜桃| 黄色一级大片看看| 美女国产视频在线观看| 亚洲欧美清纯卡通| 亚洲av不卡在线观看| 国产91av在线免费观看| 欧美性感艳星| 97热精品久久久久久| 人妻系列 视频| 免费少妇av软件| 大香蕉久久网| 亚洲国产高清在线一区二区三| 午夜福利高清视频| 免费观看a级毛片全部| 亚洲一级一片aⅴ在线观看| 亚洲国产欧美人成| 一个人观看的视频www高清免费观看| 最近2019中文字幕mv第一页| 又爽又黄a免费视频| 高清视频免费观看一区二区| 亚洲成人av在线免费| 人妻夜夜爽99麻豆av| 97人妻精品一区二区三区麻豆| 一本色道久久久久久精品综合| 精品久久国产蜜桃| 一区二区三区四区激情视频| 毛片一级片免费看久久久久| 国产美女午夜福利| 伊人久久国产一区二区| 免费黄频网站在线观看国产| 欧美高清成人免费视频www| 久久久久国产网址| 天堂中文最新版在线下载 | 国产成人免费观看mmmm| 国产亚洲最大av| 久久精品久久精品一区二区三区| 一区二区三区乱码不卡18| 秋霞在线观看毛片| 国产爽快片一区二区三区| 老女人水多毛片| 特级一级黄色大片| 91久久精品电影网| 日韩精品有码人妻一区| 国产高清有码在线观看视频| 狠狠精品人妻久久久久久综合| 白带黄色成豆腐渣| 2018国产大陆天天弄谢| 乱系列少妇在线播放| 波多野结衣巨乳人妻| 免费看a级黄色片| 久久久久久久精品精品| 丝瓜视频免费看黄片| 亚洲经典国产精华液单| 久久99热这里只有精品18| 亚洲婷婷狠狠爱综合网| 亚洲人成网站高清观看| 日本与韩国留学比较| 久久久久久久精品精品| av在线app专区| 精品熟女少妇av免费看| 欧美另类一区| 亚洲精品国产av蜜桃| 欧美区成人在线视频| av又黄又爽大尺度在线免费看| 精品一区二区三区视频在线| 亚洲国产精品999| 免费观看a级毛片全部| 欧美成人a在线观看| 国产日韩欧美亚洲二区| 国产永久视频网站| av在线天堂中文字幕| 国产v大片淫在线免费观看| 黄色日韩在线| 亚洲精品日韩在线中文字幕| 一区二区av电影网| 国产成人精品婷婷| 一区二区三区乱码不卡18| 极品教师在线视频| 97超碰精品成人国产| 交换朋友夫妻互换小说| 一本色道久久久久久精品综合| 久久久精品欧美日韩精品| 国产片特级美女逼逼视频| 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 国产精品国产三级专区第一集| 亚洲av一区综合| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 日本猛色少妇xxxxx猛交久久| 亚洲精品成人av观看孕妇| 精品少妇黑人巨大在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美成人精品一区二区| 美女国产视频在线观看| 人妻系列 视频| 白带黄色成豆腐渣| 久久久久久伊人网av| 亚洲色图综合在线观看| 中文在线观看免费www的网站| 国产女主播在线喷水免费视频网站| 亚洲国产欧美在线一区| 国产精品偷伦视频观看了| 国产黄片美女视频| 国产精品一区二区性色av| 亚洲一级一片aⅴ在线观看| 天天躁日日操中文字幕| 国产乱来视频区| 高清毛片免费看| 国产欧美日韩精品一区二区| 老司机影院毛片| 国产一区二区在线观看日韩| 亚洲精品国产色婷婷电影| 日韩免费高清中文字幕av| 超碰97精品在线观看| 成人午夜精彩视频在线观看| .国产精品久久| eeuss影院久久| 欧美性猛交╳xxx乱大交人| 日韩 亚洲 欧美在线| 中文在线观看免费www的网站| .国产精品久久| 久久精品国产自在天天线| 色哟哟·www| 嘟嘟电影网在线观看| 夜夜看夜夜爽夜夜摸| 日韩强制内射视频| 在线观看人妻少妇| 亚洲成人中文字幕在线播放| 精品人妻熟女av久视频| 欧美区成人在线视频| 国产亚洲最大av| 99热这里只有精品一区| 日韩一本色道免费dvd| 久久精品久久精品一区二区三区| 成人午夜精彩视频在线观看| av免费在线看不卡| 成年免费大片在线观看| 我的老师免费观看完整版| 欧美日韩视频高清一区二区三区二| videos熟女内射| 日日啪夜夜撸| 国产免费视频播放在线视频| 国产欧美日韩精品一区二区| av免费在线看不卡| 香蕉精品网在线| 欧美xxⅹ黑人| 欧美成人a在线观看| 久热这里只有精品99| av国产免费在线观看| 久久久色成人| 日韩,欧美,国产一区二区三区| 国产精品人妻久久久影院| 欧美丝袜亚洲另类| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| videos熟女内射| 蜜臀久久99精品久久宅男| 亚洲婷婷狠狠爱综合网| 99re6热这里在线精品视频| 免费电影在线观看免费观看| 国产精品伦人一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲aⅴ乱码一区二区在线播放| 国产探花在线观看一区二区| 九九久久精品国产亚洲av麻豆| 亚洲欧美一区二区三区国产| 免费观看在线日韩| 69av精品久久久久久| 日韩一区二区视频免费看| eeuss影院久久| 亚洲欧美中文字幕日韩二区| 九九在线视频观看精品| 久久99热这里只频精品6学生| 日韩视频在线欧美| 美女主播在线视频| 久久精品夜色国产| av福利片在线观看| 久久97久久精品| 久久久精品免费免费高清| 男女下面进入的视频免费午夜| 成年人午夜在线观看视频| 国产欧美亚洲国产| 国产黄a三级三级三级人| 亚洲不卡免费看| 日日啪夜夜撸| 熟女人妻精品中文字幕| 99视频精品全部免费 在线| 嫩草影院精品99| 在线观看人妻少妇| 国产 一区精品| 一本一本综合久久| 亚洲成色77777| 久久影院123| 亚洲av成人精品一区久久| 深夜a级毛片| 在现免费观看毛片| 国产熟女欧美一区二区| 国产免费视频播放在线视频| 一个人看视频在线观看www免费| 99久久中文字幕三级久久日本| 国产黄片视频在线免费观看| 3wmmmm亚洲av在线观看| 国产高清国产精品国产三级 | 夫妻午夜视频| 毛片女人毛片| 大又大粗又爽又黄少妇毛片口| 纵有疾风起免费观看全集完整版| 国产男女内射视频| 久久久久国产精品人妻一区二区| 亚洲欧美成人精品一区二区| 亚洲自拍偷在线| 啦啦啦中文免费视频观看日本| 蜜桃久久精品国产亚洲av| 极品少妇高潮喷水抽搐| 亚洲成色77777| 亚洲av男天堂| 亚洲av免费高清在线观看| 久久鲁丝午夜福利片| 国产一区二区亚洲精品在线观看| 97在线人人人人妻| 久久精品久久久久久噜噜老黄| 97精品久久久久久久久久精品| 熟女电影av网| 97热精品久久久久久| 一个人看视频在线观看www免费| 夫妻午夜视频| 一级爰片在线观看| 久久久久精品久久久久真实原创| 极品教师在线视频| 成人黄色视频免费在线看| 狂野欧美激情性bbbbbb| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利视频1000在线观看| 精品酒店卫生间| 国产精品女同一区二区软件| 中文资源天堂在线| 亚洲精品国产成人久久av| 免费少妇av软件| 菩萨蛮人人尽说江南好唐韦庄| 超碰97精品在线观看| 久久精品人妻少妇| 亚洲人与动物交配视频| 午夜福利高清视频| 丝袜脚勾引网站| 免费观看av网站的网址| av国产久精品久网站免费入址| 国产精品99久久久久久久久| 欧美日韩视频高清一区二区三区二| 午夜视频国产福利| 成人特级av手机在线观看| 黄色欧美视频在线观看| 我的女老师完整版在线观看| 久久久久久久久久久丰满| 成人黄色视频免费在线看| 在线免费观看不下载黄p国产| 天堂中文最新版在线下载 | 搞女人的毛片| 黄色怎么调成土黄色| 80岁老熟妇乱子伦牲交| 精品久久国产蜜桃| 99久久精品一区二区三区| 亚洲精品色激情综合| 亚洲av在线观看美女高潮| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 国产精品无大码| 亚洲国产精品成人久久小说| 免费大片黄手机在线观看| 精品午夜福利在线看| 久久精品国产a三级三级三级| 我的老师免费观看完整版| 免费观看性生交大片5| 亚洲精品一二三| 高清av免费在线| 久久久久久久亚洲中文字幕| 日日撸夜夜添| av在线老鸭窝| 亚洲精品色激情综合| 中文在线观看免费www的网站| 一区二区av电影网| 欧美极品一区二区三区四区| 成人国产麻豆网| 国产精品久久久久久久久免| 亚洲精品成人久久久久久| 三级国产精品片| 日韩av在线免费看完整版不卡| 在线免费十八禁| 国内精品美女久久久久久| 最后的刺客免费高清国语| 久久女婷五月综合色啪小说 | 色5月婷婷丁香| 久久久精品94久久精品| 亚洲人成网站在线观看播放| 精品一区在线观看国产| 午夜免费男女啪啪视频观看| 久久影院123| 欧美高清成人免费视频www| 一级毛片电影观看| 日韩在线高清观看一区二区三区| 国产精品精品国产色婷婷| 亚洲精华国产精华液的使用体验| 国产精品精品国产色婷婷| 亚洲精品一二三| 国产精品爽爽va在线观看网站| 久热这里只有精品99| 有码 亚洲区| 成人欧美大片| 精品一区在线观看国产| 国产午夜精品一二区理论片| 国产视频首页在线观看| 国产毛片在线视频| 午夜免费鲁丝| 少妇高潮的动态图| 91在线精品国自产拍蜜月| 国产精品99久久久久久久久| 国产av国产精品国产| 人妻制服诱惑在线中文字幕| 亚洲精品影视一区二区三区av| 大香蕉久久网| 成人无遮挡网站| 国产男女超爽视频在线观看| 男插女下体视频免费在线播放| 国产有黄有色有爽视频| 亚洲精品日本国产第一区| 成人毛片60女人毛片免费| 久久久久网色| 午夜日本视频在线| 女人被狂操c到高潮| 国产 一区精品| 欧美 日韩 精品 国产| 日韩中字成人| 亚洲国产精品国产精品| 国产精品国产三级国产av玫瑰| 真实男女啪啪啪动态图| 欧美精品国产亚洲| 看黄色毛片网站| 亚洲欧洲日产国产| 午夜精品国产一区二区电影 | 国产高清国产精品国产三级 | 国产成人午夜福利电影在线观看| 国产 一区精品| av国产精品久久久久影院| 少妇猛男粗大的猛烈进出视频 | 亚洲精品成人久久久久久| 99久国产av精品国产电影| 特级一级黄色大片| 精品久久久久久久久av| 直男gayav资源| 欧美日韩在线观看h| 麻豆成人午夜福利视频| 色网站视频免费| 亚洲最大成人av| 男女那种视频在线观看| freevideosex欧美| 能在线免费看毛片的网站| 精品久久久久久久人妻蜜臀av| 亚洲欧美精品专区久久| 一本久久精品| 黄色配什么色好看| 成人无遮挡网站| av.在线天堂| 亚洲av免费高清在线观看| 中文天堂在线官网| 成年免费大片在线观看| 男女下面进入的视频免费午夜| 亚洲精品一区蜜桃| 亚州av有码| 国产色婷婷99| 精品久久久久久电影网| 干丝袜人妻中文字幕| 精品亚洲乱码少妇综合久久| 国产日韩欧美亚洲二区| 午夜福利高清视频| 国产精品久久久久久av不卡| 国产伦精品一区二区三区四那| 91久久精品国产一区二区三区| 亚洲精华国产精华液的使用体验| 91在线精品国自产拍蜜月| 久久久久久久久久久免费av| 亚洲va在线va天堂va国产| 天堂中文最新版在线下载 | 欧美潮喷喷水| 青青草视频在线视频观看| 国产毛片a区久久久久| 成人无遮挡网站| 大话2 男鬼变身卡| a级毛色黄片| 欧美日韩精品成人综合77777| 国产爽快片一区二区三区| 免费电影在线观看免费观看| 国产在线男女| 国产亚洲av嫩草精品影院| 中文天堂在线官网| 又爽又黄无遮挡网站| 亚洲av成人精品一二三区| 久久精品国产a三级三级三级| 爱豆传媒免费全集在线观看| 午夜免费观看性视频| 欧美少妇被猛烈插入视频| 伦精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 色网站视频免费| 中文乱码字字幕精品一区二区三区| 神马国产精品三级电影在线观看| 在线亚洲精品国产二区图片欧美 | 一级毛片aaaaaa免费看小| 青春草亚洲视频在线观看| 大香蕉久久网| 国产精品嫩草影院av在线观看| 国产乱人视频| 国产精品蜜桃在线观看| 七月丁香在线播放| 亚洲人成网站在线观看播放| 69av精品久久久久久| 日本免费在线观看一区| 免费黄色在线免费观看| 最新中文字幕久久久久| 成人综合一区亚洲| 综合色丁香网| 在线天堂最新版资源| 亚洲人与动物交配视频| 26uuu在线亚洲综合色| 涩涩av久久男人的天堂| 精品一区在线观看国产| 成人免费观看视频高清| av女优亚洲男人天堂| 亚洲成人精品中文字幕电影| 涩涩av久久男人的天堂| 国产视频首页在线观看| 成人亚洲精品一区在线观看 | 啦啦啦啦在线视频资源| 美女国产视频在线观看| 小蜜桃在线观看免费完整版高清| 一二三四中文在线观看免费高清| 亚洲av日韩在线播放| 麻豆乱淫一区二区| 成人国产麻豆网| 青春草视频在线免费观看| 亚洲av电影在线观看一区二区三区 | 777米奇影视久久| 日日摸夜夜添夜夜添av毛片| 大片免费播放器 马上看| 亚洲精品aⅴ在线观看| 亚洲精品中文字幕在线视频 | 少妇高潮的动态图| 色播亚洲综合网| 久久久久久久久久久免费av| 人妻少妇偷人精品九色| 亚洲欧美精品专区久久| 国产精品久久久久久久久免| 久久精品国产自在天天线| 久久久久久久亚洲中文字幕| 免费大片18禁| av.在线天堂| 三级国产精品片| 精品久久久久久电影网| 国产精品国产av在线观看| 久久韩国三级中文字幕| 国内少妇人妻偷人精品xxx网站| 如何舔出高潮| 又粗又硬又长又爽又黄的视频| 国产成人一区二区在线| 一级毛片我不卡| 亚洲精品乱码久久久久久按摩| 午夜精品国产一区二区电影 | 中文在线观看免费www的网站| 精品一区在线观看国产| 夫妻性生交免费视频一级片| 少妇人妻精品综合一区二区| 日本wwww免费看| 80岁老熟妇乱子伦牲交| 亚洲高清免费不卡视频| 亚洲一区二区三区欧美精品 | 久久久久久国产a免费观看| 伊人久久精品亚洲午夜| 在线 av 中文字幕| 3wmmmm亚洲av在线观看| 成人一区二区视频在线观看| 欧美zozozo另类| 久久久久久九九精品二区国产| 成人亚洲精品av一区二区| 国产黄a三级三级三级人| 国产成人精品福利久久| 国产成人aa在线观看| 精品国产三级普通话版| 午夜视频国产福利| 欧美xxxx黑人xx丫x性爽| 午夜日本视频在线| 新久久久久国产一级毛片| 激情 狠狠 欧美| 插逼视频在线观看| 国产精品一二三区在线看| 中国三级夫妇交换| 麻豆成人av视频| 国产亚洲av嫩草精品影院| 国产一区二区三区av在线| 午夜精品一区二区三区免费看| 免费观看性生交大片5| 永久免费av网站大全| 中文字幕免费在线视频6| 国产一区亚洲一区在线观看| 伊人久久精品亚洲午夜| 下体分泌物呈黄色| 日韩不卡一区二区三区视频在线| 国产精品国产av在线观看| 午夜老司机福利剧场| 黄片无遮挡物在线观看| 国产一区有黄有色的免费视频| 日本av手机在线免费观看| 在线看a的网站| 午夜精品国产一区二区电影 | 国产成人a区在线观看| 在线观看美女被高潮喷水网站| 寂寞人妻少妇视频99o| 欧美丝袜亚洲另类| 97超视频在线观看视频| 亚洲国产日韩一区二区| 黄色视频在线播放观看不卡| 少妇人妻 视频| 亚洲欧美成人精品一区二区| 久久国内精品自在自线图片| 日本色播在线视频| 青春草视频在线免费观看| 午夜福利视频精品| 国产极品天堂在线| 国产大屁股一区二区在线视频| 日韩成人av中文字幕在线观看| 两个人的视频大全免费| 99久久人妻综合| 99九九线精品视频在线观看视频| 国产免费福利视频在线观看| 精品99又大又爽又粗少妇毛片| 精品一区在线观看国产| 肉色欧美久久久久久久蜜桃 | 精品午夜福利在线看| 少妇高潮的动态图| 日韩一区二区三区影片| 丰满人妻一区二区三区视频av| 亚洲最大成人中文| 51国产日韩欧美| 最近最新中文字幕免费大全7| 免费看光身美女| 国产精品99久久99久久久不卡 | kizo精华| 一本一本综合久久| 搡女人真爽免费视频火全软件| 免费av观看视频| 日本欧美国产在线视频| 又大又黄又爽视频免费| 熟女av电影| 精品一区二区三区视频在线| 99久久精品热视频| 在线精品无人区一区二区三 | 免费看不卡的av| 亚洲欧美清纯卡通| 少妇丰满av| 中国美白少妇内射xxxbb| 成年人午夜在线观看视频| 精品久久久久久电影网| 一个人看视频在线观看www免费| 日韩 亚洲 欧美在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 |