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

    聯(lián)合整地機(jī)勻土旋平刀輥設(shè)計(jì)與試驗(yàn)

    2023-01-16 09:35:06劉尚坤張秀花劉江濤弋景剛張晉國(guó)
    關(guān)鍵詞:刀輥刀軸壤土

    劉尚坤,劉 超,張秀花,劉江濤,弋景剛,張晉國(guó)

    ·農(nóng)業(yè)裝備工程與機(jī)械化·

    聯(lián)合整地機(jī)勻土旋平刀輥設(shè)計(jì)與試驗(yàn)

    劉尚坤,劉 超,張秀花,劉江濤,弋景剛,張晉國(guó)

    (1. 河北農(nóng)業(yè)大學(xué)機(jī)電工程學(xué)院,保定 071001;2. 河北省智慧農(nóng)業(yè)裝備技術(shù)創(chuàng)新中心,保定 071001)

    為解決聯(lián)合整地機(jī)破茬翻埋作業(yè)后地表平整度低、土壤分布不均勻的問題,該研究設(shè)計(jì)了分段反螺旋式勻土旋平刀輥,提出田口法結(jié)合離散元仿真的刀輥參數(shù)設(shè)計(jì)思路。首先定性分析并確定了影響刀輥軸向運(yùn)土和拋土的關(guān)鍵因素及水平;然后以作業(yè)后地表平整度為優(yōu)化目標(biāo),以沙壤土和黏壤土為噪聲因子,利用田口法設(shè)計(jì)試驗(yàn)并進(jìn)行離散元仿真,得出適應(yīng)不同土壤類型的最優(yōu)參數(shù)組合為:刀軸轉(zhuǎn)速300 r/min、螺旋升角70°、旋耕刀型號(hào)IT245、刀軸直徑80 mm、刀座間距62 mm、旋耕刀雙螺旋排布,此時(shí)地表平整度仿真值為14.2 mm。以最優(yōu)參數(shù)組合制造樣機(jī)并進(jìn)行田間驗(yàn)證試驗(yàn),作業(yè)后地表平整度為11.6 mm、土壤分布均勻度為92.6%、耕深穩(wěn)定性系數(shù)為93.2%、破土率為86.4%,優(yōu)于小麥種植對(duì)耕整地的要求,所設(shè)計(jì)的刀輥有效。研究結(jié)果可為聯(lián)合整地機(jī)刀輥設(shè)計(jì)提供參考。

    離散元;田口法;聯(lián)合整地機(jī);分段反螺旋刀輥;地表平整度

    0 引 言

    中國(guó)華北地區(qū)以小麥-玉米輪作種植方式為主,小麥種植土地需滿足深、細(xì)、透、實(shí)、平的整地要求[1-2]。聯(lián)合整地機(jī)通過破茬、深松、翻埋實(shí)現(xiàn)秸稈還田、打破犁底層、改變土壤通透性等功能[3-5],作業(yè)地表平整度是整地機(jī)的一個(gè)關(guān)鍵技術(shù)指標(biāo)[6],對(duì)保障小麥種植后續(xù)農(nóng)藝過程及經(jīng)濟(jì)效益具有重要意義。

    王金武等[7]利用對(duì)稱式球面圓盤合墑裝置和改進(jìn)的柵條式碎土裝置聯(lián)合作業(yè)提高聯(lián)合整地機(jī)作業(yè)后地表平整度,達(dá)到耕層上實(shí)下虛的要求;鄭侃等[8]設(shè)計(jì)了一款適用于開溝旋耕機(jī)的漸變螺旋升角勻土刀輥,提高了地表平整度和土壤軸向分布均勻度;周華等[9]設(shè)計(jì)了一款秸稈還田深松旋埋聯(lián)合耕整機(jī),通過優(yōu)化刀輥上刀具排布提升了地表平整度和其他性能;Zheng等[10]根據(jù)深松旋耕聯(lián)合作業(yè)機(jī)深松后土壤的分布特征,設(shè)計(jì)了一款長(zhǎng)短刀組合式旋耕刀輥,提高了作業(yè)地表平整度;Yu等[11]設(shè)計(jì)了一款適用于秸稈還田中耕機(jī)的雙螺旋刀輥,利用EDEM離散元仿真優(yōu)化結(jié)構(gòu)參數(shù),達(dá)到提高地表平整度目的。以上研究均沒有考慮土壤類型變化對(duì)作業(yè)效果的影響,優(yōu)化結(jié)果局限于固定土壤類型。但實(shí)際作業(yè)中,不同地塊或同一地塊的不同區(qū)域會(huì)存在土壤物理特性差異,因此,在提升地表平整度時(shí)需要考慮土壤類型對(duì)作業(yè)效果的影響。

    田口法是由田口玄一在響應(yīng)面試驗(yàn)和正交試驗(yàn)的基礎(chǔ)上提出的一種試驗(yàn)方法[12],通過增添噪聲因子設(shè)計(jì)試驗(yàn),具有提升產(chǎn)品抵御外界干擾能力的優(yōu)點(diǎn),已廣泛應(yīng)用于汽車、化工、機(jī)加工等行業(yè)[13-14]。

    本文針對(duì)聯(lián)合整地機(jī)破茬翻埋作業(yè)后地表平整度低和土壤分布不均勻的問題,設(shè)計(jì)了分段反螺旋式勻土旋平刀輥,并將田口法引入整地機(jī)勻土旋平刀輥參數(shù)優(yōu)化中,以土壤類型為噪聲因子,結(jié)合離散元仿真優(yōu)化勻土旋平刀輥工作參數(shù),根據(jù)優(yōu)化參數(shù)制造物理樣機(jī)并進(jìn)行田間試驗(yàn),驗(yàn)證分段反螺旋式勻土旋平刀輥設(shè)計(jì)的有效性和田口法試驗(yàn)設(shè)計(jì)的穩(wěn)健性。

    1 聯(lián)合整地機(jī)結(jié)構(gòu)方案與工作原理

    1.1 聯(lián)合整地機(jī)結(jié)構(gòu)方案

    聯(lián)合整地機(jī)結(jié)構(gòu)由破茬裝置、翻埋裝置、勻土旋平裝置和鎮(zhèn)壓裝置組成,其中翻埋裝置包括機(jī)具中心線上的雙翼開溝犁和相對(duì)中心線對(duì)稱安裝的4個(gè)單翼鏵式犁,結(jié)構(gòu)方案如圖1所示,整機(jī)尺寸為3 730 mm′2 538 mm′1 340 mm,工作幅寬2.1 m,破茬刀輥?zhàn)鳂I(yè)深度(1)12 cm,勻土旋平刀輥?zhàn)鳂I(yè)深度(2)15 cm,翻埋作業(yè)深度(3)30 cm。

    注:h1為破茬刀輥?zhàn)鳂I(yè)深度,mm;h2為勻土旋平刀輥?zhàn)鳂I(yè)深度,mm;h3為翻埋作業(yè)深度,mm。

    1.2 工作原理

    工作時(shí),聯(lián)合整地機(jī)與拖拉機(jī)通過三點(diǎn)懸掛連接,由拖拉機(jī)提供動(dòng)力,破茬裝置將玉米根茬打碎,然后翻埋裝置通過交錯(cuò)式雙翼開溝犁、單翼鏵式犁打破犁底層,將破碎后的玉米根茬和地表殘留秸稈埋入土壤,隨后勻土旋平裝置通過勻土旋平刀輥碎土并把土壤壅入墑溝,最后利用鎮(zhèn)壓裝置壓實(shí)表層土壤,達(dá)到單次作業(yè)即可滿足破茬翻埋、土壤細(xì)碎和地表平整、上實(shí)下虛的小麥播種整地要求。

    2 勻土旋平刀輥設(shè)計(jì)與分析

    破茬作業(yè)后,雙翼犁將機(jī)具中心土壤翻向兩側(cè),與機(jī)具中軸線對(duì)稱布置的4個(gè)單翼犁將兩側(cè)土壤翻向中心,兩側(cè)邊產(chǎn)生墑溝,由于靠近機(jī)具中軸線的2個(gè)單翼犁不能將雙翼犁翻出的土壤完全回填,使中軸線上也產(chǎn)生較淺的墑溝,翻埋后地表呈M形分布。采用普通刀輥旋平及鎮(zhèn)壓后地表仍呈M形,不能滿足小麥種植平整度要求。

    2.1 勻土旋平刀輥設(shè)計(jì)

    為實(shí)現(xiàn)旋耕后土壤均勻分布及鎮(zhèn)壓平整,借鑒螺旋輸送原理,根據(jù)翻埋后地表M形分布中兩凸起位置間距及墑溝填土需求,本文設(shè)計(jì)了一種分段反螺旋式勻土旋平刀輥。該刀輥與刮土板、齒輪箱和機(jī)架組成勻土旋平裝置。刀輥的旋耕刀排列方式及螺旋線展開如圖2所示,左右2個(gè)刀輥與齒輪箱對(duì)稱布置,每側(cè)刀軸設(shè)計(jì)為I和II兩段,分別按2條螺旋線等距布置旋耕刀,其中I段在每個(gè)位置上朝外側(cè)徑向?qū)ΨQ布置2把旋耕刀,將部分土壤壅入兩側(cè)墑溝;II段在1′、2′、3′每個(gè)位置上朝內(nèi)側(cè)徑向?qū)ΨQ布置2把旋耕刀,將部分土壤壅入中間墑溝;每側(cè)刀座數(shù)量為15,每側(cè)旋耕刀總數(shù)量為32;2段螺旋線的升角方向相反,大小相等,即1=2;2段螺旋線相鄰處1與1′的起始點(diǎn)相位差為90°,以利于勻土及載荷平衡。

    注:k1、k2為螺旋線周數(shù);D為刀軸直徑,mm;H為導(dǎo)程,mm;l1、l2為螺旋升角,(°);Z為刀座間距,mm;13、12、11,……,1′、2′、3′為刀座位置編號(hào)。

    2.2 螺旋刀輥軸向運(yùn)土影響因素分析

    為保證合適的運(yùn)土量,必須明確螺旋刀輥參數(shù),確定刀輥軸向運(yùn)土的影響因素及水平,探求最優(yōu)參數(shù)組合。

    旋耕刀在刀軸上按螺旋線分布,螺旋線軌跡為:

    式中0為螺旋線上繞動(dòng)點(diǎn)角速度,rad/s;、和分別為螺旋線上繞動(dòng)點(diǎn)坐標(biāo)值,mm;為螺旋升角,(°)。

    螺旋升角是重要的勻土效果影響因素,當(dāng)轉(zhuǎn)速一定時(shí),增加,土壤輸送速度增大[15],但過大易產(chǎn)生較大振動(dòng),而過小則會(huì)引起刀間夾土,的合理范圍為54°~85°[16],本文試驗(yàn)取54°、70°、85°三個(gè)水平。

    由式(1)可知,當(dāng)螺旋升角一定時(shí),刀軸直徑與成正比,若直徑增大則導(dǎo)程增大、運(yùn)土效果增強(qiáng)[17-18],但是刀軸直徑過大會(huì)導(dǎo)致旋耕刀與機(jī)具干涉,直徑過小會(huì)降低機(jī)械強(qiáng)度,根據(jù)機(jī)架結(jié)構(gòu)參數(shù)和刀軸力學(xué)校核,刀軸直徑取70、80和90 mm三個(gè)水平。

    螺旋刀輥軸向運(yùn)土過程中,部分土壤會(huì)從相鄰兩旋耕刀的間隙被拋到機(jī)具后側(cè),所以相鄰旋耕刀的間距也是影響軸向運(yùn)土的關(guān)鍵因素,刀座間距與旋耕刀工作幅寬存在如下關(guān)系:

    式中為刀座間距,mm;為旋耕刀工作幅寬,mm;D為余量,查《農(nóng)業(yè)機(jī)械設(shè)計(jì)手冊(cè)》為15~20 mm。

    參照國(guó)標(biāo)GB/T 5669—2008《旋耕機(jī)械刀和刀座》,刀座軸向間距取50、60和70 mm三個(gè)水平。

    此外,旋耕刀雙螺旋布置相比于單螺旋布置運(yùn)土效率增大,因此刀輥的螺旋頭數(shù)也是影響刀輥運(yùn)土性能的因素,取單螺旋和雙螺旋2個(gè)水平。

    2.3 旋耕刀拋土影響因素分析

    勻土旋平作業(yè)時(shí),旋耕刀將前側(cè)土壤打碎拋到后側(cè),土壤的落地位置會(huì)影響平整效果,故對(duì)旋耕刀拋土過程的土壤顆粒建立受力分析模型,分析旋耕刀拋土影響因素。如圖3所示,分別以土壤顆粒中心位置和刀軸中心0為原點(diǎn)建立直角坐標(biāo)系,假設(shè)土壤顆粒與旋耕刀無(wú)相對(duì)位移、土壤顆粒間無(wú)摩擦力[19],土壤顆粒即將脫離旋耕刀的瞬時(shí)有如下受力關(guān)系:

    其中:

    =(6)

    式中為土壤顆粒質(zhì)量,kg;g為重力加速度,m/s2;為土壤顆粒與旋耕刀間摩擦系數(shù)。

    注:Fx為x方向上土壤顆粒所受合力,N;G為土壤顆粒重力,N;θ為旋耕刀折彎角,(°);F1為離心力在z0軸上的分力,N;fx為x方向的摩擦力,N;Fy為y方向上土壤顆粒所受合力,N;FC為離心力,N;α為顆粒脫離時(shí)刀軸轉(zhuǎn)角,(°);fy為y方向的摩擦力,N;Fz為z方向上土壤顆粒所受合力,N;FN為旋耕刀對(duì)土壤顆粒的支持力,N;ω為刀軸回轉(zhuǎn)角速度,rad·s-1;R1為土壤顆?;剞D(zhuǎn)半徑,m。

    土壤顆粒只沿著旋耕刀表面切向被拋出,此時(shí)F=0,結(jié)合式(3)~(9)可得土壤顆粒在和方向上的合力分別為

    土壤顆粒被拋出瞬時(shí)所受到的合力為

    式中F為土壤顆粒所受到的合力,N;為土壤顆粒離開旋耕刀時(shí)的瞬時(shí)速度,m/s;為土壤顆粒離開旋耕刀時(shí)的瞬時(shí)加速度,m/s2。

    由式(10)~(12)可知,土壤顆粒的拋出速度和加速度由摩擦系數(shù)、旋耕刀折彎角、回轉(zhuǎn)半徑1、刀軸轉(zhuǎn)速?zèng)Q定,相關(guān)研究表明[20],平整度還與旋耕刀正切刃形狀、切土節(jié)距、耕深有關(guān)。

    土壤與旋耕刀的摩擦系數(shù)主要由土壤類型決定,根據(jù)DB13/T 1053-2009《山前平原區(qū)小麥-玉米減蒸降耗節(jié)水高產(chǎn)技術(shù)規(guī)程》和NY/T 1409-2007《旱地玉米機(jī)械化保護(hù)性技術(shù)規(guī)范》,小麥-玉米輪作區(qū)的土壤主要為沙壤土和黏壤土。

    旋耕刀折彎角、回轉(zhuǎn)半徑1和正切刃形狀由旋耕刀型號(hào)決定,根據(jù)GB/T 5669-2008《旋耕機(jī)械刀和刀座》,水旱田耕作常用的旋耕刀有IT225、IT245和IT260三種型號(hào)。

    刀軸轉(zhuǎn)速由式(13)計(jì)算:

    式中為旋耕刀轉(zhuǎn)速,r/min;0為機(jī)具前進(jìn)速度,m/s;0為切土區(qū)域內(nèi)圓周設(shè)置的刀片數(shù),對(duì)稱布置取2;0為切土節(jié)距,cm。

    研究表明[21],機(jī)具前進(jìn)速度為1.5 m/s時(shí),破茬裝置的破茬率達(dá)100%,根據(jù)土壤特性和《農(nóng)業(yè)機(jī)械設(shè)計(jì)手冊(cè)》,選擇切土節(jié)距分別為15、17.5和20 cm,求出對(duì)應(yīng)刀輥轉(zhuǎn)速分別為225、260和300 r/min。

    由以上分析可知,土壤類型、旋耕刀型號(hào)、刀輥轉(zhuǎn)速是影響整地機(jī)作業(yè)后地表平整度的另3個(gè)關(guān)鍵因素。

    3 基于田口法與離散元仿真的刀輥參數(shù)優(yōu)化

    在理論分析基礎(chǔ)上,借助離散元仿真獲取試驗(yàn)數(shù)據(jù),利用田口法設(shè)計(jì)試驗(yàn)并確定最佳平整度對(duì)應(yīng)的最優(yōu)參數(shù)組合,以達(dá)到減少試驗(yàn)次數(shù)、降低試驗(yàn)成本、提高地表平整度和機(jī)具穩(wěn)健性的目的。

    3.1 土槽模型的建立

    因?yàn)榉鬯榈慕斩捄透缭谕寥乐姓急容^小且大部分被翻埋土壤覆蓋,對(duì)刀輥勻土旋平作業(yè)效果影響很小,故建立僅含土壤的土槽模型分析刀輥的作業(yè)效果。沙壤土破碎后土壤顆粒保持松散狀態(tài),因此選擇Hertz Mindlin with Bonding接觸模型[22-23],相關(guān)研究表明[24-25],沙壤土的泊松比為0.35,剪切模量為1′106Pa,密度為2 550 kg/m3;土壤顆粒之間的恢復(fù)系數(shù)為0.6,靜摩擦系數(shù)為0.28,動(dòng)摩擦系數(shù)為0.11;土壤與機(jī)具元件間的恢復(fù)系數(shù)為0.6,靜摩擦系數(shù)為0.5,動(dòng)摩擦系數(shù)為0.05。

    黏壤土破碎后土壤顆粒間存在一定的黏性,因此選擇Hertz Mindlin with JKR接觸模型[26]。相關(guān)研究表明[27-28],黏壤土的泊松比為0.37,剪切模量為1.82′106Pa,密度為2315 kg/m3;土壤顆粒之間的恢復(fù)系數(shù)為0.55,靜摩擦系數(shù)為0.84,動(dòng)摩擦系數(shù)為0.1;土壤與機(jī)具元件間的恢復(fù)系數(shù)為0.3,靜摩擦系數(shù)為0.6,動(dòng)摩擦系數(shù)為0.1,JKR表面能為12.73 J/m2。

    機(jī)具材料為鋼材,其本征參數(shù)[29]為:泊松比0.3,剪切模量7.9×1010Pa,密度7 865 kg/m3。土壤與機(jī)具接觸時(shí)僅考慮基本的接觸力,所以沙壤土顆粒、黏壤土顆粒與機(jī)械元件接觸模型均采用無(wú)滑動(dòng)的Hertz Mindlin模型[30]。為了提高計(jì)算速度,沙壤土和黏壤土土壤顆粒半徑均設(shè)為10 mm,然后在離散元仿真軟件EDEM中分別建立長(zhǎng)寬高為3 000 mm′2 400 mm′400 mm的沙壤土土槽和黏壤土土槽。

    3.2 田口法試驗(yàn)設(shè)計(jì)

    在Mintab中采用田口法設(shè)計(jì)試驗(yàn),以作業(yè)后的地表平整度作為優(yōu)化目標(biāo),考慮作業(yè)土壤黏性的不可控性,將土壤類型設(shè)為噪聲因子,螺旋頭數(shù)、旋耕刀型號(hào)、刀軸轉(zhuǎn)速、螺旋升角、刀軸直徑、刀座間距為影響因素;各因素及水平如表1所示,根據(jù)表1設(shè)計(jì)L18(21′35)正交表作為田口試驗(yàn)的內(nèi)表,土壤類型沙壤土和黏壤土為田口試驗(yàn)的外表,外表對(duì)應(yīng)的2列仿真試驗(yàn)數(shù)據(jù)即為地表平整度響應(yīng)列。

    表1 試驗(yàn)因素及水平

    為保證刀輥?zhàn)鳂I(yè)前土槽土壤顆粒分布狀態(tài)相同,破茬、翻埋和勻土旋平作業(yè)獨(dú)立進(jìn)行,在Solidworks中分別建立整地機(jī)破茬裝置、翻埋裝置和勻土旋平裝置的三維模型。按1.1節(jié)中各作業(yè)部件入土深度,依次進(jìn)行破茬作業(yè)(圖4a)、翻埋作業(yè)(圖4b)和勻土旋平作業(yè)仿真(圖4c),由于勻土旋平裝置中的刮土板具有一定的刮平和鎮(zhèn)壓作用,為了更好地體現(xiàn)勻土旋平刀輥的整平效果,在三維模型中減少刮土板高度,使其處于土層上方,僅起到遮擋飛濺土壤的作用。

    圖4 離散元仿真模型

    仿真結(jié)束后,利用Clipping切片功能對(duì)勻土旋平作業(yè)后的土壤穩(wěn)定區(qū)進(jìn)行隨機(jī)采樣,每隔50 mm截取長(zhǎng)50 mm的土槽作為一個(gè)樣本,每次仿真試驗(yàn)采集3個(gè)樣本,采樣方法如圖5a所示。將每個(gè)樣本切片的切面分別導(dǎo)入CAXA中,在作業(yè)幅寬范圍內(nèi)每隔150 mm(即小麥窄行播種行距)做該處最高點(diǎn)顆粒的切線,測(cè)量此處切線與土槽底部橫線的距離,記為i(取樣點(diǎn)數(shù)=1~15),測(cè)量方法如圖5b所示,每次試驗(yàn)取3個(gè)樣本,共采集45個(gè)數(shù)據(jù),由式(14)計(jì)算均方差,表征地表平整度。

    注:Xi(i=1~15)為第i點(diǎn)土層厚度測(cè)量值,mm。

    3.3 仿真試驗(yàn)與結(jié)果分析

    按照田口法中正交表各因素組合方式(表2中因素變量1~6列)分別建立勻土旋平裝置的18個(gè)三維模型,再以破茬和翻埋作業(yè)后相同土壤顆粒分布狀態(tài)的土槽為作業(yè)對(duì)象,分別在沙壤土和黏壤土中仿真勻土旋平作業(yè),由式(14)計(jì)算每次作業(yè)的地表平整度并計(jì)入表2中沙壤土和黏壤土地表平整度響應(yīng)列。響應(yīng)列數(shù)值越小說明地表平整度越高,故選擇望小特性進(jìn)行田口試驗(yàn)分析[31],得出信噪比主效應(yīng)、均值主效應(yīng)分別如圖6a、6b所示,信噪比方差和均值方差如表3所示。

    圖6a中信噪比越大的因素水平對(duì)優(yōu)化目標(biāo)貢獻(xiàn)度越大,圖6b中均值越小的因素水平表示地表平整度越高,可知各因素對(duì)地作業(yè)后表平整度的貢獻(xiàn)度由大到小排序?yàn)?>4>2>6>5>1。

    分析表3中信噪比方差和均值方差,刀軸轉(zhuǎn)速3的顯著性指標(biāo)值分別為0.02和0.024,均已達(dá)到小于0.05的顯著水平,并且變異比值相對(duì)都最大,表明刀軸轉(zhuǎn)速對(duì)地表平整度的貢獻(xiàn)度最大;螺旋升角4的顯著性值分別為0.069和0.053,貢獻(xiàn)度次之;旋耕刀型號(hào)2對(duì)地表平整度的貢獻(xiàn)度較小;刀座間距6、刀軸直徑5、螺旋頭數(shù)1的貢獻(xiàn)度更小,對(duì)地表平整度影響相對(duì)更小。

    表2 田口試驗(yàn)設(shè)計(jì)方案及仿真結(jié)果

    圖6 信噪比主效應(yīng)和均值主效應(yīng)

    表3 信噪比方差和均值方差

    注:<0.05表示影響顯著。

    Note:<0.05 indicates significant influence.

    根據(jù)以上分析,為了達(dá)到最佳平整度、增強(qiáng)機(jī)具穩(wěn)健性,需要使信噪比最大化,根據(jù)圖6,刀軸轉(zhuǎn)速3選取3水平300 r/min,螺旋升角4和旋耕刀型號(hào)2分別選擇70°和IT245旋耕刀,刀座間距6對(duì)信噪比和均值的貢獻(xiàn)度相對(duì)較小,為了避免刀間夾土,同時(shí)為了防止旋耕刀與機(jī)具干涉,確定刀座間距為62 mm,此時(shí)每側(cè)旋耕刀位置數(shù)為15,因此圖2中兩側(cè)刀軸各安裝32把旋耕刀;刀軸直徑5的2水平和3水平貢獻(xiàn)度相似,考慮到翻埋犁作業(yè)后土壤硬度低,對(duì)刀軸強(qiáng)度要求低,從制造成本層面考慮,選取2水平80 mm;旋耕刀排列螺旋頭數(shù)1貢獻(xiàn)度相對(duì)更小,考慮旋耕刀的裝配工藝性,選取雙螺旋排列形式。

    按以上優(yōu)化參數(shù)重新建立勻土旋平裝置三維模型,然后分別導(dǎo)入破茬翻埋后的沙壤土和黏壤土土槽模型中,各重復(fù)3次仿真試驗(yàn),結(jié)果取均值,最終得到沙壤土作業(yè)的地表平整度為13.6 mm、黏壤土作業(yè)的地表平整度為14.9 mm,均小于表2中沙壤土和黏壤土相應(yīng)仿真試驗(yàn)地表平整度的最小值,表明田口法優(yōu)化刀輥參數(shù)可靠。

    為了驗(yàn)證田口法優(yōu)化刀輥參數(shù)的穩(wěn)健性,通過土壤直剪試驗(yàn)和土壤堆積試驗(yàn)(堆積角度平均值34°8′)標(biāo)定4.1節(jié)試驗(yàn)田土壤仿真參數(shù),建立試驗(yàn)田土壤參數(shù)的土槽模型,依次仿真破茬、翻埋作業(yè),重復(fù)3次結(jié)果取均值,得到地表平整度為14.2 mm,介于13.6和14.9 mm之間,說明優(yōu)化后的刀輥參數(shù)對(duì)不同土壤具有較強(qiáng)的適應(yīng)能力。

    4 田間試驗(yàn)

    4.1 試驗(yàn)條件

    為驗(yàn)證刀輥優(yōu)化設(shè)計(jì)結(jié)果的可靠性和穩(wěn)健性,采用優(yōu)化后的參數(shù)制造樣機(jī),于2021年10月在河北省石家莊市辛集市馬莊村保高豐農(nóng)場(chǎng)進(jìn)行田間試驗(yàn),試驗(yàn)田為小麥-玉米輪作田,完成玉米收獲及秸稈粉碎還田后進(jìn)行。試驗(yàn)儀器包括東方紅LG1504輪式拖拉機(jī)(額定功率110.5 kW),TJSD-750-II型土壤堅(jiān)實(shí)度測(cè)量?jī)x(范圍0~10 000 kPa,精度50 kPa),RS-TRREC-N01型土壤水份測(cè)量?jī)x(范圍0~95%RH,精度2%~3%),電子稱(量程3 kg,精度0.1 g),激光測(cè)量?jī)x(量程700 m,示值誤差±0.5 m),卷尺(量程50 m,精度0.001 m),直尺(量程100 cm,精度0.1 cm),環(huán)刀(容積100 cm3),土壤篩(篩孔直徑4 cm),水平尺(長(zhǎng)度600 mm),秒表(精度0.01 s)。經(jīng)測(cè)量,厚度在0~15 cm內(nèi)的試驗(yàn)田耕層土壤堅(jiān)實(shí)度1 159.4 kPa,密度1 257 kg/m3,平均含水率11.6%,含沙率69.8%。田間試驗(yàn)如圖7所示,在耕作秸稈堆積地段時(shí)偶有少量秸稈短暫滯留現(xiàn)象,但經(jīng)過該地段后會(huì)逐漸散落,試驗(yàn)時(shí)未安裝鎮(zhèn)壓裝置,一是為了與仿真試驗(yàn)進(jìn)行對(duì)比,二是為了避免鎮(zhèn)壓影響真實(shí)旋平效果。

    a. 聯(lián)合整地機(jī)作業(yè)試驗(yàn)a. Combined tilling machine operation testb. 作業(yè)后地表b. Surface after operation

    4.2 試驗(yàn)方法

    采用地表平整度、土壤分布均勻度、耕深穩(wěn)定性和破土率4個(gè)指標(biāo)衡量勻土旋平刀輥的綜合作業(yè)性能,其中旋耕深度穩(wěn)定性和破土率按照國(guó)標(biāo)GB/T 5668—2017《旋耕機(jī)》測(cè)定。

    4.2.1 地表平整度

    在非作業(yè)區(qū)域放置激光測(cè)量?jī)x,用水平尺調(diào)平,利用卷尺沿機(jī)具前進(jìn)方向測(cè)量作業(yè)區(qū)域長(zhǎng)度,利用直尺測(cè)出地面與激光測(cè)量?jī)x基準(zhǔn)的相對(duì)高度。作業(yè)前,選定40 m′20 m的區(qū)域作為測(cè)量區(qū),將測(cè)量區(qū)域劃分成2 m′2 m的小區(qū),每個(gè)小區(qū)內(nèi)利用五點(diǎn)采樣法測(cè)量5處高度,取均值代表該小區(qū)地面高度;作業(yè)中,根據(jù)機(jī)具作業(yè)距離和行走用時(shí)計(jì)算實(shí)際行進(jìn)速度;作業(yè)后再次測(cè)量每個(gè)小區(qū)的地面高度。計(jì)算出作業(yè)前、后各小區(qū)相對(duì)于基準(zhǔn)高度的相對(duì)高度,相對(duì)高度的均方差即為地表平整度。

    4.2.2 土壤分布均勻度

    選作業(yè)穩(wěn)定區(qū)域內(nèi)的任一作業(yè)行程,在垂直于機(jī)具前進(jìn)方向每隔2 m豎直開溝,在土壤橫斷面上用環(huán)刀和不銹鋼板每0.1 m采集一個(gè)樣本,每個(gè)斷面采樣21次,20個(gè)橫斷面共采集420個(gè)樣本,用電子稱分別測(cè)量每個(gè)樣本土壤的質(zhì)量,計(jì)算所有樣本的均值和均方差,并按式(15)計(jì)算穩(wěn)定性系數(shù)評(píng)價(jià)土壤分布均勻度。

    4.2.3 耕深穩(wěn)定性系數(shù)

    耕深穩(wěn)定性系數(shù)與土壤分布均勻度同時(shí)測(cè)量,在每個(gè)土壤橫斷面上,用直尺左右各測(cè)量一次旋耕深度,20個(gè)橫斷面共測(cè)量40個(gè)數(shù)據(jù),計(jì)算出所有測(cè)量數(shù)據(jù)的均值和均方差并代入式(15),結(jié)果代表耕深穩(wěn)定性系數(shù)。

    4.2.4 破土率

    在作業(yè)穩(wěn)定區(qū)任取10個(gè)作業(yè)行程,每個(gè)行程內(nèi)隨機(jī)選取0.5 m′0.5 m的測(cè)量區(qū),選出每個(gè)測(cè)量區(qū)旋耕層內(nèi)最長(zhǎng)邊小于4 cm的所有土塊,測(cè)量其質(zhì)量記作1,測(cè)量該區(qū)域內(nèi)所有土壤總質(zhì)量記作2,計(jì)算1/2即得單次行程的破土率,取10個(gè)行程破土率的均值代表破土率。

    4.3 試驗(yàn)結(jié)果與分析

    作業(yè)后地表相對(duì)基站高度如表4所示,根據(jù)式(14)計(jì)算表4數(shù)據(jù)的均方差即作業(yè)后地表平整度為11.6 mm,高于小麥播種整地要求,表明由田口法優(yōu)化得到的刀輥工作參數(shù)具有較高的可靠性和穩(wěn)健性,該值與試驗(yàn)田土壤仿真值14.2 mm相似度為82%,表明仿真過程具有較高的有效性,同時(shí)田間試驗(yàn)數(shù)據(jù)優(yōu)于仿真結(jié)果,分析認(rèn)為是田間作業(yè)過程中柔性刮土板起到一定的刮平和鎮(zhèn)壓作用。

    表4 作業(yè)后試驗(yàn)小區(qū)不同斷面的平均相對(duì)高度

    經(jīng)測(cè)量旋耕層土壤分布均勻度為92.6%、平均旋耕深度為12.7 cm(標(biāo)準(zhǔn)要求大于等于8 cm)、耕深穩(wěn)定性系數(shù)為93.2%(標(biāo)準(zhǔn)要求大于85%)、破土率為86.4%(標(biāo)準(zhǔn)要求大于60%),均已達(dá)到標(biāo)準(zhǔn)要求,表明勻土旋平刀輥既能保證土壤表面平整度,又能保證土壤分布均勻性、旋耕深度穩(wěn)定性和破土率等指標(biāo),所設(shè)計(jì)的勻土旋平刀輥的有效性。

    5 結(jié) 論

    1)設(shè)計(jì)了分段反螺旋式勻土旋平刀輥,分析并優(yōu)化了影響刀輥工作效果的關(guān)鍵參數(shù),提高了土壤軸向分布的均勻性,解決了聯(lián)合整地機(jī)破茬翻埋作業(yè)后,地表呈M形平整度低的問題。

    2)為滿足聯(lián)合整地機(jī)勻土旋平刀輥的穩(wěn)健性設(shè)計(jì)要求,提出了一種田口法結(jié)合離散元仿真優(yōu)化勻土旋平刀輥工作參數(shù)的新思路。得出最優(yōu)工作參數(shù)組合為:刀軸轉(zhuǎn)速300 r/min、螺旋升角70°、旋耕刀型號(hào)IT245、刀軸直徑80 mm、刀座間距62 mm、旋耕刀雙螺旋排布,此組合下仿真試驗(yàn)田土壤作業(yè),得出地表平整度為14.2 mm,表明優(yōu)化參數(shù)的刀輥在不同土壤環(huán)境中具有較高的穩(wěn)健性。

    3)田間試驗(yàn)表明:聯(lián)合整地機(jī)作業(yè)后地表平整度為11.6 mm、土壤分布均勻度92.6%、耕深穩(wěn)定性系數(shù)為93.2%、破土率為86.4%,可為小麥播種、灌溉及生長(zhǎng)提供優(yōu)質(zhì)的土壤環(huán)境;該地表平整度與仿真結(jié)果相似度為82%,證明刀輥設(shè)計(jì)參數(shù)有效且穩(wěn)定。

    [1] 董印麗,李振峰,王若倫,等.華北地區(qū)小麥、玉米兩季秸稈還田存在問題及對(duì)策研究[J]. 中國(guó)土壤與肥料,2018(1):159-163.

    Dong Yinli, Li Zhenfeng, Wang Ruolun, et al. Study on the problems and countermeasures of returning wheat and corn stalks into the soil in north China[J]. Soil and Fertilizer in China, 2018(1): 159-163.(in Chinese with English abstract)

    [2] 河北省市場(chǎng)監(jiān)督管理局. 小麥玉米節(jié)水豐產(chǎn)品種高效組合模式應(yīng)用技術(shù)規(guī)程:DB13/T 5440-2021[S]. 2021.

    [3] 張欣悅,李連豪,汪春,等.1GSZ-350型滅茬旋耕聯(lián)合整地機(jī)的設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2009,25(5):73-77.

    Zhang Xinyue, Li Lianhao, Wang Chun, et al. Design and test of 1GSZ-350 stubble-breaking and rotary tilling combined cultivating machine[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2009, 25(5): 73-77.(in Chinese with English abstract)

    [4] 車剛,張偉,萬(wàn)霖,等.基于滅茬圓盤驅(qū)動(dòng)旋耕刀多功能耕整機(jī)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(20):34-40.

    Che Gang,Zhang Wei, Wan Lin, et al. Design and experiment of multifunctional tillage machine with driven bent blade by stubble ploughing disk[J]. Transactions of Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012,28(20): 34-40.(in Chinese with English abstract)

    [5] 鄭侃,何進(jìn),李洪文,等.反旋深松聯(lián)合作業(yè)耕整機(jī)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(8):61-71.

    Zheng Kan, He Jin, Li Hongwen, et al. Design and experiment of combined tillage implement of reverse-rotary and subsoiling[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017,48(8): 61-71.(in Chinese with English abstract)

    [6] 張秀梅,夏俊芳,張居敏,等. 水旱兩用秸稈還田組合刀輥?zhàn)鳂I(yè)性能試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(9):9-15.

    Zhang Xiumei, Xia Junfang, Zhang Jumin, et al. Working performance experiment of combination blade roller for straw returning in paddy field and dry land[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(9): 9-15.(in Chinese with English abstract)

    [7] 王金武,張成亮,許春林,等.聯(lián)合整地機(jī)平整部件參數(shù)優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(2):34-37,11.

    Wang Jinwu, Zhang Chengliang, Xu Chunlin, et al. Parameters optimization on flattening component of combined cultivating implement[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(2): 34-37, 11.(in Chinese with English abstract)

    [8] 鄭侃,李宇飛,夏俊芳,等.開溝旋耕機(jī)漸變螺旋升角軸向勻土刀輥設(shè)計(jì)與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2021,52(5):63-73.

    Zheng Kan, Li Yufei, Xia Junfang, et al. Design and experiment of land leveling blade roller of ditching and rotary tiller with gradual spiral angle[J]. Transactions of the Chinese Society for Agricultural Machinery, 2021, 52(5): 63-73.(in Chinese with English abstract)

    [9] 周華,張居敏,祝英豪,等.秸稈還田深松旋埋聯(lián)合耕整機(jī)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(22):17-26.

    Zhou Hua, Zhang Jumin, Zhu Yinghao, et al. Design and experiment of combined tillage machine for subsoiling androtary burying of straw incorporated into soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(22): 17-26. (in Chinese with English abstract)

    [10] Zheng K, Mchugh A D, Li H W, et al. Design and experiment of anti-vibrating and anti-wrapping rotary components for subsoiler cum rotary tiller[J]. International Journal of Agricultural and Biological Engineering,2019, 12(4): 47-55.

    [11] Yu C Y,Liu J,Zhang J, et al. Design and optimization and experimental verification of a segmented double-helix blade roller for straw returning cultivators[J]. Journal of the Chinese Institute of Engineers, 2021, 44(3): 1-9.

    [12] 林秀雄. 田口方法實(shí)戰(zhàn)技術(shù)[M]. 深圳:海天出版社,2008.

    [13] 姜雪松,張坤,張紅麗,等.基于田口試驗(yàn)設(shè)計(jì)的TiAlN涂層膜基結(jié)合力控制參數(shù)優(yōu)化[J]. 中國(guó)電子科學(xué)研究院學(xué)報(bào),2020,15(7):698-702.

    Jiang Xuesong,Zhang Kun,Zhang Hongli, et al. Optimization of control parameters of adhesion strength between TiAlN coating and substrate based on taguchi experimental design[J]. Journal of China Electronics Research Institute, 2020, 15(7): 698-702.(in Chinese with English abstract)

    [14] 吳曉君,陳竹,周天擇,等.彈性磨具高效磨拋M300鋼曲面的接觸特征研究[J]. 機(jī)械工程學(xué)報(bào),2018,54(1):171-177.

    Wu Xiaojun,Chen Zhu,Zhou Tianze, et al. Research on processing efficiency and contact characteristics of M300steel surface grinding with elastic abrasives[J]. Journal of Mechanical Engineering, 2018, 54(1): 171-177. (in Chinese with English abstract)

    [15] 蒙賀偉,高振江,坎雜,等. 等徑變螺距奶牛精確飼喂給料裝置設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2011,27(3):103-107.

    Meng Hewei, Gao Zhenjiang, Kan Za, et al. Design and experiment on dairy cow precision-feeding device based on equal-diameter and variable-pitch[J]. Transactions of the Chinese Society of Agricultural Engineering ((Transactions of the CSAE), 2011, 27(3): 103-107.(in Chinese with English abstract)

    [16] 秦寬,丁為民,方志超,等. 復(fù)式耕整機(jī)耕深與耕寬穩(wěn)定性分析與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(9):1-8.

    Qin Kuan, Ding Weimin, Fang Zhichao, et al. Analysis and experiment of tillage depth and width stability for plowing and rotary tillage combined machine[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(9): 1-8.(in Chinese with English abstract)

    [17] 楊文武,羅錫文,王在滿,等.輪式拖拉機(jī)水田輪轍覆土裝置設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(16):26-31.

    Yang Wenwu,Luo Xiwen,Wang Zaiman, et al. Design and experiment of track filling assembly mounted on wheeled-tractor for paddy fields[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016,32(16): 26-31.(in Chinese with English abstract)

    [18] 高文英,林靜,李寶筏,等.玉米秸稈深埋還田機(jī)螺旋開溝裝置參數(shù)優(yōu)化與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2018,49(9):45-54.

    Gao Wenying,Lin Jing,Li Baofa, et al. Parameter optimization and experiment for spiral type openerdevice of maize straw deep bury and returning machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018,49(9): 45-54.(in Chinese with English abstract)

    [19] 王金峰,陳博聞,姜巖,等.水稻秸稈全量深埋還田機(jī)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2020,51(1):84-93.

    Wang Jinfeng, Chen Bowen, Jiang Yan, et al. Design and experiment on machine for rice straw full quantity deep buried into field[J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(1): 84-93. (in Chinese with English abstract)

    [20] 陳鈞,近江谷和彥,寺尾日出男.高速攝影法研究旋耕刀拋土特性[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),1994,25(3):56-60.

    Chen Jun,Ohmiya Kazuhiko,Terao Hideo.Studies on the soil-throwing property of the rotary blades with high speed photography[J]. Transactions of the Chinese Society for Agricultural Machinery, 1994, 25(3): 56-60.(in Chinese with English abstract)

    [21] Zhang X H, Wang J Z , Lao F X, et al. Simulation and optimization of working parameters of stubble breaking device in two ripening area of yumai in north China[C]//2021 ASABE Annual International Virtual Meeting. Michigan: American Society of Agricultural and Biological Engineers, 2021: 1.

    [22] 方會(huì)敏,姬長(zhǎng)英,F(xiàn)arman Ali Chandio,等. 基于離散元法的旋耕過程土壤運(yùn)動(dòng)行為分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(3):22-28.

    Fang Huimin, Ji Changying, Farman Ali Chandio, et al.Analysis of soil dynamic behavior during rotary tillagebased on distinct element method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016,47(3): 22-28.(in Chinese with English abstract)

    [23] 方會(huì)敏,姬長(zhǎng)英,Ahmed Ali Tagar,等.秸稈-土壤-旋耕刀系統(tǒng)中秸稈位移仿真分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(1):60-67.

    Fang Huimin, Ji Changying, Ahmed Ali Tagar, et al.Simulation analysis of straw movement in straw-soil-rotary blade system[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016,47(1): 60-67.(in Chinese with English abstract)

    [24] 王憲良,胡紅,王慶杰,等.基于離散元的土壤模型參數(shù)標(biāo)定方法[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(12):78-85.

    Wang Xianliang, Hu Hong, Wang Qingjie, et al. Calibration method of soil contact characteristic parameters based on dem theory[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(12): 78-85.(in Chinese with English abstract)

    [25] 石林榕,趙武云,孫偉.基于離散元的西北旱區(qū)農(nóng)田土壤顆粒接觸模型和參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(21):181-187.

    Shi Linrong, Zhao Wuyun, Sun Wei. Parameter calibration of soil particles contact model of farmland soil in northwest arid region based on discrete element method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2017, 33(21): 181-187.(in Chinese with English abstract)

    [26] Johnson K L,Kendall K, Roberts A D.Surface energy and the contact of elastic solids[J].Proceedings of the Royal Society of London, Series A: Mathematical and Physical Sciences, 1971, 324(1558): 301-313.

    [27] 祝英豪,夏俊芳,曾榮,等. 基于離散元的稻板田旋耕功耗預(yù)測(cè)模型研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2020,51(10):42-50.

    Zhu Yinghao, Xia Junfang, Zeng Rong, et al. Prediction model of rotary tillage power consumption in paddy stubble field based on discrete element method[J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(10): 42-50.(in Chinese with English abstract)

    [28] 向偉,吳明亮,呂江南,等.基于堆積試驗(yàn)的黏壤土仿真物理參數(shù)標(biāo)定[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(12):116-123.

    Xiang Wei, Wu Mingliang, Lyu Jiangnan, et al. Calibration of simulation physical parameters of clay loam based on soil accumulation test[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(12): 116-123.(in Chinese with English abstract)

    [29] Ucgul M, Fielke J M, Saunders C. 3D DEM tillage simulation: Validation of a hysteretic spring (plastic) contact model for a sweep tool operating in a cohesionless soil[J]. Soil and Tillage Research, 2014, 144: 220-227.

    [30] Thankur S C,Morrissey J P, Sun J, et al.Micromechanical analysis of cohesive granular materials using the discrete element method with an adhesive elasto-plastic contact model[J]. Granular Matter, 2014, 16(3): 383-400.

    [31] 劉春景,唐敦兵,何華,等. 基于灰色關(guān)聯(lián)和主成分分析的車削加工多目標(biāo)優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(4):293-298,292.

    Liu Chunjing,Tang Dunbing,He Hua, et al. Multi-objective optimization of turning based on grey relational and principal component analysis[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013,44(4): 293-298, 292.(in Chinese with English abstract)

    Design and tests of the soil moving and leveling blade roller for combined tilling machines

    Liu Shangkun, Liu Chao, Zhang Xiuhua, Liu Jiangtao, Yi Jinggang, Zhang Jinguo

    (1.,,071001,; 2,,071001,)

    Combined tillage can be utilized to implement two or more tillage operations in the field at the same time. The efficient and potential way can be widely used to reduce soil compaction, labor, and fuel cost. However, the current combined tilling machine cannot fully meet the harsh requirement of better surface flatness and soil distribution after stubble breaking and burying operation. In this study, a piecewise anti-helix blade roller was designed to improve the parameter robustness using Taguchi method and discrete element simulation. Firstly, the structural scheme of moving and leveling soil blade rollers was designed, according to the M-shaped surface produced by the ordinary blade roller in the previous generation prototype test. Specifically, the gearbox was put in the middle of the frame with a blade roller at each side. Each blade roller was designed as two sections with opposite and unequal lengths of spiral directions. As such, the middle and each side ridge ditches were filled for the moving and leveling of the two raised soil ridges that formed by the burying operation. Secondly, some key factors were qualitatively analyzed to optimize the work parameters of the piecewise anti-helix blade roller. The specific level of the factor was determined in the axial soil transportation and throwing off the blade roller. Taguchi method was used to design the test scheme, where the surface flatness was taken as the optimization objective, and the sandy loam and clay loam were the noise factors. The three-dimensional models were established using the test parameters after the discrete element simulation. The optimal combination of parameters was obtained suitable for the different soil types. Specifically, the cutter shaft speed and diameter, the screw angle, and the blade seat spacing were 300 r/min, 80 mm, 70°, and 62 mm, respectively. The blade roller structure was the double helix mode, with the rotary blade model of IT245. The soil operation was then simulated under the optimal combination in the farm test field. The surface flatness was achieved at 14.2mm, indicating the robustness of the optimal parameter of the blade roller. Finally, the prototype was manufactured with the optimal combination of parameters, and then verified by the field test. After that, the surface flatness, soil distribution uniformity, stability coefficient of rotary tillage depth, and soil breaking rate were 11.6 mm, 92.6%, 93.2%, and 86.4%, respectively, which were better than the requirements of wheat cultivation for the tillage and land preparation. The better performance of the blade roller was achieved under the optimal combination of parameters. The finding can provide a promising approach to designing the blade roller in a combined tilling machine.

    discrete element; Taguchi method; combined tilling machine; piecewise anti-helix blade roller; surface flatness

    10.11975/j.issn.1002-6819.2022.17.001

    S224.4

    A

    1002-6819(2022)-17-0001-09

    劉尚坤,劉超,張秀花,等. 聯(lián)合整地機(jī)勻土旋平刀輥設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(17):1-9.doi:10.11975/j.issn.1002-6819.2022.17.001 http://www.tcsae.org

    Liu Shangkun, Liu Chao, Zhang Xiuhua, et al. Design and tests of the soil moving and leveling blade roller for combined tilling machines[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(17): 1-9. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.17.001 http://www.tcsae.org

    2022-06-09

    2022-08-18

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0300907);河北省重大科技成果轉(zhuǎn)化專項(xiàng)(20287201Z);河北省重點(diǎn)研發(fā)計(jì)劃(22327206D);河北省省屬高?;究蒲袠I(yè)務(wù)費(fèi)項(xiàng)目(KY202014);河北農(nóng)業(yè)大學(xué)引進(jìn)人才科研專項(xiàng)(YJ201814)

    劉尚坤,博士,副教授,研究方向?yàn)闄C(jī)械裝備設(shè)計(jì)制造與優(yōu)化。Email:lsk1213@163.com

    猜你喜歡
    刀輥刀軸壤土
    基于響應(yīng)面分析法的輥?zhàn)邮介蛔悠茪C(jī)結(jié)構(gòu)參數(shù)的優(yōu)化
    微耕機(jī)旋耕刀輥彎刀排列優(yōu)化仿真分析*
    復(fù)雜曲面五軸數(shù)控?zé)o干涉刀軸規(guī)劃
    土壤質(zhì)地及砧木影響蘋果根際微生物功能多樣性及其碳源利用
    雙軸式旋耕滅茬機(jī)滅茬刀輥結(jié)構(gòu)設(shè)計(jì)與試驗(yàn)
    左家林場(chǎng)核桃良種基地選址調(diào)查報(bào)告
    CONTENTS
    全斷面隧道掘進(jìn)機(jī)(TBM)盤形滾刀刀軸斷裂失效分析
    瓦線橫切機(jī)刀軸偏心質(zhì)量對(duì)切刀磨損影響的研究
    基于強(qiáng)度比較的微耕機(jī)刀軸的優(yōu)化設(shè)計(jì)
    啦啦啦韩国在线观看视频| 精品一区二区三区视频在线| 好男人在线观看高清免费视频| 黄色视频,在线免费观看| 色综合色国产| 天天躁日日操中文字幕| 有码 亚洲区| 欧美xxxx黑人xx丫x性爽| 久久99热6这里只有精品| 国内精品久久久久精免费| 亚洲国产精品合色在线| 搡老岳熟女国产| 亚洲av日韩精品久久久久久密| 日韩高清综合在线| 亚洲精品乱码久久久v下载方式| 日韩在线高清观看一区二区三区 | 免费电影在线观看免费观看| 欧美xxxx黑人xx丫x性爽| 婷婷丁香在线五月| 不卡视频在线观看欧美| 久久久午夜欧美精品| 长腿黑丝高跟| 毛片女人毛片| 中文在线观看免费www的网站| 久久久久国产精品人妻aⅴ院| 精品人妻1区二区| 国产av不卡久久| www.www免费av| 色视频www国产| 成人鲁丝片一二三区免费| 成人性生交大片免费视频hd| 欧美黑人巨大hd| 欧美性猛交黑人性爽| 精品午夜福利在线看| 日韩欧美精品v在线| 神马国产精品三级电影在线观看| 日日啪夜夜撸| 久久香蕉精品热| 我的老师免费观看完整版| 男女做爰动态图高潮gif福利片| 在线观看美女被高潮喷水网站| 人妻少妇偷人精品九色| 亚洲中文日韩欧美视频| 亚洲精品一区av在线观看| 欧美bdsm另类| 久久天躁狠狠躁夜夜2o2o| 美女高潮的动态| 99久久精品一区二区三区| 不卡一级毛片| 国产男人的电影天堂91| 国产精华一区二区三区| 啦啦啦观看免费观看视频高清| 小说图片视频综合网站| 99久久精品热视频| 午夜福利在线观看吧| 亚洲av电影不卡..在线观看| 亚洲最大成人中文| eeuss影院久久| 我要搜黄色片| 老司机福利观看| 欧美人与善性xxx| 午夜免费成人在线视频| 偷拍熟女少妇极品色| 一级黄片播放器| 日日夜夜操网爽| 999久久久精品免费观看国产| 91麻豆精品激情在线观看国产| 一个人观看的视频www高清免费观看| 又爽又黄无遮挡网站| 91狼人影院| 久久欧美精品欧美久久欧美| 亚洲狠狠婷婷综合久久图片| 国产精品伦人一区二区| 亚洲av五月六月丁香网| 他把我摸到了高潮在线观看| 午夜免费激情av| 99久久成人亚洲精品观看| 国产乱人伦免费视频| 男人舔女人下体高潮全视频| videossex国产| 国产精品久久电影中文字幕| 国产亚洲精品久久久久久毛片| 一进一出好大好爽视频| 男人和女人高潮做爰伦理| 日本一本二区三区精品| ponron亚洲| 麻豆国产97在线/欧美| 久9热在线精品视频| 精品久久久噜噜| 欧美三级亚洲精品| 国产国拍精品亚洲av在线观看| 免费看美女性在线毛片视频| 两个人视频免费观看高清| 毛片一级片免费看久久久久 | 亚洲欧美日韩卡通动漫| 99久久无色码亚洲精品果冻| 亚洲av第一区精品v没综合| 国产真实乱freesex| 国产色爽女视频免费观看| 国产熟女欧美一区二区| 亚洲国产欧美人成| 色视频www国产| 国产免费男女视频| 国内少妇人妻偷人精品xxx网站| 国产午夜福利久久久久久| 国产午夜福利久久久久久| 夜夜看夜夜爽夜夜摸| 亚洲精品在线观看二区| 婷婷六月久久综合丁香| 99久久精品热视频| 性色avwww在线观看| 老师上课跳d突然被开到最大视频| av福利片在线观看| 久久99热这里只有精品18| 亚洲 国产 在线| 中文字幕熟女人妻在线| 国产精品久久久久久精品电影| 99久久精品热视频| 欧美国产日韩亚洲一区| 成人性生交大片免费视频hd| 成年人黄色毛片网站| 国产精品一区www在线观看 | 女的被弄到高潮叫床怎么办 | 成人无遮挡网站| 国产精品伦人一区二区| 日韩中文字幕欧美一区二区| 国产精品一区二区三区四区久久| 中文字幕精品亚洲无线码一区| 窝窝影院91人妻| 日韩精品青青久久久久久| 国产一区二区在线av高清观看| 久99久视频精品免费| 免费av观看视频| 无遮挡黄片免费观看| 亚洲七黄色美女视频| 亚洲久久久久久中文字幕| 三级男女做爰猛烈吃奶摸视频| 高清日韩中文字幕在线| 国产成人av教育| 尾随美女入室| 男女做爰动态图高潮gif福利片| 亚洲va在线va天堂va国产| 十八禁国产超污无遮挡网站| 午夜影院日韩av| 欧美激情在线99| 在线播放国产精品三级| 亚洲精品亚洲一区二区| 亚洲精品影视一区二区三区av| 国产大屁股一区二区在线视频| 永久网站在线| 亚洲av日韩精品久久久久久密| 国产亚洲欧美98| 欧美不卡视频在线免费观看| 精品一区二区免费观看| 非洲黑人性xxxx精品又粗又长| 国产探花极品一区二区| 深夜精品福利| 国产视频内射| 亚洲最大成人中文| 国产精品人妻久久久影院| 国产精品一区www在线观看 | 亚洲精品色激情综合| 国产精品女同一区二区软件 | 熟妇人妻久久中文字幕3abv| 九九在线视频观看精品| 神马国产精品三级电影在线观看| 欧美bdsm另类| 女生性感内裤真人,穿戴方法视频| 亚洲国产精品成人综合色| 成人av在线播放网站| 欧美极品一区二区三区四区| 国产免费一级a男人的天堂| 女生性感内裤真人,穿戴方法视频| 桃色一区二区三区在线观看| 在线天堂最新版资源| 亚洲内射少妇av| 韩国av在线不卡| 日韩,欧美,国产一区二区三区 | 亚洲人成网站高清观看| 国产午夜精品久久久久久一区二区三区 | av在线亚洲专区| avwww免费| 午夜免费激情av| 中文字幕高清在线视频| 亚洲精品成人久久久久久| 性插视频无遮挡在线免费观看| 精品一区二区三区av网在线观看| 搡老岳熟女国产| 性色avwww在线观看| 亚洲av中文字字幕乱码综合| 99久久精品热视频| 夜夜夜夜夜久久久久| 99久久九九国产精品国产免费| 久久精品国产清高在天天线| 国产亚洲欧美98| 亚洲专区国产一区二区| 琪琪午夜伦伦电影理论片6080| 中出人妻视频一区二区| 日本黄大片高清| 亚洲性夜色夜夜综合| 中文字幕免费在线视频6| 国产一级毛片七仙女欲春2| 欧美性猛交╳xxx乱大交人| 亚洲人成网站在线播| 亚洲第一区二区三区不卡| 久久久久久久午夜电影| 99久久成人亚洲精品观看| 精品人妻偷拍中文字幕| 一边摸一边抽搐一进一小说| 久久久久免费精品人妻一区二区| 色吧在线观看| 免费看日本二区| 国产精品亚洲一级av第二区| 超碰av人人做人人爽久久| av黄色大香蕉| 国产男靠女视频免费网站| 日本免费a在线| 日本欧美国产在线视频| 人人妻,人人澡人人爽秒播| 国产男靠女视频免费网站| 2021天堂中文幕一二区在线观| 国产一区二区在线av高清观看| 麻豆精品久久久久久蜜桃| 亚洲性久久影院| 一区二区三区高清视频在线| 日韩欧美精品v在线| 国产欧美日韩精品亚洲av| 色精品久久人妻99蜜桃| 极品教师在线视频| 国产精品人妻久久久久久| 国内精品美女久久久久久| 久久久久久伊人网av| 国产高清有码在线观看视频| 午夜激情福利司机影院| 午夜精品一区二区三区免费看| 一个人看视频在线观看www免费| av在线观看视频网站免费| 亚洲黑人精品在线| 韩国av一区二区三区四区| 亚洲综合色惰| 亚洲欧美精品综合久久99| 精品不卡国产一区二区三区| 日日撸夜夜添| 91久久精品国产一区二区成人| 国产69精品久久久久777片| 中文字幕高清在线视频| 亚洲国产高清在线一区二区三| 免费看a级黄色片| 日韩欧美三级三区| 久久久久国内视频| 内地一区二区视频在线| 久久久久国产精品人妻aⅴ院| 国产精品一区二区三区四区久久| 在线免费观看不下载黄p国产 | 最近最新免费中文字幕在线| 久久久久久久精品吃奶| 午夜福利视频1000在线观看| 老女人水多毛片| 全区人妻精品视频| 国产蜜桃级精品一区二区三区| 男女之事视频高清在线观看| 欧美中文日本在线观看视频| 大型黄色视频在线免费观看| 国内精品宾馆在线| 婷婷精品国产亚洲av| 少妇的逼好多水| 九九爱精品视频在线观看| 久久精品人妻少妇| 99在线视频只有这里精品首页| 色综合色国产| 色综合站精品国产| 国产成人a区在线观看| 精品久久久久久久久亚洲 | 91午夜精品亚洲一区二区三区 | 在线观看av片永久免费下载| 99久久精品国产国产毛片| 日本熟妇午夜| 俄罗斯特黄特色一大片| 天堂av国产一区二区熟女人妻| 99久久精品热视频| 如何舔出高潮| 自拍偷自拍亚洲精品老妇| 久久九九热精品免费| 成人高潮视频无遮挡免费网站| 麻豆国产97在线/欧美| 99久久成人亚洲精品观看| 国产午夜精品久久久久久一区二区三区 | 国国产精品蜜臀av免费| 中文字幕人妻熟人妻熟丝袜美| 国产精品无大码| 成人性生交大片免费视频hd| 麻豆成人av在线观看| 国产精品自产拍在线观看55亚洲| 国产成人a区在线观看| 少妇丰满av| 午夜激情欧美在线| 国产爱豆传媒在线观看| 亚洲最大成人手机在线| 欧美区成人在线视频| 免费大片18禁| 亚洲午夜理论影院| 综合色av麻豆| 欧美潮喷喷水| 国产精品一区二区三区四区久久| 能在线免费观看的黄片| 精品一区二区三区人妻视频| 99久久精品国产国产毛片| 亚洲av五月六月丁香网| 亚洲av熟女| 午夜精品一区二区三区免费看| 久久久午夜欧美精品| 国内毛片毛片毛片毛片毛片| 性欧美人与动物交配| 国产精品99久久久久久久久| 久久这里只有精品中国| www.www免费av| 观看免费一级毛片| 亚洲第一区二区三区不卡| 美女高潮喷水抽搐中文字幕| 久久午夜亚洲精品久久| 又黄又爽又刺激的免费视频.| 色哟哟哟哟哟哟| 日韩强制内射视频| xxxwww97欧美| 日韩,欧美,国产一区二区三区 | 18+在线观看网站| 国产精品人妻久久久影院| 色尼玛亚洲综合影院| 欧美日韩精品成人综合77777| 在线看三级毛片| 日韩欧美 国产精品| 最后的刺客免费高清国语| 少妇高潮的动态图| 最近中文字幕高清免费大全6 | 长腿黑丝高跟| 亚洲欧美日韩高清在线视频| 高清日韩中文字幕在线| 乱码一卡2卡4卡精品| 国内揄拍国产精品人妻在线| 中文字幕免费在线视频6| 精品久久久久久,| 一进一出抽搐动态| 国产爱豆传媒在线观看| 精品人妻熟女av久视频| 欧美一区二区精品小视频在线| 久久久久九九精品影院| 亚洲成人免费电影在线观看| a在线观看视频网站| 亚洲av五月六月丁香网| 久久6这里有精品| 国产精品亚洲美女久久久| 露出奶头的视频| 啦啦啦韩国在线观看视频| 亚洲av免费在线观看| 国内久久婷婷六月综合欲色啪| 男女下面进入的视频免费午夜| 亚洲一区高清亚洲精品| 亚洲美女黄片视频| 成人二区视频| 2021天堂中文幕一二区在线观| 18禁裸乳无遮挡免费网站照片| 久久久色成人| 国产欧美日韩一区二区精品| 中文字幕久久专区| 国产爱豆传媒在线观看| 99久久久亚洲精品蜜臀av| а√天堂www在线а√下载| 中文字幕人妻熟人妻熟丝袜美| 大型黄色视频在线免费观看| 国产伦人伦偷精品视频| 精品一区二区三区av网在线观看| 国产午夜精品久久久久久一区二区三区 | 成人一区二区视频在线观看| 美女高潮喷水抽搐中文字幕| 国产色婷婷99| 欧美黑人欧美精品刺激| 真实男女啪啪啪动态图| 国产高清视频在线观看网站| 欧美最黄视频在线播放免费| 精品久久久久久久久av| 久久久久久久久中文| 国产欧美日韩精品亚洲av| 欧美另类亚洲清纯唯美| av在线蜜桃| 制服丝袜大香蕉在线| 国产又黄又爽又无遮挡在线| 亚洲第一区二区三区不卡| 综合色av麻豆| 免费黄网站久久成人精品| 亚洲一区二区三区色噜噜| 国产精品自产拍在线观看55亚洲| 国产在线精品亚洲第一网站| 日韩一本色道免费dvd| 日本 欧美在线| 啦啦啦韩国在线观看视频| 亚洲成人中文字幕在线播放| 日本 欧美在线| 日韩欧美在线乱码| 日韩欧美一区二区三区在线观看| 日本撒尿小便嘘嘘汇集6| 日本黄大片高清| 欧美色视频一区免费| 国产老妇女一区| 国产精品精品国产色婷婷| 亚洲第一电影网av| 99国产极品粉嫩在线观看| 亚洲国产精品成人综合色| 欧美高清成人免费视频www| 九九久久精品国产亚洲av麻豆| 简卡轻食公司| 亚洲经典国产精华液单| 有码 亚洲区| 亚洲美女搞黄在线观看 | 国产一区二区在线观看日韩| 国产精品精品国产色婷婷| 亚洲av日韩精品久久久久久密| 国产人妻一区二区三区在| 久久久久久久久久成人| 美女大奶头视频| 俺也久久电影网| 韩国av在线不卡| 少妇人妻精品综合一区二区 | 国产在线男女| 琪琪午夜伦伦电影理论片6080| 久久久久久久久中文| av在线天堂中文字幕| 内射极品少妇av片p| 两个人视频免费观看高清| 成年版毛片免费区| 午夜福利视频1000在线观看| 午夜福利高清视频| 成年免费大片在线观看| 色噜噜av男人的天堂激情| 别揉我奶头 嗯啊视频| www.色视频.com| 国产色爽女视频免费观看| 国产欧美日韩精品亚洲av| 日韩一区二区视频免费看| 欧美性猛交黑人性爽| 亚洲无线观看免费| 99国产精品一区二区蜜桃av| 欧美区成人在线视频| 美女cb高潮喷水在线观看| 两个人的视频大全免费| 麻豆精品久久久久久蜜桃| 搡老岳熟女国产| 亚洲精品乱码久久久v下载方式| 丰满的人妻完整版| 日韩大尺度精品在线看网址| 久久99热这里只有精品18| 精品久久久久久成人av| 久久中文看片网| 欧美潮喷喷水| 偷拍熟女少妇极品色| 黄色一级大片看看| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区性色av| 久99久视频精品免费| 免费看a级黄色片| 校园春色视频在线观看| 亚洲av美国av| 欧美高清成人免费视频www| 一夜夜www| 在线观看66精品国产| 久久九九热精品免费| 中文字幕av成人在线电影| 最新在线观看一区二区三区| 精品久久久久久久久亚洲 | 男女边吃奶边做爰视频| 亚洲成av人片在线播放无| 狠狠狠狠99中文字幕| 女同久久另类99精品国产91| 午夜精品久久久久久毛片777| 精品人妻1区二区| 午夜福利成人在线免费观看| 亚州av有码| 99久久精品一区二区三区| 日韩欧美免费精品| 成人国产综合亚洲| 少妇裸体淫交视频免费看高清| 99精品在免费线老司机午夜| 久久婷婷人人爽人人干人人爱| 欧美丝袜亚洲另类 | 欧美成人a在线观看| 天天一区二区日本电影三级| 丰满乱子伦码专区| 最近最新中文字幕大全电影3| eeuss影院久久| 成人二区视频| 又粗又爽又猛毛片免费看| 嫁个100分男人电影在线观看| 亚洲,欧美,日韩| 国产伦人伦偷精品视频| 少妇熟女aⅴ在线视频| 免费人成视频x8x8入口观看| 波野结衣二区三区在线| 一区二区三区免费毛片| 久久精品91蜜桃| 一区二区三区免费毛片| 日本黄色视频三级网站网址| 美女cb高潮喷水在线观看| 国产精品久久久久久久电影| 欧美成人a在线观看| 免费看光身美女| 18禁在线播放成人免费| 国产精品人妻久久久久久| 床上黄色一级片| 亚洲三级黄色毛片| 亚洲无线在线观看| 老熟妇乱子伦视频在线观看| 国内精品宾馆在线| 自拍偷自拍亚洲精品老妇| 久久精品国产鲁丝片午夜精品 | 国产精品亚洲美女久久久| 午夜福利在线观看吧| 日本免费一区二区三区高清不卡| 特大巨黑吊av在线直播| 亚洲熟妇中文字幕五十中出| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| 蜜桃久久精品国产亚洲av| 欧美另类亚洲清纯唯美| 亚洲乱码一区二区免费版| 免费在线观看成人毛片| 日韩欧美精品v在线| 美女黄网站色视频| 久久久国产成人免费| 麻豆成人av在线观看| 亚洲精华国产精华精| 国产精品亚洲一级av第二区| x7x7x7水蜜桃| 免费av毛片视频| 真人做人爱边吃奶动态| 男女视频在线观看网站免费| 22中文网久久字幕| 亚洲avbb在线观看| 国产精品三级大全| 国产白丝娇喘喷水9色精品| 别揉我奶头~嗯~啊~动态视频| 97人妻精品一区二区三区麻豆| 国产综合懂色| 乱人视频在线观看| 在线天堂最新版资源| 日本撒尿小便嘘嘘汇集6| 他把我摸到了高潮在线观看| 免费av观看视频| 国模一区二区三区四区视频| 大型黄色视频在线免费观看| 亚洲图色成人| 欧美日韩中文字幕国产精品一区二区三区| 欧美不卡视频在线免费观看| 18禁黄网站禁片午夜丰满| 热99在线观看视频| 99热精品在线国产| 色播亚洲综合网| 久久久久久九九精品二区国产| 中文字幕免费在线视频6| 精品久久久久久久久av| 国产白丝娇喘喷水9色精品| 国产麻豆成人av免费视频| 婷婷六月久久综合丁香| 欧美高清成人免费视频www| 在线观看午夜福利视频| 人人妻,人人澡人人爽秒播| 日韩欧美国产在线观看| 俄罗斯特黄特色一大片| 午夜福利在线观看免费完整高清在 | 色在线成人网| 男女那种视频在线观看| 嫩草影院新地址| 不卡一级毛片| 午夜福利成人在线免费观看| 成年女人永久免费观看视频| 欧美最新免费一区二区三区| 非洲黑人性xxxx精品又粗又长| www.色视频.com| 在线观看舔阴道视频| 午夜福利视频1000在线观看| 亚洲精品日韩av片在线观看| 国产精品电影一区二区三区| 日本 欧美在线| 午夜免费成人在线视频| 男女那种视频在线观看| 日日干狠狠操夜夜爽| 久久热精品热| 精品久久久久久久久久久久久| 成人特级av手机在线观看| 精品不卡国产一区二区三区| 免费看av在线观看网站| 免费高清视频大片| 亚洲色图av天堂| aaaaa片日本免费| 国产成年人精品一区二区| 久久久久久久久中文| 精品免费久久久久久久清纯| av福利片在线观看| 毛片女人毛片| 中文字幕av成人在线电影| xxxwww97欧美| 男人舔女人下体高潮全视频| 国产精品永久免费网站| 亚洲熟妇中文字幕五十中出| 中文字幕久久专区| 露出奶头的视频| 91狼人影院| 亚洲aⅴ乱码一区二区在线播放| 黄色日韩在线| 国模一区二区三区四区视频| 国产精品久久久久久久电影| 欧美色视频一区免费| 神马国产精品三级电影在线观看| 国产精品伦人一区二区| 国产精品av视频在线免费观看| av黄色大香蕉| 亚洲内射少妇av|