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

    氣送式高速玉米精量排種器導(dǎo)流渦輪設(shè)計(jì)與試驗(yàn)

    2019-12-06 03:03:50高筱鈞賀小偉張東興
    關(guān)鍵詞:型孔種器外圈

    高筱鈞 徐 楊 賀小偉 張東興 楊 麗 崔 濤

    (1.中國農(nóng)業(yè)大學(xué)工學(xué)院, 北京 100083; 2.農(nóng)業(yè)農(nóng)村部土壤-機(jī)器-植物系統(tǒng)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 北京 100083;3.塔里木大學(xué)機(jī)械電氣化工程學(xué)院, 阿拉爾 843300)

    0 引言

    高速精量播種技術(shù)是推動(dòng)玉米產(chǎn)業(yè)可持續(xù)發(fā)展的重要途徑之一[1-2]。排種器在高速排種過程中,由于排種機(jī)構(gòu)高速回轉(zhuǎn),充種時(shí)間急劇下降,加上玉米種子形狀大小差異較大,導(dǎo)致漏播指數(shù)顯著上升。為解決上述問題,國內(nèi)外學(xué)者在現(xiàn)有重力充種的原理基礎(chǔ)上開展了提高充種性能的研究[3-11],這些研究均采用輔助手段改善充種性能,未從根本上解決因高速作業(yè)、種群堆積而造成的重力充種難問題,為此,崔濤等[12]提出一種氣力送種、氣流擾種、離心清種的高速排種方法,并設(shè)計(jì)了一種氣送式高速玉米精量排種器。

    由于排種器內(nèi)部流場復(fù)雜且雜亂無章,探討不同結(jié)構(gòu)導(dǎo)流渦輪對(duì)排種器內(nèi)部流場分布以及工作性能的影響很有必要。通過流場分析以及對(duì)比試驗(yàn)發(fā)現(xiàn),導(dǎo)流渦輪可以將排種器中心空氣引導(dǎo)至外圈,增加外圈流速流量,產(chǎn)生壓覆作用力,有效提高了排種器工作性能。近年來,CFD 數(shù)值計(jì)算在農(nóng)業(yè)機(jī)械領(lǐng)域的應(yīng)用較為廣泛。通過 CFD 技術(shù),可利用計(jì)算機(jī)分析并顯示氣流場中的現(xiàn)象,使分析結(jié)果可視化、直觀化,并能在較短時(shí)間內(nèi)預(yù)測(cè)流場[13-18]。本文通過對(duì)比不同結(jié)構(gòu)的導(dǎo)流渦輪,并采用Fluent 軟件數(shù)值模擬排種器內(nèi)流場變化,確定導(dǎo)流渦輪結(jié)構(gòu),分析對(duì)排種過程各階段的影響,通過對(duì)比試驗(yàn)驗(yàn)證安裝導(dǎo)流渦輪的必要性。為獲得在高速作業(yè)下排種器最佳工作性能,以工作速度、種子喂入量和氣送風(fēng)壓為試驗(yàn)因素,開展二次回歸正交旋轉(zhuǎn)組合試驗(yàn),以期得到排種器最佳參數(shù)組合。

    1 工作原理

    1.1 排種器組成與工作原理

    氣送式高速玉米精量排種器主要由后殼體、軸承、排種盤、型孔插件、軸、導(dǎo)流渦輪、護(hù)種板和前殼體組成,其整體結(jié)構(gòu)如圖1a所示。

    圖1 氣送式高速玉米精量排種器結(jié)構(gòu)示意圖Fig.1 Structure schematic of high speed maize precision seed metering device with seed feeding in air-assisted1.后殼體 2.軸承 3.排種盤 4.型孔插件 5.軸 6.導(dǎo)流渦輪 7.護(hù)種板 8.前殼體 9.分隔板 Ⅰ.充種區(qū) Ⅱ.攜種區(qū) Ⅲ.清種區(qū) Ⅳ.投種區(qū) Ⅴ.過渡區(qū)

    氣流從進(jìn)風(fēng)口吹入,將從喂入口進(jìn)入的種子吹入文丘里管,形成均勻有序的種子流進(jìn)入排種器[19]。排種器工作過程主要分為充種、攜種、清種、投種和過渡5個(gè)階段,逆時(shí)針轉(zhuǎn)動(dòng),如圖1b所示。軸轉(zhuǎn)動(dòng)帶動(dòng)排種盤轉(zhuǎn)動(dòng),排種盤轉(zhuǎn)動(dòng)帶動(dòng)與之連接的4個(gè)型孔插件及導(dǎo)流渦輪轉(zhuǎn)動(dòng),種子從前殼體進(jìn)種口隨氣流進(jìn)入排種器。在型孔插件以及護(hù)種板的夾持作用下,充入型孔插件的種子隨之轉(zhuǎn)動(dòng),途經(jīng)充種區(qū)和攜種區(qū)后轉(zhuǎn)動(dòng)到排種器正上方處,由于護(hù)種板在充種區(qū)和攜種區(qū)的前端傾斜部與型孔插件貼合起到囊種作用,在清種區(qū)由于其前端傾斜部分變窄,此時(shí)型孔插件內(nèi)僅能容納單粒種子,多余的種子失去了護(hù)種板的支持力,做離心運(yùn)動(dòng)脫離型孔插件,甩到前殼體內(nèi)壁面并順著分隔板滑落回排種器底部,在氣流作用下與剛進(jìn)入排種器的種子匯合,進(jìn)行再次的充種。型孔插件內(nèi)留下的單粒種子繼續(xù)隨著型孔插件轉(zhuǎn)動(dòng),直到投種區(qū)失去了護(hù)種板的包裹,種子缺少了支持力,在自身重力和氣流作用下脫離型孔插件沿著落種口排出排種器,完成排種過程(圖2)。

    圖2 排種器工作原理圖Fig.2 Schematic of metering device1.進(jìn)風(fēng)口 2.喂入口 3.文丘里管 4.玉米種子 5.分隔板 6.型孔插件 7.排種器 8.護(hù)種板 9.導(dǎo)流渦輪 10.軸

    1.2 導(dǎo)流渦輪工作原理

    導(dǎo)流渦輪固定安裝在排種盤上,隨排種盤一起轉(zhuǎn)動(dòng)。工作情況下,導(dǎo)流渦輪逆時(shí)針轉(zhuǎn)動(dòng),在相同轉(zhuǎn)速情況下排種盤內(nèi)外圈角速度一致,但外圈線速度最大,公式為

    vl=ωR

    (1)

    式中vl——線速度,m/s

    ω——角速度,rad/s

    R——半徑,m

    因此導(dǎo)流渦輪外圈將會(huì)產(chǎn)生較大的線速度,帶動(dòng)外圈氣流流動(dòng),由簡化的伯努利方程可以定性地知道,外圈氣流流速快、壓強(qiáng)小,將會(huì)產(chǎn)生較大壓差,起到很好的壓覆種子作用,公式為

    (2)

    式中p——壓強(qiáng),Pa

    ρ——流體密度,kg/m3

    g——重力加速度,取9.81 m/s2

    h——高度,m

    v——流速,m/sC——常數(shù)

    導(dǎo)流渦輪可以提高排種器內(nèi)部流場的流動(dòng)性,增大排種器外圈壓差,從而增加對(duì)種子壓覆作用,使得種子流動(dòng)均勻有序,提高充種、攜種、清種和投種性能。

    2 導(dǎo)流渦輪結(jié)構(gòu)設(shè)計(jì)

    為了探討不同導(dǎo)流渦輪結(jié)構(gòu)對(duì)排種器內(nèi)部流場以及工作性能的影響,設(shè)計(jì)了3種具有代表性的導(dǎo)流渦輪,如圖3所示,3種裝置均為中心對(duì)稱結(jié)構(gòu),采用未來800樹脂材料3D打印加工制成。

    圖3 3種結(jié)構(gòu)導(dǎo)流渦輪示意圖Fig.3 Schematic of diversion turbine structure

    對(duì)3種結(jié)構(gòu)的迎風(fēng)角α分析,由圖4可知,A型結(jié)構(gòu)迎風(fēng)角為銳角,因此具有較好的導(dǎo)流作用,但其擾流作用不明顯。B型結(jié)構(gòu)迎風(fēng)角為直角,具有較好的擾流作用,但導(dǎo)流作用不顯著。C型結(jié)構(gòu)結(jié)合上述2種結(jié)構(gòu)的特點(diǎn),將攪拌葉片設(shè)計(jì)為曲線形式,增加導(dǎo)流作用,同時(shí)2種葉片迎風(fēng)角分別為銳角和直角,兼具導(dǎo)流性和擾流性。為了得出中心空氣導(dǎo)流至外圈的運(yùn)動(dòng)軌跡,假設(shè)空氣質(zhì)點(diǎn)最佳絕對(duì)運(yùn)動(dòng)軌跡方程為|y|=vt,其中t為運(yùn)動(dòng)時(shí)間。設(shè)t=0時(shí)刻空氣質(zhì)點(diǎn)初始位置坐標(biāo)為(R′,0),排種盤逆時(shí)針轉(zhuǎn)動(dòng),當(dāng)運(yùn)動(dòng)一定時(shí)間,空氣質(zhì)點(diǎn)位置坐標(biāo)為

    圖4 迎風(fēng)角分析Fig.4 Upwindangle analysis

    (3)

    式中x——漸開線橫坐標(biāo)值,mm

    y——漸開線縱坐標(biāo)值,mm

    將空氣質(zhì)點(diǎn)絕對(duì)運(yùn)動(dòng)軌跡方程|y|=vt代入式(3)中,得

    (4)

    v=R′ω,φ=ωt,將其代入式(4)中,得出空氣相對(duì)運(yùn)動(dòng)軌跡方程

    (5)

    式中φ——轉(zhuǎn)動(dòng)角度,rad

    可見,空氣相對(duì)運(yùn)動(dòng)軌跡方程即導(dǎo)流渦輪葉片曲線方程是基圓半徑為R′的漸開線,根據(jù)排種器整體結(jié)構(gòu)尺寸要求,本設(shè)計(jì)取R′為38 mm。為了增大型孔處流速流量以及與型孔相對(duì)位置,銳角迎風(fēng)角α取40°,導(dǎo)流渦輪外圈與型孔插件上邊緣相切,直角迎風(fēng)角葉片頂端與型孔插件頂端相連接,相鄰2種葉片包圍一個(gè)型孔,由于種盤上有4個(gè)型孔,因此一共需要8個(gè)葉片,內(nèi)圓半徑r為75 mm,導(dǎo)流渦輪半徑R為150 mm,這種結(jié)構(gòu)增大了型孔處空氣流速流量,起到了增大壓覆力的作用。

    2.1 充種過程

    由于排種器工作原理的改變,將原本堆積于排種器腔體內(nèi)的種子群依靠自身重力進(jìn)行充種的方式改變?yōu)殡S氣流均勻有序進(jìn)入排種器內(nèi)的種子流,避免了種子堆積帶來的種間摩擦力大,出現(xiàn)結(jié)拱和架空等不利于充種現(xiàn)象。在充種過程中,種子之間相互碰撞,不止一顆種子被囊入型孔插件中,由于導(dǎo)流渦輪的加入,種子除了受到來自種群間的擠壓力以及與護(hù)種板、型孔插件之間的接觸力,還受到流場的壓覆作用力,其中在型孔底端的種子受力如圖5所示。圖中F0為種群擠壓力,N;F1為種子受到護(hù)種板的支持力,N;F2為種子受到型孔底部支持力,N;F3為流場壓覆力所產(chǎn)生的擠壓力,N;f為種子與幾何體之間的靜摩擦力,N。

    圖5 型孔底部種子受力圖Fig.5 Force analysis of seed at bottom

    高速工作狀態(tài)下,種子所需向心力較大,同時(shí)護(hù)種板給予型孔底部種子的支持力以及種群對(duì)型孔底部種子的擠壓力均較大,再加之導(dǎo)流渦輪高速旋轉(zhuǎn)在外圈形成較高的空氣流速,產(chǎn)生較大壓覆力,使得型孔底部的種子被牢牢壓覆于型孔內(nèi),因此在高工作速度下不會(huì)出現(xiàn)漏充的現(xiàn)象。導(dǎo)流渦輪使得排種器在原工作條件的基礎(chǔ)上增加了對(duì)種子的壓覆作用,因此避免了傳統(tǒng)排種器因種群堆積,排種盤在高速轉(zhuǎn)動(dòng)下漏充嚴(yán)重的問題。

    2.2 清種過程

    隨著高速旋轉(zhuǎn)的排種盤,型孔內(nèi)的種子到達(dá)排種盤頂端清種區(qū),由于排種盤高速旋轉(zhuǎn),種子所需向心力較大,且多來自護(hù)種板的支持力,當(dāng)種子到達(dá)清種區(qū)由于護(hù)種板寬度驟減,型孔上部多余的種子瞬間失去了來自護(hù)種板的支持力,加之導(dǎo)流渦輪作用產(chǎn)生的壓覆力,壓覆力作用效果為推出多余種子做離心運(yùn)動(dòng),更好地起到了清種效果,使得多余的種子做離心運(yùn)動(dòng),脫離型孔,其受力情況如圖6所示。

    圖6 清出種子受力圖Fig.6 Force analysis of redundant seed

    在水平方向上,多余的種子受到型孔插件的支持力以及種子之間的擠壓力,處于平衡狀態(tài),在向心方向上,由于缺少護(hù)種板的支持再加上導(dǎo)流渦輪所產(chǎn)出的壓覆力,使得種子做離心運(yùn)動(dòng),受力為

    (6)

    式中Fy——排種盤向心方向的合力,N

    Fx——排種盤切線方向的合力,N

    Fn——種子受到的支持力,N

    Fr——種子所需向心力,N

    清除了多余種子,在型孔底部僅且只有一顆種子被牢牢壓覆在型孔內(nèi),由于清種區(qū)護(hù)種板寬度蓋過型孔底部種子的重心位置,因此可以提供其足夠的向心力,不會(huì)脫離型孔,實(shí)現(xiàn)清除多余種子保證單粒率。

    2.3 投種過程

    在投種區(qū),型孔內(nèi)種子失去護(hù)種板的支持力,以及受到導(dǎo)流渦輪所產(chǎn)生的壓覆力,種子做離心運(yùn)動(dòng)。在離心運(yùn)動(dòng)過程中種子第1階段先脫離型孔,其受力運(yùn)動(dòng)情況如圖7所示。

    圖7 種子脫離受力分析圖Fig.7 Force and motion analysis of seed departure

    第2階段中,種子脫離型孔的束縛,在慣性的作用下繼續(xù)保持沿著排種盤切線方向運(yùn)動(dòng)。因?yàn)閷?dǎo)流渦輪在排種器內(nèi)外圈產(chǎn)生較大空氣流速,且在投種區(qū)種子與氣流均從落種口排出,所以種子除了受到自身重力、流場壓覆力,還受到流場中曳力的作用,曳力表達(dá)式為

    (7)

    式中FD——曳力,N

    CD——繞流阻力系數(shù)

    Q——物體在流動(dòng)方向上的投影面積,m2

    導(dǎo)流渦輪的加入使得種子在投種階段減小了水平方向的速度,并且在豎直方向做加速下落運(yùn)動(dòng),增加了投種的流暢性,降低漏播風(fēng)險(xiǎn),實(shí)現(xiàn)單粒播種,如圖8所示。

    圖8 種子下落受力圖Fig.8 Force analysis of seed dropping

    在整個(gè)投種過程中,種子除了受到自身重力以及與幾何體之間的摩擦力之外,主要受到導(dǎo)流渦輪所產(chǎn)生的壓覆力和排種器內(nèi)投種區(qū)流場中曳力的作用,壓覆力作用效果為推出型孔內(nèi)種子,做離心運(yùn)動(dòng)。其中壓覆力遠(yuǎn)遠(yuǎn)大于摩擦力,起到了較好的推種作用。種子自身的重力以及流場中的曳力有效起到快速投種的作用。

    3 流場仿真與分析

    為了探討不同結(jié)構(gòu)導(dǎo)流渦輪對(duì)排種器內(nèi)流場影響情況,采用計(jì)算流體力學(xué)(CFD)的方法進(jìn)行模擬仿真,通過流場變化分析得出較好導(dǎo)流渦輪結(jié)構(gòu)。

    3.1 流場數(shù)學(xué)模型

    氣體密度隨著壓強(qiáng)或溫度的變化而變化的性質(zhì),叫做氣體的壓縮性,它是氣體的重要屬性,通常用馬赫數(shù)來判別流動(dòng)氣體是否可壓縮[20-22]。馬赫數(shù)是流場中任意一點(diǎn)的速度v與該點(diǎn)處聲速U的比值,用Ma表示,表達(dá)式為

    (8)

    式中γ——比熱容

    M——?dú)怏w常數(shù),J/(kg·K)

    T——熱力學(xué)溫度,K

    當(dāng)Ma≤0.3時(shí),認(rèn)為流體為低速流動(dòng),密度變化忽略不計(jì),流體為不可壓縮流體;當(dāng)Ma>0.3時(shí),認(rèn)為流體為高速流動(dòng),密度變化不可忽略,考慮流體可壓縮性的影響。高速流動(dòng)也可分為3種情況:當(dāng)0.35時(shí),為超高聲速流動(dòng)。

    在本次模擬仿真中,通過試驗(yàn)測(cè)量風(fēng)速最大不超過80 m/s,常溫下氣體比熱容為1.4,氣體常數(shù)為287 J/(kg·K),絕對(duì)溫度為293 K,通過計(jì)算馬赫數(shù)小于0.3,因此氣體看作不可壓縮流體。

    根據(jù)雷諾數(shù)Re來判斷流場類型,通常情況下:當(dāng)Re<2 300時(shí),一般為層流;當(dāng)Re>4 000時(shí), 一般為湍流;當(dāng)2 300≤Re≤4 000可能是層流,也可能是湍流,與流動(dòng)環(huán)境有關(guān)。對(duì)于內(nèi)流問題,當(dāng)Re>2 300時(shí),應(yīng)用湍流模型進(jìn)行分析。雷諾數(shù)為

    (9)

    式中a——粘性系數(shù),Pa·s

    D——水力直徑,m

    標(biāo)準(zhǔn)大氣壓下,空氣密度為1.2 kg/m3,粘性系數(shù)為1.8×10-5Pa·s,通過計(jì)算排種器內(nèi)流場的雷諾數(shù)遠(yuǎn)大于2 300,因此確定流場為湍流。

    流體流動(dòng)遵循質(zhì)量守恒、動(dòng)量守恒以及能量守恒,分別對(duì)應(yīng)計(jì)算流體力學(xué)中的連續(xù)性方程、N-S方程和能量守恒方程,湍流流動(dòng)還需滿足湍流輸運(yùn)方程。

    連續(xù)性方程

    (10)

    N-S方程

    (11)

    湍流輸運(yùn)方程

    (12)

    (13)

    式中k——湍動(dòng)能

    ε——湍動(dòng)耗散率

    μ——湍流粘度

    Gk——平均速度梯度引起的湍動(dòng)能產(chǎn)生項(xiàng)

    Gb——浮力影響產(chǎn)生的湍動(dòng)能產(chǎn)生項(xiàng)

    YM——可壓縮湍流脈動(dòng)膨脹對(duì)總耗散率系數(shù)

    C1ε、C2ε、C3ε——經(jīng)驗(yàn)常數(shù)

    σk、σε——湍動(dòng)能和湍動(dòng)耗散率對(duì)應(yīng)的普朗特?cái)?shù)

    Si、Sk、Sε——源項(xiàng)

    由于本文不涉及傳質(zhì)傳熱相關(guān)內(nèi)容,因此無需考慮能量守恒方程。

    3.2 仿真方法

    對(duì)排種器內(nèi)部流場區(qū)域進(jìn)行網(wǎng)格劃分。采用ANASYS軟件下的MESHING模塊對(duì)排種器進(jìn)行網(wǎng)格劃分,劃分得網(wǎng)格總數(shù)為129 797,最大網(wǎng)格體積為9.49×10-8mm3,最小體積為9.75×10-13mm3。根據(jù)排種器工作原理,選擇排種器進(jìn)種口為壓力入口,落種口為壓力出口,如圖9所示。

    圖9 邊界條件Fig.9 Boundary conditions

    數(shù)值模擬選擇非耦合隱式求解器進(jìn)行求解,采用k-ε模型,工作參考?jí)毫c(diǎn)選擇在進(jìn)種口中心位置,絕對(duì)壓力為101 325 Pa。設(shè)置導(dǎo)流渦輪為Moving wall邊界條件,轉(zhuǎn)速設(shè)為300 r/min,其他壁面采用無滑移邊界條件,流場壓力入口邊界條件為6 kPa壓力入口。為了減少數(shù)值擴(kuò)散,計(jì)算中選取二階迎風(fēng)格式,并采用SIMPLEC算法進(jìn)行求解計(jì)算。

    3.3 結(jié)果分析

    為了研究不同結(jié)構(gòu)導(dǎo)流渦輪對(duì)排種器內(nèi)部流場的影響情況,對(duì)3種不同結(jié)構(gòu)的導(dǎo)流渦輪以及不加導(dǎo)流渦輪的排種器內(nèi)部流場進(jìn)行數(shù)值模擬仿真,探尋影響機(jī)理以及最佳結(jié)構(gòu)。

    圖10中,D型為不增加導(dǎo)流渦輪的排種器,從圖中可以看出,空氣流速較小,總壓強(qiáng)較低,與添加導(dǎo)流渦輪的相比,更加明顯??梢?,導(dǎo)流渦輪可以有效增加排種器內(nèi)部空氣的流動(dòng),大幅增加外圈型孔處空氣流速流量以及落種口空氣的流速,增大壓覆作用力,在排種過程中有助于壓覆充種、高效清種、加速投種。A型有效提高了排種器內(nèi)部流場的流動(dòng)性,具有較好的導(dǎo)流性,將中心空氣導(dǎo)向外圈,增大了外圈型種孔處空氣流速流量,增大壓覆力,與B型和C型相比,由于迎風(fēng)角度較小,擾流性不強(qiáng),外圈空氣流速增大有限。B型有效提高了排種器內(nèi)部流場的流動(dòng)性,具有較好的擾流性,增大了外圈型種孔處空氣流速,增大壓覆力,與A型和C型相比,由于不具備曲線結(jié)構(gòu),導(dǎo)流性不強(qiáng),流場中心空氣向外圍流動(dòng)效果不明顯。C型有效提高了排種器內(nèi)部流場的流動(dòng)性,由于結(jié)構(gòu)兼具較大迎風(fēng)角以及曲線結(jié)構(gòu),因此具有較好的擾流性和導(dǎo)流性,增大了外圈型種孔處空氣流量與流速,增大了壓覆力。3種結(jié)構(gòu)中,C型具有較好的適用性。模擬仿真結(jié)果與上述理論分析一致。

    圖10 腔體內(nèi)部流場速度矢量圖Fig.10 Velocity vector diagrams of flow field inside cavity

    為了定量地分析不同導(dǎo)流渦輪結(jié)構(gòu)對(duì)排種器內(nèi)部流場產(chǎn)生的影響,分別選取充種區(qū)、攜種區(qū)、清種區(qū)以及投種區(qū)型孔處動(dòng)壓為評(píng)價(jià)指標(biāo),如圖11所示。

    圖11 動(dòng)壓變化曲線Fig.11 Changing curves of dynamic pressure

    從圖11可以看出,3種導(dǎo)流渦輪在排種過程的各個(gè)階段均產(chǎn)生動(dòng)壓變化,即產(chǎn)生壓覆作用力,并且變化趨勢(shì)一致。相同工作條件下,其中導(dǎo)流渦輪C在充種、攜種、清種和投種環(huán)節(jié)所產(chǎn)出的壓覆作用力均大于其他2個(gè)。由此可見,導(dǎo)流渦輪C結(jié)構(gòu)最佳。

    4 試驗(yàn)

    4.1 試驗(yàn)條件

    試驗(yàn)所用排種器為自主研發(fā)的氣送式高速玉米精量排種器,所用檢測(cè)裝置為中國農(nóng)業(yè)大學(xué)自主研發(fā)的排種器性能檢測(cè)儀,如圖12所示。試驗(yàn)時(shí),檢測(cè)儀風(fēng)機(jī)為排種器提供風(fēng)壓,將定量排種裝置排出的種子經(jīng)由供料管吹入氣送式高速玉米精量排種器,導(dǎo)種管安裝于排種器投種口下方,用于檢測(cè)排種器性能的傳感器安裝在導(dǎo)種管中間。作業(yè)時(shí),當(dāng)種子經(jīng)過導(dǎo)種管時(shí),觸發(fā)傳感器,檢測(cè)儀通過記錄相鄰種子間的時(shí)間間隔并同時(shí)將其轉(zhuǎn)換為實(shí)際株距的方法,計(jì)算和判斷排種情況,并將計(jì)算的重播指數(shù)、漏播指數(shù)和合格指數(shù)直接輸出到顯示器上。該檢測(cè)儀的檢測(cè)結(jié)果與美國 Precision Planting研發(fā)的MeterMax檢測(cè)儀相比合格指數(shù)誤差不超過 0.7%[23-24]。

    圖12 試驗(yàn)裝置實(shí)物圖Fig.12 Experiment equipment1.排種器性能檢測(cè)儀 2.定量排種裝置 3.供料管 4.排種器 5.導(dǎo)種管

    4.2 試驗(yàn)方法

    試驗(yàn)種子為河南省農(nóng)科院種業(yè)有限公司生產(chǎn)的秋樂牌鄭單958玉米雜交種,籽粒黃色,半馬齒型,千粒質(zhì)量351 g,含水率12.5%,未分級(jí)[25]。根據(jù)GB/T 6973—2005《單粒(精密)播種機(jī)試驗(yàn)方法》的規(guī)定,試驗(yàn)指標(biāo)包括合格指數(shù)Y1=n1/N′×100%;重播指數(shù)Y2=n2/N′×100%;漏播指數(shù)Y3=n3/N′×100%。其中n1為單粒排種數(shù),n2是2粒以上排種數(shù),n3為漏排種數(shù),N′是理論排種數(shù)。根據(jù)前期試驗(yàn)研究,結(jié)合相關(guān)學(xué)者的研究成果[26-29],以影響排種性能的主要參數(shù)工作速度X1、種子喂入量X2以及氣送風(fēng)壓X3為主要試驗(yàn)因素。其中種子喂入量為試驗(yàn)開始階段前5 s的種子喂入速度,隨后喂入量與當(dāng)前工作速度下排種量保持一致,排種量與工作速度關(guān)系式為

    (14)

    式中Q——排種量,g/sS——理論株距,m

    Z——排種器工作速度,m/s

    連續(xù)記錄200粒種子的試驗(yàn)結(jié)果,每次試驗(yàn)重復(fù)3次取平均值。

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

    通過上述理論分析以及仿真分析可以看出,在高速作業(yè)下導(dǎo)流渦輪將會(huì)產(chǎn)生壓覆作用力,結(jié)合前期單因素試驗(yàn)選取工作速度取值范圍在8.3~15.8 km/h,種子喂入量取值范圍在1~2 kg/min,氣送風(fēng)壓范圍在6.5~8.5 kPa。為了找到這3個(gè)因素的最佳參數(shù),使得此排種器排種性能最佳,選擇試驗(yàn)次數(shù)少,計(jì)算方便,可以避免回歸系數(shù)間相關(guān)性的二次旋轉(zhuǎn)正交組合試驗(yàn)方法[30],試驗(yàn)因素和水平如表1所示。再根據(jù)三因素二次旋轉(zhuǎn)正交組合試驗(yàn)表進(jìn)行試驗(yàn),每組試驗(yàn)重復(fù)3次取平均值。試驗(yàn)方案與試驗(yàn)結(jié)果見表2,x1、x2、x3為因素編碼值。

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

    表2 試驗(yàn)設(shè)計(jì)方案與結(jié)果Tab.2 Experiment design and response values

    采用回歸方程方差分析法進(jìn)一步分析, 結(jié)果如表3所示。

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

    注:*表示差異顯著(P<0.05),** 表示差異極顯著(P<0.01)。

    4.4 試驗(yàn)結(jié)果回歸分析

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

    4.4.1合格指數(shù)

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

    (15)

    Y1=35.52+5.10X1+42.81X2-4.02X1X2

    (16)

    通過對(duì)式(16)回歸系數(shù)的檢驗(yàn)得出,影響排種合格指數(shù)的因素主次順序?yàn)榉N子喂入量、工作速度和氣送風(fēng)壓。

    4.4.2重播指數(shù)

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

    (17)

    (18)

    通過對(duì)式(18)回歸系數(shù)的檢驗(yàn)得出,影響重播指數(shù)的因素主次順序?yàn)榉N子喂入量、工作速度和氣送風(fēng)壓。

    4.4.3漏播指數(shù)

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

    (19)

    通過對(duì)式(20)回歸系數(shù)的檢驗(yàn)得出,影響漏播指數(shù)的因素主次順序?yàn)榉N子喂入量、氣送風(fēng)壓和工作速度。

    4.5 各因素對(duì)排種合格指數(shù)的影響

    由圖13可知:排種合格指數(shù)隨著工作速度的提高先上升后下降并保持穩(wěn)定,在10.8~13.7 km/h有最大值;隨著種子喂入量的增大先增后降,在1.2~1.5 kg/min有最大值;隨著氣送風(fēng)壓的增大先增后降,在7.5~8 kPa有最大值。

    圖13 因素交互作用對(duì)合格指數(shù)的影響Fig.13 Impact of interaction on qualified index

    由上述內(nèi)容可知,各因素以及各因素之間的交互作用對(duì)排種器工作性能影響較大,為了獲得最佳排種作業(yè)參數(shù),優(yōu)化工作參數(shù),實(shí)現(xiàn)因素間參數(shù)的合理匹配是提高合格指數(shù)的關(guān)鍵[31-33]。

    以合格指數(shù)為最終優(yōu)化目標(biāo),根據(jù)JB/T 10293—2001《單粒(精密)播種機(jī)技術(shù)條件》中的要求,播種合格指數(shù)大于等于80.0%, 重播指數(shù)小于等于15.0%,漏播指數(shù)小于等于8.0%。結(jié)合各因素邊界條件, 建立參數(shù)化數(shù)學(xué)模型,進(jìn)行優(yōu)化求解,其目標(biāo)函數(shù)和約束條件為

    Fmax=Y1-Y2-Y3

    (21)

    (22)

    Fmax為最終優(yōu)化的目標(biāo),即合格指數(shù);運(yùn)用Design-Expert 8.0.6軟件的優(yōu)化模塊,對(duì)約束目標(biāo)優(yōu)化求解,得到工作速度為9.8 km/h,種子喂入量為1.8 kg/min,氣送風(fēng)壓為8 kPa時(shí),排種合格指數(shù)最高。

    根據(jù)優(yōu)化得到的最優(yōu)參數(shù),進(jìn)行5次重復(fù)對(duì)比臺(tái)架試驗(yàn),對(duì)比結(jié)果如表4所示。對(duì)比試驗(yàn)結(jié)果表明,通過安裝導(dǎo)流渦輪可以大幅度提高合格指數(shù),降低漏播指數(shù)和重播指數(shù)。試驗(yàn)時(shí)發(fā)現(xiàn),在當(dāng)前工作條件下,裝有導(dǎo)流渦輪一組在充種、攜種過程中種子壓覆效果明顯,清種效果明顯優(yōu)于未安裝導(dǎo)流渦輪一組,投種過程中種子掉落速度明顯增加。因此通過安裝導(dǎo)流渦輪可以提高排種器的排種性能。

    表4 對(duì)比試驗(yàn)結(jié)果Tab.4 Result of contrast test %

    5 結(jié)論

    (1)為了使氣送式高速玉米精量排種器內(nèi)流場分布有序,從而提高排種器工作性能,設(shè)計(jì)了3種結(jié)構(gòu)類型的導(dǎo)流渦輪,通過理論分析得出,具有較大迎風(fēng)角、且具備曲線結(jié)構(gòu)的導(dǎo)流渦輪具有較好的擾動(dòng)性和導(dǎo)流性。

    (2)采用CFD方法模擬分析不同結(jié)構(gòu)導(dǎo)流渦輪對(duì)排種器內(nèi)部流場的影響情況,得出3種結(jié)構(gòu)均可有效提高排種器內(nèi)部空氣的流動(dòng)性,增大外圈型孔處空氣流速,增大壓覆作用力。通過對(duì)比分析排種過程中各個(gè)環(huán)節(jié)動(dòng)壓的變化,得出C型結(jié)構(gòu)效果最佳。

    (3)在排種器內(nèi)添加C型導(dǎo)流渦輪,并采用三元二次回歸正交旋轉(zhuǎn)組合試驗(yàn)方法進(jìn)行臺(tái)架試驗(yàn),建立排種性能指標(biāo)(合格指數(shù)、漏播指數(shù)、重播指數(shù))與工作速度、種子喂入量、氣送風(fēng)壓 3 個(gè)試驗(yàn)因素間的回歸數(shù)學(xué)模型,得出了影響各指標(biāo)的因素主次順序, 并運(yùn)用響應(yīng)面法分析了試驗(yàn)因素對(duì)響應(yīng)指標(biāo)的影響。采用多目標(biāo)優(yōu)化方法,確定了排種器最佳工作參數(shù)組合為工作速度為9.8 km/h、種子喂入量為1.8 kg/min、氣送風(fēng)壓為8 kPa,此時(shí)排種合格指數(shù)最高,其性能指標(biāo)為:合格指數(shù) 91.32%、漏播指數(shù) 2.83%、重播指數(shù)5.85%。對(duì)優(yōu)化結(jié)果進(jìn)行對(duì)比驗(yàn)證試驗(yàn),在相同條件下與未安裝導(dǎo)流渦輪排種器進(jìn)行對(duì)比表明,安裝導(dǎo)流渦輪可以有效提高排種器工作性能。

    猜你喜歡
    型孔種器外圈
    預(yù)切種式木薯排種機(jī)構(gòu)設(shè)計(jì)與試驗(yàn)
    深溝球軸承外圈表面凹坑缺陷分析
    哈爾濱軸承(2022年1期)2022-05-23 13:13:16
    組合型孔輪式玉米精量穴播器設(shè)計(jì)與試驗(yàn)
    基于EDEM的雙腔式棉花精量排種器排種性能仿真研究
    一種新式高調(diào)油動(dòng)機(jī)密封圈冷卻系統(tǒng)
    角接觸球軸承外圈鎖口高度自動(dòng)檢測(cè)規(guī)改進(jìn)
    哈爾濱軸承(2020年2期)2020-11-06 09:22:34
    精量排種器現(xiàn)狀及發(fā)展分析
    氣力托勺式馬鈴薯精量排種器設(shè)計(jì)
    旱作水稻鋪膜播種機(jī)排種器性能優(yōu)化試驗(yàn)
    變?nèi)萘啃涂纵喪脚欧N器實(shí)現(xiàn)了播量無級(jí)調(diào)節(jié)
    科技傳播(2011年4期)2011-08-15 00:52:53
    热99re8久久精品国产| 在线国产一区二区在线| 人人妻人人澡人人爽人人夜夜 | 成人无遮挡网站| 日日摸夜夜添夜夜爱| 狂野欧美白嫩少妇大欣赏| 久99久视频精品免费| 亚洲美女搞黄在线观看| 看免费成人av毛片| 一区福利在线观看| 午夜福利在线观看吧| 女的被弄到高潮叫床怎么办| 欧美不卡视频在线免费观看| 男女下面进入的视频免费午夜| 青青草视频在线视频观看| 51国产日韩欧美| 青春草视频在线免费观看| 国产精品久久久久久久电影| 中文精品一卡2卡3卡4更新| 午夜福利在线观看吧| 亚洲久久久久久中文字幕| 国产精品不卡视频一区二区| 日韩欧美精品v在线| 变态另类丝袜制服| 麻豆成人午夜福利视频| 久久久久久久久中文| 国内少妇人妻偷人精品xxx网站| 国产成人精品婷婷| 久久精品夜色国产| 麻豆成人av视频| 99国产极品粉嫩在线观看| 国产av不卡久久| 99久久精品国产国产毛片| 中文在线观看免费www的网站| 亚洲av二区三区四区| 看非洲黑人一级黄片| 综合色丁香网| 欧美成人精品欧美一级黄| 国产精品99久久久久久久久| 久久精品国产自在天天线| 日本免费a在线| 久久精品夜色国产| 不卡视频在线观看欧美| 男女下面进入的视频免费午夜| 99国产精品一区二区蜜桃av| 内地一区二区视频在线| 在线观看免费视频日本深夜| 日韩三级伦理在线观看| 久久九九热精品免费| 国产在线精品亚洲第一网站| 精品人妻偷拍中文字幕| 国产精品精品国产色婷婷| 热99re8久久精品国产| 亚洲人成网站在线观看播放| 日本五十路高清| 欧美性猛交╳xxx乱大交人| 非洲黑人性xxxx精品又粗又长| 午夜视频国产福利| 91aial.com中文字幕在线观看| 美女xxoo啪啪120秒动态图| 夜夜夜夜夜久久久久| 亚洲精华国产精华液的使用体验 | 91精品国产九色| 日韩国内少妇激情av| av国产免费在线观看| 日日干狠狠操夜夜爽| 免费看日本二区| 熟妇人妻久久中文字幕3abv| 婷婷六月久久综合丁香| 久久久色成人| 国产在线精品亚洲第一网站| 久久欧美精品欧美久久欧美| 亚洲精品粉嫩美女一区| 啦啦啦啦在线视频资源| 国语自产精品视频在线第100页| 热99在线观看视频| 久久久久网色| 在线国产一区二区在线| av天堂在线播放| 中国国产av一级| 久久精品久久久久久久性| 亚洲,欧美,日韩| 日本色播在线视频| 插阴视频在线观看视频| 久久久精品大字幕| 亚洲精品亚洲一区二区| 少妇裸体淫交视频免费看高清| 又爽又黄无遮挡网站| 久久精品国产亚洲av涩爱 | 精品久久久久久久久久免费视频| 91久久精品电影网| 91久久精品国产一区二区三区| 青春草国产在线视频 | 麻豆国产97在线/欧美| 成人性生交大片免费视频hd| 特级一级黄色大片| 国产精品一及| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 亚洲欧美成人综合另类久久久 | 真实男女啪啪啪动态图| 久久精品人妻少妇| av视频在线观看入口| 亚洲美女搞黄在线观看| 国内少妇人妻偷人精品xxx网站| 国产视频首页在线观看| 91精品国产九色| 久久久久久久久久久免费av| 久久人人爽人人爽人人片va| 麻豆av噜噜一区二区三区| 日韩国内少妇激情av| 亚洲欧美成人综合另类久久久 | 成人美女网站在线观看视频| 国产乱人视频| 一个人免费在线观看电影| 午夜精品国产一区二区电影 | 欧美zozozo另类| 亚洲第一电影网av| 91午夜精品亚洲一区二区三区| 免费看美女性在线毛片视频| 国产精品久久电影中文字幕| 22中文网久久字幕| 一夜夜www| 亚洲不卡免费看| 国产熟女欧美一区二区| 精品国产三级普通话版| 日韩亚洲欧美综合| 精品少妇黑人巨大在线播放 | 好男人视频免费观看在线| 国产亚洲av片在线观看秒播厂 | 国产成人午夜福利电影在线观看| 亚洲美女视频黄频| 日韩高清综合在线| 在线国产一区二区在线| 国产精品久久久久久精品电影| 给我免费播放毛片高清在线观看| 久久久久久久久久成人| www.色视频.com| 在现免费观看毛片| 麻豆成人av视频| 日韩制服骚丝袜av| 美女 人体艺术 gogo| 国产探花极品一区二区| 美女国产视频在线观看| a级毛色黄片| 三级经典国产精品| 久久久久久大精品| 色尼玛亚洲综合影院| 天堂中文最新版在线下载 | 青春草视频在线免费观看| 午夜福利在线在线| 亚洲精品久久国产高清桃花| 久久久久久久久中文| 五月伊人婷婷丁香| 老司机影院成人| 少妇的逼水好多| eeuss影院久久| 一边摸一边抽搐一进一小说| 欧美高清成人免费视频www| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 九草在线视频观看| 99在线人妻在线中文字幕| 99热这里只有是精品在线观看| 国产精品.久久久| 亚洲av免费在线观看| 精品人妻熟女av久视频| 看十八女毛片水多多多| 丰满乱子伦码专区| 亚洲国产精品国产精品| av天堂在线播放| av卡一久久| 99热精品在线国产| 男人舔女人下体高潮全视频| 国产成人aa在线观看| 成人高潮视频无遮挡免费网站| 我的女老师完整版在线观看| 欧美丝袜亚洲另类| 夜夜看夜夜爽夜夜摸| 欧美一区二区国产精品久久精品| 又爽又黄a免费视频| 国产爱豆传媒在线观看| 性插视频无遮挡在线免费观看| 免费人成在线观看视频色| 亚洲欧美精品专区久久| 亚洲国产色片| 亚洲av二区三区四区| 色哟哟·www| 久久草成人影院| 夜夜看夜夜爽夜夜摸| 乱人视频在线观看| 99热这里只有是精品50| 26uuu在线亚洲综合色| 99久国产av精品国产电影| 麻豆成人av视频| 久久中文看片网| 久久人人爽人人片av| 波野结衣二区三区在线| 深爱激情五月婷婷| 一级毛片aaaaaa免费看小| 欧美日本亚洲视频在线播放| 99在线视频只有这里精品首页| 国产伦一二天堂av在线观看| 免费电影在线观看免费观看| 如何舔出高潮| 三级国产精品欧美在线观看| 国产成人a∨麻豆精品| 免费一级毛片在线播放高清视频| 久久精品国产清高在天天线| 中出人妻视频一区二区| 男女啪啪激烈高潮av片| 最近手机中文字幕大全| 久久久久国产网址| 中文字幕熟女人妻在线| 日韩大尺度精品在线看网址| 亚洲av.av天堂| 秋霞在线观看毛片| 亚洲成人久久性| 成人午夜高清在线视频| 亚洲精华国产精华液的使用体验 | 赤兔流量卡办理| 一区福利在线观看| 干丝袜人妻中文字幕| 夜夜看夜夜爽夜夜摸| 综合色丁香网| 国产精品人妻久久久久久| 九九久久精品国产亚洲av麻豆| 午夜激情福利司机影院| 色哟哟哟哟哟哟| 免费观看人在逋| 国产在线精品亚洲第一网站| 久久久成人免费电影| 国产精品美女特级片免费视频播放器| 国产亚洲精品久久久久久毛片| 成人漫画全彩无遮挡| 真实男女啪啪啪动态图| 久久婷婷人人爽人人干人人爱| 久久久国产成人免费| 久久久久久久久久久免费av| 亚洲最大成人av| 午夜久久久久精精品| 国产三级在线视频| 99在线视频只有这里精品首页| 白带黄色成豆腐渣| 亚洲av第一区精品v没综合| 午夜精品一区二区三区免费看| av国产免费在线观看| 日韩 亚洲 欧美在线| 国产精品一区二区在线观看99 | 国产精品无大码| 午夜精品在线福利| 白带黄色成豆腐渣| 69人妻影院| 99精品在免费线老司机午夜| 久久中文看片网| 极品教师在线视频| 久久99蜜桃精品久久| 一本一本综合久久| 波野结衣二区三区在线| 亚洲aⅴ乱码一区二区在线播放| 午夜视频国产福利| 床上黄色一级片| 蜜桃久久精品国产亚洲av| 国产一级毛片七仙女欲春2| 久久精品夜夜夜夜夜久久蜜豆| 日韩,欧美,国产一区二区三区 | 国产精品国产高清国产av| 色吧在线观看| 51国产日韩欧美| 国产精品野战在线观看| 男的添女的下面高潮视频| 九九热线精品视视频播放| 国产精品久久久久久久久免| 亚洲电影在线观看av| or卡值多少钱| 中出人妻视频一区二区| 国内精品美女久久久久久| 国产成人精品婷婷| 中文字幕精品亚洲无线码一区| 中国美白少妇内射xxxbb| 久久精品国产自在天天线| 国产淫片久久久久久久久| 一级毛片aaaaaa免费看小| 亚洲国产精品sss在线观看| 亚洲高清免费不卡视频| 国产午夜精品一二区理论片| 久久国产乱子免费精品| 久久6这里有精品| 欧美高清性xxxxhd video| 国产高清不卡午夜福利| 亚洲内射少妇av| 国产单亲对白刺激| 亚洲经典国产精华液单| 亚洲欧洲日产国产| 欧美色欧美亚洲另类二区| av黄色大香蕉| 免费人成在线观看视频色| 色5月婷婷丁香| 亚洲真实伦在线观看| 综合色丁香网| 韩国av在线不卡| 91狼人影院| 亚洲国产精品sss在线观看| 亚洲欧美日韩卡通动漫| 一本精品99久久精品77| 身体一侧抽搐| 欧美性感艳星| 日本黄色片子视频| 国产免费男女视频| 午夜a级毛片| 99久国产av精品国产电影| 国产午夜精品一二区理论片| eeuss影院久久| 高清毛片免费观看视频网站| 精品久久久久久久久久免费视频| 12—13女人毛片做爰片一| 久久久精品94久久精品| 国产精品久久久久久精品电影小说 | av在线亚洲专区| 日本五十路高清| 国产高清有码在线观看视频| 最近最新中文字幕大全电影3| 日本免费一区二区三区高清不卡| 欧美性猛交╳xxx乱大交人| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久九九精品二区国产| 深爱激情五月婷婷| 国产精品久久久久久久久免| 不卡视频在线观看欧美| 精品人妻视频免费看| 免费观看人在逋| 亚洲av不卡在线观看| 国产精品国产高清国产av| avwww免费| 丝袜美腿在线中文| 久久草成人影院| h日本视频在线播放| www.色视频.com| 国产淫片久久久久久久久| 国产成人freesex在线| 久久精品国产清高在天天线| 51国产日韩欧美| 三级经典国产精品| 国产一区二区在线av高清观看| 99久国产av精品| 身体一侧抽搐| 精品免费久久久久久久清纯| 日韩欧美国产在线观看| av天堂在线播放| 在线观看午夜福利视频| 少妇的逼水好多| 九九久久精品国产亚洲av麻豆| 日韩欧美精品免费久久| 可以在线观看毛片的网站| 日本免费一区二区三区高清不卡| 乱人视频在线观看| 久久人人爽人人片av| 网址你懂的国产日韩在线| 在线观看一区二区三区| 一个人观看的视频www高清免费观看| 亚洲av成人精品一区久久| 久久人人精品亚洲av| 最近手机中文字幕大全| 久久精品国产亚洲av涩爱 | 中文字幕人妻熟人妻熟丝袜美| 大型黄色视频在线免费观看| 少妇猛男粗大的猛烈进出视频 | 欧美+亚洲+日韩+国产| 六月丁香七月| 看黄色毛片网站| 国产成人91sexporn| 亚洲成a人片在线一区二区| 国产白丝娇喘喷水9色精品| av在线播放精品| 九草在线视频观看| 99国产精品一区二区蜜桃av| 夫妻性生交免费视频一级片| 男插女下体视频免费在线播放| 全区人妻精品视频| 日韩亚洲欧美综合| 大型黄色视频在线免费观看| 国产精品一二三区在线看| 自拍偷自拍亚洲精品老妇| 啦啦啦啦在线视频资源| 色噜噜av男人的天堂激情| 女的被弄到高潮叫床怎么办| 99久久精品热视频| 国产高清三级在线| 麻豆国产97在线/欧美| 精品久久久久久久人妻蜜臀av| a级毛片免费高清观看在线播放| 国产精品野战在线观看| 国产视频首页在线观看| 尾随美女入室| 69av精品久久久久久| 国产精品.久久久| 美女内射精品一级片tv| 久久精品91蜜桃| 亚洲四区av| 在线观看免费视频日本深夜| 最近手机中文字幕大全| 亚洲欧美精品自产自拍| 亚洲最大成人手机在线| 全区人妻精品视频| 麻豆国产av国片精品| 天天躁夜夜躁狠狠久久av| 亚洲欧美中文字幕日韩二区| 一进一出抽搐gif免费好疼| 最近手机中文字幕大全| 久久精品91蜜桃| 午夜精品在线福利| 欧美+亚洲+日韩+国产| 久久99热6这里只有精品| 久久精品人妻少妇| 麻豆乱淫一区二区| 欧美日本亚洲视频在线播放| 男人狂女人下面高潮的视频| 欧美日韩精品成人综合77777| 国产v大片淫在线免费观看| 日韩欧美 国产精品| 夜夜爽天天搞| 高清在线视频一区二区三区 | 变态另类丝袜制服| 最近的中文字幕免费完整| 日本黄色视频三级网站网址| 国国产精品蜜臀av免费| 欧美色欧美亚洲另类二区| 中文资源天堂在线| 国产精品无大码| 嫩草影院入口| 国产色婷婷99| 3wmmmm亚洲av在线观看| 99在线视频只有这里精品首页| 人妻制服诱惑在线中文字幕| 少妇的逼好多水| 日本三级黄在线观看| 亚洲精品日韩在线中文字幕 | 国产午夜精品一二区理论片| 久久久久久久久大av| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 女人十人毛片免费观看3o分钟| 99精品在免费线老司机午夜| 高清日韩中文字幕在线| 校园春色视频在线观看| 久久久午夜欧美精品| 毛片一级片免费看久久久久| 国产精品野战在线观看| 国产不卡一卡二| 国产成人a∨麻豆精品| kizo精华| 18禁黄网站禁片免费观看直播| 性色avwww在线观看| 中国美白少妇内射xxxbb| 村上凉子中文字幕在线| 亚洲精品久久国产高清桃花| 国产精品99久久久久久久久| 亚洲18禁久久av| 国产91av在线免费观看| 国内少妇人妻偷人精品xxx网站| 男女视频在线观看网站免费| 看黄色毛片网站| 精品午夜福利在线看| 精品久久久噜噜| 国产视频首页在线观看| 免费无遮挡裸体视频| 给我免费播放毛片高清在线观看| 国产极品天堂在线| 欧美极品一区二区三区四区| 亚洲18禁久久av| 91精品一卡2卡3卡4卡| 一区二区三区四区激情视频 | 国产美女午夜福利| 中文欧美无线码| 直男gayav资源| 国产精品一区二区三区四区久久| 亚洲无线观看免费| 性插视频无遮挡在线免费观看| 亚洲一区二区三区色噜噜| 久久精品国产亚洲网站| 久久精品国产亚洲av香蕉五月| 亚洲三级黄色毛片| 偷拍熟女少妇极品色| 国产黄片美女视频| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 黄色配什么色好看| 亚洲经典国产精华液单| 超碰av人人做人人爽久久| 欧美高清成人免费视频www| 九色成人免费人妻av| 在线国产一区二区在线| 69人妻影院| 最近最新中文字幕大全电影3| 国产精品福利在线免费观看| 国产一级毛片在线| 一区福利在线观看| 天堂影院成人在线观看| 国产精品一区二区三区四区久久| 久久久久久久久中文| 青青草视频在线视频观看| 亚洲欧美日韩东京热| 亚洲婷婷狠狠爱综合网| 午夜福利在线观看免费完整高清在 | 一边摸一边抽搐一进一小说| 22中文网久久字幕| 青春草国产在线视频 | 国产中年淑女户外野战色| av国产免费在线观看| 欧美高清性xxxxhd video| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲美女视频黄频| 国产日韩欧美在线精品| 美女xxoo啪啪120秒动态图| 亚洲av二区三区四区| 午夜精品在线福利| 精品人妻偷拍中文字幕| 亚洲av不卡在线观看| 久久久久久久久久黄片| 国产精品不卡视频一区二区| 性插视频无遮挡在线免费观看| 国产免费一级a男人的天堂| 国产69精品久久久久777片| 菩萨蛮人人尽说江南好唐韦庄 | 久久鲁丝午夜福利片| 午夜福利视频1000在线观看| 亚洲精品亚洲一区二区| 欧美高清性xxxxhd video| 久久人人爽人人爽人人片va| 搞女人的毛片| 国产亚洲av嫩草精品影院| 日产精品乱码卡一卡2卡三| 亚洲不卡免费看| 亚洲内射少妇av| 国产伦一二天堂av在线观看| 99riav亚洲国产免费| av天堂在线播放| 色尼玛亚洲综合影院| 男女边吃奶边做爰视频| 亚洲国产精品合色在线| 国产老妇伦熟女老妇高清| 国产一区二区三区av在线 | 狂野欧美激情性xxxx在线观看| 成年版毛片免费区| 精品人妻熟女av久视频| 国产一区二区在线观看日韩| 免费观看精品视频网站| 一夜夜www| 成年女人看的毛片在线观看| 美女高潮的动态| 国产成年人精品一区二区| av在线亚洲专区| 毛片女人毛片| 国产真实伦视频高清在线观看| 最近最新中文字幕大全电影3| 一级毛片久久久久久久久女| 亚洲国产高清在线一区二区三| 麻豆乱淫一区二区| 丰满的人妻完整版| 亚洲精华国产精华液的使用体验 | 日本免费一区二区三区高清不卡| 精品99又大又爽又粗少妇毛片| 床上黄色一级片| 国产精品国产高清国产av| 久久久国产成人免费| 五月伊人婷婷丁香| 联通29元200g的流量卡| 一区二区三区四区激情视频 | 免费看日本二区| 一边亲一边摸免费视频| 国产精品一区二区三区四区免费观看| 国产美女午夜福利| 婷婷六月久久综合丁香| 午夜免费男女啪啪视频观看| 中文字幕免费在线视频6| 亚洲欧洲日产国产| 网址你懂的国产日韩在线| 国产黄色视频一区二区在线观看 | 能在线免费看毛片的网站| 日本五十路高清| 国产综合懂色| 免费观看人在逋| 一区二区三区高清视频在线| 欧美另类亚洲清纯唯美| 欧美一区二区国产精品久久精品| 女同久久另类99精品国产91| 夜夜夜夜夜久久久久| 久久久久久久久久黄片| 九草在线视频观看| 中文字幕精品亚洲无线码一区| 国产成人91sexporn| 中国美白少妇内射xxxbb| 午夜福利视频1000在线观看| 精品久久久久久久久亚洲| av专区在线播放| 美女内射精品一级片tv| av天堂在线播放| 精品久久久噜噜| 69av精品久久久久久| 欧美+亚洲+日韩+国产| 天堂影院成人在线观看| 69人妻影院| 91午夜精品亚洲一区二区三区| 亚洲国产精品国产精品| 成人漫画全彩无遮挡| 亚洲成人久久性| 久久精品国产亚洲av天美| 舔av片在线| 永久网站在线| 校园春色视频在线观看| 乱人视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久草成人影院| 综合色av麻豆| av.在线天堂|