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

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

用戶名:   密碼:         免費注冊  |   申請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
主站蜘蛛池模板: 国产日产欧美a级毛片 | 男人超碰| 精品亚洲一区二区三区四区五区 | 久久青草国产 | 欧美精品综合 | 国产精品美女久久久久久免费 | 精品国产不卡一区二区三区 | 亚洲看片 | 国产日韩一区 | 欧美 日韩 中文字幕 | 国产精品永久免费自在线观看 | 国产乱叫456 | 国产高潮在线观看 | 性开放xxxhd视频 | 人人爽日日爽 | 天天夜操 | 国产精品成人网 | 欧美黄色片免费观看 | 在线草| 男人天堂亚洲 | 欧美国产日韩一区 | 无码日韩精品一区二区免费 | 国产精品一区人伦免视频播放 | 国产不卡视频在线观看 | 精品久久久久av | 午夜免费视频网站 | 成人免费视频在线观看 | 成人精品视频在线观看 | 亚洲网站在线播放 | 亚洲成人高清 | 黄色国产 | 亚洲天堂成人在线 | 国产精品久久久久久福利一牛影视 | 成人精品鲁一区一区二区 | 精品无码久久久久久国产 | 欧美a在线 | 国产成人啪精品午夜在线观看 | 欧美另类一二三四 | 午夜影院在线观看版 | 亚洲成人影院在线观看 | 久久综合久久久 | 毛片国产 | 成人黄色一级网站 | 啪一啪操一操 | 最近免费中文字幕在线视频2 | 日韩一及片| 欧美精品一级二级 | h在线看| 欧美日韩在线精品 | 国产一区不卡视频 | 美女精品视频在线 | 精品视频一区二区 | 亚洲国产精品久久久男人的天堂 | 玖玖国产精品视频 | 国产中文字幕在线 | 男人天堂999| 亚洲一区二区在线 | 亚洲综合成人网 | 久久在线播放 | 91视频三区 | 日韩欧美国产成人一区二区 | 免费成人激情视频 | 国产高清在线 | 国产色婷婷精品综合在线播放 | 午夜激情免费在线观看 | 自拍视频在线 | 91精品久久久久久综合五月天 | 日韩av不卡在线 | 久久久久久久精 | 激情久久av一区av二区av三区 | 久草免费在线视频 | 久久国产亚洲精品 | 日韩视频精品在线 | 日韩视频精品 | 国产精品中文字幕在线播放 | 国产玖玖 | 国产99一区二区 | www中文字幕在线观看 | 国产美女精品视频免费观看 | 日韩欧美国产精品综合嫩v 国产高清av在线一区二区三区 | 91在线播 | 欧美日韩在线视频一区 | 午夜欧美一区二区三区在线播放 | 久久成人国产精品 | 在线免费中文字幕 | 999免费视频| 中文字幕av亚洲精品一部二部 | 亚洲成人一区二区 | 国产亚洲精品v | 亚洲人成网站999久久久综合 | 黄在线免费观看 | 精品在线免费视频 | 污视频在线观看免费 | 国产精品一区二区久久 | 中文字幕日韩一区二区三区 | 久久这里精品 | 久久福利 | 黄色网址av | 日韩成人在线一区 | www.精品| 99热新| 国产美女中出 | 色播99| 亚洲国产一区视频 | 欧美视频在线免费 | 色视频www在线播放国产人成 | 欧美成人精品激情在线观看 | 欧美精品三区 | 一级毛片免费看 | 午夜男人天堂 | 日韩免费高清视频 | 国产91久久精品一区二区 | 亚洲综合色自拍一区 | 亚洲精品视频免费观看 | 日本中文字幕一区二区 | 91久久精品一区 | 波多野结衣中文字幕在线视频 | 成人免费看片 | 国产成人精品免高潮在线观看 | 国产91在线观看 | 日韩一二三四 | 亚洲激情视频在线观看 | 成人免费视频网址 | 毛片网站免费在线观看 | 日本不卡高字幕在线2019 | 日韩一级免费在线观看 | 色综合色 | 91精品免费在线观看 | 欧美a在线看 | 色视频久久 | 亚洲成人免费在线观看 | 久久精品99国产精品亚洲最刺激 | 天堂综合网久久 | 一级片在线观看 | 国产精品二区三区 | 99成人在线视频 | 免费毛片网 | 欧美爱爱视频 | 成人久久 | 久久精品视频网站 | 国产日韩精品一区二区 | 国产一区二区三区不卡在线观看 | 成年人在线观看 | 久久成人综合网 | av在线天堂 | 青春草在线观看 | 日韩一区二区三区精品 | 一区二区亚洲 | 亚洲精品亚洲人成人网 | 日韩在线视频观看 | 97视频久久久 | 国产精品免费久久 | 在线99 | 99热精品在线 | 精品欧美日韩 | 草久久av| 中国妞xxxhd露脸偷拍视频 | 成人影音 | 国产一区二区三区视频在线观看 | 国产福利片在线观看 | caoporon| 成人黄页在线观看 | 国产成人免费在线 | 国产精品视频免费看 | 国产精品福利网站 | 色网在线 | 欧美日韩成人在线视频 | 成人免费高清视频 | 日日爽| 国产成人精品一区二区三区 | 国产三级视频 | 久久九九精品久久 | 日韩成人tv | 一区二区日韩视频 | 91社区在线观看 | 91精品国产一区二区三区蜜臀 | 国产一区二区综合 | 亚洲精品一区二区三区在线看 | 国产精品不卡一区 | 亚洲国产aⅴ成人精品无吗 国产精品永久在线观看 | 一级片日韩 | 成人国产在线 | 黄色一级网站 | 亚洲一区二区三区高清 | 中文字幕亚洲一区二区三区 | 国产精品视频免费看 | 在线播放黄色片网站 | 精品色区| 一区二区免费 | 成人激情视频在线免费观看 | 精品国产鲁一鲁一区二区在线观看 | 精品久久久久一区二区国产 | 亚洲精品免费观看 | 99久久国产综合精品女不卡 | 国产精品激情在线观看 | 狠狠干av | av影音 | 成人精品 | 亚洲成人在线观看视频 | 日韩成人在线观看 | 久久国产视频精品 | 久热官网 | 91福利网站在线观看 | 日本一区二区中文字幕 | 热久久这里只有精品 | 亚洲精选久久 | 欧美成人精品一区二区三区 | 亚洲久草视频 | 国产精品黄网站在线观看 | 欧美三级视频在线播放 | 日韩中文字幕一区二区高清99 | 国内久久精品视频 | 精品视频网 | 欧美伦理一区二区 | 久久国产精彩视频 | 国产欧美精品一区二区 | 国产一区二区三区在线免费观看 | 欧美男人天堂 | 一区二区精品在线 | 国产精品成人在线观看 | 免费午夜电影 | 超碰在线一区二区三区 | 久久91av | 午夜免费剧场 | 热久久这里只有精品 | 久久九| 激情综合网五月婷婷 | 欧美一区二区三区精品 | 在线观看日韩av | 二区在线视频 | 中文字幕av亚洲精品一部二部 | 色天天综合久久久久综合片 | 9999久久久久 | 欧美顶级毛片在线播放 | 午夜精品久久久久久久 | 日韩一区免费在线观看 | 成人精品一区二区三区 | 国产在线播 | 欧美日韩综合一区 | 欧美人体一区二区三区 | 国产高清一区二区三区 | 91视频网址| 日韩精品一区二区三区在线播放 | 精品一区二区三区四区五区 | 国产精品久久久久久久久久久久 | 久草视频在线资源站 | 国产中文字幕在线观看 | 青青久久 | 日日骚视频 | 天天操天操 | 亚洲成人av片 | 国产精品粉嫩白浆在线观看 | 日干夜干天天干 | 日韩特级 | 亚洲国产精品久久 | 久久综合亚洲 | 密桃av| 黄色大片免费网址 | 黄色小视频在线免费观看 | 国产亚州av | 免费的日批视频 | 天堂久久久久久 | 国产亚洲在线 | 中文一区| 亚洲欧洲日韩在线 | 一区二区三区国产 | 水卜樱一区二区av | 国产免费一区二区 | 日韩专区一区二区三区 | 成人精品鲁一区一区二区 | 麻豆精品一区二区 | 国产浪潮av色综合久久超碰 | 在线日韩视频 | 亚洲永久免费 | 国产丝袜一区二区三区免费视频 | 色玖玖综合 | 免费国产视频在线观看 | 国产精品69久久久久水密桃 | 亚洲欧洲一区二区 | 日本成人中文字幕 | 91麻豆精品国产91久久久久久 | 欧美成人激情 | 69久久99精品久久久久婷婷 | 亚洲高清免费 | 亚洲在线免费观看 | 午夜爽视频 | 欧美成在线观看 | 久久久久久成人 | 在线免费观看色视频 | 国产成人综合av | 操操操小说 | 欧美久久久久久 | 国产成人精品免费 | 毛片综合| 一区视频在线 | 一卡二卡久久 | 久热官网 | 亚洲欧美一区二区三区在线 | 成人一区av | 欧美第一视频 | 日韩在线精品视频 | 亚洲在线视频 | 久久久久久久国产 | 99久久精品一区二区 | 亚洲免费影院 | 国产精品视频免费播放 | 欧美亚洲91 | 国产伦精品一区二区三区在线 | 亚洲一区二区三区爽爽爽爽爽 | 欧美成人小视频 | 亚洲综合首页 | 日韩一级在线免费观看 | 日本天天操 | 日韩中文字幕一区二区 | 欧美精品一区在线 | 在线免费黄色 | 99久久久久 | 91观看| 亚洲一区二区三区久久久 | 91久久精品国产91久久性色tv | 欧美xxxx性| 国产精品日韩欧美 | 免费毛片网站 | 日本三级欧美三级 | 国产精品久久久久久久午夜 | 日韩成人免费视频 | 一区二区免费在线播放 | 免费中文字幕 | 中文字幕综合 | 亚洲精品日韩精品 | 在线国产视频 | av在线免费观看网站 | 亚洲欧美在线免费观看 | 欧美精品在线视频 | 成人免费久久 | 色综合社区 | 国产亲子乱弄免费视频 | 精品视频一区在线观看 | 亚洲欧美日韩另类精品一区二区三区 | av中文字幕在线观看 | 久久精品麻豆 | 亚洲a在线播放 | a在线观看 | 午夜精品久久久久久久男人的天堂 | 成人黄色电影在线观看 | 国产精品一区人伦免视频播放 | 在线观看免费av网站 | 欧美日韩精品亚洲 | 国产精品久久精品 | 亚洲一区视频在线 | 国产精品一二三 | 欧美日韩一区二区中文字幕 | 中文字幕亚洲欧美日韩在线不卡 | 欧美日韩国产综合视频 | 91久久久久久久久久久久久 | 亚洲av毛片 | 91,看片 | 国产激情免费 | av观看免费| 女人夜夜春高潮爽a∨片传媒 | 欧美一级视频免费 | 亚洲a网 | 玖玖免费 | 欧美日韩国产一区二区 | 一级片av | 99精品国产高清在线观看 | 亚洲精品一区久久久久久 | 国产精品theporn| av福利网站 | 日韩成人精品在线 | 高清一区二区三区 | 国产激情在线看 | 日本毛片视频 | 欧美黄视频在线观看 | 操人网| 精品一二区 | 亚洲国产精品精华液com | 日韩精品视频免费在线观看 | 国产成人一区二区 | 精品一二三区 | 精品免费视频 | 欧美狠狠操 | 亚洲免费在线视频 | www国产成人免费观看视频,深夜成人网 | 99精品99| 精品视频在线免费观看 | 91xx在线观看 | 亚洲精品在线观看免费 | 男人的天堂在线视频 | 中文字幕亚洲精品 | 黄色一级片看看 | 欧美成人免费 | 日韩亚洲 | 丝袜久久| 2021狠狠干| 国产精品乱码一区二区三区 | 麻豆一区二区三区 | 国产亚洲精品美女久久久久久久久久 | 亚洲成av| 国变精品美女久久久久av爽 | 久久国产精品久久久久久电车 | 91精品国产综合久久久久久丝袜 | 蜜臀91精品国产高清在线观看 | 欧美精品一级 | 成人 在线 | 国产情侣在线视频 | 色视频网站免费看 | 欧美黄色一区 | 婷五月综合| 91精品国产综合久久久蜜臀图片 | 一区二区av | 精品日韩一区二区三区 | 久久精品国产一区二区三区不卡 | 国产一区二区视频在线播放 | 国产欧美日韩 | 中文字幕三区 | 在线视频a| 伊人av在线免费观看 | 成人午夜在线视频 | 日本免费xxxx| 欧美寡妇偷汉性猛交 | 久久视频免费在线 | 国产精品久久久久久久 | 中国电影黄色一级片免费观看 | 欧美日韩亚洲一区二区 | 天天操天天曰 | 欧美日韩高清 | 成人免费视频网站 | 自拍偷拍小视频 | 欧日韩不卡在线视频 | 精品久久久久久久久久久久包黑料 | 国产精品国产精品国产专区不蜜 | 日韩成人免费 | 91精品久久 | 精品成人av | 国产精品久久久精品 | 一级做a毛片 | 亚洲一区中文字幕 | 国产96在线观看 | 国产www在线| 日韩欧洲亚洲 | 久久精品一区二区三区四区 | 91亚洲精品久久久 | 精品国产色 | 久久小草 | 黄色91在线 | 亚洲国产精品久久 | 久久黄网| 国产在线一级片 | 国产1级片 | 日韩欧美一区二区视频 | 日韩精品视频三区 | 亚洲精选久久 | 激情图区在线观看 | 古典武侠第一页久久777 | 精品久久精品 | 国产艳妇av视国产精选av一区 | 精品免费视频 | 午夜三级在线 | 国产成人av在线播放 | 久久久久久久国产精品 | 国产欧美一区二区精品忘忧草 | 夜夜爽99久久国产综合精品女不卡 | 久久一区国产 | 久久国产成人 | 欧美成人a∨高清免费观看 在线视频成人 | 一区二区三区精品 | 成人小视频在线看 | 一区二区三区日韩在线 | 美女操网站 | 成人欧美一区二区三区 | 国产中文字幕一区二区三区 | 精品少妇一区二区三区日产乱码 | 国产日韩欧美在线观看 | 国产一区二区在线看 | 久久99精品久久久久久久青青日本 | 日韩特级 | 天天舔夜夜 | 久视频在线观看 | 亚洲精品乱码久久久久久蜜桃 | 韩国av一区二区 | 日本aaaa| 成人网av | 在线免费日韩 | 欧美亚洲另类丝袜综合网动图 | 在线日韩 | 欧美日韩精品在线观看 | 亚洲国产高清高潮精品美女 | 国产精品成人国产乱一区 | 亚洲精品久久久一区二区三区 | 免费在线亚洲 | 91看片淫黄大片一级在线观看 | 狠狠干美女 | 亚洲精品视频在线播放 | 欧美男人的天堂 | 国产成人毛片 | 精品免费国产一区二区三区 | 国产精品中文字幕在线 | 欧美a网 | 精品在线二区 | 中文字幕在线看 | 久久久久久久久国产 | 麻豆毛片 | 高清三区 | 深夜av在线 | 久久亚洲综合 | 一区二区三区成人久久爱 | 欧美亚洲一 | 中文字幕亚洲精品 | 精品99在线 | 激情国产| 午夜电影网址 | 91在线精品一区二区 | 永久免费网站 | 欧美精品一区在线发布 | 国产成人久久精品一区二区三区 | 99精品亚洲 | 国产精品a一区二区三区网址 | 亚洲欧美日韩另类精品一区二区三区 | 欧洲亚洲视频 | 中文字幕一区二区三区乱码图片 | 亚洲乱码一区二区三区在线观看 | 国产一区二区三区免费视频 | 中文字幕在线看 | 国产美女一区二区 | 91文字幕巨乱亚洲香蕉 | 91麻豆精品国产91久久久久久久久 | 国产精品美女久久久久久久久久久 | 精品久久久久久久久久久久久久 | 成人片免费看 | 久久麻豆视频 | 国产不卡在线观看 | 成人xxx| 狠狠操狠狠操 | 中文字幕高清 | 日韩欧美国产精品综合嫩v 久久久久久国产精品高清 国产目拍亚洲精品99久久精品 | 四虎黄色网 | 日本日韩中文字幕 | 曰批免费视频播放免费 | 亚洲一区二区在线 | 精品无码久久久久国产 | 欧美精品在线免费观看 | 黄网站在线播放 | 伊人久久视频 | 成人网址在线观看 | 日韩精品在线观看视频 | 丁香五月亚洲综合在线 | 91精品久久久久久久 | 欧美,日韩,国产精品免费观看 | 在线无码 | 国产中文在线 | 天堂在线www| 天天综合91| 视频一区二区三区在线播放 | 欧美黄视频在线观看 | 日本高清视频一区二区三区 | 中文字幕本久久精品一区 | 亚洲一区日韩 | 五月天婷婷色综合 | 亚洲美女精品视频 | 国产一区二区视频在线观看 | 久久一道本| 成人午夜精品一区二区三区 | 日韩视频在线观看 | 中文字幕成人 | 久久国产亚洲精品 | 成人综合在线观看 | 精品国产一区二区三区久久久 | 国产毛片精品 | av中文在线 | 不卡一区二区三区四区 | 99re热精品视频 | 久久一区| 国产免费亚洲 | 电影91 | 日韩国产一区二区 | 欧美亚洲视频 | 天天操操 | 在线成人免费视频 | av网站免费在线观看 | av大全在线 | 亚洲一区二区三区四区五区中文 | 欧美一区二区视频 | 成年人综合网 | 综合精品 | 欧美激情视频久久 | 亚洲精品一区二区网址 | 国产一区二区av | 成人小视频在线观看 | 日韩爱爱视频 | 中文字幕av在线播放 | 国产在线精品一区二区三区 | 精品国产一区二区 | 国产干干干| 免费一区二区三区 | 国产成人精品一区二 | 91免费影片 | 日韩三级在线 | www.久久| 一区二区三区高清不卡 | 精品国产乱码久久久久久1区2区 | 日韩一级免费在线观看 | 日韩一区二区在线观看视频 | 日韩精品无码一区二区三区 | 精久久久 | 亚洲精品午夜aaa久久久 | 99国产精品| 欧美一级片在线观看 | 亚洲欧美综合精品久久成人 | 亚洲电影中文字幕 | 欧美aaaaa | 特级av| 国产精品日本一区二区不卡视频 | 荷兰欧美一级毛片 | 亚洲欧洲精品一区二区 |