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

    氣吸滾筒式玉米排種器充種性能仿真與試驗(yàn)優(yōu)化

    2017-07-31 20:54:59衣淑娟
    關(guān)鍵詞:供種種器滾筒

    張 昆 衣淑娟

    (黑龍江八一農(nóng)墾大學(xué)工程學(xué)院,大慶163319)

    氣吸滾筒式玉米排種器充種性能仿真與試驗(yàn)優(yōu)化

    張 昆 衣淑娟

    (黑龍江八一農(nóng)墾大學(xué)工程學(xué)院,大慶163319)

    為了提高氣吸滾筒式排種器充種性能,采用離散元分析的方法,對(duì)種層高度、振動(dòng)頻率、振動(dòng)角度分別進(jìn)行數(shù)值模擬,結(jié)果表明:在相同條件下提升種層高度,可以增長(zhǎng)充種區(qū)弧長(zhǎng),增加充種時(shí)間,降低排種器的漏充率;振動(dòng)頻率增加,種子平均法向應(yīng)力方差增大,即對(duì)種子的擾動(dòng)性增強(qiáng);合適的振動(dòng)角度可以有效提高供種高度。減小內(nèi)摩擦、增強(qiáng)種群擾動(dòng)性、提高供種高度均可有效提高排種器充種性能。為尋找最佳參數(shù)組合,以鄭單958玉米種子為播種對(duì)象,采用二次旋轉(zhuǎn)正交組合試驗(yàn)方法,對(duì)排種器進(jìn)行了排種性能試驗(yàn),建立了種層高度、振動(dòng)頻率、振動(dòng)角度3個(gè)主要因素與合格率、漏播率、重播率的數(shù)學(xué)模型,分析了各個(gè)因素及交互作用對(duì)合格率的影響規(guī)律,并進(jìn)行了參數(shù)優(yōu)化與驗(yàn)證試驗(yàn)。當(dāng)最佳參數(shù)組合為振動(dòng)角度45°,振動(dòng)頻率116~122 Hz,種層高度96~117mm時(shí),合格率大于90%,漏播率小于5%,重播率小于5%。經(jīng)試驗(yàn)驗(yàn)證,試驗(yàn)結(jié)果與分析結(jié)果基本一致。試驗(yàn)結(jié)果表明該氣吸滾筒式精密排種器對(duì)于玉米種子具有很好的播種適應(yīng)性。

    玉米;氣吸滾筒式排種器;充種性能;數(shù)值模擬;正交試驗(yàn)

    引言

    玉米機(jī)械化精量播種可以節(jié)省種子,無(wú)需間苗作業(yè),而且可使種子在田間分布均勻,播深一致,達(dá)到苗齊、苗壯,是玉米播種技術(shù)的發(fā)展方向。精量排種器作為播種機(jī)的核心部件,其排種性能成為玉米播種技術(shù)的關(guān)鍵[1-3]。

    目前,國(guó)內(nèi)外先進(jìn)的播種機(jī)普遍采用氣吸式工作原理[4-9],這主要是因?yàn)闄C(jī)械式排種器整體上存在作業(yè)速度不高、對(duì)種子形狀和尺寸要求嚴(yán)格、播種單粒率有待改善等問(wèn)題,而氣吸式排種器借助氣流作用充種、清種或攜種,與機(jī)械式排種器相比,對(duì)種子形狀和大小一致性要求不嚴(yán)格,更易于實(shí)現(xiàn)高速、精密播種。氣吸式排種器,其充種環(huán)節(jié)尤為重要,基本決定了排種器的優(yōu)劣[10]。充種過(guò)程處理不當(dāng)會(huì)造成重播或者漏播現(xiàn)象,嚴(yán)重影響排種器的工作性能,導(dǎo)致播種質(zhì)量下降。目前采用種箱振動(dòng)方式提高充種是通用的方法[11-13]。王淑銘等[14]對(duì)氣動(dòng)振動(dòng)式精密排種器振盤及種子的運(yùn)動(dòng)理論進(jìn)行了分析;陳進(jìn)等[15]采用離散元法對(duì)氣吸式精密播種機(jī)振動(dòng)種盤中水稻種群運(yùn)動(dòng)進(jìn)行研究;種群的運(yùn)動(dòng)特性對(duì)于充種性能有較大的影響。李林[16]對(duì)氣吸式排種器工作過(guò)程進(jìn)行了研究,提出了種群內(nèi)摩擦力對(duì)種子充填和所需吸力有很大影響。陳進(jìn)等[17]、祁兵等[18]利用電磁激振的方法擾動(dòng)種群,使種群產(chǎn)生“沸騰”,提高了氣吸滾筒排種器充種性能。為提高氣吸滾筒式排種器的充種性能,結(jié)合上述學(xué)者的研究成果,本文采用振動(dòng)供種的方式,研究種層高度、振動(dòng)頻率和振動(dòng)角度對(duì)充種性能的影響。

    本文以自行研制的氣吸式滾筒排種器為研究對(duì)象,借助離散元仿真軟件EDEM進(jìn)行種層高度、振動(dòng)頻率、振動(dòng)角度對(duì)種群擾動(dòng)及供種高度影響的仿真分析,并通過(guò)臺(tái)架試驗(yàn)驗(yàn)證,探尋種群擾動(dòng)和供種高度對(duì)充種性能的影響。利用二次旋轉(zhuǎn)正交組合試驗(yàn)方法,對(duì)上述因素進(jìn)行排種性能試驗(yàn),以尋找最佳參數(shù)組合來(lái)提高排種性能。

    1 排種器結(jié)構(gòu)與工作原理

    氣吸滾筒式精密排種器主要由種箱、滾筒、空心軸、絕壓輥、調(diào)節(jié)螺桿及傳動(dòng)鏈輪組成,其整體結(jié)構(gòu)如圖1所示。

    圖1 氣吸滾筒排種器結(jié)構(gòu)圖Fig.1 Structure diagram of pneumatic cylinder seed-metering device

    空心軸兩端,一端與風(fēng)機(jī)連接,另一端封閉,并與滾筒內(nèi)部負(fù)壓填充區(qū)相連通??招妮S上下開孔插入可調(diào)節(jié)的螺桿結(jié)構(gòu),受到彈簧的壓力作用使得絕壓輥與滾筒內(nèi)壁緊密接觸,從而起到隔絕負(fù)壓的作用。排種器在工作狀態(tài)下,滾筒會(huì)圍繞空心軸進(jìn)行逆時(shí)針運(yùn)動(dòng),在風(fēng)機(jī)的作用下,此時(shí)的種子受到自身重力、吸孔負(fù)壓吸力的作用充入型孔當(dāng)中,并隨滾筒轉(zhuǎn)動(dòng)。當(dāng)種子轉(zhuǎn)動(dòng)到絕壓輥下方時(shí),種子失去負(fù)壓吸力,依靠自身重力落入種床帶上,完成種子的精密排種。

    1.1 滾筒直徑

    目前國(guó)內(nèi)外氣吸滾筒式排種器的滾筒直徑多為140~260 mm,大直徑滾筒可以增加吸種孔的數(shù)量從而降低滾筒轉(zhuǎn)速,增加合格率,但是相應(yīng)也增大了負(fù)壓腔的空間,需要風(fēng)機(jī)提供的空氣流量也相應(yīng)增加,能耗必然增加,也易出現(xiàn)漏氣現(xiàn)象;綜合考慮滾筒直徑選取200mm,材料選取厚度為1 mm的不銹鋼鐵皮。

    1.2 吸種孔位置尺寸

    設(shè)計(jì)滾筒上吸種孔周向孔數(shù)時(shí),既要考慮到增加吸種孔數(shù)量有利于充種,又要考慮到兩孔之間的弧長(zhǎng)Δl不小于兩粒種子的最大尺寸。

    通過(guò)測(cè)量500粒種子得出所用玉米種子平均長(zhǎng)度9.3mm、寬度6.8mm、厚度4.6mm,其中最大尺寸不超過(guò)12.2mm,綜上兩點(diǎn)考慮以及加工方便,選取滾筒周向吸種孔的數(shù)量為8個(gè)。

    吸種孔直徑k(mm)參照經(jīng)驗(yàn)公式[20]

    式中 D——滾筒直徑,mm

    通過(guò)大量測(cè)量得出種子平均直徑為8.32 mm,形狀近似球體,由式(2)中可以得到吸種孔直徑范圍,即4.99~5.82mm。為了便于實(shí)際加工,選擇吸種孔直徑為5.5mm。

    1.3 投種裝置結(jié)構(gòu)設(shè)計(jì)

    投種裝置的主要部分是絕壓輥。絕壓輥隨著滾筒的轉(zhuǎn)動(dòng)而轉(zhuǎn)動(dòng),絕壓輥與滾筒內(nèi)壁的摩擦形式由傳統(tǒng)的滑動(dòng)摩擦改為滾動(dòng)摩擦[21],大大降低了因摩擦而損失的能耗,延長(zhǎng)了機(jī)構(gòu)的使用壽命。同時(shí),為了更好地保證絕壓輥與滾筒內(nèi)壁的緊密接觸,在絕壓輥周圍套上一層壁厚為2 mm的橡膠套,以橡膠套的微量形變使其與滾筒內(nèi)壁的接觸力加大,同時(shí)為了保證摩擦輥與滾筒內(nèi)壁的過(guò)盈配合,除了在絕壓輥外部套上橡膠套之外,還在其兩側(cè)焊接上彈簧調(diào)節(jié)機(jī)構(gòu),既保證滾筒與中心軸之間的固定,又避免了長(zhǎng)時(shí)間的摩擦使之與滾筒內(nèi)壁接觸不緊密。兩端螺母則起到了調(diào)節(jié)距離的作用,如圖2所示。

    圖2 摩擦輥與彈簧調(diào)節(jié)機(jī)構(gòu)Fig.2 Friction roller and spring adjustmentmechanism

    2 仿真模型建立

    2.1 接觸模型

    離散元法描述的是顆粒之間碰撞的過(guò)程,也就是接觸的產(chǎn)生和發(fā)生作用的過(guò)程。離散元中根據(jù)接觸方式的不同可分為硬顆粒接觸和軟顆粒接觸,軟顆粒接觸方式允許顆粒之間接觸點(diǎn)間出現(xiàn)重疊部分,并根據(jù)接觸顆粒的物理屬性和法向重疊量、切向位移計(jì)算出接觸力,而顆粒間內(nèi)摩擦力即切向力與法向力有著密切的聯(lián)系,因此本文采用軟顆粒接觸模型。考慮到種子表面無(wú)粘附力,故本文選取Hertz-Mindlin無(wú)滑移接觸模型[22-25]。

    2.2 接觸力計(jì)算

    法向力Fnij是彈簧和法向阻尼器作用在顆粒i上的彈性力和阻尼力的合力,同時(shí)顆粒i發(fā)生滑移,通過(guò)推導(dǎo)發(fā)現(xiàn)切向力Ftij與法向力Fnij關(guān)系為[26-28]

    式中 μs——靜摩擦因數(shù)

    vct——接觸點(diǎn)的滑移速度,m/s

    從式中可知,種子間切向力即為內(nèi)摩擦力,而內(nèi)摩擦力是阻礙種子間相對(duì)運(yùn)動(dòng)或者有相對(duì)運(yùn)動(dòng)趨勢(shì)的主要因素[29]。從式(3)中可以看出,影響內(nèi)摩擦力的因素有靜摩擦因數(shù)和法向力,而靜摩擦因數(shù)又由種子本身物理屬性決定,因此采用法向力Fnij來(lái)衡量種子之間的內(nèi)摩擦力[30-31]。

    2.3 顆粒模型

    采用鄭單958玉米種子進(jìn)行仿真試驗(yàn),由文獻(xiàn)[32-33]得到種子的幾何尺寸與力學(xué)特征參數(shù),如表1所示。

    表1 種子特征參數(shù)Tab.1 Seeds characteristic parameters

    對(duì)鄭單958種子進(jìn)行建模,將1 000粒種子分為大扁、小扁、大圓、小圓4類,依次按其長(zhǎng)、寬、高尺寸的正態(tài)分布均值進(jìn)行三維模型建立,并在EDEM軟件中通過(guò)多球面組合的填充方式,完成一個(gè)與種子模型外輪廓較為吻合的表面組作為一個(gè)顆粒,模擬顆粒如圖3所示[34],仿真過(guò)程中隨機(jī)生成4種形狀的種子。

    圖3 種子仿真結(jié)果Fig.3 Simulation results of seed

    2.4 排種器供種過(guò)程仿真

    種子由種箱模型內(nèi)部的顆粒工廠生成,生成的種子在自身重力的作用下分布在種箱底部。在對(duì)種箱添加激振作用下,種箱內(nèi)的種群出現(xiàn)“沸騰”狀態(tài),種間出現(xiàn)了相對(duì)的運(yùn)動(dòng)與碰撞。

    為了研究氣吸滾筒式排種器充種性能,本文對(duì)種層高度、振動(dòng)頻率和振動(dòng)角度進(jìn)行考察。通過(guò)EDEM模擬不同種層高度條件下,充種區(qū)弧長(zhǎng)對(duì)充種效果的影響。通過(guò)EDEM模擬振動(dòng)頻率條件下,種間法向力的變化對(duì)供種效果的影響。通過(guò)EDEM模擬不同振動(dòng)角度條件下,供種高度對(duì)充種效果的影響。

    3 仿真結(jié)果分析

    3.1 種層高度

    根據(jù)試驗(yàn)實(shí)際情況,模擬種箱內(nèi)種層高度分別為50、100、150mm。

    通過(guò)模擬發(fā)現(xiàn),隨著種層高度的提升,排種器漏充情況明顯減少。出現(xiàn)這樣現(xiàn)象的主要原因是由于提升種層高度,充種區(qū)弧長(zhǎng)增長(zhǎng),使得滾筒在相同轉(zhuǎn)速的情況下,充種時(shí)間增長(zhǎng),一定程度上降低了排種器的漏充率。建立各因素對(duì)于充種影響的方程組

    式中 t——滾筒轉(zhuǎn)過(guò)充種區(qū)的時(shí)間,s

    lc——充種區(qū)弧長(zhǎng),m

    s——滾筒線速度,m/s

    d——滾筒直徑,m

    n——滾筒轉(zhuǎn)速,r/min

    δ——充種區(qū)弧度,rad

    由式(4)整理得出

    式(5)表明,滾筒轉(zhuǎn)過(guò)充種區(qū)的時(shí)間t僅與滾筒轉(zhuǎn)速和充種區(qū)弧度有關(guān),驗(yàn)證了上述的分析。

    通過(guò)模擬3種不同種層高度,得出排種器漏充率變化規(guī)律如圖4所示,由圖可知,隨著種層高度增大,漏充率下降。

    圖4 種層高度與漏充率關(guān)系曲線Fig.4 Relationship curve of height of seeds layer and miss rate of seed-filling

    3.2 振動(dòng)頻率

    圖5 不同振動(dòng)頻率下平均法向應(yīng)力隨時(shí)間的變化曲線Fig.5 Changing curves of average normal stresswith time under different vibration frequencies

    根據(jù)試驗(yàn)實(shí)際情況,選取振幅為1 mm,種層高度為100mm,振動(dòng)角度與水平夾角為45°時(shí),模擬種箱振動(dòng)頻率分別為65、85、105 Hz,獲取隨時(shí)間變化種子的平均法向應(yīng)力,如圖5所示。從圖中可以看出,種子平均法向應(yīng)力隨著時(shí)間的推移出現(xiàn)無(wú)規(guī)律波動(dòng),并隨著振頻的提高,波動(dòng)幅度增大且落差也隨著增大。為了定量描述這種波動(dòng),引入方差公式

    式中 D——方差,Pa2

    N——輸出數(shù)據(jù)總數(shù),N=100

    μ——種子平均法向應(yīng)力總和,Pa

    fni——i時(shí)刻每個(gè)種子受到的平均法向應(yīng)力,Pa

    方差較大的其波動(dòng)離散程度大,也就是對(duì)種子的擾動(dòng)性較大。通過(guò)計(jì)算得出振動(dòng)頻率為65 Hz時(shí),方差為39.65 Pa2。振動(dòng)頻率為85Hz時(shí),方差為58.62 Pa2。振動(dòng)頻率為105Hz時(shí),方差為61.77 Pa2。因此提高振動(dòng)頻率可以一定程度提高種群擾動(dòng)性,利于充種[33]。

    通過(guò)模擬3種不同振動(dòng)頻率,得出排種器漏充率變化規(guī)律如圖6所示,從圖中可見(jiàn),隨振動(dòng)頻率的提高,種群擾動(dòng)性增強(qiáng),漏充率下降。

    圖6 振動(dòng)頻率與漏充率關(guān)系曲線Fig.6 Relationship curve of vibration frequency and miss rate of seed-filling

    3.3 振動(dòng)角度

    根據(jù)試驗(yàn)實(shí)際情況,選取振幅為1 mm,種層高度為100 mm,振頻為105 Hz,模擬種箱振動(dòng)角度θ分別為30°、45°、60°,如圖7所示,獲取不同振動(dòng)角度情況下種箱供種高度。

    圖7 振動(dòng)角度示意圖Fig.7 Sketch of vertical oscillation angle

    通過(guò)模擬發(fā)現(xiàn),在振動(dòng)角度為30°、45°時(shí)供種高度均有提升,而在振動(dòng)角度為60°時(shí),種子出現(xiàn)了后移的現(xiàn)象,使得供種高度出現(xiàn)了下降的情況。通過(guò)測(cè)量得出振動(dòng)角度為30°、45°、60°,對(duì)應(yīng)供種高度為105、112、90mm。振動(dòng)角度的不同影響種子水平和豎直方向的位移。振動(dòng)角度為30°,其水平向前振動(dòng)位移較大,但豎直方向位移量較小,制約其進(jìn)一步提高供種高度;振動(dòng)角度為60°,其豎直方向振動(dòng)位移較大,但水平向前位移較小,導(dǎo)致種子后移,使得供種高度不但沒(méi)有提高反而下降。振動(dòng)角度為45°,既有足夠的水平向前位移又有足夠的豎直位移,使其供種高度在三者中最佳。因此比較合適的振動(dòng)角度有利于提高供種高度從而提高充種性能,但是從模擬結(jié)果來(lái)看,相互之間的差距不大。

    通過(guò)模擬3種不同振動(dòng)角度,得出排種器漏充率變化規(guī)律如圖8所示,從圖中可見(jiàn),隨著振動(dòng)角度增加,即對(duì)應(yīng)供種高度先增大后減小,漏充率先下降后上升,同樣可以說(shuō)明適當(dāng)?shù)恼駝?dòng)角度可以提高供種高度的同時(shí)降低漏充率。

    圖8 振動(dòng)角度與漏充率關(guān)系曲線Fig.8 Relationship curve of vertical oscillation angle and miss rate of seed-filling

    4 試驗(yàn)

    4.1 試驗(yàn)條件

    試驗(yàn)在JPS-12型計(jì)算機(jī)視覺(jué)排種器性能檢測(cè)試驗(yàn)臺(tái)上進(jìn)行,所用排種器為自行設(shè)計(jì)的氣吸滾筒式排種器,激振源為FP-12型氣動(dòng)活塞振動(dòng)器,振頻檢測(cè)設(shè)備采用VC63B型測(cè)振儀,如圖9所示。

    圖9 試驗(yàn)裝置圖Fig.9 Picture of test unit

    試驗(yàn)參照GB/T 6973—2005《單粒(精密)播種機(jī)試驗(yàn)方法》,連續(xù)檢測(cè)200粒種子,選取合格率、漏播率、重播率為指標(biāo)。每次試驗(yàn)重復(fù)5次取平均值。試驗(yàn)過(guò)程中設(shè)定播種機(jī)前進(jìn)速度為3 km/h,給定株距為300mm,其他條件不變。根據(jù)

    式中 vm——播種機(jī)前進(jìn)速度,km/h

    S——株距,mm

    Z——周向吸孔數(shù),Z=8

    計(jì)算出相應(yīng)播種機(jī)前進(jìn)速度為3 km/h時(shí),滾筒轉(zhuǎn)速np為20.8 r/min。

    4.2 二次旋轉(zhuǎn)正交組合試驗(yàn)

    通過(guò)前期大量的單因素試驗(yàn)確定了種層高度的取值范圍在50~150mm,振動(dòng)頻率取值范圍在65~145 Hz,振動(dòng)角度范圍在0°~90°。為了找到這3個(gè)因素的最佳參數(shù),使得此排種器排種性能最佳,本文選擇了試驗(yàn)次數(shù)少,計(jì)算方便,可以避免回歸系數(shù)間相關(guān)性的二次旋轉(zhuǎn)正交組合試驗(yàn)方法,試驗(yàn)因素編碼如表2所示。再根據(jù)三因素二次旋轉(zhuǎn)正交組合試驗(yàn)表進(jìn)行試驗(yàn),每組試驗(yàn)重復(fù)3次取平均值。試驗(yàn)方案與試驗(yàn)結(jié)果見(jiàn)表3。

    表2 試驗(yàn)因素編碼Tab.2 Factors and levels of test

    表3 試驗(yàn)設(shè)計(jì)方案與結(jié)果Tab.3 Experim ent design and results

    4.3 回歸數(shù)學(xué)模型的建立與顯著性檢驗(yàn)

    采用Design-Expert8.0.6軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行多元回歸擬合,對(duì)試驗(yàn)結(jié)果進(jìn)行回歸分析,可以得到合格率Y1、漏播率Y2和重播率Y3的回歸方程。

    (1)合格率Y1

    通過(guò)試驗(yàn)以及對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行多元回歸擬合,得到各因素對(duì)排種合格率Y1影響的回歸模型

    回歸方程的顯著性檢驗(yàn)如表4所示。

    根據(jù)表4可知,這個(gè)模型的擬合度是極顯著的(P<0.01)。但種層高度和振動(dòng)頻率交互項(xiàng)(X1X2)的P值、種層高度和振動(dòng)角度交互項(xiàng)(X1X3)的P值以及振動(dòng)頻率和振動(dòng)角度交互項(xiàng)(X2X3)的P值均大于0.1,說(shuō)明種層高度和振動(dòng)頻率的交互項(xiàng)、種層高度和振動(dòng)角度的交互項(xiàng)以及振動(dòng)頻率和振動(dòng)角度的交互項(xiàng)對(duì)排種合格率的影響不顯著,其他各項(xiàng)的F檢驗(yàn)均極顯著或顯著,說(shuō)明相關(guān)試驗(yàn)因素對(duì)響應(yīng)值的影響存在二次關(guān)系。對(duì)于失擬項(xiàng)P=0.413,不顯著,說(shuō)明不存在其他影響指標(biāo)的主要因素。剔除不顯著因素后的回歸模型

    通過(guò)對(duì)式(9)回歸系數(shù)的檢驗(yàn)得出,影響排種合格率的因素由大到小為種層高度、振動(dòng)頻率和振動(dòng)角度。

    (2)漏播率Y2

    通過(guò)試驗(yàn)以及對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行多元回歸擬合,得到各因素對(duì)漏播率Y2影響的回歸模型

    根據(jù)表4可知,這個(gè)模型的擬合度是極顯著的(P<0.01)。但種層高度和振動(dòng)角度交互項(xiàng)(X1X3)的P值大于0.1,說(shuō)明種層高度和振動(dòng)角度的交互項(xiàng)對(duì)漏播率的影響不顯著,其他各項(xiàng)的F檢驗(yàn)均極顯著或顯著,說(shuō)明相關(guān)試驗(yàn)因素對(duì)響應(yīng)值的影響存在二次關(guān)系。對(duì)于失擬項(xiàng)P=0.110 1,不顯著,說(shuō)明不存在其他影響指標(biāo)的主要因素。剔除不顯著因素后的回歸模型為

    表4 回歸方程方差分析Tab.4 Variance analysis of regression equation

    通過(guò)對(duì)式(11)回歸系數(shù)的檢驗(yàn)得出,影響漏播率的因素由大到小為種層高度、振動(dòng)頻率和振動(dòng)角度。

    (3)重播率Y3

    通過(guò)試驗(yàn)以及對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行多元回歸擬合,得到各因素對(duì)重播率Y3影響的回歸模型

    通過(guò)對(duì)式(13)回歸系數(shù)的檢驗(yàn)得出,影響重播率的因素由大到小為振動(dòng)頻率、種層高度和振動(dòng)角度。

    4.4 最佳參數(shù)組合

    設(shè)定合格率大于90%,漏播率小于5%,重播率小于5%,振動(dòng)角度為45°,最佳參數(shù)范圍如圖10所示。

    圖10 參數(shù)分析圖Fig.10 Diagram of parameters analysis

    由圖10可知,在振動(dòng)角度為45°時(shí),種層高度與振動(dòng)頻率都處于最佳值的區(qū)域圖,其中灰色區(qū)域?yàn)樽罴呀M合區(qū)域,即振動(dòng)頻率在116~122 Hz,種層高度在96~117mm時(shí),可獲得合格率大于90%,漏播率小于5%,重播率小于5%。

    對(duì)理論結(jié)果進(jìn)行試驗(yàn)驗(yàn)證。在相同的試驗(yàn)條件下選取振動(dòng)角度為45°,振動(dòng)頻率為120 Hz,種層高度為100mm進(jìn)行3次重復(fù)驗(yàn)證試驗(yàn),得到排種器合格率平均值為92.36%,且均大于90%。漏播率平均值為3.52%,且均小于5%。重播率平均值為4.21%,且均小于5%,試驗(yàn)結(jié)果與最佳組合區(qū)域結(jié)果基本相符。

    5 結(jié)論

    (1)為探尋種層高度、種群擾動(dòng)和供種高度對(duì)排種器充種性能的影響,采用離散元分析方法,利用EDEM軟件,分別對(duì)種層高度、振動(dòng)頻率以及振動(dòng)角度進(jìn)行數(shù)值模擬,得出增長(zhǎng)充種區(qū)弧長(zhǎng)、增大種群擾動(dòng)強(qiáng)度和提高供種高度可以有效提高排種器充種性能。

    (2)采用三因素五水平二次正交旋轉(zhuǎn)組合試驗(yàn)方法進(jìn)行試驗(yàn),并對(duì)試驗(yàn)結(jié)果進(jìn)行方差分析,得出影響排種合格率的因素主次順序?yàn)榉N層高度、振動(dòng)頻率和振動(dòng)角度。

    (3)利用Design-Expert8.0.6進(jìn)行數(shù)據(jù)處理,以排種合格率、漏播率、重播率為評(píng)價(jià)指標(biāo),得出在振動(dòng)角度為45°,振動(dòng)頻率在116~122 Hz,種層高度在96~117mm時(shí),合格率可大于90%,漏播率小于5%,重播率小于5%,經(jīng)試驗(yàn)驗(yàn)證,與分析結(jié)果基本一致。

    1 許劍平,謝宇峰,陳寶昌.國(guó)外氣力式精密播種機(jī)技術(shù)現(xiàn)狀及發(fā)展趨勢(shì)[J].農(nóng)機(jī)化研究,2008,30(12):203-206.XU Jianping,XIE Yufeng,CHEN Baochang.The present technique status and developing tendency of abroad pneumatic precision drill[J].Journal of Agricultural Mechanization Research,2008,30(12):203-206.(in Chinese)

    2 李國(guó)林,宋煒,毛俐,等.國(guó)內(nèi)外幾種主要排種器的特點(diǎn)[J].農(nóng)業(yè)科技與裝備,2011(8):70-71,73.LIGuolin,SONGWei,MAO Li,etal.Features of somemain seed-metering devices athome and abroad[J].Agricultural Science and Technology and Equipment,2011(8):70-71,73.(in Chinese)

    3 張澤平,馬成林,左春檉.精播排種器及排種理論研究進(jìn)展[J].吉林工業(yè)大學(xué)學(xué)報(bào),1995,25(4):112-117.ZHANG Zeping,MA Chenglin,ZUO Chuncheng.The development of the seed-metering device for precision planter and its theoretical study[J].Journal of Jilin University of Technology,1995,25(4):112-117.(in Chinese)

    4 YAZGIA,Degirmencioglu A.Optimisation of the seed spacing uniformity performance of a vacuum-type precision seeder using response surfacemethodology[J].Biosystems Engineering,2007,97(3):347-356.

    5 BEREKETB Z.Effectof differentoperating parameters on seed holding in the single seedmetering unitof a pneumatic planter[J].Turkish Journal of Agricultural Machinery,2004,28(6):435-441.

    6 GAIKWAD B B,SIROHIN PS.Design of a low-cost pneumatic seeder for nursery plug trays[J].Biosystems Engineering,2008,99(3):322-329.

    7 周海波,馬旭,姚亞利.水稻秧盤育秧播種技術(shù)與裝備的研究現(xiàn)狀及發(fā)展趨勢(shì)[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(4):301-305.ZHOU Haibo,MA Xu,YAO Yali.Research advances and prospects in the seeding technology and equipment for tray nursing seedlings of rice[J].Transactions of the Chinese Society of Agricultural Engineering,2008,24(4):301-305.(in Chinese)

    8 廖慶喜,李繼波,覃國(guó)良.氣力式油菜精量排種器試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(8):44-48.LIAO Qingxi,LIJibo,QIN Guoliang.Experiment of pneumatic precision metering device for rapeseed[J].Transactions of the Chinese Society for Agricultural Machinery,2009,40(8):44-48.(in Chinese)

    9 張曉慧,宋建農(nóng).針吸滾筒式水稻排種器設(shè)計(jì)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(3):69-71.ZHANG Xiaohui,SONG Jiannong.Design and research on rice precision needle-like vacuum seeder[J].Transactions of the Chinese Society for Agricultural Machinery,2009,40(3):69-71.(in Chinese)

    10 楊麗,史嵩,崔濤,等.氣吸與機(jī)械輔助附種結(jié)合式玉米精量排種器[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2012,43(增刊):48-53.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.a(chǎn)spx?file_no=2012s10&flag=1.DOI:10.6041/j.issn.1000-1298.2012.S0.010.YANG Li,SHI Song,CUI Tao,et al.Air-suction corn precision metering device with mechanical supporting plate to assist carrying seed[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2012,43(Supp.):48-53.(in Chinese)

    11 胡建平,鄭賽男,劉文東.磁吸滾筒式精密排種器設(shè)計(jì)與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(3):60-63.HU Jianping,ZHENG Sainan,LIUWendong.Design and experiment of precision magnetic cylinder seeder[J].Transactions of the Chinese Society for Agricultural Machinery,2009,40(3):60-63.(in Chinese)

    12 趙湛,李耀明,陳進(jìn),等.氣吸滾筒式排種器吸種過(guò)程的動(dòng)力學(xué)分析[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(7):112-116.ZHAO Zhan,LI Yaoming,CHEN Jin,et al.Dynamic analysis of seeds pick-up process for vacuum-cylinder seeder[J].Transactions of the CSAE,2011,27(7):112-116.(in Chinese)

    13 張石平,夏靜,陳進(jìn).氣吸振動(dòng)式蔬菜穴盤育苗精密播種裝置的研究[J].農(nóng)機(jī)化研究,2007,29(8):80-83.ZHANG Shiping,XIA Jing,CHEN Jin.Study on equipment of vibration air-suction vegetable tray precise seeder[J].Journal of Agricultural Mechanization Research,2007,29(8):80-83.(in Chinese)

    14 王淑銘,魏天路,周海波.氣動(dòng)振動(dòng)式精密排種器工作參數(shù)的分析與仿真[J].農(nóng)業(yè)工程學(xué)報(bào),2010,26(增刊):56-60.WANG Shuming,WEITianlu,ZHOU Haibo.Theoretical research and simulation analysis on operating parameters of pneumatic vibration precision seed-metering[J].Transactions of the CSAE,2010,26(Supp.):56-60.(in Chinese)

    15 陳進(jìn),周韓,趙湛,等.基于EDEM的振動(dòng)種盤中水稻種群運(yùn)動(dòng)規(guī)律研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2011,42(10):79-83.CHEN Jin,ZHOU Han,ZHAO Zhan,etal.Analysis of rice seedsmotion on vibrating plate using EDEM[J].Transactions of theChinese Society for Agricultural Machinery,2011,42(10):79-83.(in Chinese)

    16 李林.氣吸式排種器理論及試驗(yàn)的初步研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),1979,10(3):56-63.LILin.A preliminary study on the theory and experimentation of the suction-type metering device for precision drill[J].Transactions of the Chinese Society for Agricultural Machinery,1979,10(3):56-63.(in Chinese)

    17 陳進(jìn),李耀明.氣吸振動(dòng)式播種試驗(yàn)臺(tái)內(nèi)種子運(yùn)動(dòng)規(guī)律的研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2002,33(1):47-50.CHEN Jin,LIYaoming.Study on seedsmovement law in sowing test stand with suction and vibration[J].Transactions of the Chinese Society for Agricultural Machinery,2002,33(1):47-50.(in Chinese)

    18 祁兵,張東興,崔濤.中央集排氣送式玉米精量排種器設(shè)計(jì)與試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(18):8-15.QI Bing,ZHANG Dongxing,CUI Tao.Design and experiment of centralized pneumatic metering device for maize[J].Transactions of the CSAE,2013,29(18):8-15.(in Chinese)

    19 張順,夏俊芳,周勇,等.氣力滾筒式水稻直播精量排種器的設(shè)計(jì)與試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(1):11-19.ZHANG Shun,XIA Junfang,ZHOU Yong,et al.Design and experiment of pneumatic cylinder-type precision direct seedmetering device for rice[J].Transactions of the CSAE,2015,31(1):11-19.(in Chinese)

    20 張國(guó)忠,羅錫文,臧英,等.水稻氣力式排種器群布吸孔吸種盤吸種精度試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(6):13-20.ZHANG Guozhong,LUO Xiwen,ZANG Ying,et al.Experiment of sucking precision of sucking plate with group holes on rice pneumatic metering device[J].Transactions of the CSAE,2013,29(6):13-20.(in Chinese)

    21 高筱鈞,周金華,賴慶輝.中草藥三七氣吸滾筒式精密排種器的設(shè)計(jì)與試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(2):20-28.GAO Xiaojun,ZHOU Jinhua,LAIQinghui.Design and experiment of pneumatic cylinder precision seed-metering device for Panax notoginseng[J].Transactions of the CSAE,2016,32(2):20-28.(in Chinese)

    22 VAN LIEDEKERKE P,TIJSKENSE,DINTWA E,etal.DEM simulations of the particle flow on a centrifugal fertilizer spreader[J].Powder Technology,2009,190(3):348-360.

    23 YOSHIYUKI S,CUNDALL P A.Three dimensional DEM simulation of bulk handling by screw conveyors[J].Journal of Engineering Mechanics,2001,127(9):864-872.

    24 胡建平,郭坤,周春健,等.磁吸滾筒式排種器種箱振動(dòng)供種仿真與試驗(yàn)[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(8):61-65.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.a(chǎn)spx?file_no=20140810&flag=1.DOI:10.6041/j.issn.1000-1298.2014.08.010.HU Jianping,GUO Kun,ZHOU Chunjian,et al.Experimentof supplying seeds in box ofmagnetic precision cylinder-seeder[J/ OL].Transactions of the Chinese Society of Agricultural Machinery,2014,45(8):61-65.(in Chinese)

    25 李帥,李曄,馬履翱,等.基于EDEM對(duì)圓形導(dǎo)料溜槽的仿真分析[J].煤礦機(jī)械,2012,33(6):71-73.LIShuai,LIYe,MA Lüao,et al.Simulation analysis of cylindrical chute based on EDEM[J].Coal Mine Machinery,2012,33(6):71-73.(in Chinese)

    26 于建群,申燕芳,牛序堂,等.組合內(nèi)窩孔精密排種器清種過(guò)程的離散元法仿真分析[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(5): 105-109.YU Jianqun,SHEN Yanfang,NIU Xutang,et al.DEM simulation and analysis of the clearing process in precision metering device with combination inner-cell[J].Transactions of the CSAE,2008,24(5):105-109.(in Chinese)

    27 王國(guó)強(qiáng),郝萬(wàn)軍,王繼新.離散單元法及其在EDEM上的實(shí)踐[M].西安:西北工業(yè)大學(xué)出版社,2010.

    28 胡國(guó)明.顆粒系統(tǒng)的離散元素法分析仿真[M].武漢:武漢理工大學(xué)出版社,2010.

    29 吳愛(ài)祥,孫業(yè)志,劉湘平.散體動(dòng)力學(xué)理論及其應(yīng)用[M].北京:冶金工業(yè)出版社,2002.

    30 滕云楠,曾慶輝,姚紅良,等.振動(dòng)摩擦機(jī)理及其非線性動(dòng)力學(xué)特性[J].摩擦學(xué)學(xué)報(bào),2010,30(2):197-202.TENG Yunnan,ZENGQinghui,YAO Hongliang,etal.Mechanism of vibration friction and its non-linear dynamic characteristics[J].Tribology,2010,30(2):197-202.(in Chinese)

    31 孔維姝,胡林,杜學(xué)能,等.用探測(cè)棒研究顆粒堆中的最大靜摩擦力[J].物理學(xué)報(bào),2007,56(4):2318-2322.KONGWeishu,HU Lin,DU Xueneng,etal.Studies ofmaximum static friction force in granular pilewith probing rod[J].Acta Physica Sinica,2007,56(4):2318-2322.(in Chinese)

    32 崔濤,劉佳,楊麗,等.基于高速攝像的玉米種子滾動(dòng)摩擦特性試驗(yàn)與仿真[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(15):34-41.CUITao,LIU Jia,YANG Li,et al.Experiment and simulation of rolling friction characteristic of corn seed based on high-speed photography[J].Transactions of the CSAE,2013,29(15):34-41.(in Chinese)

    33 劉佳.機(jī)械氣力組合式精量排種器設(shè)計(jì)研究[D].北京:中國(guó)農(nóng)業(yè)大學(xué),2013.

    34 史嵩,張東興,楊麗,等.基于EDEM軟件的氣壓組合孔式排種器充種性能模擬與驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(3): 62-69.SHISong,ZHANG Dongxing,YANG Li,et al.Simulation and verification of seed-filling performance of pneumatic-combined holesmaize precision seed-metering device based on EDEM[J].Transactions of the CSAE,2015,31(3):62-69.(in Chinese)

    Simulation and Experimental Optim ization on Filling Seeds Performance of Seed Metering Device w ith Roller of Air-suction

    ZHANG Kun YIShujuan
    (College of Engineering,Heilongjiang Bayi Agricultural University,Daqing 163319,China)

    In order to improve the filling performance of seed metering device with roller of air-suction,the numerical simulation on the height of seeds layer,the vibration frequency and the angle of vibration were carried outwith method of discrete element analysis.The results showed that the arc length in the area of filling seeds and the time of filling seeds can be increased through increasing the height of seeds layer.The leakage rate can also be reduced under the same conditions;themean value of normal stress was increased with the increase of vibration frequency,whichmeant that the disturbance to the seeds can be enhanced;the proper vibration angle can increase the height of filling seeds effectively.The filling performance of the seeds metering device can be improved effectively by reducing the internal friction force,enhancing the disturbance of the populations and increasing the heightof filling seeds.Maize seeds (Zhengdan 958)were taken as research materials in order to determine the optimal combination of parameters.The seeding performance of metering device was studied through the quadratic rotationorthogonal combination experiment.The height of seeds layer,the vibration frequency and the vibration anglewere taken asmain factors.The qualified rate,leak seeding rate and the re-seeding ratewere taken as indexes.Themathematicalmodel between three factors and the indexeswas built.The effect of each factor on the qualified rate was analyzed and optimized parameters were tested.The effect of factors on the qualified seeding rate in primary sequence was the negative pressure,forward speed and angle of seeds suction.The optimal parameters were obtained.The qualified rate was above 90%,the leak seeding ratewas lower than 5%and the re-seeding ratewas also lower than 5%,when the vibration angle was 45°,the vibration frequency was 116~122 Hz and the height of seeds layer was 96~117mm.The results of the verified experimentwere consistentwith optimal results.Conclusions can be drawn that the seed metering device with roller of air-suction had good adaptability to maize seeds.References for thedesign of themaize seedsmetering device with roller of air-suction can be provided.

    maize;roller seeder of air-suction;filling seeds performance;numerical simulation; orthogonal experiment

    S223.2+5

    A

    1000-1298(2017)07-0078-09

    2017-01-20

    2017-05-02

    “十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(2014BAD06B04-03)

    張昆(1978—),男,助理研究員,博士生,主要從事墾區(qū)旱作農(nóng)業(yè)生產(chǎn)機(jī)械化研究,E-mail:zhangkun0814@163.com

    衣淑娟(1965—),女,教授,博士生導(dǎo)師,主要從事寒區(qū)大規(guī)模水稻生產(chǎn)全程機(jī)械化研究,E-mail:yishujuan_2005@126.com

    10.6041/j.issn.1000-1298.2017.07.010

    猜你喜歡
    供種種器滾筒
    油菜旋轉(zhuǎn)盤式高速集排器螺旋供種裝置設(shè)計(jì)與試驗(yàn)
    基于振動(dòng)盤供種的混編景觀種子帶編織機(jī)設(shè)計(jì)
    我們成功了
    基于EDEM的雙腔式棉花精量排種器排種性能仿真研究
    復(fù)合剪刃滾筒飛剪開發(fā)應(yīng)用
    稻麥兼用排種器振動(dòng)供種裝置性能試驗(yàn)
    精量排種器現(xiàn)狀及發(fā)展分析
    氣力托勺式馬鈴薯精量排種器設(shè)計(jì)
    基于改進(jìn)蟻群算法的滾筒截割軌跡規(guī)劃
    傾斜拋物線型孔輪式小麥供種裝置設(shè)計(jì)與試驗(yàn)
    久久午夜福利片| 极品少妇高潮喷水抽搐| 亚洲成av片中文字幕在线观看 | 亚洲欧美成人综合另类久久久| 欧美少妇被猛烈插入视频| 亚洲 欧美一区二区三区| 女性被躁到高潮视频| 蜜桃国产av成人99| 午夜福利,免费看| 最近最新中文字幕大全免费视频 | 国产精品国产三级专区第一集| 成人国语在线视频| 精品久久久精品久久久| 亚洲av综合色区一区| 亚洲欧美色中文字幕在线| 不卡视频在线观看欧美| 国产黄频视频在线观看| 亚洲欧洲日产国产| 下体分泌物呈黄色| 精品一区二区三卡| 亚洲av日韩在线播放| 国产熟女午夜一区二区三区| 国产亚洲最大av| 久久国内精品自在自线图片| 精品99又大又爽又粗少妇毛片| 国产成人午夜福利电影在线观看| 母亲3免费完整高清在线观看 | 丰满饥渴人妻一区二区三| 青春草视频在线免费观看| 亚洲经典国产精华液单| 亚洲熟女精品中文字幕| 校园人妻丝袜中文字幕| 久久久精品区二区三区| 国产成人一区二区在线| 如日韩欧美国产精品一区二区三区| 啦啦啦在线观看免费高清www| 热re99久久国产66热| 国产精品国产三级国产专区5o| 在线观看一区二区三区激情| 免费女性裸体啪啪无遮挡网站| 亚洲欧洲国产日韩| 毛片一级片免费看久久久久| 久久精品国产a三级三级三级| 91精品三级在线观看| 精品人妻偷拍中文字幕| 日韩av免费高清视频| 咕卡用的链子| 老女人水多毛片| 美女高潮到喷水免费观看| 超色免费av| 亚洲欧洲日产国产| 亚洲精品成人av观看孕妇| 日韩视频在线欧美| 亚洲av日韩在线播放| av天堂久久9| 男人爽女人下面视频在线观看| 乱人伦中国视频| 国产在线视频一区二区| 精品国产乱码久久久久久小说| 制服人妻中文乱码| 青春草国产在线视频| 一区福利在线观看| av网站免费在线观看视频| 日韩大片免费观看网站| 亚洲精品国产av成人精品| 亚洲精品国产av蜜桃| 国产毛片在线视频| 日本av免费视频播放| 亚洲精品国产色婷婷电影| 一级毛片我不卡| 精品99又大又爽又粗少妇毛片| 捣出白浆h1v1| 美女主播在线视频| 午夜福利视频精品| 伊人久久国产一区二区| 国产有黄有色有爽视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲在久久综合| 婷婷色麻豆天堂久久| 色哟哟·www| 不卡av一区二区三区| 精品亚洲成a人片在线观看| 色94色欧美一区二区| av福利片在线| 久久午夜福利片| www.自偷自拍.com| 叶爱在线成人免费视频播放| 成人国产av品久久久| 欧美亚洲日本最大视频资源| 亚洲精品一二三| 久久女婷五月综合色啪小说| 五月开心婷婷网| 亚洲精品美女久久av网站| 久久久久国产网址| 大陆偷拍与自拍| 中文天堂在线官网| 欧美日韩视频精品一区| 久久久久久伊人网av| 亚洲精品国产一区二区精华液| 伦理电影免费视频| av国产精品久久久久影院| 一级毛片我不卡| 国产片特级美女逼逼视频| 九色亚洲精品在线播放| 在线亚洲精品国产二区图片欧美| 波多野结衣av一区二区av| 18+在线观看网站| 国产一区二区激情短视频 | 国产精品亚洲av一区麻豆 | 久久 成人 亚洲| 亚洲精品成人av观看孕妇| 久久久精品免费免费高清| 边亲边吃奶的免费视频| 国产高清国产精品国产三级| 少妇 在线观看| 少妇被粗大的猛进出69影院| 亚洲精品久久久久久婷婷小说| 丰满乱子伦码专区| 国产亚洲一区二区精品| 国产成人精品福利久久| 色婷婷av一区二区三区视频| 美女国产高潮福利片在线看| 免费在线观看完整版高清| 最新的欧美精品一区二区| 亚洲精品国产色婷婷电影| 亚洲精品成人av观看孕妇| a级片在线免费高清观看视频| 这个男人来自地球电影免费观看 | 麻豆精品久久久久久蜜桃| 18+在线观看网站| 大香蕉久久成人网| 卡戴珊不雅视频在线播放| 亚洲精品国产av蜜桃| 国产白丝娇喘喷水9色精品| 精品国产乱码久久久久久小说| 国产免费视频播放在线视频| 一级爰片在线观看| 麻豆av在线久日| 伦理电影免费视频| 狠狠精品人妻久久久久久综合| 岛国毛片在线播放| 亚洲国产av影院在线观看| 午夜91福利影院| 国产在线一区二区三区精| 90打野战视频偷拍视频| 国产成人精品婷婷| 亚洲三区欧美一区| 美女高潮到喷水免费观看| 精品久久久久久电影网| a 毛片基地| 赤兔流量卡办理| 日产精品乱码卡一卡2卡三| 在线观看国产h片| 久久久精品区二区三区| 日韩熟女老妇一区二区性免费视频| 精品一区二区免费观看| 亚洲一区中文字幕在线| 80岁老熟妇乱子伦牲交| 婷婷色综合大香蕉| 日韩免费高清中文字幕av| 日本欧美视频一区| 伊人亚洲综合成人网| 搡老乐熟女国产| 美女视频免费永久观看网站| 熟女电影av网| 免费av中文字幕在线| 黄网站色视频无遮挡免费观看| 欧美黄色片欧美黄色片| 母亲3免费完整高清在线观看 | 亚洲,一卡二卡三卡| 黄色毛片三级朝国网站| 国产精品熟女久久久久浪| 日韩一区二区三区影片| 免费黄频网站在线观看国产| 国产精品不卡视频一区二区| 午夜福利在线观看免费完整高清在| 免费观看a级毛片全部| 欧美日韩亚洲国产一区二区在线观看 | 久久精品久久精品一区二区三区| 9热在线视频观看99| 高清在线视频一区二区三区| 久久鲁丝午夜福利片| 男女下面插进去视频免费观看| 99久久中文字幕三级久久日本| 大码成人一级视频| 成人影院久久| 色婷婷久久久亚洲欧美| 高清在线视频一区二区三区| 国产成人一区二区在线| 一级片免费观看大全| 亚洲一级一片aⅴ在线观看| 91在线精品国自产拍蜜月| 黄色视频在线播放观看不卡| 涩涩av久久男人的天堂| 欧美人与性动交α欧美精品济南到 | 亚洲精品国产一区二区精华液| 性色avwww在线观看| 亚洲国产欧美日韩在线播放| 最近最新中文字幕大全免费视频 | 高清av免费在线| 建设人人有责人人尽责人人享有的| 老汉色av国产亚洲站长工具| 精品久久久久久电影网| 九色亚洲精品在线播放| 国产成人精品在线电影| 国产精品一区二区在线不卡| 极品少妇高潮喷水抽搐| 亚洲一区二区三区欧美精品| 亚洲国产最新在线播放| 亚洲第一青青草原| 色吧在线观看| 妹子高潮喷水视频| 永久网站在线| 日本爱情动作片www.在线观看| 亚洲伊人色综图| 亚洲av欧美aⅴ国产| 国产成人av激情在线播放| 亚洲伊人色综图| 伊人亚洲综合成人网| 午夜免费鲁丝| 精品亚洲乱码少妇综合久久| 最新的欧美精品一区二区| 国产1区2区3区精品| av在线老鸭窝| 波多野结衣一区麻豆| 天天躁日日躁夜夜躁夜夜| 久久狼人影院| 欧美日韩一区二区视频在线观看视频在线| 性少妇av在线| 亚洲精品国产av成人精品| 另类亚洲欧美激情| 亚洲综合色惰| 黄网站色视频无遮挡免费观看| 国产免费现黄频在线看| 亚洲国产欧美在线一区| 蜜桃国产av成人99| 伊人亚洲综合成人网| av在线播放精品| 国产精品麻豆人妻色哟哟久久| 国产精品国产三级专区第一集| 美女福利国产在线| 成人午夜精彩视频在线观看| 精品一区二区三区四区五区乱码 | 国产视频首页在线观看| 一个人免费看片子| 哪个播放器可以免费观看大片| 国语对白做爰xxxⅹ性视频网站| 人妻 亚洲 视频| 美女主播在线视频| 亚洲精品日本国产第一区| 亚洲男人天堂网一区| 最近最新中文字幕免费大全7| 亚洲精品第二区| 女性生殖器流出的白浆| 亚洲内射少妇av| 妹子高潮喷水视频| 在线看a的网站| 久久精品久久精品一区二区三区| 在线观看免费视频网站a站| 国产 一区精品| 国产在视频线精品| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 三级国产精品片| 国产精品av久久久久免费| 一级毛片电影观看| 丁香六月天网| 99香蕉大伊视频| 五月天丁香电影| 国产成人一区二区在线| 国产精品免费大片| 女人高潮潮喷娇喘18禁视频| 亚洲精品自拍成人| 99久久精品国产国产毛片| 欧美黄色片欧美黄色片| 免费久久久久久久精品成人欧美视频| 日韩制服丝袜自拍偷拍| 国产极品天堂在线| 成年人午夜在线观看视频| 免费大片黄手机在线观看| 亚洲综合色网址| 中文字幕制服av| 电影成人av| 91精品国产国语对白视频| 久久久精品免费免费高清| 日本av免费视频播放| 国产精品嫩草影院av在线观看| 久久久久国产网址| 日韩三级伦理在线观看| 国产精品久久久久久久久免| 欧美精品一区二区大全| 国产在视频线精品| 看十八女毛片水多多多| 国产成人av激情在线播放| 欧美在线黄色| 色婷婷av一区二区三区视频| 欧美成人午夜精品| 黄频高清免费视频| 久久韩国三级中文字幕| 成人午夜精彩视频在线观看| 色吧在线观看| 18禁动态无遮挡网站| a级毛片在线看网站| 毛片一级片免费看久久久久| 1024香蕉在线观看| 如日韩欧美国产精品一区二区三区| 91精品伊人久久大香线蕉| 热99久久久久精品小说推荐| 啦啦啦在线观看免费高清www| 狠狠精品人妻久久久久久综合| 亚洲激情五月婷婷啪啪| 中文字幕精品免费在线观看视频| 不卡视频在线观看欧美| 日本猛色少妇xxxxx猛交久久| 老熟女久久久| 哪个播放器可以免费观看大片| 久久久久精品性色| 欧美日韩一区二区视频在线观看视频在线| 午夜福利视频精品| 99re6热这里在线精品视频| 国产男人的电影天堂91| 精品亚洲成国产av| 国产国语露脸激情在线看| 色婷婷久久久亚洲欧美| 色94色欧美一区二区| 国产一区二区 视频在线| 精品人妻一区二区三区麻豆| 中文欧美无线码| 午夜老司机福利剧场| 国产精品不卡视频一区二区| 哪个播放器可以免费观看大片| 午夜老司机福利剧场| 欧美在线黄色| 男女啪啪激烈高潮av片| √禁漫天堂资源中文www| 欧美成人午夜精品| 捣出白浆h1v1| 极品人妻少妇av视频| 老司机影院成人| 国产av一区二区精品久久| 黄片无遮挡物在线观看| 我的亚洲天堂| 国产乱人偷精品视频| 人成视频在线观看免费观看| 久久ye,这里只有精品| 如何舔出高潮| 夫妻性生交免费视频一级片| 可以免费在线观看a视频的电影网站 | 一级片'在线观看视频| 免费看av在线观看网站| 1024视频免费在线观看| 色婷婷久久久亚洲欧美| 精品国产露脸久久av麻豆| 天美传媒精品一区二区| 亚洲男人天堂网一区| 国产亚洲午夜精品一区二区久久| 青春草视频在线免费观看| 桃花免费在线播放| 亚洲精品一区蜜桃| 国产国语露脸激情在线看| 热re99久久精品国产66热6| 母亲3免费完整高清在线观看 | 国产欧美亚洲国产| 国产精品女同一区二区软件| 另类亚洲欧美激情| 免费av中文字幕在线| 久久国产精品男人的天堂亚洲| 一区二区av电影网| 婷婷色av中文字幕| 大香蕉久久成人网| 一区二区三区激情视频| 久久99精品国语久久久| 韩国av在线不卡| 母亲3免费完整高清在线观看 | 国产精品一区二区在线观看99| 国产日韩一区二区三区精品不卡| 少妇被粗大的猛进出69影院| 欧美另类一区| 日韩欧美一区视频在线观看| 国产男女内射视频| 久久婷婷青草| 男女边吃奶边做爰视频| 日韩欧美一区视频在线观看| √禁漫天堂资源中文www| 国产免费又黄又爽又色| 久久久久国产一级毛片高清牌| 极品少妇高潮喷水抽搐| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品古装| 人成视频在线观看免费观看| 久久久久精品性色| 母亲3免费完整高清在线观看 | 夫妻性生交免费视频一级片| 久久精品国产综合久久久| 国产一区二区激情短视频 | 80岁老熟妇乱子伦牲交| 国产有黄有色有爽视频| 免费不卡的大黄色大毛片视频在线观看| 国产精品熟女久久久久浪| 亚洲第一青青草原| 日韩制服骚丝袜av| 在线免费观看不下载黄p国产| 亚洲久久久国产精品| 日韩中文字幕视频在线看片| 日本欧美国产在线视频| 深夜精品福利| 国产一区二区在线观看av| 国产无遮挡羞羞视频在线观看| 美女脱内裤让男人舔精品视频| 黄片播放在线免费| 欧美 日韩 精品 国产| 不卡视频在线观看欧美| 日韩中文字幕视频在线看片| 婷婷色av中文字幕| www.av在线官网国产| 五月伊人婷婷丁香| av在线播放精品| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 久久综合国产亚洲精品| 日本91视频免费播放| 国产精品一国产av| 亚洲综合精品二区| 十八禁网站网址无遮挡| 天天操日日干夜夜撸| 80岁老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 卡戴珊不雅视频在线播放| 日韩制服丝袜自拍偷拍| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 五月天丁香电影| 国产探花极品一区二区| 国产精品秋霞免费鲁丝片| 一边摸一边做爽爽视频免费| 99国产综合亚洲精品| 久久人人爽av亚洲精品天堂| 久久久久久久亚洲中文字幕| 一区在线观看完整版| 久久久久精品性色| 久久女婷五月综合色啪小说| 成人漫画全彩无遮挡| 国产深夜福利视频在线观看| 国产极品天堂在线| 婷婷色综合www| 亚洲第一青青草原| 啦啦啦在线观看免费高清www| 久久久久久人人人人人| 熟女电影av网| 日韩成人av中文字幕在线观看| 99re6热这里在线精品视频| xxxhd国产人妻xxx| 亚洲中文av在线| 在线观看免费高清a一片| 欧美最新免费一区二区三区| 久久精品国产a三级三级三级| 亚洲精品视频女| 久久人人爽人人片av| 亚洲av欧美aⅴ国产| 久久久久网色| freevideosex欧美| 少妇熟女欧美另类| 亚洲熟女精品中文字幕| 九九爱精品视频在线观看| 少妇人妻久久综合中文| 国产成人精品一,二区| 午夜激情久久久久久久| 80岁老熟妇乱子伦牲交| 黄片播放在线免费| 亚洲中文av在线| 久久精品国产综合久久久| 人体艺术视频欧美日本| 2021少妇久久久久久久久久久| 国产女主播在线喷水免费视频网站| 老司机影院成人| 久久久久久久久久久久大奶| 一区二区三区激情视频| 欧美bdsm另类| videos熟女内射| 亚洲av福利一区| 亚洲精品久久成人aⅴ小说| 久久久精品94久久精品| 亚洲在久久综合| 午夜福利视频精品| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 精品人妻偷拍中文字幕| 精品国产超薄肉色丝袜足j| 超色免费av| 晚上一个人看的免费电影| 久久这里只有精品19| 精品亚洲成国产av| 成人国产av品久久久| 丝袜喷水一区| 有码 亚洲区| 国产精品不卡视频一区二区| 久久97久久精品| 老司机影院毛片| 99久久精品国产国产毛片| av.在线天堂| a级毛片黄视频| 婷婷色麻豆天堂久久| 欧美亚洲日本最大视频资源| 亚洲欧洲国产日韩| 女人高潮潮喷娇喘18禁视频| 成年美女黄网站色视频大全免费| 亚洲欧美中文字幕日韩二区| 丝袜人妻中文字幕| 婷婷色av中文字幕| 亚洲综合色惰| 日韩一区二区视频免费看| 成人18禁高潮啪啪吃奶动态图| 久久鲁丝午夜福利片| 日韩欧美精品免费久久| 精品福利永久在线观看| 超碰97精品在线观看| 久久久久精品久久久久真实原创| tube8黄色片| 激情五月婷婷亚洲| 久久精品国产亚洲av涩爱| 在线观看免费高清a一片| 免费观看av网站的网址| 97精品久久久久久久久久精品| 如日韩欧美国产精品一区二区三区| 日韩制服丝袜自拍偷拍| 一级片'在线观看视频| 亚洲,欧美精品.| 国产成人a∨麻豆精品| 欧美激情高清一区二区三区 | 人妻 亚洲 视频| 天美传媒精品一区二区| 在线天堂最新版资源| 九色亚洲精品在线播放| 日本午夜av视频| 国产亚洲av片在线观看秒播厂| 高清av免费在线| 伦精品一区二区三区| 午夜日韩欧美国产| 国产成人精品久久二区二区91 | 男人舔女人的私密视频| 日韩三级伦理在线观看| 考比视频在线观看| 制服人妻中文乱码| 91精品三级在线观看| 少妇猛男粗大的猛烈进出视频| 欧美日韩一区二区视频在线观看视频在线| 免费黄色在线免费观看| 街头女战士在线观看网站| 久久精品国产a三级三级三级| 国产男女内射视频| 免费看不卡的av| 美女福利国产在线| 欧美精品一区二区大全| 飞空精品影院首页| 久久久久视频综合| 少妇的逼水好多| 久久国产亚洲av麻豆专区| 久久99一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲av.av天堂| 另类亚洲欧美激情| 丝袜美足系列| 天天躁夜夜躁狠狠久久av| 侵犯人妻中文字幕一二三四区| 亚洲色图 男人天堂 中文字幕| 97在线人人人人妻| 国产精品欧美亚洲77777| 在线精品无人区一区二区三| 丰满少妇做爰视频| 边亲边吃奶的免费视频| 亚洲三区欧美一区| 最新中文字幕久久久久| 成年动漫av网址| 亚洲,一卡二卡三卡| 满18在线观看网站| 久热久热在线精品观看| 18禁裸乳无遮挡动漫免费视频| 尾随美女入室| 大片免费播放器 马上看| 久久久久久久久久人人人人人人| 水蜜桃什么品种好| 免费高清在线观看视频在线观看| 欧美亚洲日本最大视频资源| 欧美日韩精品成人综合77777| 日日啪夜夜爽| 国产男女内射视频| 国产精品蜜桃在线观看| 少妇的丰满在线观看| 深夜精品福利| av电影中文网址| 大香蕉久久成人网| 一个人免费看片子| 国产一区二区 视频在线| 国产成人精品在线电影| videossex国产| 老汉色av国产亚洲站长工具| 国产成人av激情在线播放| 亚洲精品一二三| 青春草亚洲视频在线观看| 日韩在线高清观看一区二区三区| 欧美精品av麻豆av| 宅男免费午夜| 观看美女的网站| 伦理电影大哥的女人| 好男人视频免费观看在线| 叶爱在线成人免费视频播放| 不卡av一区二区三区| 在线看a的网站| 丝袜脚勾引网站| 精品人妻偷拍中文字幕| 国产乱来视频区| 亚洲美女黄色视频免费看| 国产精品蜜桃在线观看| 亚洲欧美成人综合另类久久久| 97在线人人人人妻|