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

    改進(jìn)多目標(biāo)PSO算法優(yōu)化機(jī)器人軌跡跟蹤模糊PID控制器

    2021-07-03 05:45:54蔣清澤王宏濤
    應(yīng)用科技 2021年3期
    關(guān)鍵詞:軌跡向量粒子

    蔣清澤,王宏濤

    南京航空航天大學(xué) 機(jī)電學(xué)院,江蘇 南京 210016

    關(guān)節(jié)機(jī)器人是一個(gè)強(qiáng)耦合、高非線性的復(fù)雜系統(tǒng),其軌跡跟蹤控制是調(diào)整每個(gè)關(guān)節(jié)電機(jī)輸出的控制扭矩,使得關(guān)節(jié)角度能夠達(dá)到期望值[1]。關(guān)節(jié)機(jī)器人軌跡跟蹤常用的控制方法主要有基于機(jī)器人模型的運(yùn)動(dòng)控制方法、PID控制、滑模控制和模糊控制等。其中PID控制技術(shù)成熟且易于實(shí)現(xiàn),被廣泛應(yīng)用于機(jī)器人軌跡跟蹤控制中,但PID算法難以實(shí)現(xiàn)非線性系統(tǒng)的高精度控制[2]。因此,有學(xué)者提出利用模糊控制無(wú)須依靠精確數(shù)學(xué)模型便能對(duì)非線性系統(tǒng)良好逼近的特性,采用模糊PID控制實(shí)現(xiàn)機(jī)器人軌跡跟蹤控制[3]。然而,在實(shí)際應(yīng)用中,控制器輸出力矩具有飽和性,如何設(shè)計(jì)模糊PID控制器使其在有限的輸出力矩下達(dá)到理想的控制效果是需要專業(yè)人員的經(jīng)驗(yàn),存在設(shè)計(jì)不當(dāng)難以獲得最佳效果的不足。均衡評(píng)估控制力矩和軌跡跟蹤控制效果這2個(gè)優(yōu)化目標(biāo)可以通過多目標(biāo)優(yōu)化算法對(duì)模糊PID控制的隸屬函數(shù)和模糊規(guī)則進(jìn)行尋優(yōu)調(diào)整來(lái)實(shí)現(xiàn)[4?5]。粒子群算法(PSO)是一種模擬鳥群覓食的啟發(fā)式算法優(yōu)化技術(shù),已成功應(yīng)用于多種優(yōu)化問題[6]。已有文獻(xiàn)證明PSO在單目標(biāo)優(yōu)化上具有較快的收斂速度和較好的優(yōu)化結(jié)果,在多目標(biāo)優(yōu)化問題上具有較大優(yōu)勢(shì)[7]。Mahmoodabadi等[8]用多目標(biāo)粒子群算法對(duì)非線性系統(tǒng)的模糊控制器進(jìn)行優(yōu)化,提高了模糊控制器對(duì)非線性系統(tǒng)的控制精度;進(jìn)一步以控制偏差和控制器輸出力矩作為2個(gè)優(yōu)化目標(biāo),采用多目標(biāo)粒子群算法對(duì)雙足機(jī)器人軌跡跟蹤滑??刂七M(jìn)行了優(yōu)化,驗(yàn)證了多目標(biāo)粒子群算法在多目標(biāo)優(yōu)化中的優(yōu)越性[9]。

    本文以控制器輸出力矩和軌跡跟蹤偏差為2個(gè)優(yōu)化目標(biāo),實(shí)現(xiàn)改進(jìn)多目標(biāo)PSO算法優(yōu)化模糊控制器的隸屬函數(shù)和模糊規(guī)則優(yōu)化調(diào)整。首先設(shè)計(jì)了機(jī)器人軌跡跟蹤模糊PID控制器;然后根據(jù)模糊PID控制器的多目標(biāo)優(yōu)化問題,設(shè)計(jì)改進(jìn)多目標(biāo)PSO算法;最后分別使用改進(jìn)多目標(biāo)PSO和基本多目標(biāo)PSO應(yīng)用于模糊PID控制的多目標(biāo)優(yōu)化問題,對(duì)比分析了2種優(yōu)化算法的優(yōu)化結(jié)果。

    1 機(jī)器人軌跡跟蹤模糊PID控制

    機(jī)器人軌跡跟蹤模糊PID控制系統(tǒng)如圖1所示,圖中qd和q分別表示六關(guān)節(jié)機(jī)器人各關(guān)節(jié)期望角度和實(shí)際角度;e表示qd與q的偏差,e=qd?q;ec表示偏差變化率,ec=de/dt;圖1虛線部分表示模糊控制器,其中e和ec是模糊控制器的輸入變量,ΔKP、ΔKI、ΔKD表 示PID調(diào)整參數(shù),是模糊控制器的輸出變量;τ表示模糊PID控制器的輸出控制力矩??刂七^程為:首先模糊控制器根據(jù)e和ec的值通過模糊化、模糊推理和解模糊得到PID調(diào)整參數(shù);然后實(shí)現(xiàn)PID控制器3個(gè)參數(shù)的在線自適應(yīng)調(diào)整,如式(1),根據(jù)式(2)輸出控制力矩;最后將控制力矩作用于機(jī)器人,得到實(shí)際關(guān)節(jié)角度,并反饋給比較器。

    圖1 機(jī)器人軌跡跟蹤模糊PID控制系統(tǒng)

    式中:KP0、KI0、KD0為 PID參數(shù)的初值;t為時(shí)間。

    1.1 模糊PID控制器設(shè)計(jì)

    模糊PID控制器中的模糊控制器輸入變量偏差e和偏差變化率ec以及輸出變量ΔKP、ΔKI、ΔKD的 模糊集個(gè)數(shù)均設(shè)計(jì)為7個(gè),其模糊集取為{NL,NM, NS, ZE, PS, PM, PL},即{負(fù)大, 負(fù)中, 負(fù)小, 零,正小, 正中, 正大}。模糊集隸屬函數(shù)選取三角形隸屬函數(shù),根據(jù)隸屬函數(shù)的特性,某變量所有的模糊集的隸屬函數(shù)關(guān)于Y軸對(duì)稱,且將某一個(gè)三角隸屬函數(shù)的“左腳”和“右腳”分別設(shè)置為相鄰左右2個(gè)三角隸屬函數(shù)的頂點(diǎn)橫坐標(biāo),如圖2所示。故在設(shè)計(jì)每一個(gè)變量的隸屬函數(shù)時(shí),需要確定其論域,以及Y軸某一側(cè)3個(gè)參數(shù)X1、X2、X3。模糊PID控制器所有變量的隸屬函數(shù)具體參數(shù)如表1所示。

    圖2 三角形隸屬函數(shù)

    表1 模糊PID控制器輸入輸出變量隸屬函數(shù)設(shè)計(jì)

    2個(gè)輸入變量e和ec的模糊集個(gè)數(shù)為7,可得輸出變量ΔKP、ΔKI、ΔKD的每一個(gè)變量所對(duì)應(yīng)的模糊規(guī)則均為49條,每條模糊規(guī)則的形式為

    本文設(shè)計(jì)的ΔKP、ΔKI、ΔKD的 模糊規(guī)則如表2所示,數(shù)字1~7分別代表NL、NM、NS、ZE、PS、PM、PL。

    表2 ΔKP/Δ KI/Δ KD的模糊規(guī)則表

    1.2 機(jī)器人動(dòng)力學(xué)模型

    忽略關(guān)節(jié)摩擦和末端負(fù)載,六關(guān)節(jié)機(jī)器人動(dòng)力學(xué)方程可以寫成如下形式

    為了對(duì)六關(guān)節(jié)機(jī)器人進(jìn)行仿真,根據(jù)上述建立的動(dòng)力學(xué)方程式(3)可計(jì)算方程中關(guān)節(jié)角加速度

    由t時(shí)刻的關(guān)節(jié)角加速度,假設(shè)時(shí)間步長(zhǎng)Δt內(nèi)不變化, 在初始條件為的情況下,對(duì)當(dāng)前t時(shí)刻關(guān)節(jié)信息qt和進(jìn)行數(shù)值積分, 計(jì)算出下一時(shí)刻t+Δt的角速度和關(guān)節(jié)角qt+Δt為

    2 模糊PID控制的多目標(biāo)PSO優(yōu)化

    2.1 模糊PID控制的多目標(biāo)問題

    優(yōu)化機(jī)器人軌跡跟蹤模糊PID控制,需要將待優(yōu)化參數(shù)進(jìn)行編碼形成一個(gè)決策向量。粒子群算法中的決策向量為某個(gè)粒子i的位置向量xi,具體表達(dá)為

    式中:D為優(yōu)化問題中待優(yōu)化參數(shù)的個(gè)數(shù)。

    PSO優(yōu)化模糊PID控制的粒子位置向量xi中的元素所代表的參數(shù)及參數(shù)取值范圍如表3所示,其中元素xi,1?xi,15分別表示輸入變量e、ec和輸出變量ΔKP、ΔKI、ΔKD的隸屬函數(shù)中心頂點(diǎn)值,每個(gè)變量需指定3個(gè)中心頂點(diǎn)值,故5個(gè)變量需要15個(gè)表示模糊PID控制器隸屬函數(shù)的設(shè)計(jì)參數(shù);模糊PID控制器總共有3個(gè)規(guī)則表,每個(gè)規(guī)則表有49條模糊規(guī)則結(jié)果需要設(shè)計(jì),故需要3×49=147個(gè)元素表示,即xi,16?xi,64、xi,65?xi,113、xi,114?xi,162分別表示ΔKP、ΔKI、ΔKD規(guī)則的49條模糊規(guī)則。故總共需要162位元素基來(lái)表示,D=162。

    表3 粒子位置向量參數(shù)

    模糊PID控制有2個(gè)優(yōu)化目標(biāo)函數(shù),分別與軌跡跟蹤偏差和控制器輸出力矩有關(guān),如式(7)所示,則目標(biāo)向量為f=(f1,f2)。

    式中:e(t)為t時(shí)刻的機(jī)器人笛卡爾空間仿真軌跡與期望軌跡的偏差;τn(t)為t時(shí) 刻模糊PID控制器作用在第n個(gè) 關(guān)節(jié)的輸出力矩;T為運(yùn)行模糊PID控制器完成一次軌跡跟蹤控制仿真的總時(shí)間。

    模糊PID控制器多目標(biāo)優(yōu)化時(shí),最小化控制器輸出力矩可能導(dǎo)致軌跡跟蹤誤差的提高,因此需要在這2個(gè)互相沖突的目標(biāo)中找到權(quán)衡的解,而這樣的解不唯一。2個(gè)粒子的位置向量x1和x2,若對(duì)于?k=1,2,有fk(x1)≤fk(x2),且對(duì)于?k=1,2有fk(x1)<fk(x2),稱x1支配x2;否則稱x1與x2為非支配關(guān)系。多目標(biāo)優(yōu)化算法的任務(wù)是找到所有不被其他粒子位置向量支配的粒子所對(duì)應(yīng)的目標(biāo)向量f=(f1,f2)的集合,這些目標(biāo)向量互為非支配關(guān)系,其集合稱為非支配解集。

    2.2 改進(jìn)多目標(biāo)PSO算法

    粒子群算法是通過各個(gè)粒子位置的不斷調(diào)整,使得粒子不斷向最優(yōu)位置移動(dòng)。算法中每個(gè)粒子均具有形如[xi,1(k),xi,2(k),···,xi,D(k)]的位置向量xi(k)和形如[vi,1(k),vi,2(k),···,vi,D(k)]的速度向量vi(k),其中,i=1,2,···,S是粒子在粒子群中的序號(hào),k為迭代次數(shù),S是粒子群的個(gè)數(shù)。

    改進(jìn)多目標(biāo)PSO算法優(yōu)化機(jī)器人軌跡跟蹤控制模糊PID控制器的流程如圖3所示,算法主要步驟如下。

    圖3 多目標(biāo)粒子群算法優(yōu)化模糊PID控制器流程

    1)生成初始粒子種群:根據(jù)表3,隨機(jī)生成S個(gè)形如[xi,1(0),xi,2(0),···,xi,162(0)]的粒子位置向量xi(0),且初始種群為{x1(0),x2(0),···,xS(0)},設(shè)置初始種群中所有粒子的速度向量vi(0)=0。

    2)更新粒子群:設(shè)置第k代粒子群為當(dāng)前粒子群。

    3)更新隸屬函數(shù)和模糊規(guī)則:讀取當(dāng)前粒子群中第j(j=1,2,···,S)個(gè)粒子的位置向量,根據(jù)編碼xi,1(k)?xi,10(k)生成新的輸入變量e、ec和輸出變量ΔKP、ΔKI、ΔKD的隸屬函數(shù),更新隸屬函數(shù);根據(jù)編碼xi,11(k)?xi,63(k)產(chǎn)生新的輸出變量ΔKP、ΔKI、ΔKD的 模糊規(guī)則,更新模糊規(guī)則。

    4)計(jì)算目標(biāo)函數(shù)值:運(yùn)行模糊PID控制器完成軌跡跟蹤控制仿真,獲得仿真控制中t時(shí)刻機(jī)器人仿真軌跡與期望軌跡的偏差e(t)和各關(guān)節(jié)控制器的輸出ui(t),根據(jù)式(7)計(jì)算目標(biāo)向量f=(f1,f2)。

    5)判斷是否遍歷S個(gè)粒子:若粒子數(shù)量j不滿足j>S,則j+1,進(jìn)入步驟3);若滿足j>S,則進(jìn)入步驟6)。

    6)更新個(gè)體檔案和全局檔案:個(gè)體檔案和全局檔案分別是存儲(chǔ)個(gè)體最佳位置和全局最優(yōu)粒子的非支配解集。每個(gè)粒子均有自己的個(gè)體檔案,而全局檔案為所有粒子共享。個(gè)體檔案和全局檔案是為步驟8中選擇pbest和gbest做準(zhǔn)備。若當(dāng)代粒子群中的粒子與檔案中的解集為非支配關(guān)系,則將該粒子添加到檔案中。

    7)判斷是否終止優(yōu)化算法:判斷粒子迭代次數(shù)k是否達(dá)到設(shè)定值,若到達(dá)最大值,則結(jié)束優(yōu)化算法,返回當(dāng)前全局檔案中的粒子位置向量;若未到達(dá),則進(jìn)入步驟8)。

    8)選擇pbest和gbest:需從個(gè)體檔案和全局檔案中的非支配解集中分別選擇出來(lái)一個(gè)確定的值。選擇pbest和gbest的具體方法如下,首先對(duì)于每一個(gè)粒子均被賦予一個(gè)值σi,且對(duì)于個(gè)人檔案和全局檔案中的每個(gè)粒子也被賦予一個(gè)值σpi和σgi,目標(biāo)函數(shù)值為f=(f1,f2)的粒子值σ被定義為

    然后將個(gè)體檔案和全局檔案中粒子σpi和σgi最接近粒子i的σi的檔案成員作為該粒子的向?qū)best和gbest。

    9)計(jì)算速度向量所需參數(shù):速度向量所需參數(shù)分別為第k代粒子的慣性權(quán)重ω(k)、正加速度系數(shù)c1(k)和c2(k),其具體計(jì)算公式如式(9)—式(11)所示

    式中:nk為算法設(shè)置的最大迭代次數(shù);ω(0)和ω(nk)分別為初始慣性權(quán)重和最終慣性權(quán)重;ω(k)為第k代的慣性權(quán)重,且ω(0)≥ω(nk)。

    式中:c1(0)和c1(nk)分別為c1的初值與終值;c2(0)和c2(nk)分別為c2的初值與終值,且c1(0)≥c1(nk),c2(0)≤c2(nk)。

    10)計(jì)算粒子速度向量:下一次迭代的速度向量vi(k+1)中第d維速度計(jì)算如式(12)所示

    式中:Vmax,d為第d維上允許的最大速度;計(jì)算如式(13)所示

    式中:xi,d(k)和vi,d(k)分別為在k時(shí)刻粒子i第d維上的位置和速度,d=1,2,···,D;ω為慣性權(quán)重系數(shù);pbesti,d(k)為粒子i從個(gè)體檔案中選擇出來(lái)的pbest的第d維,gbestd(k)為在k時(shí)刻當(dāng)前全局檔案中選擇出來(lái)的gbest的第d維;系數(shù)r1、r2則為在區(qū)間[0,1]上服從均勻分布U(0,1)的隨機(jī)數(shù)。

    11)更新粒子位置:根據(jù)式(14)更新粒子位置向量xi(k+1):

    12)粒子位置變異:根據(jù)式(15)對(duì)新粒子群中的粒子位置進(jìn)行變異,進(jìn)入步驟2)。

    式中:xi,d為粒子i在第d維上的值;xd,max和xd,min分別為粒子在第d維上值的上界與下界;系數(shù)r為在區(qū)間[0,1]上服從均勻分布U(0,1)的隨機(jī)數(shù)。

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

    在機(jī)器人軌跡跟蹤控制仿真中,設(shè)笛卡爾空間的期望軌跡為一段螺旋線,其參數(shù)方程如下

    式中:t為時(shí)間;Tc為周期時(shí)長(zhǎng);x、y、z分別為笛卡爾空間的坐標(biāo)值。

    六關(guān)節(jié)機(jī)器人的D-H參數(shù)如表4所示。根據(jù)機(jī)器人運(yùn)動(dòng)學(xué)可將笛卡爾空間的期望軌跡轉(zhuǎn)換成機(jī)器人關(guān)節(jié)空間各個(gè)關(guān)節(jié)的期望角度軌跡,如圖4所示。本文搭建的機(jī)器人軌跡跟蹤模糊PID控制的MATLAB仿真系統(tǒng)如圖5所示,仿真步長(zhǎng)Δt=0.01 s,仿真總時(shí)長(zhǎng)T=10 s。運(yùn)用該控制器可以得到各個(gè)時(shí)刻各關(guān)節(jié)角度的仿真值和輸出力矩值,用于計(jì)算優(yōu)化算法中式(7)的2個(gè)目標(biāo)函數(shù)值f1、f2,按照?qǐng)D3的流程圖,設(shè)置表5的優(yōu)化算法參數(shù),運(yùn)行上述設(shè)計(jì)的控制仿真系統(tǒng),完成關(guān)節(jié)機(jī)器人軌跡跟蹤模糊PID控制器的改進(jìn)多目標(biāo)PSO優(yōu)化。

    表4 關(guān)節(jié)機(jī)器人D-H參數(shù)表

    圖4 關(guān)節(jié)空間期望角度軌跡

    圖5 機(jī)器人軌跡跟蹤模糊PID控制的MATLAB仿真系統(tǒng)

    表5 NSGA II優(yōu)化算法參數(shù)

    分別采用帶精英策略的非支配排序遺傳算法(NSGA II)、多目標(biāo)基本粒子群算法(MOPSO)和改進(jìn)多目標(biāo)粒子群算法(Modify MOPSO)進(jìn)行軌跡跟蹤模糊PID控制器優(yōu)化。NSGA II優(yōu)化算法的參數(shù)設(shè)置如表5所示,2種粒子群優(yōu)化算法的參數(shù)設(shè)置如表6所示。通過3種優(yōu)化算法優(yōu)化后得到的非支配解集如圖6所示,解集中各個(gè)解的2個(gè)目標(biāo)函數(shù)值如表7所示。從圖6可以看出,Modify PSO的解集分布較廣且均勻,而通過MOPSO得到的部分解集在某一個(gè)區(qū)域比較集中,且涵蓋的范圍較小,由NSGA II所得到的解集雖然分布較廣,但分布不均勻。采用能夠同時(shí)評(píng)價(jià)解集收斂性和多樣性的超體積指標(biāo)來(lái)比較3種優(yōu)化算法獲得的非支配解集,超體積越大說(shuō)明算法綜合性能越好。設(shè)置超體積的計(jì)算參考點(diǎn)為(36.5,14 900),NSGA II、MOPSO和Modify MOPSO的超體積分別為885.4、3 369.1和7 247.6,可見Modify PSO獲得的非支配解集最優(yōu)。

    表6 MOPSO和Modify MOPSO優(yōu)化算法參數(shù)

    表7 目標(biāo)函數(shù)值

    圖6 多目標(biāo)優(yōu)化結(jié)果

    綜上,改進(jìn)多目標(biāo)PSO算法應(yīng)用于關(guān)節(jié)機(jī)器人軌跡跟蹤模糊PID控制器的多目標(biāo)優(yōu)化問題,可以得到比基本多目標(biāo)PSO算法和NSGA II更優(yōu)的非支配解集,驗(yàn)證了該算法在機(jī)器人軌跡跟蹤模糊PID控制優(yōu)化上有效性和優(yōu)越性。

    4 結(jié)論

    本文設(shè)計(jì)了一個(gè)模糊PID控制器應(yīng)用于六關(guān)節(jié)機(jī)器人軌跡跟蹤控制,以控制器輸出力矩和軌跡跟蹤偏差為2個(gè)優(yōu)化目標(biāo),實(shí)現(xiàn)了改進(jìn)多目標(biāo)PSO算法優(yōu)化模糊控制器的隸屬函數(shù)和模糊規(guī)則優(yōu)化調(diào)整。對(duì)比改進(jìn)多目標(biāo)PSO和基本多目標(biāo)PSO算法、遺傳算法應(yīng)用于關(guān)節(jié)機(jī)器人軌跡跟蹤模糊PID控制器的多目標(biāo)問題,結(jié)果表明,改進(jìn)PSO獲得的非支配解集更好,驗(yàn)證了該算法應(yīng)用于機(jī)器人軌跡跟蹤模糊PID控制優(yōu)化的有效性和優(yōu)越性。進(jìn)一步的研究工作是繼續(xù)探索多目標(biāo)PSO算法與其他優(yōu)化算法(如蟻群算法、神經(jīng)網(wǎng)絡(luò)算法等)實(shí)現(xiàn)機(jī)器人軌跡跟蹤控制器優(yōu)化的效果,進(jìn)行對(duì)比分析。

    猜你喜歡
    軌跡向量粒子
    向量的分解
    聚焦“向量與三角”創(chuàng)新題
    軌跡
    軌跡
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    軌跡
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    進(jìn)化的軌跡(一)——進(jìn)化,無(wú)盡的適應(yīng)
    向量垂直在解析幾何中的應(yīng)用
    向量五種“變身” 玩轉(zhuǎn)圓錐曲線
    欧美黑人精品巨大| 一边摸一边做爽爽视频免费| 脱女人内裤的视频| 99国产精品一区二区三区| 亚洲精品第二区| 亚洲精品中文字幕在线视频| 久久狼人影院| 视频区欧美日本亚洲| 丝袜美腿诱惑在线| 欧美老熟妇乱子伦牲交| 亚洲三区欧美一区| 一级毛片女人18水好多| 99香蕉大伊视频| 亚洲人成77777在线视频| 精品久久久久久久毛片微露脸 | 国产一区有黄有色的免费视频| 黄色 视频免费看| 国产精品久久久久久精品古装| 999精品在线视频| 国产精品1区2区在线观看. | 久久天躁狠狠躁夜夜2o2o| 欧美日韩视频精品一区| 菩萨蛮人人尽说江南好唐韦庄| 一本一本久久a久久精品综合妖精| 男女午夜视频在线观看| 99国产精品一区二区三区| 男人爽女人下面视频在线观看| 亚洲精品在线美女| 考比视频在线观看| 久久久久久久久久久久大奶| 在线av久久热| 91精品国产国语对白视频| 国产高清videossex| 久久精品国产亚洲av高清一级| 亚洲激情五月婷婷啪啪| 免费不卡黄色视频| 亚洲精品久久午夜乱码| 中文字幕制服av| 日韩精品免费视频一区二区三区| 91av网站免费观看| 最黄视频免费看| 久久99热这里只频精品6学生| 热99国产精品久久久久久7| 999精品在线视频| 国产淫语在线视频| 美女脱内裤让男人舔精品视频| 操出白浆在线播放| 久久综合国产亚洲精品| 国产免费现黄频在线看| 国产精品一区二区在线观看99| 久久九九热精品免费| 免费在线观看黄色视频的| 啦啦啦 在线观看视频| 国产无遮挡羞羞视频在线观看| 午夜免费鲁丝| 国产男女内射视频| 成人国产av品久久久| 国产老妇伦熟女老妇高清| 91精品伊人久久大香线蕉| 9色porny在线观看| 亚洲国产中文字幕在线视频| 亚洲免费av在线视频| 最黄视频免费看| 如日韩欧美国产精品一区二区三区| 超色免费av| 男女边摸边吃奶| 少妇被粗大的猛进出69影院| 亚洲中文字幕日韩| 日本91视频免费播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美成人综合另类久久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美清纯卡通| 精品亚洲乱码少妇综合久久| 国产欧美亚洲国产| 色精品久久人妻99蜜桃| 国产精品免费大片| 精品国产国语对白av| 少妇粗大呻吟视频| 亚洲视频免费观看视频| 国产成人精品久久二区二区91| 91国产中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 黄色毛片三级朝国网站| 人人妻,人人澡人人爽秒播| 黄片小视频在线播放| 大型av网站在线播放| 久久久精品国产亚洲av高清涩受| 久久久久久久精品精品| 女人高潮潮喷娇喘18禁视频| 91成年电影在线观看| 一区在线观看完整版| 免费av中文字幕在线| 亚洲国产欧美在线一区| 我的亚洲天堂| 脱女人内裤的视频| 国产精品欧美亚洲77777| 亚洲成国产人片在线观看| 各种免费的搞黄视频| 日本vs欧美在线观看视频| 国产福利在线免费观看视频| 色94色欧美一区二区| 欧美黄色片欧美黄色片| 深夜精品福利| 国产伦理片在线播放av一区| 一区福利在线观看| 国产又爽黄色视频| 日韩欧美一区视频在线观看| 中文字幕人妻丝袜制服| 国产成人免费观看mmmm| 亚洲国产欧美网| 中文字幕制服av| 动漫黄色视频在线观看| avwww免费| 久久精品国产亚洲av香蕉五月 | 婷婷色av中文字幕| 日本精品一区二区三区蜜桃| 俄罗斯特黄特色一大片| 日本a在线网址| videosex国产| 国产日韩欧美在线精品| 一区二区三区激情视频| 亚洲精品第二区| 欧美日韩精品网址| 精品福利观看| 9191精品国产免费久久| 精品少妇内射三级| 人人妻人人添人人爽欧美一区卜| 天堂8中文在线网| 国产精品一区二区在线不卡| 国产男人的电影天堂91| av又黄又爽大尺度在线免费看| 国产精品免费大片| 超碰成人久久| 欧美变态另类bdsm刘玥| 在线观看免费日韩欧美大片| 久久久精品94久久精品| 天堂8中文在线网| 欧美性长视频在线观看| 一区福利在线观看| 亚洲精品日韩在线中文字幕| 亚洲av片天天在线观看| 伊人久久大香线蕉亚洲五| 黄色视频不卡| 欧美日韩一级在线毛片| 99国产精品一区二区三区| 免费黄频网站在线观看国产| 国产免费现黄频在线看| 大码成人一级视频| 女人爽到高潮嗷嗷叫在线视频| 好男人电影高清在线观看| 久久av网站| 日韩大码丰满熟妇| av国产精品久久久久影院| 免费在线观看影片大全网站| 热99国产精品久久久久久7| 纵有疾风起免费观看全集完整版| kizo精华| 亚洲精品在线美女| av线在线观看网站| 日韩精品免费视频一区二区三区| 黄片小视频在线播放| 美女高潮到喷水免费观看| 久久国产精品男人的天堂亚洲| 99re6热这里在线精品视频| 午夜精品国产一区二区电影| 侵犯人妻中文字幕一二三四区| 激情视频va一区二区三区| 午夜福利影视在线免费观看| 手机成人av网站| 精品久久久精品久久久| 精品少妇一区二区三区视频日本电影| 免费在线观看影片大全网站| 少妇粗大呻吟视频| 日韩欧美一区视频在线观看| 欧美黄色片欧美黄色片| 亚洲中文av在线| 国产成人精品无人区| 国产色视频综合| 一区二区av电影网| 日本一区二区免费在线视频| 日韩一卡2卡3卡4卡2021年| 久久久久国产一级毛片高清牌| 亚洲成国产人片在线观看| 久久午夜综合久久蜜桃| 精品少妇久久久久久888优播| 日日摸夜夜添夜夜添小说| 国产1区2区3区精品| 欧美日韩国产mv在线观看视频| 欧美激情久久久久久爽电影 | 18禁国产床啪视频网站| 最近最新免费中文字幕在线| 男女国产视频网站| 成人av一区二区三区在线看 | 国产成人一区二区三区免费视频网站| 亚洲精品一卡2卡三卡4卡5卡 | 日本黄色日本黄色录像| 99久久国产精品久久久| 日韩免费高清中文字幕av| 国产深夜福利视频在线观看| 午夜福利一区二区在线看| 美女主播在线视频| 亚洲午夜精品一区,二区,三区| 在线天堂中文资源库| 婷婷丁香在线五月| 亚洲三区欧美一区| 在线看a的网站| 一级毛片女人18水好多| 欧美精品高潮呻吟av久久| 国产成人精品无人区| 黑人猛操日本美女一级片| 在线 av 中文字幕| 亚洲欧美色中文字幕在线| 中文欧美无线码| 亚洲色图综合在线观看| 一个人免费看片子| 亚洲激情五月婷婷啪啪| 窝窝影院91人妻| 99热全是精品| 永久免费av网站大全| 亚洲欧洲精品一区二区精品久久久| 日韩一区二区三区影片| 亚洲第一青青草原| 汤姆久久久久久久影院中文字幕| 99九九在线精品视频| tube8黄色片| 电影成人av| 一本—道久久a久久精品蜜桃钙片| 色老头精品视频在线观看| 日韩欧美一区二区三区在线观看 | 国产日韩欧美视频二区| 久久久水蜜桃国产精品网| 国产精品秋霞免费鲁丝片| 我要看黄色一级片免费的| 男人添女人高潮全过程视频| 亚洲激情五月婷婷啪啪| 欧美人与性动交α欧美软件| 宅男免费午夜| 精品久久久久久电影网| 一级黄色大片毛片| 国产伦人伦偷精品视频| 男女国产视频网站| 欧美+亚洲+日韩+国产| 大陆偷拍与自拍| 成年女人毛片免费观看观看9 | 久久久久久亚洲精品国产蜜桃av| 欧美日韩精品网址| 高清视频免费观看一区二区| 精品国产一区二区久久| 久久久久久久精品精品| 新久久久久国产一级毛片| 狠狠狠狠99中文字幕| 欧美久久黑人一区二区| 999久久久国产精品视频| 国产三级黄色录像| 精品国产一区二区久久| 9热在线视频观看99| 午夜免费观看性视频| 国产男女超爽视频在线观看| 亚洲国产av新网站| 老司机亚洲免费影院| 亚洲av成人不卡在线观看播放网 | 亚洲成国产人片在线观看| 黄片播放在线免费| 欧美 日韩 精品 国产| 久久天堂一区二区三区四区| 成人三级做爰电影| av在线老鸭窝| 精品久久久久久电影网| 伊人久久大香线蕉亚洲五| 欧美日韩成人在线一区二区| 亚洲专区国产一区二区| 在线看a的网站| 9191精品国产免费久久| videosex国产| 丁香六月欧美| 在线天堂中文资源库| 欧美精品av麻豆av| 午夜福利乱码中文字幕| 国产亚洲欧美精品永久| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 国产成人欧美| 国产成人精品无人区| 日韩大片免费观看网站| 精品人妻1区二区| 男人操女人黄网站| 丰满饥渴人妻一区二区三| 高清黄色对白视频在线免费看| 91麻豆av在线| 18禁观看日本| 久久人妻熟女aⅴ| av国产精品久久久久影院| 在线观看一区二区三区激情| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区 | 亚洲欧美一区二区三区黑人| 黑人操中国人逼视频| 一区二区日韩欧美中文字幕| 少妇人妻久久综合中文| 久久久久国产精品人妻一区二区| 美女大奶头黄色视频| 亚洲精品久久午夜乱码| 成年动漫av网址| 久9热在线精品视频| 亚洲国产欧美网| 日韩欧美国产一区二区入口| 大香蕉久久成人网| 午夜福利影视在线免费观看| 亚洲视频免费观看视频| 亚洲国产av影院在线观看| 一本大道久久a久久精品| 深夜精品福利| 成年女人毛片免费观看观看9 | 午夜老司机福利片| 亚洲精品中文字幕在线视频| 亚洲色图 男人天堂 中文字幕| 一级黄色大片毛片| 久久精品成人免费网站| 久久精品国产综合久久久| 热99久久久久精品小说推荐| 极品人妻少妇av视频| 少妇粗大呻吟视频| 男女边摸边吃奶| 王馨瑶露胸无遮挡在线观看| 宅男免费午夜| 国产精品 欧美亚洲| 99久久人妻综合| 亚洲av国产av综合av卡| 咕卡用的链子| 99精国产麻豆久久婷婷| 青青草视频在线视频观看| 日本91视频免费播放| 嫁个100分男人电影在线观看| 国产成+人综合+亚洲专区| 亚洲成av片中文字幕在线观看| 少妇的丰满在线观看| 国产精品99久久99久久久不卡| 久久av网站| 爱豆传媒免费全集在线观看| 侵犯人妻中文字幕一二三四区| 国产精品国产av在线观看| 亚洲天堂av无毛| 亚洲三区欧美一区| 99国产精品一区二区蜜桃av | 热99久久久久精品小说推荐| 男女国产视频网站| 亚洲欧美激情在线| 人妻一区二区av| 欧美人与性动交α欧美精品济南到| 亚洲一卡2卡3卡4卡5卡精品中文| 50天的宝宝边吃奶边哭怎么回事| 国产精品国产三级国产专区5o| 色老头精品视频在线观看| 国产日韩欧美视频二区| 丝袜喷水一区| 人妻一区二区av| 黄片大片在线免费观看| 久久人人爽av亚洲精品天堂| 国产一级毛片在线| 久久女婷五月综合色啪小说| 99香蕉大伊视频| 午夜91福利影院| 国产老妇伦熟女老妇高清| 精品亚洲成国产av| 男女床上黄色一级片免费看| 操美女的视频在线观看| 国产成人欧美| 嫁个100分男人电影在线观看| 99九九在线精品视频| 考比视频在线观看| 男女国产视频网站| 国产精品久久久av美女十八| 午夜老司机福利片| 麻豆乱淫一区二区| 国产成人影院久久av| 超碰成人久久| 国产av精品麻豆| 搡老乐熟女国产| 日韩电影二区| 成人国产av品久久久| 狠狠狠狠99中文字幕| 91大片在线观看| 午夜影院在线不卡| 国产av精品麻豆| 欧美另类亚洲清纯唯美| av在线app专区| 亚洲激情五月婷婷啪啪| 国产精品99久久99久久久不卡| 欧美激情高清一区二区三区| 亚洲av片天天在线观看| 亚洲精品国产精品久久久不卡| 亚洲精品av麻豆狂野| 欧美午夜高清在线| 亚洲欧洲日产国产| 在线观看一区二区三区激情| 99久久国产精品久久久| 国产亚洲av片在线观看秒播厂| 亚洲一卡2卡3卡4卡5卡精品中文| 天堂俺去俺来也www色官网| 国产免费福利视频在线观看| 一本综合久久免费| 国产亚洲欧美在线一区二区| 久久99热这里只频精品6学生| 日韩欧美一区二区三区在线观看 | 黄色视频在线播放观看不卡| 日本a在线网址| 国产精品 欧美亚洲| 他把我摸到了高潮在线观看 | 视频区图区小说| 飞空精品影院首页| 两性夫妻黄色片| 国产精品99久久99久久久不卡| 国产一区二区激情短视频 | 国产成人免费无遮挡视频| 亚洲精品一二三| 一区在线观看完整版| 飞空精品影院首页| 国产一区二区三区av在线| 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频| 免费一级毛片在线播放高清视频 | av电影中文网址| 19禁男女啪啪无遮挡网站| 91大片在线观看| 亚洲av成人不卡在线观看播放网 | 18禁观看日本| 亚洲熟女精品中文字幕| 韩国高清视频一区二区三区| 中文字幕人妻熟女乱码| 国产精品久久久久久精品古装| 欧美日韩福利视频一区二区| 午夜视频精品福利| 日本wwww免费看| 久久久久久人人人人人| 操美女的视频在线观看| 涩涩av久久男人的天堂| 极品人妻少妇av视频| 亚洲 国产 在线| 十分钟在线观看高清视频www| 成人国产一区最新在线观看| 国产国语露脸激情在线看| 亚洲精品中文字幕在线视频| 亚洲精品国产一区二区精华液| av电影中文网址| 亚洲视频免费观看视频| √禁漫天堂资源中文www| 免费少妇av软件| 热99re8久久精品国产| 欧美日韩av久久| 啦啦啦 在线观看视频| 汤姆久久久久久久影院中文字幕| 午夜福利在线观看吧| 91精品国产国语对白视频| 欧美黑人精品巨大| 极品人妻少妇av视频| 日韩 亚洲 欧美在线| 美女高潮到喷水免费观看| 十八禁网站免费在线| 熟女少妇亚洲综合色aaa.| 乱人伦中国视频| 久久香蕉激情| 欧美日韩av久久| 麻豆av在线久日| 黄色视频在线播放观看不卡| 欧美精品一区二区大全| 热99国产精品久久久久久7| 精品国产一区二区三区久久久樱花| 国产av国产精品国产| 一二三四社区在线视频社区8| 国产成人啪精品午夜网站| 99国产极品粉嫩在线观看| 99香蕉大伊视频| 色婷婷久久久亚洲欧美| 日韩人妻精品一区2区三区| 午夜激情av网站| 91国产中文字幕| 自线自在国产av| 日韩大片免费观看网站| 欧美亚洲日本最大视频资源| 成人国语在线视频| 国产91精品成人一区二区三区 | 老司机影院成人| 超色免费av| 亚洲av成人不卡在线观看播放网 | 精品人妻熟女毛片av久久网站| 99九九在线精品视频| 美女大奶头黄色视频| 老司机靠b影院| 操出白浆在线播放| 久久狼人影院| 亚洲 欧美一区二区三区| 中文字幕人妻丝袜制服| 天堂俺去俺来也www色官网| 日韩制服骚丝袜av| 亚洲精品久久成人aⅴ小说| 免费在线观看日本一区| 肉色欧美久久久久久久蜜桃| 建设人人有责人人尽责人人享有的| 久久久精品国产亚洲av高清涩受| avwww免费| 脱女人内裤的视频| 在线看a的网站| 亚洲精品国产区一区二| 中文字幕色久视频| 99国产精品免费福利视频| 亚洲精品国产av成人精品| 丝袜人妻中文字幕| 精品人妻一区二区三区麻豆| 正在播放国产对白刺激| 久久国产精品大桥未久av| 欧美日韩av久久| 精品熟女少妇八av免费久了| 精品国产一区二区三区四区第35| 免费在线观看日本一区| 大型av网站在线播放| 亚洲国产欧美日韩在线播放| 狂野欧美激情性xxxx| 亚洲男人天堂网一区| 搡老熟女国产l中国老女人| 色精品久久人妻99蜜桃| 在线观看www视频免费| 国产97色在线日韩免费| 亚洲男人天堂网一区| 亚洲欧美清纯卡通| 中文欧美无线码| 天天躁狠狠躁夜夜躁狠狠躁| 两个人看的免费小视频| 欧美激情高清一区二区三区| 999久久久精品免费观看国产| 夜夜夜夜夜久久久久| 亚洲,欧美精品.| 亚洲中文日韩欧美视频| 视频在线观看一区二区三区| 精品人妻1区二区| 欧美大码av| 亚洲性夜色夜夜综合| 秋霞在线观看毛片| 成人国产av品久久久| 高清黄色对白视频在线免费看| 婷婷成人精品国产| 777米奇影视久久| 久久久水蜜桃国产精品网| 手机成人av网站| 黄色片一级片一级黄色片| 精品一品国产午夜福利视频| 啦啦啦在线免费观看视频4| 欧美精品一区二区大全| 日韩欧美国产一区二区入口| 国产高清国产精品国产三级| 精品国产乱码久久久久久男人| 在线观看www视频免费| 亚洲第一欧美日韩一区二区三区 | 国产免费福利视频在线观看| 女人被躁到高潮嗷嗷叫费观| 老司机福利观看| www.av在线官网国产| 水蜜桃什么品种好| 久久久精品国产亚洲av高清涩受| 蜜桃在线观看..| 咕卡用的链子| 国产一区二区三区av在线| 日韩视频在线欧美| 中文字幕av电影在线播放| 亚洲黑人精品在线| 日韩制服丝袜自拍偷拍| www.熟女人妻精品国产| 成人免费观看视频高清| 又紧又爽又黄一区二区| 成人国产av品久久久| 亚洲激情五月婷婷啪啪| av在线app专区| 我的亚洲天堂| 啪啪无遮挡十八禁网站| 好男人电影高清在线观看| 欧美精品av麻豆av| 久久人妻熟女aⅴ| 亚洲国产成人一精品久久久| 欧美中文综合在线视频| 大陆偷拍与自拍| 欧美亚洲日本最大视频资源| 亚洲国产中文字幕在线视频| 亚洲性夜色夜夜综合| 亚洲视频免费观看视频| 一本大道久久a久久精品| 久9热在线精品视频| 两个人免费观看高清视频| 午夜91福利影院| 制服诱惑二区| 日韩欧美一区视频在线观看| 一区二区三区乱码不卡18| 久久九九热精品免费| 日本欧美视频一区| 亚洲中文日韩欧美视频| 国产黄频视频在线观看| 男女无遮挡免费网站观看| 99久久99久久久精品蜜桃| 日本av免费视频播放| 久久亚洲精品不卡| 国产在视频线精品| av福利片在线| 欧美老熟妇乱子伦牲交| 欧美精品啪啪一区二区三区 | 免费在线观看影片大全网站| 欧美少妇被猛烈插入视频| 免费观看人在逋| 亚洲国产欧美日韩在线播放| 男女国产视频网站| 成年美女黄网站色视频大全免费| 91字幕亚洲| 久久热在线av| 高清黄色对白视频在线免费看| 乱人伦中国视频| 午夜免费鲁丝| 日韩人妻精品一区2区三区|