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

    振動(dòng)式水稻精密播種裝置機(jī)理分析與試驗(yàn)

    2018-07-05 11:38:14鹿芳媛邢緒坡李宏偉郭林杰
    關(guān)鍵詞:供種雜交稻定量

    鹿芳媛 馬 旭 齊 龍 邢緒坡 李宏偉 郭林杰

    (華南農(nóng)業(yè)大學(xué)工程學(xué)院, 廣州 510642)

    0 引言

    水稻是我國主要糧食作物,其種植機(jī)械化是我國水稻生產(chǎn)機(jī)械化的薄弱環(huán)節(jié),也是我國水稻生產(chǎn)全程機(jī)械化發(fā)展中的瓶頸[1],提高水稻秧苗育秧播種機(jī)械化水平是實(shí)現(xiàn)我國水稻種植機(jī)械化的關(guān)鍵環(huán)節(jié)之一。

    目前,水稻機(jī)械化種植方式主要包括機(jī)栽植和機(jī)直播。美國、意大利、澳大利亞等發(fā)達(dá)國家采用機(jī)械化直播方式種植水稻,對秧盤育秧播種機(jī)械的研制較少[2-3];亞洲的水稻秧盤育秧設(shè)備較多,其中日本和韓國的技術(shù)水平最高,如日本久保田、洋馬、井關(guān)等株式會(huì)社生產(chǎn)的育秧播種流水線,采用機(jī)械式排種器,設(shè)備自動(dòng)化程度高、工藝精湛,但主要適用于常規(guī)稻4~8 粒/格的撒播或3~6 粒/穴的穴播作業(yè),不適用雜交稻1~3粒/穴的精密播種需要[4]。

    我國水稻機(jī)械化種植主要采用育秧移栽。水稻精密播種器的研發(fā)始于20世紀(jì)70年代,早期研究以機(jī)械式排種器為主,80年代初開始研究氣力式排種器[5],隨著水稻秧盤育秧流水線的研制,90年代起振動(dòng)式播種流水線發(fā)展起來[6]。可以看出,我國對于水稻精密排種器的研究,按其播種方式分為機(jī)械式、振動(dòng)式和氣力式。傳統(tǒng)機(jī)械式播種器結(jié)構(gòu)簡單、適應(yīng)性強(qiáng),主要適用于水稻大播量條播和撒播,較難滿足雜交稻精密播種需求[3]。氣力式水稻精密播種器作業(yè)效率高,對種子損傷小,能夠?qū)崿F(xiàn)單粒精量點(diǎn)播[7],但由于水稻主要采用芽種播種育秧,帶芽播種時(shí)種芽的脫落易造成吸孔堵塞而影響播種性能,且吸孔堵塞問題較難解決。

    振動(dòng)式水稻精密播種是一種有效的精密播種方式,可分為機(jī)械振動(dòng)式、電磁振動(dòng)式和氣力振動(dòng)式。張學(xué)義、楊堅(jiān)、李志偉等[6,8-9]改進(jìn)電磁振動(dòng)式播種器用于水稻育秧播種,可實(shí)現(xiàn)3~5 粒/穴的播種要求。張斌等[10]提出利用導(dǎo)向板使稻種在振動(dòng)式排種器上沿胚胎方向有序排列,實(shí)現(xiàn)了超級(jí)雜交稻定向播種。齊龍等[11]采用交叉導(dǎo)流式振動(dòng)種箱配合螺旋勺式槽輪定量供種,與氣動(dòng)振動(dòng)勻種機(jī)構(gòu)相結(jié)合,研制了一種振動(dòng)式水稻精密播種裝置,能夠滿足雜交稻和常規(guī)稻種子播種育秧要求,但該裝置的結(jié)構(gòu)與工作參數(shù)還需優(yōu)化,以保證穩(wěn)定的工作性能。

    為此,本文將解析交叉導(dǎo)流式振動(dòng)種箱的工作機(jī)理,利用離散元法[12-15]模擬分析與試驗(yàn)探究該機(jī)構(gòu)的供種性能,分析振動(dòng)勻種機(jī)構(gòu)的不同種槽板對水稻種子導(dǎo)向作用和流動(dòng)速度(以下簡稱流速)的影響,探究不同類型種槽板對播種性能的影響,并通過播種性能試驗(yàn)優(yōu)化種槽板結(jié)構(gòu)參數(shù)。

    1 振動(dòng)式精密播種裝置工作原理

    振動(dòng)式水稻精密播種裝置是水稻秧盤育秧精密播種流水線的關(guān)鍵技術(shù)部件,其播種質(zhì)量直接影響育秧流水線的作業(yè)性能和效率。該播種裝置主要由兩部分組成:定量供種機(jī)構(gòu)和振動(dòng)勻種機(jī)構(gòu),其結(jié)構(gòu)示意圖如圖1所示。

    圖1 振動(dòng)式水稻精密播種裝置結(jié)構(gòu)示意圖Fig.1 Sketch map of vibrating precision sowing devise for rice1.種箱 2.種箱振動(dòng)板 3.電磁振動(dòng)器 4.調(diào)節(jié)門 5.清種毛刷 6.螺旋勺式排種輪 7.支架 8.氣動(dòng)振動(dòng)器 9.振盤振動(dòng)臂 10.振動(dòng)器安裝板 11.種槽板 12.T型板 13.調(diào)速電動(dòng)機(jī) 14.第二導(dǎo)種板 15.第一導(dǎo)種板

    定量供種機(jī)構(gòu)主要包括種箱、第一導(dǎo)種板、種箱振動(dòng)板、第二導(dǎo)種板、電磁振動(dòng)器、調(diào)節(jié)門、清種毛刷、螺旋勺式排種輪(以下簡稱排種輪)和調(diào)速電動(dòng)機(jī);振動(dòng)勻種機(jī)構(gòu)包括氣動(dòng)振動(dòng)器、振動(dòng)器安裝板、振盤振動(dòng)臂、振動(dòng)種盤、支架等部件。其工作原理為:定量供種機(jī)構(gòu)采用交叉導(dǎo)流式振動(dòng)種箱(由圖1中零件1、2、3、14、15構(gòu)成),使水稻種子在排種輪上方的充填區(qū)形成“Z”型連續(xù)種流,然后隨排種輪轉(zhuǎn)動(dòng)定量地排出種箱,落入振動(dòng)勻種機(jī)構(gòu)的振動(dòng)種盤中。在振動(dòng)勻種機(jī)構(gòu)的振動(dòng)作用下,水稻種子在振動(dòng)種盤上部T型板中進(jìn)行篩分勻種,在振動(dòng)種盤下部種槽板的導(dǎo)向作用下有序排隊(duì),形成連續(xù)、均勻的種子流,最終精量地落入育秧流水線輸送帶上的秧盤中,完成精少量播種。

    定量供種機(jī)構(gòu)連續(xù)定量的供種性能和振動(dòng)勻種機(jī)構(gòu)均勻穩(wěn)定的播種性能,是實(shí)現(xiàn)水稻精密播種的重要條件,二者缺一不可。因此,分別探究定量供種機(jī)構(gòu)的供種機(jī)理與振動(dòng)勻種機(jī)構(gòu)的播種機(jī)理,分析影響裝置播種性能的因素,使播種合格指數(shù)提高到92%以上,空穴指數(shù)控制在2%以下。

    1.1 定量供種機(jī)構(gòu)的供種機(jī)理

    在播種過程中隨著排種輪的轉(zhuǎn)動(dòng)種量不斷減少,常規(guī)種箱內(nèi)靠近排種輪出流口的種子流速較其他位置快,中部種子流速較種箱內(nèi)壁快,排種輪中部與兩側(cè)充填特性不一致,導(dǎo)致供種量不穩(wěn)定。為此,李志偉、馬旭、周海波等[6,16-17]提出了多層導(dǎo)種板結(jié)構(gòu)的種箱,提高了種箱的供種性能,實(shí)現(xiàn)了均勻定量供種。

    本文研究的定量供種機(jī)構(gòu)交叉導(dǎo)流式振動(dòng)種箱,采用第一導(dǎo)種板、種箱振動(dòng)板和第二導(dǎo)種板交疊安裝形成的交叉導(dǎo)流結(jié)構(gòu),利用多層導(dǎo)種板分擔(dān)了種箱內(nèi)種子的大部分重力,避免了種子直接堆積在排種輪上,而在排種輪上方充填區(qū)形成“Z”型連續(xù)種流,有效提高了充種性能,且減少了多余種子在清種毛刷部位旋轉(zhuǎn)翻滾產(chǎn)生的渦流與碎種現(xiàn)象。此外,由于水稻種子自身的物理特性,外殼存在芒刺,種子間的摩擦因數(shù)較大,無序的分布狀態(tài)使其在種箱內(nèi)的流動(dòng)性差。同時(shí),催芽后的水稻種子抵抗損傷的能力很弱,種子間較大的摩擦力及排種輪的強(qiáng)制帶動(dòng)作用會(huì)加劇傷種。因此,通過電磁振動(dòng)器連接振動(dòng)板,在種箱內(nèi)側(cè)產(chǎn)生周期性的小幅振動(dòng),加速種箱內(nèi)側(cè)種子的流動(dòng)性,并使排種輪充填區(qū)的種子分布狀態(tài)較為“蓬松”,減少種子的機(jī)械損傷,優(yōu)化供種性能。

    1.2 振動(dòng)勻種機(jī)構(gòu)的勻種機(jī)理

    振動(dòng)種盤的振動(dòng)性能直接決定了水稻種子的運(yùn)動(dòng)特性以及機(jī)構(gòu)的播種性能。為分析振動(dòng)種盤的運(yùn)動(dòng)特性,將振動(dòng)勻種機(jī)構(gòu)分為3部分:支架、氣動(dòng)振動(dòng)器及受迫振動(dòng)系統(tǒng)(包括振盤振動(dòng)臂和振動(dòng)種盤)。氣動(dòng)振動(dòng)器對受迫振動(dòng)系統(tǒng)產(chǎn)生激振力,在激振作用下振動(dòng)種盤繞振盤振動(dòng)臂與機(jī)架的固定點(diǎn)產(chǎn)生周期性的往復(fù)擺動(dòng)。為探明該往復(fù)運(yùn)動(dòng)特性,對受迫振動(dòng)系統(tǒng)進(jìn)行動(dòng)力學(xué)分析。

    1.2.1受迫振動(dòng)系統(tǒng)動(dòng)力學(xué)分析

    通過振動(dòng)特性測試分析可知,振動(dòng)種盤產(chǎn)生的擺角始終小于0.1°,對于該系統(tǒng)來說振動(dòng)種盤在鉛垂方向上產(chǎn)生的位移可忽略不計(jì),位移主要產(chǎn)生于水平方向。因此,可將振盤振動(dòng)臂的作用力在鉛垂方向上等效為提拉振動(dòng)種盤的拉力,在水平方向上等效為彈性力與阻尼力,則振動(dòng)種盤在水平方向上做單自由度簡諧受迫振動(dòng)[18],其力學(xué)模型如圖2所示。為便于分析,將受迫振動(dòng)系統(tǒng)簡化為質(zhì)量為M的等效物塊。

    圖2 單自由度受迫振動(dòng)系統(tǒng)力學(xué)模型Fig.2 Mechanical model of single degree of freedom forced vibration system

    氣動(dòng)振動(dòng)器產(chǎn)生的激振力運(yùn)動(dòng)方程為

    F(t)=Hsin(ωt)

    (1)

    式中H——激振力的振幅,mm

    ω——激振力的角頻率,rad/s

    以平衡位置O點(diǎn)為坐標(biāo)原點(diǎn),Y軸鉛垂向上為正,則等效物塊沿鉛垂方向的運(yùn)動(dòng)微分方程為

    (2)

    其中

    δ=α+γ

    式中Fl——振動(dòng)系統(tǒng)受到的支架拉力,N

    α——振動(dòng)器安裝板夾角,(°)

    γ——槽板與水平面夾角,(°)

    Fl(t)=Mg-Hsin(ωt)cosδ

    (3)

    表明振動(dòng)勻種機(jī)構(gòu)在工作過程中,振盤支架以及振盤振動(dòng)臂所受到的拉力隨激振力的變化而發(fā)生周期性改變。

    對等效物塊水平方向進(jìn)行受力分析,X水平向左為正,則等效物塊的運(yùn)動(dòng)微分方程為

    (4)

    式中c——等效阻尼系數(shù)

    k——等效彈簧系數(shù)

    其解[19]為

    X(t)=X1(t)+X2(t)=
    Ae-ntsin(ωdt+Ψ)+Bsin(ωt-φ)sinδ

    (5)

    其中

    式中n——系統(tǒng)的衰減系數(shù)

    ωd——衰減振動(dòng)的角頻率,rad/s

    Ψ——自由振動(dòng)的相位差,(°)

    ωn——系統(tǒng)的固有角頻率,rad/s

    A——初始振幅,mm

    B——受迫振動(dòng)的振幅,mm

    φ——受迫振動(dòng)的相位差,(°)

    可以看出振動(dòng)種盤的運(yùn)動(dòng)由兩部分組成:角頻率為ωd的衰減振動(dòng)和角頻率為ω的受迫振動(dòng)。在受迫振動(dòng)系統(tǒng)的衰減振動(dòng)中,由于阻尼的存在,衰減振動(dòng)會(huì)逐漸消失,在衰減振動(dòng)完全消失之前,系統(tǒng)的振動(dòng)過程稱為暫態(tài)響應(yīng)。衰減振動(dòng)消失后,系統(tǒng)振動(dòng)的過程為穩(wěn)定的等幅受迫振動(dòng),該過程稱為振動(dòng)系統(tǒng)的穩(wěn)態(tài)響應(yīng)。當(dāng)受迫振動(dòng)系統(tǒng)達(dá)到穩(wěn)態(tài)響應(yīng)階段,其運(yùn)動(dòng)方程為

    X(t)=X2(t)=Bsin(ωt-φ)sinδ

    可以看出穩(wěn)態(tài)過程即為簡諧振動(dòng),其頻率與激振力頻率相等,相位角與激振力的相位角相差φ。通過上述分析,在使用振動(dòng)勻種機(jī)構(gòu)時(shí),需開機(jī)運(yùn)行等待衰減振動(dòng)的振幅逐漸減小,在受迫振動(dòng)系統(tǒng)振動(dòng)狀態(tài)穩(wěn)定后再進(jìn)行調(diào)試播種。

    1.2.2種子在振動(dòng)種盤上的運(yùn)動(dòng)理論分析

    振動(dòng)勻種機(jī)構(gòu)工作時(shí),水稻種子大部分呈松散均勻的狀態(tài)分布在振動(dòng)種盤中。以單粒水稻種子為研究對象,探究種子在振動(dòng)種盤上部T型板和下部槽板上的運(yùn)動(dòng)特性[11]。分別對種子在T型板和種槽板上的受力進(jìn)行分析(振動(dòng)種盤結(jié)構(gòu)見圖1),如圖3所示。

    圖3 種子在振動(dòng)種盤上的受力圖Fig.3 Force diagrams of seed on vibration plate

    基于上述單自由度振動(dòng)系統(tǒng)分析,設(shè)激振力產(chǎn)生的振動(dòng)幅值為h(mm),則振動(dòng)種盤的位移為

    s=hsin(ωt-φ)sinδ

    由圖3a分析得到種子在T型板上的下滑條件為

    (6)

    式中μ——種子在振動(dòng)種盤上的滑動(dòng)摩擦角,(°)

    β——T型板與水平面夾角,(°)

    種子在T型板上的跳動(dòng)條件為

    (7)

    由圖3b得到種子在槽板上的下滑條件為

    (8)

    式中θ——種槽板夾角,(°)

    種子在槽板上的跳動(dòng)條件為

    (9)

    其中,當(dāng)0γ,則滿足種子在槽板上的下滑或跳動(dòng)條件時(shí),必定滿足種子在T型板上的下滑或跳動(dòng)條件。此外,由式(8)可知,種子在種槽板上的下滑行走條件與種槽板夾角θ有關(guān),下面將通過離散元法探究種槽板夾角θ對種子運(yùn)動(dòng)狀態(tài)的影響。

    2 振動(dòng)式精密播種裝置仿真分析

    利用離散元仿真軟件EDEM,分別對振動(dòng)式精密播種裝置的定量供種機(jī)構(gòu)和振動(dòng)勻種機(jī)構(gòu)進(jìn)行播種性能仿真。

    2.1 定量供種機(jī)構(gòu)的供種過程仿真

    為探究供種時(shí)交叉導(dǎo)流式種箱內(nèi)水稻種子的流動(dòng)特性,以及振動(dòng)對供種性能的影響,利用離散元仿真軟件EDEM模擬定量供種機(jī)構(gòu)的工作過程,以種箱振動(dòng)板有、無振動(dòng)作用進(jìn)行仿真對比試驗(yàn)。通過三維建模軟件Solidworks創(chuàng)建定量供種機(jī)構(gòu)的三維實(shí)體模型并導(dǎo)入EDEM軟件中,創(chuàng)建水稻芽種顆粒,由23個(gè)基礎(chǔ)顆粒球組成[20]。設(shè)置水稻泊松比為0.25,剪切模量為108 MPa,密度為1 098 kg/m3;不銹鋼板(種箱及導(dǎo)種板材料)泊松比為0.30,剪切模量為79 000 MPa,密度為7 850 kg/m3;尼龍1010(螺旋勺式排種輪材料)泊松比為0.28,剪切模量為1 070 MPa,密度為1 050 kg/m3。水稻種子含水率約37.3%,種子顆粒與各材料間的接觸參數(shù)見表1。根據(jù)定量供種機(jī)構(gòu)實(shí)際工作狀態(tài)設(shè)置模型運(yùn)動(dòng)參數(shù),排種輪模型運(yùn)動(dòng)設(shè)置為線性旋轉(zhuǎn),轉(zhuǎn)速為7.0 r/min。使用的電磁振動(dòng)器型號(hào)為DH180-2-1,振動(dòng)頻率為50 Hz,通過振動(dòng)測試分析得到種箱振動(dòng)板在電磁振動(dòng)器的作用下產(chǎn)生簡諧振動(dòng),且振幅可調(diào)范圍為0~0.02 mm,因此設(shè)置振動(dòng)板運(yùn)動(dòng)為正弦平移,頻率為50 Hz,根據(jù)播量要求振幅選用0.015 mm。

    表1 仿真所需接觸參數(shù)Tab.1 Contact parameters of simulation

    2.1.1種子在種箱內(nèi)的流動(dòng)特性分析

    以振動(dòng)板的振動(dòng)作用作為試驗(yàn)因素,設(shè)計(jì)兩組供種過程模擬試驗(yàn),一組設(shè)置振動(dòng),另一組不設(shè)置振動(dòng),其余參數(shù)設(shè)置均相同,仿真時(shí)間為85 s。仿真完成后,對種箱內(nèi)的種子顆粒進(jìn)行染色分層,由下至上將種子分為10層,以最底層為第1種層,以最頂層為第10種層。分別選取兩組仿真試驗(yàn)過程中的4個(gè)時(shí)刻進(jìn)行分析,如圖4所示,左圖為種箱振動(dòng)板有振動(dòng)作用,右圖為無振動(dòng)。

    圖4 交叉導(dǎo)流式振動(dòng)種箱供種過程仿真Fig.4 Feeding process simulation of cross-flow vibration seeds box

    圖4a為初始時(shí)刻水稻種子自然填充于種箱內(nèi)第一導(dǎo)種板下部、種箱振動(dòng)板下部及第二導(dǎo)種板之間,形成了“Z”型連續(xù)種流,緩解了大量種子直接堆積在排種輪上方產(chǎn)生的正壓力。圖4b為排種輪轉(zhuǎn)動(dòng)了3 s時(shí),在此過程中,左圖供種機(jī)構(gòu)在振動(dòng)作用下第1~4種層種子按種層順序排出種箱,第5~10種層種子整體流動(dòng)平穩(wěn)下落;右圖供種機(jī)構(gòu)無振動(dòng)時(shí)第1~4種層種子隨排種輪的轉(zhuǎn)動(dòng)相互混合,亂序排出種箱,第5~7種層種子開始呈現(xiàn)漏斗狀流動(dòng)[21],第8~10種層平穩(wěn)下落。圖4c、4d分別為排種輪轉(zhuǎn)動(dòng)了10 s和22 s時(shí),在此過程中,左圖種箱右側(cè)由于振動(dòng)作用種子流速較快,故種子整體高度向右下方傾斜,呈現(xiàn)擴(kuò)散流動(dòng),流速較平穩(wěn);右圖種箱內(nèi)漏斗狀明顯,各層種子均形成了漏斗流,流速不穩(wěn)定。

    觀察兩組仿真試驗(yàn)整個(gè)供種過程,與導(dǎo)種板內(nèi)壁接觸的種子均形成了不流動(dòng)區(qū)。當(dāng)右側(cè)振動(dòng)板無振動(dòng)時(shí),種箱左右兩側(cè)皆存在不流動(dòng)區(qū),因此種子形成漏斗流動(dòng),速率不穩(wěn)定,種子密度不均勻,影響排種輪的充填性能。而當(dāng)右側(cè)振動(dòng)板振動(dòng)時(shí),振動(dòng)作用提高了右側(cè)種子流速,僅左側(cè)導(dǎo)種板上存在不流動(dòng)區(qū),緩解了漏斗流,促進(jìn)種子擴(kuò)散流動(dòng),穩(wěn)定了種子整體的流速和密度,提高了供種穩(wěn)定性。

    根據(jù)兩組試驗(yàn)仿真完成時(shí)間,有振動(dòng)時(shí)供種試驗(yàn)種箱內(nèi)的種子在78 s內(nèi)被排空,無振動(dòng)時(shí)則在82 s內(nèi)被排空,表明在振動(dòng)作用下排種輪的供種頻率(單位時(shí)間內(nèi)排種輪的供種質(zhì)量,g/s)較高。為確定振動(dòng)對供種頻率的提高比率,對兩組仿真試驗(yàn)的供種質(zhì)量進(jìn)行統(tǒng)計(jì),繪制定量供種機(jī)構(gòu)的供種質(zhì)量隨時(shí)間變化曲線對比,如圖5所示。對圖5中2條曲線做線性擬合,得到種箱振動(dòng)板有振動(dòng)時(shí)排種輪的供種質(zhì)量

    Q1=10.177t

    (10)

    圖5 定量供種機(jī)構(gòu)有無振動(dòng)時(shí)供種質(zhì)量對比Fig.5 Feeding weight comparison of quantitative-seeds feeding device with or without vibration

    Q2=9.851t

    (11)

    2.1.2種子在充填區(qū)分布狀態(tài)與流速分析

    為研究振動(dòng)對排種輪充填區(qū)內(nèi)的種子分布“蓬松”程度和種子流速的影響,在排種輪的強(qiáng)制層和帶動(dòng)層添加網(wǎng)格進(jìn)行分析,劃分強(qiáng)制層、帶動(dòng)層Ⅰ、帶動(dòng)層Ⅱ、帶動(dòng)層Ⅲ,如圖6所示。

    隨機(jī)選擇仿真過程的4個(gè)時(shí)刻,分別統(tǒng)計(jì)強(qiáng)制層、帶動(dòng)層Ⅰ、帶動(dòng)層Ⅱ、帶動(dòng)層Ⅲ網(wǎng)格內(nèi)的種子數(shù)量,統(tǒng)計(jì)結(jié)果如表2所示。通過表2中各網(wǎng)格內(nèi)的顆粒數(shù)量對比可以看出,振動(dòng)板有振動(dòng)比無振動(dòng)時(shí)強(qiáng)制層內(nèi)的種子數(shù)量多,而帶動(dòng)層(Ⅰ、Ⅱ、Ⅲ)內(nèi)的種子數(shù)量少,表明種箱內(nèi)側(cè)的振動(dòng)作用促進(jìn)排種輪充填區(qū)的種子分布狀態(tài)更為“蓬松”;振動(dòng)使強(qiáng)制層的種子獲得較快的流速,單位時(shí)間內(nèi)充填進(jìn)入排種輪凹槽內(nèi)的種子較無振動(dòng)時(shí)多,因而排種輪的供種頻率較高。此外,振動(dòng)時(shí)的顆粒數(shù)量變異系數(shù)小于無振動(dòng),表明振動(dòng)作用下排種輪的充種量更穩(wěn)定。

    圖6 排種輪強(qiáng)制層與帶動(dòng)層網(wǎng)格劃分Fig.6 Meshing for compulsion layers and driven layer of wheal

    表2 各網(wǎng)格區(qū)域顆粒數(shù)量和變異系數(shù)Tab.2 Particle number of grid regions

    圖7 不同區(qū)域顆粒流速對比曲線Fig.7 Particles velocity comparison curves of different regions

    對帶動(dòng)層和強(qiáng)制層的顆粒速度進(jìn)行統(tǒng)計(jì),繪制仿真3~21 s過程中種子的流速曲線,如圖7所示,進(jìn)一步分析振動(dòng)對充種區(qū)種子流速的影響。通過圖7a種子流速變化曲線對比,振動(dòng)作用對排種輪強(qiáng)制層的種子流速大小影響不明顯,而計(jì)算和比較兩組種子流速平均值時(shí)則有振動(dòng)略微大于無振動(dòng)時(shí)。由圖7b可知,有振動(dòng)時(shí)種子流速曲線整體高于無振動(dòng)時(shí),表明振動(dòng)作用加速了帶動(dòng)層Ⅰ區(qū)域的種子流速,從而促進(jìn)種子充填進(jìn)入排種輪的凹槽內(nèi),提高排種輪的供種頻率。通過圖7c可以看出,對于帶動(dòng)層Ⅱ的種子,在無振動(dòng)時(shí)速度變化劇烈,而在有振動(dòng)時(shí)速度穩(wěn)定,表明振動(dòng)作用有效地穩(wěn)定了該層種子流速。圖7d為帶動(dòng)層Ⅲ種子流速,從兩條速度對比曲線來看,種子流動(dòng)速度變化差異不明顯,通過趨勢線對比表明有振動(dòng)時(shí)種子流速變化稍平穩(wěn)。

    上述分析以定量供種機(jī)構(gòu)大播種量供種(適用于常規(guī)稻,通常播量標(biāo)準(zhǔn)秧盤80~120 g/盤,即3~7粒/穴或格)為例進(jìn)行探究。針對低播量供種(適用于雜交稻,通常播量標(biāo)準(zhǔn)秧盤40~60 g/盤,即1~4粒/穴)以相同的方法進(jìn)行模擬分析,得到種子在種箱內(nèi)交叉導(dǎo)流式結(jié)構(gòu)上的流動(dòng)特性以及在排種輪充填區(qū)的分布狀態(tài)和流速,在有無振動(dòng)作用對比仿真試驗(yàn)中表現(xiàn)出的規(guī)律與常規(guī)稻供種時(shí)基本一致,雜交稻低播量時(shí)有振動(dòng)比無振動(dòng)時(shí)排種輪的供種頻率高5.41%。

    綜上采用離散元法對定量供種機(jī)構(gòu)在有振動(dòng)與無振動(dòng)條件下,分別進(jìn)行常規(guī)稻及雜交稻供種過程模擬,從種箱內(nèi)的種子流速變化及分布“蓬松程度”等方面分析,驗(yàn)證了交叉導(dǎo)流式結(jié)構(gòu)能夠?yàn)榕欧N輪提供連續(xù)穩(wěn)定的種子流,揭示了振動(dòng)作用穩(wěn)定種子流速和提高機(jī)構(gòu)供種頻率的機(jī)理,表明該定量供種機(jī)構(gòu)較優(yōu)的供種性能能夠滿足精密播種需求。

    2.2 振動(dòng)勻種機(jī)構(gòu)的播種性能仿真

    振動(dòng)勻種機(jī)構(gòu)的種槽板對種子的限位作用決定了種子的有序排隊(duì)狀態(tài),進(jìn)而影響條播或穴播時(shí)種子流與秧盤的對行或?qū)ρㄐ阅堋榱诉_(dá)到較好的播種性能,播種時(shí)到達(dá)種槽出流口的種子最好垂直落入該種槽對應(yīng)的秧盤行或秧穴中,尤其當(dāng)雜交稻低播量播種時(shí)需要種子在種槽板上保持沿長軸方向有序排隊(duì),同時(shí)避免劇烈跳動(dòng)或擺動(dòng),防止掉落過程中發(fā)生彈跳落入相鄰的秧盤行或秧穴中。

    為探究不同類型與參數(shù)的種槽板對種子流速、排隊(duì)狀態(tài)及播種性能的影響,設(shè)計(jì)了9組種槽板,即不同種槽夾角θ的V型種槽板:V-60°、V-75°、V-90°、V-105°、V-120°、V-135°、V-150°,和不同半徑的U型種槽板:U1(半徑為6.30 mm)、U2(半徑為7.00 mm),進(jìn)行水稻精少量播種性能模擬試驗(yàn)。旨在篩選種槽板因素水平,為裝置優(yōu)化試驗(yàn)提供理論依據(jù),減少試驗(yàn)次數(shù),節(jié)省時(shí)間和成本。

    仿真所需材料的物性參數(shù)及接觸參數(shù)參照2.1節(jié)進(jìn)行設(shè)置,秧盤運(yùn)動(dòng)設(shè)置為勻速行走,速度為0.083 m/s(播種流水線生產(chǎn)率500盤/h),振動(dòng)種盤運(yùn)動(dòng)參考文獻(xiàn)[14]進(jìn)行設(shè)置,仿真時(shí)間為12 s。

    2.2.1不同種槽板對種子流速影響

    統(tǒng)計(jì)在播種過程穩(wěn)定后4~10 s內(nèi),9組仿真試驗(yàn)種槽板上的種子流動(dòng)速度,每組取100個(gè)數(shù)據(jù),繪制種子速度箱線圖,如圖8所示。通過觀察箱線圖中每組數(shù)據(jù)箱體和線段的長短,可得到數(shù)據(jù)的四分位距大小,以及正常值的分布范圍,從而直觀地判斷不同種槽板對種子流速及穩(wěn)定性的影響。

    圖8 不同種槽板上種子速度箱線圖Fig.8 Box plot of seeds velocity on different striper tanks

    從圖8中各組箱線圖的高度分布來看,V型種槽板夾角θ從60°到150°,各組箱線的高度、速度最大值和中位數(shù)隨θ增大而下降,表明隨著種槽板夾角θ的增大,種子在種槽板上的流動(dòng)速度整體呈減小的趨勢;對于U型種槽板,U1型種槽板上的種子流速略大于U2型種槽板,表明種槽半徑較大時(shí)種子流速較慢。從各組箱線長度和異常值分布情況可以看出,隨著種槽板夾角θ的增大,種子速度分布范圍逐漸集中,種子速度異常值逐漸減少,且在種槽板夾角θ大于等于120°時(shí),種子流速分布范圍相近,表明隨著種槽板夾角θ的增大,種子流速趨于穩(wěn)定,且在θ達(dá)到120°時(shí)種子流速穩(wěn)定性隨槽板夾角的增大變化不明顯。兩種U型種槽板上種子速度異常值少,表明種子在U型槽上的流速較V型槽穩(wěn)定。

    2.2.2不同種槽板對種子排隊(duì)狀態(tài)影響

    觀察種子在種槽板上的排隊(duì)狀態(tài)及到達(dá)種槽板出流口落入秧盤內(nèi)的掉落狀態(tài),分析播種過程中不同種槽板對種子的限位作用和對行、對穴性能影響。圖9為9組種槽板上種子的排隊(duì)狀態(tài),從圖中可以看出,對于V型種槽板,當(dāng)種槽板夾角θ較小(60°、75°和90°)時(shí),種槽板對種子的限位作用強(qiáng),此時(shí)種槽內(nèi)的種子可沿其長軸方向有序排隊(duì),并形成較均勻的單粒到1.5粒種子流;隨著種槽板夾角θ增大,種槽板對種子的限位作用逐漸減弱,當(dāng)θ為105°和120°時(shí),種子與種槽的接觸面積增大,種槽內(nèi)的種子可前后交疊排列形成較均勻的1.5粒到2粒種子流;當(dāng)θ為135°和150°時(shí),種槽板夾角過大,此時(shí)種槽板對種子的限位作用弱,無法引導(dǎo)種子在槽板內(nèi)有序排隊(duì)形成均勻的種子流。對于U型槽板,種子與種槽板的接觸面積大,種子排列較不規(guī)律,當(dāng)U型槽半徑較小(U1半徑為6.30 mm)時(shí)可形成較均勻的1.5粒到2粒種子流,當(dāng)U型槽半徑較大(U2半徑為7.00 mm)時(shí)形成無序的2粒種子流。

    圖9 種子在不同種槽板上的排隊(duì)狀態(tài)Fig.9 Queue conditions of seeds on different striper tanks

    此外,通過觀察種子在不同種槽板出流口落入秧盤時(shí)的形態(tài),得到V-60°、V-75°和V-90°種槽板中,由于種槽板的限位作用好,種子能夠垂直掉落到下方的秧盤中;V-105°和V-120°種槽板限位性一般,但由于種子流速較穩(wěn)定,種子隨槽板振動(dòng)作用前后擺動(dòng)行走,也能夠垂直掉落至秧盤中,較準(zhǔn)確地實(shí)現(xiàn)對行對穴;而V-135°、V-150°、U1和U2型種槽板,由于限位性差,種子隨槽板振動(dòng)作用產(chǎn)生小幅度的左右擺動(dòng),掉落時(shí)易發(fā)生翻滾,對行對穴性能較差。

    2.2.3不同種槽板對播種性能影響的仿真分析

    為進(jìn)一步確定播種性能較優(yōu)的種槽板,統(tǒng)計(jì)不同種槽板的播種合格指數(shù)與空穴指數(shù)等性能指標(biāo)。每組仿真試驗(yàn)分別播種3個(gè)秧盤,統(tǒng)計(jì)秧盤中每穴內(nèi)的種子數(shù)量,評(píng)價(jià)指標(biāo)為

    (12)

    式中Y1——播種合格指數(shù)

    Y2——空穴指數(shù)

    N1——合格穴數(shù),合格區(qū)間為1~4 粒/穴

    N2——空穴數(shù),無種子為空穴

    N——總秧穴數(shù)

    統(tǒng)計(jì)結(jié)果見表3。

    表3 不同種槽板播種性能仿真試驗(yàn)結(jié)果Tab.3 Simulation result of sowing performance on different striper tanks

    從表3試驗(yàn)結(jié)果可以看出,V-60°和V-75°種槽板的播種空穴指數(shù)高,不能滿足水稻精密播種要求,其原因是在種槽板上形成的是單粒種子流,存在斷流的位置播入秧盤中即形成了空穴;V-135°、V-150°和U2型種槽板空穴指數(shù)高,其原因是槽板的限位性差,不能形成均勻的種子流;V-90°、V-105°、V-120°和U1型種槽板的播種合格指數(shù)均超過90%,空穴指數(shù)較低,平均每穴粒數(shù)滿足播種要求,且標(biāo)準(zhǔn)差較小,表明較穩(wěn)定的種子流速和連續(xù)均勻的種子流是實(shí)現(xiàn)水稻精少量播種的重要條件。

    綜上對種子在不同種槽板上的流速及排隊(duì)狀態(tài)分析,以播種合格指數(shù)和空穴指數(shù)作為評(píng)判播種性能的指標(biāo),優(yōu)選了播種效果較好的V-90°、V-105°、V-120°和U1型種槽板,為進(jìn)行機(jī)構(gòu)的播種性能優(yōu)化試驗(yàn)提供了參考。

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

    分別進(jìn)行定量供種機(jī)構(gòu)在種箱振動(dòng)板有振動(dòng)與無振動(dòng)條件下的供種頻率對比試驗(yàn),以及不同種槽板播種性能試驗(yàn),以驗(yàn)證仿真的準(zhǔn)確性,確定最優(yōu)種槽板能夠提高精密播種裝置的性能。試驗(yàn)地點(diǎn)在華南農(nóng)業(yè)大學(xué)農(nóng)業(yè)工程實(shí)驗(yàn)室,采用研制的振動(dòng)式精密播種裝置,安裝在水稻秧盤育秧精密播種流水線上,裝置如圖10所示。試驗(yàn)所用水稻品種:雜交稻品種為培雜泰豐,種子含水率約為27.5%,千粒質(zhì)量為28.20 g;常規(guī)稻品種為華航38號(hào),種子含水率約為27.38%,千粒質(zhì)量為27.59 g。兩個(gè)水稻品種均進(jìn)行催芽處理,種子破胸露白后再進(jìn)行播種試驗(yàn)。

    圖10 振動(dòng)式水稻精密播種裝置實(shí)物圖Fig.10 Vibration rice precision sowing device

    3.1 定量供種機(jī)構(gòu)供種頻率對比試驗(yàn)

    為探明振動(dòng)作用對定量供種機(jī)構(gòu)供種穩(wěn)定性和供種頻率的影響,以種箱振動(dòng)板的振動(dòng)作用作為試驗(yàn)因素,對供種機(jī)構(gòu)有、無振動(dòng)時(shí)的供種頻率進(jìn)行測定分析。有振動(dòng)時(shí)開啟電磁振動(dòng)器并調(diào)節(jié)振動(dòng)板產(chǎn)生0.015 mm的振幅(可通過振動(dòng)測試得到),無振動(dòng)時(shí)關(guān)閉電磁振動(dòng)器,分別進(jìn)行定量供種機(jī)構(gòu)常規(guī)稻與雜交稻供種試驗(yàn)。常規(guī)稻供種時(shí)排種輪轉(zhuǎn)速為7.0 r/min,雜交稻低播量時(shí)轉(zhuǎn)速為3.3 r/min。

    供種頻率測定方法為:在機(jī)構(gòu)工作穩(wěn)定后,有、無振動(dòng)條件下分別接取排種輪轉(zhuǎn)動(dòng)60 s所排出的種子質(zhì)量,使用精度0.01 g的電子秤進(jìn)行稱量,計(jì)算排種輪的供種頻率,每組試驗(yàn)重復(fù)6次,求取平均供種頻率;供種穩(wěn)定性的測定方法為:接取排種輪轉(zhuǎn)動(dòng)5 s所排出的種子質(zhì)量并稱量,有、無振動(dòng)條件下各重復(fù)5次,求取供種頻率變異系數(shù)。試驗(yàn)結(jié)果見表4。

    表4 定量供種機(jī)構(gòu)供種頻率對比試驗(yàn)結(jié)果Tab.4 Feeding frequency comparison test result of quantitative-seeds supply device

    由表4可知,常規(guī)稻供種時(shí)電磁振動(dòng)作用使供種頻率提高3.64%,雜交稻低播量時(shí)提高5.52%,與2.1.2節(jié)仿真結(jié)果進(jìn)行對比,仿真相對誤差分別為9.07%和1.99%,表明仿真結(jié)果精度較高。常規(guī)稻供種時(shí)有、無振動(dòng)條件下的供種頻率變異系數(shù)均較低播量供種時(shí)小,表明常規(guī)稻供種時(shí)排種輪單位時(shí)間內(nèi)排出的稻種量更加穩(wěn)定;在振動(dòng)板有振動(dòng)作用時(shí),可提高常規(guī)稻供種穩(wěn)定性2.47%,雜交稻供種穩(wěn)定性提高1.33%,表明振動(dòng)板振動(dòng)對常規(guī)稻供種穩(wěn)定性提高作用更明顯。

    3.2 不同種槽板播種性能試驗(yàn)

    以種槽板類型、播種量、氣壓作為試驗(yàn)因素進(jìn)行全因素播種性能試驗(yàn)。仿真試驗(yàn)優(yōu)選的種槽板類型因素(A)水平為:V-90°、V-105°、V-120°和U1型種槽板;播種量因素(B)水平為:常規(guī)稻(80 g/盤)和雜交稻(40 g/盤);常用氣壓(C)水平為:0.24、0.26、0.28 MPa,共進(jìn)行24組試驗(yàn)。

    水稻秧盤育秧精密播種流水線播種生產(chǎn)率為500盤/h,秧盤為14×29缽體毯狀軟盤。每組試驗(yàn)播種6盤,每盤隨機(jī)選取2個(gè)區(qū)域,區(qū)域1橫向完整2行(14格×2行)和區(qū)域2橫向6行、縱向7列(6格×7格),統(tǒng)計(jì)并計(jì)算每組的播種合格指數(shù)、空穴指數(shù)、平均每穴粒數(shù)和標(biāo)準(zhǔn)差。雜交稻低播量播種合格區(qū)間為1~4粒/穴,常規(guī)稻大播量播種合格區(qū)間為3~7粒/穴。分別計(jì)算4組種槽板在播種雜交稻和常規(guī)稻時(shí)的性能指標(biāo)平均值k,試驗(yàn)結(jié)果見表5,播種效果如圖11所示。

    通過表5播種性能試驗(yàn)結(jié)果比較,雜交稻低播量播種時(shí)和常規(guī)稻大播量播種時(shí),試驗(yàn)因素A在k3水平下的播種平均合格指數(shù)最高、平均空穴指數(shù)最低,平均每穴粒數(shù)滿足播種要求且標(biāo)準(zhǔn)差較小,表明兩種播量下V-120°種槽板的播種性能均最優(yōu),與2.2.2節(jié)仿真試驗(yàn)優(yōu)選結(jié)果一致。同時(shí),采用該種槽板進(jìn)行低播量播種,當(dāng)氣壓為0.26 MPa時(shí)播種合格指數(shù)為93.91%,空穴指數(shù)為0.94%;常規(guī)稻播種時(shí)當(dāng)氣壓為0.28 MPa,播種合格指數(shù)為96.10%,空穴指數(shù)為0,滿足常規(guī)稻與雜交稻精少量播種育秧技術(shù)要求,優(yōu)化了振動(dòng)式水稻精密播種裝置的播種性能。

    表5 播種性能試驗(yàn)結(jié)果Tab.5 Test result of sowing performance

    圖11 播種效果Fig.11 Seeding effect

    4 結(jié)論

    (1)采用離散元法進(jìn)行了有無振動(dòng)時(shí)定量供種機(jī)構(gòu)的工作過程仿真研究,得出在振動(dòng)條件下進(jìn)行供種,可提高供種頻率和穩(wěn)定性。通過分析種箱內(nèi)不同種層的種子流動(dòng)特性,驗(yàn)證了交叉導(dǎo)流式振動(dòng)種箱可為排種輪提供連續(xù)穩(wěn)定的種子流;通過對比分析充填區(qū)種子的流速與“蓬松”程度得知,振動(dòng)作用能夠有效地穩(wěn)定種子流速和分布密度,提高供種頻率。

    (2)進(jìn)行了不同類型與參數(shù)種槽板播種性能仿真試驗(yàn),得到了種子在不同種槽板上的流速規(guī)律。通過分析不同種槽板對種子流速、有序排列及播種性能的影響可知,種子流速隨著種槽板夾角θ增大而減小并趨于穩(wěn)定,且當(dāng)θ達(dá)到120°時(shí),流速穩(wěn)定性隨θ增大的變化不明顯;以播種合格指數(shù)和空穴指數(shù)為指標(biāo),篩選出了V-90°、V-105°、V-120°和U1型4組種槽板的播種效果較好,為裝置優(yōu)化試驗(yàn)減少了試驗(yàn)次數(shù)。

    (3)進(jìn)行了定量供種機(jī)構(gòu)的供種頻率試驗(yàn)和精密播種裝置性能優(yōu)化試驗(yàn),確定了播種效果最佳的種槽板為V-120°種槽板。通過定量供種機(jī)構(gòu)試驗(yàn),得到常規(guī)稻、雜交稻供種時(shí)電磁振動(dòng)作用分別可提高供種頻率3.64%和5.52%,驗(yàn)證了定量供種機(jī)構(gòu)在有無振動(dòng)作用下的仿真結(jié)果誤差較小;通過精密播種裝置性能優(yōu)化試驗(yàn),得到采用V-120°種槽板,以0.26 MPa氣壓進(jìn)行低播量播種時(shí),合格指數(shù)為93.91%,空穴指數(shù)為0.94%,以0.28 MPa氣壓進(jìn)行常規(guī)稻播種時(shí),合格指數(shù)為96.10%,空穴指數(shù)為0,滿足常規(guī)稻與雜交稻精密播種育秧技術(shù)要求,為水稻精密播種裝置的研究提供了基礎(chǔ)。

    1 李澤華,馬旭,李秀昊,等.水稻栽植機(jī)械化技術(shù)研究進(jìn)展[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2018,49(5):1-20. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20180501&flag=1. DOI:10.6041/j.issn.1000-1298.2018.05.001.

    LI Zehua,MA Xu,LI Xiuhao,et al.Research progress of rice transplanting mechanization[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2018,49(5):1-20. (in Chinese)

    2 YOO S N, CHOI Y S, SUH S R. Development of a precision seed metering device for direct seeding of rice [J]. Journal of Biosystems Engineering, 2005,30(5):261-267.

    3 李洪昌,高芳,趙湛,等.國內(nèi)外精密排種器研究現(xiàn)狀與發(fā)展趨勢[J] .中國農(nóng)機(jī)化學(xué)報(bào),2014,35(2): 12-16,56.

    LI Hongchang, GAO Fang, ZHAO Zhan, et al. Domestic and overseas research status and development trend of precision seed-metering device [J]. Journal of Chinese Agricultural Mechanization, 2014, 35(2): 12-16,56. (in Chinese)

    4 李耀明,徐立章,向忠平,等.日本水稻種植機(jī)械化技術(shù)的最新研究進(jìn)展[J].農(nóng)業(yè)工程學(xué)報(bào),2005,21(11):182-185.

    LI Yaoming, XU Lizhang, XIANG Zhongping, et al. Research advances of rice planting mechanization in Japan [J]. Transactions of the CSAE, 2005, 21(11): 182-185. (in Chinese)

    5 孫勇飛,吳崇友,張文毅,等.水稻育秧播種機(jī)的發(fā)展概況與趨勢[J]. 農(nóng)機(jī)化研究,2013,35(12):210-215.

    SUN Yongfei, WU Chongyou, ZHANG Wenyi, et al. Development situation and tendency of rice seedling planter [J]. Journal of Agricultural Mechanization Research, 2013,35(12):210-215. (in Chinese)

    6 李志偉,邵耀堅(jiān).電磁振動(dòng)式水稻穴盤精量播種機(jī)的設(shè)計(jì)與試驗(yàn)[J] .農(nóng)業(yè)機(jī)械學(xué)報(bào),2000,31(5): 32-34.

    LI Zhiwei, SHAO Yaojian. Study and test of electromagnetic vibrating type rice seeder for hill seedling nursery box [J]. Transactions of the Chinese Society for Agricultural Machinery, 2000, 31(5): 32-34. (in Chinese)

    7 張石平,陳進(jìn),李耀明,等.振動(dòng)氣吸式穴盤精播裝置振動(dòng)條件理論分析與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2008,39(7): 56-59,79.

    ZHANG Shiping, CHEN Jin, LI Yaoming, et al. Theoretical analysis and experiment on vibration conditions for vibrational air-suction tray precision seeding machine [J]. Transactions of the Chinese Society for Agricultural Machinery, 2008, 39(7): 56-59,79. (in Chinese)

    8 張學(xué)義,邵耀堅(jiān).電磁振動(dòng)排種器種子運(yùn)動(dòng)分析及試驗(yàn)研究[J] .農(nóng)業(yè)工程學(xué)報(bào),1996,12(1): 81-86.

    ZHANG Xueyi, SHAO Yaojian. Theoretical analysis and experiment on the kinematics of feeding device under electromagnetic vibration [J]. Transactions of the CSAE, 1996, 12(1): 81-86. (in Chinese)

    9 楊望,楊堅(jiān),周曉蓉,等.電磁振動(dòng)排種器振動(dòng)系統(tǒng)參數(shù)的確定[J] .農(nóng)機(jī)化研究,2009,31(2): 47-50,54.

    YANG Wang, YANG Jian, ZHOU Xiaorong, et al. Parameters study on the vibrant system electro magnetic vibrated seeding apparatus [J]. Journal of Agricultural Mechanization Research, 2009, 31(2): 47-50,54. (in Chinese)

    10 張斌,俞亞新,趙勻.超級(jí)稻稻種定向播種的原理研究[J].浙江理工大學(xué)學(xué)報(bào),2008,25(4): 454-456.

    ZHANG Bin, YU Yaxin, ZHAO Yun. Study on principle of the rice seeds embryo oriented seeding [J]. Journal of Zhejiang Sci-Tech University, 2008, 25(4): 454-456. (in Chinese)

    11 齊龍,譚祖庭,馬旭,等.氣動(dòng)振動(dòng)式勻種裝置工作參數(shù)的優(yōu)化及試驗(yàn)[J].吉林大學(xué)學(xué)報(bào):工學(xué)版,2014,44(6): 1684-1691.

    QI Long, TAN Zuting, MA Xu, et al. Optimization and test of operational parameters of pneumatic vibration uniform-seeds device [J]. Journal of Jilin University: Engineering and Technology Edition, 2014, 44(6): 1684-1691. (in Chinese)

    12 ZHU H P, ZHOU Z Y, YANG R Y, et al. Discrete particle simulation of particulate systems: a review of major applications and findings [J] . Chemical Engineering Science, 2008, 63(23): 5728-5770.

    13 YU A B. Discrete element method: an effective method for particle scale research of particulate matter [J] . Engineering Computations, 2004, 21(2-4): 205-214.

    14 鹿芳媛,馬旭,齊龍,等.基于離散元法的雜交稻振動(dòng)勻種裝置參數(shù)優(yōu)化與試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(10): 17-25.

    LU Fangyuan, MA Xu, QI Long, et al. Parameter optimization and experiment of vibration seed-uniforming device for hybrid rice based on discrete element method [J]. Transactions of the CSAE, 2016, 32(10): 17-25. (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,100.

    CHEN Jin, ZHOU Han, ZHAO Zhan, et al. Analysis on rice seeds motion on vibrating plate using EDEM [J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(10): 79-83, 100. (in Chinese)

    16 玉大略,馬旭,周海波,等.振動(dòng)流動(dòng)式水稻秧盤育秧播種器的試驗(yàn)研究[J].農(nóng)機(jī)化研究,2010,32(2): 137-139.

    YU Dalue, MA Xu, ZHOU Haibo, et al. Experimental research of vibration flow type rice seedlings seeder [J]. Journal of Agricultural Mechanization Research, 2010, 32(2): 137-139. (in Chinese)

    17 周海波,梁秋艷,魏天路,等.雙級(jí)振動(dòng)精密排種器外槽輪式定量供種裝置設(shè)計(jì)與試驗(yàn)[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2016, 47(增刊): 57-61,83.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=2016s009&flag=1.DOI:10.6041/j.issn.1000-1298.2016.S0.009.

    ZHOU Haibo, LIANG Qiuyan, WEI Tianlu, et al. Design and experiment of quantitative seed supply device with fluted roller used for double-vibrating precision seed meter [J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(Supp.): 57-61,83. (in Chinese)

    18 胡繼云,廉振紅,邵洪濤.單自由度簡諧振動(dòng)機(jī)械彈性支承剛度和激振力的設(shè)計(jì)[J].河南工業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2006,27(4): 11-14.

    HU Jiyun, LIAN Zhenhong, SHAO Hongtao. The design of elastic supporting stiffness and exciting force for simple harmonic vibration machine with single-degree-of-freedom [J]. Journal of Henan University of Technology: Natural Science Edition, 2006, 27(4): 11-14. (in Chinese)

    19 劉習(xí)軍,賈啟芬.工程振動(dòng)理論與測試技術(shù)[M].北京:高等教育出版社,2004: 39-41.

    20 鹿芳媛,馬旭,譚穗妍,等.水稻芽種離散元主要接觸參數(shù)仿真標(biāo)定與試驗(yàn)[J/OL]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2018,49(2):93-99. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20180212&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2018.02.012.

    LU Fangyuan, MA Xu, TAN Suiyan, et al.Simulative calibration and experiment of rice bud seeds discrete elements main contact parameters [J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2018,49(2):93-99. (in Chinese)

    21 陳長冰.基于整體流型的粉體料倉設(shè)計(jì)分析[J].化工設(shè)備與管道,2006,43(3): 34-38.

    CHEN Changbing. Design and analysis of silos with mass flow [J]. Process Equipment and Piping, 2006, 43(3): 34-38. (in Chinese)

    猜你喜歡
    供種雜交稻定量
    油菜旋轉(zhuǎn)盤式高速集排器螺旋供種裝置設(shè)計(jì)與試驗(yàn)
    基于振動(dòng)盤供種的混編景觀種子帶編織機(jī)設(shè)計(jì)
    稻麥兼用排種器振動(dòng)供種裝置性能試驗(yàn)
    顯微定量法鑒別林下山參和園參
    當(dāng)歸和歐當(dāng)歸的定性與定量鑒別
    中成藥(2018年12期)2018-12-29 12:25:44
    傾斜拋物線型孔輪式小麥供種裝置設(shè)計(jì)與試驗(yàn)
    10 種中藥制劑中柴胡的定量測定
    中成藥(2017年6期)2017-06-13 07:30:35
    樹脂包膜緩釋肥在機(jī)插雜交稻上的肥效初探
    慢性HBV感染不同狀態(tài)下HBsAg定量的臨床意義
    石門縣超級(jí)雜交稻“種三產(chǎn)四”高產(chǎn)栽培技術(shù)
    作物研究(2014年6期)2014-03-01 03:39:11
    成年人免费黄色播放视频| 成人国语在线视频| 国产xxxxx性猛交| 国产高清videossex| 啦啦啦 在线观看视频| 丁香欧美五月| 黄网站色视频无遮挡免费观看| 侵犯人妻中文字幕一二三四区| 国产乱人伦免费视频| 日韩精品中文字幕看吧| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久久久久久大奶| 18禁美女被吸乳视频| 亚洲一区二区三区欧美精品| 如日韩欧美国产精品一区二区三区| 成人18禁在线播放| 中文亚洲av片在线观看爽| 男女下面进入的视频免费午夜 | 久久欧美精品欧美久久欧美| 国产av精品麻豆| 国产欧美日韩一区二区精品| 色综合婷婷激情| 黄色 视频免费看| 日本黄色视频三级网站网址| 99在线视频只有这里精品首页| 交换朋友夫妻互换小说| 18禁美女被吸乳视频| 精品久久久久久久毛片微露脸| 叶爱在线成人免费视频播放| 日本免费a在线| 一个人观看的视频www高清免费观看 | 免费不卡黄色视频| 亚洲av片天天在线观看| 一本综合久久免费| 欧美日韩av久久| 麻豆久久精品国产亚洲av | 怎么达到女性高潮| 久久人人97超碰香蕉20202| 老汉色av国产亚洲站长工具| 精品国产一区二区久久| 欧美 亚洲 国产 日韩一| 亚洲少妇的诱惑av| 欧美日韩乱码在线| www.精华液| 日韩精品中文字幕看吧| 午夜免费观看网址| 成人手机av| 淫妇啪啪啪对白视频| 精品一区二区三区视频在线观看免费 | 久久久久久久久久久久大奶| 丰满迷人的少妇在线观看| 一边摸一边做爽爽视频免费| 日本精品一区二区三区蜜桃| 日本一区二区免费在线视频| 在线观看舔阴道视频| 99久久人妻综合| 午夜福利在线免费观看网站| 老司机午夜福利在线观看视频| 高清欧美精品videossex| 亚洲av成人不卡在线观看播放网| 色综合站精品国产| bbb黄色大片| 亚洲中文日韩欧美视频| 亚洲自拍偷在线| 亚洲人成77777在线视频| 很黄的视频免费| 亚洲少妇的诱惑av| 一级黄色大片毛片| 亚洲午夜精品一区,二区,三区| 久久精品亚洲熟妇少妇任你| 亚洲熟女毛片儿| 欧美一区二区精品小视频在线| 久久久国产欧美日韩av| 国产精品影院久久| 精品一区二区三卡| 亚洲激情在线av| 天天添夜夜摸| 亚洲av熟女| 中文字幕人妻熟女乱码| 精品欧美一区二区三区在线| 亚洲色图综合在线观看| 欧美精品亚洲一区二区| 亚洲美女黄片视频| 久久精品91无色码中文字幕| 黄色怎么调成土黄色| 欧美亚洲日本最大视频资源| 五月开心婷婷网| 无人区码免费观看不卡| 免费在线观看视频国产中文字幕亚洲| 欧美午夜高清在线| 99久久国产精品久久久| 国产亚洲精品综合一区在线观看 | 日韩高清综合在线| 视频在线观看一区二区三区| 久久精品成人免费网站| 男人的好看免费观看在线视频 | 亚洲情色 制服丝袜| 丰满迷人的少妇在线观看| 午夜日韩欧美国产| 久久国产精品人妻蜜桃| 亚洲色图综合在线观看| 色尼玛亚洲综合影院| 亚洲欧美日韩高清在线视频| 不卡av一区二区三区| 久久久久久免费高清国产稀缺| 在线天堂中文资源库| 国产成人av教育| 国产熟女xx| 国产在线观看jvid| 麻豆av在线久日| 久久热在线av| 亚洲全国av大片| 亚洲av电影在线进入| 少妇的丰满在线观看| 国产精品乱码一区二三区的特点 | 啪啪无遮挡十八禁网站| 我的亚洲天堂| av中文乱码字幕在线| 欧美激情高清一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲精品在线观看二区| 亚洲黑人精品在线| 看免费av毛片| 欧美日本亚洲视频在线播放| 丝袜人妻中文字幕| 欧美日韩国产mv在线观看视频| 丁香欧美五月| 超色免费av| 日韩大尺度精品在线看网址 | 国产av精品麻豆| 成人三级做爰电影| 亚洲精品中文字幕在线视频| 琪琪午夜伦伦电影理论片6080| 两人在一起打扑克的视频| 女性被躁到高潮视频| 国产精品国产av在线观看| 久久性视频一级片| 午夜福利影视在线免费观看| 91麻豆精品激情在线观看国产 | 天天影视国产精品| 波多野结衣一区麻豆| 欧美日韩国产mv在线观看视频| 国产又色又爽无遮挡免费看| 国产欧美日韩综合在线一区二区| 久久久久国产精品人妻aⅴ院| 1024香蕉在线观看| aaaaa片日本免费| 不卡av一区二区三区| 欧美午夜高清在线| 亚洲一区高清亚洲精品| 国产无遮挡羞羞视频在线观看| 久久精品人人爽人人爽视色| 亚洲avbb在线观看| 久久久久久人人人人人| 亚洲五月天丁香| xxxhd国产人妻xxx| 亚洲片人在线观看| 操出白浆在线播放| 水蜜桃什么品种好| 国产成人精品在线电影| 精品久久久久久成人av| 国产成人欧美| 精品免费久久久久久久清纯| 啦啦啦在线免费观看视频4| 天天躁狠狠躁夜夜躁狠狠躁| 久久性视频一级片| 又黄又爽又免费观看的视频| 亚洲,欧美精品.| 国产成人一区二区三区免费视频网站| 国产精品免费视频内射| 国产精品九九99| 欧美黄色淫秽网站| 悠悠久久av| 免费观看人在逋| 国产精品一区二区在线不卡| 日韩成人在线观看一区二区三区| 欧美色视频一区免费| 香蕉久久夜色| 精品国产美女av久久久久小说| 女警被强在线播放| 亚洲五月色婷婷综合| 久久久久久久午夜电影 | 免费高清视频大片| 久久伊人香网站| 久久久久久人人人人人| 这个男人来自地球电影免费观看| 欧美乱色亚洲激情| 国产黄a三级三级三级人| 国产亚洲欧美98| 一二三四在线观看免费中文在| 欧美激情极品国产一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 97碰自拍视频| a级毛片在线看网站| 在线观看免费高清a一片| 日韩 欧美 亚洲 中文字幕| 18禁美女被吸乳视频| 免费久久久久久久精品成人欧美视频| av天堂在线播放| av网站免费在线观看视频| xxxhd国产人妻xxx| 免费av毛片视频| 日韩一卡2卡3卡4卡2021年| 777久久人妻少妇嫩草av网站| 欧美日本中文国产一区发布| svipshipincom国产片| 天天躁夜夜躁狠狠躁躁| 嫁个100分男人电影在线观看| 亚洲av片天天在线观看| 亚洲精品美女久久久久99蜜臀| 午夜老司机福利片| 亚洲欧美精品综合一区二区三区| 国产一区二区三区综合在线观看| 正在播放国产对白刺激| 日本欧美视频一区| 美女 人体艺术 gogo| 又紧又爽又黄一区二区| 国产成人精品久久二区二区91| 看免费av毛片| av视频免费观看在线观看| 精品无人区乱码1区二区| 亚洲专区国产一区二区| 在线看a的网站| 一级毛片精品| 在线永久观看黄色视频| 一级片'在线观看视频| 老司机靠b影院| 久久久久久久精品吃奶| 丁香六月欧美| a级片在线免费高清观看视频| 亚洲欧美激情综合另类| 亚洲精品成人av观看孕妇| 天天添夜夜摸| 看免费av毛片| 国产欧美日韩一区二区三| 女人精品久久久久毛片| www.自偷自拍.com| 两性夫妻黄色片| 欧美激情久久久久久爽电影 | 精品一区二区三区四区五区乱码| 99久久久亚洲精品蜜臀av| 交换朋友夫妻互换小说| 欧美 亚洲 国产 日韩一| 一二三四社区在线视频社区8| 熟女少妇亚洲综合色aaa.| 男人舔女人的私密视频| 男人操女人黄网站| 亚洲 欧美 日韩 在线 免费| 在线看a的网站| 99国产精品99久久久久| 亚洲精品美女久久av网站| 亚洲avbb在线观看| 亚洲美女黄片视频| 极品人妻少妇av视频| 777久久人妻少妇嫩草av网站| 亚洲av电影在线进入| 国产精品一区二区免费欧美| 欧美黑人精品巨大| 久久久久国产一级毛片高清牌| 欧美日韩亚洲高清精品| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 久久亚洲真实| 午夜91福利影院| 欧美激情 高清一区二区三区| 久久精品国产99精品国产亚洲性色 | 一级毛片精品| 国产精品免费一区二区三区在线| www.熟女人妻精品国产| 国产一区在线观看成人免费| 黄色片一级片一级黄色片| 国产亚洲精品久久久久久毛片| а√天堂www在线а√下载| 久久 成人 亚洲| 国产精品免费一区二区三区在线| 黄色丝袜av网址大全| 国产三级在线视频| 午夜精品在线福利| 看免费av毛片| 日韩成人在线观看一区二区三区| 亚洲一区二区三区欧美精品| 又黄又爽又免费观看的视频| 亚洲av第一区精品v没综合| 91精品国产国语对白视频| 国产精品免费一区二区三区在线| 国产av一区在线观看免费| 国产精品久久久av美女十八| 久久久国产成人精品二区 | 狂野欧美激情性xxxx| 好男人电影高清在线观看| 黄片大片在线免费观看| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 亚洲欧美一区二区三区久久| 久久久精品欧美日韩精品| av国产精品久久久久影院| 一进一出好大好爽视频| 亚洲精品国产色婷婷电影| 国产xxxxx性猛交| e午夜精品久久久久久久| 纯流量卡能插随身wifi吗| 久久久久九九精品影院| 亚洲专区国产一区二区| 久久天躁狠狠躁夜夜2o2o| 欧美日韩瑟瑟在线播放| 日日干狠狠操夜夜爽| 日韩欧美在线二视频| 亚洲成人免费av在线播放| 99精品久久久久人妻精品| 亚洲 欧美 日韩 在线 免费| 大香蕉久久成人网| 亚洲中文字幕日韩| 久久久久九九精品影院| 亚洲中文字幕日韩| 久久精品亚洲av国产电影网| 成人特级黄色片久久久久久久| 夜夜爽天天搞| 岛国在线观看网站| 免费人成视频x8x8入口观看| 女性被躁到高潮视频| 欧美激情高清一区二区三区| 高清欧美精品videossex| 免费高清在线观看日韩| 91九色精品人成在线观看| 久久香蕉国产精品| 一级毛片精品| 美女高潮到喷水免费观看| 亚洲在线自拍视频| 人人妻人人爽人人添夜夜欢视频| 丰满迷人的少妇在线观看| 在线永久观看黄色视频| 男女高潮啪啪啪动态图| 久久精品国产清高在天天线| 国产高清videossex| 精品一品国产午夜福利视频| 午夜久久久在线观看| 久久狼人影院| 高清av免费在线| 天天躁狠狠躁夜夜躁狠狠躁| 大码成人一级视频| 精品国产一区二区三区四区第35| 亚洲avbb在线观看| 国产免费av片在线观看野外av| 成人亚洲精品av一区二区 | 丰满迷人的少妇在线观看| 久久天堂一区二区三区四区| 久久人妻av系列| 18禁国产床啪视频网站| 久久人妻av系列| 宅男免费午夜| 人人澡人人妻人| 欧美激情 高清一区二区三区| 正在播放国产对白刺激| a级片在线免费高清观看视频| 搡老乐熟女国产| 午夜免费观看网址| 国产成人av激情在线播放| 99热只有精品国产| 国产99久久九九免费精品| 亚洲成av片中文字幕在线观看| 黑人欧美特级aaaaaa片| 欧美黑人欧美精品刺激| 亚洲色图综合在线观看| 午夜福利免费观看在线| 99久久99久久久精品蜜桃| 亚洲人成77777在线视频| 亚洲av片天天在线观看| 99精国产麻豆久久婷婷| 久久精品亚洲熟妇少妇任你| 侵犯人妻中文字幕一二三四区| 电影成人av| 亚洲欧美精品综合一区二区三区| 亚洲人成电影观看| 精品国产超薄肉色丝袜足j| 日韩大码丰满熟妇| 色哟哟哟哟哟哟| a在线观看视频网站| 亚洲成人免费电影在线观看| 欧美成狂野欧美在线观看| 动漫黄色视频在线观看| 久9热在线精品视频| 国产精品综合久久久久久久免费 | 久久久国产一区二区| 精品国产一区二区三区四区第35| 久久久国产成人精品二区 | 丰满迷人的少妇在线观看| 99热只有精品国产| 窝窝影院91人妻| 国产亚洲欧美在线一区二区| 9191精品国产免费久久| 51午夜福利影视在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品一区二区三区在线| 亚洲av五月六月丁香网| 午夜免费激情av| 不卡一级毛片| 热99re8久久精品国产| 日韩精品青青久久久久久| 性少妇av在线| 国产在线观看jvid| 男女高潮啪啪啪动态图| www.999成人在线观看| 日韩大尺度精品在线看网址 | 国产麻豆69| 人人妻,人人澡人人爽秒播| 久久久久久亚洲精品国产蜜桃av| 免费少妇av软件| 欧美日韩亚洲国产一区二区在线观看| 免费av毛片视频| 国产在线精品亚洲第一网站| 天堂√8在线中文| 麻豆一二三区av精品| 亚洲av电影在线进入| netflix在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲男人的天堂狠狠| 欧美黑人精品巨大| 精品午夜福利视频在线观看一区| 欧美+亚洲+日韩+国产| 亚洲九九香蕉| 精品久久久久久久毛片微露脸| 中文字幕最新亚洲高清| 久久欧美精品欧美久久欧美| 国产色视频综合| 亚洲国产欧美网| 88av欧美| 国产av一区在线观看免费| 国产精品久久久av美女十八| bbb黄色大片| 亚洲精品在线美女| 精品福利永久在线观看| 亚洲中文日韩欧美视频| 亚洲熟妇中文字幕五十中出 | 九色亚洲精品在线播放| 大陆偷拍与自拍| 亚洲精品久久午夜乱码| av视频免费观看在线观看| 最新在线观看一区二区三区| 超碰97精品在线观看| 国产精品 国内视频| 99国产精品99久久久久| 亚洲视频免费观看视频| 美女 人体艺术 gogo| 757午夜福利合集在线观看| 90打野战视频偷拍视频| 午夜久久久在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩福利视频一区二区| 午夜激情av网站| 俄罗斯特黄特色一大片| 50天的宝宝边吃奶边哭怎么回事| 欧美成人免费av一区二区三区| 精品一区二区三区四区五区乱码| 亚洲五月色婷婷综合| 露出奶头的视频| 一级毛片女人18水好多| 欧美老熟妇乱子伦牲交| 大码成人一级视频| 日韩成人在线观看一区二区三区| 亚洲成a人片在线一区二区| 亚洲国产欧美日韩在线播放| 中亚洲国语对白在线视频| 色老头精品视频在线观看| 最近最新中文字幕大全免费视频| 91在线观看av| 嫩草影院精品99| 亚洲色图综合在线观看| 日本精品一区二区三区蜜桃| 91大片在线观看| 亚洲成人免费av在线播放| 精品一区二区三区四区五区乱码| 少妇粗大呻吟视频| 国产欧美日韩一区二区精品| 999久久久精品免费观看国产| 亚洲人成电影观看| 侵犯人妻中文字幕一二三四区| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影 | av有码第一页| 脱女人内裤的视频| 久久精品影院6| 怎么达到女性高潮| 亚洲av成人av| 午夜免费鲁丝| 免费在线观看完整版高清| 正在播放国产对白刺激| 精品人妻在线不人妻| 亚洲男人的天堂狠狠| av在线播放免费不卡| 高清av免费在线| av在线播放免费不卡| 黄片播放在线免费| 丝袜美腿诱惑在线| 亚洲人成网站在线播放欧美日韩| 亚洲人成伊人成综合网2020| 亚洲精品久久成人aⅴ小说| 人人妻人人添人人爽欧美一区卜| 大陆偷拍与自拍| 一级a爱视频在线免费观看| 精品熟女少妇八av免费久了| 99re在线观看精品视频| 一级毛片女人18水好多| 国产有黄有色有爽视频| 黄网站色视频无遮挡免费观看| 久久这里只有精品19| 亚洲精品久久成人aⅴ小说| 精品国内亚洲2022精品成人| 日日干狠狠操夜夜爽| 动漫黄色视频在线观看| 女人被狂操c到高潮| 欧美老熟妇乱子伦牲交| 欧美乱码精品一区二区三区| 99riav亚洲国产免费| 99久久人妻综合| 色尼玛亚洲综合影院| 色综合站精品国产| 亚洲精品一卡2卡三卡4卡5卡| 99精品久久久久人妻精品| 丁香欧美五月| 一级毛片女人18水好多| 大型av网站在线播放| 成人av一区二区三区在线看| 高清av免费在线| 啦啦啦 在线观看视频| www日本在线高清视频| 久久精品亚洲av国产电影网| 久久国产乱子伦精品免费另类| 精品第一国产精品| 久久人妻福利社区极品人妻图片| 黄网站色视频无遮挡免费观看| 久久这里只有精品19| 亚洲精品久久成人aⅴ小说| 免费一级毛片在线播放高清视频 | 丰满人妻熟妇乱又伦精品不卡| 成人免费观看视频高清| 久热爱精品视频在线9| 精品无人区乱码1区二区| 另类亚洲欧美激情| 日日爽夜夜爽网站| 亚洲黑人精品在线| 黑人巨大精品欧美一区二区蜜桃| 男女高潮啪啪啪动态图| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 人人妻人人爽人人添夜夜欢视频| 不卡av一区二区三区| 麻豆国产av国片精品| 久久狼人影院| 成人精品一区二区免费| 这个男人来自地球电影免费观看| 涩涩av久久男人的天堂| 欧美久久黑人一区二区| 亚洲精品中文字幕一二三四区| 亚洲成人久久性| 美女扒开内裤让男人捅视频| 视频区欧美日本亚洲| 久久国产精品影院| 操美女的视频在线观看| 亚洲成av片中文字幕在线观看| 老司机亚洲免费影院| 国产熟女午夜一区二区三区| 精品一区二区三卡| 一区二区三区国产精品乱码| 又大又爽又粗| 亚洲国产精品sss在线观看 | 久久中文看片网| 免费在线观看视频国产中文字幕亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 成人三级黄色视频| 久久青草综合色| 电影成人av| 精品人妻1区二区| 两个人看的免费小视频| 最近最新中文字幕大全免费视频| 999久久久国产精品视频| 国产av一区在线观看免费| 丝袜美腿诱惑在线| 国产成人欧美| 免费av中文字幕在线| 俄罗斯特黄特色一大片| 久久婷婷成人综合色麻豆| 熟女少妇亚洲综合色aaa.| 一个人免费在线观看的高清视频| 日本黄色视频三级网站网址| 精品午夜福利视频在线观看一区| 成人免费观看视频高清| 黄网站色视频无遮挡免费观看| 99久久综合精品五月天人人| 美女高潮到喷水免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲成国产人片在线观看| 99在线人妻在线中文字幕| 欧美日韩av久久| 97超级碰碰碰精品色视频在线观看| 日本精品一区二区三区蜜桃| 亚洲av五月六月丁香网| 操美女的视频在线观看| 80岁老熟妇乱子伦牲交| 久久青草综合色| 日韩欧美免费精品| 一二三四在线观看免费中文在| 一级作爱视频免费观看| 久久久国产欧美日韩av| 高清在线国产一区| 午夜亚洲福利在线播放| ponron亚洲| 成人18禁高潮啪啪吃奶动态图| 91成年电影在线观看| 高清欧美精品videossex| 超色免费av| 国产亚洲精品第一综合不卡| 看免费av毛片| 99riav亚洲国产免费| 精品国产一区二区三区四区第35| 亚洲国产欧美网|