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

    復式耕整機耕深與耕寬穩(wěn)定性分析與試驗

    2016-12-19 08:54:27丁為民方志超杜濤濤趙思琪
    農(nóng)業(yè)工程學報 2016年9期
    關(guān)鍵詞:耕深復式耕作

    秦 寬,丁為民※,方志超,杜濤濤,趙思琪,王 朕

    (1. 南京農(nóng)業(yè)大學工學院,南京 210031; 2. 江蘇省智能化農(nóng)業(yè)裝備重點實驗室,南京 210031)

    復式耕整機耕深與耕寬穩(wěn)定性分析與試驗

    秦 寬1,2,丁為民1,2※,方志超1,2,杜濤濤1,趙思琪1,王 朕1

    (1. 南京農(nóng)業(yè)大學工學院,南京 210031; 2. 江蘇省智能化農(nóng)業(yè)裝備重點實驗室,南京 210031)

    針對設(shè)計的復式耕整機出現(xiàn)的耕作穩(wěn)定性問題,結(jié)合復式耕整機整體結(jié)構(gòu)及工作原理,從牽引、水平面內(nèi)受力、機器振動3個角度分析影響耕作穩(wěn)定性因素,確定影響耕作穩(wěn)定性關(guān)鍵因素為牽引角、犁體配置斜角、旋耕刀升角。以牽引角、犁體配置斜角、旋耕刀升角為試驗因素,以工況耕深穩(wěn)定系數(shù)和工況耕寬穩(wěn)定系數(shù)為性能評價指標進行二次正交旋轉(zhuǎn)組合試驗。正交試驗結(jié)果表明:試驗因素對評價指標影響程度從高到低皆為:犁體配置斜角、牽引角、旋耕升角,當各影響因素分別取值為17.3°、27.8°、72.6°時,工況耕深穩(wěn)定系數(shù)和工況耕寬穩(wěn)定系數(shù)分別為91.8%、93.4%。以影響因素最優(yōu)參數(shù)組合為基礎(chǔ)進行的驗證試驗結(jié)果表明:試驗后工況耕深穩(wěn)定系數(shù)和工況耕寬穩(wěn)定系數(shù)為91.5%、93.1%,與軟件分析結(jié)果基本一致,且其他耕作指標均達到農(nóng)藝要求。該研究可為復式整地機械的耕作穩(wěn)定性研究提供技術(shù)參考。

    農(nóng)業(yè)機械;優(yōu)化;穩(wěn)定性;犁旋作業(yè);復式;耕深;耕寬

    0 引言

    農(nóng)作物秸稈機械化耕整還田既能減少秸稈露天焚燒現(xiàn)象,防止空氣污染,保護環(huán)境,又有利于增加土壤肥力,改善土壤結(jié)構(gòu),是發(fā)展生態(tài)農(nóng)業(yè)和實現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展的重要措施[1-2]。目前,復式作業(yè)機械已成為農(nóng)業(yè)機械發(fā)展趨勢,針對秸稈粉碎還田問題,復式耕整機能夠?qū)崿F(xiàn)鏵式犁翻耕與旋耕機碎壟耕整作業(yè),推動沙土及黏土地區(qū)秸稈還田步伐[3]。復式耕整機作為犁耕、旋耕復式作業(yè)機械,質(zhì)量大、結(jié)構(gòu)復雜、工作阻力大,因此機器作業(yè)過程中易出現(xiàn)側(cè)傾、旋轉(zhuǎn)現(xiàn)象或趨勢,影響機器耕作穩(wěn)定性[4]。

    復式耕整機作業(yè)后,要求耕后溝底與地表平整,耕深穩(wěn)定,耕作層膨松均勻,秸稈粉碎率高,入土均勻性好[5]。其中復式耕整機耕作穩(wěn)定性直接影響整地質(zhì)量,影響因素包括機組結(jié)構(gòu)布置、機具懸掛狀態(tài)、土壤物理特性、機器作業(yè)速度、耕前地表平整度等。因此,如何解決耕作穩(wěn)定性問題成為眾多學者的長期研究目標。Ryan 等[6]研究了美國西北部小麥留茬地整地機器的作業(yè)平衡性問題,提出耕作機械與牽引機械之間距離應不小于1 200 mm。孫松林等[7]從入土過程和正常耕作兩部分對犁體進行受力分析,確定了犁體支持面與地平面夾角對耕作穩(wěn)定性的作用關(guān)系。柳克令[8]對懸掛犁耕作穩(wěn)定性進行力學、運動學分析,提出了瞬心位置對耕深穩(wěn)定性的影響和以此為基礎(chǔ)的調(diào)節(jié)原則。熊一兵[9]從牽引的角度闡述手扶拖拉機組縱向與橫向平衡重心對耕作穩(wěn)定性的影響。但國內(nèi)外對復式整地機械耕作穩(wěn)定性的研究較少。本文圍繞如何提高復式耕整機耕作穩(wěn)定性,重點研究影響耕作穩(wěn)定性的關(guān)鍵因素耕深與耕寬,以期得到提高耕作穩(wěn)定性的最優(yōu)參數(shù)組合,為優(yōu)化復式耕整機復式耕作質(zhì)量提供技術(shù)參考。

    1 總體結(jié)構(gòu)及工作原理

    復式耕整機整體結(jié)構(gòu)如圖1所示,主要由前置犁耕機構(gòu)。后置旋耕機構(gòu)。傳動機構(gòu)、旋耕機架、犁耕機架、懸掛機架等機構(gòu)組成。犁耕機構(gòu)犁鏵數(shù)為4個,犁耕幅寬1 400 mm。旋耕機構(gòu)動力來自于拖拉機尾部輸出,通過萬向軸與減速箱相配合的傳動系統(tǒng)傳輸于旋耕軸,旋耕機構(gòu)采用反轉(zhuǎn)作業(yè)方式,旋耕刀采用彎刀,彎刀按雙頭螺旋線排列,旋向相反,升角相同,刀軸旋轉(zhuǎn)一周內(nèi)總有一把旋耕刀入土,旋耕幅寬1 800 mm。本機器采用三點式懸掛機構(gòu),可通過調(diào)節(jié)拖拉機的三點懸掛機構(gòu)牽引角度與位置,實現(xiàn)耕深調(diào)節(jié)。機器主要參數(shù)為:外形尺寸1 650 mm×1 940 mm×1 140 mm、轉(zhuǎn)彎半徑≤1 700 mm、耕寬1 700~2 000 mm、耕深180~240 mm、拖拉機配套動力≥75 kW。

    工作時,由拖拉機牽引復式耕整機進行犁翻旋耕復式作業(yè),犁耕機構(gòu)在牽引力的作用下首先對土壤進行翻耕作業(yè),旋耕機構(gòu)在驅(qū)動力作用下進而對土壤進行反轉(zhuǎn)旋耕作業(yè)。復式耕整機一次性作業(yè)夠能完成秸稈粉碎、還田、碎土、覆土、平整地表等多道工序,從而減少作業(yè)次數(shù)、提高工作效率。

    圖1 復式耕整機結(jié)構(gòu)圖Fig.1 Structure schematic of plowing and rotary tillage combined machine

    2 影響耕作穩(wěn)定性關(guān)鍵部件設(shè)計

    影響耕作穩(wěn)定性因素多樣而復雜[10],包括懸掛方式、牽引點位置、土壤阻力、機架振動狀態(tài)、機具前進速度等[11-14]。復式耕整機作為復式整地機械,機具質(zhì)量大,結(jié)構(gòu)復雜,受力項多,因此需要確定影響其耕作穩(wěn)定性的關(guān)鍵因素。

    2.1 機器牽引角設(shè)計

    2.1.1 機器質(zhì)心位置的確定

    復式耕整機由前置犁耕機構(gòu)與后置旋耕機構(gòu)組成,旋耕機構(gòu)除了具有旋耕作業(yè)功能外,還具有平衡犁耕機構(gòu)配重作用[15],復式耕整機總質(zhì)量為1 050 kg,后置旋耕機構(gòu)質(zhì)量為700 kg,使機器重心置后,防止發(fā)生翹尾現(xiàn)象[16]。圖2為縱垂面牽引線調(diào)節(jié)示意圖,犁耕機架長度L1為1 250 mm,旋耕機構(gòu)長度L2為700 mm,犁耕機構(gòu)與旋耕機構(gòu)質(zhì)量均勻分布,因此按質(zhì)量分步推算,重心位置O應與機架起始位置距離L3為1 287.5 mm。

    圖2 縱垂面牽引線調(diào)節(jié)示意圖Fig.2 Schematic diagram of draft line adjustment in vertical

    2.1.2 牽引角的確定

    復式耕整機犁耕機構(gòu)在作業(yè)過程中所有阻力包括犁體曲面所受土壤阻力與犁側(cè)板與溝壁摩擦力[17],2種作用力在縱垂面內(nèi)分力之和用Rxz表示,如圖2所示,由于犁體在縱垂面內(nèi)從前往后依次等距排開,因此將力Rxz對犁體作用點簡化至四犁體長度幾何中點A點;旋耕機構(gòu)在作業(yè)過程所受阻力為旋耕刀切削土壤阻力[18],此作用力在縱垂面內(nèi)分力用Pxz表示,Pxz對旋耕刀作用點簡化至旋耕刀端點B。過A點與B點連線與過質(zhì)心垂線相交與點C,牽引線π1e一端通過機器縱垂面內(nèi)瞬心π1點,另一端點e點希望與C點相接近,e點若相對于C點向上偏移距離過大,牽引線上移,牽引角α(牽引線與水平線夾角)減小,則可能造成犁耕深度加深,旋耕機構(gòu)上翹,整機動力消耗增加[19];e點若相對于C點向下偏移距離過大,牽引線下移,牽引角α增大,則可能造成犁耕機構(gòu)前端上翹,致使犁耕機構(gòu)入土變淺,旋耕機構(gòu)入土過深[20]。因此,為保證機器在縱垂面內(nèi)的耕作穩(wěn)定性,牽引角α的取值范圍應在15°~30°。

    2.2 犁體配置斜角設(shè)計

    復式耕整機關(guān)鍵結(jié)構(gòu)直接影響機器受力平衡,從而影響機器耕作穩(wěn)定性。復式耕整機在作業(yè)時,會受到拖拉機牽引力、土壤阻力及自身重力作用[21]。在牽引力滿足工作需求的前提下,機器在縱垂面內(nèi)受力基本可以達到平衡狀態(tài),而機器在水平面內(nèi)受力由于無土壤約束,自由度大,容易由于受力不均產(chǎn)生傾斜、旋轉(zhuǎn)現(xiàn)象,對耕作穩(wěn)定性造成影響[22],因此設(shè)計機具結(jié)構(gòu)時重點考慮水平面內(nèi)受力平衡。如圖3所示,為機器水平面受力投影,四犁體所受土壤阻力在水平面內(nèi)投影為Rxy,犁體所受土壤阻力以力場形式作用于犁體曲面[23],因此將犁體曲面所受土壤阻力簡化至犁體曲面中心,垂直于犁體曲面中心切線且與滾垡方向相反;四犁體犁側(cè)板與溝壁摩擦力在水平面內(nèi)投影為Qxy,作用在犁側(cè)板中心點上,沿犁側(cè)板方向且與犁體前進方向相反,由于四犁體所受土壤阻力、犁側(cè)板與溝壁板摩擦力基本相同,因此將Rxy、Qxy簡化至犁體中心點連線的中點E。牽引力在水平面內(nèi)投影為Fxy,作用在犁耕機架前端中心點D。旋耕刀削土阻力在水平面內(nèi)投影為Pxz,由于旋耕刀左、右對稱排列,因此在y軸方向阻力基本可以相互抵消,旋耕刀削土阻力方向與x軸一致,作用在旋耕機架中心位置的旋耕刀端點M。根據(jù)平面力系平衡條件,在y軸上列平衡方程,得

    對質(zhì)心位置取矩得

    式中Fy為機具水平面內(nèi)y軸方向受力,N;Fxy為牽引力水平面內(nèi)投影,N;Rxy為犁體所受土壤阻力水平面內(nèi)投影,N;MO2為水平面內(nèi)質(zhì)心位置力矩,N·m;Qxy為犁側(cè)板與溝壁摩擦力在水平面內(nèi)投影,N;Pxy為旋耕刀削土阻力在水平面內(nèi)投影,N;x1為E點至y軸距離,mm;y1為E點至x軸距離,mm;y2為D點至x軸距離,mm;y3為M點至x軸距離,mm;β為犁體配置斜角,(°)。

    由式(1)、式(2)可知,犁體配置斜角β的角度直接影響機器橫向受力與力偶平衡,從而影響機器的耕作穩(wěn)定性,因此為保證耕作穩(wěn)定性,β范圍應在23°~30°。

    圖3 機器水平面受力投影Fig.3 Force projection of machine in level surface

    2.3 旋耕刀升角設(shè)計

    復式作耕整作業(yè)時會受到外部載荷激勵與自身系統(tǒng)內(nèi)部激勵作用,使整機在作業(yè)過程中出現(xiàn)振動現(xiàn)象,從而影響整個機器的耕作穩(wěn)定性。復式耕整機犁耕機構(gòu)與旋耕作業(yè)時都會產(chǎn)生振動現(xiàn)象,犁耕機構(gòu)振動主要由犁體所受土壤阻力引起,旋耕機構(gòu)振動主要由旋耕刀所受土壤阻力與萬向軸、主減速器、副減速器等傳動機構(gòu)嚙合高速旋轉(zhuǎn)引起,此外復式耕整機振動還受到拖拉機機體振動的影晌[24]。在上述所有引起振動的激勵因素中,旋耕刀所受土壤阻力是引起機器振動的最主要激勵。

    圖4為土壤對旋耕刀作用力示意圖,旋耕刀所受土壤阻力Q可用三分力表示,如式(3)所示。

    式中Q為旋耕刀所受土壤阻力,N;i,j,k為單位矢量;Hs,Rv,Ph分別為旋耕刀所受土壤阻力的橫向水平分力(x軸)、縱向水平分力(y軸)和垂直水平分力(z軸)。復式耕整機旋耕刀安裝時繞刀輥以雙頭螺旋線形式左右彎刀交替排列,旋向相反,升角相同,此時旋耕刀作業(yè)瞬時橫向水平力

    式中n1,n2分別為該瞬時同時工作的左、右彎刀片數(shù);Hsi為單個旋耕刀所受土壤阻力橫向水平分力。由于旋耕刀左、右彎刀對稱排列,因此橫向水平力總體可相互抵消,但旋耕刀片是沿螺旋線等角度(升角)間隔排列,各個瞬間交變負荷總是存在的,這種交變負荷是引起機組振動的主要激勵[25]。當升角過大,則單位工作周期內(nèi)旋耕刀入土次數(shù)少,容易對機器產(chǎn)生低頻激勵,而帶有旋耕機構(gòu)的復式耕作機械對低頻激勵比較敏感,容易產(chǎn)生較大振動;升角過小,則作業(yè)時容易夾土堵塞[26],因此同一螺旋線上相鄰旋耕刀升角范圍為54°~85°。

    圖4 土壤對旋耕刀作用力示意圖Fig.4 Acting force diagram of soil to rotary blade

    3 試驗及結(jié)果分析

    3.1 試驗設(shè)備

    試驗機具為委托南通世創(chuàng)公司加工生產(chǎn)的1LFG-140型犁翻旋耕復式作業(yè)耕整機。試驗器材包括:皮尺(1~100 m,0.001m)、直尺(1~600 mm,1 mm)、水平儀(上海鈺誠電子有限公司,iLevel 5,精度:0.029°)、耕深尺(精度:0.01 mm)、土壤堅實度儀(浙江托普儀器有限公司,TJSD-750,±0.5‰FS)、水分測試儀(上海婉源電子科技有限公司,SK-100,0.01%)。

    3.2 試驗方法

    2015年6月14日在江蘇省常州市金壇區(qū)沙湖村對復式耕整機穩(wěn)定性進行試驗,平均作業(yè)時速為1.43 m/s,配套動力為久保田854。試驗田塊特性參數(shù)如表1所示,試驗現(xiàn)場如圖5所示。

    表1 試驗田塊特性參數(shù)Table 1 Characteristics of experimental field environment

    圖5 試驗現(xiàn)場Fig.5 Test of field

    試驗方法參照GB/T14225-2008《鏵式犁 試驗方法》與GB/T5668-2008《旋耕機械 試驗方法》。具體方案如下。

    1)耕深及工況耕深穩(wěn)定系數(shù)試驗方法。沿機組前進方向每隔2 m左、右兩側(cè)各取一點,每個行程測量20個點,用耕深尺測量每個測量點耕深度,共測3個行程。耕深按照式(5)、式(6)計算;工況耕深穩(wěn)定系數(shù)按式(5)-式(10)計算,計算時,一個行程中左右兩測量點各算一個單獨行程。

    式中aj為第j個行程的耕深平均值,cm;aji為第j個行程中的第i個點的耕深值,cm;nj為第j個行程中的測定點數(shù);U為工況的耕深穩(wěn)定系數(shù),%;a為工況耕深平均值,cm;N為同一工況中的行程數(shù);Sj為第j個行程耕深標準差,cm;S為工況的耕深標準差,cm;V為工況的耕深變異系數(shù),%。

    2)耕寬及工況耕寬穩(wěn)定系數(shù)試驗方法。沿垂直機組運動方向測定2個相鄰行程溝墻之間的水平距離即為耕寬,在測定耕深的相應處進行測量耕寬,分別計算工況耕寬和工況耕寬穩(wěn)定性系數(shù),計算方法同耕深。

    3)地表平整度試驗方法。沿垂直于機組前進方向,在地表最高點為基準取一長度為機器幅寬的水平基準線,均分成10等分,測定各等分點至地表的距離,一個行程測3組,共測3個行程。按式(4)計算其標準差,并以標準差的平均值表示平整度。

    4)碎土率試驗方法。在已耕地上測定0.5 m×0.5 m面積內(nèi)的全耕層土壤,土塊大小按其最長邊分為小于4 cm、4~8 cm和大于8 cm三級。并以小于4 cm的土塊質(zhì)量占總質(zhì)量的百分比為碎土率,每一行程測定一點,共測3個行程。

    5)植被覆蓋率試驗方法。在測區(qū)內(nèi)對角線上取5點,每點按1 m2面積貼地面剪下露出地表的植被,稱其質(zhì)量,并計算出5點的平均值,每個行程測量1點,共測3個行程,按式(11)計算植被覆蓋率。式中Fb為植被覆蓋率,%;Wq為耕前植被平均值,g;Wh為耕后植被平均值,g。

    3.3 試驗設(shè)計

    為探究復式耕整機關(guān)鍵參數(shù)對耕作穩(wěn)定性的影響,進行三因素二水平二次正交旋轉(zhuǎn)組合試驗,試驗因素為牽引角α、犁體配置斜角β、旋耕刀升角γ,復式耕整機作業(yè)后耕深穩(wěn)定性與耕寬穩(wěn)定性最能夠體現(xiàn)機器耕作穩(wěn)定性情況,因此選擇工況耕深穩(wěn)定系數(shù)R1與工況耕寬穩(wěn)定系數(shù)R2作為試驗指標,試驗后用Design-Expert軟件對數(shù)據(jù)進行處理,建立回歸方程與優(yōu)化模型,得到試驗因素對試驗指標影響的主次關(guān)系與最優(yōu)組合,設(shè)計因素水平編碼表如表2所示。

    表2 因素水平編碼表Table 2 Coding with factors and levels

    3.4 多因素試驗結(jié)果與分析

    根據(jù)二次正交旋轉(zhuǎn)組合試驗進行試驗,以α,β,γ取值為試驗因素,以單幅寬左右耕深差R1、工況耕深穩(wěn)定系數(shù)R2為相應試驗指標,試驗結(jié)果如表3所示。

    表3 試驗方案和試驗結(jié)果Table 3 Protocols and results

    3.4.1 方差分析

    對回歸模型中各項回歸系數(shù)進行F檢驗和方差分析,工況耕深穩(wěn)定系數(shù)R1與工況耕寬穩(wěn)定系數(shù)R2方差分析結(jié)果如表4所示。

    表4 工況耕深穩(wěn)定系數(shù)與工況耕寬穩(wěn)定系數(shù)的二次項模型方差分析Table 4 Anova of quadratic model for working condition tillage depth stability factor and working condition tillage width stability factor

    對表4中的數(shù)據(jù)進行二次多元回歸擬合,選用二次項模型建立工況耕深穩(wěn)定系數(shù)R1、工況耕寬穩(wěn)定系數(shù)R2與各個影響因素之間的回歸模型,去除其中不顯著項后,得到R1、R2對牽引角、犁體配置斜角、旋耕刀升角的二次多元回歸方程為式(12)、式(13)。

    式中R1為工況耕深穩(wěn)定系數(shù),%;R2為工況耕寬穩(wěn)定系數(shù),%;α為牽引角,(°);β為犁體配置斜角,(°);γ為旋耕刀升角,(°)。

    由表4可知,目標函數(shù)R1、R2的模型失擬項P值分別為0.4619、0.512,均大于0.05,說明無失擬因素存在,可以用上述回歸方程代替試驗真實點對試驗結(jié)果進行分析。

    由表4方差分析可知,工況耕深穩(wěn)定系數(shù)R1、工況耕寬穩(wěn)定系數(shù)R2模型顯著性P值分別為0.0092與0.0083,均小于0.05,說明該模型具有統(tǒng)計學意義。對于目標函數(shù)R1,因素β、γ、αβ、β2、γ2非常顯著,因素βγ顯著;對于目標函數(shù)R2,因素β、γ、βγ非常顯著,因素α、αβ、β2、γ2顯著。表中F值表示各個影響因素對試驗指標的影響,F(xiàn)值越大對試驗指標影響越大,由表4可知,各個試驗因素對工況耕深穩(wěn)定系數(shù)R1、工況耕寬穩(wěn)定系數(shù)R2的影響程度從大到小皆依次是:犁體配置斜角β、牽引角α、旋耕刀升角γ。

    3.4.2 響應曲面分析

    固定對試驗指標影響最小的旋耕刀升角γ為零水平,令其值為69.5°,考察牽引角α、犁體配置斜角β兩因素對工況耕深穩(wěn)定系數(shù)與工況耕寬穩(wěn)定系數(shù)的影響規(guī)律,得到對應的響應曲面,如圖6所示。當旋耕刀升角γ為69.5°,牽引角α為23.7°,犁體配置斜角β為26.8°時工況耕深穩(wěn)定系數(shù)存在一個最優(yōu)值,為91.7%;當旋耕刀升角γ為69.5°,牽引角α為25.2°,犁體配置斜角β為27.3°時工況耕寬穩(wěn)定系數(shù)存在一個最優(yōu)值,為93.2%。

    3.4.3 最佳參數(shù)組合的確定

    通過二次正交旋轉(zhuǎn)組合試驗,需要確定影響耕作穩(wěn)定性因素參數(shù)的最優(yōu)組合,從而提高工況耕深穩(wěn)定系數(shù)、工況耕寬穩(wěn)定系數(shù)等耕作穩(wěn)定性指標,達到提高耕作穩(wěn)定性的目的。根據(jù)農(nóng)藝要求,保證工況耕寬穩(wěn)定系數(shù)在90%以上,工況耕深穩(wěn)定系數(shù)越大越好。因此確定目標函數(shù)為式(14),約束函數(shù)為式(15)。

    式中F(α,β,γ)為目標函數(shù);s.t.G(α,β,γ)為約束函數(shù)。

    根據(jù)目標函數(shù)與約束函數(shù)模型,利用Design-Expert軟件對回歸方程式(12)、式(13)進行優(yōu)化,得到影響復式耕整機工況耕深穩(wěn)定系數(shù)、工況耕寬穩(wěn)定系數(shù)因素的最優(yōu)參數(shù)組合。當影響因素參數(shù)組合為牽引角α為17.3°、犁體配置角β為27.8°、旋耕升角γ為72.6°,此時工況耕深穩(wěn)定系數(shù)為91.8%,工況耕寬穩(wěn)定系數(shù)為93.4%,此組為最優(yōu)參數(shù)組合。

    圖6 牽引角和犁體配置斜角對工況耕深穩(wěn)定系數(shù)與工況耕寬穩(wěn)定系數(shù)的響應曲面Fig.6 Response surface showing effects of angle of traction between installing bevel angle of plough to working condition tillage depth stability factor and working condition tillage width stability factor

    3.5 驗證試驗

    根據(jù)優(yōu)化分析得到的最優(yōu)參數(shù)組合,于2015年6月16日在江蘇省常州市金壇區(qū)沙湖村進行驗證試驗。根據(jù)二次正交旋轉(zhuǎn)組合試驗后的優(yōu)化結(jié)果,對機器選取牽引角α為17.3°、犁體配置角β為27.8°、旋耕升角γ為72.6°的最優(yōu)參數(shù)組合,進行驗證試驗。試驗指標除工況耕深穩(wěn)定系數(shù)與工況耕寬穩(wěn)定系數(shù)外,同時考察耕深、耕寬、地表平整度、碎土率、秸稈覆蓋率5項耕作指標,以檢驗分析所得最優(yōu)參數(shù)組合對機器耕作穩(wěn)定性的全面影響。共進行5次重復試驗,試驗后取平均值,試驗結(jié)構(gòu)計算方法見3.2節(jié)。對試驗數(shù)據(jù)進行處理與分析,驗證試驗結(jié)果如表5所示。

    表5 驗證試驗結(jié)果Table 5 Test result of verification

    驗證試驗得到工況耕深穩(wěn)定系數(shù)為91.5%,工況耕寬穩(wěn)定性系數(shù)為93.1%,與軟件分析得到的工況耕深穩(wěn)定系數(shù)91.8%,工況耕寬穩(wěn)定系數(shù)93.4%,相差較小,說明軟件優(yōu)化參數(shù)具有準確性與可行性。復式耕整機其他耕作指標均達到農(nóng)藝要求,說明最優(yōu)參數(shù)組合下機器的耕作質(zhì)量可以達到播種前整地的農(nóng)業(yè)標準。

    4 結(jié)論與討論

    1)本文對復式耕整機耕作穩(wěn)定性關(guān)鍵部件進行研究,從牽引、水平面內(nèi)受力、機組振動3個角度分析出影響耕作穩(wěn)定性的3個關(guān)鍵因素為牽引角、犁體配置斜角、旋耕刀升角。

    2)通過三因素二水平二次正交旋轉(zhuǎn)組合試驗,得出各影響因素對工況耕深穩(wěn)定系數(shù)和工況耕寬穩(wěn)定系數(shù)影響程度從高到低皆為:犁體配置斜角、牽引角、旋耕刀升角。優(yōu)化分析得出牽引角為17.3°、犁體配置角為27.8°、旋耕升角為72.6°為最優(yōu)參數(shù)組合,此時,工況耕深穩(wěn)定系數(shù)為91.8%,工況耕寬穩(wěn)定系數(shù)為93.4%。

    3)驗證試驗表明,試驗后工況耕深穩(wěn)定系數(shù)和工況耕寬穩(wěn)定系數(shù)分別為91.5%與93.1%,與軟件分析結(jié)果基本一致,且復式耕整機其他耕作標均達到農(nóng)藝要求。

    本文從牽引、機器水平面內(nèi)受力、振動3個關(guān)鍵方向研究復式耕整機耕作穩(wěn)定性,而實際作業(yè)過程中,機器前進速度、機器縱垂面內(nèi)受力,土壤環(huán)境等其他多種因素均會對耕作穩(wěn)定性產(chǎn)生影響,因此對于復式耕整機耕作穩(wěn)定性更全面的研究,有待于進一步展開。

    [1] 趙亞麗,薛志偉,郭海斌,等. 耕作方式與秸稈還田對土壤呼吸的影響及機理[J]. 農(nóng)業(yè)工程學報,2014,30(19):155-165. Zhao Yali, Xue Zhiwei, Guo Haibin, et al. Effects of tillage and crop residue management on soil respiration and its mechanism[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(19): 155-165. (in Chinese with English abstract)

    [2] 劉世平,聶新濤,張洪程,等. 稻麥兩熟條件下不同土壤耕作方式與秸稈還田效用分析[J]. 農(nóng)業(yè)工程學報,2006,22(7):48-51. Liu Shiping, Nie Xintao, Zhang Hongcheng, et al. Effects of tillage and straw returning on soil fertility and grain yield in a wheat-rice double cropping system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2006, 22(7): 48-51. (in Chinese with English abstract)

    [3] 沈丹波,繆明,丁煒. 稻麥秸稈犁翻旋耕復式作業(yè)耕整機的研制[J]. 農(nóng)業(yè)裝備技術(shù),2014,40(4):14-16.

    [4] 張青松,汲文峰,廖宜濤,等. 油菜直播機鏵式開畦溝前犁曲面分析與阻力特性試驗[J]. 農(nóng)業(yè)機械學報,2014,45(2):130-135. Zhang Qingsong, Ji Wenfeng, Liao Yitao, et al. Surface analysis and resistance characteristics experiment on ditch Plow ahead of direct rapeseed seeder[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(2): 130-135. (in Chinese with English abstract)

    [5] 汲文峰,賈洪雷,佟金. 旋耕-碎茬仿生刀片田間作業(yè)性能的試驗研究[J]. 農(nóng)業(yè)工程學報,2012,28(12):24-30. Ji Wenfeng, Jia Honglei, Tong Jin. Experiment on working performance of bionic blade for soil-rototilling and stubble-breaking[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(12): 24-30. (in Chinese with English abstract)

    [6] Ryan W H, Stephen S J, Arron H C. Wheat cultivar performance and stability between no-till and conventional tillage systems in the pacific Northwest of the United States[J]. Sustainability, 2013, 5: 882-895.

    [7] 孫松林,吳明亮,謝方平,等. 水田耕整機犁耕作業(yè)穩(wěn)定性分析[J]. 農(nóng)業(yè)機械學報,2007,38(12):217-221.

    [8] 熊一兵. 手扶拖拉機耕深穩(wěn)定性分析[J]. 農(nóng)業(yè)機械學報,1966,9(3):244-247.

    [9] 柳克令. 懸掛犁耕深穩(wěn)定性研究的初步評述[J]. 農(nóng)業(yè)機械學報,1966,9(3):275-281.

    [10] 張欣悅,趙大勇,許春林,等. 1GMMZ-280/4型壟作組合式滅茬旋耕整地機[J]. 農(nóng)機化研究,2012,34(6):23-26. Zhang Xinyue, Zhao Dayong, Xu Chunlin, et al. The development of 1GMMZ-280/4 type ridge culture combined-type stubble ploughing, spin tillage soil preparation machine[J]. Journal of Agricultural Mechanization Research, 2012, 34(6): 23-26. (in Chinese with English abstract)

    [11] 張居敏,周勇,夏俊芳,等. 旋耕埋草機螺旋橫刀的數(shù)學建模與參數(shù)分析[J]. 農(nóng)業(yè)工程學報,2013,29(1):18-25. Zhang Jumin, Zhou Yong, Xia Junfang, et al. Mathematical modeling and analysis of helical blade for stubble burying rotary tiller[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(1): 18-25. (in Chinese with English abstract)

    [12] Onwualu A P,Watts K C. Draught and vertical forces obtained from dynamic soil cutting by plane tillage tools[J]. Soil & Tillage Research, 1998, 48(4): 239-253.

    [13] 張居敏,賀小偉,夏俊芳,等. 高茬秸稈還田耕整機功耗檢測系統(tǒng)設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報,2014,30(18):38-46. Zhang Jumin, He Xiaowei, Xia Junfang, et al. Design and field experiment of power consumption measurement system for high stubble returning and tillage machine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(18): 38-46. (in Chinese with English abstract)

    [14] 趙大勇,李連豪,許春林,等. 1ZQHF-350/5型前后分置懸掛式聯(lián)合整地機[J]. 農(nóng)業(yè)機械學報,2014,45(10):91-96. Zhao Dayong, Li Lianhao, Xu Chunlin, et al. 1ZQHF-350/5 hang combined cultivating machine with front-stubblebreaking, post-subsoil and rotary-tilling equipment[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(10): 91-96. (in Chinese with English abstract)

    [15] 張秀梅,張居敏,夏俊芳,等. 水旱兩用秸稈還田耕整機關(guān)鍵部件設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報,2015,31(11):10-16. Zhang Xiumei, Zhang Jumin, Xia Junfang, et al. Design and experiment on critical component of cultivator for straw returning in paddy field and dry land[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(11): 10-16. (in Chinese with English abstract)

    [16] 車剛,張偉,萬霖,等. 基于滅茬圓盤驅(qū)動旋耕刀多功能耕整機設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報,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 the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(20): 34-40. (in Chinese with English abstract)

    [17] 車剛,張偉,梁遠,等. 3ZFC-7 型全方位復式中耕機的設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報,2011,27(1):130-135. Che Gang, Zhang Wei, Liang Yuan, et al. Design and experiment of the 3ZFC-7 omni-bearing duplex type cultivator[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(1): 130-135. (in Chinese with English abstract)

    [18] 羅金海,孫佳民,楊莉,等. 9BQS-3.0型氣吹式松土播種復式作業(yè)機設(shè)計與試驗[J]. 農(nóng)業(yè)機械學報,2013,44(增刊1):52-66. Luo Jinhai, Sun Jiamin, Yang Li, et al. Design and experiment of type 9BQS-3.0 pneumatic scarifying and sowing compound operation machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(Supp.1): 52-66. (in Chinese with English abstract)

    [19] 蔣金琳,龔麗農(nóng),王東偉,等. 免耕播種機雙刀盤有支撐切茬破茬裝置的研制與試驗[J]. 農(nóng)業(yè)工程學報,2012,28(21):17-22. Jiang Jinlin, Gong Linong, Wang Dongwei, et al. Development and experiment for driving double coulters anti-blockage device of no-till planter[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(21): 17-22. (in Chinese with English abstract)

    [20] 謝斌,李皓,朱忠祥,等. 基于傾角傳感器的拖拉機懸掛機組耕深自動測量方法[J]. 農(nóng)業(yè)工程學報,2013,29(4):15-21. Xie Bin, Li Hao, Zhu Zhongxiang, et al. Measuring tillage depth for tractor implement automatic using inclinometer[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(4): 15-21. (in Chinese with English abstract)

    [21] 邱進,吳明亮,官春云,等. 動定刀同軸水稻秸稈切碎還田裝置結(jié)構(gòu)設(shè)計與試驗[J]. 農(nóng)業(yè)工程學報,2015,31(10):11-19. Qiu Jin, Wu Mingliang, Guan Chunyun, et al. Design and experiment of chopping device with dynamic fixed knife coaxial for rice straw[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(10): 11-19. (in Chinese with English abstract)

    [22] 陳玉侖,丁為民,方志超,等. 全喂入式聯(lián)合收割機碎草脫粒裝置的改進設(shè)計[J]. 農(nóng)業(yè)工程學報,2013,29(24):11-16. Chen Yulun, Ding Weimin, Fang Zhichao, et al. Improved design of straw-cutting type threshing mechanism of full-feeding combine harvester[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(24): 11-16. (in Chinese with English abstract)

    [23] Radite P, Hermawan W, Mulyana F, et al. Experimental investigation on the application of vibration to reduce draft requirement of subsoiler[J]. International Agricultural Engineering Journal, 2010, 19(1): 31-38.

    [24] 賈洪雷,陳忠亮,郭紅,等. 旋耕碎茬工作機理研究和通用刀輥的設(shè)計[J]. 農(nóng)業(yè)機械學報,2000,31(4):29-32. Jia Honglei, Chen Zhongliang, Guo Hong, et al. Study on working principle of rotary tillage and stubble cutting and design of universal knife roller[J]. Transactions of the Chinese Society for Agricultural Machinery, 2000, 31(4): 29-32. (in Chinese with English abstract)

    [25] Siemens M C, Wilkins D E, Correa R F. Development and evaluation of a residue management wheel for hoe-type no-till drills[J]. Transactions of the ASAE, 2004, 47(2): 397-404.

    [26] 余水生. 水田高茬秸稈還田耕整機的研制[D]. 武漢:華中農(nóng)業(yè)大學,2012. Yu Shuisheng. Research and Trial-manufacture of Straw Returning and Tillage Machine for High Stubble Paddy Field[D]. Wuhan: Huazhong Agricultural University, 2012. (in Chinese with English abstract)

    Analysis and experiment of tillage depth and width stability for plowing and rotary tillage combined machine

    Qin Kuan1,2, Ding Weimin1,2※, Fang Zhichao1,2, Du Taotao1, Zhao Siqi1, Wang Zhen1
    (1. College of Engineering, Nanjing Agricultural University, Nanjing 210031, China; 2. Jiangsu Key Laboratory for Intelligent Agricultural Equipment, Nanjing 210031, China)

    The plowing and rotary tillage combined machine developed in this study is a duplex-operation scarification machine designed with the combination of plough mechanism and rotary tillage mechanism in the front. Such a design allows the cultivator to accomplish multiple tasks simultaneously, such as plow tillage, rotary tillage, straw chopping for mulching, soil pulverization, soil covering, and surface leveling. To investigate the stability of the plowing and rotary tillage combined machine, the factors such as cultivator tillage stability were examined from 3 different aspects: traction, force analysis in the horizontal plane, and vibration analysis. The center of mass of the machine could be determined on the longitudinal vertical plane of the machine, where the traction line passing through the instantaneous center of rotation intersected with the vertical line passing through the center of mass. The traction angle (15°-30°), i.e., the angle between traction line and horizontal line , was a crucial factor affecting tillage stability. The forces acting on the machine on the horizontal plane during operation were analyzed to obtain the plow tilt angle (23°-30°) from the equilibrium equations of the plough, rotary blade, and traction forces on the horizontal plane; the plow tilt angle affected tillage stability by directly impacting the force balance on the horizontal plane. From an analysis of the vibrational excitation during machine operation, it was determined that alternating load in rotary blade operations was the main source of machine vibration excitation. The lift angle of the rotary blade (54°-85°) affected the alternating load, therefore influencing tillage stability. The traction angle, plow tilt angle, and lift angle of the rotary blade were used as the experimental variables, and the stability coefficient under tillage depth-based working conditions and the stability coefficient under tillage width-based working conditions were used as the experimental indicators in a 3-factor/2-level quadratic orthogonal rotating combinatorial test to determine the optimal parameter combination of the influencing factors. Subsequently, the regression equations with the stability coefficient as the objective functions as well as the response surface for the stability coefficient could be obtained by analyzing the test results using the Design-Expert software. The variance analysis showed that among the 3 variables in the test, the plow tilt angle had the greatest influence on the tillage depth and tillage width stability coefficients, while the lift angle of the rotary blade had the least influence on the coefficients. The response surface analysis showed that with a fixed lift angle of the rotary blade of 69.5°, an optimal tillage depth stability coefficient of 91.7% could be obtained with traction angle and plow tilt angle of 23.7° and 26.8°, respectively. An optimal tillage width stability coefficient of 93.2% could be achieved with traction angle and plow tilt angle of 25.2° and 27.3°, respectively. The optimal parameter combination was traction angle of 17.3°, plow tilt angle of 27.8°, and rotary blade’s lift angle of 72.6°. This optimal combination could achieve the working condition tillage depth stability coefficient of 91.8% and the tillage width stability coefficient of 93.4%, respectively. The validation experiments showed that with the optimal parameter combination of the influencing factors, the working condition tillage depth stability coefficient and the tillage width stability coefficient were 91.5% and 93.1%, respectively; these results were consistent with the ones obtained via software analysis. Other tillage performance indicators of the machine, such as tillage depth, tillage width, surface leveling degree, pulverization rate, and straw coverage rate which were respectively 1.87, 1.98, 21.20, 90.30% and 90.70%, all met the agronomic requirements.

    agricultural machinery; optimization; stability; plough and rotary tillage; recombination; tilling depth; tilling width

    10.11975/j.issn.1002-6819.2016.09.001

    S222.4

    A

    1002-6819(2016)-09-0001-08

    秦 寬,丁為民,方志超,杜濤濤,趙思琪,王 朕.復式耕整機耕深與耕寬穩(wěn)定性分析與試驗[J]. 農(nóng)業(yè)工程學報,2016,32(9):1-8.

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

    Qin Kuan, Ding Weimin, Fang Zhichao, Du Taotao, Zhao Siqi, Wang Zhen. 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) doi:10.11975/j.issn.1002-6819.2016.09.001 http://www.tcsae.org

    2015-10-25

    2016-03-05

    國家科技支撐計劃項目資助(2013BAD08B04)

    秦 寬,男,安徽蚌埠人,博士生,主要研究方向為農(nóng)業(yè)機械化裝備研究。南京 南京農(nóng)業(yè)大學工學院,210031。Email:qinkuan_njau@163.com

    ※通信作者:丁為民,男,安徽合肥人,教授,博士生導師,主要從事農(nóng)業(yè)機械化裝備研究。南京 南京農(nóng)業(yè)大學工學院,20031。Email:wmding@njau.edu.cn

    猜你喜歡
    耕深復式耕作
    拖拉機多重模糊PID變論域耕深調(diào)節(jié)研究
    各式各樣的復式條形統(tǒng)計圖
    基于卡爾曼濾波融合算法的深松耕深檢測裝置研究
    略談“垂直互動”復式教學的作用
    甘肅教育(2020年4期)2020-09-11 07:42:04
    懸掛式深松機耕深監(jiān)測系統(tǒng)的設(shè)計與試驗
    價值工程(2019年31期)2019-12-04 04:11:05
    澳門復式住宅
    線性擬合與Kalman預測法修正耕深測量誤差
    耕作深度對紫色土坡地旋耕機耕作侵蝕的影響
    玉米保護性耕作的技術(shù)要領(lǐng)
    草地耕作技術(shù)在澳大利亞的應用
    土壤與作物(2015年3期)2015-12-08 00:46:58
    国产欧美另类精品又又久久亚洲欧美| 美女主播在线视频| 久久久欧美国产精品| 97精品久久久久久久久久精品| 国产精品三级大全| 人妻夜夜爽99麻豆av| 国产av一区二区精品久久| 在线亚洲精品国产二区图片欧美 | 日本黄大片高清| 精品人妻一区二区三区麻豆| 国产在视频线精品| 搡老乐熟女国产| 这个男人来自地球电影免费观看 | 欧美另类一区| 黄片无遮挡物在线观看| av福利片在线观看| 成年人免费黄色播放视频 | 免费少妇av软件| 久久精品久久久久久久性| 一区二区三区免费毛片| 欧美精品国产亚洲| freevideosex欧美| 久久婷婷青草| 超碰97精品在线观看| 久久久国产一区二区| 国产精品一区二区在线观看99| 高清在线视频一区二区三区| av女优亚洲男人天堂| 人人澡人人妻人| 亚洲伊人久久精品综合| 22中文网久久字幕| 国产精品一区二区三区四区免费观看| 又爽又黄a免费视频| 国产69精品久久久久777片| 秋霞在线观看毛片| 一本久久精品| av在线播放精品| 亚洲情色 制服丝袜| av一本久久久久| 国产精品人妻久久久影院| 又大又黄又爽视频免费| 一区二区三区精品91| 高清午夜精品一区二区三区| 在线精品无人区一区二区三| av黄色大香蕉| 久热久热在线精品观看| 少妇熟女欧美另类| 最近中文字幕高清免费大全6| 亚洲精品色激情综合| 美女国产视频在线观看| 在线观看免费日韩欧美大片 | 中文资源天堂在线| 国产日韩欧美在线精品| 亚洲精品国产av蜜桃| 日韩免费高清中文字幕av| 日日爽夜夜爽网站| 成人影院久久| 在线观看三级黄色| 欧美bdsm另类| 国产永久视频网站| 精品国产一区二区三区久久久樱花| 99久久精品一区二区三区| 精品人妻一区二区三区麻豆| 成人二区视频| 精品一区二区三区视频在线| 亚洲高清免费不卡视频| 久久精品久久精品一区二区三区| 少妇人妻久久综合中文| 中文欧美无线码| av有码第一页| 伊人久久国产一区二区| 久热久热在线精品观看| 国产精品国产三级国产专区5o| 久久青草综合色| 亚洲精品国产色婷婷电影| 一级毛片 在线播放| av在线观看视频网站免费| 日韩av不卡免费在线播放| 你懂的网址亚洲精品在线观看| 日韩成人伦理影院| 欧美精品高潮呻吟av久久| 日本黄色片子视频| 国产中年淑女户外野战色| 国产精品一区二区在线不卡| 亚洲av.av天堂| 91aial.com中文字幕在线观看| 美女大奶头黄色视频| 亚洲国产精品一区二区三区在线| 成人免费观看视频高清| 日韩在线高清观看一区二区三区| tube8黄色片| 极品少妇高潮喷水抽搐| 中文字幕制服av| 51国产日韩欧美| 亚洲国产av新网站| 人体艺术视频欧美日本| 国产黄片美女视频| 在线精品无人区一区二区三| 亚洲性久久影院| 久久久久国产网址| 在线精品无人区一区二区三| 久久久精品免费免费高清| 日韩精品免费视频一区二区三区 | 国产91av在线免费观看| 韩国av在线不卡| 热re99久久国产66热| 免费大片18禁| 欧美精品国产亚洲| 街头女战士在线观看网站| 永久网站在线| 国产高清不卡午夜福利| 伦精品一区二区三区| 免费看不卡的av| 国产成人精品一,二区| 黑人巨大精品欧美一区二区蜜桃 | 欧美三级亚洲精品| 边亲边吃奶的免费视频| 久久国产精品大桥未久av | 免费不卡的大黄色大毛片视频在线观看| 一区二区av电影网| 成人国产麻豆网| 免费观看a级毛片全部| av在线观看视频网站免费| 亚洲精品日韩在线中文字幕| 美女大奶头黄色视频| 少妇高潮的动态图| 亚洲,欧美,日韩| 人人妻人人看人人澡| 国产色婷婷99| 精品亚洲成国产av| 国产色爽女视频免费观看| av线在线观看网站| 国产乱人偷精品视频| tube8黄色片| 大又大粗又爽又黄少妇毛片口| 精品一品国产午夜福利视频| 一级毛片aaaaaa免费看小| av在线app专区| 午夜影院在线不卡| 久久久a久久爽久久v久久| 亚洲精品乱码久久久v下载方式| 国产91av在线免费观看| 久久久久久久大尺度免费视频| 久久精品国产亚洲网站| 欧美日韩av久久| 香蕉精品网在线| 久久久久久久大尺度免费视频| 18禁动态无遮挡网站| 中文字幕制服av| 国产精品蜜桃在线观看| 国产精品蜜桃在线观看| 欧美成人精品欧美一级黄| 亚洲不卡免费看| 亚洲三级黄色毛片| 99久久精品一区二区三区| 久久这里有精品视频免费| 亚洲久久久国产精品| 精品人妻熟女毛片av久久网站| 不卡视频在线观看欧美| 国产毛片在线视频| 18+在线观看网站| 一个人免费看片子| 亚洲国产毛片av蜜桃av| 男人添女人高潮全过程视频| 免费少妇av软件| 色5月婷婷丁香| 少妇高潮的动态图| 免费观看在线日韩| 亚洲av国产av综合av卡| av在线观看视频网站免费| 亚洲精品日韩av片在线观看| 国内揄拍国产精品人妻在线| 99热网站在线观看| 婷婷色综合大香蕉| 久久女婷五月综合色啪小说| a级片在线免费高清观看视频| 午夜免费观看性视频| 伊人久久精品亚洲午夜| 精品国产一区二区三区久久久樱花| 欧美日韩一区二区视频在线观看视频在线| 国产女主播在线喷水免费视频网站| 纵有疾风起免费观看全集完整版| 国产亚洲5aaaaa淫片| 婷婷色综合大香蕉| 久久精品久久久久久噜噜老黄| 十分钟在线观看高清视频www | 91精品国产国语对白视频| 日韩免费高清中文字幕av| 在线 av 中文字幕| 99热这里只有是精品在线观看| 中文在线观看免费www的网站| 国产一区有黄有色的免费视频| 国产亚洲av片在线观看秒播厂| 亚洲国产精品国产精品| 欧美精品高潮呻吟av久久| 一个人免费看片子| 纯流量卡能插随身wifi吗| 久久99热这里只频精品6学生| 狂野欧美激情性xxxx在线观看| 一区二区三区乱码不卡18| 精品一区二区三区视频在线| 国产精品伦人一区二区| 国产色爽女视频免费观看| 色5月婷婷丁香| 亚洲欧美精品自产自拍| 亚洲精品久久久久久婷婷小说| 18禁裸乳无遮挡动漫免费视频| 搡老乐熟女国产| 好男人视频免费观看在线| 亚洲国产色片| 亚洲不卡免费看| 国产精品不卡视频一区二区| 国产熟女欧美一区二区| 91aial.com中文字幕在线观看| 国产精品免费大片| 精品少妇内射三级| 在线免费观看不下载黄p国产| 又爽又黄a免费视频| 国产伦在线观看视频一区| 高清黄色对白视频在线免费看 | 亚洲欧美清纯卡通| 国产片特级美女逼逼视频| 日韩一区二区三区影片| 中文字幕亚洲精品专区| 国产一区有黄有色的免费视频| 午夜福利视频精品| 3wmmmm亚洲av在线观看| 九色成人免费人妻av| 激情五月婷婷亚洲| 亚洲情色 制服丝袜| 国产成人免费无遮挡视频| 观看美女的网站| 18禁在线播放成人免费| 久久精品久久久久久噜噜老黄| 国产熟女午夜一区二区三区 | 久久99热6这里只有精品| 超碰97精品在线观看| 成年av动漫网址| 亚洲国产成人一精品久久久| 乱人伦中国视频| 国产高清三级在线| 在线免费观看不下载黄p国产| 亚洲欧美日韩另类电影网站| 热re99久久精品国产66热6| 国产一级毛片在线| 内地一区二区视频在线| 国产精品一二三区在线看| 亚洲国产色片| 一区二区三区乱码不卡18| 一区二区三区乱码不卡18| 亚洲欧美中文字幕日韩二区| 国产精品一区二区在线不卡| 97在线视频观看| 成人国产麻豆网| 亚洲欧洲国产日韩| 国产一区二区在线观看av| 午夜福利网站1000一区二区三区| 国产精品99久久久久久久久| 在线精品无人区一区二区三| 国产精品一区二区在线观看99| 人人妻人人澡人人爽人人夜夜| 久久人妻熟女aⅴ| 边亲边吃奶的免费视频| 亚洲在久久综合| 人人妻人人爽人人添夜夜欢视频 | 久久久国产一区二区| 亚洲精品国产成人久久av| 色视频在线一区二区三区| 欧美日韩视频高清一区二区三区二| 免费观看性生交大片5| 美女福利国产在线| 性色avwww在线观看| 国产综合精华液| 高清不卡的av网站| av卡一久久| 国产淫片久久久久久久久| 啦啦啦啦在线视频资源| 精品亚洲乱码少妇综合久久| 亚洲第一av免费看| 91精品国产九色| 青春草视频在线免费观看| 欧美最新免费一区二区三区| 亚洲精品乱码久久久久久按摩| 日韩不卡一区二区三区视频在线| 成人特级av手机在线观看| 天美传媒精品一区二区| 简卡轻食公司| 超碰97精品在线观看| 极品教师在线视频| 国产无遮挡羞羞视频在线观看| 久久精品夜色国产| 国产无遮挡羞羞视频在线观看| 国产精品久久久久久精品古装| 国产色爽女视频免费观看| 观看免费一级毛片| 免费少妇av软件| videos熟女内射| 久久人人爽人人爽人人片va| 欧美日韩一区二区视频在线观看视频在线| 久久精品久久久久久噜噜老黄| 91精品国产国语对白视频| 欧美丝袜亚洲另类| 亚洲丝袜综合中文字幕| 亚洲欧美日韩东京热| 亚洲精品亚洲一区二区| av免费观看日本| 免费久久久久久久精品成人欧美视频 | 国产一区亚洲一区在线观看| 成人影院久久| 国产高清有码在线观看视频| 黄色毛片三级朝国网站 | av又黄又爽大尺度在线免费看| 亚洲精品乱久久久久久| 午夜激情久久久久久久| 国产黄片视频在线免费观看| 大片免费播放器 马上看| 日韩欧美精品免费久久| 色婷婷久久久亚洲欧美| 久久精品国产亚洲网站| 日韩视频在线欧美| 午夜福利,免费看| 寂寞人妻少妇视频99o| 精品久久久久久久久亚洲| 中国三级夫妇交换| 久久久久人妻精品一区果冻| 2018国产大陆天天弄谢| 精品久久国产蜜桃| 亚洲情色 制服丝袜| 成年人午夜在线观看视频| 久久精品国产亚洲av天美| 日韩精品免费视频一区二区三区 | 免费大片18禁| 美女福利国产在线| 国产成人aa在线观看| 日日啪夜夜爽| 在线播放无遮挡| 成人亚洲欧美一区二区av| 日韩av在线免费看完整版不卡| 免费观看在线日韩| 亚洲国产av新网站| 成人18禁高潮啪啪吃奶动态图 | 精品亚洲成a人片在线观看| 精品一品国产午夜福利视频| 岛国毛片在线播放| 我的老师免费观看完整版| av播播在线观看一区| 最黄视频免费看| 99久久中文字幕三级久久日本| 久久ye,这里只有精品| 黄色毛片三级朝国网站 | 色视频在线一区二区三区| 久久久久网色| 亚洲精品乱码久久久v下载方式| 久久鲁丝午夜福利片| 伦理电影大哥的女人| 国产精品成人在线| 亚洲国产精品专区欧美| 亚洲一级一片aⅴ在线观看| 视频中文字幕在线观看| 一级毛片我不卡| 婷婷色综合www| 免费黄色在线免费观看| 久久6这里有精品| 好男人视频免费观看在线| 妹子高潮喷水视频| 精品亚洲成a人片在线观看| 成年av动漫网址| 另类亚洲欧美激情| 国产永久视频网站| 99久久精品热视频| 久久国产亚洲av麻豆专区| 亚洲第一区二区三区不卡| 欧美日韩在线观看h| 国产精品人妻久久久影院| 天美传媒精品一区二区| 老司机亚洲免费影院| 久久久久久人妻| 99久久精品一区二区三区| 亚洲经典国产精华液单| 高清黄色对白视频在线免费看 | 美女脱内裤让男人舔精品视频| 男人舔奶头视频| 视频中文字幕在线观看| 韩国av在线不卡| 国产伦在线观看视频一区| 欧美日韩在线观看h| 久久青草综合色| 伊人久久国产一区二区| 欧美xxⅹ黑人| 五月玫瑰六月丁香| 日本av手机在线免费观看| 韩国av在线不卡| 国产亚洲av片在线观看秒播厂| 中国美白少妇内射xxxbb| 国产又色又爽无遮挡免| 亚洲av欧美aⅴ国产| 我要看黄色一级片免费的| 人妻系列 视频| 热99国产精品久久久久久7| 99久久精品国产国产毛片| 丝袜在线中文字幕| 日韩精品免费视频一区二区三区 | 九色成人免费人妻av| 免费av中文字幕在线| 久久久久人妻精品一区果冻| 大香蕉97超碰在线| 观看免费一级毛片| 人妻制服诱惑在线中文字幕| 国产精品无大码| 18禁在线播放成人免费| 国产淫片久久久久久久久| 精品一区二区三卡| 免费大片黄手机在线观看| 中文天堂在线官网| 日韩伦理黄色片| 最近中文字幕高清免费大全6| 人妻系列 视频| 精品一区二区三卡| 欧美最新免费一区二区三区| 如日韩欧美国产精品一区二区三区 | 国产女主播在线喷水免费视频网站| 我的女老师完整版在线观看| 日日摸夜夜添夜夜爱| 美女福利国产在线| 啦啦啦啦在线视频资源| 又黄又爽又刺激的免费视频.| 日本欧美视频一区| av有码第一页| 一级,二级,三级黄色视频| 国产男人的电影天堂91| 91精品伊人久久大香线蕉| 午夜91福利影院| 成人18禁高潮啪啪吃奶动态图 | 午夜精品国产一区二区电影| 少妇的逼水好多| 精品国产露脸久久av麻豆| 一本色道久久久久久精品综合| 中国国产av一级| 秋霞在线观看毛片| av视频免费观看在线观看| 日韩欧美 国产精品| www.色视频.com| 亚洲国产最新在线播放| 免费黄色在线免费观看| 亚洲人成网站在线播| 午夜视频国产福利| 一级毛片黄色毛片免费观看视频| 日韩不卡一区二区三区视频在线| 偷拍熟女少妇极品色| 99热全是精品| 久热这里只有精品99| 免费av不卡在线播放| 亚洲欧美成人综合另类久久久| 在线观看免费高清a一片| 欧美成人午夜免费资源| 色94色欧美一区二区| 一区二区av电影网| 亚洲婷婷狠狠爱综合网| 99久久人妻综合| 一级爰片在线观看| 三级经典国产精品| 人人妻人人澡人人看| 久久久久久久精品精品| 91精品国产九色| 少妇被粗大的猛进出69影院 | 久久韩国三级中文字幕| 51国产日韩欧美| 亚洲人与动物交配视频| 一个人免费看片子| 久久热精品热| 国产精品一区二区在线不卡| 另类精品久久| 国产男女超爽视频在线观看| 精品国产一区二区久久| 啦啦啦视频在线资源免费观看| 麻豆乱淫一区二区| 一级毛片aaaaaa免费看小| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 老司机影院成人| 午夜免费男女啪啪视频观看| 最新中文字幕久久久久| 国产 精品1| 成人免费观看视频高清| 亚洲自偷自拍三级| 中文字幕精品免费在线观看视频 | 黄色视频在线播放观看不卡| 国产伦精品一区二区三区四那| 伊人亚洲综合成人网| √禁漫天堂资源中文www| 久久久久国产精品人妻一区二区| 亚洲欧美成人精品一区二区| av在线观看视频网站免费| 日日摸夜夜添夜夜爱| 免费黄网站久久成人精品| 国产精品.久久久| 99热网站在线观看| 精品午夜福利在线看| 国产精品嫩草影院av在线观看| 亚洲丝袜综合中文字幕| 日本与韩国留学比较| 最近中文字幕高清免费大全6| 国产乱来视频区| 看免费成人av毛片| 午夜久久久在线观看| 99热这里只有是精品在线观看| 日韩欧美一区视频在线观看 | 成人国产av品久久久| 日韩精品有码人妻一区| 国内少妇人妻偷人精品xxx网站| 国产男女内射视频| 在线 av 中文字幕| 99久久精品一区二区三区| 熟女电影av网| 国产国拍精品亚洲av在线观看| 国产在视频线精品| 国产欧美日韩精品一区二区| 国产成人freesex在线| 免费av中文字幕在线| 久久久久久久国产电影| 精品酒店卫生间| 亚洲伊人久久精品综合| 亚洲综合精品二区| kizo精华| 美女国产视频在线观看| 国产黄片视频在线免费观看| 91成人精品电影| 成人毛片a级毛片在线播放| 欧美性感艳星| 免费黄网站久久成人精品| 中文字幕人妻丝袜制服| 97在线视频观看| 免费在线观看成人毛片| 久久女婷五月综合色啪小说| 日韩伦理黄色片| 久久久久久久久大av| 在线看a的网站| 9色porny在线观看| 亚洲精品国产av蜜桃| av.在线天堂| 久久人人爽人人爽人人片va| 麻豆乱淫一区二区| 看免费成人av毛片| 免费不卡的大黄色大毛片视频在线观看| 偷拍熟女少妇极品色| 久久精品国产亚洲网站| 91精品一卡2卡3卡4卡| av福利片在线| 国产片特级美女逼逼视频| 久久精品国产亚洲av涩爱| 十分钟在线观看高清视频www | 欧美精品一区二区大全| 久久人人爽人人爽人人片va| 久久精品久久久久久久性| 免费大片黄手机在线观看| 中文精品一卡2卡3卡4更新| 美女视频免费永久观看网站| av福利片在线观看| 精品亚洲成a人片在线观看| 亚洲精华国产精华液的使用体验| 欧美日韩在线观看h| 久久6这里有精品| 在线观看人妻少妇| 街头女战士在线观看网站| 国模一区二区三区四区视频| 国产又色又爽无遮挡免| 寂寞人妻少妇视频99o| 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 欧美高清成人免费视频www| 日韩不卡一区二区三区视频在线| 久久人妻熟女aⅴ| 亚洲婷婷狠狠爱综合网| 国产高清不卡午夜福利| 免费av中文字幕在线| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看 | 黑人高潮一二区| 国产黄色视频一区二区在线观看| 日韩精品免费视频一区二区三区 | 人妻系列 视频| 麻豆成人av视频| 大片电影免费在线观看免费| 大香蕉97超碰在线| 亚洲三级黄色毛片| 五月伊人婷婷丁香| 亚洲欧美清纯卡通| 亚洲无线观看免费| 18禁在线播放成人免费| 九九爱精品视频在线观看| 夫妻性生交免费视频一级片| 日本午夜av视频| 色婷婷久久久亚洲欧美| 三上悠亚av全集在线观看 | 久久久亚洲精品成人影院| 午夜91福利影院| 亚洲,欧美,日韩| 在线观看三级黄色| 男人狂女人下面高潮的视频| 肉色欧美久久久久久久蜜桃| 国精品久久久久久国模美| 国产亚洲av片在线观看秒播厂| 男的添女的下面高潮视频| 毛片一级片免费看久久久久| 大又大粗又爽又黄少妇毛片口| 国产免费一级a男人的天堂| 黄色毛片三级朝国网站 | 国产精品一区二区在线观看99| √禁漫天堂资源中文www| 久久久国产欧美日韩av| 久久久久久久久久久久大奶| 91久久精品国产一区二区成人| 美女中出高潮动态图| 亚洲av成人精品一二三区| 少妇 在线观看|