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

    基于改進雞群優(yōu)化算法的機械臂軌跡規(guī)劃*

    2024-01-03 01:50:10李長勇
    組合機床與自動化加工技術 2023年12期
    關鍵詞:樣條公雞母雞

    馬 龍,李長勇,司 慶

    (新疆大學機械工程學院,烏魯木齊 830017)

    0 引言

    隨著工業(yè)的發(fā)展,機械臂在工業(yè)領域應用的越來越廣泛[1]。機械臂的軌跡規(guī)劃是機械臂控制領域的研究重點,機械臂軌跡規(guī)劃分為笛卡爾空間軌跡規(guī)劃和關節(jié)空間規(guī)劃,笛卡爾空間軌跡規(guī)劃是對機械臂末端執(zhí)行器的軌跡規(guī)劃,而關節(jié)空間軌跡規(guī)劃是對各個關節(jié)角度進行規(guī)劃。目前進行軌跡優(yōu)化的主要是對關節(jié)空間進行軌跡規(guī)劃,目的包括以下4種:縮短運動時間,降低運動能量,減少關節(jié)的運動沖擊,以及綜合考慮多目的[2]。

    當前諸多學者對機械臂的軌跡規(guī)劃進行了研究,劉貝貝等[3]為了減少沖擊和振動對機械臂造成的損害,采用一種5次多項式結合均勻5次B樣條的來進行軌跡規(guī)劃,得到了較為光滑的關節(jié)角度、角速度、角加速度關于時間的變化曲線,減少了運動的沖擊,但是對時間沒有優(yōu)化,無法提升機械臂工作效率;郝晗等[4]進行了遺傳算法的初始種群生成方式、變異算子和交叉算子的構造優(yōu)化研究,使用3次B樣條曲線進行軌跡插值,將改善后的遺傳算法運用到軌跡優(yōu)化中,仿真結果證明該方法有效的縮短了機械臂運行時間,但未考慮機械臂運動中沖擊造成振動與磨損的問題;袁旭華等[5]通過增加控制點對三次B樣條曲線進行完善,并給出了優(yōu)化搜索策略蜂群算法的解決途徑,提高了曲線擬合準確性,并減少了規(guī)劃時間,但是機械臂在軌跡起點的初始加加速度不是零且曲線不平滑,會造成機械臂機運動沖擊過大;孫颙琰等[6]在傳統(tǒng)蝙蝠算法中的局部搜索中加入動態(tài)干擾系數(shù),并改進局部搜索和全局搜索策略,使用5次多項式進行軌跡插值,將時間、能耗和沖擊設為優(yōu)化目標,使用改進算法對機械臂運動進行軌跡優(yōu)化,有效的提高了機械臂的運行效率,但是求解過程復雜,計算耗時長。

    目前,機械臂的軌跡優(yōu)化主要把時間最優(yōu)或沖擊最優(yōu)作為單一優(yōu)化目標,機械臂在工廠工作中,要長時間搬運重物,為了提高工作效率、減少關節(jié)磨損,綜合考慮機械臂運動時間和沖擊,進行六自由度機械臂軌跡優(yōu)化研究。利用5次B樣條在關節(jié)空間構造插值曲線,并加入雞群優(yōu)化算法。將鯨魚算法與蜉蝣算法融入到雞群優(yōu)化算法中,提高算法收斂速度;引入萊維飛行策略,進而提升算法全局搜索能力;引入變異操作,防止陷入局部最優(yōu)。通過改進的雞群優(yōu)化算法對時間和沖擊進行優(yōu)化,規(guī)劃出關節(jié)軌跡,達到縮短時間減小沖擊的目標。

    1 機械臂運動學分析

    為了研究機械臂軌跡規(guī)劃,對機械臂進行運動學分析必不可少,本文中以埃夫特ER3A-C60型機械臂為研究對象,如圖1所示,其擁有6個關節(jié)并均為旋轉副,能承受較大載荷,可以搬運較重的零件。對機械臂關節(jié)建立連桿固連坐標系,如圖2所示。

    用標準D-H參數(shù)法對該機械臂進行模型建立,該機械機器臂D-H參數(shù)如表1所示。

    表1 D-H參數(shù)

    使用MATLAB機器人工具箱建立模型,世界坐標系與機械臂的基坐標系相互重合,減少計算量,機械臂在零位時的模型如圖3所示,借助搭建的模型,可以驗證機械臂正、逆運動學正確性。

    圖3 ER3A-C60型機器人模型

    1.1 機械臂的運動學分析

    (1)

    式中:c為cos,s為sin。

    (2)

    (3)

    式中:[n,o,a]3×3表示機械臂末端相對于基坐標的姿態(tài)矩陣,[Px,Py,Pz]T表示機械臂末端相對于基坐標的位置矩陣。

    1.2 軌跡規(guī)劃模型

    機械臂在工作過程中執(zhí)行器末端在笛卡爾空間會通過一系列軌跡點,{oi},i=0,1,…,m,其中o0代表機械臂末端起始軌跡點,om代表機械臂末端終止軌跡點,通過機械臂逆運動學求解,可以推導出軌跡點對應的關節(jié)空間位置為{θij},j=1,2,…,n,其中n是機械臂關節(jié)的個數(shù)[5]。優(yōu)化目標定義為:

    (4)

    (5)

    式中:F1代表機械臂運動總時間,ti代表從軌跡點oi到oi+1的時間,i=1,2,…,m-1,F2代表機械臂關節(jié)沖擊。

    歸一化目標函數(shù)如下:

    Ff=KtF1+KjF2+KpFp

    (6)

    式中:Kt為時間系數(shù),Kj為沖擊系數(shù),Kp為懲罰項系數(shù),設定Kt=0.4,Kj=0.6,Kp=5;Fp是懲罰項,初始值是0,當關節(jié)超過各約束都會使其值增加1。

    機械臂運動過程中,需要滿足以下約束函數(shù):

    (7)

    1.3 B樣條插值構造

    B樣條曲線具有較高的靈活性,可以進行局部微小的修改操作,而且各個關節(jié)需要滿足關節(jié)角度、角速度、角加速度和角加加速度的約束,使用五次B樣條曲線對機械臂進行軌跡規(guī)劃插值擬合[7-8],得到角度、角速度、角加速度、角加加速度都平滑連續(xù)的函數(shù)曲線。B樣條曲線方程為:

    (8)

    (9)

    式中:u是變量,di是控制頂點,規(guī)定0/0=0;Bi,k(t)是k次B樣條基函數(shù),區(qū)間u∈[ui,ui+k+1]定義節(jié)點矢量U=[u0,u1,u,…,u2k],對于一點ui最多有k+1個非0的k次B樣條Br,k(u),其中r∈[i-1,i],其他k次B樣條基函數(shù)都為0。因此B樣條曲線可另外表示為:

    (10)

    將插值曲線的定義域轉換成規(guī)范定義域,節(jié)點矢量為:

    (11)

    對時間節(jié)點歸一化,得到:

    (12)

    得到滿足插值條件的n+1個方程:

    (13)

    式中:ui+k∈[uk,un+k],i=0,1,…,n。

    使用切矢邊條件約束可以得到k-1個條件,5次B樣條曲線可以得到:

    (14)

    式中:vs和ve是關節(jié)啟停速度,as和ae是關節(jié)啟停加速度。根據(jù)德布爾遞推公式求出各階矢導:

    (15)

    (16)

    根據(jù)上述公式可得個n+5個方程,以矩陣形式表示:

    Ajdj=pj

    (17)

    (18)

    通過式(17)可反推求解出控制點:

    (19)

    根據(jù)式(12)、式(19)、式(9)求出的時間節(jié)點向量、控制頂點向量和B樣條基函數(shù);由式(10)得到關節(jié)在t∈[t0,tn]的軌跡曲線,由式(15)可以得出各階導曲線,可在求解最優(yōu)軌跡規(guī)劃時構造約束函數(shù)。

    2 雞群優(yōu)化算法

    常用的軌跡優(yōu)化算法有遺傳算法,粒子群算法等,但上述算法針對多目標優(yōu)化存在求解精度低,后期迭代速度慢等問題,而雞群優(yōu)化算法不需要對問題進行特殊處理,比如不需要引入特殊的編碼方式、特殊的交配和突變操作等,有著收斂速度快,收斂精度高的優(yōu)點,因此本文選用雞群優(yōu)化算法來進行機械臂的軌跡優(yōu)化。

    2.1 原始雞群優(yōu)化算法(CSO)

    雞群優(yōu)化算法中種群分為3類:公雞、母雞和小雞,將雞群按適應度值排序,靠前的NR個是公雞,最后的NC個是小雞,剩下的是母雞;將雞群分為NR個組,保證每個組中都有一個公雞,母雞隨機分配到各個組中,與公雞確立跟隨關系;隨機選取NH個母雞,小雞隨機分配到母雞身邊,與母雞確立母子關系。且不同的種群的位置更新方法不同。若第i只雞用i表示,d表示解的維度,表示解中的第j維,G表示當前迭代輪數(shù),則表示在第G輪迭代時,第i只雞在第j維上的位置[9]。公雞的位置更新方法如式(20)、式(21)所示。

    (20)

    (21)

    式中:N(0,σ2)表示均值為0、方差為σ2的正態(tài)分布;k代表不同于公雞i的另一只公雞k,f為每只公雞的最優(yōu)適應度值,ε為計算機中最小的常數(shù),其作用是防止公式中的分母為0造成無法計算。

    母雞跟隨自己組中的公雞尋找食物,母雞位置更新方法如式(22)~式(24)所示。

    (22)

    (23)

    S2=exp(fr2-fi)

    (24)

    式中:R表示區(qū)間[0,1]上的隨機數(shù),abs()表示取絕對值運算,r1為母雞i所在組內的公雞,r2為雞群優(yōu)化中除了母雞i和與其同組的公雞r1之外的另一只公雞或母雞。

    小雞跟隨自己的母雞進行覓食,小雞的位置更新方法如式(25)所示。

    (25)

    式中:m表示小雞i的母親,F表示小雞跟隨母雞進行覓食的調節(jié)參數(shù),由于小雞群優(yōu)化體中的每個個體也存在差異性,因此F為區(qū)間[0,2]上的隨機數(shù)。

    2.2 改進的雞群優(yōu)化算法(ICSO)

    公雞的位置更新方式呈現(xiàn)正態(tài)分布,其隨意性造成公雞錯過全局最優(yōu)解,本文將鯨魚算法[10]的包圍機制和蜉蝣算法[11]的速度機制融入公雞的位置更新方式中,增強全局搜索能力,ICSO中的公雞的位置更新方法為:

    (26)

    (27)

    S=2Ab2-A

    (28)

    (29)

    (30)

    基于萊維的運動方式是非破壞覓食環(huán)境中覓食者最佳搜索方法[12],為了使母雞和小雞可以更精準的模仿雞的覓食,引入萊維飛行策略,進而提升算法全局搜索能力。ICSO算法中的母雞的位置更新方法為:

    (31)

    (32)

    (33)

    式中:L(λ)表示萊維飛行產(chǎn)生的步長,β為常數(shù),取值1.5;μ、v為均值為0、方差為1的正態(tài)分布;Γ為Gamma函數(shù)[13]。

    小雞跟隨母雞,為了增強其全局搜索能力,將小雞跟隨母雞進行覓食的調節(jié)參數(shù)替換成萊維飛行步長,ICSO算法中的小雞位置更新方法為:

    (34)

    雞群迭代前期,種群差異較大,但是在迭代后期,個體之間差異程度降低,導致種群多樣性降低,致使過早收斂,陷入局部最優(yōu)。為了優(yōu)化雞群優(yōu)化算法這個缺陷,本文引入隨機差分變異操作[14],其定義為:

    X(G+1)=Q(pbest-X(G))+Q(X′(G)-X(G))

    (35)

    式中:Q為差分變異算子,pbest為當前最優(yōu)個體,X′(G)為群體中隨機選出的個體,X(G)為進行變異操作的個體,其作用是對差分值進行縮放。

    Q取較大值時,可以讓算法搜索范圍變廣,但是會導致收斂速度減慢;Q取值較小的時候,雖然搜索速度加快,但是會導致陷入局部最優(yōu)。為了避免上述問題,對差分變異算子進行改進,使算法前期搜索大范圍,隨著迭代進行加快搜索速度,進后的定義如下:

    (36)

    式中:M為最大迭代次數(shù),G為當前迭代次數(shù),rand為[0,1]之間的隨機數(shù)。

    雞群迭代前期,不需要對太多個體進行變異操作,因此,設定前期隨機選取變異個體少,隨著迭代進行,種群差異變多,增加變異個體,對變異數(shù)量的設定為:

    (37)

    式中:min()為取其中最小值,int表示取整操作,Hm為變異數(shù)量上限,定為種群規(guī)模的20%。

    2.3 ICSO算法步驟

    算法基本步驟為:

    步驟1:對算法參數(shù)初始化,設定雞群大小N、最大迭代次數(shù)M、個體位置維度、公雞母雞小雞的比例;

    步驟2:對雞群進行初始化,將雞群按適應度值排序,確定公雞、母雞與小雞,并進行分組,確定母雞與小雞之間的對應關系;

    步驟3:開始迭代,判斷是否更新分組,是否跟新集群中的關系,需要則更新;否則公雞、母雞與小雞各自按照式(26)、式(31)、式(34)對自己的位置進行更新,同時計算更新后的適應度值;

    步驟4:比較新舊位置適應度值,保留小適應度值的位置;

    步驟5:根據(jù)式(35)對種群中隨機個體進行差分變異操作;

    步驟6:判斷是否滿足迭代停止條件,若是滿足。那么算法結束,輸出最優(yōu)解;否則回到步驟3。

    2.4 性能測試

    為了驗證ICSO的有效性,使用MATLAB R2018A平臺對ICSO進行性能測試仿真試驗,同時添加CSO、WOA兩種算法進行對比。設定群數(shù)量M為30,最大迭代次數(shù)為500次,所有測試函數(shù)如式(38)~(43)所示,實驗結果如圖4所示。

    f1=max{|xi|,1≤i≤n},[-100,100]n

    (38)

    (39)

    (40)

    (41)

    (42)

    (43)

    式中:f1為單峰測試函數(shù),用來衡量算法求解精度;f2、f3為多峰測試函數(shù),用來度量算法全局搜索能力;f4、f5、f6為固定維測試函數(shù),用來測試算法的求解復雜多維能力。

    由圖4a可知,在求解f1時隨著迭代次數(shù)的增加,ICSO比WOA和CSO達到更低的函數(shù)值,說明ICSO具有更好的收斂精度。由圖4b和圖4c可知,在求解f2、f3時,起始階段ICSO就開始快速收斂,并能跳出局部最優(yōu),求解f2時ICSO盡管收斂速度比WOA慢,但是能在140次左右達到最優(yōu)值,且收斂值比WOA低,而CSO收斂速度慢而且收斂值差;求解f3時盡管3種算法收斂值一樣,但ICSO在前期能多次跳出局部最優(yōu)值,最快到達全局最優(yōu),體現(xiàn)了ICSO良好的全局尋優(yōu)能力。由圖4d~圖4f可知,在求解f4時,ICSO在10次就完成收斂;在求解f5和f6時,ICSO收斂值最低,且不斷跳出局部最優(yōu)解,驗證了ICSO的有效性。

    圖4 測試函數(shù)仿真結果

    3 實驗與分析

    為了驗證本文提出的ICSO在機械臂軌跡規(guī)劃中的有效性,以埃夫特ER3A-C60型機械臂為研究對象,使用MATLAB中機器人工具箱對其建模,如圖3所示,可通過其中Teach模塊控制機械臂運動。機械臂在笛卡爾空間內經(jīng)過一條軌跡,其中包括起始點和終止點的8個軌跡點對應的各個關節(jié)的角度信息如表2所示,將這些關節(jié)角度信息記錄到電腦中。

    表2 軌跡點關節(jié)角度

    在MATLAB中采用五次B樣條,CSO優(yōu)化和ICSO優(yōu)化分別對起始點到目標點的共8個軌跡點的的運動軌跡規(guī)劃仿真對比實驗。智能算法迭代次數(shù)設置為200次,雞群數(shù)量設置為30,公雞、母雞、小雞的比例為1:2:2,以式(6)中F1、F2為優(yōu)化指標,得出不同方案下的各關節(jié)角度、角速度、角加速度和沖擊隨時間的變化曲線如圖5所示,CSO和ICSO算法軌跡規(guī)劃收斂曲線如圖6所示。

    (a) 關節(jié)1

    圖6 軌跡規(guī)劃收斂曲線

    由圖5可知,ICSO收斂值小于CSO,且在迭代20次、40次時均可跳出局部最優(yōu)。由圖6可知,不同方法下規(guī)劃的曲線都是連續(xù)光滑,但能看出經(jīng)過CSO優(yōu)化和ICSO優(yōu)化,機械臂運動時間減少,沖擊曲線起伏變小。具體的,各段軌跡運行時間、總運行時間如表3所示,不同方案下關節(jié)最大沖擊如表4所示。

    表3 不同方案機械臂運動時間 (s)

    表4 不同方案下關節(jié)最大沖擊 ((°)·s-3)

    由表3可知,CSO優(yōu)化后比B樣條軌跡規(guī)劃的曲線時間減少9%,而ICSO優(yōu)化后時間減少12%。由表4可知,在經(jīng)過CSO優(yōu)化后,比B樣條軌跡規(guī)劃,關節(jié)1沖擊減少42%,關節(jié)2沖擊減少66%,關節(jié)3沖擊減少12%,關節(jié)4沖擊減少65%,關節(jié)5沖擊減少67%;而在ICSO優(yōu)化后,比B樣條軌跡規(guī)劃,關節(jié)1沖擊減少63%,關節(jié)2沖擊減少78%,關節(jié)3沖擊減少25%,關節(jié)4沖擊減少73%,關節(jié)5沖擊減少74%。使用CSO和ICSO軌跡優(yōu)化后關節(jié)最大沖擊減小,但ICSO效果更優(yōu)。因此本文提出的ICSO算法不僅降低了機械臂運行時間,還降低了機械臂運動時的沖擊,延長了機械臂使用壽命。

    4 結論

    基于改進雞群優(yōu)化算法的機械臂軌跡規(guī)劃,該算法參考鯨魚算法、蜉蝣算法來改良公雞位置更新方式,引入萊維飛行策略優(yōu)點來改良母雞和小雞的位置更新方式,提升全局搜索能力,并引入變異操作,提高了種群的復雜多樣性,得到收斂速率快,全局搜索能力強大的雞群優(yōu)化算法。

    將ICSO應用于機械臂的軌跡優(yōu)化中,并與5次B樣條軌跡規(guī)劃、CSO軌跡優(yōu)化結果對比,實驗結果表明,基于ICSO的軌跡優(yōu)化,在機械臂運動中,減少了路徑運動時間,且各個關節(jié)最大沖擊,大大減少,為機械臂軌跡優(yōu)化領域的研究提供一定的借鑒。

    猜你喜歡
    樣條公雞母雞
    一元五次B樣條擬插值研究
    母雞
    兩只公雞
    快樂語文(2021年36期)2022-01-18 05:48:40
    母雞下蛋
    母雞
    三次參數(shù)樣條在機床高速高精加工中的應用
    三次樣條和二次刪除相輔助的WASD神經(jīng)網(wǎng)絡與日本人口預測
    軟件(2017年6期)2017-09-23 20:56:27
    基于樣條函數(shù)的高精度電子秤設計
    說話的公雞
    小母雞回來了
    幼兒園(2016年6期)2016-04-26 10:57:10
    x7x7x7水蜜桃| 久99久视频精品免费| 美女cb高潮喷水在线观看| 午夜福利视频1000在线观看| 九色国产91popny在线| 少妇高潮的动态图| 欧美一区二区精品小视频在线| 香蕉av资源在线| 免费观看人在逋| 在线免费观看不下载黄p国产 | 成人亚洲精品av一区二区| 在线国产一区二区在线| 亚洲专区中文字幕在线| 九九热线精品视视频播放| 成人午夜高清在线视频| 最近在线观看免费完整版| 级片在线观看| 亚洲色图av天堂| 此物有八面人人有两片| 一个人免费在线观看的高清视频| 久久久国产成人免费| 色av中文字幕| 99国产精品一区二区三区| 波多野结衣高清无吗| 天堂√8在线中文| 在线观看午夜福利视频| 日本撒尿小便嘘嘘汇集6| 亚洲人成网站在线播| netflix在线观看网站| 国产精品久久久久久久电影| 亚洲av免费在线观看| 日韩av在线大香蕉| 午夜福利成人在线免费观看| 伊人久久精品亚洲午夜| 少妇丰满av| 一级毛片久久久久久久久女| 两人在一起打扑克的视频| 看黄色毛片网站| 九九在线视频观看精品| 99热只有精品国产| 国产精华一区二区三区| 亚洲欧美日韩高清在线视频| 免费在线观看影片大全网站| 欧美黑人欧美精品刺激| 91久久精品国产一区二区成人| 一区二区三区免费毛片| 国产精品久久电影中文字幕| 国产探花在线观看一区二区| 日本黄大片高清| 国产精品一及| 午夜精品在线福利| 美女被艹到高潮喷水动态| 成年女人看的毛片在线观看| 国产成年人精品一区二区| 亚洲精华国产精华精| 99热这里只有精品一区| 国产老妇女一区| 久久这里只有精品中国| 一边摸一边抽搐一进一小说| 99热这里只有是精品50| 91字幕亚洲| 成年女人永久免费观看视频| 可以在线观看的亚洲视频| 十八禁国产超污无遮挡网站| 欧美三级亚洲精品| 亚洲av五月六月丁香网| 丁香欧美五月| 夜夜爽天天搞| 欧美高清性xxxxhd video| 亚洲欧美日韩东京热| 久久国产乱子伦精品免费另类| 欧美一区二区国产精品久久精品| 成人无遮挡网站| 国产单亲对白刺激| 国产黄片美女视频| 国产精品一及| 欧美日韩福利视频一区二区| 可以在线观看毛片的网站| 国产精品一区二区性色av| 久久精品国产清高在天天线| 麻豆成人午夜福利视频| 小蜜桃在线观看免费完整版高清| 亚洲成人精品中文字幕电影| 夜夜爽天天搞| 人妻丰满熟妇av一区二区三区| 国产午夜福利久久久久久| 欧美黑人欧美精品刺激| 无人区码免费观看不卡| 91在线观看av| 午夜福利18| 亚洲性夜色夜夜综合| 狂野欧美白嫩少妇大欣赏| 美女被艹到高潮喷水动态| 日韩精品青青久久久久久| 国产三级黄色录像| 少妇人妻精品综合一区二区 | 国产伦精品一区二区三区视频9| 国产精品三级大全| 亚洲av熟女| 日韩有码中文字幕| 亚洲av不卡在线观看| 看黄色毛片网站| 久久久久亚洲av毛片大全| 九九热线精品视视频播放| 久久99热6这里只有精品| 内射极品少妇av片p| www日本黄色视频网| 国产高清有码在线观看视频| www.熟女人妻精品国产| 午夜福利在线观看免费完整高清在 | 亚洲 欧美 日韩 在线 免费| 亚洲av电影在线进入| 午夜激情福利司机影院| 国产熟女xx| 日本五十路高清| av在线天堂中文字幕| 久久久国产成人精品二区| 国产精品99久久久久久久久| 丰满的人妻完整版| 精品人妻熟女av久视频| 欧美激情久久久久久爽电影| 波多野结衣高清作品| 国产人妻一区二区三区在| 中文字幕高清在线视频| 免费看a级黄色片| 国产国拍精品亚洲av在线观看| 成人av在线播放网站| 色综合婷婷激情| 亚洲国产精品999在线| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品久久男人天堂| 大型黄色视频在线免费观看| 丝袜美腿在线中文| 怎么达到女性高潮| 亚洲国产精品sss在线观看| 久久久久久久精品吃奶| 日韩大尺度精品在线看网址| 亚洲欧美日韩无卡精品| 午夜免费激情av| 久久国产精品人妻蜜桃| av在线观看视频网站免费| 日本免费a在线| 成年女人毛片免费观看观看9| 人妻制服诱惑在线中文字幕| 国产精品一区二区三区四区久久| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久久电影| 制服丝袜大香蕉在线| 内地一区二区视频在线| 在线观看av片永久免费下载| 丝袜美腿在线中文| 成人欧美大片| 亚洲精品日韩av片在线观看| 精品欧美国产一区二区三| 小蜜桃在线观看免费完整版高清| 国产高清激情床上av| 免费观看的影片在线观看| 国产av麻豆久久久久久久| 最好的美女福利视频网| 中国美女看黄片| 国产成人a区在线观看| 国产精品影院久久| 日本 欧美在线| 日韩大尺度精品在线看网址| 日韩中字成人| 熟女人妻精品中文字幕| 白带黄色成豆腐渣| 丰满乱子伦码专区| 久久伊人香网站| 成年版毛片免费区| 色哟哟·www| 久久亚洲真实| 国产伦精品一区二区三区四那| 最好的美女福利视频网| 在现免费观看毛片| 久久久久国产精品人妻aⅴ院| 在线免费观看的www视频| 国产成人aa在线观看| 91狼人影院| 久久这里只有精品中国| 美女 人体艺术 gogo| 天天一区二区日本电影三级| 国产三级黄色录像| 午夜久久久久精精品| 久久热精品热| 啪啪无遮挡十八禁网站| 国产欧美日韩精品一区二区| 欧美在线黄色| 别揉我奶头~嗯~啊~动态视频| 国产伦在线观看视频一区| 九九久久精品国产亚洲av麻豆| 麻豆国产av国片精品| 熟妇人妻久久中文字幕3abv| 哪里可以看免费的av片| 色尼玛亚洲综合影院| 三级国产精品欧美在线观看| а√天堂www在线а√下载| 中文字幕人妻熟人妻熟丝袜美| 大型黄色视频在线免费观看| 老司机午夜十八禁免费视频| 欧美成人a在线观看| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 一个人观看的视频www高清免费观看| 午夜精品一区二区三区免费看| 国产精品久久久久久亚洲av鲁大| www.www免费av| av专区在线播放| 日韩欧美精品v在线| 国产亚洲av嫩草精品影院| 成人一区二区视频在线观看| 少妇被粗大猛烈的视频| 国产精品美女特级片免费视频播放器| 国产精品1区2区在线观看.| 国产国拍精品亚洲av在线观看| 极品教师在线视频| 嫩草影视91久久| 97热精品久久久久久| 99久久精品一区二区三区| 少妇被粗大猛烈的视频| 97超级碰碰碰精品色视频在线观看| 三级男女做爰猛烈吃奶摸视频| 久久久久久久午夜电影| 久久国产精品人妻蜜桃| 一级a爱片免费观看的视频| 久久久色成人| 国产白丝娇喘喷水9色精品| 波多野结衣高清作品| 老司机福利观看| 欧美激情国产日韩精品一区| www日本黄色视频网| 女生性感内裤真人,穿戴方法视频| 亚洲国产精品sss在线观看| 亚洲精品影视一区二区三区av| 国产综合懂色| 成人高潮视频无遮挡免费网站| 久久久久久久久大av| 午夜老司机福利剧场| 老司机福利观看| 午夜福利在线观看免费完整高清在 | 久久亚洲精品不卡| 午夜精品久久久久久毛片777| 欧美日韩瑟瑟在线播放| 亚洲,欧美精品.| 国产三级黄色录像| 老司机福利观看| 日韩中文字幕欧美一区二区| 中文字幕熟女人妻在线| 亚洲av电影不卡..在线观看| 国产午夜福利久久久久久| 美女xxoo啪啪120秒动态图 | av在线观看视频网站免费| 日韩av在线大香蕉| 久久国产精品人妻蜜桃| 午夜福利在线观看免费完整高清在 | 欧美bdsm另类| 我要搜黄色片| 精品久久久久久成人av| h日本视频在线播放| 日本 av在线| 日本a在线网址| 精品99又大又爽又粗少妇毛片 | 午夜视频国产福利| 亚洲成人精品中文字幕电影| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久久久久久久| 国产精品久久久久久久电影| 长腿黑丝高跟| 能在线免费观看的黄片| 99久久99久久久精品蜜桃| 最好的美女福利视频网| 最近中文字幕高清免费大全6 | 日韩欧美国产在线观看| 少妇被粗大猛烈的视频| 久久草成人影院| 国产三级中文精品| 久久九九热精品免费| 黄色配什么色好看| 日本a在线网址| 国产精品人妻久久久久久| 国产精品一区二区三区四区久久| av在线老鸭窝| 一进一出好大好爽视频| 国产精品永久免费网站| 欧美bdsm另类| 色av中文字幕| 久久精品国产自在天天线| 在线免费观看不下载黄p国产 | 国产三级黄色录像| 窝窝影院91人妻| 精品乱码久久久久久99久播| 国产精品亚洲av一区麻豆| 91久久精品电影网| 欧美成人a在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 我要搜黄色片| 欧美乱色亚洲激情| 色综合站精品国产| av专区在线播放| 国产成年人精品一区二区| 亚洲第一电影网av| 搡女人真爽免费视频火全软件 | 69av精品久久久久久| 国产午夜精品论理片| 国产成人福利小说| 丁香六月欧美| 国产成人福利小说| av专区在线播放| 波多野结衣高清作品| 精品人妻一区二区三区麻豆 | 成年版毛片免费区| 免费人成视频x8x8入口观看| 国产精品不卡视频一区二区 | 蜜桃久久精品国产亚洲av| 村上凉子中文字幕在线| 俺也久久电影网| 日本熟妇午夜| 日韩精品青青久久久久久| 丝袜美腿在线中文| 国产白丝娇喘喷水9色精品| 天美传媒精品一区二区| 中文字幕高清在线视频| 露出奶头的视频| 免费在线观看日本一区| 亚洲av免费在线观看| 精品人妻一区二区三区麻豆 | 小说图片视频综合网站| 夜夜爽天天搞| 麻豆成人av在线观看| 最近中文字幕高清免费大全6 | 久久香蕉精品热| 狠狠狠狠99中文字幕| 国产一区二区激情短视频| 婷婷色综合大香蕉| 我要看日韩黄色一级片| 国产成人av教育| 伊人久久精品亚洲午夜| 神马国产精品三级电影在线观看| 欧美激情国产日韩精品一区| 日韩中字成人| a在线观看视频网站| 女人被狂操c到高潮| 变态另类成人亚洲欧美熟女| 床上黄色一级片| 欧美极品一区二区三区四区| 日本黄色视频三级网站网址| 成人特级av手机在线观看| 久久久久精品国产欧美久久久| 精品人妻偷拍中文字幕| 夜夜爽天天搞| 赤兔流量卡办理| 亚洲片人在线观看| 久久久色成人| 欧美精品啪啪一区二区三区| 国产高清视频在线观看网站| 成人亚洲精品av一区二区| 亚洲专区中文字幕在线| 嫩草影院入口| 午夜免费男女啪啪视频观看 | 人妻久久中文字幕网| 每晚都被弄得嗷嗷叫到高潮| 日韩av在线大香蕉| 黄色日韩在线| 欧美黑人巨大hd| 又黄又爽又免费观看的视频| 一个人观看的视频www高清免费观看| 久久久久国产精品人妻aⅴ院| 51午夜福利影视在线观看| 国产伦一二天堂av在线观看| 日韩欧美免费精品| 成人高潮视频无遮挡免费网站| 最近中文字幕高清免费大全6 | 久久香蕉精品热| 久久久久九九精品影院| 一区二区三区激情视频| 麻豆国产av国片精品| 12—13女人毛片做爰片一| 午夜精品久久久久久毛片777| 一进一出抽搐gif免费好疼| 在现免费观看毛片| 又粗又爽又猛毛片免费看| 中文字幕精品亚洲无线码一区| 午夜福利高清视频| 精品无人区乱码1区二区| 97碰自拍视频| 亚洲成a人片在线一区二区| 国产三级在线视频| 亚洲av不卡在线观看| 久久精品夜夜夜夜夜久久蜜豆| 特大巨黑吊av在线直播| 在线观看一区二区三区| 亚洲av不卡在线观看| 精品国产三级普通话版| 最近最新中文字幕大全电影3| 亚洲激情在线av| 国产综合懂色| 国产成人福利小说| 在线观看舔阴道视频| 在线观看美女被高潮喷水网站 | 国产一区二区三区在线臀色熟女| 搞女人的毛片| 亚洲国产精品成人综合色| 免费看光身美女| 一区二区三区四区激情视频 | 1024手机看黄色片| 久久欧美精品欧美久久欧美| 最新在线观看一区二区三区| 亚洲欧美日韩无卡精品| 国产久久久一区二区三区| 国产伦精品一区二区三区视频9| 精品人妻一区二区三区麻豆 | 国产一级毛片七仙女欲春2| 一级作爱视频免费观看| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | АⅤ资源中文在线天堂| 91在线观看av| 日韩av在线大香蕉| 成年女人看的毛片在线观看| 日韩有码中文字幕| 精品人妻视频免费看| 精品午夜福利视频在线观看一区| 欧美午夜高清在线| 美女被艹到高潮喷水动态| 舔av片在线| 特级一级黄色大片| 99久久无色码亚洲精品果冻| 亚洲成人免费电影在线观看| 能在线免费观看的黄片| 精品久久久久久,| 亚洲av不卡在线观看| 丰满的人妻完整版| 国产欧美日韩精品一区二区| 欧美性猛交╳xxx乱大交人| 欧美精品啪啪一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产精品98久久久久久宅男小说| 亚洲精品久久国产高清桃花| 麻豆久久精品国产亚洲av| 欧美日韩福利视频一区二区| 人人妻人人澡欧美一区二区| 国产精品久久久久久久久免 | 成年免费大片在线观看| 国产高清激情床上av| 天堂影院成人在线观看| 国产精品日韩av在线免费观看| 免费一级毛片在线播放高清视频| 亚洲av电影在线进入| 欧美激情国产日韩精品一区| 97热精品久久久久久| 我的老师免费观看完整版| 国内精品久久久久久久电影| 午夜影院日韩av| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| bbb黄色大片| 国产熟女xx| 国产色婷婷99| av天堂在线播放| 国产伦人伦偷精品视频| 欧美成人一区二区免费高清观看| 波多野结衣巨乳人妻| 亚洲熟妇中文字幕五十中出| 国产 一区 欧美 日韩| 九九在线视频观看精品| 欧美激情在线99| 国产又黄又爽又无遮挡在线| www日本黄色视频网| 麻豆av噜噜一区二区三区| 91久久精品电影网| 高清毛片免费观看视频网站| 伊人久久精品亚洲午夜| 欧美一区二区亚洲| 99精品在免费线老司机午夜| 香蕉av资源在线| x7x7x7水蜜桃| 无遮挡黄片免费观看| 欧美黄色淫秽网站| 99久久久亚洲精品蜜臀av| 午夜福利视频1000在线观看| 三级国产精品欧美在线观看| 免费看美女性在线毛片视频| 国产色婷婷99| 琪琪午夜伦伦电影理论片6080| 久久久久久久亚洲中文字幕 | 久久久久久国产a免费观看| 黄色丝袜av网址大全| 午夜福利在线在线| 久久亚洲真实| av福利片在线观看| 一区二区三区免费毛片| 欧美最新免费一区二区三区 | 亚洲人成伊人成综合网2020| 国产高清视频在线播放一区| 一区二区三区免费毛片| 亚洲中文字幕日韩| 久久久久国内视频| 亚洲经典国产精华液单 | 久久久久国产精品人妻aⅴ院| 日本a在线网址| 成人国产综合亚洲| 18美女黄网站色大片免费观看| 日日摸夜夜添夜夜添av毛片 | 精品久久国产蜜桃| 高潮久久久久久久久久久不卡| 亚洲av五月六月丁香网| 免费大片18禁| 亚洲午夜理论影院| 欧美一区二区亚洲| 欧美绝顶高潮抽搐喷水| 91在线精品国自产拍蜜月| 黄色配什么色好看| www.www免费av| 亚洲成人久久性| 热99re8久久精品国产| 午夜免费成人在线视频| 很黄的视频免费| 日韩精品青青久久久久久| 最近视频中文字幕2019在线8| 丁香六月欧美| 欧美性猛交黑人性爽| 国产成人a区在线观看| 免费电影在线观看免费观看| 老司机午夜十八禁免费视频| 精品人妻1区二区| aaaaa片日本免费| 波多野结衣高清无吗| 又爽又黄无遮挡网站| 2021天堂中文幕一二区在线观| 国产真实乱freesex| 欧美日韩黄片免| 久久九九热精品免费| 日本黄色片子视频| 精品欧美国产一区二区三| 精品国产亚洲在线| 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线观看二区| 99久久精品国产亚洲精品| 欧美不卡视频在线免费观看| 女生性感内裤真人,穿戴方法视频| 中文字幕久久专区| 久久99热这里只有精品18| 乱码一卡2卡4卡精品| 男女下面进入的视频免费午夜| 极品教师在线视频| 热99在线观看视频| 一本精品99久久精品77| 国产伦精品一区二区三区视频9| 一级av片app| 久久国产乱子免费精品| 欧美黄色淫秽网站| 亚洲真实伦在线观看| 亚洲七黄色美女视频| 欧美乱色亚洲激情| a级毛片a级免费在线| 黄色一级大片看看| 欧美黑人巨大hd| 精品国产亚洲在线| 久99久视频精品免费| 每晚都被弄得嗷嗷叫到高潮| 99热6这里只有精品| 久久久久久九九精品二区国产| 亚洲七黄色美女视频| 中文字幕高清在线视频| 国产精品三级大全| 免费av观看视频| 亚洲av第一区精品v没综合| 看片在线看免费视频| a级毛片免费高清观看在线播放| 中文字幕高清在线视频| 性色av乱码一区二区三区2| 日韩精品中文字幕看吧| 嫁个100分男人电影在线观看| 国产亚洲精品久久久久久毛片| 亚洲欧美日韩卡通动漫| 国产伦一二天堂av在线观看| 老鸭窝网址在线观看| 中文资源天堂在线| 日韩欧美免费精品| 久久久久国产精品人妻aⅴ院| 午夜亚洲福利在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 两个人视频免费观看高清| 女同久久另类99精品国产91| 精品人妻熟女av久视频| 国产三级中文精品| 51国产日韩欧美| 精品人妻视频免费看| 在线播放国产精品三级| 日韩欧美三级三区| 9191精品国产免费久久| 女人十人毛片免费观看3o分钟| 亚洲av二区三区四区| 国模一区二区三区四区视频| 两个人的视频大全免费| 午夜激情福利司机影院| 欧美+亚洲+日韩+国产| 欧美丝袜亚洲另类 | av福利片在线观看| 成人午夜高清在线视频| 熟女电影av网| 美女 人体艺术 gogo| 亚洲精品成人久久久久久| 日本一本二区三区精品| 欧美一区二区国产精品久久精品| 免费看光身美女| 国产精品久久视频播放| 身体一侧抽搐| 亚洲欧美日韩高清专用| 欧美又色又爽又黄视频| 色av中文字幕| 国产主播在线观看一区二区| 久久国产精品人妻蜜桃|