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

    錐盤排種器離心推送結(jié)構(gòu)設(shè)計與充種機理分析

    2018-04-02 05:55:48李沐桐李天宇官曉東趙廣闊周福君
    農(nóng)業(yè)機械學(xué)報 2018年3期
    關(guān)鍵詞:漏播種器種間

    李沐桐 李天宇 官曉東 趙廣闊 周福君

    (東北農(nóng)業(yè)大學(xué)工程學(xué)院, 哈爾濱 150030)

    0 引言

    精密排種技術(shù)是國內(nèi)外研究人員廣泛關(guān)注的熱點,由于機械式排種器整體上存在作業(yè)速度不高、對種子形狀和尺寸要求嚴格、播種單粒率有待改善等問題[1],近年來,相關(guān)領(lǐng)域?qū)W者開始對氣力式精密排種器進行改進與研究。通過目前的研究發(fā)現(xiàn),氣力式排種器在適應(yīng)性、效率和均勻性上皆優(yōu)于傳統(tǒng)機械式排種器,然而,其密封結(jié)構(gòu)復(fù)雜,在氣流腔出現(xiàn)縫隙導(dǎo)致氣壓流失等問題出現(xiàn)時,極易發(fā)生漏播現(xiàn)象;同時在田間非平整土壤的高速作業(yè)中,若氣流壓力不足,震動加劇會造成種粒脫離吸盤產(chǎn)生嚴重的重播、漏播等現(xiàn)象;同時,其氣壓由拖拉機動力輸出軸帶動的單離心風(fēng)機產(chǎn)生,隨播種作業(yè)幅寬的增加能耗增加明顯。

    目前,國內(nèi)有關(guān)排種器充排方式的研究報道較少,相關(guān)研究多為機理分析與優(yōu)化改進等[2-5],在國外,2007年德國Amazone公司的新產(chǎn)品DEX氣力集排式播種機作業(yè)速度可達15 km/h[6],其采用正壓滾筒集中充種排種后通過氣送投種,工作幅寬可達9 m,雖然作業(yè)能耗高,但符合國外地廣人稀的大面積作業(yè)模式;在此基礎(chǔ)上,祁兵[7]提出了擾種氣流結(jié)構(gòu)以提高種粒群流動性,進一步對其結(jié)構(gòu)參數(shù)與性能進行了初步研究和優(yōu)化試驗。

    綜合上述問題,本文提出一種立式離心錐盤排種器。相比傳統(tǒng)精播技術(shù)方式:籽粒靜態(tài)堆積—取種—清種—排種—導(dǎo)種,本方式的動作流程為:籽粒群進入錐向離心盤后被賦予持續(xù)動能并具有與錐盤相同周向角速度,通過與排種盤摩擦、動錐盤推送、區(qū)域波紋振動作用,依次有序地以相同姿態(tài)排出,從而在大幅度提高作業(yè)速度的同時,提升排種質(zhì)量穩(wěn)定性,滿足精密播種要求。

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

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

    結(jié)構(gòu)設(shè)計思路是使籽粒群進入充種腔內(nèi)被賦予與充種槽相同角速度的動量,縮小籽粒與充種槽的相對切向線速度,利用錐形結(jié)構(gòu)周轉(zhuǎn)運動時的離心作用,提高籽粒在法向的填充能力,進而改善排種能力。離心錐盤排種器主要由排種管、導(dǎo)種室、靜止盤、動錐盤、傳動軸、彈性銷、殼體、鏈輪和開口銷部件組成,如圖1所示。

    圖1 離心錐盤排種器結(jié)構(gòu)圖Fig.1 Structure diagrams of centrifugal cone grain metering device1.排種管 2.導(dǎo)種室 3.靜止盤 4.動錐盤 5.傳動軸 6.殼體 7.鏈輪 8.彈性銷 9.銷軸 10.種槽 11.推送片 12.導(dǎo)種口 13.波紋震動凹面 14.護種板

    如圖1a所示,按其結(jié)構(gòu)特點,以靜止盤配裝基準,其左側(cè)依次固裝導(dǎo)種室與排種管,其右側(cè)依次扣裝動錐盤、殼體,其中殼體與靜止盤固裝,動錐盤可轉(zhuǎn)動地同軸配裝于靜止盤與殼體之間,動錐盤內(nèi)表面與靜止盤右側(cè)表面呈常接觸閉合狀態(tài)。

    1.2 工作原理

    如圖1b所示,作業(yè)時,外輸入動力通過鏈輪驅(qū)動動錐盤相對于靜止盤逆時針轉(zhuǎn)動,籽粒從導(dǎo)種室進入到動錐盤與靜止盤之間的腔室內(nèi),在充種區(qū)域內(nèi),由動錐盤內(nèi)壁的推送片帶動產(chǎn)生的離心力、籽粒群法向壓力和自身重力三者作用下,籽粒沿法向方向逐一充入動錐盤內(nèi)均布的種槽內(nèi),在經(jīng)過清種波紋凹面振動后,籽粒通過護種板并沿切向方向在左下方所設(shè)投送區(qū)域投入至排種管內(nèi)并排出,依次循環(huán)。

    2 充種力學(xué)分析及相關(guān)參數(shù)設(shè)定

    2.1 無推送結(jié)構(gòu)(靜-動)充種機理

    為進一步探究該型排種器的技術(shù)特點,首先選取無推送片的充種腔室作為對照分析,為清晰區(qū)別于有推送片的離心排種器,按籽粒充種運動狀態(tài),將無推送結(jié)構(gòu)視為靜-動充種機理,其充種力學(xué)分析是進一步研究有推送結(jié)構(gòu)(動-動)充種機理的理論基礎(chǔ)。

    對于錐盤排種器,在籽粒處于下端充種區(qū)域時,其充種利用種間垂直壓應(yīng)力σ1、壁面水平壓應(yīng)力σ2及重力應(yīng)力σG實現(xiàn)法向復(fù)合充填力σ,受力情況如圖2所示。

    圖2 無推送片結(jié)構(gòu)充種機理示意圖Fig.2 Filling stress diagram of cone disc seed-metering device

    2.2 徑向力學(xué)模型

    決定充種性能的關(guān)鍵在于最大化提高充填力,離心錐盤排種器對籽粒充種作用力共3個,即籽粒受上層種群的堆積和擠壓引起的正壓力、籽粒自身自重和錐盤內(nèi)壁面產(chǎn)生的反作用力。3個力復(fù)合即為離心錐盤排種器的充填力。

    由于籽粒在波紋凹面清種區(qū)域內(nèi)主要受自身重力完成清種動作,需留出一定落種空間,本文設(shè)定種層高度為60%,錐盤內(nèi)壁面的反作用力由種間水平側(cè)壓應(yīng)力與自身重力提供,其水平側(cè)壓應(yīng)力分布規(guī)律與物料力學(xué)性質(zhì)及物料堆積深度有關(guān),為簡便計算,將排種器充種腔模型簡化為深倉,并假設(shè)任何水平面上的垂直壓力均為一常數(shù)[8-9],可得

    (1)

    式中dy——微小種層長度,m

    dσ1——種層增加dy厚度時垂直壓應(yīng)力增加量,N

    γ——籽粒重度,玉米為12.25 kN/m3

    fs——物料和內(nèi)壁間靜態(tài)滑動摩擦因數(shù)

    Rh——液力半徑,取0.008 5 m[10]

    當y=0時,σ1=0,積分可得出

    (2)

    式中σ2——種間水平側(cè)壓應(yīng)力,即錐盤內(nèi)壁面產(chǎn)生的反作用力,kPa

    壓力比為

    (3)

    式中φ——籽粒內(nèi)摩擦角,(°)

    玉米φ約為25°[10-11](含水率12%~14%),則K為0.405。

    在安全含水率條件下,玉米籽粒與鋁板之間的靜態(tài)滑動摩擦角φs為31°,則

    fs=tanφs

    (4)

    種槽充填處與種面的距離為

    y=0.6(D-L)

    (5)

    式中D——動錐盤最大外徑,依據(jù)現(xiàn)有排種器尺寸及前期機架空間分布,初步設(shè)計尺寸為215 mm

    L——種槽長度,通過前期針對楔形玉米籽粒的尺寸采集與統(tǒng)計,本文設(shè)計值為14 mm

    0.6為排種器種腔內(nèi)種層高度占總體高度的比值。綜上分析,對于位于充種區(qū)域近最低點的籽粒,復(fù)合充填力σ可表示為

    σ=σGcosθ+σ1cosθ+σ2sinθ

    (6)

    其中

    σG=G/(WT)

    式中σG——籽粒重力應(yīng)力,kPa

    其中W與T分別為楔形籽粒平均寬度與厚度,取W為8.8 mm,T為5.1 mm(本文所研究充種姿態(tài)均為籽粒側(cè)臥于種槽內(nèi),其充種側(cè)臥姿態(tài)宏觀上呈較顯著概率)。

    將式(2)、(4)、(5)代入式(6)計算化簡可得

    σ(θ)=0.055cosθ+0.422(1-e-3.48)cosθ+
    0.171(1-e-3.48)sinθ

    (7)

    應(yīng)用Matlab軟件繪制出籽粒復(fù)合充填應(yīng)力σ(θ)和動錐盤錐角θ的關(guān)系曲線,如圖3所示。

    圖3 充填應(yīng)力函數(shù)曲線Fig.3 Filling force function line

    由圖3可知,錐角區(qū)域[0°,23°]為填充應(yīng)力較大區(qū)間,在17.2°時復(fù)合充填應(yīng)力σ(θ)達到最大,為0.472 kPa。隨著錐角θ的增大,種間垂直壓應(yīng)力和籽粒自身重力的有效分力比例減小,壁面反作用力有效分力比例增大,由于結(jié)構(gòu)特點,水平方向的應(yīng)力在充種過程中的作用并不大,同時錐角的增大會導(dǎo)致容腔寬度增加,可知錐角在[10°,220°]范圍內(nèi)變化不大,為便于加工,初步設(shè)定錐角為20°。

    2.3 周向力學(xué)模型

    上述分析是基于籽粒處在種腔最底端時錐角對充填力的影響,而真實情況下,隨著排種器內(nèi)動錐盤的轉(zhuǎn)動,種槽由最低端的豎直位置轉(zhuǎn)動到水平位置,在此過程中處于不同角度位置的籽粒由于有效種層高度逐漸降低,且重力方向與充填方向夾角逐漸增大,導(dǎo)致種間垂直壓應(yīng)力σ1和重力應(yīng)力σG均逐步減小,同時,逐漸增加的種間水平側(cè)壓應(yīng)力σ2對充填力影響并不起關(guān)鍵作用,將造成籽粒的法向充填應(yīng)力σn逐漸減小,使排種器充種區(qū)域縮小至種腔底端附近,充填性能降低,很大程度上限制了排種器充種頻率和工作速度。為改善上述問題,首先需要分析籽粒的有效充填區(qū)域和充填應(yīng)力在不同角度位置的變化規(guī)律,如圖4所示。

    圖4 不同角度位置籽粒法向充填應(yīng)力分布Fig.4 Distribution of grain filling force at different angles

    設(shè)籽粒位置角β為籽粒重心和圓心連線與鉛垂線夾角,由圖4可知,當β=0°時,籽粒法向充填應(yīng)力最大,其主要由重力應(yīng)力σG和種間垂直壓應(yīng)力σ1提供;然后隨籽粒位置角β逐漸增大,種層有效高度降低使中間垂直壓應(yīng)力σ1減小,同時,重力應(yīng)力對于法向充填應(yīng)力σn的貢獻量逐漸降低,使種間水平側(cè)壓應(yīng)力σ2逐步轉(zhuǎn)化為法向充填應(yīng)力的主要作用應(yīng)力,然而,由上述分析可知,種間水平側(cè)壓應(yīng)力σ2不僅隨種間垂直壓應(yīng)力σ1減小而減小,且其數(shù)值上也遠小于種間垂直壓應(yīng)力,因此,參考上節(jié)分析,種槽充填處與種面的距離可整合為

    (8)

    籽粒充種在區(qū)域內(nèi)的充填應(yīng)力數(shù)學(xué)模型可進一步修正為(θ=20°)

    σ(β)=0.052cosβ+0.397[1-e-(3.19cosβ-0.232)]cosβ+
    0.058 5[1-e-(3.19cosβ-0.232)]sinβ

    (9)

    通過Matlab軟件繪制修正后的復(fù)合填充力σ(β)和籽粒位置角β的關(guān)系曲線,如圖5所示。

    圖5 充填應(yīng)力修正函數(shù)曲線Fig.5 Filling force updated function line

    由圖5可知,在無推送片的充種條件下,最適宜籽粒充填的區(qū)域較小,即在籽粒位置角[-5°,16.5°]區(qū)域內(nèi),此角度區(qū)域為籽粒復(fù)合充填應(yīng)力較高范圍,最大值為0.42 kPa,且隨籽粒位置角增大,復(fù)合填充應(yīng)力逐漸減小。綜上分析,對于無推送(靜-動)充種方式,在排種器工作轉(zhuǎn)速提高時,動錐盤與靜態(tài)籽粒間較大的相對線速度在較小的有效充種區(qū)域內(nèi)并不適合高頻率的充種動作。若要提高充種效率,不僅要提高充填應(yīng)力,也要擴大有效充種區(qū)域,最大程度減小動錐盤與預(yù)充填籽粒的相對切向線速度。

    2.4 推送結(jié)構(gòu)(動-動)充種機理

    綜合以上分析,為提高充填應(yīng)力并擴大有效作用區(qū)域,在動錐盤內(nèi)壁設(shè)置推送結(jié)構(gòu),使進入種腔內(nèi)的籽粒群被持續(xù)勻速轉(zhuǎn)動的推送片賦予相同角速度,在充種區(qū)域近種槽位置的籽粒主要由離心壓應(yīng)力σC、重力應(yīng)力σG、種間垂直壓應(yīng)力σ1和種間水平側(cè)壓應(yīng)力σ2的共同作用作為復(fù)合充填力σ。由于不同結(jié)構(gòu)形狀的推送片對種層壓力方向不同,復(fù)合充填力σ的方向也由上方種層壓力傳遞決定,故選取規(guī)則條形推送片作為籽粒充種受力的初步分析對象,如圖6所示。

    圖6 籽粒在推送結(jié)構(gòu)下的法向充填應(yīng)力Fig.6 Normal filling stress of grain with push blade

    由于籽粒群在運動過程中產(chǎn)生的摩擦、碰撞和黏合力非常復(fù)雜,為方便分析,假設(shè)籽粒上層種群之間相互作用力為等效傳遞,將上層隨推送結(jié)構(gòu)運動的種群視為整體質(zhì)點P,其繞軸心O勻速轉(zhuǎn)動,其合力為外層籽粒所提供的方向指向軸心O的向心拉應(yīng)力σN;而對于下層待充填的籽粒而言,法向復(fù)合充填應(yīng)力σ′n主要由持續(xù)的離心壓應(yīng)力σC(σC=σN)、自身重力應(yīng)力σG、種間垂直壓應(yīng)力σ1和種間水平側(cè)壓應(yīng)力σ2共同提供。為研究推送作用下的復(fù)合充填力,需先分析影響σC的因素,進而確定相關(guān)參數(shù)以計算復(fù)合充填力。

    處于勻速圓周運動的上方運動種層主要為相鄰兩推送片之間的環(huán)形部分,其動態(tài)種層質(zhì)點質(zhì)心與軸心O的距離r為

    (10)

    式中RTE——推送片外徑,mm

    RTI——推送片內(nèi)徑,為保證籽粒流動順利且重心傳動裝配方便,本文取40 mm

    α——相鄰兩推送片內(nèi)壁面夾角,(°)

    結(jié)合式(10)并代入σC=Δmω2r/ΔS可得

    (11)

    ΔS=απRTE(TC-ΔT)/180

    (12)

    (13)

    式中Δm——環(huán)形區(qū)域內(nèi)種層質(zhì)量,kg

    ΔS——下層待充填種層內(nèi)環(huán)面面積,m2

    ρg——玉米籽粒密度(安全含水率),本文取1.197 kg/m3

    TC——錐盤種腔內(nèi)壁厚度,取24 mm

    ΔD——錐盤內(nèi)推送片側(cè)向直徑差,本文設(shè)定30 mm

    ΔT——錐盤側(cè)向無效厚度,本文為12 mm

    ω——工作角速度,ω<8.9 rad/s(需保證重力清種ωr

    由式(11)、(13)可看出,RTE越大,σC越大,然而排種器工作過程中既要保證上方有足夠空間清種,又要避免推送片和護板對籽粒造成的擠壓、夾傷和撞擊,在本文中設(shè)定RTE為80 mm。

    由式(11)可知,隨工作轉(zhuǎn)速的增大,離心拉應(yīng)力呈幾何級增長,而相鄰?fù)扑推瑠A角在0°~360°范圍內(nèi),離心拉應(yīng)力呈先增加后降低規(guī)律變化,當夾角約180°時,離心拉應(yīng)力接近最大值。分析以上結(jié)果可知,如果夾角過小,環(huán)面內(nèi)勻速運動種層質(zhì)量變小,同時種層質(zhì)心距軸心半徑r較小,所以離心拉應(yīng)力σC有效作用較小,因此,相鄰兩推送片夾角α設(shè)定為180°,推送片個數(shù)為2個,為發(fā)揮離心拉應(yīng)力的最大效率,需要計算排種器工作角速度的極限取值。

    2.5 充填可靠性分析

    為保證籽粒及時、完整的充填,需對目標籽粒相對于種槽的法向充填時間tn與種槽相對目標籽粒的周向運動時間tτ進行計算比較,以確定種槽寬度及工作轉(zhuǎn)速范圍。如圖7所示,目標籽粒在充填過程中具有與種槽相同的角速度ω,其籽粒法向充填時間為

    (14)

    其中

    an=gcosβ+ω2r+a1+a2

    (15)

    式中S——種槽長度,本文取16 mm

    an——籽粒法向合加速度,m/s2

    r——動態(tài)種層質(zhì)點質(zhì)心與軸心O的距離,由上述分析取0.124 mm

    a1——種間垂直壓力加速度,m/s2

    a2——種間水平側(cè)壓力加速度,m/s2

    種槽相對于籽粒的周向運動時間為

    (16)

    式中β0——種槽內(nèi)兩壁面延長線分別與軸心O連線的夾角,本文取8°

    為實現(xiàn)目標籽粒及時充填且達到完全充填效果,差值T需滿足

    T=tn-tτ≤0

    (17)

    將時間差T作為目標函數(shù),籽粒位置角β和排種器工作角速度ω作為變量,采用ORIGIN軟件生成曲面變化圖,如圖8所示。

    圖8 工作角速度與籽粒位置角對時間差T的影響Fig.8 Influence of grain position angle and working angle velocity on T value

    由圖8可看出,籽粒位置角β在范圍(-30°,60°)內(nèi)對時間差T的影響較小,而工作角速度ω對T影響區(qū)間較大,隨著工作角速度逐漸增大,T由負值逐漸增加,當工作角速度大于2.9 rad/s時,T值大于零,即tn>tτ,在此工況下,籽粒已不能保證完全充填入種槽內(nèi),由此設(shè)定排種器工作轉(zhuǎn)速最大值為27.7 r/min。

    2.6 力學(xué)模型

    通過上述分析,籽粒在有推送裝置的充填應(yīng)力數(shù)學(xué)模型最終確定為(α=180°)

    σ′(β)=0.052cosβ+0.124ω2+
    0.397[1-e-(3.12cosβ-0.41)]cosβ+
    0.058 5[1-e-(3.12cosβ-0.41)]sinβ

    (18)

    (19)

    式中y″——上方有效種層高度,mm

    由于兩推送結(jié)構(gòu)對稱排布,當一端推送片運動到種腔上端時,種層高度已不能通過上述公式計算,因此,結(jié)合設(shè)定的最終推送結(jié)構(gòu)位置特點,并保證籽粒在內(nèi)部有良好的流動性,設(shè)定種層有效高度為50%的種腔高度,以便簡化計算。

    將工作角速度初步設(shè)定5個水平,分別為2.30、2.51、2.72、2.93、3.14 rad/s,對應(yīng)工作轉(zhuǎn)速為22、24、26、28、30 r/min,分別代入式(18)并通過Matlab繪制出籽粒在推送裝置下的復(fù)合充填應(yīng)力變化曲線,如圖9所示。

    圖9 推送結(jié)構(gòu)下復(fù)合充填應(yīng)力的變化曲線Fig.9 Changing curves of filling stress with push structure

    由圖9可知,在設(shè)置離心推送裝置后籽粒復(fù)合充填應(yīng)力有了明顯增加,在所選工作轉(zhuǎn)速范圍內(nèi),籽粒位置角約8°時,復(fù)合充填應(yīng)力均達到最大值;隨籽粒位置角的增大,復(fù)合充填應(yīng)力減小的幅度較小,當工作轉(zhuǎn)速較低時(n=22 r/min),充填應(yīng)力最大值約1.0 kPa;當工作轉(zhuǎn)速較高時(n=30 r/min),充填應(yīng)力最大值約為1.65 kPa,與上述無推送裝置(靜-動)充種方式相比(圖5),安裝離心推送裝置后的充填應(yīng)力提高了近4倍,同時,相比(靜-動)方式的有效充填角度區(qū)域[-5°,16°],其有效充種角度區(qū)域有了明顯的擴大。綜上理論分析,推送裝置的增加可大幅度提高充種性能。

    3 EDEM虛擬仿真

    采用EDEM(Engineering discrete element method)軟件對離心錐盤排種器進行虛擬仿真。依據(jù)上述分析可知,有推送裝置(動-動)離心充種方式相比無推送裝置(靜-動)充種方式在充種效率上有顯著優(yōu)勢,故選取以上兩種結(jié)構(gòu)作為對照,進行虛擬仿真及單體臺架試驗,以充分驗證離心錐盤排種器在設(shè)置推送結(jié)構(gòu)后的充種性能。

    3.1 幾何模型建立

    為便于仿真模擬及計算,將籽粒運動過程中與接觸無關(guān)的部件去除,應(yīng)用Solidworks對排種器進行三維建模,將.igs格式文件導(dǎo)入EDEM軟件中,添加旋轉(zhuǎn)方向和轉(zhuǎn)速,創(chuàng)建虛擬工廠作為籽粒生成的虛擬面,其幾何體模型如圖10所示。

    圖10 帶推送結(jié)構(gòu)離心錐盤排種器仿真模型Fig.10 Simulation model of centrifugal cone seed metering device with push structure

    3.2 定義顆粒模型

    選取國內(nèi)廣泛種植的紫甜糯六號玉米籽粒,其楔狀外形及尺寸與大部分糧食玉米相近,如圖11a所示。隨機挑選100顆籽粒,測量記錄籽粒的長、寬、高平均值,在EDEM軟件中通過多球面組合方式填充,模擬籽粒狀態(tài)如圖11b所示。

    圖11 籽粒外形與顆粒模型Fig.11 Grain shape and particle model

    3.3 虛擬試驗方案

    在仿真過程中,設(shè)定顆粒與顆粒、顆粒與壁面之間均采用Hertz模型,固定時間步長為Rayleigh時間步長的15%,仿真時間為20 s,網(wǎng)格尺寸為顆粒尺寸的0.5倍,所涉及相關(guān)物理參數(shù)如表1所示。虛擬試驗選取排種器工作轉(zhuǎn)速為變量,分別對有推送結(jié)構(gòu)和無推送結(jié)構(gòu)的離心錐盤排種器的工作過程進行仿真,同時,在后處理過程中取充種區(qū)域內(nèi)任意時間的隨機顆粒作為標記檢測點,識別顆粒在充填后與種槽壁面碰撞的瞬時合外力數(shù)值,作為近似充填力以便比較計算。以漏播指數(shù)M為評價指標[12],其工作轉(zhuǎn)速與作業(yè)速度的關(guān)系式為

    (20)

    式中n——排種器工作轉(zhuǎn)速,r/min

    v——作業(yè)速度,km/h

    l——理論株距,本文設(shè)定為240 mm

    MK——動錐盤內(nèi)種槽數(shù)量,可由上文分析確定為30個

    作業(yè)速度取值域為7~13 km/h,對應(yīng)排種器工作轉(zhuǎn)速為16.20~30.09 r/min。仿真過程截圖如圖12所示。

    表1 全局變量參數(shù)設(shè)置Tab.1 Pre-treatment parameters setting

    圖12 排種器仿真過程截圖Fig.12 Screenshot of seed-metering device simulation process

    4 臺架試驗

    為進一步檢驗理論分析與虛擬仿真結(jié)果,以離心錐盤排種器單體為載體,分別配裝帶推送裝置和無推送裝置的動錐盤進行單體臺架驗證試驗。試驗地點在黑龍江省農(nóng)業(yè)機械工程科學(xué)研究院排種器實驗室,玉米籽粒選取紫甜糯6號,所用設(shè)備為JPS-12型排種器性能檢測試驗臺。

    為便于記錄試驗結(jié)果,在試驗過程中卸載排種管,實時觀測投種區(qū)域豁口處的種槽狀態(tài),成功充填的籽粒從種槽投出后落至涂有油層的種床帶上;種槽內(nèi)空漏無籽粒時記錄為一次漏播。漏播指數(shù)計算公式為

    (21)

    式中N——理論充填數(shù)n0——漏播次數(shù)

    5 結(jié)果分析

    在不同工作轉(zhuǎn)速條件下,離心錐盤排種器帶推送結(jié)構(gòu)與無推送結(jié)構(gòu)的漏播指數(shù)仿真值與試驗值對比如表2所示。

    表2 排種器不同工作轉(zhuǎn)速下漏播指數(shù)的仿真值與試驗值Tab.2 Simulated and experimental values of missing rate at different rotational speeds of seed metering device %

    在EDEM仿真后處理結(jié)果中,測定標記顆粒的位置角(約45°),將角度代入式(18)中得出理論充填應(yīng)力;同時提取標記顆粒與種槽壁面碰撞的瞬時合外力,并轉(zhuǎn)化為充填應(yīng)力,繪制出轉(zhuǎn)速在16.2~30.1 r/min內(nèi)的仿真值與理論值的對比圖,如圖13所示。

    圖13 充填應(yīng)力對比Fig.13 Comparison of filling stress

    由表2結(jié)果分析可知,有推送結(jié)構(gòu)與無推送結(jié)構(gòu)的離心錐盤排種器的仿真值與試驗值基本相符,誤差小于4.8%,其產(chǎn)生原因可能是由實際工作中的轉(zhuǎn)速不穩(wěn)定、籽粒尺寸差異造成,且在實際試驗中由于先前使用有機玻璃的破碎情況較多,故本次試驗觀測方式無法直觀記錄清種區(qū)域的籽粒實時狀態(tài),期間發(fā)生的籽粒卡滯、壓迫等問題不能避免。

    當排種器工作轉(zhuǎn)速在16.2~25.5 r/min時,帶推送結(jié)構(gòu)的漏播指數(shù)接近3%,且基本保持穩(wěn)定,當轉(zhuǎn)速大于25.5 r/min時,漏播指數(shù)有明顯升高;在相同轉(zhuǎn)速范圍(16.2~30.1 r/min)內(nèi),無推送結(jié)構(gòu)的漏播指數(shù)隨工作轉(zhuǎn)速的增加而持續(xù)升高,當轉(zhuǎn)速大于18.5 r/min時,仿真值與試驗值均較高,在轉(zhuǎn)速為27.8 r/min時,無推送結(jié)構(gòu)的漏播指數(shù)達到7.99%,接近GB/T 6973—2005《單粒(精密)播種機試驗方法》要求的漏播指數(shù)(小于8%)。由此可見,離心推送結(jié)構(gòu)與無推送結(jié)構(gòu)相比,在充種性能上有明顯改善。

    由圖13可知,推送充種方式的充填應(yīng)力仿真值整體低于理論值,但總體變化趨勢基本一致。其原因可能由于顆粒在充填過程中存在顆粒間無規(guī)律的摩擦、碰撞、擠壓等復(fù)雜接觸,且在運動種槽內(nèi)壁面的過程中存在能量損失導(dǎo)致。而當工作轉(zhuǎn)速超過25.5 r/min時,仿真值的增長趨勢趨于平緩,與理論值的差值逐漸增大,其原因可能由于隨著轉(zhuǎn)速的提高,顆粒與內(nèi)壁存在相對滑動,所需充填時間大于種槽周向運動時間,即不能及時完成顆粒的完整充填,此時顆粒位于種槽內(nèi)的部分與位于種腔的部分產(chǎn)生了相對運動,將造成顆粒與種槽在切向方向連續(xù)的撞擊與震動,導(dǎo)致離心作用不能持續(xù)增加。總體而言,各轉(zhuǎn)速下的充填應(yīng)力接近預(yù)測值,其數(shù)學(xué)模型可用于充種機理研究。

    6 性能驗證試驗

    為進一步驗證帶離心推送結(jié)構(gòu)的錐盤排種器的單體作業(yè)性能,在上述工況下,選取合格指數(shù)、漏播指數(shù)、重播指數(shù)和破碎率為指標,分別對每組進行3次試驗(圖14),其結(jié)果如表3所示。

    圖14 臺架驗證試驗Fig.14 Verification test of bench1.試驗臺架 2.排種器 3.傳動鏈輪 4.驅(qū)動電動機 5.種床帶

    由表3可知,排種器漏播指數(shù)與上述試驗結(jié)果基本一致,重播指數(shù)隨轉(zhuǎn)速的增加表現(xiàn)出降低趨勢,對于合格指數(shù)而言,隨著轉(zhuǎn)速增加呈現(xiàn)出先升高后降低的規(guī)律,在轉(zhuǎn)速分別為20.8、23.1、25.5 r/min的3組試驗中,合格指數(shù)分別為90.2%、91.0%和90.5%;而對于破碎率而言,隨轉(zhuǎn)速的增加,呈現(xiàn)先降低后升高的規(guī)律,造成此結(jié)果的原因可能是在低轉(zhuǎn)速過程中,種槽內(nèi)充入多余的籽粒,在工作過程中的內(nèi)部結(jié)構(gòu)與籽粒碰撞劇烈;而在高轉(zhuǎn)速情況下,內(nèi)部結(jié)構(gòu)與籽粒碰撞頻率提高,其轉(zhuǎn)速大于27.8 r/min時,破碎率達到了0.4%。

    表3 帶離心推送結(jié)構(gòu)的錐盤排種器工作性能試驗結(jié)果Tab.3 Experiment results of working performance of cone disc seeding device with a centrifugal push structure

    7 結(jié)論

    (1)傳統(tǒng)機械式排種器重力充種過程中,籽粒位置角的增大導(dǎo)致充填力直線下降是限制充種效率的關(guān)鍵因素。因此,增加充填力種類和擴大有效充種區(qū)域是提高排種器充種性能的重要途徑。

    (2)相比無推送結(jié)構(gòu)的立式錐盤排種器,在種腔內(nèi)設(shè)置推送結(jié)構(gòu)可充分利用上層種群的離心壓力提高充種性能。通過臺架對比試驗可知,帶推送結(jié)構(gòu)的動錐盤轉(zhuǎn)速可達30.1 r/min,單體作業(yè)速度為13 km/h,在此范圍內(nèi)平均漏播指數(shù)約3.93%,較無推送結(jié)構(gòu)充種方式的平均漏播指數(shù)降低了2.52%。同時,參考文獻[12],在相同試驗條件下,作業(yè)速度試驗值比指夾式排種器提高了3 km/h,且在相同作業(yè)速度(9 km/h)條件下,漏播指數(shù)降低了3.94%。由此可知,離心推送結(jié)構(gòu)有效提高了充填能力,并有效提高了充種頻率,從而提升了作業(yè)速度,且滿足精播技術(shù)要求中對漏播指數(shù)的要求。

    (3)性能驗證試驗表明,當工作轉(zhuǎn)速在20.8~25.5 r/min范圍內(nèi)時,排種器各項性能及指標表現(xiàn)較優(yōu)狀態(tài),而破碎率變化區(qū)間較大,仍需對內(nèi)部材料進一步改善。

    (4)通過對籽粒在立式錐盤排種器中充種機理的分析,研究了影響充填應(yīng)力和有效充填區(qū)域的因素,確定了關(guān)鍵結(jié)構(gòu)與核心參數(shù),其試驗也主要針對充種性能進行檢驗。然而,排種器動錐盤的種槽尺寸、清種裝置結(jié)構(gòu)和參數(shù)、投種區(qū)域位置和排種管最佳正壓值等因素是影響重播指數(shù)、破碎率和粒距變異系數(shù)等指標的重要因素,仍需進一步優(yōu)化相關(guān)參數(shù)、改善材料來提高錐盤排種器整體可靠性。

    1楊麗,顏丙新,張東興,等.玉米精密播種技術(shù)研究進展[J/OL].農(nóng)業(yè)機械學(xué)報,2016,47(11):38-48.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20161106&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.11.006.

    YANG Li,YAN Bingxin,ZHANG Dongxing,et al.Research progress on precision planting technology of maize[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2016,47(11):38-48.(in Chinese)

    2張昆,衣淑娟.氣吸滾筒式玉米排種器充種性能仿真與試驗優(yōu)化[J/OL].農(nóng)業(yè)機械學(xué)報,2017,48(7):78-86.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20170710&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2017.07.010.

    ZHANG Kun,YI Shujuan.Simulation and experimental optimization on filling seeds performance of seed metering device with roller of air-suction[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(7): 78-86.(in Chinese)

    3劉佳,崔濤,張東興,等.氣吹式精密排種器工作壓力試驗研究[J].農(nóng)業(yè)工程學(xué)報,2011,27(12):18-22.

    LIU Jia,CUI Tao,ZHANG Dongxing, et al.Experimental study on pressure of air-blowing precision seed-metering device[J].Transactions of the CSAE,2011,27(12):18-22.(in Chinese)

    4YU H, ZHAO X, LIU Z, et al. Study on negative-pressure precision millet seed-metering device[J]. Advance Journal of Food Science & Technology, 2015, 7(2):139-143.

    5賀俊林,裘祖榮.新型氣壓式精密排種器的試驗研究[J].農(nóng)業(yè)工程學(xué)報,2001,17(2):80-83.

    HE Junlin,QIU Zurong.Experimental study on a new type seed meter[J].Transactions of the CSAE,2001,17(2):80-83.(in Chinese)

    6Werke H Dreyer GmbH&Co.EDX precision airplanter overview [EB/OL].(2016-10-4)[2016-10-4]. http:∥www.amazone.net/1262. asp.

    7祁兵. 中央集排氣送式精量排種器設(shè)計與試驗研究[D].北京.中國農(nóng)業(yè)大學(xué), 2014.

    QI Bing.Design and experiment of acentralized metering device [D].Bejing: China Agricultural University, 2014. (in Chinese)

    8劉宏新,徐曉萌,郭麗峰,等.具有復(fù)合充填力的立式淺盆型排種器充種機理[J].農(nóng)業(yè)工程學(xué)報,2014,30(21):9-16.

    LIU Hongxin,XU Xiaomeng,GUO Lifeng,et al.Research on seed-filling mechanism of vertical shallow basin type seed-metering device with composite filling force[J].Transactions of the CSAE,2014,30(21):9-16.(in Chinese)

    9劉宏新,劉俊孝,唐師法,等.對置斜盤高速精密大豆排種器設(shè)計與充種機理分析[J].農(nóng)業(yè)工程學(xué)報,2016,32(20):24-31.

    LIU Hongxin,LIU Junxiao, TANG Shifa,et al.Design on opposed inclined-plate high-speed precision seed-metering device and its working mechanism analysis[J].Transactions of the CSAE,2016,32(20):24-31. (in Chinese)

    10王云霞,梁志杰,張東興,等.基于離散元的玉米種子顆粒模型種間接觸參數(shù)標定[J].農(nóng)業(yè)工程學(xué)報,2016,32(22):36-42.

    WANG Yunxia,LIANG Zhijie,ZHANG Dongxing,et al.Calibration method of contact characteristic parameters for corn seeds based on EDEM[J] .Transactions of the CSAE,2016,32(22):36-42.(in Chinese)

    11崔濤,劉佳,楊麗,等.基于高速攝像的玉米種子滾動摩擦特性試驗與仿真[J].農(nóng)業(yè)工程學(xué)報,2013,29(15):34-41.

    CUI Tao,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)

    12王金武,唐漢,王奇,等.基于EDEM軟件的指夾式精量排種器排種性能數(shù)值模擬與試驗[J].農(nóng)業(yè)工程學(xué)報,2015,31(21):43-50.

    WANG Jinwu,TANG Han,WANG Qi,et al.Numerical simulation and experiment on seeding performance of pickup finger precision seed-metering device based on EDEM[J]. Transactions of the CSAE,2015,31(21):43-50.(in Chinese)

    13WANG Jinwu, TANG Han, WANG Jinfeng,et al. Measurement and analysis of restitution coefficient between maize seed and soil based on high-speed photography[J]. International Journal of Agricultural and Biological Engineering, 2017,10(3):102-114.

    14WANG Jinwu, TANG Han, WANG Jinfeng,et al. Optimization design and experiment on ripple surface type pickup finger of precision maize seed metering device[J]. International Journal of Agricultural and Biological Engineering, 2017,10(1):61-71.

    15周福君,李沐桐,李天宇.離心錐盤式高速排種單體:201710334326.1 [P].2017-05-12.

    16于建群,王剛,心男,等.型孔輪式排種器工作過程與性能仿真[J].農(nóng)業(yè)機械學(xué)報,2011,42(12):83-87,101.

    YU Jianqun,WANG Gang,XIN Nan,et al.Simulation analysis of working process and performance of cell wheel metering device[J].Transactions of the Chinese Society for Agricultural Machinery,2011,42(12):83-87,101. (in Chinese)

    17廖慶喜, 張朋玲, 廖宜濤,等. 基于EDEM的離心式排種器排種性能數(shù)值模擬[J/OL]. 農(nóng)業(yè)機械學(xué)報, 2014, 45(2):109-114.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20140219&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.02.019.

    LIAO Qingxi, ZHANG Pengling, LIAO Yitao, et al. Numerical simulation on seeding performance of centrifugal rapeseed metering device based on EDEM[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(2):109-114.(in Chinese)

    18劉宏新, 徐曉萌, 劉俊孝,等. 利用高速攝像及仿真分析立式淺盆型排種器工作特性[J]. 農(nóng)業(yè)工程學(xué)報, 2016, 32(2):13-19.

    LIU Hongxin, XU Xiaomeng, LIU Junxiao, et al. Working characteristics of vertical shallow-basin type seed-metering device based on high-speed photography and virtual simulation[J]. Transactions of the CSAE, 2016,32(2):13-19.(in Chinese)

    19YU J J, LIAO Y Y, CONG J L, et al. Simulation analysis and match experiment on negative and positive pressures of pneumatic precision metering device for rapeseed[J]. International Journal of Agricultural and Biological Engineering, 2014, 7(3): 1-12.

    20ZHANG G Z, ZANG Y, LUO X W,et al. Design and indoor simulated experiment of pneumatic rice seed metering deivce[J]. International Journal of Agricultural and Biological Engineering, 2015,8(4): 10-18.

    21YANG S, ZHANG S M. Design and parameter optimization of flexible comb-type grass seed metering device[J]. International Journal of Agricultural and Biological Engineering, 2015,8(1): 9-16.

    猜你喜歡
    漏播種器種間
    玉米免耕播種漏播補償裝置研究
    干旱條件對鬼針草和醉魚草種間相互作用及生長的影響
    植物研究(2023年5期)2023-09-09 08:01:22
    三峽庫區(qū)支流花溪河浮游植物種間關(guān)聯(lián)及影響因子分析
    基于EDEM的雙腔式棉花精量排種器排種性能仿真研究
    籽粒作物漏播補償方法研究現(xiàn)狀分析
    馬鈴薯播種機漏播檢測預(yù)警系統(tǒng)設(shè)計
    精量排種器現(xiàn)狀及發(fā)展分析
    氣力托勺式馬鈴薯精量排種器設(shè)計
    江蘇省宜興市茶園秋季雜草種間生態(tài)關(guān)系及群落分類
    型孔輪式水稻排種器排種監(jiān)控系統(tǒng)設(shè)計
    中文字幕熟女人妻在线| 精品日产1卡2卡| 久久久精品大字幕| 亚洲国产高清在线一区二区三| 欧美在线黄色| 亚洲av熟女| 免费在线观看成人毛片| 国产欧美日韩精品一区二区| 人妻久久中文字幕网| 国产v大片淫在线免费观看| 免费av观看视频| 99久久久亚洲精品蜜臀av| 日本 av在线| 亚洲自偷自拍三级| 噜噜噜噜噜久久久久久91| 最近最新免费中文字幕在线| 成人av在线播放网站| 国产麻豆成人av免费视频| 中文字幕免费在线视频6| 国内精品久久久久精免费| 白带黄色成豆腐渣| 人人妻人人澡欧美一区二区| 国产免费av片在线观看野外av| 亚洲av成人av| 色播亚洲综合网| 好男人在线观看高清免费视频| 非洲黑人性xxxx精品又粗又长| 欧美黄色淫秽网站| 国产色婷婷99| 日韩人妻高清精品专区| 国产精品精品国产色婷婷| 成人性生交大片免费视频hd| 日韩欧美在线乱码| 久久国产乱子伦精品免费另类| 国产精品爽爽va在线观看网站| 精品人妻熟女av久视频| 小蜜桃在线观看免费完整版高清| 99国产极品粉嫩在线观看| 国产亚洲欧美在线一区二区| 久久99热6这里只有精品| av专区在线播放| 国产一级毛片七仙女欲春2| 搡女人真爽免费视频火全软件 | av在线观看视频网站免费| 无遮挡黄片免费观看| 亚洲精品影视一区二区三区av| 嫩草影院入口| 桃红色精品国产亚洲av| 欧美丝袜亚洲另类 | 亚洲国产欧洲综合997久久,| 亚洲片人在线观看| 亚洲精品色激情综合| 99国产综合亚洲精品| 男人的好看免费观看在线视频| 中文字幕免费在线视频6| 天堂av国产一区二区熟女人妻| 国产老妇女一区| 日韩欧美国产一区二区入口| 丰满的人妻完整版| 我要搜黄色片| 色精品久久人妻99蜜桃| 国产成人欧美在线观看| 亚洲三级黄色毛片| av天堂中文字幕网| 无遮挡黄片免费观看| 露出奶头的视频| 久久久色成人| 久久久久精品国产欧美久久久| 精品久久久久久久久久久久久| 看黄色毛片网站| 精品一区二区三区视频在线观看免费| 久久人人精品亚洲av| 精品午夜福利在线看| 69人妻影院| 女人被狂操c到高潮| 久久人妻av系列| 特级一级黄色大片| 成人毛片a级毛片在线播放| 亚洲激情在线av| 亚洲av第一区精品v没综合| 真人一进一出gif抽搐免费| 中文在线观看免费www的网站| 国产精品亚洲美女久久久| 白带黄色成豆腐渣| 国语自产精品视频在线第100页| 亚洲电影在线观看av| 国产爱豆传媒在线观看| 欧美日本亚洲视频在线播放| 亚州av有码| 午夜亚洲福利在线播放| 男人和女人高潮做爰伦理| 国产69精品久久久久777片| 最近最新中文字幕大全电影3| 久久久久免费精品人妻一区二区| av天堂中文字幕网| 日本五十路高清| 久久精品国产清高在天天线| 听说在线观看完整版免费高清| 禁无遮挡网站| 色综合欧美亚洲国产小说| 国产av麻豆久久久久久久| 国产精品98久久久久久宅男小说| 免费电影在线观看免费观看| 老司机午夜十八禁免费视频| 欧美日韩瑟瑟在线播放| 欧美高清成人免费视频www| 亚洲精品在线观看二区| 性色av乱码一区二区三区2| 精品99又大又爽又粗少妇毛片 | 深夜a级毛片| 日日干狠狠操夜夜爽| 亚洲美女搞黄在线观看 | 国产精品亚洲美女久久久| 美女免费视频网站| 午夜福利视频1000在线观看| 免费av观看视频| 少妇人妻一区二区三区视频| 欧美xxxx黑人xx丫x性爽| 国产高清三级在线| 亚洲 国产 在线| 成人特级av手机在线观看| 最近最新免费中文字幕在线| 精品人妻熟女av久视频| 久久久色成人| 91麻豆av在线| 亚洲三级黄色毛片| 日本熟妇午夜| 亚洲片人在线观看| 成人av在线播放网站| 国产精品不卡视频一区二区 | 99热这里只有精品一区| 麻豆久久精品国产亚洲av| 成人特级av手机在线观看| 老司机午夜十八禁免费视频| 亚洲av中文字字幕乱码综合| 国产免费男女视频| 99视频精品全部免费 在线| 国产精品久久久久久精品电影| 国产一区二区三区在线臀色熟女| 欧美xxxx黑人xx丫x性爽| 国产欧美日韩一区二区精品| 欧美成人免费av一区二区三区| 免费av观看视频| 少妇高潮的动态图| 人妻夜夜爽99麻豆av| 搞女人的毛片| 精品久久久久久成人av| 99久久精品一区二区三区| 成人国产综合亚洲| 老熟妇仑乱视频hdxx| 久久精品国产清高在天天线| 狠狠狠狠99中文字幕| 国产中年淑女户外野战色| 国产精品电影一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产三级在线视频| 丰满乱子伦码专区| 1024手机看黄色片| 免费黄网站久久成人精品 | 一a级毛片在线观看| 亚洲av一区综合| 亚洲人成电影免费在线| 国产探花在线观看一区二区| 在线观看免费视频日本深夜| 无遮挡黄片免费观看| 日本成人三级电影网站| 午夜精品一区二区三区免费看| 天堂影院成人在线观看| 好看av亚洲va欧美ⅴa在| 全区人妻精品视频| 国产高清视频在线播放一区| 国产av一区在线观看免费| 99热精品在线国产| 老司机深夜福利视频在线观看| 一区二区三区高清视频在线| 色综合亚洲欧美另类图片| 亚洲性夜色夜夜综合| 精品无人区乱码1区二区| 男女那种视频在线观看| av天堂中文字幕网| 久久精品久久久久久噜噜老黄 | 亚洲国产高清在线一区二区三| 国产欧美日韩一区二区三| 午夜福利在线观看吧| 亚洲久久久久久中文字幕| av天堂中文字幕网| 免费看美女性在线毛片视频| 啦啦啦观看免费观看视频高清| 国产三级在线视频| 欧美色视频一区免费| 一个人看视频在线观看www免费| 丰满人妻一区二区三区视频av| 中文字幕精品亚洲无线码一区| 黄色配什么色好看| 哪里可以看免费的av片| 国产av一区在线观看免费| 又粗又爽又猛毛片免费看| 欧美在线一区亚洲| 91在线观看av| 亚州av有码| 成人鲁丝片一二三区免费| 在线天堂最新版资源| 乱人视频在线观看| 精品国内亚洲2022精品成人| www.色视频.com| 18禁黄网站禁片午夜丰满| 国产精品国产高清国产av| 亚洲精品粉嫩美女一区| 在线播放国产精品三级| 日韩成人在线观看一区二区三区| av女优亚洲男人天堂| 又黄又爽又免费观看的视频| 国产精品自产拍在线观看55亚洲| 精品一区二区三区av网在线观看| 久久精品国产自在天天线| 久久久国产成人免费| 成人国产综合亚洲| 欧美性猛交╳xxx乱大交人| 在线播放国产精品三级| 97人妻精品一区二区三区麻豆| 一本久久中文字幕| 搡老熟女国产l中国老女人| 午夜a级毛片| 可以在线观看的亚洲视频| 亚洲色图av天堂| 亚洲国产欧美人成| 国产成人a区在线观看| 精品久久久久久,| 日韩亚洲欧美综合| av福利片在线观看| 国产精品精品国产色婷婷| 国产 一区 欧美 日韩| 精品久久久久久,| 久久久色成人| 亚洲五月天丁香| 少妇的逼水好多| 久久久久久久久大av| 亚洲人成伊人成综合网2020| 午夜两性在线视频| 欧美zozozo另类| 赤兔流量卡办理| 国产精品一及| 精品99又大又爽又粗少妇毛片 | 综合色av麻豆| 亚洲三级黄色毛片| 精品乱码久久久久久99久播| 欧洲精品卡2卡3卡4卡5卡区| 日韩成人在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 中文亚洲av片在线观看爽| 国产精品女同一区二区软件 | 国产一区二区三区视频了| 日日干狠狠操夜夜爽| 看片在线看免费视频| 男女做爰动态图高潮gif福利片| 国产中年淑女户外野战色| 亚洲 欧美 日韩 在线 免费| 天天躁日日操中文字幕| 国产精品99久久久久久久久| 国产欧美日韩精品亚洲av| 日韩欧美在线二视频| 亚洲真实伦在线观看| 欧美成狂野欧美在线观看| 成人高潮视频无遮挡免费网站| 麻豆国产97在线/欧美| 国产成人啪精品午夜网站| 此物有八面人人有两片| 久99久视频精品免费| 九九在线视频观看精品| 色精品久久人妻99蜜桃| 伦理电影大哥的女人| 听说在线观看完整版免费高清| 人妻丰满熟妇av一区二区三区| 国产主播在线观看一区二区| 色av中文字幕| 国产午夜精品久久久久久一区二区三区 | 国产精品野战在线观看| 精品午夜福利在线看| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区在线观看日韩| 欧美中文日本在线观看视频| 国产野战对白在线观看| 国内精品久久久久久久电影| 好男人电影高清在线观看| 美女cb高潮喷水在线观看| 亚洲av五月六月丁香网| 两个人视频免费观看高清| 日韩欧美在线乱码| 淫秽高清视频在线观看| 小蜜桃在线观看免费完整版高清| 色哟哟·www| 国产在视频线在精品| 91av网一区二区| 美女 人体艺术 gogo| 夜夜看夜夜爽夜夜摸| 女人十人毛片免费观看3o分钟| 中文字幕av在线有码专区| 欧美精品国产亚洲| 亚洲欧美日韩高清在线视频| 精品乱码久久久久久99久播| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 最新在线观看一区二区三区| 免费av不卡在线播放| 我的女老师完整版在线观看| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 国产午夜福利久久久久久| a级毛片a级免费在线| 极品教师在线视频| 亚洲 欧美 日韩 在线 免费| or卡值多少钱| 欧美高清性xxxxhd video| 久久精品91蜜桃| 亚洲av日韩精品久久久久久密| 久久草成人影院| 国产乱人伦免费视频| АⅤ资源中文在线天堂| 免费在线观看亚洲国产| 99视频精品全部免费 在线| a级毛片免费高清观看在线播放| 91麻豆av在线| 亚洲一区二区三区色噜噜| 国产黄色小视频在线观看| 午夜日韩欧美国产| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添小说| 亚洲自偷自拍三级| 高潮久久久久久久久久久不卡| 欧美午夜高清在线| 国产精品久久电影中文字幕| 婷婷六月久久综合丁香| 黄色丝袜av网址大全| 日本 av在线| 特大巨黑吊av在线直播| 丁香欧美五月| 成人亚洲精品av一区二区| 少妇人妻一区二区三区视频| 91麻豆av在线| 最新中文字幕久久久久| 国产v大片淫在线免费观看| 2021天堂中文幕一二区在线观| 在线免费观看不下载黄p国产 | 免费电影在线观看免费观看| 一本精品99久久精品77| 丰满人妻一区二区三区视频av| 国产成人啪精品午夜网站| 久久精品国产自在天天线| 蜜桃亚洲精品一区二区三区| 成人三级黄色视频| 成熟少妇高潮喷水视频| 草草在线视频免费看| 蜜桃久久精品国产亚洲av| 少妇丰满av| 欧美成人免费av一区二区三区| 少妇人妻一区二区三区视频| 精品日产1卡2卡| 97人妻精品一区二区三区麻豆| 性欧美人与动物交配| 村上凉子中文字幕在线| 91麻豆av在线| 搡女人真爽免费视频火全软件 | 日韩av在线大香蕉| 免费无遮挡裸体视频| 久久久久久久精品吃奶| 亚洲专区国产一区二区| 嫁个100分男人电影在线观看| 黄色日韩在线| 桃红色精品国产亚洲av| 我的老师免费观看完整版| 国产亚洲av嫩草精品影院| 欧美日韩综合久久久久久 | 久久久久国内视频| 日本黄色视频三级网站网址| 精品欧美国产一区二区三| 精品一区二区三区视频在线观看免费| 久久久国产成人免费| 12—13女人毛片做爰片一| 成人av在线播放网站| 国内少妇人妻偷人精品xxx网站| 久9热在线精品视频| 成人毛片a级毛片在线播放| 日韩有码中文字幕| 人人妻,人人澡人人爽秒播| 免费在线观看成人毛片| 91av网一区二区| 精品日产1卡2卡| avwww免费| 久久精品国产亚洲av天美| 九九久久精品国产亚洲av麻豆| 国产av一区在线观看免费| 欧美午夜高清在线| 中文字幕久久专区| 90打野战视频偷拍视频| 欧美绝顶高潮抽搐喷水| 99久久无色码亚洲精品果冻| 在线a可以看的网站| 久久久久九九精品影院| 悠悠久久av| 一进一出抽搐动态| 午夜福利18| 亚洲欧美精品综合久久99| 久久精品国产自在天天线| 日韩欧美在线二视频| 他把我摸到了高潮在线观看| 天堂网av新在线| 美女高潮喷水抽搐中文字幕| 亚洲av中文字字幕乱码综合| 国产精品女同一区二区软件 | 少妇被粗大猛烈的视频| 精品人妻一区二区三区麻豆 | 午夜福利在线观看免费完整高清在 | 国产成+人综合+亚洲专区| 男女之事视频高清在线观看| 国产精品1区2区在线观看.| 99riav亚洲国产免费| 欧美日韩黄片免| 可以在线观看的亚洲视频| 久久久久亚洲av毛片大全| 色视频www国产| 婷婷色综合大香蕉| www.熟女人妻精品国产| 哪里可以看免费的av片| 国产熟女xx| 日本a在线网址| 国产黄片美女视频| 青草久久国产| 亚洲av成人不卡在线观看播放网| 久久草成人影院| 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 亚洲国产欧美人成| 久久精品国产亚洲av涩爱 | 国产极品精品免费视频能看的| 最近在线观看免费完整版| 成人国产一区最新在线观看| 亚洲一区高清亚洲精品| 搡老岳熟女国产| 国产久久久一区二区三区| 宅男免费午夜| 成人精品一区二区免费| 狠狠狠狠99中文字幕| 在线观看66精品国产| 亚洲av免费高清在线观看| 大型黄色视频在线免费观看| 国产精品伦人一区二区| 一本精品99久久精品77| 高清在线国产一区| 精品一区二区三区视频在线观看免费| 69人妻影院| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 最近最新中文字幕大全电影3| 国产高清激情床上av| 中出人妻视频一区二区| 久久欧美精品欧美久久欧美| 国产真实乱freesex| 身体一侧抽搐| 俺也久久电影网| 免费av不卡在线播放| 欧美成人性av电影在线观看| 色5月婷婷丁香| 亚洲av日韩精品久久久久久密| 日韩大尺度精品在线看网址| 免费电影在线观看免费观看| 国内久久婷婷六月综合欲色啪| 男女做爰动态图高潮gif福利片| 国产成人a区在线观看| 午夜福利免费观看在线| 亚洲av五月六月丁香网| 亚洲av熟女| 在线十欧美十亚洲十日本专区| 国产91精品成人一区二区三区| 国产三级黄色录像| 直男gayav资源| 深夜a级毛片| 97超级碰碰碰精品色视频在线观看| 久久人人精品亚洲av| 久久伊人香网站| 国产av不卡久久| 在线观看舔阴道视频| 男女视频在线观看网站免费| 国产精品综合久久久久久久免费| 欧美黄色淫秽网站| 国产精品,欧美在线| 午夜视频国产福利| 亚洲欧美日韩卡通动漫| av在线蜜桃| 欧美丝袜亚洲另类 | 国产av不卡久久| 成人永久免费在线观看视频| 俺也久久电影网| 此物有八面人人有两片| 亚洲精品乱码久久久v下载方式| 欧美精品国产亚洲| 一个人看视频在线观看www免费| 国产成人影院久久av| 精品人妻偷拍中文字幕| 中文字幕免费在线视频6| www.999成人在线观看| 9191精品国产免费久久| 毛片一级片免费看久久久久 | 亚洲国产精品999在线| 国产综合懂色| 欧美性猛交╳xxx乱大交人| .国产精品久久| 欧美日韩福利视频一区二区| 精品久久久久久久久亚洲 | 国产免费一级a男人的天堂| 两人在一起打扑克的视频| 亚洲美女视频黄频| a级一级毛片免费在线观看| 最近中文字幕高清免费大全6 | 日韩精品青青久久久久久| 免费黄网站久久成人精品 | 亚洲最大成人中文| 午夜精品在线福利| 亚洲成人中文字幕在线播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产69精品久久久久777片| 国产精品久久久久久精品电影| 赤兔流量卡办理| 丝袜美腿在线中文| 免费在线观看影片大全网站| 两人在一起打扑克的视频| 亚洲激情在线av| netflix在线观看网站| 午夜福利在线观看吧| 不卡一级毛片| 久久久国产成人免费| 三级男女做爰猛烈吃奶摸视频| 亚洲精品在线美女| 青草久久国产| 一级黄色大片毛片| 国产精华一区二区三区| 中文亚洲av片在线观看爽| 欧美成人a在线观看| 亚洲成人免费电影在线观看| 国产精品嫩草影院av在线观看 | www.熟女人妻精品国产| 香蕉av资源在线| 日韩欧美在线二视频| 97超视频在线观看视频| 日本一本二区三区精品| 国产欧美日韩精品亚洲av| 亚洲av中文字字幕乱码综合| 舔av片在线| 韩国av一区二区三区四区| 性色av乱码一区二区三区2| 激情在线观看视频在线高清| 99riav亚洲国产免费| 色哟哟哟哟哟哟| 51国产日韩欧美| 国产真实伦视频高清在线观看 | 怎么达到女性高潮| 女生性感内裤真人,穿戴方法视频| 欧美日本视频| 欧美成人性av电影在线观看| 日本成人三级电影网站| 毛片女人毛片| 中文字幕av成人在线电影| 最新在线观看一区二区三区| 欧美bdsm另类| 国产 一区 欧美 日韩| 国产精品三级大全| 亚洲无线观看免费| 夜夜爽天天搞| 久久精品91蜜桃| 18禁黄网站禁片午夜丰满| 日本五十路高清| 精品午夜福利视频在线观看一区| 久久婷婷人人爽人人干人人爱| 真实男女啪啪啪动态图| 欧美黑人欧美精品刺激| 免费一级毛片在线播放高清视频| 一个人免费在线观看电影| 看十八女毛片水多多多| 给我免费播放毛片高清在线观看| 国产爱豆传媒在线观看| 又爽又黄无遮挡网站| 国产欧美日韩一区二区三| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 麻豆国产97在线/欧美| av黄色大香蕉| 久久久成人免费电影| 国产高清三级在线| 亚洲精品色激情综合| 免费观看精品视频网站| АⅤ资源中文在线天堂| 啦啦啦韩国在线观看视频| 天堂√8在线中文| 国产在线男女| 一级毛片久久久久久久久女| 黄色配什么色好看| 久久久久国内视频| 日本黄大片高清| 亚洲人成伊人成综合网2020| 村上凉子中文字幕在线| 人妻久久中文字幕网| 午夜福利欧美成人| 国产大屁股一区二区在线视频| 床上黄色一级片| 国产熟女xx| 欧美黑人巨大hd| 亚洲第一区二区三区不卡| 可以在线观看的亚洲视频| 国产精品久久久久久亚洲av鲁大| 国内精品久久久久久久电影| 97人妻精品一区二区三区麻豆| 又爽又黄a免费视频| 成年女人永久免费观看视频| 内射极品少妇av片p| 日韩欧美精品免费久久 | 亚洲最大成人手机在线| 97碰自拍视频|