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

    大部件對接中基于單位四元數(shù)的iGPS測量位姿比對研究

    2014-03-28 10:39:58徐碧菡孫涪龍邢宏文
    圖學學報 2014年5期
    關(guān)鍵詞:測量點位姿點位

    徐碧菡, 孫涪龍, 趙 罡, 邢宏文

    (1.北京航空航天大學中法工程師學院,北京 100191;2.北京航空航天大學機械工程及自動化學院,北京 100191;3.上海飛機制造有限公司,上海 200436)

    iGPS (indoor GPS)是一種用于大尺寸空間測量與定位的新型三坐標工業(yè)測量系統(tǒng),其工作原理類似于GPS,利用三角測量原理建立三維坐標系,不同的是采用紅外激光代替了衛(wèi)星(微波)信號。該系統(tǒng)具有測量范圍大、精度高、實時性好、系統(tǒng)安裝方便等特點[1],并且支持多用戶,即可以使用多個傳感器來完成不同的測量任務,如部件定位、指導機器人工作及同步跟蹤[2]等,具有廣闊的發(fā)展前景。目前國內(nèi)各大航空航天企業(yè)已經(jīng)開始引進iGPS,但對其在飛機裝配過程中的應用仍然處于探索階段。

    在大尺寸部件數(shù)字化裝配的過程中,運用iGPS進行合理布局,可較為便捷地測得部件關(guān)鍵點的位置信息,此時最為關(guān)鍵的問題即為對所測數(shù)據(jù)進行處理并解算出部件的位姿。但由于iGPS實測數(shù)據(jù)所采用的坐標系統(tǒng)為測量系統(tǒng)定義的局部三維直角坐標系,與理論模型所采用的設(shè)計坐標系進行統(tǒng)一后,實測數(shù)據(jù)點與理論模型并不一致,因此針對iGPS所測數(shù)據(jù)與理論模型進行位姿比對研究是進一步利用iGPS輔助部件位姿調(diào)整的重要前提。

    在測量數(shù)據(jù)與理論模型進行比對的過程中,需要首先實現(xiàn)剛體姿態(tài)解算。目前用于剛體姿態(tài)解算的方法主要有:奇異值分解法、三點法、最小二乘法、七參數(shù)法、單位四元數(shù)法等[3-5]。其中,奇異值分解法雖然求解結(jié)果較為準確,然而其應用前提為假設(shè)質(zhì)點經(jīng)過坐標轉(zhuǎn)換后完全重合,實際在大部件對接中,由于測量誤差的存在,質(zhì)點經(jīng)過坐標轉(zhuǎn)換后不完全重合,從而影響精度;三點法只適用于不共線的3個點,不適用于iGPS同時測量多個關(guān)鍵點的使用條件;最小二乘法精度受初始值影響大[3];七參數(shù)法在實際使用過程中涉及針對較小旋轉(zhuǎn)角度的線性化[4],因此若旋轉(zhuǎn)角度較大,就會產(chǎn)生較大的誤差,不適宜用于大部件對接過程中的剛體位姿描述;單位四元數(shù)法描述剛體的旋轉(zhuǎn)與平移變換直觀明了,可針對多點進行計算[5],并且具有無奇點、不涉及大量三角函數(shù)計算等優(yōu)點,具備作為大部件對接過程中剛體姿態(tài)描述方法的條件。

    本文根據(jù)iGPS本身測量數(shù)據(jù)的特點,對測量數(shù)據(jù)的不同形式進行分析與選擇,并根據(jù)多點同時測量、大旋轉(zhuǎn)角等測量要求,提出一種基于單位四元數(shù)的測量數(shù)據(jù)與理論模型位姿比對計算方法,最后通過實驗對該算法進行了分析與驗證。

    1 大部件對接中的iGPS測量

    iGPS由激光發(fā)射器、傳感器、手持探頭、接收器電路和系統(tǒng)軟件組成。測量時,發(fā)射器產(chǎn)生兩個激光平面在工作區(qū)域內(nèi)旋轉(zhuǎn),轉(zhuǎn)速約為3000 r/min。接收器根據(jù)其所能接收到的激光信號,對水平角與垂直角進行測量[6]。

    iGPS的測量數(shù)據(jù)以Frame為單位。Frame與坐標系類似,每個Frame可包含6個自由度信息,即3個方位自由度信息與3個位置自由度信息。但根據(jù)測量前用戶對Frame中所包含傳感器個數(shù)及相對位置的不同設(shè)定,F(xiàn)rame的自由度信息也會不同。例如,僅包含由單個傳感器返回數(shù)據(jù)的Frame,只包含點位信息,即3個位置自由度信息;包含兩個相距固定距離的傳感器所返回數(shù)據(jù)的Frame具有5個自由度信息;而包含4個相對位置固定的傳感器(且其中至少3個傳感器不共線)所返回數(shù)據(jù)的Frame,則包含6個自由度信息,即由iGPS底層數(shù)據(jù)處理軟件Surveyor擬合得出的Frame整體的位置與姿態(tài)信息。

    使用iGPS進行對大部件上多個測量點進行數(shù)據(jù)測量時,根據(jù)Frame的不同數(shù)據(jù)結(jié)構(gòu),有兩種數(shù)據(jù)形式可以選擇,分別為點位測量與位姿測量。點位測量僅包含3個位置自由度信息。當對多個測量點分別布置單個傳感器,并將每個傳感器所測數(shù)據(jù)設(shè)定為單獨的Frame時,得到的測量數(shù)據(jù)是每個傳感器幾何中心的位置坐標。而位姿測量則需將Frame設(shè)置為由多個相互位置固定的傳感器(且至少3個傳感器不共線)及其所測數(shù)據(jù)構(gòu)成的整體,得到的數(shù)據(jù)是由Surveyor經(jīng)過數(shù)據(jù)融合計算得出的Frame坐標系原點位置和三坐標軸信息,以及各個點在Frame坐標系下的坐標,但無法獲取每個傳感器所測的原始數(shù)據(jù)。

    在實際測量過程中,點位測量與位姿測量分別有各自的優(yōu)點與不足。點位測量形式的優(yōu)點為可直接采用測量點的原始數(shù)據(jù)進行計算,并根據(jù)實際情況有針對性地設(shè)計位姿標定算法;缺點為計算過程可能較為復雜。位姿測量形式的優(yōu)點為可以直接方便地獲取對應Frame的位置與姿態(tài),即測量目標的位置與姿態(tài);缺點為該位姿標定過程由Surveyor自動完成,計算過程涉及多次數(shù)據(jù)融合,用戶難以針對實際測量情況控制位姿標定過程,誤差分析較難進行。

    因此,本文為避免由系統(tǒng)自身多次數(shù)據(jù)融合所帶來的誤差,采用點位測量形式來配置Frame數(shù)據(jù)結(jié)構(gòu),即測量各點的位置坐標信息(3自由度),并以此為基礎(chǔ)設(shè)計適用于飛機對接過程的模型位姿比對算法,為后續(xù)的部件對接提供輸入值。

    2 基于單位四元數(shù)的iGPS測量數(shù)據(jù)位姿比對算法

    將所測點位數(shù)據(jù)導入理論模型設(shè)計軟件后,iGPS測量坐標系與理論模型的全局坐標系重合,便可獲知在同一坐標系下,兩對接部件各自的理論位姿與實際位姿。此時應首先求解部件實測點位信息與理論模型對應點位信息之間的對應關(guān)系,即求解兩組一一對應點之間的坐標轉(zhuǎn)換關(guān)系。根據(jù)前文對目前幾種常用姿態(tài)解算方法的分析,本文采用單位四元數(shù)法描述部件姿態(tài)。

    2.1 基于單位四元數(shù)的坐標轉(zhuǎn)換方法

    設(shè)初始點集為A,目標點集為E,其中目標點與參考點數(shù)目相同、一一對應。若每個點集中點數(shù)為k,則可記A={ai},E={ei},i=1,2,…,k。計算兩組一一對應點的坐標轉(zhuǎn)換關(guān)系,即為求解滿足:

    的最佳旋轉(zhuǎn)矩陣R與平移向量T。

    根據(jù)剛體定點轉(zhuǎn)動理論中的歐拉定理:剛體繞固定點的任一位移,都可由繞過此點的某一軸轉(zhuǎn)過一個角度得到。假設(shè)坐標系OXYZ下該旋轉(zhuǎn)矩陣對應剛體在三維空間中繞單位轉(zhuǎn)軸向量v=(ax,ay,az)旋轉(zhuǎn)θ角,則基于單位四元數(shù)的坐標轉(zhuǎn)換方法思想[6]如下:

    假設(shè)旋轉(zhuǎn)轉(zhuǎn)換向量為單位四元數(shù)

    其中各分量分別為[7]:

    且各分量滿足:

    則根據(jù)單位四元數(shù)可得到3×3的旋轉(zhuǎn)矩陣R(qR)[8-9]:

    設(shè)平移向量為:

    則完整的坐標變換向量即為:

    考慮到飛機部件尺寸大剛性差的特點,難以做到坐標轉(zhuǎn)換后與模型完全重合,根據(jù)式(1),結(jié)合最小二乘原理設(shè)計優(yōu)化函數(shù)如下:

    至此求解對應點集間的最佳坐標轉(zhuǎn)換向量問題轉(zhuǎn)化為求使得f(q)具有最小值的坐標轉(zhuǎn)換向量q。

    分別求出初始點集A與目標點集E的重心:

    由集合A與集合E構(gòu)成協(xié)方差矩陣:

    然后利用反對稱矩陣:

    的輪換對稱項來構(gòu)成列向量:

    并進一步構(gòu)造用于求解最優(yōu)旋轉(zhuǎn)變換向量的矩陣Q(ΣAE)[8]:

    其中,I3為3×3的單位矩陣,tr(ΣAE)為矩陣ΣAE的跡。對應Q(ΣAE)取最大特征值的特征向量qR=[qqqq]T即為最優(yōu)的旋轉(zhuǎn)變換向量。因此,最優(yōu)的平移向量為:

    至此,旋轉(zhuǎn)矩陣與平移向量求解完畢,即可實現(xiàn)坐標轉(zhuǎn)換。旋轉(zhuǎn)矩陣所對應的旋轉(zhuǎn)軸v與旋轉(zhuǎn)角度值θ可結(jié)合式(3)計算得出。

    2.2 iGPS測量數(shù)據(jù)與理論模型比對方法

    假設(shè)理論模型中,給定兩對接部件M、N,其中M為固定部件,N為移動部件。M與N分別對應于實際部件M′與N′。測量時分別在M′與N′上選定并布置m與n個測量點(m≥3,n≥3),m與n之間無相互約束關(guān)系,可視實際測量需求而定。

    首先考慮固定部件,由于其理論模型M在設(shè)計坐標系下所有坐標點的信息已知,則可獲取設(shè)計坐標系下m個測量點的理論坐標信息;實測部件M′上所布置的m個測量點的坐標信息通過iGPS實測獲得。此時便有了兩組一一對應的坐標值,可采用2.1節(jié)中的基于單位四元數(shù)的轉(zhuǎn)換方法計算其轉(zhuǎn)換關(guān)系Q。此處涉及的轉(zhuǎn)換為將理論模型M上的m個點轉(zhuǎn)換到實測部件位置。

    其次考慮移動部件,確定在實際部件N′上的n個測量點后,通過iGPS實測獲取相對應的n個坐標點信息。此時由于根據(jù)第一步中求出的轉(zhuǎn)換關(guān)系Q,即可對移動部件理論模型N上的n個坐標點進行同樣的坐標變換,從而求出實際測量場地中移動部件N′的目標位置N′′。

    至此確定了移動部件的實際位置N′與目標位置N′′,即可再一次運用2.1節(jié)中所介紹的基于單位四元數(shù)的坐標轉(zhuǎn)換方法,計算出移動部件到達目標位置所需的坐標轉(zhuǎn)換關(guān)系P,從而為后續(xù)調(diào)姿步驟提供輸入數(shù)據(jù)。

    該方法的具體實現(xiàn)流程如圖1所示。

    圖1 iGPS測量數(shù)據(jù)與理論模型比對總體方法計算流程

    3 實驗與分析

    3.1 驗證實驗方案

    為驗證本文提出的基于單位四元數(shù)的iGPS測量數(shù)據(jù)位姿比對算法的精確度,設(shè)計驗證實驗。整體思路為,首先根據(jù)兩對接部件實測點位信息與理論點位信息計算相關(guān)坐標轉(zhuǎn)換關(guān)系;接著在理論模型中兩部件對接面上取3個既屬于固定部件又屬于移動部件的點,用于驗證,理想情況即為:計算得到的移動部件目標位姿上的這3個點,與固定部件實測位姿上的這3個點基本重合;于是對這3個點分別應用第一步擬合得出的坐標變換關(guān)系,得出固定部件上此3點的實際位置信息,和移動部件上此3點的目標位置信息,計算這兩組位置信息之間的誤差,考察2節(jié)中所提出算法的精度。

    具體實驗過程如下:

    (1) 首先為實際固定部件與移動部件分別選取m個與n個測量點;

    (2) 根據(jù)2.1節(jié)中坐標轉(zhuǎn)換方法,由固定部件實測點位信息M′與其理論點位信息M,計算M向M′的轉(zhuǎn)換關(guān)系Q;

    (3) 根據(jù)式(3)計算Q中單位四元數(shù)對應的旋轉(zhuǎn)軸vQ與旋轉(zhuǎn)角度θQ;

    (4) 根據(jù)2.1節(jié)中坐標轉(zhuǎn)換方法,由移動部件實測點位信息N'與其理論點位信息N,計算由N向N′的轉(zhuǎn)換關(guān)系S;

    (5) 根據(jù)式(3)計算S中單位四元數(shù)對應的旋轉(zhuǎn)軸vS與旋轉(zhuǎn)角度θS;

    (6) 對移動部件理論點位信息N應用轉(zhuǎn)換關(guān)系Q,得到移動部件目標點位信息N′′;

    (7) 根據(jù)2.1節(jié)中坐標轉(zhuǎn)換方法,由移動部件實測點位信息N′與移動部件目標點位信息N′′,計算固定部件由N′向N′′的轉(zhuǎn)換關(guān)系P;

    (8) 根據(jù)式(3)計算P中單位四元數(shù)對應的旋轉(zhuǎn)軸vP與旋轉(zhuǎn)角度θP;

    (9) 對移動部件實測點位信息N′應用轉(zhuǎn)換關(guān)系P,得到移動部件最終點位信息N′′′;

    (10) 在理論模型中兩部件對接面上取3個驗證點,則這3個點既屬于固定部件,又屬于移動部件。理想條件下對接后固定部件和移動部件上的這3個點應當重合。將此3點位置信息記為MNTest;

    (11) 對固定部件對接面上的3個點應用變換關(guān)系Q,得到固定部件驗證點的實際位置信息MTest;

    (12) 對移動部件對接面上的3個點應用變換關(guān)系S,使其到達移動部件實測位置N′附近,得出移動部件驗證點的實際位置信息NTest;

    (13) 對NTest的3個點應用變換關(guān)系P,得到移動部件驗證點的目標位置信息OTest;

    (14) 對比MTest與OTest兩組驗證點的偏差。計算偏差平均值

    與標準差σTest。

    3.2 實驗數(shù)據(jù)

    下面根據(jù)3.1節(jié)中的實驗方案進行3組實驗。

    為實際固定部件與移動部件分別選取4個測量點,即m=n=4。(實際應用過程中,應根據(jù)實際測量點的選取情況確定m與n的取值。)

    固定部件理論模型上4個對應點位信息M如表1所示。

    移動部件理論模型上4個對應點位信息N如表2所示。

    表2 移動部件理論模型對應點位信息N

    實驗結(jié)果如表3所示。計算過程中涉及的各坐標轉(zhuǎn)換關(guān)系所對應的旋轉(zhuǎn)軸、旋轉(zhuǎn)角度、旋轉(zhuǎn)矩陣以及平移向量結(jié)果也均在表中呈現(xiàn)。

    表3 驗證實驗結(jié)果

    (待續(xù))

    (續(xù)表)

    (待續(xù))

    (續(xù)表)

    3.3 結(jié)果分析

    由實驗結(jié)果可以看出,在第I組與第Ⅲ組實驗中,固定部件的實際點位信息不變,僅改變移動部件的實際點位信息,因此移動部件由理論模型點位信息向?qū)崪y點位信息轉(zhuǎn)動角度以及由實測點位信息向目標點位信息轉(zhuǎn)動角度都不同,平移向量也不同,但最終驗證點之間的誤差都小于10-11數(shù)量級。

    第Ⅱ組實驗與第I組相比,則將固定部件與移動部件的實際位置同時改變,但可以看出,無論是大旋轉(zhuǎn)角還是小旋轉(zhuǎn)角,坐標轉(zhuǎn)換的誤差始終小于10-11數(shù)量級,標準差也小于10-11數(shù)量級。

    4 結(jié) 論

    本文針對大部件對接過程中iGPS多點同時測量、對接旋轉(zhuǎn)角大的特點,提出了基于單位四元數(shù)的iGPS測量數(shù)據(jù)位姿比對算法,用于部件對接過程中測量數(shù)據(jù)與理論模型的位姿比對。該算法描述位姿變換簡潔直觀,計算量小,效率高。誤差計算與驗證實驗結(jié)果表明,該算法對旋轉(zhuǎn)角大小和平移量無特殊要求,且在不同條件下均有很高的轉(zhuǎn)換精度,能有效適用于iGPS的測量位姿比對計算。

    [1]呂景亮,張春富,唐文彥,郝欣宇.Indoor GPS 工業(yè)測量系統(tǒng)自標定技術(shù)研究[J].計量學報,2011,32(1):12-15.

    [2]Norman A R,Sch?nberg A,Gorlach I A,Schmitt R.Validation of iGPS as an external measurement system for cooperative robot positioning [J].The International Journal of Advanced Manufacturing Technology,2013,64(1-4): 427-446.

    [3]羅 芳,鄒 方,周萬勇.飛機大部件對接中的位姿計算方法[J].航空制造技術(shù),2011,(3): 91-94.

    [4]陳 宇,白征東.基于非線性最小二乘算法的空間坐標轉(zhuǎn)換[J].大地測量與地球動力學,2010,30(2):129-132.

    [5]Horn B K P.Closed-form solution of absolute orientation using unit quaternions [J].JOSA A,1987,4(4): 629-642.

    [6]吳曉峰,張國雄.室內(nèi)GPS測量系統(tǒng)及其在飛機裝配中的應用[J].航空精密制造技術(shù),2006,42(5): 1-5.

    [7]Yoon D,Narduzzi M,Shen J.Computing rotation minimizing frames using quaternions [J].Computer-Aided Design and Applications,2012,9(5):679-690.

    [8]Besl P J,McKay N D.A Method for registration of 3-D shapes [J].IEEE Transactions on Pattern Analysis and Machine Intelligence,1992,14(2): 239-256.

    [9]滕志遠,張愛武.單位四元素法在激光點云坐標轉(zhuǎn)換中的應用[J].測繪通報,2010(11): 7-10.

    猜你喜歡
    測量點位姿點位
    1年期國債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    飛機部件數(shù)字化調(diào)姿定位測量點的優(yōu)選與構(gòu)造算法
    航空學報(2022年5期)2022-07-04 02:24:32
    機器人快速示教方法及示教點位姿變換的研究
    淺析沖壓件測量點的規(guī)劃
    模具制造(2019年10期)2020-01-06 09:13:08
    基于CAD模型的三坐標測量機測量點分布規(guī)劃
    PM2.5空中探測器的設(shè)計
    基于共面直線迭代加權(quán)最小二乘的相機位姿估計
    基于CAD模型的單目六自由度位姿測量
    垂直面內(nèi)建立基線的特殊點位高程測量法
    小型四旋翼飛行器位姿建模及其仿真
    中文字幕免费在线视频6| 国内精品宾馆在线| 中国国产av一级| 精品一区二区免费观看| 亚洲欧美日韩无卡精品| 国产爱豆传媒在线观看| av.在线天堂| 久久久久国产网址| 在线国产一区二区在线| 级片在线观看| 老司机午夜福利在线观看视频| 日韩 亚洲 欧美在线| 国内精品宾馆在线| 插阴视频在线观看视频| 91久久精品国产一区二区成人| 18禁在线播放成人免费| 亚洲性夜色夜夜综合| 亚洲av成人av| 99热只有精品国产| 搡女人真爽免费视频火全软件 | 波野结衣二区三区在线| 久久久久久久久久成人| 亚洲中文日韩欧美视频| 联通29元200g的流量卡| 亚洲人成网站高清观看| 男女那种视频在线观看| 男插女下体视频免费在线播放| 日本a在线网址| 免费看美女性在线毛片视频| 国产精品日韩av在线免费观看| 国产成人91sexporn| 久久精品夜夜夜夜夜久久蜜豆| 精品人妻视频免费看| 亚洲自拍偷在线| 久久久午夜欧美精品| 国产爱豆传媒在线观看| ponron亚洲| 国产爱豆传媒在线观看| 日韩一本色道免费dvd| 国产午夜福利久久久久久| 亚洲中文字幕日韩| 成年版毛片免费区| 亚洲国产精品合色在线| 伊人久久精品亚洲午夜| 久久鲁丝午夜福利片| 国产伦在线观看视频一区| 伊人久久精品亚洲午夜| 久久久久久久久久成人| 中文资源天堂在线| 久久久久免费精品人妻一区二区| 久久久欧美国产精品| 美女被艹到高潮喷水动态| 啦啦啦韩国在线观看视频| 亚洲在线观看片| 男女那种视频在线观看| 久久久久久久亚洲中文字幕| 国产综合懂色| 直男gayav资源| 国内久久婷婷六月综合欲色啪| 欧美高清性xxxxhd video| 国产精品久久电影中文字幕| 欧美另类亚洲清纯唯美| 人妻久久中文字幕网| 老司机福利观看| 蜜桃久久精品国产亚洲av| 午夜福利成人在线免费观看| 蜜桃久久精品国产亚洲av| 尤物成人国产欧美一区二区三区| 亚洲国产精品国产精品| 亚洲av成人精品一区久久| 又爽又黄a免费视频| 国产 一区 欧美 日韩| 成人美女网站在线观看视频| 国产激情偷乱视频一区二区| 久久久欧美国产精品| 成人永久免费在线观看视频| 一个人看视频在线观看www免费| 国产探花极品一区二区| 熟妇人妻久久中文字幕3abv| 欧美+日韩+精品| 搡老熟女国产l中国老女人| 久久久久久久久久成人| 黄色一级大片看看| 成人亚洲欧美一区二区av| 国产午夜福利久久久久久| 99久久久亚洲精品蜜臀av| 在线观看av片永久免费下载| av天堂在线播放| 精品久久久久久久久亚洲| 91av网一区二区| 日韩精品有码人妻一区| 成人国产麻豆网| 精品一区二区免费观看| 国产成年人精品一区二区| 精品久久久久久久久av| 九九热线精品视视频播放| 中文字幕熟女人妻在线| 最近视频中文字幕2019在线8| 美女xxoo啪啪120秒动态图| 在线观看美女被高潮喷水网站| 欧美日韩乱码在线| av在线播放精品| 久久午夜福利片| 国内精品宾馆在线| 狂野欧美白嫩少妇大欣赏| 国产69精品久久久久777片| 如何舔出高潮| 亚洲国产精品国产精品| 成人综合一区亚洲| 国产在线男女| 国产探花极品一区二区| 三级男女做爰猛烈吃奶摸视频| 永久网站在线| 午夜福利高清视频| 色在线成人网| 婷婷六月久久综合丁香| 色吧在线观看| 97在线视频观看| or卡值多少钱| 最近中文字幕高清免费大全6| 性色avwww在线观看| 欧美xxxx性猛交bbbb| 国产精品伦人一区二区| 精品久久久久久久久亚洲| 一区福利在线观看| 波多野结衣高清无吗| 亚洲av中文av极速乱| 91av网一区二区| 久久99热这里只有精品18| 国产私拍福利视频在线观看| 亚洲精品粉嫩美女一区| 国产一区二区三区在线臀色熟女| 欧美一区二区亚洲| 日韩av在线大香蕉| 国产亚洲av嫩草精品影院| 精品无人区乱码1区二区| 亚洲久久久久久中文字幕| 黄色欧美视频在线观看| 成人精品一区二区免费| 久久精品国产亚洲av涩爱 | 内射极品少妇av片p| 亚洲精品一卡2卡三卡4卡5卡| 国产极品精品免费视频能看的| 天天躁夜夜躁狠狠久久av| 精品久久久久久久久亚洲| 国产高清视频在线观看网站| 亚洲最大成人av| 天天躁日日操中文字幕| 久久国内精品自在自线图片| 乱系列少妇在线播放| 成人二区视频| 精品久久久久久久久久久久久| 大又大粗又爽又黄少妇毛片口| 97超碰精品成人国产| 精品福利观看| 不卡一级毛片| 亚洲av中文av极速乱| 国产高清激情床上av| 国产成人福利小说| 亚洲七黄色美女视频| 国产亚洲精品久久久com| 97超级碰碰碰精品色视频在线观看| 成人亚洲精品av一区二区| 成人美女网站在线观看视频| 熟妇人妻久久中文字幕3abv| 国产精品亚洲一级av第二区| 午夜福利在线在线| 久久久久久久久大av| 99在线人妻在线中文字幕| 欧美日韩乱码在线| 俺也久久电影网| 99热这里只有是精品在线观看| 午夜久久久久精精品| 国产亚洲精品综合一区在线观看| 波野结衣二区三区在线| 久久热精品热| 国产蜜桃级精品一区二区三区| 亚洲国产欧美人成| 99久久无色码亚洲精品果冻| 久久精品夜色国产| 亚洲国产精品久久男人天堂| 小蜜桃在线观看免费完整版高清| 亚洲专区国产一区二区| 日韩一本色道免费dvd| 久久久成人免费电影| or卡值多少钱| 精品欧美国产一区二区三| 欧美3d第一页| 观看免费一级毛片| 日韩,欧美,国产一区二区三区 | 97热精品久久久久久| 美女 人体艺术 gogo| 中文字幕人妻熟人妻熟丝袜美| 99久久精品一区二区三区| 内射极品少妇av片p| 1000部很黄的大片| 亚洲美女视频黄频| 国产黄a三级三级三级人| 日韩av在线大香蕉| 三级毛片av免费| 国产乱人偷精品视频| 在线播放无遮挡| 免费观看在线日韩| 午夜影院日韩av| 国产男人的电影天堂91| 日日摸夜夜添夜夜爱| 别揉我奶头~嗯~啊~动态视频| 晚上一个人看的免费电影| 国产伦精品一区二区三区四那| 欧美不卡视频在线免费观看| 中文字幕av成人在线电影| 国产私拍福利视频在线观看| 亚洲av一区综合| 国产真实乱freesex| 亚洲美女搞黄在线观看 | 日韩欧美三级三区| 日日干狠狠操夜夜爽| 中文字幕精品亚洲无线码一区| 欧洲精品卡2卡3卡4卡5卡区| 免费一级毛片在线播放高清视频| 亚洲国产高清在线一区二区三| 亚洲美女黄片视频| 亚洲国产欧美人成| 99riav亚洲国产免费| 极品教师在线视频| 国产欧美日韩精品一区二区| 亚洲专区国产一区二区| 欧美日韩在线观看h| 亚洲熟妇中文字幕五十中出| 免费人成视频x8x8入口观看| 精品无人区乱码1区二区| 久久久a久久爽久久v久久| 日韩国内少妇激情av| 国内少妇人妻偷人精品xxx网站| 少妇丰满av| 偷拍熟女少妇极品色| 最近在线观看免费完整版| 国产精品一区www在线观看| 中出人妻视频一区二区| 少妇的逼水好多| 午夜福利在线观看免费完整高清在 | 欧美性感艳星| 99久国产av精品国产电影| 久久精品国产鲁丝片午夜精品| 国产精品国产高清国产av| 日本与韩国留学比较| 国产精品一区二区三区四区久久| 一边摸一边抽搐一进一小说| 女人被狂操c到高潮| 波多野结衣巨乳人妻| 一区二区三区高清视频在线| 色综合亚洲欧美另类图片| 欧美日韩乱码在线| 99热全是精品| 看片在线看免费视频| 在线观看免费视频日本深夜| 男女下面进入的视频免费午夜| 国产精品女同一区二区软件| 最近2019中文字幕mv第一页| 国产精品一区二区免费欧美| 亚洲经典国产精华液单| 精品久久久久久久久亚洲| 国产精品爽爽va在线观看网站| 99久国产av精品国产电影| 国产aⅴ精品一区二区三区波| 一个人看视频在线观看www免费| 日日啪夜夜撸| 亚洲无线在线观看| 亚洲最大成人手机在线| 久久久久精品国产欧美久久久| 婷婷亚洲欧美| 91av网一区二区| 日本爱情动作片www.在线观看 | 国内少妇人妻偷人精品xxx网站| 国产私拍福利视频在线观看| 亚洲婷婷狠狠爱综合网| a级毛色黄片| 18+在线观看网站| 亚洲精品亚洲一区二区| 久久精品影院6| 男人的好看免费观看在线视频| 国内久久婷婷六月综合欲色啪| 少妇高潮的动态图| 一级a爱片免费观看的视频| 舔av片在线| 欧美xxxx性猛交bbbb| 免费观看精品视频网站| 我的老师免费观看完整版| 亚洲自偷自拍三级| 俺也久久电影网| 一本精品99久久精品77| 成年版毛片免费区| 日本 av在线| 日日摸夜夜添夜夜爱| 晚上一个人看的免费电影| 欧美人与善性xxx| 亚洲自拍偷在线| 免费看美女性在线毛片视频| 国产69精品久久久久777片| 中文资源天堂在线| 免费看a级黄色片| 国产精品精品国产色婷婷| 美女 人体艺术 gogo| 18+在线观看网站| 97碰自拍视频| 国国产精品蜜臀av免费| 别揉我奶头~嗯~啊~动态视频| 成人美女网站在线观看视频| 国内精品美女久久久久久| АⅤ资源中文在线天堂| 国产片特级美女逼逼视频| 国产色婷婷99| 国产高清视频在线播放一区| 亚洲精品色激情综合| 高清毛片免费观看视频网站| 国产综合懂色| 日本黄色视频三级网站网址| 日日撸夜夜添| 插阴视频在线观看视频| 久久久久久久亚洲中文字幕| 国产免费一级a男人的天堂| 99在线人妻在线中文字幕| 日韩亚洲欧美综合| 亚洲人成网站高清观看| 三级国产精品欧美在线观看| 麻豆国产97在线/欧美| 久99久视频精品免费| 久久久欧美国产精品| 嫩草影视91久久| 嫩草影院精品99| 日本黄色片子视频| 亚洲无线观看免费| 精品人妻视频免费看| 成人午夜高清在线视频| 高清毛片免费看| 欧美区成人在线视频| 国产精品人妻久久久影院| 麻豆国产av国片精品| av专区在线播放| 午夜福利在线在线| 国产精品女同一区二区软件| 欧美中文日本在线观看视频| 内射极品少妇av片p| 3wmmmm亚洲av在线观看| 亚洲精品影视一区二区三区av| av黄色大香蕉| 不卡一级毛片| 国产一区二区亚洲精品在线观看| 在线免费观看不下载黄p国产| 亚洲欧美成人精品一区二区| 亚洲精品日韩av片在线观看| 免费搜索国产男女视频| 精品人妻一区二区三区麻豆 | 亚洲国产精品sss在线观看| 亚洲aⅴ乱码一区二区在线播放| 99久久成人亚洲精品观看| 国产一级毛片七仙女欲春2| 国产一区二区激情短视频| 欧美色欧美亚洲另类二区| 国产探花极品一区二区| 亚洲欧美成人综合另类久久久 | 国产熟女欧美一区二区| 亚洲熟妇中文字幕五十中出| 在线免费观看不下载黄p国产| 久久久久免费精品人妻一区二区| 搡老岳熟女国产| 欧美xxxx黑人xx丫x性爽| 精品久久久久久久久久免费视频| 在线观看av片永久免费下载| 乱人视频在线观看| 亚洲欧美精品综合久久99| 波野结衣二区三区在线| 亚洲av美国av| 熟女人妻精品中文字幕| 日韩国内少妇激情av| 色吧在线观看| 国产精品一区www在线观看| 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 夜夜爽天天搞| 精品一区二区三区av网在线观看| 午夜福利高清视频| 国产黄色小视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 免费在线观看影片大全网站| 露出奶头的视频| 午夜福利成人在线免费观看| 高清午夜精品一区二区三区 | 99久久中文字幕三级久久日本| 别揉我奶头~嗯~啊~动态视频| 伦精品一区二区三区| 丰满的人妻完整版| 99在线视频只有这里精品首页| 欧美高清成人免费视频www| 成人三级黄色视频| 亚洲高清免费不卡视频| 亚洲人成网站高清观看| 久久久久久九九精品二区国产| 久久综合国产亚洲精品| 亚洲第一电影网av| 老熟妇仑乱视频hdxx| 嫩草影院新地址| 天堂影院成人在线观看| 老熟妇仑乱视频hdxx| 我要搜黄色片| 国产色爽女视频免费观看| 日韩精品青青久久久久久| 免费人成视频x8x8入口观看| 国产精品一区二区免费欧美| 此物有八面人人有两片| 精品免费久久久久久久清纯| 国产精品永久免费网站| 91久久精品电影网| 俄罗斯特黄特色一大片| 尤物成人国产欧美一区二区三区| 日韩人妻高清精品专区| 两个人视频免费观看高清| 天堂影院成人在线观看| 波多野结衣高清无吗| 一个人观看的视频www高清免费观看| 精品国内亚洲2022精品成人| 亚洲乱码一区二区免费版| 亚洲高清免费不卡视频| 亚洲美女搞黄在线观看 | av视频在线观看入口| 亚洲在线自拍视频| 午夜福利在线在线| 性插视频无遮挡在线免费观看| 一级毛片aaaaaa免费看小| 精品人妻视频免费看| 日本欧美国产在线视频| 亚洲中文日韩欧美视频| 日本 av在线| 欧美xxxx性猛交bbbb| 看十八女毛片水多多多| 免费高清视频大片| 欧美另类亚洲清纯唯美| 日本免费一区二区三区高清不卡| 最近中文字幕高清免费大全6| 精品久久久噜噜| 亚洲欧美日韩无卡精品| 99久久精品国产国产毛片| 欧美+日韩+精品| 十八禁网站免费在线| 欧美中文日本在线观看视频| 国产极品精品免费视频能看的| 最近中文字幕高清免费大全6| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕 | 综合色丁香网| 十八禁国产超污无遮挡网站| 午夜福利视频1000在线观看| 少妇熟女欧美另类| 99视频精品全部免费 在线| 乱人视频在线观看| 精品一区二区三区av网在线观看| 午夜影院日韩av| 久久国内精品自在自线图片| 久久鲁丝午夜福利片| 99在线人妻在线中文字幕| 精品99又大又爽又粗少妇毛片| 午夜福利高清视频| av卡一久久| 搡老熟女国产l中国老女人| 久久久久久伊人网av| 成人二区视频| av在线蜜桃| 久久韩国三级中文字幕| 赤兔流量卡办理| 我的女老师完整版在线观看| 性欧美人与动物交配| 寂寞人妻少妇视频99o| 免费搜索国产男女视频| 在线国产一区二区在线| 91麻豆精品激情在线观看国产| 国产单亲对白刺激| 国产免费一级a男人的天堂| 91久久精品国产一区二区成人| 看黄色毛片网站| 丰满人妻一区二区三区视频av| 亚洲精品456在线播放app| 欧美一区二区精品小视频在线| 欧美成人一区二区免费高清观看| 久久亚洲国产成人精品v| 日韩三级伦理在线观看| 特大巨黑吊av在线直播| 亚洲精品日韩在线中文字幕 | videossex国产| 中文资源天堂在线| 长腿黑丝高跟| 久久久久精品国产欧美久久久| 国产一区二区三区av在线 | www日本黄色视频网| 男人和女人高潮做爰伦理| 成年版毛片免费区| 欧美性猛交╳xxx乱大交人| 大又大粗又爽又黄少妇毛片口| 国产高清激情床上av| 成年免费大片在线观看| 麻豆久久精品国产亚洲av| 18禁在线无遮挡免费观看视频 | 亚洲精品粉嫩美女一区| 国产在视频线在精品| 精品久久久久久久久亚洲| 一级黄片播放器| 人妻夜夜爽99麻豆av| 中国美女看黄片| 亚洲国产精品合色在线| 九九热线精品视视频播放| 大又大粗又爽又黄少妇毛片口| 国产精品av视频在线免费观看| 亚洲精品456在线播放app| 国产高潮美女av| 日本欧美国产在线视频| 看十八女毛片水多多多| 免费大片18禁| 免费av观看视频| 国产精品三级大全| 精品久久久噜噜| 婷婷六月久久综合丁香| 校园春色视频在线观看| 国产精品无大码| 一进一出抽搐动态| 日韩欧美 国产精品| 久久九九热精品免费| 变态另类成人亚洲欧美熟女| 91在线精品国自产拍蜜月| 久久亚洲国产成人精品v| 国产精品三级大全| 午夜精品一区二区三区免费看| 日韩欧美一区二区三区在线观看| 中文字幕免费在线视频6| 亚洲自拍偷在线| 三级毛片av免费| 日本免费一区二区三区高清不卡| 精品一区二区三区视频在线| 日本黄色视频三级网站网址| 一区二区三区免费毛片| 免费看光身美女| 天美传媒精品一区二区| 中文资源天堂在线| 亚洲精品国产成人久久av| 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| av在线老鸭窝| 久久久久久久久久黄片| 又爽又黄a免费视频| 国产不卡一卡二| 校园春色视频在线观看| 五月伊人婷婷丁香| 老司机午夜福利在线观看视频| 久久人人精品亚洲av| 秋霞在线观看毛片| 国产三级中文精品| 简卡轻食公司| 国产高清视频在线播放一区| 午夜福利视频1000在线观看| 亚洲国产高清在线一区二区三| 最后的刺客免费高清国语| 国产老妇女一区| 亚洲精品久久国产高清桃花| 老熟妇乱子伦视频在线观看| 91久久精品国产一区二区成人| 欧美激情久久久久久爽电影| 狠狠狠狠99中文字幕| 不卡视频在线观看欧美| 国产黄色小视频在线观看| 亚洲精华国产精华液的使用体验 | 色哟哟·www| 美女被艹到高潮喷水动态| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 精品久久久久久久久av| 99久久精品一区二区三区| 淫秽高清视频在线观看| 亚洲成av人片在线播放无| 嫩草影院精品99| 色5月婷婷丁香| 18禁在线播放成人免费| av在线亚洲专区| 成人漫画全彩无遮挡| 日日摸夜夜添夜夜添av毛片| 最近视频中文字幕2019在线8| 久久久久免费精品人妻一区二区| 美女xxoo啪啪120秒动态图| 国产伦在线观看视频一区| 午夜亚洲福利在线播放| 伊人久久精品亚洲午夜| 日本a在线网址| 一个人看视频在线观看www免费| 男女之事视频高清在线观看| 国产伦一二天堂av在线观看| 国产黄色小视频在线观看| 国产成人freesex在线 | 亚洲欧美成人精品一区二区| 我要看日韩黄色一级片| 国产伦在线观看视频一区| 少妇裸体淫交视频免费看高清| 色5月婷婷丁香| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 卡戴珊不雅视频在线播放| 成年版毛片免费区| 亚洲精品影视一区二区三区av| 在线观看66精品国产| 日韩欧美精品免费久久| 一级黄色大片毛片| 久久九九热精品免费| 国产高清视频在线播放一区| 非洲黑人性xxxx精品又粗又长| 亚洲美女黄片视频| 成年av动漫网址| 日日撸夜夜添| 亚洲欧美日韩高清在线视频| 久久精品影院6|