谁有免费黄色网址_黄色电影免费提供_国产综合一区二区_欧美精品在线视频_免费观看日韩毛片_一级在线观看

  為機床工具企業提供深度市場分析                     

用戶名:   密碼:         免費注冊  |   申請VIP  |  

English  |   German  |   Japanese  |   添加收藏  |  
沈陽機床

車床 銑床 鉆床 數控系統 加工中心 鍛壓機床 刨插拉床 螺紋加工機床 齒輪加工機床
磨床 鏜床 刀具 功能部件 配件附件 檢驗測量 機床電器 特種加工 機器人

機器人

電工電力 工程機械 航空航天 汽車 模具
儀器儀表 通用機械 軌道交通 船舶

搜索
熱門關鍵字:

數控機床

 | 數控車床 | 數控系統 | 滾齒機 | 數控銑床 | 銑刀 | 主軸 | 立式加工中心 | 機器人
您現在的位置:機器人> 技術前沿>工業機器人位姿誤差空間 IDSW 插值補償方法研究
工業機器人位姿誤差空間 IDSW 插值補償方法研究
2017-7-11  來源:江蘇無錫  作者:陳宵燕 張秋菊 孫沂琳 陳海衛
        在焊接、激光切割、自動化裝配等實際生產和應用中,機器人精度是衡量其性能的重要指標[1]。 一般而言,工業機器人具有較高的重復定位精度,而絕對定位精度較差,前者可達到0.1 mm人,而后者一般2~3 mm,無法滿足某些行業的生產要求。
 
        目前提高絕對定位精度的方法主要分為兩種[2]:誤差預防法和誤差補償法。 前者主要是通過提高機器人設計和加工精度,這種方法不僅成本高而且不易實現。 后者也可稱為機器人標定,該方法通過使用先進的測量技術并且成功辨識出所建機器人模型中各項參數的準確值,然后通過修改機器人控制器中的參數或者增加一定的控制算法以達到提高定位精度的目的。 機器人標定是離線編程技術實用化的關鍵技術之一,其最后一步且關鍵一步為誤差補償[3],現階段誤差補償方法主要分為:基于模型的方法[4-7]和基于插補的方法[8-9]。

        基于模型的補償方法的基本原理是:在假設主要誤差源為關節軸傾斜和關節偏移等幾何誤差的情況下,用機器人運動學模型來描述機器人位姿誤差。校準中要進行機器人模型參數對測量姿態誤差的擬合,用校準后的機器人模型預測各個位姿的誤差。基于模型的補償方法主要分為兩種[10]:
 
        1) 關節空間補償。 利用標定后的結果直接在關節空間修正,對于某可達位姿,重新計算其關節值,來達到提高機器人位姿精度的目的;

        2) 微分誤差補償。該方法一般是針對基于微分變換思想建立機器人誤差模型,識別幾何參數名義值與真實值之間的微小偏差,并把其補償控制器的名義參數上,達到提高機器人精度的目的。

        基于模型的方法缺點是需建立復雜的運動學誤差模型,并需修改各項運動學參數從而造成機器人運動學逆解的困難,計算復雜且難以實現。

        對于基于插補思想的誤差補償法,需要對機器人末端在機器人工作空間中的特定網格點的位姿進行測量,目前多數采用在工作空間用雙線性插補和多項式擬合等插補技術來擬合相應于不同定位誤差的關節補償指令。 此外,文獻[9] 提出一種空間插值法,該方法在前期將機器人工作空間劃分成的網格,對網格點進行位置誤差測量,然后采用空間插值來補償定位誤差。 但這種方法存在很多缺陷,目前只在沒有考慮姿態變化對位置誤差的影響下對位置誤差進行了補償,且沒有對姿態誤差進行補償,因此存在很大的局限性。

        隨著高精度測量工具的不斷發展及普及應用,如FARO 激光跟蹤儀、LEICA 激光跟蹤儀等,使得基于插補的補償方法變得更為便捷、可靠,同時,完善并優化基于插補的補償方法也變得更為重要和緊迫。 針對空間插值補償方法,對機器人位姿誤差與其末端位置之間相關性進行探討,提出基于均勻數據場的空間IDSW 插值算法的機器人位姿誤差補償方法,并設計對比試驗,利用 MATLAB 進行仿真,結果證明該方法對機器人位姿精度提高的有效性。
 
        1.機器人運動學誤差模型

        DH 模型常用來對機器人空間的位姿進行表示和建模,這個建模方法由 Denavit 和 Hartenberg 兩人創新性提出的。 該建模方法優點在于其非常簡單且通用性非常好,可對任何構型的機器人的連桿和關節進行建模,方法為
 
=Rot(θ)Trans(0,0,)Trans(,0,0)Rot(α)        (1) 

        對于六自由度的串聯機器人,其末端的位置和姿態的坐標變換為
 
  

        經典的 DH 模型存在參數不連續的缺點:假如相鄰關節軸線實際上不平行而理論上卻以平行來建模,則這兩個軸線間的相對位置誤差的建模將不能滿足“小誤差建模” 條件,造成參數突變。 “小誤差模型”定義如下:對于連桿坐標變換的微小位姿誤差能通過連桿參數的微小偏差對其進行建模[11]。

        為了解決這種問題,Hayati 等[12]在經典 DH 運動學模型的基礎上,增加了一個新的繞y 軸的旋轉參數 βi,該模型被稱之為 MDH 模型。 該模型的相鄰連桿坐標變換的數學表達式為
 
Rot(θ)Trans(0,0,)Trans(,0,0)Rot(α)Rot(β)         (3)

        機器人定位精度誤差大體可以分為兩類:運動學(幾何)參數誤差和動力學(非幾何)參數誤差,具體包括幾何參數的偏差、關節偏移、連桿和關節的柔性、齒輪系結構、齒輪間隙、環境溫度等等。 運動學參數誤差主要由制造和裝配的誤差造成的,所導致的誤差占總誤差的 80%左右[13]。 所建立的機器人運動學誤差模型,主要考慮運動學參數誤差對機器人位姿誤差的影響。 除了以上運動學參數存在的誤差,在建立位姿誤差模型時,測量機器人末端實際位姿將引入機器人基坐標系相對于機器人構造的基坐標系之間的誤差 Δσ。 采用微分法[11]建立機器人運動學誤差模型,對式(3)進行全微分,得
 
θΔθΔΔαΔαβΔβ            ( 4 )
 
        式中: Δθi、Δdi、Δαi、Δai、Δβi分別表示連桿長度偏差、連桿偏置、扭角偏差、關節角偏差和相鄰平行軸的角度偏差。

        由于該方法較成熟,在此對過程不做贅述,具體可參見文獻[11],最終所建立的運動學誤差模型為
 
E = J·ε(5)
 
        式中:E 為位姿誤差;J 為辨識雅克比矩陣;ε 為運動學參數誤差。
 
  

        2  定位位姿誤差模型

        2.1  機器人位姿誤差與末端位置相關性
        空間插值是指通過已知的數據點或已知的已劃為各個相對小一些的區域內的數據點,計算出相關的其他未知點或相關區域內的任意點的方法[14]。 研究的空間IDSW 插值補償法是基于笛卡爾坐標系下的空間插值,即選擇一定數量的機器人末端位置作為測量點,測量這些點所對應的位姿誤差,根據這些已知點所對應的位姿誤差來計算相關區域內其它任意位置的位姿誤差。 由此可知,機器人位置和機器人位姿誤差兩者之間的高相關程度是本文插值補償法有效的先決條件。

        Pearson 相關(也稱為積差相關或積矩相關) 是英國統計學家皮爾遜于 20 世紀提出的一種計算直線相關的方法。 Pearson 相關系數用來衡量兩個數據集合是否在一條線上面。 其計算公式為
 
  
 
        式中相關系數 r 的值介于-1 與+1 之間。
 
        相關系數性質如下:
        1) 當r >0 時,表示兩變量正相關,r <0 時,兩變量為負相關;
        2) 當| r | =1 時,表示兩變量為完全線性相關,即為函數關系;
        3) 當r =0 時,表示兩變量間無線性相關關系。
 
        一般可按三級劃分:| r | <0.4 為低度線性相關;0.4 ≤| r | <0.7 為顯著性相關;0.7 ≤| r |<1 為高度線性相關。 低度線性相關又 可 分為 弱 相 關(0.2~0.4)和極弱相關或無相關(0.0~0.2)。

        利用 Pearson 相關系數來衡量機器人位置和機器人位姿誤差之間的相關性程度。以 Staubli TX60L機器人為例,根據第 1 節所述方法建立運動學誤差模型,可知機器人的位置誤差和姿態誤差都是機器人關節角 θi的函數,由于姿態誤差數量級較位置誤差小,所以將位姿誤差分為位置誤差和姿態誤差,并分別進行分析。
        利用 MATLAB 仿真功能,設計實驗 1,其步驟如下:

        1) Staubli TX60L 機器人名義運動參數及設定的運動參數偏差值如表 1 所示,利用蒙特卡洛法隨機選取 101 組 6 個關節角度值,記錄其所對應的101 組機器人末端位置值。

表 1  Staubli TX60L 機器人運動參數名義值和偏差值
 
  

         2) 將這 101 組關節值代入式(5),計算其對應的位置誤差 dΔ和位姿誤差 δΔ。

        3) 計算第 1 組數據所對應的位置誤差dΔ和位置值分別與其它 100 組對應的置誤差dΔ和末端位置值的位置誤差差值 ΔdΔ和位置值差值 dz,利用MATLAB 得到 ΔdΔ和 dz的兩者的相關性(用散點圖表示,見圖 1)和 Pearson 相關系數r1;
 
 圖 1  兩點的位置誤差差值與dz相關性

        4) 計算第 1 組數據所對應的姿態誤差δΔ和位置值分別與其它 100 組對應的位姿誤差δΔ和末端位置值的姿態誤差差值 ΔδΔ和位置值差值 dz,利用MATLAB 得到 ΔδΔ和 dz的兩者之間的相關性(用散點圖表示,見圖 2)和 Pearson 相關系數r2;
 
 
圖 2  兩點的姿態誤差差值與dz相關性
 

        利用 MATLAB 計算得到 ΔdΔ和 dz的 Pearson 相關系數 r1= 0.015 7, ΔδΔ和 dz的 Pearson 相關系數r2= -0.079 2。 圖 1 所描述的兩點所對應的位置誤差之間的差值與位置值間的差值不具有明顯的關系,根據計算所得相關系數r1可知兩變量的線性關系幾乎為無相關。 圖 2 所示兩點所對應的姿態誤差之間的差值與位置值間的差值也不具有明顯關系,同時根據所得的相關系數r2可知其線性關系也幾乎為無相關。 而對機器人來說,末端執行器的位姿由位置和姿態共同決定,所以在研究機器人位置和機器人位姿誤差兩者之間的相關性時,存在一個可控變量———姿態。 因此,設計實驗 2,其步驟如下:

        1) 利用蒙特卡洛法選取 Staubli TX60L 機器人工作空間 101 組機器人末端位置值,記錄其所對應的101 組機器人位置值,任意選取某一姿態,在該姿態下,對上述 101 個位置經運動學逆解得到 101 組6 個關節角度值;

        2) 同實驗 1 的步驟 2);

        3) 同實驗 1 的步驟 3),亦利用 MATLAB 得到ΔdΔ和 dz的兩者之間的相關性(用散點圖表示,見圖3 所示)和 Pearson 相關系數r3;
 
 
 
圖 3  兩點的位置誤差差值與dz相關性
 
       
        4) 同實驗 1 的步驟 4),亦利用 MATLAB 得到ΔδΔ和 dz的兩者之間的關系圖(用散點圖表示,如圖4 所示)和 Pearson 相關系數r4。
 
圖 4  兩點的姿態誤差差值與dz相關性
 

        利用 MATLAB 計算得到 ΔdΔ和 dz的 Pearson 相關系數 r3為 0.73,ΔδΔ和 dz的 Pearson 相關系數r4為 0.98。 比較圖 1 和圖 3,并根據所求的相關系數r1和 r3,可知通過控制機器人末端姿態這一變量后,兩點對應的位置誤差間的差值與位置值間的差值之間呈高度線性相關;分析比較圖 2 和圖 4,并根據所求的相關系數 r2和 r4,可知在機器人末端姿態不變的情況下,兩點對應的姿態誤差間的差值與位置值間的差值之間呈高度線性相關。 綜上所述,在機器人的工作空間內,若機器人末端姿態保持一致,則兩點所對應的位置值的差值 dz越大,所對應的位置誤差差值 ΔdΔ和姿態誤差差值 ΔδΔ也越大,且兩者關系均呈高度線性相關。 而此時,所提出的空間 IDSW插值補償法對機器人位姿誤差補償有效的前提條件成立。

        2.2  姿態變化因素下的定位位姿誤差模型

        根據 2.1 節的結論可知,為使空間 IDSW 插值補償方法能對機器人位姿誤差補償在任意姿態下均有效,必須考慮姿態變化對位姿誤差的影響。 所以將重點研究如何建立固定位置時姿態變化因素下的定位位姿誤差模型。

        機器人運動姿態的有兩種常用描述方法:機械人末端的運動方向和繞軸 x,y,z 的旋轉序列(歐拉角)。 前者是關于機器人運動學模型計算得到接近矢量 a,方向矢量n,法線矢量o,即式(2) 機器人末端坐標系變換矩陣的旋轉矩陣對應的 3 個列矢量;后者是笛卡爾坐標空間下,對機器人位置控制時運動姿態的 3 個重要參數,一般工業機器人的姿態表示法采用的是 x_y_z 歐拉角,即歐拉角Φ,θ,φ。 這種描述在基坐標系{B} 的方位的法則如下:基坐標系{B} 與參考系{A} 相同,首先使{B} 繞xB軸轉 Φ 角,然后繞yB軸轉 θ 角,最后繞zB軸轉 φ 角,如圖 5 所示。
 
 
圖 5 x?y?z 歐拉角法則
 

        根據右乘的原則,可以求得與之等價的旋轉矩陣,即(圖片)
由此可知,給定x-y-z 歐拉角 Φ,θ,φ, 便可確定對應的旋轉矩陣。當給定旋轉矩陣,也可求出等價的x-y-z 歐拉角Φ,θ,φ。 令
 
 
        由上述可知,繞zB軸轉 φ 角只影響方向矢量 n和法線矢量 o,與接近矢量a 無關,同時經研究可知旋轉角 φ 對位姿誤差的影響呈三角函數關系。 因此,對固定位置時姿態變化的定位位姿誤差模型確立分為兩步:正弦函數相位角偏移誤差建模和矢量夾角加權平均誤差建模。
 
 
圖 6  定位下所選取的 9 個姿態 

        首先對機器人姿態操作空間以接近矢量劃分,令 φ =0 只考慮歐拉角 Φ 和 θ,選取 9 個具有代表性的姿態,如圖 6 所示。 其次就是考慮歐拉角φ 對機器人位姿誤差的影響,對以上9 個姿態分別測量3 個不同歐拉角 φ 時的位姿誤差,建立正弦函數,即
 
                                   A = B*sin(φ - w)- C        (12)

        由此得到任意歐拉角φ時這9個姿態所對應的位姿誤差,這就是正弦函數相位角偏移誤差。

        對機器人的其他操作姿態所對應的位姿誤差可利用矢量夾角加權平均求得。 根據如圖 6 所示的已知 9 個姿態的每個姿態均有相對應的接近矢量、方向矢量和法線矢量,與所求姿態所對應的接近矢量、方向矢量和法線矢量分別進行夾角算,選取夾角量最小的 3個姿態并根據第一步獲得該 3 個姿態對應的位姿誤差來進行量夾角加權平均計算,得到該姿態下的位姿誤差。 綜合上述兩步就是建立的在姿態變化因素下的定位位姿誤差模型為
 
E′ =i =× λ   ( 1 3 )
 
        式中:E(i) 為 9個代表性姿態經第一步獲得所求姿態下對應的位姿誤差;λi為矢量夾角加權值。

        3  空間IDSW 插值補償模型

        對機器人位姿誤差補償是基于空間插值法,空間插值方法很多。 依據空間插值的基本假設和數學本質可分類為:幾何方、統計方法、空間統計方法、函數方法、隨機模擬方法、物理模型模擬方法和綜合方法[15]。 幾何方法包括反距離加權法、泰森多邊形法等;統計方法包括趨勢面、多元回歸法和多面函數法等;空間統計方法包括克里金法、最小二乘配置和楊赤中法等;函數方法包括樣條函數、雙線性內插方法[16]。 所提出的空間 IDSW 插值補償法選擇的是反距離加權法來進行機器人位姿誤差插值補償,該算法是以距離作為權重的滑動平均加權插值法[17],它是根據離散點生成規則網格數據最常見的空間內插方法,計算過程簡單且運算速度快。

        根據第 2.1 節的研究結果可知,當機器人末端姿態一致時,其末端任意兩點所對應的位置誤差間的差值 ΔdΔ和姿態誤差間的差值 ΔδΔ分別與該兩點位置間的差值 dz具有高度的線性相關性。 根據第2.2 節所建立的定位位姿誤差模型,可以將所測量點(即已知點) 的姿態所對應的位姿誤差轉換到需要預測的位姿的姿態所對應的位姿誤差,使得機器人末端姿態保持一致的條件成立,滿足了第 2.1 節所提出的空間插值補償法有效性的先決條件,從而為空間 IDSW 插值算法在機器人位姿誤差補償上奠定了理論基礎和應用基礎。

        根據采集點的分布可以采集數據分為規則采集點數據集點和離散采集點數據集。 前者可使反距離加權算法具有更好的逼近程度,本文將對機器人操作空間的規則采集點數據集稱為均勻數據場。 因此,將機器人的工作空間劃分為均勻網格,而對于某一目標點 M,使用快速搜索法得到距離M 點最近的8 個包絡點,如圖 7 所示。
 
 
圖 7  基于反距離加權算法的位姿補償原理 

        利用反距離加權法[18-19],可得位置誤差和姿態誤差的補償模型為:
 
Δi =Δ× ρ(14)
 
δΔi =δΔ× ρ(15 )

        式中:d*Δ和 δ*Δ為目標點的預測位置誤差和姿態誤差;dΔ(i)和δΔ(i)為目標點周圍 8 個包絡點的位置誤差和姿態誤差;ρi為 8 個包絡點的權重,權重的公式為
 
 ρj =(16)
 
        式中:di或 dj為目標點與周圍 8 個包絡點之間的距離;k 為冪指數。
 
         在式(16)中,k 值越大,插補結果越具有平滑,當取 k= 2 時, 通常稱為反距離 平 方 加 權 ( 簡 稱IDSW)算法。 因此,在基于空間插值的反距離平方加權算法下,提出機器人位姿誤差空間 IDSW 插值補償方法,上述所建立的補償模型即為基于均勻數據場的空間 IDSW 插值補償模型。

        4  仿真驗證

        為驗證空間 IDSW 插值算法對機器人位姿誤差補償的有效性,將利用目前較為成熟的微分誤差補償法與本文所提出的方法進行對比試驗。 微分誤差補償法即基于位姿誤差模型,將辨識得到的參數偏差來修正名義參數,來達到誤差補償效果。

        在本文第 2 節中,以 Staubli TX60L 機器人為例已建立了運動學位姿誤差模型,然后對辨識雅克比矩陣進行 QR 分解,剔除兩個無法辨識的參數:dZB和 δZB, 再使用最小二乘法對其余 28 個參數進行辨識,結果如表 2 所示。 最后將該辨識得到的參數偏差用于基于微分誤差補償法補償實驗。
 
表 2  最小二乘法參數辨識結果
 
 
        基于空間 IDSW 插值算法的機器人位姿誤差補償方法首先需要建立合適的空間網格,將 StaubliTX60L 機器人的某局部工作空間劃分為均勻網格,令網格邊長為 100 mm,以機器人末端的笛卡爾坐標系來劃分,X 軸向工作范圍為(300 mm~500 mm),Y軸向工作范圍為( -400 mm ~ 400 mm),Z 軸向工作范圍為( -100 mm ~ 200 mm),具有 48 個網格,108個網格頂點。仿真實驗如下:利用 MATLAB 建立基于 MDH模型的 Staubli TX60L 機器人運動學模型,并編寫空間 IDSW 插值算法和遞歸搜索算法;選取 3 個不同姿態,分別在這 3 個姿態下,在網格空間中任取 3 個點,共 9 個選取的驗證點如表 3 所示。 利用表 1 預定義的運動參數的偏差值,修改運動學模型的參數,分別計算這 9 個點所對應的位置誤差和姿態誤差,結果分別見圖 8和圖 9。
 
表 3  選取的驗證點
 
 
 
 
圖 8  補償前位置誤差 
 
圖 9  補償前姿態誤差 
 
        分別計算 108 網格頂點在如上 3 個姿態下所對應的位置誤差和姿態誤差,所得誤為采集點數據集,用以建立空間 IDSW 插值算法所需的均勻數據場;分別根據微分誤差補償法(圖中簡稱微補法)和空間 IDSW 插值算法(圖中簡稱 IDSW 法)對機器人位姿誤差進行補償,其結果如圖 10 和圖 11 所示。
 
 
圖 10  微補法和 IDSW 法補償后的位置誤差 1 和 2
 
  
 圖 11  微補法和 IDSW 法補償后的姿態誤差 1 和 2

        根據 MATLAB 仿真實驗結果,在 3 個不同姿態下, 微 分 誤 差 補 償 法 補 償 后 位 置 精 度 平 均 提84.8%,83.1%,59.4%,姿態精度平提高 95.9%,96.6%,94.6%;而經反距離平方比加權算法補償后位置精度平均提高 98.0%,98.5%,93.5%,姿態精度平均提高 98.7%,98.5%,97.3%。 所以采用微分誤差補償法和空間 IDSW 插值算法對位姿誤差都具有較好的補償效果。 分析圖 10 和圖 11 可知,空間IDSW 插值算法對位置補償效果優于微分誤差補償法,同時姿態補償效果也優于微分誤差補償法。 因此,所提出的補償方法補償效果更好,且不需要建立復雜的機器人位姿誤差模型,不需要進行繁瑣的參數辨識過程。 綜上所述,經 MATLAB 的對比實驗仿真驗證,在只考慮運動學參數偏差引起的位姿誤差情況下,所提出的插值算法對機器人位姿誤差補償是有效的。

        5  結論

        在探討機器人誤差補償的現有方法之上,研究基于空間插值的補償方法,并提出機器人末端位置與機器人末端位姿誤差的相關性程度是決定插補補償方法是否有效的充分條件。 首先利用 Staubli TX60L 的DH 參數建立運動學模型,并預定義各項運動參數的偏差值,使用 MATLAB 建模仿真分析,得到當機器人末端姿態不變時,任意兩點對應的位置誤差之間的差值和姿態誤差之間的差值分別與兩點間的位置差值具有高度線性關系。 然而機器人末端任意姿態時,兩點的位置誤差的差值和姿態誤差的差值分別與兩點間的位置差值之間幾乎無關。 因此,提出了一種姿態變化因素下的定位位姿誤差模型。 最終提出機器人位姿誤差空間 IDSW 插值補償方法,并由 MATLAB 仿真驗證了該補償方法的有效性。

        所提出的位姿誤差補償方法最大的優點是不需要建立機器人位姿誤差模型,也不需要辨識機器人復雜、冗余的各項參數,在測量上也只需要有限的采集點,尤其是對在機器人操作時末端執行器具有任意姿態且對姿態要求高的機器人應用領域具有較好的應用前景,
該方法具有計算過程簡單、速度快、補償效果好等優點。 基于空間 IDSW 插值算法是一種空間內插方法,只能對網格空間內的數據進行補償,所以需探索擴展方法,研究具有外推能力的插值方法將是后續的研究工作。
 
    投稿箱:
        如果您有機床行業、企業相關新聞稿件發表,或進行資訊合作,歡迎聯系本網編輯部, 郵箱:skjcsc@vip.sina.com
主站蜘蛛池模板: 欧美专区中文字幕 | 天天舔天天爽 | 亚洲精品久久久一区二区三区 | 久二影院| 欧美精品成人一区二区三区四区 | 成人av影院 | 欧美国产视频 | 亚洲一区二区三区日韩 | 一区小视频 | 一区二区三区国产好 | 久久国产亚洲 | 免费看的毛片 | 亚洲综合无码一区二区 | 天天干国产 | 欧美精品一区二区在线观看 | 99re6在线视频精品免费 | 羞羞羞网站 | 中文字幕在线网址 | 日韩精品在线一区 | av电影手机版 | 欧美黄色一区 | 天天干干干干 | 91精品久久久久久久久久入口 | av网站网址 | 你懂的在线视频播放 | 国产一区二区三区免费 | 色av色av色av | 日韩欧美在线观看一区二区三区 | 日韩国产欧美视频 | 国产成人av电影 | 亚洲一区二区三区 | 久久视频国产 | aaa在线观看 | 综合久久久| 久久精品久久久久久久久久久久久 | 欧美日韩一区二区中文字幕 | 欧美日韩高清在线一区 | 久久综合91 | 久久久久久久久久影院 | 欧美一区二区在线播放 | 久草视频在线看 | 亚洲精品二三区 | 久久黄网 | 中文字幕在线免费播放 | 可以在线观看的黄色 | 日韩中文一区 | 欧美一级视频在线观看 | 午夜激情视频在线观看 | 国产精品人人做人人爽人人添 | 精品乱子伦一区二区三区 | 激情一区| a在线看 | 日韩中文字幕在线观看 | 久久三区| 欧美不卡视频 | 国产成人av在线 | 亚洲视频在线观看 | 91资源在线观看 | 国产精品第一国产精品 | 精品视频在线观看一区二区三区 | 日韩成人免费av | 午夜精品久久 | 国产一区二区三区免费看 | 色精品| 在线国产一区 | 四虎影院观看 | 毛片网站免费在线观看 | 久久夜精| 韩国精品一区二区三区 | 免费国产一区二区 | 狠狠插狠狠操 | 欧美日韩中文 | 国产三区在线成人av | 亚洲久久 | 久久精品国产99国产精品 | 久久国产精品精品国产 | 亚洲精品成人 | a级在线 | 久久精品不卡 | 我看一级毛片 | 91精品国产一区二区三区四区在线 | 精品国产乱码久久久久久闺蜜 | 欧美日一区 | 免费av电影观看 | 国产亚洲精品美女久久久久久久久久 | 精品久久中文 | 午夜老湿影院 | 91亚洲日本| 久久综合久 | 国产一区中文字幕 | 欧美一级片免费看 | 国产精品美女久久久久久久网站 | 久久高清国产 | 狠狠入ady亚洲精品经典电影 | 国产精品久久久久久久久久免费 | 欧美狠狠操 | 91久久国产综合久久蜜月精品 | 久久精品无码一区二区日韩av | 一级毛片电影 | 91精品国产综合久久久久久软件 | 中文一二区 | 久久精品亚洲精品 | 色视频网站在线观看 | 午夜精品久久久久久 | 国产偷自视频区视频 | 日韩在线观看中文字幕 | 欧美视频网站 | 日韩日b视频 | 中文字幕在线三区 | 百性阁综合另类 | 成人福利在线 | 99中文字幕 | 国产精品99一区二区三区 | 久久精品视频免费看 | 一区二区影视 | 午夜私人影院 | 在线免费国产视频 | 国产精品久久视频 | av国产精品| 欧美日韩在线精品 | 久久国产日韩 | 一区二区精品视频在线观看 | 亚洲动漫在线观看 | 国产成人福利在线 | 在线成人av | 91精品国产综合久久国产大片 | 99视频精品在线 | 免费视频一区二区 | 精品久久久久久久久久久 | 日韩高清中文字幕 | 久久av网| 亚洲综合欧美日韩 | 亚洲国产情侣自拍 | 欧洲av在线 | 五月婷婷丁香 | 久久亚洲综合 | 久久综合九九 | 亚洲精品电影在线一区 | 成人一区二区三区在线观看 | 四虎影院最新地址 | 伊人超碰 | 日韩在线视频一区 | 九九九视频精品 | www.成人在线视频 | 国产免费久久 | 日韩高清国产一区在线 | 日本福利网站 | 日本成人久久 | 成人精品鲁一区一区二区 | 国产精品久久电影观看 | 中文字幕av一区二区 | 精品久久影院 | 亚洲精品3区 | 精品久久久久久久人人人人传媒 | 狠狠操操 | 国产精品久久久精品 | 99在线观看 | 天天干天天操 | 美女黄网| 免费观看一级特黄欧美大片 | 久久久久久91香蕉国产 | 亚洲精品日韩在线 | 午夜精品福利网 | 色综合久久一区二区三区 | 精品国产一区一区二区三亚瑟 | 精品国产一区二区三区久久久 | 欧美日韩国产在线观看 | 精品一区二区三区三区 | 欧美成人三区 | 二区中文字幕 | 性开放xxxhd视频 | 最新日韩精品在线观看 | 成人激情视频在线观看 | 欧美区视频 | 黄色片免费看 | 欧美久草 | 亚洲狠狠爱一区二区三区 | 国产欧美精品一区二区三区 | 日韩亚洲精品在线观看 | 狠狠狠色丁香婷婷综合久久五月 | 国产成人精品在线 | 日韩精品视频免费专区在线播放 | 91大片| 午夜一区二区三区 | 日韩久久精品一区二区 | 狠狠干狠狠操 | 色综合99 | 亚洲欧美国产精品久久久久 | 在线观看毛片网站 | 日本三级2018 | av福利网站 | 一级一片免费视频 | 久草色视频在线观看 | 日韩爱爱视频 | 日韩免费网站 | 成人精品免费视频 | 99精品免费在线 | 国产中文一区 | 欧美性猛交一区二区三区精品 | 免费成人在线观看视频 | 免费视频久久 | 久久久国产精品 | 日本在线观看 | 亚洲一区日韩 | 日韩成人精品在线 | 欧美精品1区 | 国产欧美精品区一区二区三区 | 国产成人精品一区二区三区视频 | 99久久婷婷国产精品综合 | 亚洲电影一区二区三区 | 久久黄色片 | 久久精品一区二区三区不卡牛牛 | 欧美国产精品一区二区三区 | 一区二区三区中文字幕 | 黄色免费观看 | 日韩免费av | 草比网站 | 91av导航 | 国产精品久久久久久一区二区三区 | 久久精品久久久久电影 | 亚洲热在线视频 | 亚洲精美视频 | 精品久久一区二区 | 在线中文字幕观看 | 五月香婷婷 | 成人高清视频免费观看 | 欧美午夜一区二区三区免费大片 | 最新免费视频 | 日韩精品一区二区三区中文在线 | 青青草在线视频免费观看 | 国产成人精品久久二区二区91 | 欧洲国产伦久久久久久久 | 中文字幕在线免费看 | 欧美一区二区三区爽大粗免费 | 北条麻妃国产九九九精品小说 | 一区二区在线电影 | 亚洲在线一区二区 | 国产精品永久久久久久久久久 | 欧美日韩在线免费观看 | 91色在线观看 | 免费成人在线网站 | 99精品九九 | 欧洲亚洲视频 | 国产精品777一区二区 | 中文字幕在线精品 | 久久久精品国产 | 成人在线看片 | 亚洲精品在线播放 | 亚洲久草视频 | 日韩欧美国产精品 | 亚洲欧美网址 | 国产成人黄色 | 精品国产99 | 中文二区| 日韩一区二区中文字幕 | av在线一区二区 | 亚洲视频中文字幕 | 欧美日韩国产一区二区三区 | 精品一区视频 | 婷婷国产在线观看 | 天天操夜夜操免费视频 | 一区二区在线看 | 亚洲成人精品久久久 | 人人澡人人草 | 狠狠狠色丁香婷婷综合久久五月 | 久久久精 | 欧美精品自拍 | 中文字幕_第2页_高清免费在线 | 国产超碰在线 | 欧美精品成人一区二区在线 | 久久久一区二区三区 | 国产日韩欧美 | 国产视频二 | 精品久久久久久国产 | 国产精品一区av | 9色porny自拍视频一区二区 | www久久99| 日本不卡免费新一二三区 | 中文字幕第100页 | www亚洲成人| 黄色一级大片在线免费看产 | 美女黄网| 国产精品乱码一区二区三区 | 日韩中文字幕无码一区二区三区 | 久久亚洲一区二区三区四区 | 在线观看国产视频 | 日本精品在线观看 | 三级成人在线 | 狠狠的日 | 久久久影视 | 精品在线二区 | 中文字幕 在线观看 | 免费大片在线观看网站 | 国产精品欧美日韩 | 欧美黄视频在线观看 | 亚洲一区中文字幕 | 久久毛片| 动漫泳衣美女 | 成人一区视频 | 欧美精品一区二区三区四区五区 | 欧美久久久久久久久久久久久久 | 亚洲一区二区三区在线 | 伊人网视频在线观看 | 欧美涩涩网| 伊人av在线免费观看 | 国产精品免费视频一区 | 日韩欧美国产一区二区 | 色呦呦网站在线观看 | 91成人免费视频 | 精品一区二区在线观看 | av一区二区三区 | 欧美日韩中文 | 日韩一区电影 | 久久先锋 | 看亚洲a级一级毛片 | 6080yy精品一区二区三区 | 久久亚洲综合 | 视频精品一区 | 日韩色在线 | 毛片网站在线 | 日本免费三片免费观看 | 欧美亚洲91 | 特级淫片女子高清视频在线观看 | 日韩在线精品强乱中文字幕 | 亚洲电影在线观看 | 有码在线 | 日韩超碰 | 亚洲成人av | 国产一级免费在线 | 久久久精品 | 午夜精品一区二区三区在线 | 超级乱淫片国语对白免费视频 | 观看av | 国产精品毛片久久久久久 | 久久久久久9 | 久久免费视频观看 | 97久久久 | 久久久久国产精品视频 | 国产精品乱码一二三区的特点 | 天天操综 | 一区二区av在线 | 日本成人中文字幕 | 亚洲一级毛片 | 国产一区二区日韩 | 一级毛片免费观看 | 欧一区二区| 久久性 | 国产日韩欧美激情 | 国产高清在线a视频大全 | 成人精品在线视频 | 日本一区二区不卡视频 | 日日精品| 91社区在线高清 | 国产精品福利在线观看 | 91成人在线免费视频 | 久久国产精品免费一区二区三区 | 自拍视频在线观看免费 | 91精品久久久久久久久中文字幕 | www.日韩.com | 91亚色| 亚洲 欧美 精品 | 国产精品99久久 | 久久男人天堂 | 国产在线一区二区 | 午夜精品成人一区二区 | 日韩视频中文字幕在线观看 | 日韩成人免费中文字幕 | 欧洲美女性开放视频 | 亚洲va欧美va天堂v国产综合 | 亚洲一区二区三区四区五区中文 | 国产成人精品一区二区在线 | 麻豆专区一区二区三区四区五区 | 呦呦在线观看 | 精品国产伦一区二区三区观看说明 | 一区国产精品 | 亚洲欧美另类久久久精品2019 | 国产成人一区二区三区影院在线 | 久久99精品久久久久久园产越南 | 日韩亚洲一区二区 | 日本不卡高字幕在线2019 | 精品午夜久久 | 国产视频亚洲 | 三级无遮挡污在线观看 | 欧美一级片毛片免费观看视频 | 国产伦精品一区二区三区照片91 | 视频一区二区三区在线观看 | 在线免费观看羞羞视频 | 中文字幕在线观看1 | 91精品国产人妻国产毛片在线 | 国产精品大片在线观看 | 亚洲天天 | 久久这里只有精品首页 | 欧美中文字幕在线 | 美女福利网站 | 日韩中文字幕一区二区 | av片在线观看网站 | 国产亚洲精品久久久久久久 | 成人二区 | 国产成人精品久久 | 伊人超碰在线 | 99久久99热这里只有精品 | 午夜精品网站 | 亚洲精品一区二三区不卡 | 亚洲一区在线日韩在线深爱 | 亚洲av毛片一级二级在线 | 欧美激情视频一区二区三区在线播放 | 色婷婷av一区二区三区软件 | 日韩成人免费在线 | 中文字幕在线影院 | 国产成人精品无人区一区 | 一级片在线播放 | 国产精品亚洲第一 | aⅴ色国产 欧美 | 亚洲欧美一区二区三区在线 | 亚洲精品久久久久久一区二区 | 日韩视频在线观看 | 天天澡天天狠天天天做 | 国产一区国产二区在线观看 | 欧美一级视频免费 | 午夜精品网站 | 国产真实精品久久二三区 | 亚洲黄页 | 日本中文字幕一区 | 亚洲免费人成在线视频观看 | 中文日韩在线 | 免费一区 | 亚洲三级视频 | 性xxxxxxxxx18欧美 | 超碰在线影院 | 日韩专区在线播放 | 久二影院| 国产高清av在线一区二区三区 | 欧美 日韩 亚洲 一区 | 久久精品99国产精品亚洲最刺激 | 99久热在线精品视频观看 | 精品国产乱码久久久久久久软件 | 日本中文一区二区 | 99久久精品一区二区成人 | 亚洲久久| 久久久久成人精品 | 欧美啊v | 欧美日韩视频在线播放 | 成人免费淫片视频观 | 色综合久久久久 | 亚洲精品视频在线免费播放 | 成人午夜影院 | 日本一本视频 | 91视频在线免费观看 | av一道本| 国产精品大片在线观看 | 欧美成人在线网站 | 亚洲网站在线播放 | 成人一区二区三区 | 香蕉久久av一区二区三区 | 亚洲精品一| 免费的一级视频 | 亚洲www.| 久久新| 国产成人精品一区二区在线 | 国产视频一区二区 | 欧美综合久久 | 香蕉久久夜色精品国产使用方法 | 国产日韩欧美一区 | 婷婷成人免费视频 | 国产精品一区二区久久久久 | 91综合网 | 日本一区二区三区四区 | 天天夜夜操 | 免费看特级毛片 | 九色网址 | 日韩久久久久 | 91视频一88av| 日本成人中文字幕在线观看 | 国产黄色在线免费看 | 在线视频国产一区 | 一区在线视频 | 成人在线一区二区三区 | 成人欧美一区二区三区白人 | 久久久久成人精品 | 国产视频久久久 | 99精品欧美一区二区三区综合在线 | 看一级黄色大片 | 在线观看国产视频 | 国产成人看片 | 精品国产欧美一区二区 | 91精品国产综合久久久久久漫画 | 国产精品久久av | 国产日韩一区二区三区 | 精品一区二区久久久久久久网站 | 国产精品久久久久aaaa九色 | 大陆一级毛片免费视频观看 | 亚洲精品久久久久久下一站 | 男人天堂亚洲天堂 | 国产精品一区二区三区在线 | 九九热免费看 | 久久久久久久久久久久久久久久久久久 | 成人二区 | 一区二区在线视频 | 日韩一区二区免费视频 | 亚洲免费在线视频 | 亚洲精品乱码久久久久久蜜桃不卡 | 成人国产精品免费网站 | 极品久久| 色视频免费在线观看 | 久久精品在线 | 国产精品国产三级国产aⅴ入口 | 国产午夜手机精彩视频 | 亚洲综合色自拍一区 | 精品成人免费一区二区在线播放 | 日韩欧美理论片 | 国产日韩免费视频 | 日日夜夜天天 | 欧美日韩亚洲一区二区 | 成年人网站在线免费看 | 日本福利网站 | 日本黄色一级电影 | 9191在线| 亚洲精品wwww | 精品在线一区二区三区 | 中国一级毛片 | 成人午夜剧场 | 精品一区二区av | 久久久精品一区 | 亚洲男人天堂网 | av色伊人久久综合一区二区 | 91夜夜操 | 中文字幕第一页在线视频 | 在线观看国精产品二区1819 | 国产三级久久久久 | 成人伊人| 日韩中文字幕三区 | 在线亚洲不卡 | 亚洲免费在线看 | 中文字幕乱码一区二区三区 | 免费成人在线视频网站 | 国产一级视频 | 香蕉av在线 | 日韩精品在线播放 | 在线第一页 | 久久亚洲高清 | 久久精品一 | 午夜国产影院 | 免费看的毛片 | 久久精品亚洲精品国产欧美 | 91麻豆精品国产91久久久久 | 毛片视频观看 | 久久中文视频 | 黄理论视频 | 91视频18| 欧美精品一区二区在线观看 | 精品人成| 三级av网站 | 欧美偷偷操 | 久草久| 亚洲精品视频在线观看网站 | 男女视频在线观看 | 插插插干干干 | 日韩激情一区二区三区 | 免费看片国产 | 激情久久av一区av二区av三区 | 羞羞视频在线网站观看 | 国产福利在线观看 | 欧美精品日韩 | 日本韩国欧美一区 | 九九视频在线 | 综合久久综合 | 国产欧美日韩一区二区三区 | 国产成人在线一区 | 中文字幕在线三区 | 中文字幕1区 | 久草免费在线 | 国产精品一卡二卡 | 黄网站色大毛片 | 一区二区三区高清 | 波多野吉衣网站 | 在线视频中文字幕 | 日本特黄特色aaa大片免费 | 日韩在线资源 | 91精品国产九九九久久久亚洲 | 成人在线片 | 午夜视频在线 | 97超碰在线免费 | 日韩一区二区三区在线 | 成人影院www在线观看 | 日韩中文视频 | 国产精品久久免费观看spa | 欧美日本韩国一区二区三区 | 久久久精 | 日韩在线中文 | 久久久久久久久久久网站 | 农村妇女毛片精品久久久 | 精品久久久久久亚洲综合网 | 亚洲人人| 久久线视频 | 99热手机在线观看 | 最新超碰 | 欧美在线国产 | 亚洲国内精品 | 在线观看欧美一区二区三区 | 国产免费观看一区二区三区 | 中文字幕在线观看一区二区三区 | 日韩免费网站 | 国产精品欧美一区二区三区 | 亚洲国产高清在线 | 日韩视频―中文字幕 | 国产男女做爰免费网站 | 国产精品久久久久久久久 | 亚洲免费视频一区 | 精品亚洲一区二区三区四区五区 | 国产欧美精品一区二区三区 | 国产在线精品一区 | 正在播放国产精品 | 日日插日日操 | 久草.com| 欧美美女黄色网 |