• <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
    国产深夜福利视频在线观看| 激情在线观看视频在线高清| 亚洲一区高清亚洲精品| 电影成人av| 另类亚洲欧美激情| 在线观看免费视频网站a站| 欧美老熟妇乱子伦牲交| 亚洲精华国产精华精| 在线观看日韩欧美| 亚洲在线自拍视频| 一区二区三区激情视频| 在线观看免费视频网站a站| 亚洲中文av在线| 夜夜爽天天搞| 国产成人一区二区三区免费视频网站| 亚洲情色 制服丝袜| 两性午夜刺激爽爽歪歪视频在线观看 | 免费一级毛片在线播放高清视频 | 久久精品aⅴ一区二区三区四区| 淫秽高清视频在线观看| 国产av又大| 欧美中文综合在线视频| 色婷婷久久久亚洲欧美| 亚洲成人精品中文字幕电影 | 成熟少妇高潮喷水视频| 国产精品久久久久成人av| 悠悠久久av| 少妇被粗大的猛进出69影院| а√天堂www在线а√下载| 欧美日韩视频精品一区| 人妻丰满熟妇av一区二区三区| 女性被躁到高潮视频| 亚洲精品在线美女| av视频免费观看在线观看| 国产午夜精品久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 久久午夜亚洲精品久久| 国产成人av教育| 国产精品秋霞免费鲁丝片| 久久草成人影院| 久9热在线精品视频| 日本三级黄在线观看| 一本综合久久免费| 日韩中文字幕欧美一区二区| 男女床上黄色一级片免费看| 国产男靠女视频免费网站| 成人国产一区最新在线观看| 日韩欧美一区二区三区在线观看| 亚洲精华国产精华精| 中出人妻视频一区二区| 91老司机精品| 黄色 视频免费看| 国产精品久久久av美女十八| 俄罗斯特黄特色一大片| 久久精品亚洲av国产电影网| xxx96com| 色老头精品视频在线观看| 亚洲国产精品sss在线观看 | 一进一出抽搐动态| 伦理电影免费视频| 日本黄色视频三级网站网址| 亚洲一区二区三区色噜噜 | 99久久久亚洲精品蜜臀av| 久久99热6这里只有精品| 一区二区三区免费毛片| 99精品久久久久人妻精品| 成年人黄色毛片网站| 他把我摸到了高潮在线观看| 精品人妻偷拍中文字幕| 国产色爽女视频免费观看| 国产精品亚洲av一区麻豆| 两人在一起打扑克的视频| 别揉我奶头~嗯~啊~动态视频| 精品无人区乱码1区二区| 国产精品美女特级片免费视频播放器| 久久久久久久精品吃奶| 国产成人福利小说| 嫩草影院精品99| 舔av片在线| 国产一区二区激情短视频| 欧美+亚洲+日韩+国产| 人妻夜夜爽99麻豆av| 久久人人精品亚洲av| 日韩欧美精品v在线| 国产精品,欧美在线| 亚洲成人精品中文字幕电影| 亚洲精品成人久久久久久| 国产精品综合久久久久久久免费| 免费在线观看影片大全网站| 国产精品久久久久久久电影| 成人特级黄色片久久久久久久| 嫩草影视91久久| 久久精品91蜜桃| 欧美日韩综合久久久久久 | 久久婷婷人人爽人人干人人爱| 亚洲国产精品sss在线观看| 神马国产精品三级电影在线观看| 国产视频一区二区在线看| 国产在线男女| 男女那种视频在线观看| 免费在线观看亚洲国产| 一进一出抽搐动态| 精品午夜福利视频在线观看一区| 九九久久精品国产亚洲av麻豆| 欧美日韩黄片免| 又紧又爽又黄一区二区| 一a级毛片在线观看| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久久免 | 中文字幕av成人在线电影| 国产欧美日韩精品亚洲av| 久久亚洲真实| 欧美中文日本在线观看视频| 日本熟妇午夜| 亚洲av成人不卡在线观看播放网| 极品教师在线免费播放| 日本撒尿小便嘘嘘汇集6| 91在线观看av| 精品人妻偷拍中文字幕| 中文字幕人成人乱码亚洲影| 国产成人影院久久av| 蜜桃亚洲精品一区二区三区| 激情在线观看视频在线高清| 精品人妻视频免费看| 亚洲欧美日韩卡通动漫| 国产69精品久久久久777片| 国内久久婷婷六月综合欲色啪| 在线观看66精品国产| 日本与韩国留学比较| 69人妻影院| 午夜激情福利司机影院| 国产 一区 欧美 日韩| 国内精品美女久久久久久| 亚洲国产精品999在线| 日日干狠狠操夜夜爽| 51国产日韩欧美| 亚洲成av人片在线播放无| 国产真实伦视频高清在线观看 | 在线观看午夜福利视频| 精品日产1卡2卡| h日本视频在线播放| 免费看a级黄色片| 亚洲av美国av| 很黄的视频免费| 天堂网av新在线| 欧美一区二区精品小视频在线| 男女做爰动态图高潮gif福利片| 国产精品女同一区二区软件 | 免费在线观看亚洲国产| 午夜视频国产福利| 午夜福利欧美成人| 亚洲av熟女| 在线播放国产精品三级| 激情在线观看视频在线高清| 亚洲成人中文字幕在线播放| а√天堂www在线а√下载| 亚洲人成伊人成综合网2020| h日本视频在线播放| 搡老岳熟女国产| 国模一区二区三区四区视频| 亚洲精品456在线播放app | 欧美日韩综合久久久久久 | 国产午夜福利久久久久久| 国产黄a三级三级三级人| 日韩高清综合在线| 久久精品国产99精品国产亚洲性色| 99热精品在线国产| 亚洲成人精品中文字幕电影| 两个人视频免费观看高清| 日本在线视频免费播放| 网址你懂的国产日韩在线| 日本五十路高清| 成人午夜高清在线视频| 亚洲精品乱码久久久v下载方式| 国产视频一区二区在线看| 国产不卡一卡二| 偷拍熟女少妇极品色| 亚洲国产精品久久男人天堂| 一个人看的www免费观看视频| 丰满人妻一区二区三区视频av| 色精品久久人妻99蜜桃| 午夜久久久久精精品| 日韩中文字幕欧美一区二区| 日本 av在线| 1024手机看黄色片| 免费观看精品视频网站| 欧美bdsm另类| 欧美一区二区亚洲| av黄色大香蕉| 亚洲一区二区三区色噜噜| 99国产综合亚洲精品| 精品99又大又爽又粗少妇毛片 | 欧美丝袜亚洲另类 | 午夜福利高清视频| 免费搜索国产男女视频| 亚洲天堂国产精品一区在线| 最好的美女福利视频网| 中文字幕人成人乱码亚洲影| 国产精品女同一区二区软件 | 久久久久久九九精品二区国产| 亚洲av免费在线观看| 精华霜和精华液先用哪个| 欧美性猛交╳xxx乱大交人| 天堂动漫精品| 嫩草影院入口| 男插女下体视频免费在线播放| 男人狂女人下面高潮的视频| 在线观看舔阴道视频| 熟女电影av网| 亚洲精品在线美女| 亚洲乱码一区二区免费版| 91字幕亚洲| 99久久九九国产精品国产免费| .国产精品久久| 国产精品久久久久久精品电影| 嫁个100分男人电影在线观看| 国产伦精品一区二区三区四那| 人人妻,人人澡人人爽秒播| 国产三级黄色录像| 狂野欧美白嫩少妇大欣赏| 麻豆一二三区av精品| 日韩欧美国产一区二区入口| 欧美精品国产亚洲| 别揉我奶头~嗯~啊~动态视频| 国产色婷婷99| 亚洲专区国产一区二区| 高清日韩中文字幕在线| 99久久无色码亚洲精品果冻| 久久精品国产99精品国产亚洲性色| 午夜免费男女啪啪视频观看 | 国内毛片毛片毛片毛片毛片| 国产精品伦人一区二区| 国产在线精品亚洲第一网站| 色5月婷婷丁香| 国产爱豆传媒在线观看| 亚洲自偷自拍三级| 黄色视频,在线免费观看| 91午夜精品亚洲一区二区三区 | 欧美绝顶高潮抽搐喷水| 真人做人爱边吃奶动态| 欧美日韩黄片免| АⅤ资源中文在线天堂| 久99久视频精品免费| 窝窝影院91人妻| 永久网站在线| 亚洲成人中文字幕在线播放| 两个人视频免费观看高清| 国产精品伦人一区二区| 国产av麻豆久久久久久久| 国产黄片美女视频| 日韩欧美精品v在线| 国产精品,欧美在线| 美女xxoo啪啪120秒动态图 | 国产精品永久免费网站| 国产探花在线观看一区二区| 好看av亚洲va欧美ⅴa在| 日韩亚洲欧美综合| av福利片在线观看| 看片在线看免费视频| 亚洲欧美日韩东京热| av在线蜜桃| 精品久久久久久,| 亚洲成人中文字幕在线播放| 有码 亚洲区| 俺也久久电影网| 两人在一起打扑克的视频| 亚洲av中文字字幕乱码综合| 在线十欧美十亚洲十日本专区| 麻豆久久精品国产亚洲av| 好男人电影高清在线观看| 欧美日韩综合久久久久久 | 丝袜美腿在线中文| 岛国在线免费视频观看| 午夜福利成人在线免费观看| 内射极品少妇av片p| 成人av一区二区三区在线看| 美女免费视频网站| 欧美午夜高清在线| 偷拍熟女少妇极品色| 精品人妻熟女av久视频| 国产色爽女视频免费观看| 又紧又爽又黄一区二区| www.色视频.com| 免费搜索国产男女视频| 亚洲男人的天堂狠狠| 成人高潮视频无遮挡免费网站| 亚洲成人精品中文字幕电影| 国产精品一区二区三区四区久久| 国产精品女同一区二区软件 | 亚洲人与动物交配视频| 日本成人三级电影网站| 欧美激情在线99| 中文字幕高清在线视频| 哪里可以看免费的av片| 日本一本二区三区精品| 免费电影在线观看免费观看| 日韩欧美在线二视频| 亚洲精品在线观看二区| 国产淫片久久久久久久久 | 久久精品久久久久久噜噜老黄 | 97超视频在线观看视频| bbb黄色大片| 成人精品一区二区免费| av欧美777| 成人无遮挡网站| 99久国产av精品| h日本视频在线播放| 色在线成人网| 日日干狠狠操夜夜爽| 直男gayav资源| 在线天堂最新版资源| 天堂动漫精品| 蜜桃亚洲精品一区二区三区| 身体一侧抽搐| www.999成人在线观看| 精品人妻熟女av久视频| 少妇的逼水好多| 又粗又爽又猛毛片免费看| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 免费av不卡在线播放| 久久精品久久久久久噜噜老黄 | av在线天堂中文字幕| 黄色丝袜av网址大全| 99久久久亚洲精品蜜臀av| 美女cb高潮喷水在线观看| 国产精品不卡视频一区二区 | 亚洲aⅴ乱码一区二区在线播放| 人人妻人人看人人澡| 久久亚洲精品不卡| 99精品久久久久人妻精品| 看十八女毛片水多多多| 两个人视频免费观看高清| 国产精品电影一区二区三区| 国产免费av片在线观看野外av| 亚洲人成伊人成综合网2020| 变态另类成人亚洲欧美熟女| 99久久无色码亚洲精品果冻| 看黄色毛片网站| 热99在线观看视频| 校园春色视频在线观看| 欧美bdsm另类| 欧美中文日本在线观看视频| 身体一侧抽搐| 亚洲五月天丁香| 国产午夜精品论理片| 亚洲自偷自拍三级| 91麻豆av在线| 成人三级黄色视频| 国产美女午夜福利| 中文亚洲av片在线观看爽| 欧美成人a在线观看| 伦理电影大哥的女人| 午夜亚洲福利在线播放| 中文字幕高清在线视频| 黄色视频,在线免费观看| 久久久久久久亚洲中文字幕 | x7x7x7水蜜桃| 久久精品91蜜桃| 俺也久久电影网| 内射极品少妇av片p| 国产精品一区二区三区四区久久| 日本在线视频免费播放| 精品人妻视频免费看| 桃色一区二区三区在线观看| 国产乱人视频| 午夜激情福利司机影院| 精品一区二区三区视频在线观看免费| 男女之事视频高清在线观看| 成人国产综合亚洲| 色播亚洲综合网| 亚洲av五月六月丁香网| 久久亚洲真实| 国产亚洲精品久久久久久毛片| 精品久久久久久久末码| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣高清作品| 欧美黄色淫秽网站| 精品人妻视频免费看| 18禁在线播放成人免费| 欧美黄色淫秽网站| 757午夜福利合集在线观看| 成人特级黄色片久久久久久久| 国产一区二区在线观看日韩| 男插女下体视频免费在线播放| av专区在线播放| 国产一区二区在线观看日韩| 久久久久九九精品影院| 少妇的逼水好多| 日本 av在线| 国产精品久久视频播放| 国产精品一区二区免费欧美| 久久精品91蜜桃| 久久久久精品国产欧美久久久| 悠悠久久av| 免费大片18禁| 日本在线视频免费播放| 久久精品人妻少妇| 91在线精品国自产拍蜜月| 国产伦精品一区二区三区四那| 日韩欧美一区二区三区在线观看| 99热这里只有是精品在线观看 | 少妇裸体淫交视频免费看高清| 窝窝影院91人妻| 一级a爱片免费观看的视频| 亚洲精品456在线播放app | 亚洲欧美精品综合久久99| 麻豆久久精品国产亚洲av| 欧美成人一区二区免费高清观看| 在线国产一区二区在线| 最近视频中文字幕2019在线8| 国产综合懂色| 有码 亚洲区| 乱码一卡2卡4卡精品| 男女床上黄色一级片免费看| 啦啦啦韩国在线观看视频| 精品久久久久久久久av| 亚洲三级黄色毛片| 久久久久久久午夜电影| 欧美一级a爱片免费观看看| 亚洲片人在线观看| 亚洲欧美激情综合另类| 日韩免费av在线播放| 欧美性感艳星| 99热这里只有是精品50| 美女黄网站色视频| 成熟少妇高潮喷水视频| 搡女人真爽免费视频火全软件 | 99热只有精品国产| 国产av麻豆久久久久久久| 久久久久久久午夜电影| 精品国内亚洲2022精品成人| 久久99热这里只有精品18| 乱码一卡2卡4卡精品| 一二三四社区在线视频社区8| 真人做人爱边吃奶动态| 51午夜福利影视在线观看| 国产精品亚洲美女久久久| 午夜亚洲福利在线播放| 国产成+人综合+亚洲专区| 一夜夜www| 国产精品人妻久久久久久| 久久人人爽人人爽人人片va | 在线观看一区二区三区| 国产精品一区二区免费欧美| 综合色av麻豆| 亚洲最大成人中文| 成人午夜高清在线视频| 久久精品91蜜桃| 亚洲精品久久国产高清桃花| 每晚都被弄得嗷嗷叫到高潮| 一级黄片播放器| 两性午夜刺激爽爽歪歪视频在线观看| 天堂动漫精品| 亚洲人成网站高清观看| 国产精品日韩av在线免费观看| 观看免费一级毛片| 淫秽高清视频在线观看| 国产探花极品一区二区| 男人的好看免费观看在线视频| 身体一侧抽搐| 欧美黄色淫秽网站| 久久久久久久久大av| 国产精品电影一区二区三区| 亚洲18禁久久av| 美女被艹到高潮喷水动态| 免费大片18禁| 欧美日韩瑟瑟在线播放| 禁无遮挡网站| 亚洲精品一区av在线观看| 欧美色欧美亚洲另类二区| avwww免费| 中文字幕免费在线视频6| 人人妻人人澡欧美一区二区| 三级毛片av免费| 婷婷色综合大香蕉| 极品教师在线视频| 2021天堂中文幕一二区在线观| 午夜福利成人在线免费观看| 男人舔女人下体高潮全视频| 国产v大片淫在线免费观看| 我的女老师完整版在线观看| 国产色爽女视频免费观看| 3wmmmm亚洲av在线观看| 国产综合懂色| 亚洲av免费在线观看| 99久国产av精品| 欧美黄色片欧美黄色片| 午夜精品一区二区三区免费看| 香蕉av资源在线| 日本三级黄在线观看| 99久久精品热视频| 日韩欧美在线二视频| 哪里可以看免费的av片| 日日夜夜操网爽| 久久久成人免费电影| 精品久久久久久久久久久久久| 国产一区二区亚洲精品在线观看| 欧美高清成人免费视频www| 免费搜索国产男女视频| 免费无遮挡裸体视频| a级毛片a级免费在线| 国内精品一区二区在线观看| 亚洲黑人精品在线| 欧美乱色亚洲激情| 免费观看的影片在线观看| 成人特级黄色片久久久久久久| 高清毛片免费观看视频网站| 久久精品人妻少妇| 99在线视频只有这里精品首页| 深夜a级毛片| 中文字幕人妻熟人妻熟丝袜美| 久久精品综合一区二区三区| 国产成+人综合+亚洲专区| 嫩草影院入口| 欧美日韩综合久久久久久 | 欧美成人a在线观看| 亚洲色图av天堂| 午夜福利视频1000在线观看| 一区二区三区激情视频| 欧美色视频一区免费| 丁香欧美五月| 两人在一起打扑克的视频| 欧美日本亚洲视频在线播放| 午夜福利成人在线免费观看| 国产蜜桃级精品一区二区三区| 波多野结衣巨乳人妻| 我要看日韩黄色一级片| 熟女电影av网| 精品久久久久久久末码| av天堂中文字幕网| 亚洲欧美日韩高清专用| 97碰自拍视频| 亚洲精品乱码久久久v下载方式| 国产欧美日韩一区二区三| 丰满人妻熟妇乱又伦精品不卡| 一二三四社区在线视频社区8| 看片在线看免费视频| 欧美日韩中文字幕国产精品一区二区三区| 日韩国内少妇激情av| 午夜老司机福利剧场| 此物有八面人人有两片| 天堂影院成人在线观看| 国产精品女同一区二区软件 | 男女做爰动态图高潮gif福利片| 麻豆一二三区av精品| 亚洲五月天丁香| 久久欧美精品欧美久久欧美| 中文字幕人妻熟人妻熟丝袜美| 国产精品爽爽va在线观看网站| 亚洲av二区三区四区| 国产美女午夜福利| 国产精品久久久久久精品电影| 在线播放无遮挡| 国产精华一区二区三区| 搡老岳熟女国产| 欧美潮喷喷水| 中亚洲国语对白在线视频| 亚洲 国产 在线| 精品久久久久久久久久久久久| av在线天堂中文字幕| 久久精品91蜜桃| 亚洲精品久久国产高清桃花| 亚洲美女黄片视频| 男女那种视频在线观看| 欧美一区二区国产精品久久精品| 免费无遮挡裸体视频| 青草久久国产| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 国产亚洲欧美在线一区二区| 麻豆成人av在线观看| 国产亚洲av嫩草精品影院| 五月伊人婷婷丁香| 日韩欧美在线二视频| 久久久精品欧美日韩精品| 日本黄色视频三级网站网址| 99精品久久久久人妻精品| 国产伦精品一区二区三区视频9| 久久国产精品人妻蜜桃| 偷拍熟女少妇极品色| 婷婷精品国产亚洲av| 两人在一起打扑克的视频| 欧美bdsm另类| 精品人妻一区二区三区麻豆 | 身体一侧抽搐| 欧美另类亚洲清纯唯美| or卡值多少钱| 不卡一级毛片| 欧洲精品卡2卡3卡4卡5卡区| 噜噜噜噜噜久久久久久91| 亚洲精品在线美女| 精品久久久久久成人av| 亚洲性夜色夜夜综合| 日韩中文字幕欧美一区二区| 亚洲欧美精品综合久久99| 有码 亚洲区| 国产午夜精品论理片| 青草久久国产| 国产精品影院久久| 99热这里只有是精品50| 国产亚洲欧美在线一区二区| 亚洲成人免费电影在线观看| 精品人妻熟女av久视频| 国产真实伦视频高清在线观看 | 国产精品电影一区二区三区| 听说在线观看完整版免费高清| 亚洲无线观看免费| 国产真实乱freesex| 91av网一区二区| 看黄色毛片网站| 午夜精品在线福利| 国产蜜桃级精品一区二区三区| 亚洲三级黄色毛片| 亚洲五月天丁香|