《電子技術(shù)應(yīng)用》
您所在的位置:首頁 > 嵌入式技術(shù) > 設(shè)計應(yīng)用 > 帶激光筆的機器人主動視覺手眼標定
帶激光筆的機器人主動視覺手眼標定
2015年微型機與應(yīng)用第17期
黃朝興1,陳 丹1,唐旭晟2
(1.福州大學(xué) 電氣工程與自動化學(xué)院 先進控制技術(shù)研究中心,福建 福州 350108; 2.福州大學(xué) 機械工程及自動化學(xué)院,福建 福州 350108)
摘要: 以固高GRB-400機器人和攝像機組成手眼系統(tǒng),在手眼關(guān)系旋轉(zhuǎn)矩陣的標定方面,分析了基于主動視覺的標定方法。為實現(xiàn)手眼關(guān)系平移向量的標定,提出以固定于機械臂末端的激光筆來獲取工件平臺上特征點的基坐標,并結(jié)合已標定的旋轉(zhuǎn)矩陣來標定平移向量。最后,從圖像求取多個特征點之間的距離并與實際值進行誤差比較,平面特征點間的長度測量誤差在±0.8 mm之間,表明手眼標定精度較高,可滿足機器人進行工件定位與自動抓取的要求。
Abstract:
Key words :

  摘  要: 以固高GRB-400機器人和攝像機組成手眼系統(tǒng),在手眼關(guān)系旋轉(zhuǎn)矩陣的標定方面,分析了基于主動視覺的標定方法。為實現(xiàn)手眼關(guān)系平移向量的標定,提出以固定于機械臂末端的激光筆來獲取工件平臺上特征點的基坐標,并結(jié)合已標定的旋轉(zhuǎn)矩陣來標定平移向量。最后,從圖像求取多個特征點之間的距離并與實際值進行誤差比較,平面特征點間的長度測量誤差在±0.8 mm之間,表明手眼標定精度較高,可滿足機器人進行工件定位與自動抓取的要求。

  關(guān)鍵詞: 機器人;手眼系統(tǒng);主動視覺;激光筆;標定

0 引言

  當今,機器人視覺系統(tǒng)以其精度高、成本效率高、連續(xù)性、靈活性等特點被廣泛應(yīng)用于電子原件裝配、焊接工程、零配件尺寸檢查等領(lǐng)域。手眼標定是機器視覺系統(tǒng)中的關(guān)鍵技術(shù),手眼標定的精度將直接影響系統(tǒng)的定位精度。因此,為提高系統(tǒng)的定位精度,需要大量的實驗來實現(xiàn)手眼標定,以減小最終誤差。在手眼標定方面許多學(xué)者做了努力,并取得了一定的成果,如Ma[1]提出了主動視覺的方法,通過攝像機的正交平移運動,利用極點建立約束方程組,求解攝像機內(nèi)參數(shù)及手眼關(guān)系;楊廣林[2]等給出一種手眼標定方法,僅需場景中的兩個特征點,通過控制攝像機做兩次平移運動和一次旋轉(zhuǎn)運動即可實現(xiàn)手眼標定。此外,針對一些機器人如SCARA機器人,其在做旋轉(zhuǎn)運動時僅能繞z軸旋轉(zhuǎn),對于一些手眼標定方法該類機器人無法實現(xiàn),如Shiu[3]和Tsai[4]的手眼標定方法均要求至少需要旋轉(zhuǎn)軸不平行的兩組運動,才能唯一確定手眼矩陣的各分量。本文以固高GRB-400機器人和攝像機組成的手眼系統(tǒng)為研究對象,采用基于主動視覺的手眼標定方法對手眼關(guān)系的旋轉(zhuǎn)矩陣進行標定,借助于激光筆來獲取工件平臺上特征點的世界坐標,然后根據(jù)特征點的世界坐標和已標定的旋轉(zhuǎn)矩陣來求取手眼關(guān)系的平移向量。

1 固高GRB-400機器人

001.jpg

  四自由度GRB-400機器人是典型的SCARA機器人,具有3個旋轉(zhuǎn)關(guān)節(jié)和1個移動關(guān)節(jié),采用D-H法建立坐標系,如圖1所示。從控制器中獲取關(guān)節(jié)變量?茲1、?茲2、?茲4、d3,由圖1可求得各連桿的變換矩陣,從而得到機械臂末端的位姿參數(shù):

  1.png

  其中,9`NN$K1]EL{]O0IGN{JS)~9.jpg分別表示基坐標系與機械臂末端坐標系之間的旋轉(zhuǎn)矩陣和平移向量。

  易知其中的旋轉(zhuǎn)矩陣{AFAA@5U@KA0%IM{2P)@6)U.jpg具有如下形式:

  2.png

2 手眼關(guān)系的標定

  2.1 手眼系統(tǒng)中的坐標變換

002.jpg

  手眼關(guān)系示意圖如圖2所示,Qw-xyz、Qh-xyz和Qc-xyz分別表示機器人基坐標系、機械臂末端坐標系和攝像機坐標系,9`NN$K1]EL{]O0IGN{JS)~9.jpg分別表示基坐標系與機械臂末端坐標系之間的旋轉(zhuǎn)矩陣和平移向量,Rhc和Thc分別表示機械臂末端坐標系與攝像機坐標系之間的旋轉(zhuǎn)矩陣和平移向量。設(shè)P為空間中的一點,Pw、Ph和Pc分別為點P在坐標系Qw-xyz、Qh-xyz和Qc-xyz中的坐標,則P在手眼系統(tǒng)中各坐標系之間的關(guān)系如下:

  34.png

  2.2 手眼關(guān)系中旋轉(zhuǎn)矩陣的標定

  在手眼系統(tǒng)中,當控制機械手做平移運動時,攝像機也隨著機械手做平移運動。參考文獻[1]給出了確定攝像機坐標系平移方向的方法。

  當控制機械手做一次平移運動時,設(shè)其平移向量為  kFM[4~0W2GI57K[VFC4ZG584.jpg1,相應(yīng)的攝像機平移向量為kKMTASD)]WSY$MQ6V]FEHJSP.jpg1(KMTASD)]WSY$MQ6V]FEHJSP.jpg1、FM[4~0W2GI57K[VFC4ZG584.jpg1均為單位向量),設(shè)P點在機械臂坐標系平移前后的坐標分別為Ph和Ph+kFM[4~0W2GI57K[VFC4ZG584.jpg1,在攝像機坐標系平移前后的坐標分別為Pc和Pc+kKMTASD)]WSY$MQ6V]FEHJSP.jpg1,根據(jù)式(3)可得:

 56.png

  類似地,控制機械手沿不同方向做兩次平移運動,設(shè)其平移向量分別為kFM[4~0W2GI57K[VFC4ZG584.jpg2、kFM[4~0W2GI57K[VFC4ZG584.jpg3,相應(yīng)的攝像機平移向量分別為kKMTASD)]WSY$MQ6V]FEHJSP.jpg2、kKMTASD)]WSY$MQ6V]FEHJSP.jpg3(KMTASD)]WSY$MQ6V]FEHJSP.jpgi、FM[4~0W2GI57K[VFC4ZG584.jpgi均為單位向量,i=2,3),同樣可推得:

  78.png

  將式(6)~(8)合成矩陣形式有:

  9.png

  由此可得手眼關(guān)系的旋轉(zhuǎn)矩陣:

  10.png

  在實際標定過程中,可以只通過控制機械手沿相互正交的兩個方向平移,然后采用施密特正交化方法對得到的兩個攝像機平移向量進行正交化,第三個平移向量KMTASD)]WSY$MQ6V]FEHJSP.jpg3、FM[4~0W2GI57K[VFC4ZG584.jpg3由前兩項叉乘得出,即KMTASD)]WSY$MQ6V]FEHJSP.jpg3=KMTASD)]WSY$MQ6V]FEHJSP.jpgKMTASD)]WSY$MQ6V]FEHJSP.jpg2,FM[4~0W2GI57K[VFC4ZG584.jpg3=FM[4~0W2GI57K[VFC4ZG584.jpgFM[4~0W2GI57K[VFC4ZG584.jpg2,由此可求得正交的旋轉(zhuǎn)矩陣%8P~KFO%`4}ZV8RJQO~JY77.jpg

  2.3 手眼關(guān)系中平移向量的標定

  手眼關(guān)系的平移向量是通過以固定于機械臂末端的激光筆來獲取工件平臺上特征點的基坐標,并結(jié)合已標定的旋轉(zhuǎn)矩陣來標定的。

  2.3.1 基于激光筆求工件平臺上特征點的基坐標

  基于激光筆求特征點基坐標的示意圖如圖3所示,工件平臺與基坐標的z軸垂直(設(shè)工件平臺在z軸的坐標位置已知),P為工件平臺上位置固定的特征點,設(shè)其基坐標表示為(wxp,wyp,wzp)T,激光筆安裝于機械臂末端,即Qh-xyz坐標系的A處,且使激光射線平行于Qh-xyz的z軸。可通過控制機械手上下移動觀察激光射線是否打在工件平臺上的同一點來判斷激光射線是否平行于Qh-xyz的z軸,若不是則調(diào)整激光筆直到激光射線與Qh-xyz的z軸相平行。易知當激光射線與Qh-xyz的z軸相互平行時,激光射線上的點在Qh-xyz上具有相同的x、y軸坐標值且不隨機械手的移動而改變,設(shè)分別為hxp、hyp。

003.jpg

  首先,通過調(diào)整機器人的四個關(guān)節(jié)狀態(tài)使激光射線打在P點上(如圖3(a)所示),設(shè)此時P點在Qh-xyz中的坐標為(wxp,wyp,wzp)T,從控制器中獲取當前機器人的關(guān)節(jié)變量分別為[~KWIAGAZN$O07%N)[AP[FY.png,根據(jù)D-H法可求得機械臂末端的位姿參數(shù)HOX{N6F}K{G[C[%V`WV_3ET.png,易知旋轉(zhuǎn)矩陣{R[5XOS1PW7{FY0KCD`UD75.jpg具有如式(2)的形式,設(shè)為:

 1112.png 

  接著,保持機械臂末端與工件平臺的距離不變,即機器人的移動關(guān)節(jié)參數(shù)值仍為d3,改變機器人的其他三個旋轉(zhuǎn)關(guān)節(jié),且仍使激光射線打在P點上(如圖3(b)所示),顯然此時的P點在Qh-xyz中的坐標仍為(wxp,wyp,wzp)T,從控制器中獲取當前機器人的關(guān)節(jié)變量分別為[~KWIAGAZN$O07%N)[AP[FY.png,同理可求得XQ4047IT(10AV4I5R}9ST@L.jpg,且旋轉(zhuǎn)矩陣7N@K(%449G(R}QG}ZJ3~]0U.png具有式(2)的形式:

  13.jpg

  由式(12)、(13)可解得:

  14.png

  在基坐標系下,對于工件平臺上的特征點的z軸坐標值相等且為已知,因此只需求取x、y軸的坐標值。將求得的(hxp,hyp)T帶入式(12)即可求得P點在基坐標系下的x、y軸的坐標值。為減小(hxp,hyp)T的誤差,可采用對不同特征點重復(fù)多次實驗求平均的方法。

  2.3.2 平移向量的標定

  根據(jù)手眼系統(tǒng)中的坐標變換關(guān)系,由式(3)、(4)可得:

 1519.jpg

  標定過程:控制機械手使工件平臺上已知基坐標的特征點置于場景中,從控制器中獲取當前機器人的關(guān)節(jié)變量求得9`NN$K1]EL{]O0IGN{JS)~9.jpg,從圖像上獲取特征點的圖像像素坐標(u,v),zc可由參考文獻[2]給出的方法求得,6])K%3AGN@_4W]`7H)ZL)QN.jpg、M均已標定求出,則由式(19)即可求得手眼關(guān)系的平移向量4HYH~4WWP@IFZ0[AN29NL$V.png。通過重復(fù)多次實驗求平均,可減少誤差。

3 實驗結(jié)果及分析

  本文的手眼系統(tǒng)的硬件部分主要有:四自由度GRB-400機器人、機器人控制柜、控制計算機、CCD攝像頭、圖像采集卡及其他輔助設(shè)備。攝像機安置于機械臂手爪的旋轉(zhuǎn)關(guān)節(jié)上,其采集的圖像尺寸為2 048×       1 536。

  3.1 手眼關(guān)系旋轉(zhuǎn)矩陣的標定實驗

  控制機械手沿機械臂末端坐標系的Zh+(“+”表示坐標軸的正方向,“-”表示坐標軸的負方向)方向移動    10 mm,根據(jù)移動前后圖像特征對應(yīng)點求取攝像機的平移向量kKMTASD)]WSY$MQ6V]FEHJSP.jpg1,再控制機械手沿機械臂末端坐標系的Yh+方向移動10 mm,同樣可求得攝像機的平移向量kKMTASD)]WSY$MQ6V]FEHJSP.jpg2,利用Schmidt正交化,令KMTASD)]WSY$MQ6V]FEHJSP.jpg3=KMTASD)]WSY$MQ6V]FEHJSP.jpgKMTASD)]WSY$MQ6V]FEHJSP.jpg2,同時利用機械手平移的方向向量由式(10)求得手眼關(guān)系的旋轉(zhuǎn)矩陣%8P~KFO%`4}ZV8RJQO~JY77.jpg。同理控制機械手沿機械臂末端坐標系的Zh-、Yh+、Zh+、Yh-、Zh-、Yh-方向移動10 mm,每兩次移動就可確定一個%8P~KFO%`4}ZV8RJQO~JY77.jpg。旋轉(zhuǎn)矩陣的多次實驗結(jié)果數(shù)據(jù)如表1所示,其中最大偏差為旋轉(zhuǎn)矩陣的計算值與平均值之間的偏差矩陣中的最大值。從表1可看出旋轉(zhuǎn)矩陣的偏差是穩(wěn)定的,取平均值作為%8P~KFO%`4}ZV8RJQO~JY77.jpg的標定結(jié)果。即:

  19+.png

006.jpg

  3.2 手眼關(guān)系平移向量的標定實驗

  在求平移向量之前,需通過激光筆來獲取工件平臺上特征點的基坐標,因此先對激光筆在機械臂末端坐標系x、y軸的坐標位置hxp、hyp進行標定。通過多次標定得到的數(shù)據(jù)結(jié)果如表2所示,其中偏差為(hxp,hyp)T的實驗計算值與平均值之間的差值。

007.jpg

  由表2可以看出,實驗得出的hxp、hyp比較穩(wěn)定,取平均值作為hxp、hyp的標定結(jié)果。即:

  hxphyp=-6.950 2-46.765 7

  在求得hxp、hyp后,控制機器人使激光射線打在工件平臺上基坐標待求的特征點P上,從控制器中獲取當前機器人的關(guān)節(jié)變量求得機械臂末端的位姿參數(shù),由式(12)即可求得P點的基坐標。

  接著對手眼關(guān)系平移向量進行標定:控制機械手使工件平臺上基坐標已知的4個特征點置于場景中,從控制器中獲取當前機器人的關(guān)節(jié)變量求得9`NN$K1]EL{]O0IGN{JS)~9.jpg,從圖像上獲取這些特征點的圖像像素坐標,這些特征點的zc值按參考文獻[2]的方法另外求取,%8P~KFO%`4}ZV8RJQO~JY77.jpg、M均已知,由式(19)即可求得手眼關(guān)系的平移向量4HYH~4WWP@IFZ0[AN29NL$V.png。平移向量的標定結(jié)果數(shù)據(jù)如表3所示,其中最大偏差為平移向量的計算值與平均值之間的偏差向量中的最大值。

008.jpg

  由表3可看出,平移向量的最大偏差小于0.5 mm,取平均值為其標定結(jié)果。即:

  19++.jpg

  3.3 標定結(jié)果測試

004.jpg

  選取工件平臺上的4個點作為測試點,且將特征點置于不同的場景中,如圖4所示,通過求取特征點的基坐標以求得這4個特征點兩兩之間的距離,實驗結(jié)果數(shù)據(jù)如表4所示,其中實際值為由游標卡尺直接測量所得特征點之間的間距。

009.jpg

001.jpg

  由表4可得不同場景中特征點間距計算值與實際值的誤差曲線,如圖5所示。從圖5可以看出,特征點間距計算值與實際值的誤差在±0.8 mm之間,說明所標定的手眼關(guān)系的旋轉(zhuǎn)矩陣%8P~KFO%`4}ZV8RJQO~JY77.jpg和平移向量4HYH~4WWP@IFZ0[AN29NL$V.png精度較高,可滿足機器人的工件定位與自動抓取的要求。

005.jpg

010.jpg

  表5列出了本文方法和其他文獻方法的手眼標定精度。從表5可以看出,參考文獻[5]的標定精度最高,達到了0.5 mm,但其標定過程要求機器人的外臂坐標系到機器人基坐標系的旋轉(zhuǎn)矩陣以及末端執(zhí)行器的高度保持不變,具有一定的局限性。而參考文獻[1]和參考文獻[6]的手眼標定方法在求解平移向量時均需要機器人做旋轉(zhuǎn)運動,并且不能僅僅只繞某一軸旋轉(zhuǎn),若僅繞某一軸旋轉(zhuǎn)將無法求出平移向量。在標定精度方面,參考文獻[6]的標定精度與本文方法相當,而參考文獻[1]和參考文獻[5]的標定精度均低于本文方法。

4 結(jié)論

  本文提出了通過固定在機械臂末端的激光筆來獲取工件平臺上特征點的基坐標,以此來標定手眼關(guān)系的平移向量;通過場景中的兩個特征點,通過精確控制機械臂做2次特定的平移運動即可標定出手眼關(guān)系的旋轉(zhuǎn)矩陣。實驗結(jié)果表明,本文方法的標定精度較高(±0.8 mm),可滿足一般機器人工件自動抓取、跟蹤控制等作業(yè)的使用要求。此外,本文提出的方法對于平移向量的標定,只需要求機器人末端手臂垂直于工件平臺,并使固定在該臂上的激光筆的射線平行于該臂坐標系的z軸,這對于大多數(shù)的機器人都是易于實現(xiàn)的,因此本文的方法精度高,適用范圍廣。

  參考文獻

  [1] MA S D. A self-calibration technique for active vision systems [J]. IEEE Transactionsons on Robotics and Automation,1996,12(1):114-120.

  [2] 楊廣林,孔令富,王潔.一種新的機器人手眼關(guān)系標定方法[J].機器人,2006,28(4):400-405.

  [3] SHIU Y C, AHMAD S. Calibration of wrist-mounted robotic sensors by solving homogeneous transform equations of the form AX=XB[J]. IEEEE Transactions on Robotics & Automation, 1989,5(1):16-29.

  [4] TSAI R Y, LENZ R K. A new technique for fully autonomous and efficient 3D robotics hand/eye calibration[J]. IEEEE Transactions on Robotics & Automation, 1989,5(3):345-358.

  [5] 戚芳,席鋒.簡便高精度的機器人手眼視覺標定方法[J].傳感器與微系統(tǒng),2006,25(1):27-29.

  [6] 張黎爍,趙志梅.基于主動視覺的結(jié)構(gòu)光手眼系統(tǒng)自標定方法[J].科學(xué)技術(shù)與工程,2014,14(9):202-206.


此內(nèi)容為AET網(wǎng)站原創(chuàng),未經(jīng)授權(quán)禁止轉(zhuǎn)載。
亚洲一区二区欧美_亚洲丝袜一区_99re亚洲国产精品_日韩亚洲一区二区
99综合在线| 久久人91精品久久久久久不卡| 亚洲欧美日韩精品一区二区 | 一区二区三区四区五区精品| 欧美日本在线播放| 欧美r片在线| 欧美不卡一区| 亚洲一区二区少妇| 一区二区三区国产在线观看| 亚洲一区二区3| 在线视频欧美精品| 激情综合视频| 欧美午夜性色大片在线观看| 久久久人成影片一区二区三区| 欧美一区二区高清| 新片速递亚洲合集欧美合集| 日韩视频欧美视频| 日韩一区二区电影网| 午夜在线视频一区二区区别| 亚洲黄一区二区三区| 国产精品影音先锋| 欧美精品在线看| 久久免费精品视频| 久久久7777| 久久综合色婷婷| 免费国产一区二区| 欧美大香线蕉线伊人久久国产精品| 欧美插天视频在线播放| 欧美黄色免费网站| 欧美日韩中文字幕在线视频| 欧美视频在线观看一区| 裸体一区二区| 欧美激情亚洲一区| 欧美日韩一区二区三区在线视频| 乱码第一页成人| 嫩模写真一区二区三区三州| 久久精品成人欧美大片古装| 久久久久九九视频| 免费永久网站黄欧美| 欧美日产一区二区三区在线观看 | 国产欧美一区二区三区另类精品| 国产亚洲精品一区二区| 欧美视频在线免费看| 国产精品久久久久国产a级| 国产欧美日韩视频| 国产精品老牛| 国产一区视频观看| 国产女主播一区二区三区| 国产综合网站| 亚洲国产合集| 亚洲成人在线观看视频| 亚洲人成在线观看| 最新成人av在线| 亚洲午夜久久久久久尤物| 亚洲精品中文字幕在线| 亚洲丶国产丶欧美一区二区三区 | 久久免费一区| 欧美专区在线观看一区| 午夜精品福利在线| 久久久综合网| 欧美日韩精品系列| 国产中文一区| 一本高清dvd不卡在线观看| 欧美在线视频一区二区| 日韩视频在线免费观看| 亚洲日本中文字幕免费在线不卡| 国产精品99久久久久久有的能看 | 在线免费观看成人网| 一区二区免费看| 日韩网站在线观看| 新狼窝色av性久久久久久| 香蕉久久夜色精品国产| 亚洲精品乱码久久久久| 亚洲精品乱码久久久久久黑人| 亚洲国产精品视频一区| 亚洲在线视频观看| 老**午夜毛片一区二区三区| 国产精品国产三级国产普通话99| 在线观看亚洲一区| 亚洲欧美成人精品| 欧美一区二区精品在线| 久久riav二区三区| 亚洲午夜国产一区99re久久| 亚洲欧美日韩视频二区| 美女黄色成人网| 国产免费成人av| 一本久道久久久| 亚洲精品日韩一| 亚洲综合视频在线| 欧美国产日韩一区二区| 欧美日韩综合不卡| 国产精品―色哟哟| 国产亚洲欧美另类中文| 日韩香蕉视频| 亚洲人体1000| 狼人天天伊人久久| 国产一区二区三区自拍| 亚洲欧美99| 亚洲欧美日韩精品久久久久| 欧美在线影院| 国产精品久久久久久久久久久久久久 | 亚洲视频播放| 欧美大片免费观看| 好吊一区二区三区| 亚洲精品综合精品自拍| 亚洲欧美成人综合| 久久国产主播| 欧美在线|欧美| 国产精品蜜臀在线观看| 9i看片成人免费高清| 99视频精品免费观看| 久久av资源网| 国产精品综合视频| 亚洲性视频h| 亚洲欧美综合v| 免费视频一区| 国产精品激情电影| 日韩亚洲视频| 久久精品视频免费观看| 亚洲一级黄色片| 欧美体内谢she精2性欧美| 国产亚洲欧美在线| 午夜欧美理论片| 一区二区三区导航| 欧美日韩国产在线观看| 亚洲伦伦在线| 在线视频精品一区| 欧美视频在线免费看| 一区二区三区四区精品| 亚洲一区激情| 欧美国产亚洲视频| 国产日韩欧美| 一区二区三区国产| 亚洲欧美高清| 国产欧美日韩综合一区在线播放| 亚洲欧美日韩一区在线| 久久精品国产精品亚洲| 国产一区视频在线观看免费| 亚洲国产精品第一区二区三区| 美女精品自拍一二三四| 亚洲激情黄色| 亚洲一区二区成人| 国产精品天天摸av网| 亚洲免费观看高清在线观看 | 亚洲欧美国产精品专区久久| 欧美高清不卡在线| 国产私拍一区| 久久精品国产2020观看福利| 午夜电影亚洲| 欧美屁股在线| 亚洲高清精品中出| 欧美自拍偷拍| 欧美阿v一级看视频| 亚洲精品国产精品国自产在线| 亚洲午夜久久久久久久久电影院 | 在线亚洲+欧美+日本专区| 亚洲免费电影在线观看| 欧美日韩亚洲一区二区三区| 亚洲国产一区二区三区高清| 中国成人亚色综合网站| 欧美人与禽性xxxxx杂性| 在线观看日韩av| 99精品99| 欧美日韩国产一级片| 亚洲图片在线| 另类欧美日韩国产在线| 99精品热视频| 久久久久久久综合| 亚洲美女91| 一区二区三区免费看| 欧美精品日韩一本| 亚洲一区二区三| 美女91精品| 亚洲一区二区视频在线| 久久一区二区精品| 精品91久久久久| 亚洲国产第一页| 美女网站在线免费欧美精品| 日韩午夜av电影| 久久久欧美精品sm网站| 好男人免费精品视频| 亚洲丰满少妇videoshd| 蜜臀久久99精品久久久久久9 | 久久精品国产2020观看福利| 亚洲人成7777| 国产精品99久久久久久久女警| 国产欧美日本| 久久国产精品99久久久久久老狼| 欧美激情一区二区三区蜜桃视频| 亚洲午夜羞羞片| 午夜亚洲性色福利视频| 亚洲电影在线| 久久精品国产一区二区三| 日韩视频亚洲视频| 另类天堂av| 亚洲欧美国产77777| 欧美精品激情| 亚洲电影自拍| 欧美精品www在线观看| 亚洲精品社区|