亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于卡爾曼濾波PID控制的精量排種器優(yōu)化設(shè)計(jì)與試驗(yàn)*

        2022-10-27 04:01:50劉思幸柴巖柳天虹李爽繆宏韓非
        關(guān)鍵詞:合格率

        劉思幸,柴巖,柳天虹,李爽,繆宏,韓非

        (1. 揚(yáng)州大學(xué)機(jī)械工程學(xué)院,江蘇揚(yáng)州,225127; 2. 揚(yáng)州大學(xué)信息工程學(xué)院,江蘇揚(yáng)州,225127)

        0 引言

        中國(guó)是農(nóng)業(yè)大國(guó),蔬菜產(chǎn)業(yè)是我國(guó)除糧食作物外栽培面積最廣、經(jīng)濟(jì)地位最重要的作物。而設(shè)施蔬菜種植打破了季節(jié)、氣候等條件的限制,保證了蔬菜的產(chǎn)量供應(yīng),提高了蔬菜產(chǎn)品質(zhì)量,有利于蔬菜生產(chǎn)的產(chǎn)業(yè)化[1-2]。

        設(shè)施蔬菜產(chǎn)業(yè)的快速發(fā)展,對(duì)其配套的播種裝備提出了新的要求。對(duì)于不同的設(shè)施蔬菜種植,其對(duì)應(yīng)的株距各不相同[3]。另外,為了有效利用土地,使蔬菜增產(chǎn)增效,設(shè)施蔬菜往往采用套種模式,不同蔬菜的株距、密度不一,需要每一行都有獨(dú)立的排種器根據(jù)不同株距調(diào)整轉(zhuǎn)速[4-6]。在播種機(jī)上設(shè)置合理的排種裝置,對(duì)解決設(shè)施蔬菜播種過(guò)程中的株距不一致問(wèn)題,具有極其重大的意義[7-9]。傳統(tǒng)排種驅(qū)動(dòng)方式容易由于地輪打滑造成漏播現(xiàn)象,降低株距合格率,采用步進(jìn)電機(jī)直接驅(qū)動(dòng)排種軸為該問(wèn)題提供了解決辦法,但是電機(jī)控制容易受到噪聲干擾,使靈敏度和精度降低,無(wú)法實(shí)現(xiàn)精量播種。對(duì)此,在控制系統(tǒng)中加入濾波算法,減少噪聲影響,成為播種機(jī)精量播種的重要研究方向[10-15]。隨著我國(guó)農(nóng)業(yè)機(jī)械化水平的不斷提升,對(duì)于小粒徑的設(shè)施蔬菜精量播種機(jī)的研究逐漸增多。解彬彬[16]、張靜[17]依據(jù)測(cè)速傳感器和角編碼器采集數(shù)據(jù),通過(guò)模糊算法自動(dòng)調(diào)整排種器轉(zhuǎn)速,實(shí)現(xiàn)精量排種。王冰媛[18]提出了小粒徑種子流體噴灑播種方法,設(shè)計(jì)了噴播式小粒種子流體播種機(jī),滿足蔬菜種子播種要求。孫靜鑫等[19-20]結(jié)合了環(huán)槽式和窩眼輪式兩種排種器優(yōu)點(diǎn),利用排種器的護(hù)種板輔助完成排種工作。李金鳳等[21-23]采用負(fù)壓內(nèi)側(cè)吸種,正壓吹氣投種的方式,解決播種機(jī)傳動(dòng)系統(tǒng)和氣流管路復(fù)雜、布局困難的問(wèn)題。國(guó)外發(fā)達(dá)國(guó)家的設(shè)施蔬菜精量播種機(jī)主要采用氣力式排種器進(jìn)行播種,技術(shù)成熟且裝備完善,具有較高的智能化、精量化和通用化水平[24-25]。綜上,國(guó)外上述播種機(jī)體型龐大,價(jià)格高昂,配件供應(yīng)不及時(shí),且土壤種植條件、氣候、降水等均與國(guó)內(nèi)存在明顯差異,并不適應(yīng)中國(guó)的設(shè)施蔬菜種植模式;而國(guó)內(nèi)學(xué)者對(duì)精量播種排種器結(jié)構(gòu)設(shè)計(jì)進(jìn)行了深入研究,但對(duì)排種裝置控制系統(tǒng)的無(wú)噪聲干擾方面的研究較少。

        本文針對(duì)我國(guó)設(shè)施機(jī)械的現(xiàn)狀和設(shè)施蔬菜種子精量播種過(guò)程中漏種、株距不一致、株距合格率低等問(wèn)題,設(shè)計(jì)了一種基于卡爾曼濾波PID控制方法的精量排種器裝置。并進(jìn)行旋轉(zhuǎn)正交臺(tái)架試驗(yàn),得出該排種裝置精量排種作業(yè)的最優(yōu)參數(shù)組合,驗(yàn)證所設(shè)計(jì)排種裝置的性能是否滿足設(shè)施蔬菜精量種植要求。

        1 結(jié)構(gòu)和工作原理

        1.1 排種器結(jié)構(gòu)

        如圖1、圖2所示,播種機(jī)排種裝置主要由排種步進(jìn)電機(jī)、排種器、種箱、測(cè)速地輪和編碼器5部分組成。播種機(jī)在機(jī)架側(cè)面設(shè)有步進(jìn)電機(jī)安裝架,安裝架通過(guò)螺栓固定在機(jī)架上,步進(jìn)電機(jī)通過(guò)聯(lián)軸器與排種軸連接,步進(jìn)電機(jī)外側(cè)裝有編碼器,測(cè)速地輪置于機(jī)架前方,用于檢測(cè)牽引機(jī)前進(jìn)速度,輪架上設(shè)有拉伸彈簧,使地輪保持貼地,并具有減震功能。

        (a) 排種器關(guān)鍵部件 (b) 地輪附件

        圖2 排種部件結(jié)構(gòu)圖

        1.2 工作原理

        播種機(jī)工作時(shí),首先由整地輪對(duì)播種之前的作業(yè)平面進(jìn)行平整,隨后由監(jiān)測(cè)地輪上的脈沖采集裝置采集到地輪轉(zhuǎn)動(dòng)生成的脈沖信號(hào),結(jié)合地輪直徑計(jì)算出播種機(jī)前進(jìn)速度,之后控制器通過(guò)計(jì)算采集到的播種機(jī)行進(jìn)的速度信號(hào),針對(duì)不同設(shè)施蔬菜株距要求,驅(qū)動(dòng)排種步進(jìn)電機(jī)工作,進(jìn)而驅(qū)動(dòng)排種輪的轉(zhuǎn)動(dòng),并加以反饋控制,從而完成整個(gè)排種作業(yè)。

        2 控制系統(tǒng)設(shè)計(jì)

        2.1 總體設(shè)計(jì)

        排種控制系統(tǒng)由主控制器、電平轉(zhuǎn)換單元、執(zhí)行單片機(jī)、步進(jìn)電機(jī)驅(qū)動(dòng)器、步進(jìn)電機(jī)、編碼器和排種軸組成。主控制器和執(zhí)行器通過(guò)RS-485信號(hào)進(jìn)行排種指令的發(fā)送和執(zhí)行的反饋,主控制器發(fā)出排種數(shù)字信號(hào),執(zhí)行器將接收到的信號(hào)轉(zhuǎn)為脈沖信號(hào)下發(fā)給驅(qū)動(dòng)器,進(jìn)而驅(qū)動(dòng)步進(jìn)電機(jī)帶動(dòng)排種軸轉(zhuǎn)動(dòng)相應(yīng)角度,設(shè)施種子沿排種盤下落,完成一次排種作業(yè)。其控制系統(tǒng)總體框圖如圖3所示,CAN總線通訊和步進(jìn)電機(jī)驅(qū)動(dòng)電路設(shè)計(jì)圖如圖4所示。

        圖3 控制系統(tǒng)總體框圖

        (a) CAN總線通訊電路

        (b) 步進(jìn)電機(jī)驅(qū)動(dòng)電路

        2.2 控制原理

        設(shè)施蔬菜精量播種需要對(duì)排種器排種進(jìn)行精準(zhǔn)控制,而在實(shí)際播種過(guò)程中,地塊不平整、打滑、磁場(chǎng)等因素嚴(yán)重影響排種器排種精度,引入一種卡爾曼濾波PID控制方法,對(duì)步進(jìn)電機(jī)進(jìn)行實(shí)時(shí)反饋控制。播種過(guò)程中,控制器根據(jù)用戶設(shè)定株距值并結(jié)合測(cè)速地輪采集的實(shí)時(shí)速度數(shù)據(jù),計(jì)算出當(dāng)前電機(jī)所需脈沖數(shù),與編碼器采集的實(shí)時(shí)電機(jī)轉(zhuǎn)速進(jìn)行分析和判斷,對(duì)于產(chǎn)生的數(shù)據(jù)偏差通過(guò)控制方法及時(shí)修正,從而產(chǎn)生當(dāng)前應(yīng)該輸出的控制信號(hào),經(jīng)過(guò)執(zhí)行機(jī)構(gòu)施加到步進(jìn)電機(jī)上,實(shí)現(xiàn)精量排種的目的。

        輸入脈沖數(shù)與理論株距計(jì)算關(guān)系如下:株距調(diào)節(jié)依靠步進(jìn)電機(jī)獨(dú)立驅(qū)動(dòng)實(shí)現(xiàn),設(shè)理論設(shè)定株距為s(cm),排種器轉(zhuǎn)速即步進(jìn)電機(jī)轉(zhuǎn)速為n1(r/min),排種器每轉(zhuǎn)排出種子數(shù)為k,則在時(shí)間t(s)內(nèi)播種機(jī)前進(jìn)的距離

        L=n1tks

        (1)

        因此要滿足在時(shí)間t內(nèi)排出種子的株距要求,牽引機(jī)需要前進(jìn)同樣的距離L(m)。已知牽引機(jī)前進(jìn)速度由測(cè)速地輪獲得,具體為由安裝在地輪上的編碼器采集到時(shí)間t內(nèi)的脈沖數(shù)計(jì)算而得,設(shè)編碼器的分辨率為P(即為編碼器旋轉(zhuǎn)一周輸出的脈沖數(shù)),時(shí)間t內(nèi)編碼器輸出的脈沖數(shù)為ΔP,則編碼器測(cè)得地輪的轉(zhuǎn)速

        (2)

        設(shè)測(cè)速地輪的直徑為d(mm),可得時(shí)間t內(nèi)地輪走過(guò)的距離LZ(m),有

        LZ=πdnt

        (3)

        可知LZ即牽引機(jī)前進(jìn)的距離,即LZ=L,將式(1)、式(3)聯(lián)立,可得

        πdn=n1ks

        (4)

        進(jìn)而可得牽引機(jī)前進(jìn)速度V(m/s),有

        (5)

        步進(jìn)電機(jī)轉(zhuǎn)動(dòng)靠輸入脈沖信號(hào)實(shí)現(xiàn),每輸入一個(gè)脈沖信號(hào)給電機(jī),轉(zhuǎn)子轉(zhuǎn)過(guò)一個(gè)角度(步距角),設(shè)步距角為β,輸入脈沖數(shù)為kc,則步進(jìn)電機(jī)轉(zhuǎn)速

        (6)

        將式(5)、式(6)聯(lián)立,可知要滿足理論設(shè)定株距的要求,需要向步進(jìn)電機(jī)輸入的脈沖數(shù)

        (7)

        2.3 控制算法實(shí)現(xiàn)

        傳統(tǒng)PID算法存在參數(shù)整定不良、超調(diào)大、適應(yīng)性差、控制性能不佳、干擾后調(diào)整時(shí)間長(zhǎng)、穩(wěn)態(tài)誤差大等問(wèn)題,不適用于精量播種系統(tǒng)調(diào)節(jié),并且由于噪聲的影響,控制器的參數(shù)難以自動(dòng)調(diào)節(jié),不能達(dá)到理想的控制效果。在傳統(tǒng)PID控制系統(tǒng)中添加極小極快無(wú)量綱高斯白噪聲后,并選取其調(diào)節(jié)參數(shù)分別為10、5、1,其MATLAB中Simulink控制模型搭建和仿真結(jié)果如圖5所示,可以看出,此時(shí)由于微分項(xiàng)的影響,在系統(tǒng)達(dá)到穩(wěn)定值后處于高頻波動(dòng)狀態(tài),不利于電機(jī)反饋控制,難以實(shí)現(xiàn)精量排種。

        (a) Simulink控制模型

        卡爾曼濾波是一種高效率的遞歸濾波器,它能夠從一系列的不完全及包含噪聲的測(cè)量中,估計(jì)動(dòng)態(tài)系統(tǒng)的狀態(tài)。其可以在貝葉斯濾波的基礎(chǔ)上做出如下假設(shè),其中,式(8)為預(yù)測(cè)方程,式(9)為觀測(cè)方程。

        XK=FXK-1+QK

        (8)

        YK=hXK+RK

        (9)

        式中:QK——預(yù)測(cè)噪聲,QK服從高斯分布;

        RK——觀測(cè)噪聲,RK服從高斯分布。

        根據(jù)貝葉斯濾波的預(yù)測(cè)步和更新步方程可計(jì)算出卡爾曼濾波的預(yù)測(cè)步和更新步方程為

        (10)

        K——卡爾曼增益。

        對(duì)傳統(tǒng)PID控制器進(jìn)行卡爾曼濾波融合后,在MATLAB中編寫控制程序并在Simulink中搭建控制模型完成仿真測(cè)試,其控制仿真模型如圖6所示,系統(tǒng)響應(yīng)結(jié)果如圖7所示,由結(jié)果可以得出,采用卡爾曼濾波后,對(duì)于相同量級(jí)的高斯白噪聲干擾,其信號(hào)的諧波幅值明顯降低,穩(wěn)定性提高,趨近于理想信號(hào)。顯然,卡爾曼濾波的引入,對(duì)干擾起到了抑制與降低,表明該系統(tǒng)有著較好的噪聲抑制能力,這對(duì)于電機(jī)精確控制及排種器精量播種的實(shí)現(xiàn)起到了關(guān)鍵作用。

        圖6 控制仿真模型

        (a) 融合卡爾曼濾波前

        3 排種試驗(yàn)

        3.1 評(píng)價(jià)指標(biāo)

        試驗(yàn)選取設(shè)施蔬菜種子中具有代表性的丸?;S瓜種子,千粒質(zhì)量為28 g,含水率為7%。根據(jù)GB/T 6973—2005《單粒(精密)播種機(jī)試驗(yàn)方法》規(guī)定的試驗(yàn)方法和指標(biāo),選取株距合格率、重播率、漏播率和株距變異系數(shù)作為排種性能評(píng)價(jià)指標(biāo),進(jìn)行設(shè)施蔬菜智能精量播種機(jī)排種試驗(yàn)評(píng)價(jià)。

        1) 株距合格率

        (11)

        式中:n1——株距合格數(shù);

        N′——區(qū)間數(shù)。

        2) 重播率

        (12)

        式中:n2——重播數(shù)。

        3) 漏播率

        (13)

        式中:n0——漏播數(shù)。

        4) 株距變異系數(shù)

        (14)

        (15)

        C=σ×100%

        (16)

        式中:C——株距變異系數(shù),%;

        σ——標(biāo)準(zhǔn)差。

        3.2 試驗(yàn)方案與結(jié)果分析

        3.2.1 單因素試驗(yàn)

        通過(guò)實(shí)地調(diào)試,確定排種盤轉(zhuǎn)速和播種機(jī)行走速度是影響設(shè)施蔬菜種子株距水平的主要因素。搭建排種試驗(yàn)臺(tái)架,如圖8所示。

        圖8 臺(tái)架試驗(yàn)圖

        為了得到理想穩(wěn)定的株距水平,在排種試驗(yàn)臺(tái)上對(duì)排種盤轉(zhuǎn)速和作業(yè)速度進(jìn)行單因素試驗(yàn),確定排種盤轉(zhuǎn)速20 r/min、行走速度2 km/h 為排種器單因素試驗(yàn)初步優(yōu)選值。

        3.2.2 正交旋轉(zhuǎn)多因素試驗(yàn)

        為尋求排種盤轉(zhuǎn)速和作業(yè)速度兩個(gè)因素對(duì)于株距水平影響的最優(yōu)參數(shù)組合,以株距合格率、重播率、漏播率和株距變異系數(shù)作為評(píng)價(jià)指標(biāo),進(jìn)行兩因素五水平二次回歸正交旋轉(zhuǎn)組合試驗(yàn)。根據(jù)單因素試驗(yàn)結(jié)果確定試驗(yàn)因素水平設(shè)置,試驗(yàn)因素水平編碼表如表1所示,試驗(yàn)方案及試驗(yàn)結(jié)果如表2所示。

        表1 試驗(yàn)因素水平Tab. 1 Experimental factors and levels

        運(yùn)用Design Expert 8.0.6軟件對(duì)試驗(yàn)結(jié)果進(jìn)行方差分析,并進(jìn)行多元回歸擬合,得到株距合格率Y1、重播率Y2、漏播率Y3和株距變異系數(shù)Y4各試驗(yàn)指標(biāo)的回歸方程。其方差分析結(jié)果如表3~表6所示。

        Y1=94.91-0.57X1-1.59X2+0.13X1X2-

        0.31X12+0.38X22

        (17)

        Y2=2.47-0.20X1-1.20X2+0.019X1X2-

        0.12X12-0.30X22

        (18)

        Y3=2.88+0.21X1+1.40X2-0.15X1X2-

        0.20X12-0.077X22

        (19)

        Y4=14.43+0.77X1+2.10X2-0.59X1X2-

        0.42X12-0.55X22

        (20)

        由表3~表6的回歸模型方差分析結(jié)果可知,四個(gè)回歸模型的P值均小于0.000 1,表明回歸模型高度顯著;在試驗(yàn)因素中,由排種盤轉(zhuǎn)速和行走速度的P值可判斷行走速度相比于排種盤轉(zhuǎn)速對(duì)4個(gè)試驗(yàn)指標(biāo)的影響更為顯著,株距合格率、重播率、漏播率和株距變異系數(shù)回歸模型中的P值均小于0.05,表明4個(gè)回歸項(xiàng)在回歸模型中交互影響顯著。對(duì)于重播率Y2,回歸項(xiàng)X22也影響極顯著;對(duì)于株距變異系數(shù)Y4,回歸項(xiàng)X1X2也影響顯著。四模型的決定系數(shù)R2均接近1,變異系數(shù)分別為0.83%、0.96%、0.95%、0.87%,說(shuō)明該株距合格率、重播率、漏播率和株距變異系數(shù)擬合回歸模型具有較高的可靠性。

        表3 株距合格率方差分析Tab. 3 Variance analysis of qualified rate of plant spacing

        表4 重播率方差分析Tab. 4 Variance analysis of multiple rate

        表5 漏播率方差分析Tab. 5 Variance analysis of missing rate

        表6 株距變異系數(shù)方差分析Tab. 6 Variance analysis of coefficient of

        經(jīng)過(guò)數(shù)據(jù)處理和分析,通過(guò)Design Expert 8.0.6軟件得到各因素對(duì)于試驗(yàn)指標(biāo)的響應(yīng)曲面。排種盤轉(zhuǎn)速和行走速度對(duì)株距合格率、重播率、漏播率和株距變異系數(shù)的響應(yīng)曲面分別如圖9~圖12所示,行走速度由1 km/h提高到3 km/h過(guò)程中,排種盤轉(zhuǎn)速越快,株距合格率和重播率越低;當(dāng)行走速度由快變慢時(shí),株距合格率和重播率隨著排種盤轉(zhuǎn)速的減小而提高。當(dāng)排種盤轉(zhuǎn)速由10 r/min增加到30 r/min過(guò)程中,行走速度越高,漏播率和株距變異系數(shù)越高;當(dāng)排種盤轉(zhuǎn)速由高變低時(shí),漏播率和株距變異系數(shù)隨著行走速度的減小而減小。

        圖9 株距合格率指數(shù)的雙因素響應(yīng)曲面

        圖10 重播率指數(shù)的雙因素響應(yīng)曲面

        圖11 漏播率指數(shù)的雙因素響應(yīng)曲面

        圖12 株距變異系數(shù)指數(shù)的雙因素響應(yīng)曲面

        隨著行走速度的提高,株距合格率和重播率呈下降趨勢(shì),而漏播率和株距變異系數(shù)呈上升趨勢(shì);在同一設(shè)定速度下,隨著排種盤轉(zhuǎn)速的提高,株距合格率和重播率呈下降趨勢(shì),而漏播率和株距變異系數(shù)呈上升趨勢(shì)。排種盤轉(zhuǎn)速一致時(shí),行走速度減小,重播率的增幅較大;行進(jìn)速度增大,漏播率的增幅相對(duì)較大。基于交互效應(yīng)分析可知,排種盤轉(zhuǎn)速和行走速度同時(shí)減小時(shí),株距合格率指數(shù)、漏播率指數(shù)和株距變異系數(shù)指數(shù)顯優(yōu);排種盤轉(zhuǎn)速和行走速度同時(shí)增大時(shí),重播率指數(shù)顯優(yōu)。

        3.3 參數(shù)優(yōu)化

        為了得到各試驗(yàn)因素的最佳水平組合,運(yùn)用Design Expert 8.0.6軟件中的Optimization功能對(duì)回歸模型進(jìn)行求解,以株距合格率最大、重播率最小、漏播率最小、株距變異系數(shù)最小為條件,求解模型得到的排種盤轉(zhuǎn)速為10 r/min,行走速度為1.6 km/h時(shí),滿足上述條件,其株距合格率為95.9%,重播率為2.9%,漏播率為1.9%,株距變異系數(shù)為12.1%,均符合標(biāo)準(zhǔn)規(guī)定。

        4 結(jié)論

        1) 針對(duì)設(shè)施蔬菜排種器存在的漏播、重播、精度低等問(wèn)題,對(duì)排種器裝置進(jìn)行了結(jié)構(gòu)和控制系統(tǒng)優(yōu)化設(shè)計(jì),利用步進(jìn)電機(jī)直接驅(qū)動(dòng)排種器進(jìn)行排種作業(yè),通過(guò)控制系統(tǒng)實(shí)現(xiàn)精量播種,實(shí)驗(yàn)得出影響設(shè)施蔬菜精量播種的主要因素為排種盤的轉(zhuǎn)速和機(jī)具的行走速度。

        2) 對(duì)傳統(tǒng)PID控制器進(jìn)行優(yōu)化設(shè)計(jì),融合卡爾曼濾波,仿真結(jié)果表明,其可以有效減少噪聲對(duì)電機(jī)的干擾,對(duì)提升系統(tǒng)控制過(guò)程的平穩(wěn)性起到重要作用。

        3) 通過(guò)分析排種盤轉(zhuǎn)速和行走輪行走速度對(duì)排種質(zhì)量的單因素試驗(yàn),得到理想排種質(zhì)量的各參數(shù)合理變化范圍:排種盤轉(zhuǎn)速為10~30 r/min,行走輪行走速度為1~3 km/h。

        4) 采用兩因素五水平正交試驗(yàn),建立各個(gè)因素與評(píng)價(jià)指標(biāo)間的回歸模型,得到影響株距合格率、重播率、漏播率和株距變異系數(shù)的因素為:行走速度>排種盤轉(zhuǎn)速。確定最佳參數(shù)組合:當(dāng)排種盤轉(zhuǎn)速為10 r/min,行進(jìn)速度為1.6 km/h時(shí),株距合格率為95.9%,重播率為2.9%,漏播率為1.9%,株距變異系數(shù)為12.1%,滿足設(shè)施蔬菜的精量播種要求。

        猜你喜歡
        合格率
        湖北化肥質(zhì)量監(jiān)督抽查合格率96%
        2020年畜禽產(chǎn)品獸藥殘留抽檢樣品合格率99.49%
        2020年農(nóng)產(chǎn)品質(zhì)量安全例行監(jiān)測(cè)合格率達(dá)97.8%
        提高肉雞出欄合格率的綜合措施
        提高某傘形罩加工合格率工藝探討
        提高鋁電解電容器合格率的多元化工藝技術(shù)
        電子制作(2018年12期)2018-08-01 00:47:52
        淺談如何提高遙測(cè)狀態(tài)估計(jì)合格率
        淺談如何提高電網(wǎng)狀態(tài)估計(jì)遙測(cè)估計(jì)合格率
        提高重癥醫(yī)學(xué)科床頭抬高合格率
        提高精神科住院患者衛(wèi)生合格率
        国产91在线|亚洲| 亚洲在中文字幕乱码熟女| 亚洲熟少妇一区二区三区| 视频女同久久久一区二区| 丰满精品人妻一区二区| 秋霞在线视频| 天天噜日日噜狠狠噜免费| 男人j进女人j啪啪无遮挡| 日韩精品一区二区亚洲av| 国产鲁鲁视频在线播放| 无码中文av有码中文av| 在线一区二区三区免费视频观看| 亚洲中文字幕久久精品一区| 久久婷婷五月综合97色直播| 天天躁狠狠躁狠狠躁夜夜躁| 中字无码av电影在线观看网站| 亚洲中文字幕精品久久久久久直播| 中文天堂一区二区三区| 精品高清一区二区三区人妖| 亚洲av中文无码乱人伦在线咪咕| 免费操逼视频| 日韩人妻精品无码一区二区三区| 亚洲成年网站在线777| 91麻豆精品一区二区三区| 成人国产精品一区二区八戒网| 国产精品久久久久久久久绿色| 中文字幕人妻av一区二区| 成人免费xxxxx在线视频| 美腿丝袜美腿国产在线| 蜜臀人妻精品一区二区免费| 国产情侣真实露脸在线| 狠狠躁夜夜躁人人爽天天天天97| 四虎无码精品a∨在线观看| 国产精品自拍网站在线| 欧美亚洲一区二区三区| 乱子伦视频在线看| 国产小车还是日产的好 | 视频在线观看一区二区三区| 久久伊人少妇熟女大香线蕉| 久久久精品人妻一区亚美研究所| 一区二区三区不卡免费av|