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

    基于Kriging模型插值的孔位修正策略

    2020-12-01 12:22:36石循磊張繼文劉順濤陳懇
    航空學(xué)報 2020年9期
    關(guān)鍵詞:制孔孔位插值

    石循磊,張繼文,劉順濤,陳懇,*

    1. 清華大學(xué) 機械工程系,北京 100084 2. 航空工業(yè)成都飛機工業(yè)(集團)有限責(zé)任公司,成都 610091

    飛機普遍采用鉚接或螺栓連接的形式實現(xiàn)零部件的裝配,疊層制孔和連接任務(wù)占據(jù)了大部分裝配任務(wù)[1]?,F(xiàn)代飛機一般具有復(fù)雜的氣動外形,且零件剛度較低,裝配時容易產(chǎn)生較大的變形和位置偏移[2-3],因此難以像普通加工機床一樣,通過離線編程和在線對刀的形式一次性完成制孔作業(yè)任務(wù)。為了適應(yīng)零部件的變形和位置移動,保證制孔孔位的邊距準(zhǔn)確性,在機器制孔之前通過人工方式或特殊方式制備基準(zhǔn)孔(或稱預(yù)連接孔),再通過激光跟蹤儀或者其他視覺測量系統(tǒng)檢測這些基準(zhǔn)孔的實際位置,將其與理論位置相比較,計算出基準(zhǔn)孔的孔位偏差[4-5]。依據(jù)已制的基準(zhǔn)孔位偏差,通過插值算法推算出待制孔的位置偏差,再對自動制孔程序進行校正,調(diào)整實際制孔位置,實現(xiàn)孔位補償?shù)墓δ堋?/p>

    在自動孔位尋找校正過程中,通過已制基準(zhǔn)孔位推斷待制孔孔位偏差的補償算法顯得尤為重要。然而,飛機零件的弱剛性以及復(fù)雜外形特性為孔位推算帶來了多個難點:其一,裝配過程中,工件的變形并非簡單的剛性位置移動或旋轉(zhuǎn),通常還伴隨有較大幅度的彈性扭轉(zhuǎn)和彎曲[6],因此無法基于一系列基準(zhǔn)孔采用裝配體的剛體變換形式獲得待制孔孔位;其二,基準(zhǔn)孔通常位于不同的零件裝配件上,例如飛機的隔框縱梁長桁結(jié)構(gòu)交替出現(xiàn),基準(zhǔn)孔也隨之交替出現(xiàn),若基準(zhǔn)孔位與零件綁定,則為鑒別和記錄帶來繁雜的標(biāo)記過程,因而極易出現(xiàn)錯誤,進而導(dǎo)致整個孔位校正失效。

    目前大多數(shù)工程化解決方案采用線性或三次樣條插值的方法校正兩個或多個基準(zhǔn)孔之間的孔位偏差。這是一種近似的解決方案,飛機的復(fù)雜曲面及變形特征很難服從這一假設(shè),特別在一些曲率變化較大的位置,發(fā)生孔位校正錯誤的概率大大增加[7]。為了保證校正的局部特性,不得不增加基準(zhǔn)孔的數(shù)量,并嚴(yán)重依賴工藝實驗和人工經(jīng)驗選擇基準(zhǔn)孔位[8]。同時線性插值、三次樣條插值需要選定插值節(jié)點,然而,對于三維空間位置而言,很難區(qū)分樣值點的內(nèi)外關(guān)系,也難以確定插值節(jié)點,而且在外插值情況下容易出現(xiàn)較大誤差[9]。

    在理論研究方面,Zhu等[10]通過基準(zhǔn)孔構(gòu)建了雙線性誤差平面,實現(xiàn)了對制孔區(qū)域的孔位進行線性插值補償。周叔陽等[11]將模擬退火算法應(yīng)用于預(yù)連接孔的排布問題,并抽象為數(shù)學(xué)模型,給出了預(yù)連接孔排布優(yōu)化算法及流程。董輝躍等[12]從飛機曲面構(gòu)造原理角度,提出了一種基于誤差Coons曲面函數(shù)的補償方法。王龍飛等[13]通過視覺輔助和手眼標(biāo)定,計算出基準(zhǔn)孔對應(yīng)的機器人驅(qū)動坐標(biāo)誤差變換矩陣,進而補償待制孔誤差。王青等[14]提出一種基于孔邊距偏差修正約束孔位和利用Shepard插值修正其余孔位的策略,實現(xiàn)了較好的孔位修正效果,但是缺少對孔位預(yù)測誤差的分析以及基準(zhǔn)孔的補制策略。本文提出一種采用Kriging模型插值算法的孔位偏差校正算法。該方法利用Kriging插值模型,依據(jù)基準(zhǔn)孔位偏差推算待制孔偏差和預(yù)測標(biāo)準(zhǔn)差。算法依賴孔位估計的均方誤差評估基準(zhǔn)孔的數(shù)量和分布質(zhì)量,進而可以在基準(zhǔn)孔實際選擇和布置時,依賴反復(fù)的迭代和計算,指導(dǎo)增添新的基準(zhǔn)孔,從而提升孔位偏差估計的精度。

    1 孔位補償原理

    根據(jù)基準(zhǔn)孔位偏差推算待制孔位置偏差的原理如圖1所示。

    圖1 由基準(zhǔn)孔位偏差推算待制孔位置偏差的原理Fig.1 Principle of estimating position deviation of hole to be made from reference hole

    已經(jīng)依據(jù)理論位置及工件裝配實際狀況制備出N個基準(zhǔn)孔?;鶞?zhǔn)孔實際位置與CAD模型中的理論位置存在一定偏差。

    (1)

    設(shè)CAD模型中的一個待制孔理論位置坐標(biāo)為(u*,v*,w*)T,則補償目標(biāo)是推算出輸出實際制孔位置與待制孔理論位置的偏差值:

    (2)

    上述過程本質(zhì)是分別建立了孔位偏差的3個空間坐標(biāo)分量Δu、Δv、Δw與孔位理論位置[u,v,w]T之間的響應(yīng)面模型。建立這種響應(yīng)面關(guān)系的一種方法是通過材料力學(xué)、彈性力學(xué)等復(fù)雜的理論推導(dǎo)出這一關(guān)系,但難度過大,而本文方法通過多次實驗的點以及一種特殊的函數(shù)關(guān)系模擬出這種關(guān)系,建立Kriging模型。

    2 Kriging模型建立

    Kriging模型屬于響應(yīng)面模型,其本質(zhì)是一種數(shù)據(jù)擬合的過程,將函數(shù)視作一種隨機過程,通過有限的數(shù)據(jù)擬合出這一隨機過程的概率分布規(guī)律,從而得到函數(shù),計算出在空間坐標(biāo)變化時函數(shù)數(shù)值的變化規(guī)律,進而計算出未知點的數(shù)值。Kriging模型具有良好的非線性函數(shù)擬合特性,不僅提供了未知點的預(yù)測值,還提供了預(yù)測標(biāo)準(zhǔn)差[15]。

    Kriging模型的基本表達式為

    y(x)=fT(x)β+μ(x)

    (3)

    式中:f(x)為q維基函數(shù)向量;β為q維系數(shù)向量;μ(x)為隨機變量,每給出一個采樣點(或稱自變量)x,在f(x)、β和μ(x)已知的條件下,即可計算出其響應(yīng)函數(shù)值y(x)[16]。β和μ(x)未知時,當(dāng)具有一系列采樣點x(1),x(2),…,x(N)及對應(yīng)的y(1),y(2),…,y(N),即可估計出其參數(shù),進而求得其他位置x*的函數(shù)值y(x*)。

    如果將μ(x)當(dāng)做滿足正態(tài)分布、與位置無關(guān)的隨機變量μ,則可以通過多次實驗消除。然而,在Kriging模型內(nèi),μ(x)是一個依賴于x的隨機變量,且兩個樣值點x(i)、x(j)距離越近,其相關(guān)性越強。針對全局范圍內(nèi)的采樣點x,滿足數(shù)學(xué)期望E[μ(x)]=0,同時,對于任意兩個采樣點x(i)、x(j),滿足協(xié)方差關(guān)系:

    (4)

    式中:σ2為采樣點對應(yīng)的隨機變量的方差;R(θ,p,x(i),x(j))為關(guān)聯(lián)函數(shù),體現(xiàn)了x(i)、x(j)兩個采樣點間的相關(guān)系數(shù):

    (5)

    式中:θh為x向量中h分量的相關(guān)強度;ph為x向量中h分量的相關(guān)形式[17]。與擬合回歸不同,Kriging模型主要通過μ(x)描述y(x)的響應(yīng)值,并通常將fT(x)設(shè)置為簡單的1T常數(shù)。該模型假設(shè)μ(x)和y(x)僅與狀態(tài)x有關(guān),因而不存在多項式擬合中的擬合誤差。

    在fT(x)=1T和p預(yù)先指定的條件下,當(dāng)給出多組采樣點時,通過最大似然估計法,估算參數(shù){β,σ2,θ},從而得到模型估計結(jié)果:

    (6)

    式(6)顯示Kriging模型估計的結(jié)果是一系列關(guān)聯(lián)函數(shù)的線性組合所表達的連續(xù)可微函數(shù)式:

    (7)

    式中:ai為最大似然假設(shè)所估計的N個常數(shù)[19]。

    依據(jù)R(θ,x*,x(i))的定義(設(shè)ph=2時):

    (8)

    3 孔位偏差估計及改進策略

    針對孔位偏差問題,距離待制孔距離越大的基準(zhǔn)孔位偏差對當(dāng)前待制孔偏差的影響越小,而距離待制孔距離越小的基準(zhǔn)孔孔位偏差則對當(dāng)前待制孔的影響較大,這種特性能夠有效描述待制孔與基準(zhǔn)孔之間的偏差關(guān)聯(lián)關(guān)系。如圖2所示。

    根據(jù)Kriging模型的定義及參數(shù)估計方法,當(dāng)x*=x(i)時,即取已知實驗點時,由響應(yīng)面函數(shù)計算可得:

    (9)

    式(9)表明由式(6)計算出的基準(zhǔn)孔x(i)處的偏差值仍然是實測值y(i),且其估計偏差為0。因此,Kriging模型所建立的是對已知實驗點的一種特殊的插值,而相比于線性插值或樣條插值,孔位偏差不再限于內(nèi)插值,而是可以實現(xiàn)一定程度的外插值。同時對于一個待制孔的誤差,其不僅取決于兩端的基準(zhǔn)孔,還取決于與其臨近的若干個孔的偏差,這與一般裝配的誤差特性一致。

    圖2 待制孔與基準(zhǔn)孔之間的偏差關(guān)聯(lián)關(guān)系Fig.2 Deviation correlation between hole to be made and reference hole

    (10)

    (11)

    因此,采用Kriging模型插值推算孔位偏差,能夠依賴孔位估計的均方誤差評估基準(zhǔn)孔的數(shù)量和分布質(zhì)量,進而可以在基準(zhǔn)孔實際選擇和布置時,依賴反復(fù)的迭代和計算,指導(dǎo)增添新的基準(zhǔn)孔,從而提升孔位偏差估計的精度。此外,可以通過式(7)采用一種連續(xù)可微的函數(shù)形式模擬孔位偏差,從而繪制出孔位誤差分布圖,有利于裝配質(zhì)量以及基準(zhǔn)孔制備的迭代改進。算法的總體結(jié)構(gòu)如圖3所示。

    圖3 算法總體結(jié)構(gòu)Fig.3 Overall structure of algorithm

    4 實驗與分析

    飛機孔位偏差來源可分為待制孔零件的變形和平移旋轉(zhuǎn)兩種,實驗中分別驗證兩種偏差來源形式下本文算法的有效性。

    針對零件變形帶來的孔位偏差,采用有限元仿真實驗來驗證算法的有效性。利用有限元軟件Abaqus建模,如圖4所示,將蒙皮簡化為邊長約4.8 m、厚度3 mm的梯形曲面工件,連接于骨架上,并施加載荷。坐標(biāo)系y軸定義為與曲面蒙皮3個相距較遠頂點組成的平面垂直的方向,即y軸近似與蒙皮表面垂直,而x和z軸在該平面內(nèi)。

    圖4 有限元仿真載荷示意圖Fig.4 Load diagram of finite element simulation

    有限元仿真結(jié)果如圖5所示,圖中標(biāo)出了25個基準(zhǔn)孔,在蒙皮上另取50個點作為待制孔的位置。將有限元仿真結(jié)果中基準(zhǔn)孔的位置作為實際孔位,可得基準(zhǔn)孔的位置偏差(位移)Δu、Δv、Δw作為Kriging模型的因變量,以基準(zhǔn)孔理論位置為自變量,可建立3個Kriging模型。利用建立好的模型估計待制孔的應(yīng)制孔位置,結(jié)果如圖6所示。

    圖6中,eu、ev、ew是Kriging模型估計的孔位偏差值(位移)的誤差,定義為

    (12)

    從圖6中可以看出,變形導(dǎo)致的x和z方向位移量較小,孔位估計誤差均小于0.1 mm。y方向位移達到10 mm量級,可認(rèn)為是實際裝配時壁板的極限變形,因此主要分析y方向?qū)嶒灲Y(jié)果。

    針對零件整體平移及旋轉(zhuǎn)帶來的孔位偏差,偏移前后孔位坐標(biāo)有Q=R0Q0+T0

    (13)

    式中:Q0為孔位理論坐標(biāo);Q為偏移后坐標(biāo);R0為旋轉(zhuǎn)矩陣;T0為平移變換矩陣,因此孔位偏差量為Q-Q0。

    在實際裝配過程中,蒙皮壁板放置在飛機骨架上,受骨架約束,可認(rèn)為壁板僅繞圖8中y軸旋轉(zhuǎn)。驗證實驗中取旋轉(zhuǎn)角度為1.15°,可得旋轉(zhuǎn)矩陣:

    圖5 有限元仿真位移結(jié)果云圖Fig.5 Displacement nephogram of finite element simulation

    (14)

    同樣采用圖5中的25個基準(zhǔn)孔建立Kriging模型,50個待制孔的x方向和z方向最大位移均達10 mm量級,而孔位估計誤差如圖9所示,均小于0.05 mm。因此對于零件整體旋轉(zhuǎn)平移導(dǎo)致的孔位偏差,該算法能很好地估計偏移后的位置。

    利用上述方法在某型飛機數(shù)字化裝配現(xiàn)場實驗的結(jié)果見表1,實物實驗中基準(zhǔn)孔的偏移量均小于1 mm,該方法估計出的待制孔的偏差量與測量值相差很小。

    考慮該算法的計算效率,實驗使用的計算機CPU為i7-8700k,利用R語言編程,由25個基準(zhǔn)孔建立Kriging模型并估計50個待制孔位所需運算時間為0.214 5 s,50個基準(zhǔn)孔估計100個待制孔耗時0.306 2 s,90個基準(zhǔn)孔估計150個待制孔耗時0.480 7 s,而實際生產(chǎn)中孔位估計并無實時性要求,所以計算效率完全能達到要求。與線性插值等傳統(tǒng)方法相比較,該算法除了能給出待制孔的預(yù)測標(biāo)準(zhǔn)誤差外,另一個顯著優(yōu)勢體現(xiàn)在曲面區(qū)域外插值的情況。例如在上述有限元仿真實驗中,一個待制孔的坐標(biāo)為(-1 592.7, 75.251, -1 008.3),該待制孔附近的兩個基準(zhǔn)孔坐標(biāo)為(-1 692.1, 66.179, -1 334.5)和(-1 690.7, 66.312, -814.05),以這兩個點為插值節(jié)點,線性插值得待制孔的y方向偏差值為154.78 mm,而仿真結(jié)果為2.990 3 mm,本文算法估計結(jié)果為2.847 1 mm,此時線性插值出現(xiàn)很大誤差,而該算法仍有較好效果。

    圖6 變形導(dǎo)致孔位偏差的實驗結(jié)果Fig.6 Experimental results of hole displacement caused by deformation

    圖7 增添基準(zhǔn)孔后實驗結(jié)果Fig.7 Experimental results after adding reference holes

    圖8 實驗件坐標(biāo)系定義Fig.8 Definition of experimental coordinate system

    多組實驗表明本文方法能較準(zhǔn)確、高效地通過基準(zhǔn)孔信息估計待制孔的位移,并且給出有實際參考意義的待制孔位置的預(yù)測標(biāo)準(zhǔn)誤差,用以指導(dǎo)排布或增添基準(zhǔn)孔,最后通過在合適位置增添基準(zhǔn)孔能有效提高待制孔位移估計的準(zhǔn)確度。

    圖9 平移旋轉(zhuǎn)導(dǎo)致孔位偏差的實驗結(jié)果Fig.9 Experimental results of hole displacement caused by translation and rotation

    表1 實物實驗待制孔位估計結(jié)果

    5 結(jié) 論

    1) 提出了一種基于基準(zhǔn)孔的孔位校正方法,該方法利用Kriging模型,依據(jù)基準(zhǔn)孔位偏差推算出待制孔偏差和預(yù)測標(biāo)準(zhǔn)誤差。

    2)相比于線性或三次樣條插值,該方法適用范圍更大,且能給出孔位估計的標(biāo)準(zhǔn)誤差,以評估基準(zhǔn)孔的數(shù)量和分布質(zhì)量,進而指導(dǎo)增添新的基準(zhǔn)孔。

    3)利用有限元仿真實驗,驗證了在零件變形及平移旋轉(zhuǎn)導(dǎo)致孔位偏差情況下本文方法的有效性,在合理增添基準(zhǔn)孔后,孔位估計的誤差小于0.3 mm,滿足航空裝配要求。

    猜你喜歡
    制孔孔位插值
    增強現(xiàn)實眼鏡輔助的線纜連接器裝配狀態(tài)智能檢錯方法
    柔性導(dǎo)軌制孔系統(tǒng)研究
    難加工材料(CFRP/Ti)疊層自適應(yīng)制孔研究
    機器人制孔系統(tǒng)與制孔工藝參數(shù)優(yōu)化方法研究
    主貯油罐地腳螺栓有限元受力分析
    山西建筑(2019年12期)2019-08-05 09:24:50
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    基于Twincat的移動機器人制孔系統(tǒng)
    一種改進FFT多譜線插值諧波分析方法
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    高厚徑比鉆孔孔位精度提升研究
    亚洲av免费高清在线观看| 成人综合一区亚洲| 成人亚洲精品一区在线观看| 久久99蜜桃精品久久| 少妇的丰满在线观看| 丝袜美足系列| 韩国精品一区二区三区 | av在线播放精品| 国产 精品1| 国产在线一区二区三区精| 久久精品国产亚洲av涩爱| 亚洲国产精品一区二区三区在线| 最近最新中文字幕免费大全7| 亚洲 欧美一区二区三区| 亚洲精品第二区| 日本91视频免费播放| 国产xxxxx性猛交| 免费看不卡的av| 黄色配什么色好看| 欧美日韩综合久久久久久| 精品熟女少妇av免费看| 91午夜精品亚洲一区二区三区| 亚洲精品日韩在线中文字幕| 男女啪啪激烈高潮av片| 亚洲精品乱码久久久久久按摩| 这个男人来自地球电影免费观看 | 日韩制服丝袜自拍偷拍| videos熟女内射| 久久狼人影院| 国产成人一区二区在线| 性高湖久久久久久久久免费观看| 五月开心婷婷网| 国产一区二区激情短视频 | 在线亚洲精品国产二区图片欧美| 侵犯人妻中文字幕一二三四区| 亚洲精品成人av观看孕妇| 韩国av在线不卡| 国产一级毛片在线| 成人二区视频| 如何舔出高潮| 日日啪夜夜爽| 超碰97精品在线观看| 18禁动态无遮挡网站| av黄色大香蕉| 精品国产一区二区三区四区第35| 国产探花极品一区二区| 看十八女毛片水多多多| 大话2 男鬼变身卡| 99久久中文字幕三级久久日本| av卡一久久| 亚洲精品日本国产第一区| 中文字幕亚洲精品专区| 亚洲少妇的诱惑av| 秋霞伦理黄片| a级片在线免费高清观看视频| 欧美+日韩+精品| 国产成人精品福利久久| av.在线天堂| 亚洲一区二区三区欧美精品| 26uuu在线亚洲综合色| av片东京热男人的天堂| 国产精品三级大全| 久久青草综合色| 好男人视频免费观看在线| 51国产日韩欧美| 久久婷婷青草| 免费人妻精品一区二区三区视频| 精品一区在线观看国产| 免费观看性生交大片5| 9191精品国产免费久久| av不卡在线播放| 少妇人妻精品综合一区二区| 欧美 日韩 精品 国产| 久久精品aⅴ一区二区三区四区 | 婷婷成人精品国产| 日本91视频免费播放| 热99久久久久精品小说推荐| 亚洲成色77777| 亚洲成人av在线免费| 国产精品 国内视频| 日本av免费视频播放| 国产综合精华液| 啦啦啦啦在线视频资源| 日韩熟女老妇一区二区性免费视频| 欧美日韩精品成人综合77777| 亚洲精品一二三| 色5月婷婷丁香| 国产女主播在线喷水免费视频网站| 如日韩欧美国产精品一区二区三区| 亚洲综合色网址| 国产精品一区www在线观看| 午夜免费观看性视频| 亚洲精品自拍成人| 黄片播放在线免费| 韩国精品一区二区三区 | 最近2019中文字幕mv第一页| 午夜福利,免费看| 日本与韩国留学比较| 最近中文字幕高清免费大全6| 黄色 视频免费看| 欧美+日韩+精品| 久久精品国产自在天天线| 久久 成人 亚洲| 中文字幕亚洲精品专区| 天堂中文最新版在线下载| 一边亲一边摸免费视频| 免费黄色在线免费观看| 精品人妻熟女毛片av久久网站| 五月伊人婷婷丁香| 高清欧美精品videossex| 日本欧美视频一区| 视频在线观看一区二区三区| 黑人欧美特级aaaaaa片| 久久精品国产亚洲av天美| 久久久国产精品麻豆| 人人妻人人澡人人看| 亚洲高清免费不卡视频| 亚洲精品久久成人aⅴ小说| 1024视频免费在线观看| 综合色丁香网| 久久久久久久国产电影| 亚洲av日韩在线播放| 激情五月婷婷亚洲| 天堂俺去俺来也www色官网| 国产毛片在线视频| 美女国产视频在线观看| www日本在线高清视频| 秋霞伦理黄片| 亚洲欧美成人综合另类久久久| 亚洲欧美日韩卡通动漫| 精品国产一区二区三区四区第35| 国产精品成人在线| 91在线精品国自产拍蜜月| 精品99又大又爽又粗少妇毛片| 国产精品一区二区在线观看99| 成人国产av品久久久| 青春草亚洲视频在线观看| 韩国av在线不卡| 久久鲁丝午夜福利片| 菩萨蛮人人尽说江南好唐韦庄| 热99国产精品久久久久久7| 亚洲精品久久成人aⅴ小说| 欧美日韩成人在线一区二区| 99精国产麻豆久久婷婷| www.色视频.com| 曰老女人黄片| 女性生殖器流出的白浆| 国产国拍精品亚洲av在线观看| 欧美xxxx性猛交bbbb| 精品一区二区三卡| 国产黄频视频在线观看| videossex国产| 亚洲丝袜综合中文字幕| 亚洲熟女精品中文字幕| 欧美人与性动交α欧美软件 | 黄色毛片三级朝国网站| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片 | 国产一区二区三区综合在线观看 | 久久国产精品男人的天堂亚洲 | 天天躁夜夜躁狠狠躁躁| 精品熟女少妇av免费看| 免费观看性生交大片5| videos熟女内射| 成人免费观看视频高清| 少妇人妻 视频| 欧美日韩国产mv在线观看视频| 日韩一区二区视频免费看| 精品熟女少妇av免费看| 免费不卡的大黄色大毛片视频在线观看| 在现免费观看毛片| h视频一区二区三区| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 啦啦啦啦在线视频资源| 久久99一区二区三区| 制服人妻中文乱码| 黑人巨大精品欧美一区二区蜜桃 | 999精品在线视频| 黄色配什么色好看| 黑人猛操日本美女一级片| 晚上一个人看的免费电影| videosex国产| 国产日韩一区二区三区精品不卡| 久久精品国产亚洲av天美| 纵有疾风起免费观看全集完整版| 91精品伊人久久大香线蕉| 国产乱人偷精品视频| 国产av码专区亚洲av| 啦啦啦在线观看免费高清www| 侵犯人妻中文字幕一二三四区| 少妇高潮的动态图| 欧美精品av麻豆av| 高清毛片免费看| 国产精品久久久久久久电影| 综合色丁香网| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕 | 欧美激情极品国产一区二区三区 | 少妇被粗大的猛进出69影院 | 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜制服| 中文字幕av电影在线播放| 国产xxxxx性猛交| 午夜av观看不卡| 一区二区三区乱码不卡18| 成人午夜精彩视频在线观看| 色吧在线观看| 国产69精品久久久久777片| 国产男女超爽视频在线观看| 日韩成人av中文字幕在线观看| av线在线观看网站| 欧美精品一区二区免费开放| 中文字幕精品免费在线观看视频 | 国精品久久久久久国模美| 亚洲国产精品国产精品| 婷婷色麻豆天堂久久| 卡戴珊不雅视频在线播放| 在线观看免费视频网站a站| 看非洲黑人一级黄片| 亚洲国产精品成人久久小说| 亚洲精品美女久久av网站| 国产国拍精品亚洲av在线观看| 999精品在线视频| 国产免费福利视频在线观看| 国产精品久久久久久av不卡| 亚洲性久久影院| 99热网站在线观看| 少妇熟女欧美另类| 欧美xxⅹ黑人| 国产免费一区二区三区四区乱码| 18+在线观看网站| 极品人妻少妇av视频| 欧美成人午夜精品| 欧美日韩一区二区视频在线观看视频在线| 国产激情久久老熟女| 欧美日韩亚洲高清精品| 国产永久视频网站| 18禁在线无遮挡免费观看视频| 久久这里有精品视频免费| 精品国产一区二区久久| 亚洲五月色婷婷综合| 九九在线视频观看精品| 亚洲欧洲精品一区二区精品久久久 | 欧美xxⅹ黑人| 99九九在线精品视频| 免费看不卡的av| 18禁裸乳无遮挡动漫免费视频| 欧美亚洲日本最大视频资源| 精品人妻偷拍中文字幕| 婷婷色综合www| 亚洲精品美女久久av网站| 色哟哟·www| 内地一区二区视频在线| 亚洲天堂av无毛| 人妻少妇偷人精品九色| 久久久久国产网址| 欧美少妇被猛烈插入视频| 亚洲国产成人一精品久久久| 大香蕉久久成人网| 亚洲欧美清纯卡通| 久久久久国产网址| 观看美女的网站| 亚洲丝袜综合中文字幕| av又黄又爽大尺度在线免费看| 久久99热6这里只有精品| 一级毛片黄色毛片免费观看视频| 少妇猛男粗大的猛烈进出视频| h视频一区二区三区| av片东京热男人的天堂| 内地一区二区视频在线| 汤姆久久久久久久影院中文字幕| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 久久精品久久精品一区二区三区| 午夜视频国产福利| 亚洲av福利一区| 亚洲图色成人| 国产精品偷伦视频观看了| 国产老妇伦熟女老妇高清| 女人久久www免费人成看片| 激情五月婷婷亚洲| 久久这里有精品视频免费| 美女国产视频在线观看| 在线亚洲精品国产二区图片欧美| 亚洲美女黄色视频免费看| av女优亚洲男人天堂| 欧美日韩综合久久久久久| 最近中文字幕高清免费大全6| 少妇被粗大的猛进出69影院 | 韩国精品一区二区三区 | 一个人免费看片子| 精品国产乱码久久久久久小说| 一边摸一边做爽爽视频免费| 成人毛片60女人毛片免费| 精品久久国产蜜桃| 亚洲国产成人一精品久久久| 69精品国产乱码久久久| 在线观看www视频免费| 中文天堂在线官网| 亚洲国产精品成人久久小说| 国产福利在线免费观看视频| 大香蕉97超碰在线| 制服丝袜香蕉在线| 又黄又粗又硬又大视频| 最近中文字幕2019免费版| 观看美女的网站| 免费看不卡的av| 国产国拍精品亚洲av在线观看| 97在线人人人人妻| 日韩制服骚丝袜av| 免费黄网站久久成人精品| 男人操女人黄网站| 黄色 视频免费看| 蜜臀久久99精品久久宅男| 亚洲欧美中文字幕日韩二区| 大话2 男鬼变身卡| 女性被躁到高潮视频| 亚洲av电影在线进入| 国精品久久久久久国模美| 人妻一区二区av| 多毛熟女@视频| 日韩视频在线欧美| 一本—道久久a久久精品蜜桃钙片| 中文字幕亚洲精品专区| 日韩伦理黄色片| 国产精品久久久久久久久免| 99香蕉大伊视频| 国产日韩欧美视频二区| 一区二区三区精品91| 免费黄色在线免费观看| 2022亚洲国产成人精品| 久久99热6这里只有精品| 亚洲精品乱码久久久久久按摩| 下体分泌物呈黄色| 永久免费av网站大全| 一区二区三区乱码不卡18| 人人妻人人澡人人爽人人夜夜| 人人妻人人添人人爽欧美一区卜| 91精品三级在线观看| 欧美性感艳星| 亚洲国产精品999| 免费久久久久久久精品成人欧美视频 | www.av在线官网国产| xxxhd国产人妻xxx| 欧美激情国产日韩精品一区| 美女xxoo啪啪120秒动态图| 80岁老熟妇乱子伦牲交| 久久ye,这里只有精品| 久久婷婷青草| 欧美少妇被猛烈插入视频| 午夜免费观看性视频| 一边亲一边摸免费视频| 两个人免费观看高清视频| 午夜久久久在线观看| 国产成人精品一,二区| 亚洲国产av影院在线观看| 久久人人97超碰香蕉20202| 男人操女人黄网站| 久久久久久久国产电影| 国产一区二区在线观看av| 黑人巨大精品欧美一区二区蜜桃 | 国产日韩欧美在线精品| 国产国拍精品亚洲av在线观看| 国产色婷婷99| xxx大片免费视频| 国产色爽女视频免费观看| 一本色道久久久久久精品综合| 精品人妻一区二区三区麻豆| 国产精品人妻久久久影院| videos熟女内射| 亚洲国产精品专区欧美| 国产乱人偷精品视频| 日韩制服骚丝袜av| av在线老鸭窝| 亚洲内射少妇av| 日韩一区二区视频免费看| 男女边吃奶边做爰视频| 日韩成人伦理影院| 天堂8中文在线网| 欧美老熟妇乱子伦牲交| 久久人人爽av亚洲精品天堂| 99国产精品免费福利视频| 亚洲美女视频黄频| 国产精品久久久久久av不卡| 不卡视频在线观看欧美| 校园人妻丝袜中文字幕| 在线天堂中文资源库| 久久国产亚洲av麻豆专区| 黄色一级大片看看| 在线观看免费高清a一片| 国产一区二区激情短视频 | 欧美丝袜亚洲另类| 26uuu在线亚洲综合色| 18禁观看日本| 看非洲黑人一级黄片| 女的被弄到高潮叫床怎么办| 欧美日韩成人在线一区二区| 超碰97精品在线观看| 哪个播放器可以免费观看大片| 欧美成人精品欧美一级黄| 26uuu在线亚洲综合色| 高清毛片免费看| 亚洲av男天堂| 91成人精品电影| xxx大片免费视频| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级| 日韩一区二区视频免费看| 激情五月婷婷亚洲| 蜜桃在线观看..| 青春草亚洲视频在线观看| 在线精品无人区一区二区三| 精品少妇内射三级| 国产成人精品久久久久久| 亚洲国产精品成人久久小说| 蜜臀久久99精品久久宅男| 亚洲精品中文字幕在线视频| 国产一区二区在线观看日韩| 国产乱人偷精品视频| 亚洲精品乱久久久久久| 考比视频在线观看| 校园人妻丝袜中文字幕| 免费av不卡在线播放| 夜夜骑夜夜射夜夜干| 男女免费视频国产| 51国产日韩欧美| 中文欧美无线码| 天天躁夜夜躁狠狠躁躁| 亚洲精品色激情综合| 精品国产一区二区三区四区第35| 国产精品一区二区在线不卡| 九草在线视频观看| 亚洲国产最新在线播放| 满18在线观看网站| 侵犯人妻中文字幕一二三四区| 街头女战士在线观看网站| 婷婷成人精品国产| 国产成人91sexporn| 一个人免费看片子| 精品一区在线观看国产| 久久久久国产网址| 又黄又粗又硬又大视频| 亚洲精品自拍成人| 在线亚洲精品国产二区图片欧美| 一本色道久久久久久精品综合| 免费人成在线观看视频色| 中文字幕最新亚洲高清| av网站免费在线观看视频| 春色校园在线视频观看| 成人毛片a级毛片在线播放| 免费黄色在线免费观看| 少妇被粗大猛烈的视频| videos熟女内射| 一边亲一边摸免费视频| 久久久久久久久久成人| 美女内射精品一级片tv| 寂寞人妻少妇视频99o| 亚洲精品第二区| 欧美人与性动交α欧美精品济南到 | 亚洲欧美一区二区三区国产| 国产精品成人在线| 国产乱人偷精品视频| 一边亲一边摸免费视频| 亚洲在久久综合| 美女福利国产在线| 亚洲欧美精品自产自拍| 久久99热6这里只有精品| 欧美亚洲 丝袜 人妻 在线| 亚洲国产精品成人久久小说| 亚洲精品自拍成人| 国产爽快片一区二区三区| 免费观看a级毛片全部| 黄网站色视频无遮挡免费观看| av在线老鸭窝| 99热国产这里只有精品6| 国产精品蜜桃在线观看| 最近最新中文字幕免费大全7| 日韩欧美一区视频在线观看| 观看美女的网站| 亚洲人成77777在线视频| 精品一区二区三区视频在线| 国产精品偷伦视频观看了| av线在线观看网站| 欧美国产精品一级二级三级| 国产在线免费精品| 中文乱码字字幕精品一区二区三区| 欧美日韩成人在线一区二区| 亚洲av在线观看美女高潮| 国产亚洲欧美精品永久| 亚洲性久久影院| 美女内射精品一级片tv| 成人毛片a级毛片在线播放| 亚洲精品中文字幕在线视频| 18禁裸乳无遮挡动漫免费视频| 午夜免费男女啪啪视频观看| 热re99久久国产66热| 免费观看在线日韩| 狂野欧美激情性bbbbbb| 免费播放大片免费观看视频在线观看| 成人综合一区亚洲| 最近中文字幕高清免费大全6| 一级毛片黄色毛片免费观看视频| 秋霞在线观看毛片| 亚洲国产av新网站| 久久毛片免费看一区二区三区| 一二三四中文在线观看免费高清| 最近最新中文字幕大全免费视频 | 国产麻豆69| 日韩精品免费视频一区二区三区 | 80岁老熟妇乱子伦牲交| 国产免费视频播放在线视频| 免费女性裸体啪啪无遮挡网站| 97精品久久久久久久久久精品| 最黄视频免费看| 少妇精品久久久久久久| 国产精品一区二区在线不卡| 成人亚洲欧美一区二区av| 一二三四在线观看免费中文在 | 免费看av在线观看网站| 各种免费的搞黄视频| 蜜桃在线观看..| 王馨瑶露胸无遮挡在线观看| 最近最新中文字幕大全免费视频 | 日韩三级伦理在线观看| 亚洲成av片中文字幕在线观看 | 51国产日韩欧美| 日韩不卡一区二区三区视频在线| 精品人妻在线不人妻| 色视频在线一区二区三区| 精品国产一区二区三区久久久樱花| 熟女人妻精品中文字幕| 免费人妻精品一区二区三区视频| 91成人精品电影| 欧美性感艳星| 天堂中文最新版在线下载| 热re99久久精品国产66热6| 少妇人妻 视频| 国产午夜精品一二区理论片| av国产久精品久网站免费入址| 欧美精品一区二区大全| 国产爽快片一区二区三区| 男女边吃奶边做爰视频| 亚洲av电影在线进入| 久久人人爽人人片av| av片东京热男人的天堂| 色5月婷婷丁香| 欧美人与性动交α欧美软件 | 久久久a久久爽久久v久久| 搡老乐熟女国产| 黑丝袜美女国产一区| 国产乱来视频区| 少妇人妻久久综合中文| 免费看不卡的av| 一本久久精品| 成人黄色视频免费在线看| 男人添女人高潮全过程视频| 卡戴珊不雅视频在线播放| 日韩一本色道免费dvd| 久久精品国产鲁丝片午夜精品| 97在线视频观看| 成人亚洲精品一区在线观看| 久久人人爽人人片av| 十分钟在线观看高清视频www| 高清视频免费观看一区二区| 99久久精品国产国产毛片| 成人黄色视频免费在线看| 国产又爽黄色视频| 国产成人免费无遮挡视频| 大话2 男鬼变身卡| 国产爽快片一区二区三区| 国产亚洲最大av| 大陆偷拍与自拍| 婷婷色av中文字幕| 九色成人免费人妻av| 亚洲经典国产精华液单| 国产乱来视频区| 国产一区亚洲一区在线观看| 亚洲精品美女久久久久99蜜臀 | 只有这里有精品99| 亚洲第一av免费看| 久久精品久久精品一区二区三区| 国产永久视频网站| 涩涩av久久男人的天堂| 欧美性感艳星| 男女下面插进去视频免费观看 | 国产精品三级大全| 亚洲国产精品999| 在现免费观看毛片| 又大又黄又爽视频免费| 亚洲精品一二三| 老司机亚洲免费影院| 国产国语露脸激情在线看| 成人午夜精彩视频在线观看| av电影中文网址| 最黄视频免费看| 国产不卡av网站在线观看| 一本久久精品| 中文字幕人妻熟女乱码| 久久久欧美国产精品| 亚洲综合精品二区| 中文字幕人妻熟女乱码| 国产又色又爽无遮挡免| 最近手机中文字幕大全| 免费在线观看黄色视频的| 亚洲精品一二三| 女性被躁到高潮视频| 精品一区二区三区四区五区乱码 | 国产精品人妻久久久影院| 欧美变态另类bdsm刘玥| 黄片播放在线免费| 你懂的网址亚洲精品在线观看| 亚洲精品456在线播放app| 男女国产视频网站| 日本免费在线观看一区|