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

    軸流氣吸式排種器的設(shè)計(jì)與試驗(yàn)

    2023-12-11 21:57:58張宇乾李貴榮張成堯楊文彩劉峰張海東
    江蘇農(nóng)業(yè)科學(xué) 2023年16期
    關(guān)鍵詞:正交試驗(yàn)數(shù)值模擬

    張宇乾 李貴榮 張成堯 楊文彩 劉峰 張海東

    摘要:針對(duì)傳統(tǒng)氣吸式排種器負(fù)壓沿程損失大、結(jié)構(gòu)復(fù)雜等問題,基于軸流風(fēng)機(jī)原理及流體力學(xué)相關(guān)理論設(shè)計(jì)了一種不需要配備真空泵及配套管道設(shè)施的軸流氣吸式排種器,簡(jiǎn)化了排種器結(jié)構(gòu)。建立排種器型孔處種子的力學(xué)模型;以大葉香菜種子為排種對(duì)象,基于三維建模軟件、CFD相關(guān)軟件進(jìn)行結(jié)構(gòu)設(shè)計(jì)與數(shù)值模擬,確定了排種器關(guān)鍵部件的結(jié)構(gòu)參數(shù);搭建排種器試驗(yàn)臺(tái)架,以扇葉轉(zhuǎn)速、排種盤轉(zhuǎn)速、型孔數(shù)量為影響因子,以種子吸附率為性能指標(biāo),進(jìn)行三因子三水平二次回歸正交試驗(yàn),利用Design-Expert軟件建立各影響因子與性能指標(biāo)之間的回歸模型,分析了各個(gè)因子對(duì)排種性能的影響規(guī)律。采用主目標(biāo)函數(shù)法進(jìn)行了優(yōu)化,確定了最佳參數(shù)組合為扇葉轉(zhuǎn)速1 914 r/min、排種盤轉(zhuǎn)速 11 r/min、吸種型孔數(shù)量18個(gè),此時(shí)理論種子吸附率為92.81%。以此參數(shù)組合在試驗(yàn)臺(tái)上進(jìn)行試驗(yàn)驗(yàn)證,結(jié)果顯示種子吸附率平均值為92.22%,與理論結(jié)果基本相符。該排種器可以滿足蔬菜精密播種對(duì)排種器的性能要求。

    關(guān)鍵詞:軸流氣吸;排種器;計(jì)算流體動(dòng)力學(xué)(CFD);數(shù)值模擬;正交試驗(yàn)

    中圖分類號(hào):S223.2+3文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1002-1302(2023)16-0179-08

    收稿日期:2022-11-28

    基金項(xiàng)目:云南省重大科技專項(xiàng)(編號(hào):202102AE090042-06、2018ZC001-4)。

    作者簡(jiǎn)介:張宇乾(1995—),男,浙江紹興人,碩士研究生,研究方向?yàn)檗r(nóng)業(yè)裝備設(shè)計(jì)制造。E-mail:598220645@qq.com。

    通信作者:張海東,副教授,研究方向?yàn)檗r(nóng)業(yè)裝備設(shè)計(jì)制造、農(nóng)產(chǎn)品品質(zhì)無損檢測(cè)。E-mail:zhd_74@126.com。

    精密排種器是決定播種機(jī)播種性能的核心部件[1,直接影響播種機(jī)的能耗、播種效率、播種均勻度、傷種率等性能指標(biāo),最終影響作物的產(chǎn)量、質(zhì)量和成本2。精密排種器主要有機(jī)械式和氣力式兩大類3。機(jī)械式排種器利用重力等實(shí)現(xiàn)充種,并通過機(jī)械裝置完成吸種、排種、清種過程。常見的機(jī)械式排種器有型孔盤式、窩眼輪式、型孔帶式、指夾式等[4。氣力式排種器通常由動(dòng)力輸出軸或者液壓系統(tǒng)帶動(dòng)風(fēng)機(jī)工作,由風(fēng)機(jī)產(chǎn)生真空吸力(壓力)對(duì)種子進(jìn)行控制,使種子完成預(yù)期運(yùn)動(dòng)。氣力式排種器按照工作原理可分為氣吸式、氣壓式、氣吹式[5。其中,氣吸式排種器是精密播種機(jī)上常用的排種器之一,這種排種器具有適用范圍廣、作業(yè)效率高、播量精準(zhǔn)、不損傷種子等優(yōu)點(diǎn)6-8,但也存在一些不足之處。風(fēng)機(jī)產(chǎn)生的負(fù)壓必須通過專門的管道送至排種器負(fù)壓室,由于管道較長,易形成壓力損失,帶來不必要的能耗。且排種器設(shè)置正壓區(qū)和負(fù)壓區(qū),這2個(gè)區(qū)域由吸種盤隔開,整體結(jié)構(gòu)較復(fù)雜,裝配、維修保養(yǎng)均較困難。

    本研究擬探索設(shè)計(jì)一種無需負(fù)壓輸送管道,且結(jié)構(gòu)較為簡(jiǎn)單的氣吸式排種器,以期解決現(xiàn)有氣吸式排種器易負(fù)壓沿程損失大、結(jié)構(gòu)復(fù)雜等不足之處,也為氣吸式排種器的設(shè)計(jì)提供理論參考。

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

    區(qū)別于傳統(tǒng)的氣吸式排種器,本研究設(shè)計(jì)的軸流氣吸式排種器不需要配備真空泵及配套管道設(shè)施,也無需在排種器內(nèi)加工出正壓區(qū)、負(fù)壓區(qū)及各種復(fù)雜的孔道等,大大簡(jiǎn)化了排種器結(jié)構(gòu)。該軸流氣吸式排種器主要由傳動(dòng)軸、軸承固定盤、扇葉、氣室、排種盤、種箱、擋壓板等組成。

    排種器結(jié)構(gòu)設(shè)計(jì)如圖1所示。擋圈與氣室通過螺釘聯(lián)接,其中心轂孔用來安置軸承,傳動(dòng)軸由擋圈和聯(lián)軸器支承,其主要作用是將電機(jī)輸出軸的動(dòng)力與運(yùn)動(dòng)傳遞給扇葉;排種盤由聯(lián)軸器連接減速電機(jī)驅(qū)動(dòng)。排種盤一側(cè)轂孔安裝軸承,傳動(dòng)軸與該軸承內(nèi)圈緊聯(lián)接,從而實(shí)現(xiàn)扇葉與排種盤的差速旋轉(zhuǎn);擋壓板固定在負(fù)壓腔氣室,與排種盤最底部相接觸。工作時(shí)扇葉高速旋轉(zhuǎn),將氣體沿軸向從氣室一側(cè)排出,氣室另一側(cè)將產(chǎn)生負(fù)壓,即排種盤型孔處將擁有吸力,此時(shí)排種盤旋轉(zhuǎn),完成吸種、攜種。當(dāng)排種盤上的吸種孔攜帶種子轉(zhuǎn)到最低處時(shí),受擋壓板作用,吸種型孔失去負(fù)壓,種子在重力作用下完成排種過程。

    2 型孔處種子力學(xué)模型

    將種子視為近似球體,忽略播種過程中的振動(dòng),排種器工作時(shí),種子在型孔上受到吸力Fi、種子所受支持力N等。以種子所受摩擦力為x軸正向、所受離心力為y軸正向、所受吸附力為z軸正向,對(duì)排種型孔處種子進(jìn)行分析(圖2)。其中Q為吸種作用點(diǎn)。

    如圖2所示,建立種子平衡方程如公式(1)所示。

    式中:G為種子所受重力,N;N為種子所受支持力,N;J為種子所受離心力,N;Ff為種子在排種盤上所受的摩擦力,N;ω為排種盤旋轉(zhuǎn)角速度,rad/s;θ為種子所受離心力與重力之間的夾角,°;Fr為G、J、Ff的合力,N;NX為支持力x軸分力,N;NY為支持力y軸分力,N;Nz為支持力z軸分力,N;Fi為種子所受負(fù)壓吸力,N;α為N與Nxy夾角,°。

    由式(1)可得到種子在xy平面所受的支持力為:

    所受G、J、Ff的合力Fr為:

    在吸種過程中,吸種型孔處的種子受力平衡,則:

    在吸種型孔處的負(fù)壓大小為:

    式中:P為負(fù)壓大小,Pa;A為型孔截面面積,mm2。

    從式(5)中可初步看出吸種型孔負(fù)壓與排種盤轉(zhuǎn)速n、排種盤直徑D、型孔直徑d、種子所受重力G(與質(zhì)量m有關(guān))、種子與排種盤間的摩擦系數(shù)有關(guān),因此排種盤設(shè)計(jì)過程中需考慮以上因素。

    3 關(guān)鍵部件設(shè)計(jì)與仿真

    本研究提出了一種有別于傳統(tǒng)氣吸式排種器的結(jié)構(gòu),而排種盤、氣室、扇葉是本設(shè)計(jì)的關(guān)鍵部位,其尺寸參數(shù)直接影響排種性能,因此基于流體力學(xué)相關(guān)理論、三維建模軟件、CFD相關(guān)軟件進(jìn)行結(jié)構(gòu)設(shè)計(jì)和仿真試驗(yàn)。

    3.1 排種盤

    3.1.1 排種盤結(jié)構(gòu)尺寸 排種盤的整體結(jié)構(gòu)尺寸會(huì)直接影響氣室、扇葉的結(jié)構(gòu)尺寸,從而直接影響到整個(gè)排種器的尺寸,所以排種盤的尺寸參數(shù)是決定排種器整體結(jié)構(gòu)尺寸的關(guān)鍵因子之一[9。通過相關(guān)資料可知,現(xiàn)有垂直圓盤氣吸式排種器的種盤直徑一般為140~260 mm[10,吸種型孔一般設(shè)置在距離排種盤邊緣15~20 mm處,綜合考慮,取排種盤直徑為200 mm,型孔中心與種盤中心距離為85 mm。

    3.1.2 型孔設(shè)計(jì)與仿真分析 (1)型孔大小確定。吸種型孔直徑大小取決于種子的幾何特性,以大葉香菜種子為研究對(duì)象,隨機(jī)選取500粒顆粒飽滿的種子,測(cè)得平均粒徑為3.14 mm。吸種型孔入口直徑d與種子粒徑的關(guān)系[11為:

    d=(0.6~0.7)dk。(6)

    式中:dk為種子粒徑,mm。

    根據(jù)式(6)可得吸種型孔入口直徑d取值為1.88~2.20 mm,根據(jù)上述型孔處種子受力分析可知,入口直徑越大則吸附種子所需的負(fù)壓就越小,因此取d=2.20 mm。

    (2)型孔形狀確定。吸種型孔結(jié)構(gòu)設(shè)計(jì)對(duì)于排種盤的吸種、投種效能至關(guān)重要。根據(jù)流量連續(xù)性方程、伯努利方程設(shè)計(jì)3種變截面的吸種型孔,分別為直筒型、倒角型、錐角型;其中直筒型作為對(duì)照組,3種型孔的入口直徑和軸向長度相等。對(duì)幾種型孔進(jìn)行數(shù)值模擬。使用ANSYS MESHING對(duì)其進(jìn)行網(wǎng)格劃分,將網(wǎng)格數(shù)據(jù)平移至FLUENT進(jìn)行仿真分析。在FLUENT中設(shè)置湍流模型為realizable k-模型。根據(jù)湍流模型的壁面函數(shù)要求,將y+控制在30~300之間[12。參考相關(guān)文獻(xiàn)并根據(jù)前期研究,設(shè)置型孔出口的氣流速度為10 m/s;設(shè)置算法為coupled;控制方程采用有限體積法進(jìn)行離散化;設(shè)置迭代步數(shù)為300步;隨著迭代步數(shù)的增加,殘差逐漸降低至0.000 1,且監(jiān)測(cè)的物理量波動(dòng)越來越小并趨于穩(wěn)定,表示此時(shí)計(jì)算已達(dá)到收斂[13-14。創(chuàng)建中間截面,并顯示結(jié)果,數(shù)值模擬結(jié)果如表1所示。

    從表1可以看出,在型孔入口直徑與軸向尺寸相同情況下,直筒型型孔中氣流速度最大值為11.46 m/s,且由于氣流在穿過型孔時(shí),型孔各橫截面的面積相同,故氣流無加速過程,產(chǎn)生的負(fù)壓也較小,最大負(fù)壓為76.42 Pa;倒角型、錐角型入口速度明顯增加,且最大值分別達(dá)到了40.44、42.88 m/s;型孔入口氣流速度明顯大于型孔出口,導(dǎo)致型孔入口處負(fù)壓急劇升高,因此產(chǎn)生負(fù)壓效果好,其中倒角型入口負(fù)壓又遠(yuǎn)大于錐角型。綜上分析,根據(jù)流量連續(xù)性方程、伯努利方程設(shè)計(jì)的型孔結(jié)構(gòu)產(chǎn)生負(fù)壓效果良好,且倒角型負(fù)壓效果要優(yōu)于錐角型,因此選用倒角型作為排種盤吸種型孔。

    (3)型孔數(shù)量確定 根據(jù)農(nóng)藝要求香菜種子播種株距為3~5 cm,設(shè)計(jì)型孔數(shù)量時(shí),排種器型孔數(shù)、排種盤轉(zhuǎn)速、播種機(jī)作業(yè)速度及株距之間的關(guān)系為:

    式中:Z為排種器型孔數(shù);v為播種機(jī)作業(yè)速度,km/h;n為排種盤轉(zhuǎn)速,r/min;l為播種株距,m。

    根據(jù)式(7)可知,型孔數(shù)量與排種盤轉(zhuǎn)速成反比,轉(zhuǎn)速過高會(huì)導(dǎo)致型孔數(shù)量過少;而研究表明,排種器轉(zhuǎn)速n=10~30 r/min時(shí),排種效果較好[15,因此取排種盤轉(zhuǎn)速n=10 r/min。參考相似株距播種機(jī)作業(yè)速度v=0.45 km/h[16,則根據(jù)公式(7)可得排種型孔數(shù)量Z的取值為15~25,暫取Z=20。

    3.2 氣室

    3.2.1 氣室結(jié)構(gòu)尺寸 本研究所設(shè)計(jì)氣室從結(jié)構(gòu)上分為扇葉腔室、負(fù)壓腔室這2個(gè)腔室。其中扇葉腔室長90 mm、外徑210 mm、壁厚3 mm;負(fù)壓腔室長60 mm、外徑196 mm、壁厚3 mm。扇葉腔室內(nèi)安裝扇葉,通過扇葉的高速旋轉(zhuǎn)使氣流沿軸向流動(dòng),并在負(fù)壓腔室內(nèi)產(chǎn)生負(fù)壓,從而使排種盤型孔處擁有吸種負(fù)壓;負(fù)壓腔室與排種盤、擋壓板相連,完成吸種、攜種、排種過程。其結(jié)構(gòu)如圖3-a所示,該氣室為喇叭型結(jié)構(gòu),由直筒型結(jié)構(gòu)(圖3-b)通過仿真分析后優(yōu)化得到。

    由圖3可知,喇叭型氣室和直筒型氣室的主要區(qū)別在于負(fù)壓腔室的變化。相對(duì)于直筒型氣室而言,喇叭型氣室的負(fù)壓腔室與扇葉腔室之間的截面變化較大,此設(shè)計(jì)可大幅降低扇葉腔室對(duì)負(fù)壓腔室的氣流擾動(dòng)作用。

    3.2.2 氣室流場(chǎng)仿真分析 為了驗(yàn)證氣室結(jié)構(gòu)設(shè)計(jì)的可行性與可靠性,使用Space Claim建立2種氣室的仿真模型,并導(dǎo)入FLUENT中進(jìn)行數(shù)值模擬,設(shè)置扇葉轉(zhuǎn)速為1 800 r/min,計(jì)算結(jié)束后在CFD-Post中觀察其中間截面壓力和速度,以型孔處的速度和負(fù)壓值作為評(píng)價(jià)標(biāo)準(zhǔn),調(diào)整壓強(qiáng)范圍。

    直筒型結(jié)構(gòu)數(shù)值模擬結(jié)果如圖4所示。圖4壓強(qiáng)和速度云圖顯示負(fù)壓腔室邊緣顏色為橘黃色,中部顏色為藍(lán)色。表示此時(shí)負(fù)壓腔室壁面處壓強(qiáng)為正。這是由于左、右兩部分內(nèi)徑相近,沿軸向無明顯截面變化,左右兩室相接相通,因此型孔處受氣流擾動(dòng)導(dǎo)致負(fù)壓效果差,甚至無法獲得負(fù)壓。

    喇叭型數(shù)值模擬結(jié)果如圖5所示。圖5壓強(qiáng)云圖顯示負(fù)壓腔室整體顏色為淡藍(lán)色,速度云圖顯示負(fù)壓腔室整體顏色為藍(lán)色,而型孔處速度為淡藍(lán)色,表示此時(shí)負(fù)壓腔室壁面處壓強(qiáng)為負(fù)壓且型孔入口處氣流速度增加。這是由于喇叭型氣室左、右兩部分沿軸向截面產(chǎn)生變化,因此右側(cè)氣室基本不受氣流擾動(dòng),而型孔由于自身結(jié)構(gòu)特點(diǎn)導(dǎo)致入口處流速增加,使得負(fù)壓在不損失的前提下進(jìn)一步提升。

    2種結(jié)構(gòu)的型孔入口處速度和負(fù)壓均值如表2所示。由表2可知,喇叭型結(jié)構(gòu)的型孔入口負(fù)壓為120.57 Pa,直筒型結(jié)構(gòu)為16.07 Pa;喇叭型結(jié)構(gòu)的型孔入口速度為14.02 m/s,直筒型結(jié)構(gòu)為4.83 m/s。因此喇叭型結(jié)構(gòu)產(chǎn)生負(fù)壓的效果遠(yuǎn)優(yōu)于直筒型,所以將喇叭型結(jié)構(gòu)作為排種器氣室結(jié)構(gòu)。

    3.3 扇葉

    3.3.1 扇葉結(jié)構(gòu)尺寸 氣室扇葉腔室內(nèi)徑為 207 mm,因此扇葉直徑需小于此尺寸。對(duì)市場(chǎng)上常見的幾種軸流風(fēng)機(jī)進(jìn)行風(fēng)量和靜壓對(duì)比,最終以德喜機(jī)電HZ200型軸流風(fēng)機(jī)的扇葉葉型作為參考進(jìn)行設(shè)計(jì)。設(shè)計(jì)扇葉直徑為187 mm,葉片和壁面厚度2.5 mm,扇葉結(jié)構(gòu)示意如圖6所示。扇葉輪轂形狀為半球形,此設(shè)計(jì)可以提高葉片根部的強(qiáng)度。

    3.3.2 扇葉流固耦合仿真分析 為了確定扇葉轉(zhuǎn)速范圍,對(duì)扇葉進(jìn)行仿真試驗(yàn)。以型孔入口平均壓強(qiáng)和扇葉變形量作為指標(biāo)對(duì)扇葉進(jìn)行流固耦合仿真分析。首先建立控制方程組,采用N-S方程、realizable k-模型、coupled算法進(jìn)行計(jì)算。連續(xù)性方程:

    N-S方程:

    realizable k-ε模型湍動(dòng)能和耗散率方程分別如下:

    單向流固耦合分析方法[30如下:(1)Fluent流場(chǎng)分析;(2)將流場(chǎng)分析結(jié)果傳遞給結(jié)構(gòu)場(chǎng);(3)耦合分析并輸出結(jié)果。流固耦合分析流程如圖7所示。

    使用MESHING對(duì)模型進(jìn)行命名設(shè)置和網(wǎng)格劃分,然后導(dǎo)入FLUENT中,使用多重參考系(MRF)方法進(jìn)行流場(chǎng)模擬,模擬結(jié)束后,計(jì)算型孔入口平均負(fù)壓并將流場(chǎng)數(shù)據(jù)傳輸至Static Structure。在Engineering Data中添加Plastic(ABS)材料。在Static Structure中進(jìn)行材料設(shè)置,添加重力和設(shè)置固定支撐等操作并進(jìn)行扇葉網(wǎng)格劃分。

    對(duì)扇葉轉(zhuǎn)速進(jìn)行流固耦合仿真試驗(yàn),試驗(yàn)結(jié)果如圖8所示。

    由圖9可見,型孔入口平均負(fù)壓和扇葉變形量隨著扇葉轉(zhuǎn)速增加而增大,呈正相關(guān)性。其中當(dāng)扇葉轉(zhuǎn)速為600 r/min時(shí),型孔入口負(fù)壓最小,為12.75 Pa;當(dāng)扇葉轉(zhuǎn)速為2 700 r/min時(shí),型孔入口負(fù)壓最大,為259.38 Pa。其中當(dāng)扇葉轉(zhuǎn)速為2 700 r/min時(shí),扇葉變形情況如圖10所示。

    由圖10可知,扇葉中部顏色為藍(lán)色,表示此處變形最?。簧热~從根部到邊緣顏色逐漸從黃綠變?yōu)槌燃t色,表示變形量逐漸增加,到達(dá)邊緣時(shí)最大變形約為0.10 mm。由此可知扇葉邊緣為扇葉結(jié)構(gòu)的薄弱之處,若轉(zhuǎn)速太高,則扇葉邊緣將會(huì)發(fā)生變形,從而擾動(dòng)氣流導(dǎo)致排風(fēng)減少,則負(fù)壓腔室型孔處產(chǎn)生負(fù)壓效果將變差,因此扇葉轉(zhuǎn)速不宜過高。

    根據(jù)上述型孔處種子受力分析并忽略摩擦力,初步估計(jì)吸種負(fù)壓應(yīng)>26 Pa。當(dāng)扇葉轉(zhuǎn)速過低時(shí),產(chǎn)生的負(fù)壓難以吸附種子;而當(dāng)轉(zhuǎn)速過高時(shí),產(chǎn)生的振動(dòng)噪聲、功耗以及葉片變形都會(huì)較大。綜合考慮選擇扇葉轉(zhuǎn)速范圍為1 200~2 100 r/min。

    4 排種器排種性能臺(tái)架試驗(yàn)

    4.1 試驗(yàn)材料及儀器設(shè)備

    本試驗(yàn)地點(diǎn)為云南省昆明市云南農(nóng)業(yè)大學(xué)機(jī)電工程實(shí)驗(yàn)教學(xué)中心,采用大葉香菜種子作為排種器排種對(duì)象。

    本試驗(yàn)的排種器采用3D打印技術(shù)加工,使用材料為ABS白色樹脂;試驗(yàn)臺(tái)采用鋁合金型材進(jìn)行搭建,選用普斯調(diào)速電機(jī)和感應(yīng)調(diào)速電機(jī)(4IK25GN-C)分別驅(qū)動(dòng)扇葉和排種盤,所搭建的試驗(yàn)臺(tái)如圖11所示。

    4.2 試驗(yàn)設(shè)計(jì)與方法

    根據(jù)前期研究確立扇葉轉(zhuǎn)速范圍為1 200~2 100 r/min;排種盤轉(zhuǎn)速范圍為10~30 r/min;型孔數(shù)量為15~25個(gè)。以種子吸附率Y1為響應(yīng)值,其中Y1=吸附種子孔數(shù)/總孔數(shù);以扇葉轉(zhuǎn)速X1、排種盤轉(zhuǎn)速X2、型孔數(shù)量X3為試驗(yàn)因子。每組試驗(yàn)重復(fù)3次,試驗(yàn)結(jié)果取3次平均值。采用二次回歸正交設(shè)計(jì)試驗(yàn),試驗(yàn)因子水平編碼如表3所示。

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

    試驗(yàn)結(jié)果如表4所示,模型顯著性檢驗(yàn)如表5所示。

    種子之間的機(jī)械物理特性(如直徑、質(zhì)量)存在差異,且存在種箱高度改變、振動(dòng)情況不同等因素,因此同一因素水平下種子吸附率存在波動(dòng)。

    對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合分析,以種子吸附率Y為響應(yīng)函數(shù),以各影響因子的水平編碼值為自變量,建立回歸數(shù)學(xué)模型。通過對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合,得到各因子水平對(duì)排種器種子吸附率影響的數(shù)學(xué)回歸模型,如公式(12)所示。

    Y1=83.21+10.05X1-5.48X2-1.99X3+0.73X1X2+0.73X1X3-0.011X2X3-4.32X21-1.61X22-1.04X23。(12)

    根據(jù)表5可知,氣吸式排種器的種子吸附率回歸模型的P值<0.000 1,故影響極顯著;失擬項(xiàng)的P值為0.125 3(P>0.05),故影響不顯著。因此可以確定,在一定的參數(shù)范圍內(nèi),本試驗(yàn)種子吸附率的回歸模型與實(shí)際的播種情況擬合度較高。由扇葉轉(zhuǎn)速、排種盤轉(zhuǎn)速、型孔數(shù)量的P值可判斷,扇葉轉(zhuǎn)速、排種盤轉(zhuǎn)速對(duì)種子吸附率都有極顯著影響,型孔數(shù)量對(duì)種子吸附率無顯著影響。進(jìn)一步可知回歸項(xiàng)X1、X2、X21的P值均<0.01,影響極顯著;回歸項(xiàng)X3、X1X2、X1X3、X2X3、X22、X23的P值均>0.05,影響不顯著。在剔除不顯著因子回歸項(xiàng)后,種子吸附率的回歸模型可以表示為:

    Y=83.21+10.05X1-5.48X2-4.32X21。(13)

    為更加直觀地分析各影響因子與排種性能之間的關(guān)系,根據(jù)回歸模型分析結(jié)果,使用Design-Expert10對(duì)數(shù)據(jù)進(jìn)行處理,可以得到扇葉轉(zhuǎn)速X1、排種盤轉(zhuǎn)速X2、型孔數(shù)量X3、對(duì)排種器排種過程種子吸附率的影響。任意選取其中1個(gè)因子固定其水平,根據(jù)響應(yīng)面圖及等值線圖,分別分析另外2個(gè)因子對(duì)種子吸附率的影響,響應(yīng)面和等值線如圖12所示。

    圖12-a是型孔數(shù)量為20個(gè)時(shí),扇葉轉(zhuǎn)速和排種盤轉(zhuǎn)速對(duì)排種器種子吸附率的響應(yīng)面圖。可知,當(dāng)扇葉轉(zhuǎn)速一定時(shí),種子吸附率隨排種盤轉(zhuǎn)速增大而減小,是因?yàn)樾涂淄ㄟ^充種區(qū)域的時(shí)間減少,種子被吸附的機(jī)會(huì)也隨之減小。當(dāng)排種盤轉(zhuǎn)速一定時(shí),種子吸附率隨扇葉轉(zhuǎn)速增大而增大,是因?yàn)樾涂兹肟谔幍呢?fù)壓增加,吸種能力得到提高。

    圖12-b是排種盤轉(zhuǎn)速為20 r/min時(shí),扇葉轉(zhuǎn)速和型孔數(shù)量對(duì)排種器種子吸附率的響應(yīng)面圖??芍?dāng)扇葉轉(zhuǎn)速一定時(shí),隨著型孔數(shù)量的增多,型孔入口處的負(fù)壓變化不大,因此種子吸附率變化很小,說明型孔數(shù)量對(duì)排種器種子吸附率影響不顯著。當(dāng)型孔數(shù)量一定時(shí),隨著扇葉轉(zhuǎn)速增大,型孔入口處負(fù)壓增加,吸種能力得到提高,因此種子吸附率增大。

    圖12-c是扇葉轉(zhuǎn)速為1 650 r/min時(shí),排種盤轉(zhuǎn)速和型孔數(shù)量對(duì)排種器種子吸附率的響應(yīng)面圖??芍?,當(dāng)排種盤轉(zhuǎn)速一定時(shí),種子吸附率隨型孔數(shù)量增加而變化很小,是因?yàn)樾涂兹肟谔幍呢?fù)壓變化不大。當(dāng)型孔數(shù)量一定時(shí),種子吸附率隨排種盤轉(zhuǎn)速的增大而減小。

    4.4 參數(shù)優(yōu)化與驗(yàn)證

    通過上述分析,為了得到吸種最佳參數(shù)組合,采用多重響應(yīng)方法中的主目標(biāo)函數(shù)對(duì)影響因子扇葉轉(zhuǎn)速、排種盤轉(zhuǎn)速、型孔數(shù)量進(jìn)行優(yōu)化,以種子吸附率最大值為性能指標(biāo)函數(shù),根據(jù)上述響應(yīng)面分析和回歸方程,建立其數(shù)學(xué)模型:

    運(yùn)用Design-Expert 10軟件中進(jìn)行優(yōu)化求解,得出影響因子的最佳參數(shù)組合:扇葉轉(zhuǎn)速為 1 914 r/min,排種盤轉(zhuǎn)速為11 r/min,型孔數(shù)量為18個(gè),目標(biāo)函數(shù)預(yù)測(cè)值為:種子吸附率92.81%。為了驗(yàn)證最佳參數(shù)組合,在上述最優(yōu)參數(shù)組合試驗(yàn)條件下,進(jìn)行5次重復(fù)試驗(yàn),試驗(yàn)結(jié)果顯示種子吸附率平均值為92.22%,可見優(yōu)化結(jié)果與試驗(yàn)結(jié)果基本一致。

    5 結(jié)論

    本研究針對(duì)傳統(tǒng)氣吸式排種器能耗大、結(jié)構(gòu)復(fù)雜的問題,依據(jù)軸流風(fēng)機(jī)原理、流量連續(xù)性方程、伯努利方程設(shè)計(jì)了軸流氣吸式排種器,對(duì)其進(jìn)行了數(shù)值模擬和臺(tái)架試驗(yàn),現(xiàn)得出主要結(jié)論如下:

    (1)對(duì)排種器吸種過程中種子進(jìn)行受力分析,結(jié)果表明,吸種型孔負(fù)壓與排種盤轉(zhuǎn)速、排種盤直徑、型孔直徑、種子質(zhì)量、種子與排種盤間的摩擦系數(shù)有關(guān)。

    (2)研制了一種軸流氣吸式排種器:針對(duì)排種盤結(jié)構(gòu)尺寸、型孔參數(shù)、氣室與扇葉結(jié)構(gòu)進(jìn)行理論分析,并結(jié)合三維建模軟件、數(shù)值模擬軟件以及排種器吸種過程受力分析,確定了各部件結(jié)構(gòu)的形狀和參數(shù)。

    (3)對(duì)所設(shè)計(jì)的排種器進(jìn)行加工,搭建排種器試驗(yàn)臺(tái)。以扇葉轉(zhuǎn)速、排種盤轉(zhuǎn)速、吸種型孔數(shù)量作為影響因子,以種子吸附率作為評(píng)價(jià)指標(biāo),對(duì)排種器開展3因子5水平的排種性能試驗(yàn)設(shè)計(jì),得到回歸模型。結(jié)果得到理論最優(yōu)組合為扇葉轉(zhuǎn)速 1 914 r/min、排種盤轉(zhuǎn)速11 r/min、吸種型孔數(shù)量18個(gè),并以此參數(shù)組合在試驗(yàn)臺(tái)上進(jìn)行試驗(yàn)驗(yàn)證。試驗(yàn)結(jié)果與理論結(jié)果基本相符,期待可為后續(xù)排種器設(shè)計(jì)提供參考依據(jù)。

    參考文獻(xiàn):

    [1]張 昆,衣淑娟,劉海軍,等. 氣吸滾筒式玉米植質(zhì)缽盤精密播種裝置優(yōu)化設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2021,52(4):62-69.

    [2]張恒榜,趙 健,胡廣銳,等. 窩眼輪式小麥排種器參數(shù)優(yōu)化試驗(yàn)研究[J]. 農(nóng)機(jī)化研究,2020,42(9):139-144.

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

    [4]袁文勝,吳崇友,金誠謙. 異形孔窩眼輪式油菜排種器設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(5):72-75.

    [5]高麗萍,施彬彬,廖慶喜,等. 正負(fù)氣壓組合油菜精量排種器錐孔盤排種性能[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(6):22-33.

    [6]李玉環(huán),楊 麗,張東興,等. 氣吸式玉米高速精量排種器投種性能分析與結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(8):1-11.

    [7]孫 晗,郭廣川,唐 帥,等. 蔬菜育苗播種流水線壓穴滾筒裝置改進(jìn)與控制系統(tǒng)設(shè)計(jì)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(22):41-48.

    [8]張 靜.氣吸滾筒式的小粒扁平種子精量播種機(jī)理研究[D]. 廣州:華南農(nóng)業(yè)大學(xué),2017.

    [9]倪向東,徐國杰,王 琦,等. 氣吸滾筒陣列式棉花精密排種器設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(12):58-67.

    [10]陳海濤,王洪飛,王業(yè)成,等. 三葉式自動(dòng)清換種大豆育種氣吸排種器設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2020,51(12):75-85.

    [11]Cao X L,Ma X,Li H W,et al. Design and experiments of pneumatic roller type precision seed-metering device for rapeseed plug seedlings[J]. Transactions of the Chinese Society of Agricultural Engineering,2021,37(2):51-60.

    [12]Zhu D,Xiao R F,Yao Z F,et al. Optimization design for reducing the axial force of a vaned mixed-flow pump[J]. Engineering Applications of Computational Fluid Mechanics,2020,14(1):882-896.

    [13]Han D D,Zhang D X,Yang L,et al. Effects of seeds morphology on pressure distribution in the cone hole of air-blowing seed-metering device based on CFD[C]//Spokane,Washington.St.Joseph,MI:American Society of Agricultural and Biological Engineers,2017.

    [14]Liu L,Du Y F,Wang L. Design and experiment of pneumatic roller type seed-metering device with small particle size[C]. 2020 ASABE Annual International Virtual Meeting.

    [15]楊文彩,闞成龍,張效偉,等. 氣吸式三七育苗播種精密排種器設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2021,52(6):95-105.

    [16]尹文慶,趙 璐,李 驊,等. 氣力槽輪組合式蔬菜精密排種器吸嘴型孔設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(4):68-76,136.

    猜你喜歡
    正交試驗(yàn)數(shù)值模擬
    傾斜式子棉清理機(jī)工作質(zhì)量影響因素分析
    復(fù)方萬年青膠囊神經(jīng)保護(hù)作用成分的水提取工藝研究
    基于MADYMO的航空座椅約束系統(tǒng)優(yōu)化設(shè)計(jì)
    科技視界(2016年18期)2016-11-03 21:44:44
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    祁山藥醇提物提取工藝研究
    科技視界(2016年20期)2016-09-29 11:36:29
    一種基于液壓緩沖的減震管卡設(shè)計(jì)與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    少妇的逼水好多| www.熟女人妻精品国产 | 18禁观看日本| 亚洲久久久国产精品| 国产成人精品久久久久久| 精品久久久久久电影网| 卡戴珊不雅视频在线播放| 国产日韩欧美视频二区| kizo精华| 久久久久久久久久久久大奶| 一级a做视频免费观看| 一本大道久久a久久精品| 日韩制服丝袜自拍偷拍| 国产欧美亚洲国产| 欧美bdsm另类| 好男人视频免费观看在线| 黄色 视频免费看| 男人添女人高潮全过程视频| 边亲边吃奶的免费视频| 国产精品免费大片| 少妇人妻精品综合一区二区| 国产午夜精品一二区理论片| 18在线观看网站| 国产男女内射视频| xxxhd国产人妻xxx| 最近的中文字幕免费完整| 丝袜喷水一区| 你懂的网址亚洲精品在线观看| 国产精品秋霞免费鲁丝片| 久久99热这里只频精品6学生| 久久久精品区二区三区| 看免费成人av毛片| 国产极品粉嫩免费观看在线| 久久热在线av| 欧美丝袜亚洲另类| 欧美激情国产日韩精品一区| 午夜av观看不卡| 久久午夜综合久久蜜桃| 亚洲精品中文字幕在线视频| 欧美激情国产日韩精品一区| 精品国产一区二区久久| 国产免费视频播放在线视频| 欧美激情 高清一区二区三区| 国产在线免费精品| 哪个播放器可以免费观看大片| 国产视频首页在线观看| 99久久精品国产国产毛片| 亚洲国产精品国产精品| 成人亚洲精品一区在线观看| freevideosex欧美| 久久久久人妻精品一区果冻| a级毛色黄片| 国产又色又爽无遮挡免| 亚洲天堂av无毛| 天美传媒精品一区二区| 乱人伦中国视频| 色94色欧美一区二区| 又大又黄又爽视频免费| 午夜福利影视在线免费观看| 成年女人在线观看亚洲视频| 国产精品 国内视频| 一级毛片电影观看| 国产精品久久久久久久久免| 岛国毛片在线播放| 色吧在线观看| 国产黄色视频一区二区在线观看| 国产不卡av网站在线观看| 国产精品 国内视频| 国产激情久久老熟女| 一本大道久久a久久精品| 国产色爽女视频免费观看| 欧美日韩一区二区视频在线观看视频在线| 国产亚洲欧美精品永久| 99热6这里只有精品| 午夜影院在线不卡| 满18在线观看网站| 国产精品三级大全| 国产精品久久久久久精品古装| kizo精华| 精品久久久久久电影网| 春色校园在线视频观看| 午夜视频国产福利| 一区二区三区精品91| 亚洲欧洲日产国产| 99热6这里只有精品| kizo精华| 日韩中字成人| 最近中文字幕高清免费大全6| 午夜精品国产一区二区电影| 99精国产麻豆久久婷婷| 亚洲精品色激情综合| 夜夜骑夜夜射夜夜干| 一二三四在线观看免费中文在 | 22中文网久久字幕| 免费不卡的大黄色大毛片视频在线观看| 欧美 亚洲 国产 日韩一| 久久国产精品男人的天堂亚洲 | 中文字幕亚洲精品专区| 肉色欧美久久久久久久蜜桃| 在线观看三级黄色| 欧美日韩av久久| 久久毛片免费看一区二区三区| 日韩av不卡免费在线播放| 国产视频首页在线观看| 99久久人妻综合| 国产有黄有色有爽视频| 精品一区在线观看国产| 三上悠亚av全集在线观看| 老司机影院毛片| 日韩大片免费观看网站| 一级黄片播放器| 国产成人一区二区在线| 日韩精品有码人妻一区| 最后的刺客免费高清国语| 国产日韩欧美在线精品| 国产黄色免费在线视频| 啦啦啦中文免费视频观看日本| 久久久久久久国产电影| 美女大奶头黄色视频| 国产高清国产精品国产三级| 日韩制服丝袜自拍偷拍| 国产免费一级a男人的天堂| 欧美精品av麻豆av| 国产深夜福利视频在线观看| 亚洲欧美一区二区三区黑人 | 久久人人爽人人片av| 亚洲精华国产精华液的使用体验| 91国产中文字幕| 国精品久久久久久国模美| 亚洲av在线观看美女高潮| 免费观看av网站的网址| 久久影院123| 欧美 日韩 精品 国产| 女性生殖器流出的白浆| 精品久久国产蜜桃| 久久久久久久国产电影| 日韩一本色道免费dvd| 成人漫画全彩无遮挡| 欧美日韩亚洲高清精品| 国产福利在线免费观看视频| 黄片播放在线免费| 少妇精品久久久久久久| 22中文网久久字幕| 99九九在线精品视频| 男人添女人高潮全过程视频| 久久精品久久久久久噜噜老黄| 咕卡用的链子| 久久亚洲国产成人精品v| 日韩中字成人| 狂野欧美激情性bbbbbb| 国产免费现黄频在线看| 97在线视频观看| 亚洲成人一二三区av| 欧美日韩综合久久久久久| 久久久国产一区二区| 国产69精品久久久久777片| 国产白丝娇喘喷水9色精品| 国产欧美另类精品又又久久亚洲欧美| 免费看av在线观看网站| 美女脱内裤让男人舔精品视频| 国产成人精品婷婷| 丝袜在线中文字幕| 男女下面插进去视频免费观看 | 国产男女超爽视频在线观看| 亚洲精品国产色婷婷电影| 国产精品久久久久久久久免| 国产av码专区亚洲av| 国产男人的电影天堂91| 欧美 日韩 精品 国产| 国产精品久久久久久久电影| 久久热在线av| 巨乳人妻的诱惑在线观看| 亚洲精品,欧美精品| 王馨瑶露胸无遮挡在线观看| 成人午夜精彩视频在线观看| 精品国产一区二区久久| 99热这里只有是精品在线观看| 国产成人精品在线电影| 男女国产视频网站| 老女人水多毛片| 欧美精品国产亚洲| 亚洲精品国产av成人精品| 精品人妻在线不人妻| 免费观看a级毛片全部| 日韩欧美精品免费久久| 国产成人一区二区在线| 精品午夜福利在线看| 亚洲成人一二三区av| 你懂的网址亚洲精品在线观看| 久久午夜综合久久蜜桃| 欧美97在线视频| 日本av免费视频播放| 欧美人与性动交α欧美软件 | 9191精品国产免费久久| 免费看不卡的av| 中文字幕人妻熟女乱码| 九色亚洲精品在线播放| 波野结衣二区三区在线| 亚洲国产精品一区三区| 亚洲欧美色中文字幕在线| 街头女战士在线观看网站| 日韩免费高清中文字幕av| 国产午夜精品一二区理论片| 极品人妻少妇av视频| 中文乱码字字幕精品一区二区三区| 99re6热这里在线精品视频| 亚洲精品乱久久久久久| 中文字幕制服av| 最近最新中文字幕大全免费视频 | 中文字幕人妻丝袜制服| 亚洲av电影在线观看一区二区三区| 久久精品国产a三级三级三级| 成人毛片a级毛片在线播放| 校园人妻丝袜中文字幕| 亚洲人成网站在线观看播放| 美女国产视频在线观看| 三上悠亚av全集在线观看| 国产免费一级a男人的天堂| 亚洲情色 制服丝袜| 久久久久久人妻| 一级,二级,三级黄色视频| 亚洲 欧美一区二区三区| av在线app专区| 亚洲av中文av极速乱| videossex国产| 亚洲精华国产精华液的使用体验| 午夜激情久久久久久久| 97在线人人人人妻| 精品国产乱码久久久久久小说| 国产精品久久久久久精品电影小说| 国产欧美亚洲国产| 伦理电影免费视频| 黑人巨大精品欧美一区二区蜜桃 | 老熟女久久久| 日本色播在线视频| 秋霞在线观看毛片| 国产免费一级a男人的天堂| 一二三四中文在线观看免费高清| 久久婷婷青草| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品一区蜜桃| freevideosex欧美| 亚洲精品日本国产第一区| 成人免费观看视频高清| 黄片播放在线免费| 亚洲欧美成人精品一区二区| 九色成人免费人妻av| 如何舔出高潮| 91午夜精品亚洲一区二区三区| 99久久综合免费| 18禁国产床啪视频网站| 国产亚洲午夜精品一区二区久久| 少妇精品久久久久久久| av片东京热男人的天堂| 中文字幕免费在线视频6| 久久久亚洲精品成人影院| 大片免费播放器 马上看| 久久av网站| 一级片'在线观看视频| 精品99又大又爽又粗少妇毛片| 一级片免费观看大全| 国产精品久久久久久av不卡| 国产一区二区三区综合在线观看 | 久久久久久久久久成人| av天堂久久9| 国产xxxxx性猛交| 精品国产一区二区久久| 精品99又大又爽又粗少妇毛片| 亚洲精品久久久久久婷婷小说| 99久久精品国产国产毛片| 亚洲中文av在线| 国产精品一区二区在线不卡| 男女无遮挡免费网站观看| 亚洲国产av影院在线观看| 欧美亚洲日本最大视频资源| 黄片无遮挡物在线观看| 亚洲精品成人av观看孕妇| 肉色欧美久久久久久久蜜桃| 精品久久久精品久久久| 亚洲av电影在线进入| 校园人妻丝袜中文字幕| 人妻一区二区av| 欧美精品国产亚洲| 国产日韩欧美在线精品| 久久久久视频综合| 国产在线一区二区三区精| 男人操女人黄网站| 丰满饥渴人妻一区二区三| 国产福利在线免费观看视频| 欧美国产精品va在线观看不卡| 国产精品人妻久久久久久| 中国三级夫妇交换| 亚洲少妇的诱惑av| 少妇被粗大猛烈的视频| 亚洲精品456在线播放app| freevideosex欧美| 99热全是精品| 亚洲精品久久久久久婷婷小说| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99| 亚洲国产精品一区三区| 亚洲,欧美精品.| 国产激情久久老熟女| 国产一区有黄有色的免费视频| 免费在线观看黄色视频的| 在线观看www视频免费| 久久久亚洲精品成人影院| 色婷婷久久久亚洲欧美| 狠狠婷婷综合久久久久久88av| 国产1区2区3区精品| 国产毛片在线视频| 久久ye,这里只有精品| 国产免费又黄又爽又色| 777米奇影视久久| 精品亚洲乱码少妇综合久久| 日韩不卡一区二区三区视频在线| 亚洲美女搞黄在线观看| 国产欧美日韩综合在线一区二区| 久热久热在线精品观看| 一级,二级,三级黄色视频| av女优亚洲男人天堂| 免费女性裸体啪啪无遮挡网站| 美女内射精品一级片tv| 成人亚洲欧美一区二区av| 亚洲,一卡二卡三卡| 免费在线观看完整版高清| 国产亚洲午夜精品一区二区久久| 久热这里只有精品99| 啦啦啦在线观看免费高清www| av黄色大香蕉| 精品久久久久久电影网| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 欧美 日韩 精品 国产| 香蕉国产在线看| 国产日韩欧美视频二区| 成人手机av| 在线免费观看不下载黄p国产| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 免费观看av网站的网址| 五月天丁香电影| 亚洲欧美精品自产自拍| 亚洲国产色片| 天天躁夜夜躁狠狠久久av| 五月天丁香电影| 国产成人精品久久久久久| 我要看黄色一级片免费的| 波野结衣二区三区在线| 午夜老司机福利剧场| 亚洲成人av在线免费| 美女内射精品一级片tv| 国产av码专区亚洲av| 性色avwww在线观看| 亚洲av电影在线进入| 毛片一级片免费看久久久久| 欧美最新免费一区二区三区| 一区二区三区乱码不卡18| 黑丝袜美女国产一区| 亚洲成人手机| 成人亚洲欧美一区二区av| 精品卡一卡二卡四卡免费| 久久女婷五月综合色啪小说| 美国免费a级毛片| 人人妻人人澡人人爽人人夜夜| 国产av国产精品国产| 精品人妻熟女毛片av久久网站| 色视频在线一区二区三区| 91精品三级在线观看| 久久青草综合色| 亚洲精品美女久久av网站| 久久青草综合色| 在线天堂中文资源库| 国产av国产精品国产| 国产成人精品福利久久| 中文字幕av电影在线播放| 欧美激情国产日韩精品一区| 日本爱情动作片www.在线观看| 免费看光身美女| 精品少妇久久久久久888优播| 亚洲久久久国产精品| xxxhd国产人妻xxx| 大香蕉久久网| 亚洲婷婷狠狠爱综合网| 日本-黄色视频高清免费观看| 大香蕉97超碰在线| 国产女主播在线喷水免费视频网站| av女优亚洲男人天堂| 久久久久久久久久久久大奶| 麻豆精品久久久久久蜜桃| 免费人妻精品一区二区三区视频| 老司机亚洲免费影院| 满18在线观看网站| 秋霞在线观看毛片| 如日韩欧美国产精品一区二区三区| 妹子高潮喷水视频| 精品国产乱码久久久久久小说| 免费黄网站久久成人精品| 亚洲熟女精品中文字幕| 飞空精品影院首页| 搡女人真爽免费视频火全软件| 男人爽女人下面视频在线观看| 人体艺术视频欧美日本| 欧美 亚洲 国产 日韩一| 国产成人精品无人区| 欧美国产精品va在线观看不卡| 男女免费视频国产| 亚洲人与动物交配视频| 欧美日韩视频精品一区| 精品99又大又爽又粗少妇毛片| 日韩制服骚丝袜av| 少妇精品久久久久久久| 免费av中文字幕在线| 久久这里有精品视频免费| 欧美精品国产亚洲| 国产午夜精品一二区理论片| 午夜免费鲁丝| 国产高清不卡午夜福利| 午夜老司机福利剧场| 99国产综合亚洲精品| 中文字幕精品免费在线观看视频 | 国产亚洲一区二区精品| 美女主播在线视频| 精品亚洲成a人片在线观看| 最近最新中文字幕免费大全7| 又黄又爽又刺激的免费视频.| 少妇的逼好多水| 国产不卡av网站在线观看| 校园人妻丝袜中文字幕| 亚洲成人av在线免费| 国产精品熟女久久久久浪| 日韩人妻精品一区2区三区| 久久人人爽人人片av| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区蜜桃 | 免费在线观看完整版高清| 免费观看性生交大片5| 亚洲av电影在线进入| 亚洲四区av| av网站免费在线观看视频| 国产精品一国产av| 午夜福利视频精品| 在线看a的网站| www日本在线高清视频| 母亲3免费完整高清在线观看 | 97在线人人人人妻| 色哟哟·www| 考比视频在线观看| 国精品久久久久久国模美| 国产av精品麻豆| 9热在线视频观看99| 大香蕉久久成人网| 亚洲av在线观看美女高潮| 国产成人免费观看mmmm| 亚洲丝袜综合中文字幕| 免费久久久久久久精品成人欧美视频 | 欧美精品高潮呻吟av久久| 亚洲av福利一区| 午夜免费男女啪啪视频观看| av国产精品久久久久影院| 国产av精品麻豆| 人妻少妇偷人精品九色| 69精品国产乱码久久久| 丝袜喷水一区| 青春草亚洲视频在线观看| 99国产精品免费福利视频| av国产久精品久网站免费入址| 男女免费视频国产| 妹子高潮喷水视频| 99久久中文字幕三级久久日本| 在线观看免费视频网站a站| 欧美少妇被猛烈插入视频| 久久鲁丝午夜福利片| 亚洲成人手机| 纵有疾风起免费观看全集完整版| 99热全是精品| 精品久久久久久电影网| 精品一品国产午夜福利视频| 夫妻午夜视频| www.av在线官网国产| 亚洲精品av麻豆狂野| 免费av不卡在线播放| 大香蕉97超碰在线| 国产成人精品无人区| 丰满饥渴人妻一区二区三| 日韩制服骚丝袜av| 少妇的逼水好多| 久久 成人 亚洲| 国产有黄有色有爽视频| 麻豆精品久久久久久蜜桃| 一区二区av电影网| 久久婷婷青草| 日日摸夜夜添夜夜爱| 人体艺术视频欧美日本| 麻豆乱淫一区二区| 亚洲av成人精品一二三区| 国产成人精品一,二区| 两个人看的免费小视频| 国产男女超爽视频在线观看| 久久久久国产精品人妻一区二区| 免费播放大片免费观看视频在线观看| 成年女人在线观看亚洲视频| 国产av国产精品国产| 大话2 男鬼变身卡| 飞空精品影院首页| 又粗又硬又长又爽又黄的视频| 中国美白少妇内射xxxbb| 嫩草影院入口| 亚洲精品自拍成人| 国产乱人偷精品视频| 日本猛色少妇xxxxx猛交久久| 99re6热这里在线精品视频| 一本大道久久a久久精品| 超色免费av| av在线app专区| 国产精品一区二区在线不卡| 精品国产一区二区三区久久久樱花| 亚洲精品第二区| 精品一区二区免费观看| 赤兔流量卡办理| 丝袜美足系列| 天天躁夜夜躁狠狠久久av| 在线观看www视频免费| 欧美bdsm另类| 高清在线视频一区二区三区| 亚洲欧美成人综合另类久久久| 婷婷色av中文字幕| 一区二区三区精品91| 欧美成人午夜免费资源| av一本久久久久| 亚洲中文av在线| 黄色视频在线播放观看不卡| 国产极品粉嫩免费观看在线| 久热这里只有精品99| 毛片一级片免费看久久久久| 91午夜精品亚洲一区二区三区| 韩国精品一区二区三区 | 欧美激情极品国产一区二区三区 | 国产 一区精品| 国产xxxxx性猛交| 国产av国产精品国产| 边亲边吃奶的免费视频| 免费黄网站久久成人精品| 男女高潮啪啪啪动态图| 99精国产麻豆久久婷婷| 在线观看www视频免费| 精品久久蜜臀av无| 黄色配什么色好看| 九色亚洲精品在线播放| 少妇人妻久久综合中文| 97超碰精品成人国产| 精品国产国语对白av| 亚洲精品久久久久久婷婷小说| 精品亚洲成a人片在线观看| 欧美性感艳星| 肉色欧美久久久久久久蜜桃| 亚洲精品456在线播放app| 2018国产大陆天天弄谢| 国产欧美亚洲国产| 久久久久国产精品人妻一区二区| 亚洲av男天堂| 尾随美女入室| 亚洲欧美清纯卡通| 精品人妻偷拍中文字幕| 亚洲国产欧美日韩在线播放| 精品久久久久久电影网| 人人妻人人爽人人添夜夜欢视频| 五月天丁香电影| 视频区图区小说| 亚洲精品一二三| av不卡在线播放| 看非洲黑人一级黄片| 欧美变态另类bdsm刘玥| 亚洲综合精品二区| 欧美人与性动交α欧美精品济南到 | 黑人欧美特级aaaaaa片| 欧美最新免费一区二区三区| a级毛片黄视频| 国产一区二区三区av在线| 精品亚洲乱码少妇综合久久| 精品久久国产蜜桃| 久久久久国产网址| 一本久久精品| 亚洲av欧美aⅴ国产| 久久99热6这里只有精品| 欧美日韩精品成人综合77777| 国产亚洲精品第一综合不卡 | 99香蕉大伊视频| 欧美日韩亚洲高清精品| 七月丁香在线播放| 国产又色又爽无遮挡免| 精品少妇内射三级| 熟女av电影| 欧美日韩国产mv在线观看视频| 在线看a的网站| 建设人人有责人人尽责人人享有的| 日韩中字成人| 美女视频免费永久观看网站| 女的被弄到高潮叫床怎么办| 亚洲国产av新网站| 午夜影院在线不卡| 午夜久久久在线观看| 18禁国产床啪视频网站| 国产1区2区3区精品| 国产成人精品久久久久久| 最近的中文字幕免费完整| 久久99蜜桃精品久久| 国产69精品久久久久777片| 婷婷色综合大香蕉| 亚洲国产精品国产精品| 男女边吃奶边做爰视频| 国产在线视频一区二区| av国产久精品久网站免费入址| 男的添女的下面高潮视频| 999精品在线视频|