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

    改進(jìn)適應(yīng)度函數(shù)的CMA-ES 算法在機(jī)器人逆運(yùn)動(dòng)學(xué)中的應(yīng)用

    2022-05-11 07:31:36譚薪興薛晨慷于權(quán)偉
    關(guān)鍵詞:運(yùn)動(dòng)學(xué)逆向適應(yīng)度

    譚薪興, 李 光, 薛晨慷, 易 靜, 于權(quán)偉

    (湖南工業(yè)大學(xué) 機(jī)械工程學(xué)院, 湖南 株洲 412007)

    0 引 言

    機(jī)器人運(yùn)動(dòng)學(xué)包括正向運(yùn)動(dòng)學(xué)和逆向運(yùn)動(dòng)學(xué),正向運(yùn)動(dòng)學(xué)是由一組關(guān)節(jié)角度確定機(jī)器人末端執(zhí)行器一個(gè)姿態(tài),逆向運(yùn)動(dòng)學(xué)是由目標(biāo)姿態(tài)求解各關(guān)節(jié)角度,逆向運(yùn)動(dòng)學(xué)求解復(fù)雜且存在多解。

    機(jī)器人逆向運(yùn)動(dòng)學(xué)求解是進(jìn)行機(jī)器人軌跡規(guī)劃、路徑規(guī)劃、避障等研究的基礎(chǔ)。 很多學(xué)者在該領(lǐng)域做了大量研究,有許多逆向運(yùn)動(dòng)學(xué)的求解理論和方法,傳統(tǒng)的方法有代數(shù)法、幾何法、數(shù)值法等。代數(shù)法主要以消元的方式將高維方程組轉(zhuǎn)化為低維方程組,還需要大量的三角代換,計(jì)算過(guò)程十分復(fù)雜,文獻(xiàn)[1]提出基于分離重構(gòu)技術(shù),降低了求解復(fù)雜度,得到了逆運(yùn)動(dòng)學(xué)全部解;幾何法針對(duì)機(jī)器人的特殊結(jié)構(gòu)進(jìn)行簡(jiǎn)化再求解,一般無(wú)法單獨(dú)使用甚至無(wú)法使用;數(shù)值法可以精確搜索到一組逆解,求解速度與初始值相關(guān),甚至存在無(wú)解,需要的運(yùn)算量大,不適合實(shí)時(shí)性要求較高的場(chǎng)合。

    近年來(lái),利用智能算法求解機(jī)器人逆向運(yùn)動(dòng)學(xué)的方法取得一定的成果。 神經(jīng)網(wǎng)絡(luò)算法的計(jì)算代價(jià)高,計(jì)算能力很大程度取決于數(shù)據(jù)的大小,同時(shí)也取決于網(wǎng)絡(luò)的深度和復(fù)雜程度;遺傳算法全局搜索能力強(qiáng),局部搜索能力較弱,容易誤入局部最優(yōu)解;蟻群算法和魚(yú)群算法參數(shù)設(shè)置復(fù)雜,如果參數(shù)設(shè)置不當(dāng),容易偏離最優(yōu)解;粒子群算法操作簡(jiǎn)單,收斂速度快,在眾多領(lǐng)域得到廣泛應(yīng)用,但存在提前收斂、維數(shù)災(zāi)難等問(wèn)題。

    自適應(yīng)協(xié)方差矩陣進(jìn)化策略(Covariance MatrixAdaptation Evolution Strategy, CMA-ES)是在進(jìn)化策略(Evolution Strategy ,ES)算法的基礎(chǔ)上發(fā)展而來(lái)的無(wú)約束優(yōu)化算法,在全局優(yōu)化、多峰優(yōu)化、多目標(biāo)優(yōu)化、大規(guī)模優(yōu)化和結(jié)構(gòu)性工程等領(lǐng)域得到了大量應(yīng)用,該算法的搜索性能較高,并且能夠達(dá)到比較高的搜索精度。 文獻(xiàn)[6]提出了改進(jìn)的CMA-ES 算法,在保留原算法的優(yōu)點(diǎn)上,可以實(shí)時(shí)且高精度地求出逆向運(yùn)動(dòng)學(xué)解;文獻(xiàn)[7]引入1/5 成功理論,調(diào)整CMA-ES 算法的規(guī)模和步長(zhǎng),提高了算法求解效率;上述智能算法提高了機(jī)器人逆向運(yùn)動(dòng)學(xué)求解的效率,但是忽略了機(jī)器人在實(shí)際運(yùn)動(dòng)過(guò)程中的平穩(wěn)性和能耗。

    本文提出了一種新的適應(yīng)度函數(shù)用于CMA-ES算法求逆運(yùn)動(dòng)學(xué)解:與原算法相比,該算法將各關(guān)節(jié)運(yùn)動(dòng)范圍作為約束條件,通過(guò)加權(quán)最小二乘法和位姿誤差建立新的適應(yīng)度函數(shù),求解出唯一的一組運(yùn)動(dòng)學(xué)逆解。 將該算法應(yīng)用于REBot-V-6R-6500 型六自由度機(jī)器人,仿真結(jié)果表明,該方法可以得到機(jī)器人高精度和平穩(wěn)的逆運(yùn)動(dòng)學(xué)解。

    1 CMA-ES 算法

    1.1 算法原理

    CMA-ES 算法是一種進(jìn)化策略類(lèi)算法。 經(jīng)典ES 算法尋找最優(yōu)解主要依靠突變,但是如何調(diào)整突變的方向沒(méi)有成熟的理論支撐,需要依據(jù)經(jīng)驗(yàn)設(shè)置,因此該方法會(huì)導(dǎo)致無(wú)效的突變。 為克服經(jīng)典ES 的不足,CMA-ES 采用正態(tài)分布(,) 生成一定數(shù)量的搜索種群,(,) 表示搜索分布均值為,以協(xié)方差矩陣為的多元正態(tài)分布,且,其中,的列向量由的特征向量正交基組成,為 對(duì)角矩陣,化簡(jiǎn)多元正態(tài)分布(,) 得式(1)。

    其中,“~”表示服從相同的分布。

    由式(1)的逆序可知,確定等概率密度圓球面,通過(guò)、的變化可以確定的橢球面;球面的分布尺度由中的對(duì)角元素決定,主軸的方向由確定,變化、實(shí)現(xiàn)球面的分布旋轉(zhuǎn)。

    1.2 算法步驟

    CAM-ES 算法流程如圖1 所示。

    圖1 CAM-ES 算法流程圖Fig.1 CAM-ES algorithm flow chart

    突變。 突變的目的是使個(gè)體之間產(chǎn)生差異,CMA-ES 算法的突變過(guò)程是以當(dāng)代的均值為中心,產(chǎn)生下一代種群,式(2)。

    選擇和重組。 經(jīng)典ES 算法的選擇策略有兩種:一種是() 策略,適應(yīng)度排名由上一代的父代與當(dāng)前子代一起參與,同時(shí)選擇優(yōu)秀的個(gè)體;另一種是(,) 策略,即父代不參與當(dāng)前代的適應(yīng)度排名競(jìng)爭(zhēng),優(yōu)秀個(gè)體只從當(dāng)前子代中選擇。 CMA -ES 采用的是第二種選擇策略,對(duì)個(gè)子代進(jìn)行適應(yīng)度評(píng)價(jià),根據(jù)適應(yīng)能力由高到低排序,選取適應(yīng)能力強(qiáng)的前個(gè)個(gè)體,作為下一代種群的父代,更新策略參數(shù)、,以此來(lái)傳遞優(yōu)良的基因。

    (1) 均值移動(dòng):下一代起始搜索點(diǎn),即新的均值是從樣本,…,中選擇出來(lái)的個(gè)最優(yōu)子群的加權(quán)平均值,式(3) 和式(4)。

    (2)協(xié)方差矩陣自適應(yīng)調(diào)整(CMA),式(5) ~式(9)。

    其中,cp的更新學(xué)習(xí)速率;h為赫維賽德階躍()函數(shù),用于調(diào)整‖p的過(guò)大增長(zhǎng);μ為方差有效選擇質(zhì)量,且1 <μ≤,和c分別為的“秩1” 和“秩” 的更新學(xué)習(xí)速率。

    (3)全局步長(zhǎng)控制(CSA),式(10)和式(11)。

    其中,可看作步長(zhǎng)變化的伸縮因子;cp的更新學(xué)習(xí)率;d為接近于1 的阻尼系數(shù);‖(0,)‖為歸一化進(jìn)化路徑在隨機(jī)選擇下的期望長(zhǎng)度。

    判斷是否達(dá)到最大迭代次數(shù)或精度,若是,則停止,輸出最優(yōu)解和最優(yōu)決策向量,否則返回Step 2。

    2 逆運(yùn)動(dòng)學(xué)計(jì)算

    2.1 正向運(yùn)動(dòng)學(xué)數(shù)學(xué)模型

    采用標(biāo)準(zhǔn)D-H 法建立的REBot-V-6R-6500型六自由度機(jī)器人結(jié)構(gòu)簡(jiǎn)圖,如圖2 所示。 目前,普遍認(rèn)為D-H 法是對(duì)機(jī)器人運(yùn)動(dòng)學(xué)建模最好的方法,可以確定機(jī)器人的關(guān)節(jié)參數(shù)和關(guān)節(jié)變量,與機(jī)器人結(jié)構(gòu)順序和復(fù)雜程度無(wú)關(guān),通過(guò)齊次變換矩陣描述兩個(gè)相鄰坐標(biāo)系間的空間位姿關(guān)系。

    圖2 REBot-V-6R-6500 結(jié)構(gòu)簡(jiǎn)圖Fig.2 Structure diagram of REBot-V-6R-6500

    首先,建立連桿坐標(biāo)系,通過(guò)4 次旋轉(zhuǎn)平移變換,根據(jù)DH 參數(shù),第連桿的參數(shù)定義為: a為連桿長(zhǎng)度;α為相鄰兩關(guān)節(jié)軸線的夾角;d為第根連桿和第1 連桿的偏置距離;θ為第連桿的旋轉(zhuǎn)角。 變換過(guò)程:繞z 軸旋轉(zhuǎn)θ角,沿z 軸平移d,沿x軸平移a, 繞x軸旋轉(zhuǎn)α,可以得到機(jī)器人連桿坐標(biāo)系相對(duì)于連桿1 坐標(biāo)系的變換矩陣(矩陣)即式(12), REBot-V-6R-6500 機(jī)器人關(guān)節(jié)參數(shù),見(jiàn)表1。

    表1 REBot-V-6R-6500 機(jī)器人關(guān)節(jié)參數(shù)Tab.1 Joint parameters of REBot-V-6R-6500 robot

    根據(jù)式(12)和表(1)中的參數(shù),得到該型號(hào)六自由度機(jī)器人正向運(yùn)動(dòng)學(xué)數(shù)學(xué)模型式(13),通過(guò)正向運(yùn)動(dòng)學(xué)控制機(jī)器人末端執(zhí)行器到達(dá)指定位置。

    2.2 適應(yīng)度函數(shù)

    用機(jī)器人末端執(zhí)行器的接近矢量和位置矢量建立適應(yīng)度函數(shù)和,式(14)。 適應(yīng)度函數(shù)為位置誤差范數(shù),表示目標(biāo)位置P和實(shí)際位置P之間的誤差;適應(yīng)度函數(shù)為姿態(tài)誤差范數(shù),表示目標(biāo)位置a和實(shí)際位置a之間的誤差;適應(yīng)度函數(shù)表示位置誤差范數(shù)和姿態(tài)誤差范數(shù)之和。 系數(shù)和是為方便調(diào)整求解過(guò)程中位置和姿態(tài)誤差而設(shè)置的權(quán)重,以便求得最合適的逆解,本文中1,1。

    六自由度機(jī)器人逆運(yùn)動(dòng)學(xué)存在多解,為了獲得唯一的逆解,本文基于加權(quán)最小二乘法,結(jié)合“最佳柔順性”規(guī)則,通過(guò)附加約束建立目標(biāo)函數(shù),使用CMA-ES 算法求解適應(yīng)度函數(shù)最小值,獲得最優(yōu)解,最佳柔順性的含義為式(15):

    其中,θ() -θ(1) 表示關(guān)節(jié)角與其自身前一個(gè)位置關(guān)節(jié)角的差值,使得關(guān)節(jié)運(yùn)動(dòng)更平滑,減少機(jī)器人運(yùn)動(dòng)路徑距離;系數(shù)是為了保證算法運(yùn)行時(shí)不產(chǎn)生局部收斂甚至發(fā)散等情況,平衡“最佳柔順性”準(zhǔn)則對(duì)整個(gè)適應(yīng)度函數(shù)優(yōu)化結(jié)果的影響。 本文中0.001;系數(shù)w遵循“少運(yùn)動(dòng)上臂,多運(yùn)動(dòng)下臂”的選取原則。 該原則的優(yōu)點(diǎn)是運(yùn)動(dòng)時(shí)更為平穩(wěn),減少能源消耗,提高工作效率等。 在本研究中,加權(quán)系數(shù)w設(shè)計(jì)為[3 2 1.5 1 1 1]。

    3 仿真與結(jié)果分析

    3.1 求解工作空間

    機(jī)器人的關(guān)節(jié)角取值范圍參考該型號(hào)機(jī)器人說(shuō)明書(shū),見(jiàn)表2。 求解機(jī)器人工作空間方法有:解析法、幾何法、數(shù)值法等,但這幾種計(jì)算復(fù)雜、效率低,本文在機(jī)器人各關(guān)節(jié)允許的角度值限制內(nèi),以蒙特卡洛法為基礎(chǔ),用隨機(jī)函數(shù)給各關(guān)節(jié)角賦值,得到10 000組關(guān)節(jié)角,計(jì)算機(jī)器人可達(dá)工作空間;利用正向運(yùn)動(dòng)學(xué)公式(13)解出末端執(zhí)行器在笛卡爾空間中離散點(diǎn)圖,從而得到機(jī)器人三維工作空間,如圖3所示。

    表2 關(guān)節(jié)角取值范圍Tab.2 Range of joint angle

    圖3 三維工作空間Fig.3 Three-dimensional workspace

    3.2 參數(shù)設(shè)置

    CMA-ES 算法的種群數(shù)量100,優(yōu)秀個(gè)體50,初始步長(zhǎng)05,初始均值在0~1 內(nèi)隨機(jī)產(chǎn)生,求解過(guò)程中,前一個(gè)逆向運(yùn)動(dòng)學(xué)解作為當(dāng)前初始點(diǎn)均值,初始均值縮放系數(shù)0.65,算法的停止條件為進(jìn)化次數(shù)200 或者適應(yīng)度函數(shù)≤10。為了保證算法穩(wěn)定運(yùn)行,兩種算法的適應(yīng)度函數(shù)中,權(quán)值系數(shù)1,1,0001。

    3.3 點(diǎn)到點(diǎn)運(yùn)動(dòng)仿真

    仿真的計(jì)算機(jī)配置:操作系統(tǒng)為64 位Window 10 企業(yè)版,處理器為AMD Ryzen 5 2600X 3.60 GHz,內(nèi)存8.0 GB。

    機(jī)器人初始狀態(tài)時(shí)選取各關(guān)節(jié)角度為[0 0 0 0 0 0],目標(biāo)位置選取工作空間內(nèi)的一組關(guān)節(jié)角,根據(jù)正向運(yùn)動(dòng)學(xué)公式(13)求出末端執(zhí)行器的位姿矩陣, 分別使用適應(yīng)度函數(shù)和, 在CMA-ES 算法上獨(dú)立運(yùn)行1 000 次,求逆運(yùn)動(dòng)學(xué)的解。 隨機(jī)選取機(jī)器人的一組關(guān)節(jié)角度為[0.412 0.365 0.254 0.121 0.454 0.235]。

    兩種適應(yīng)度函數(shù)的算法單獨(dú)運(yùn)行1 000 次的結(jié)果見(jiàn)表3,結(jié)果表明改進(jìn)后適應(yīng)度函數(shù)代表位置誤差的適應(yīng)度函數(shù)的平均值比原適應(yīng)度函數(shù)中的平均值小9 個(gè)數(shù)量級(jí);任意選取1 000 組結(jié)果中的一組的位置誤差結(jié)果見(jiàn)表4,可以看出改進(jìn)適應(yīng)度函數(shù)的位置誤差值比未改進(jìn)的適應(yīng)度函數(shù)求得的位置精度高10數(shù)量級(jí)。

    表3 CMA-ES 獨(dú)立運(yùn)行1 000 次的適應(yīng)度函數(shù)f1的值Tab.3 The value of the fitness function f1 for CMA-ES running independently for 1 000 times

    表4 點(diǎn)對(duì)點(diǎn)運(yùn)動(dòng)位置誤差Tab.4 Position error of point-to-point motion

    3.4 連續(xù)軌跡跟蹤仿真

    在工作空間內(nèi)的軌跡方程為式(16):

    其中,0≤≤2π,步長(zhǎng)為π/60(rad),長(zhǎng)度單位為m,為了滿(mǎn)足實(shí)際工作需求,保證末端執(zhí)行器關(guān)節(jié)軸在工作過(guò)程中始終平行于世界坐標(biāo)系軸正方向,運(yùn)動(dòng)軌跡的每一個(gè)點(diǎn)的矢量姿態(tài)均為[1;0;0]。

    使用原適應(yīng)度函數(shù)和改進(jìn)適應(yīng)度函數(shù)對(duì)選取的空間內(nèi)軌跡方程(16)進(jìn)行連續(xù)軌跡跟蹤,適應(yīng)度函數(shù)和中代表位置誤差的適應(yīng)度函數(shù)的結(jié)果見(jiàn)表5,連續(xù)求解過(guò)程中適應(yīng)度函數(shù)值的波動(dòng)情況如圖4 所示。 仿真結(jié)果表明,改進(jìn)的適應(yīng)度函數(shù)中的平均值比原適應(yīng)度函數(shù)中的平均值少13個(gè)數(shù)量級(jí),同時(shí)對(duì)軌跡上所有點(diǎn)進(jìn)行軌跡跟蹤,軌跡跟蹤以及軌跡仿真如圖5 所示,可以看到能完全實(shí)現(xiàn)對(duì)給定曲線的跟蹤,且運(yùn)動(dòng)過(guò)程平穩(wěn)。

    表5 連續(xù)軌跡跟蹤的適應(yīng)度函數(shù)f1的值Tab.5 The value of the fitness function f1 for continuous trajectory tracking

    圖4 適應(yīng)度函數(shù)f1的波動(dòng)Fig.4 Fluctuation of fitness function f1

    圖5 軌跡仿真與跟蹤Fig.5 Trajectory simulation and tracking

    兩種適應(yīng)度函數(shù)使用CMA-ES 算法求解軌跡的關(guān)節(jié)角度變化如圖6 所示。 結(jié)果表明,在保持原適應(yīng)度函數(shù)的特點(diǎn)上,基于加權(quán)最小二乘法“最佳柔順性”規(guī)則改進(jìn)的適應(yīng)度函數(shù)求解的結(jié)果唯一且更光滑,關(guān)節(jié)運(yùn)動(dòng)距離更小,能源消耗更少,運(yùn)行過(guò)程平穩(wěn)。

    圖6 關(guān)節(jié)角變化Fig.6 Changes of joint angle

    在兩種適應(yīng)度函數(shù)求解下,笛卡爾空間中目標(biāo)位置(,,) 和仿真位置(,,) 之間的誤差變化,如表6、表7 和圖7 表示。 可以看出,兩者求解的位置精度相差較大,采用新適應(yīng)度函數(shù)求解得到的平均位置誤差值穩(wěn)定在10m 數(shù)量級(jí),最小值誤差為0 m,且代表位置誤差穩(wěn)定度的標(biāo)準(zhǔn)差值也穩(wěn)定在10m 數(shù)量級(jí),在軌跡跟蹤精度和軌跡跟蹤穩(wěn)定度方面明顯優(yōu)于使用原適應(yīng)度函數(shù)解得的結(jié)果。

    表6 原適應(yīng)度函數(shù)軌跡跟蹤的位置誤差Tab.6 Position error of original fitness function for trajectory tracking

    表7 改進(jìn)適應(yīng)度函數(shù)軌跡跟蹤的位置誤差Tab.7 Position error of improved fitness function for trajectory tracking

    圖7 空間位置誤差波動(dòng)Fig.7 Fluctuation of spatial position error

    4 結(jié)束語(yǔ)

    本文提出改進(jìn)的適應(yīng)度函數(shù),將加權(quán)最小二乘法的“最佳柔順性”規(guī)則和位姿誤差結(jié)合,形成改進(jìn)的適應(yīng)度函數(shù),使用CMA-ES 算法求逆運(yùn)動(dòng)學(xué)解,使其解具有唯一性;在單點(diǎn)求得的逆解中, 改進(jìn)適應(yīng)度函數(shù)中代表位置誤差的適應(yīng)度函數(shù)的平均值比原適應(yīng)度函數(shù)中的平均值少10 個(gè)數(shù)量級(jí);在連續(xù)的軌跡求解中,改進(jìn)適應(yīng)度函數(shù)中的平均值比原適應(yīng)度函數(shù)逆解的結(jié)果提升13 個(gè)數(shù)量級(jí),且使用改進(jìn)適應(yīng)度函數(shù)求解得到的平均位置誤差值穩(wěn)定在10m。

    基于CMA-ES 算法將機(jī)器人關(guān)節(jié)角范圍作為算法約束條件,改進(jìn)的適應(yīng)度函數(shù)求解精度更高,且解出的各關(guān)節(jié)角位移平滑,路徑更短,遵循“多運(yùn)動(dòng)下臂,少運(yùn)動(dòng)上臂” 的原則,使得機(jī)器人能耗最少,提高了工業(yè)機(jī)器人工作效率。

    猜你喜歡
    運(yùn)動(dòng)學(xué)逆向適應(yīng)度
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    逆向而行
    基于MATLAB的6R機(jī)器人逆運(yùn)動(dòng)學(xué)求解分析
    基于D-H法的5-DOF串并聯(lián)機(jī)床運(yùn)動(dòng)學(xué)分析
    逆向解答
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    基于運(yùn)動(dòng)學(xué)原理的LBI解模糊算法
    雙足機(jī)器人運(yùn)動(dòng)學(xué)分析與仿真
    逆向工程技術(shù)及應(yīng)用
    少數(shù)民族大學(xué)生文化適應(yīng)度調(diào)查
    中文资源天堂在线| 99热这里只有是精品在线观看| 午夜福利视频1000在线观看| 精品久久久久久久久久久久久| 久久精品国产鲁丝片午夜精品| 麻豆乱淫一区二区| 亚洲一级一片aⅴ在线观看| 在线观看一区二区三区| 亚洲18禁久久av| 国产黄色小视频在线观看| 激情 狠狠 欧美| 丰满少妇做爰视频| 蜜桃亚洲精品一区二区三区| 18+在线观看网站| 一级黄色大片毛片| 亚洲国产精品sss在线观看| 国产精品国产三级国产av玫瑰| 国产国拍精品亚洲av在线观看| 欧美日本亚洲视频在线播放| 成年女人永久免费观看视频| 欧美xxxx黑人xx丫x性爽| 精品人妻熟女av久视频| 尤物成人国产欧美一区二区三区| 成人午夜精彩视频在线观看| 精品一区二区三区人妻视频| 日本wwww免费看| 黑人高潮一二区| 欧美一区二区亚洲| 18禁在线无遮挡免费观看视频| 久99久视频精品免费| 欧美丝袜亚洲另类| 午夜福利在线观看吧| 亚洲一级一片aⅴ在线观看| 国产三级在线视频| 午夜日本视频在线| 亚洲精品久久久久久婷婷小说 | 成人毛片60女人毛片免费| 可以在线观看毛片的网站| 国产精品久久久久久久久免| 亚洲av一区综合| 高清日韩中文字幕在线| 久久这里只有精品中国| 精品人妻偷拍中文字幕| 国产免费福利视频在线观看| 老司机影院毛片| 亚洲性久久影院| 国产精品野战在线观看| av在线老鸭窝| 免费黄色在线免费观看| 日韩av在线大香蕉| 亚洲av.av天堂| 26uuu在线亚洲综合色| 五月伊人婷婷丁香| 亚洲真实伦在线观看| 高清午夜精品一区二区三区| 国产又黄又爽又无遮挡在线| 国产精品久久久久久久久免| 国产精品久久久久久久久免| 99热全是精品| 亚洲国产成人一精品久久久| kizo精华| 亚洲精品久久久久久婷婷小说 | 国语对白做爰xxxⅹ性视频网站| 三级毛片av免费| 国产乱人偷精品视频| 麻豆成人av视频| 边亲边吃奶的免费视频| 五月玫瑰六月丁香| 狠狠狠狠99中文字幕| 国产精品久久视频播放| 国产又黄又爽又无遮挡在线| 亚洲电影在线观看av| 国产精品福利在线免费观看| 色哟哟·www| 一区二区三区四区激情视频| a级毛色黄片| 国产探花在线观看一区二区| 国产午夜精品论理片| 国产毛片a区久久久久| 久99久视频精品免费| 九九在线视频观看精品| 亚洲性久久影院| 赤兔流量卡办理| 欧美丝袜亚洲另类| 久久久色成人| 国产在视频线精品| 韩国av在线不卡| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区免费观看| 久久精品国产自在天天线| 久久久久久久午夜电影| 国语对白做爰xxxⅹ性视频网站| 18禁动态无遮挡网站| 级片在线观看| 永久免费av网站大全| 草草在线视频免费看| 日本免费a在线| videossex国产| 男人狂女人下面高潮的视频| av在线亚洲专区| 成人三级黄色视频| 丰满乱子伦码专区| 欧美xxxx性猛交bbbb| 亚洲av成人精品一区久久| 人妻制服诱惑在线中文字幕| 色综合亚洲欧美另类图片| 国产亚洲最大av| 美女cb高潮喷水在线观看| av视频在线观看入口| 能在线免费看毛片的网站| 一卡2卡三卡四卡精品乱码亚洲| av在线亚洲专区| 日本黄色视频三级网站网址| 日日撸夜夜添| 午夜老司机福利剧场| 蜜桃久久精品国产亚洲av| 高清视频免费观看一区二区 | 日韩欧美国产在线观看| 亚洲成人av在线免费| 久久国产乱子免费精品| 天美传媒精品一区二区| 久久精品国产亚洲网站| 中文在线观看免费www的网站| 网址你懂的国产日韩在线| 久久国产乱子免费精品| 一区二区三区四区激情视频| 国产精品久久电影中文字幕| 在线a可以看的网站| 精品久久久久久久久亚洲| 欧美性猛交╳xxx乱大交人| 99在线视频只有这里精品首页| 成人av在线播放网站| 99久久成人亚洲精品观看| 国产毛片a区久久久久| 欧美又色又爽又黄视频| 国产黄a三级三级三级人| 欧美色视频一区免费| 级片在线观看| 中文亚洲av片在线观看爽| 久久久久精品久久久久真实原创| 少妇猛男粗大的猛烈进出视频 | 国产午夜精品久久久久久一区二区三区| 国产精品一区二区三区四区久久| 国产精品美女特级片免费视频播放器| 成人综合一区亚洲| 插逼视频在线观看| 在线免费观看的www视频| 久久午夜福利片| 成年av动漫网址| 日韩在线高清观看一区二区三区| 免费看a级黄色片| 国产视频首页在线观看| 干丝袜人妻中文字幕| 国产成人a区在线观看| 寂寞人妻少妇视频99o| 日日撸夜夜添| 国产片特级美女逼逼视频| av黄色大香蕉| 人人妻人人澡人人爽人人夜夜 | 桃色一区二区三区在线观看| 国产一区有黄有色的免费视频 | 国产精品野战在线观看| 亚洲激情五月婷婷啪啪| 色综合站精品国产| 亚洲欧美日韩高清专用| 国产在线男女| 日韩人妻高清精品专区| 国产高清国产精品国产三级 | 国产成人一区二区在线| 国产三级在线视频| 国产精品嫩草影院av在线观看| 观看免费一级毛片| 日本黄色视频三级网站网址| 国产亚洲91精品色在线| 国产成人aa在线观看| 日韩国内少妇激情av| 人人妻人人澡欧美一区二区| 最近中文字幕高清免费大全6| 超碰97精品在线观看| 国产毛片a区久久久久| 免费不卡的大黄色大毛片视频在线观看 | 又粗又硬又长又爽又黄的视频| 69av精品久久久久久| 成人二区视频| 日韩国内少妇激情av| 国产午夜精品久久久久久一区二区三区| 国产成人精品婷婷| 久久久国产成人免费| 亚洲不卡免费看| 看片在线看免费视频| videos熟女内射| 亚洲图色成人| 亚洲久久久久久中文字幕| 韩国高清视频一区二区三区| 国产亚洲精品久久久com| 最后的刺客免费高清国语| 99久久精品国产国产毛片| 国产人妻一区二区三区在| 久久午夜福利片| 成人av在线播放网站| 干丝袜人妻中文字幕| 亚洲国产成人一精品久久久| av线在线观看网站| 欧美成人午夜免费资源| 看黄色毛片网站| 国产老妇女一区| 欧美日韩综合久久久久久| 亚洲成人中文字幕在线播放| 一夜夜www| 国产精品国产三级国产专区5o | 日韩精品有码人妻一区| 欧美成人a在线观看| 精品久久久久久成人av| 不卡视频在线观看欧美| 欧美日韩精品成人综合77777| 一个人免费在线观看电影| 干丝袜人妻中文字幕| 日韩,欧美,国产一区二区三区 | 久热久热在线精品观看| 三级男女做爰猛烈吃奶摸视频| 国产伦精品一区二区三区视频9| 天美传媒精品一区二区| 成人欧美大片| 日韩av不卡免费在线播放| 99九九线精品视频在线观看视频| 老司机福利观看| 夫妻性生交免费视频一级片| 国产黄a三级三级三级人| 能在线免费观看的黄片| 网址你懂的国产日韩在线| 国产私拍福利视频在线观看| 三级毛片av免费| 国产成人一区二区在线| 自拍偷自拍亚洲精品老妇| 免费观看性生交大片5| 国产在视频线精品| 精品国内亚洲2022精品成人| 波多野结衣高清无吗| 淫秽高清视频在线观看| 精品久久久久久成人av| 神马国产精品三级电影在线观看| 女人十人毛片免费观看3o分钟| 国产精品永久免费网站| 网址你懂的国产日韩在线| 久久久久久久久大av| 99久国产av精品国产电影| 国产欧美日韩精品一区二区| 久久精品久久久久久噜噜老黄 | 亚洲成av人片在线播放无| 亚洲自拍偷在线| 亚洲欧美日韩卡通动漫| 真实男女啪啪啪动态图| 久久久午夜欧美精品| 日韩av在线大香蕉| 老女人水多毛片| 国产免费一级a男人的天堂| 国产亚洲最大av| 晚上一个人看的免费电影| 成人亚洲欧美一区二区av| 国产精品人妻久久久影院| 性色avwww在线观看| 日韩成人伦理影院| 亚洲av男天堂| 久久精品夜色国产| 老司机影院毛片| 亚洲精品色激情综合| 99久国产av精品| 99热精品在线国产| 国产一区二区在线av高清观看| 在线观看一区二区三区| 国产男人的电影天堂91| 亚洲av熟女| 在线免费观看不下载黄p国产| 99久久人妻综合| 特级一级黄色大片| 国产成人精品久久久久久| 久久99蜜桃精品久久| 天堂√8在线中文| 精品不卡国产一区二区三区| 精品一区二区三区人妻视频| av在线蜜桃| 久久久久久久久大av| 少妇猛男粗大的猛烈进出视频 | videos熟女内射| 精品人妻一区二区三区麻豆| 免费无遮挡裸体视频| 亚洲欧美精品专区久久| 国产精品精品国产色婷婷| 精品久久久久久久末码| 欧美成人免费av一区二区三区| 人人妻人人看人人澡| 两个人的视频大全免费| 日日撸夜夜添| 嫩草影院入口| 91精品国产九色| 国产成人一区二区在线| 亚洲av日韩在线播放| 男人舔奶头视频| 国产一区二区在线av高清观看| 国内精品宾馆在线| 男女那种视频在线观看| 色5月婷婷丁香| 一区二区三区高清视频在线| 波多野结衣高清无吗| 国产成人福利小说| 免费观看的影片在线观看| av黄色大香蕉| 白带黄色成豆腐渣| 午夜a级毛片| 99国产精品一区二区蜜桃av| 国产精品伦人一区二区| 在线播放无遮挡| 精品酒店卫生间| 日本与韩国留学比较| 91久久精品国产一区二区三区| 大香蕉97超碰在线| 九九爱精品视频在线观看| 少妇高潮的动态图| 亚洲精品日韩av片在线观看| 波多野结衣巨乳人妻| 日韩一区二区视频免费看| 亚洲中文字幕日韩| 久久久久久久久久黄片| 久久6这里有精品| 美女国产视频在线观看| 中文乱码字字幕精品一区二区三区 | 三级经典国产精品| 亚洲成人精品中文字幕电影| 淫秽高清视频在线观看| 午夜激情福利司机影院| 亚洲成人av在线免费| 亚洲av男天堂| 一区二区三区乱码不卡18| 级片在线观看| 国产老妇女一区| 一区二区三区免费毛片| 夜夜爽夜夜爽视频| 欧美三级亚洲精品| 高清日韩中文字幕在线| 日韩精品有码人妻一区| 国产精品熟女久久久久浪| 国产av一区在线观看免费| 国产女主播在线喷水免费视频网站 | 成人无遮挡网站| 精品久久久久久久久久久久久| 白带黄色成豆腐渣| 国内揄拍国产精品人妻在线| 天美传媒精品一区二区| 免费看日本二区| 亚洲在久久综合| 亚洲真实伦在线观看| 美女xxoo啪啪120秒动态图| 亚洲欧美中文字幕日韩二区| 国产精品永久免费网站| 亚洲,欧美,日韩| 成人三级黄色视频| 国产真实伦视频高清在线观看| 中文亚洲av片在线观看爽| 亚洲欧美成人综合另类久久久 | 免费av不卡在线播放| 一区二区三区高清视频在线| 一个人免费在线观看电影| 麻豆乱淫一区二区| 搞女人的毛片| 日韩精品有码人妻一区| 我要搜黄色片| 精品一区二区三区人妻视频| 亚洲综合色惰| 亚洲高清免费不卡视频| 亚洲av成人精品一区久久| 亚洲av二区三区四区| 亚洲精品久久久久久婷婷小说 | 美女国产视频在线观看| 亚洲最大成人手机在线| 天堂网av新在线| 欧美高清成人免费视频www| av在线观看视频网站免费| 在线观看美女被高潮喷水网站| 国产av一区在线观看免费| 欧美bdsm另类| 一级黄片播放器| 黄色日韩在线| 校园人妻丝袜中文字幕| 在线播放国产精品三级| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧洲日产国产| 老司机福利观看| 永久免费av网站大全| 校园人妻丝袜中文字幕| 亚洲精品456在线播放app| 国产探花极品一区二区| 免费播放大片免费观看视频在线观看 | 久久这里有精品视频免费| 国内少妇人妻偷人精品xxx网站| 久久久精品94久久精品| 精品无人区乱码1区二区| 免费av观看视频| ponron亚洲| 视频中文字幕在线观看| av黄色大香蕉| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 午夜激情欧美在线| 国产成人午夜福利电影在线观看| 人体艺术视频欧美日本| 精品久久久久久电影网 | 青春草亚洲视频在线观看| 亚洲国产精品成人久久小说| 日韩视频在线欧美| 久久婷婷人人爽人人干人人爱| 亚洲精品影视一区二区三区av| 久久久a久久爽久久v久久| 免费观看在线日韩| 久久精品久久精品一区二区三区| 超碰97精品在线观看| 久久久久国产网址| av播播在线观看一区| 亚洲综合精品二区| 成人欧美大片| 国产精品福利在线免费观看| 一夜夜www| 成人二区视频| 国产精品乱码一区二三区的特点| 日韩精品有码人妻一区| 青春草亚洲视频在线观看| 韩国av在线不卡| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧洲国产日韩| 免费大片18禁| 色视频www国产| 国产伦理片在线播放av一区| 精品久久久久久久久亚洲| 一边摸一边抽搐一进一小说| 久久久久久久亚洲中文字幕| 欧美又色又爽又黄视频| 久久久久久久久大av| 国产精品久久久久久精品电影小说 | 成人漫画全彩无遮挡| 99久久精品一区二区三区| 亚洲一级一片aⅴ在线观看| 国产久久久一区二区三区| 2022亚洲国产成人精品| 免费看日本二区| 国产三级在线视频| 少妇熟女欧美另类| av播播在线观看一区| 麻豆乱淫一区二区| 在线观看av片永久免费下载| 韩国av在线不卡| 午夜福利在线在线| 在现免费观看毛片| 身体一侧抽搐| 一本一本综合久久| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| 国产成人精品久久久久久| 久久久久九九精品影院| 免费无遮挡裸体视频| av在线播放精品| 特级一级黄色大片| 综合色av麻豆| 在线播放无遮挡| 最近的中文字幕免费完整| 熟女电影av网| 成人无遮挡网站| 一个人看视频在线观看www免费| 97在线视频观看| 国产成人一区二区在线| 国产精品精品国产色婷婷| 亚洲人成网站在线观看播放| 九色成人免费人妻av| 欧美成人精品欧美一级黄| 97人妻精品一区二区三区麻豆| 99久久人妻综合| 国内揄拍国产精品人妻在线| 内地一区二区视频在线| 26uuu在线亚洲综合色| 中文在线观看免费www的网站| 欧美日韩精品成人综合77777| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区在线观看99 | 国内精品一区二区在线观看| 有码 亚洲区| 青春草视频在线免费观看| 五月玫瑰六月丁香| 免费观看精品视频网站| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 91在线精品国自产拍蜜月| 欧美激情国产日韩精品一区| 一级二级三级毛片免费看| 国产伦在线观看视频一区| 欧美一区二区亚洲| 日日摸夜夜添夜夜添av毛片| 久久久久久国产a免费观看| 欧美日韩国产亚洲二区| 毛片一级片免费看久久久久| 成年版毛片免费区| 少妇猛男粗大的猛烈进出视频 | 成人特级av手机在线观看| 一个人看视频在线观看www免费| 69人妻影院| www.av在线官网国产| 99国产精品一区二区蜜桃av| 国产免费男女视频| 久久久久九九精品影院| 天堂中文最新版在线下载 | 村上凉子中文字幕在线| 熟女人妻精品中文字幕| 亚洲真实伦在线观看| 国产69精品久久久久777片| 国产精品无大码| 麻豆av噜噜一区二区三区| 色综合亚洲欧美另类图片| 国产av一区在线观看免费| 欧美日韩精品成人综合77777| 久99久视频精品免费| 日本色播在线视频| 欧美另类亚洲清纯唯美| 色噜噜av男人的天堂激情| 亚洲人成网站在线播| 国产精品,欧美在线| 国产亚洲av片在线观看秒播厂 | 午夜久久久久精精品| 26uuu在线亚洲综合色| 午夜激情福利司机影院| 热99在线观看视频| 老师上课跳d突然被开到最大视频| 国产亚洲午夜精品一区二区久久 | 国模一区二区三区四区视频| 久久久亚洲精品成人影院| 亚洲av不卡在线观看| 精品人妻一区二区三区麻豆| 26uuu在线亚洲综合色| 嘟嘟电影网在线观看| 黄片无遮挡物在线观看| 在线免费观看不下载黄p国产| 又粗又爽又猛毛片免费看| 成人av在线播放网站| 欧美三级亚洲精品| 春色校园在线视频观看| 国产毛片a区久久久久| 免费搜索国产男女视频| 少妇熟女aⅴ在线视频| 最近中文字幕高清免费大全6| 小蜜桃在线观看免费完整版高清| 床上黄色一级片| 精品国内亚洲2022精品成人| 亚洲国产欧美在线一区| 中文字幕久久专区| 我要看日韩黄色一级片| 国产老妇女一区| 韩国高清视频一区二区三区| 最近最新中文字幕免费大全7| 伦精品一区二区三区| 日韩高清综合在线| 亚洲欧美精品专区久久| 真实男女啪啪啪动态图| 18禁在线播放成人免费| 国产亚洲精品av在线| av视频在线观看入口| 两性午夜刺激爽爽歪歪视频在线观看| 在线a可以看的网站| 国模一区二区三区四区视频| 中文天堂在线官网| 国产乱来视频区| 日韩欧美精品v在线| 国产久久久一区二区三区| 热99在线观看视频| 国产 一区 欧美 日韩| 亚洲欧洲国产日韩| 99热这里只有精品一区| 久久欧美精品欧美久久欧美| 久久精品国产亚洲网站| 自拍偷自拍亚洲精品老妇| 国产69精品久久久久777片| 亚洲欧美日韩无卡精品| 麻豆av噜噜一区二区三区| 亚洲精品乱码久久久久久按摩| 国产视频内射| 只有这里有精品99| 又爽又黄无遮挡网站| 日韩精品青青久久久久久| 免费黄色在线免费观看| 日本免费一区二区三区高清不卡| 亚洲av中文av极速乱| 建设人人有责人人尽责人人享有的 | av专区在线播放| 91精品伊人久久大香线蕉| 亚洲精品久久久久久婷婷小说 | 久久亚洲精品不卡| 欧美日本视频| 欧美色视频一区免费| 久久久久精品久久久久真实原创| 免费av观看视频| 99久久精品国产国产毛片| 男人舔女人下体高潮全视频| 国产激情偷乱视频一区二区| 午夜老司机福利剧场| 国产午夜精品久久久久久一区二区三区| 国产探花在线观看一区二区| 久久午夜福利片| 中文欧美无线码| 久久这里有精品视频免费| 一本久久精品| 26uuu在线亚洲综合色| 国产午夜精品一二区理论片| 18禁裸乳无遮挡免费网站照片| 亚洲精品日韩av片在线观看| 秋霞伦理黄片| 亚洲欧美清纯卡通| 最新中文字幕久久久久| 人人妻人人看人人澡| 精品欧美国产一区二区三| 少妇被粗大猛烈的视频|