• <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| 国产成人精品婷婷| 欧美日韩一区二区视频在线观看视频在线| 亚洲av成人精品一区久久| 国产精品久久久久久精品古装| 国产淫片久久久久久久久| 亚洲av福利一区| 22中文网久久字幕| 人妻一区二区av| 亚洲国产精品国产精品| 国产av国产精品国产| 免费久久久久久久精品成人欧美视频 | 伦理电影免费视频| 成人黄色视频免费在线看| 你懂的网址亚洲精品在线观看| 国产高清三级在线| 少妇熟女欧美另类| videossex国产| 免费不卡的大黄色大毛片视频在线观看| 一级,二级,三级黄色视频| 国产高清不卡午夜福利| 久久久a久久爽久久v久久| 欧美精品亚洲一区二区| 狂野欧美激情性xxxx在线观看| 日韩中字成人| 欧美激情国产日韩精品一区| 亚洲精品中文字幕在线视频 | 天美传媒精品一区二区| 成人国产av品久久久| 久久精品国产鲁丝片午夜精品| 亚洲色图综合在线观看| 一区二区三区精品91| 青春草视频在线免费观看| 亚洲人与动物交配视频| 亚洲激情五月婷婷啪啪| 秋霞伦理黄片| 少妇熟女欧美另类| 亚洲情色 制服丝袜| 久久精品国产自在天天线| 人妻制服诱惑在线中文字幕| 乱人伦中国视频| 欧美日韩国产mv在线观看视频| 日本-黄色视频高清免费观看| 久久久久久人妻| 18禁在线无遮挡免费观看视频| 少妇人妻一区二区三区视频| 黄色配什么色好看| 免费看日本二区| 男女国产视频网站| 免费看日本二区| 一级,二级,三级黄色视频| 精品少妇黑人巨大在线播放| 免费看日本二区| 成人无遮挡网站| 日日啪夜夜撸| 特大巨黑吊av在线直播| 下体分泌物呈黄色| 好男人视频免费观看在线| 男女边吃奶边做爰视频| 亚洲av欧美aⅴ国产| 五月开心婷婷网| 中文字幕久久专区| 韩国高清视频一区二区三区| 精品少妇内射三级| 免费久久久久久久精品成人欧美视频 | 久久久久网色| 卡戴珊不雅视频在线播放| 一级片'在线观看视频| 亚洲av成人精品一二三区| 纵有疾风起免费观看全集完整版| www.av在线官网国产| 日本黄色片子视频| 一本大道久久a久久精品| 亚洲av电影在线观看一区二区三区| 男女边摸边吃奶| 少妇被粗大的猛进出69影院 | 亚洲av电影在线观看一区二区三区| 大码成人一级视频| 国产极品粉嫩免费观看在线 | 免费黄色在线免费观看| 97超视频在线观看视频| 黄片无遮挡物在线观看| 国产精品一区二区在线观看99| 高清在线视频一区二区三区| 极品少妇高潮喷水抽搐| 久久久久国产精品人妻一区二区| 久久97久久精品| xxx大片免费视频| 精品国产乱码久久久久久小说| 亚洲精品视频女| 国产熟女欧美一区二区| 亚洲,欧美,日韩| 国产精品麻豆人妻色哟哟久久| 国产片特级美女逼逼视频| 插阴视频在线观看视频| 久久午夜福利片| 啦啦啦在线观看免费高清www| 哪个播放器可以免费观看大片| 亚洲国产精品一区二区三区在线| .国产精品久久| 久久亚洲国产成人精品v| 高清毛片免费看| 秋霞在线观看毛片| 久久久久网色| 国产黄片视频在线免费观看| 亚洲av福利一区| 国产极品天堂在线| 乱人伦中国视频| 青春草国产在线视频| 天堂8中文在线网| 成人国产麻豆网| 97超视频在线观看视频| 精华霜和精华液先用哪个| 一边亲一边摸免费视频| 久久午夜综合久久蜜桃| 亚洲国产精品一区三区| 国产伦精品一区二区三区视频9| 亚洲精品亚洲一区二区| 成人亚洲精品一区在线观看| 亚洲不卡免费看| 午夜福利视频精品| av又黄又爽大尺度在线免费看| 国产亚洲5aaaaa淫片| 亚洲中文av在线| 精品视频人人做人人爽| 欧美 日韩 精品 国产| 国产精品国产三级国产av玫瑰| 免费大片黄手机在线观看| 综合色丁香网| 99热网站在线观看| 久热这里只有精品99| 国产日韩欧美视频二区| 中国国产av一级| 久久99蜜桃精品久久| 最近中文字幕2019免费版| 日本欧美视频一区| 国产69精品久久久久777片| 午夜福利视频精品| 夫妻性生交免费视频一级片| 日韩欧美一区视频在线观看 | 91成人精品电影| 国产一区亚洲一区在线观看| 午夜福利在线观看免费完整高清在| 在线看a的网站| 秋霞伦理黄片| 久久久久久久久久久丰满| 国产免费一区二区三区四区乱码| 人人澡人人妻人| 日日撸夜夜添| 国产成人免费观看mmmm| 久久精品国产自在天天线| 91aial.com中文字幕在线观看| 欧美精品一区二区大全| 熟女电影av网| 伦理电影大哥的女人| 中文字幕久久专区| 亚洲av不卡在线观看| 91精品国产九色| 97精品久久久久久久久久精品| 国产黄色视频一区二区在线观看| 黄色视频在线播放观看不卡| 麻豆乱淫一区二区| 在线观看美女被高潮喷水网站| 日韩强制内射视频| 国产一区亚洲一区在线观看| 又粗又硬又长又爽又黄的视频| 亚洲精品乱久久久久久| 亚洲精品乱久久久久久| 日韩三级伦理在线观看| 亚洲美女视频黄频| 2018国产大陆天天弄谢| 国产伦精品一区二区三区视频9| 日本午夜av视频| 国产欧美日韩综合在线一区二区 | 老司机影院毛片| 国产伦在线观看视频一区| 精品少妇黑人巨大在线播放| 亚洲国产精品999| 久久精品国产亚洲av涩爱| 午夜视频国产福利| 久久婷婷青草| 涩涩av久久男人的天堂| 亚洲精品中文字幕在线视频 | 国产黄片视频在线免费观看| 少妇猛男粗大的猛烈进出视频| 黄色毛片三级朝国网站 | 精品一区二区三卡| 精品一品国产午夜福利视频| 黄色欧美视频在线观看| 一个人免费看片子| 精品亚洲乱码少妇综合久久| 特大巨黑吊av在线直播| 国产乱来视频区| 少妇熟女欧美另类| 人妻人人澡人人爽人人| 国语对白做爰xxxⅹ性视频网站| 少妇猛男粗大的猛烈进出视频| 国产av国产精品国产| 9色porny在线观看| 我的女老师完整版在线观看| 久久精品久久久久久噜噜老黄| 欧美日韩精品成人综合77777| 久久99蜜桃精品久久| 久久热精品热| 边亲边吃奶的免费视频| 永久网站在线| 99视频精品全部免费 在线| 伊人久久精品亚洲午夜| 简卡轻食公司| 欧美日韩国产mv在线观看视频| 国产av国产精品国产| 久久99热这里只频精品6学生| 亚洲精品国产色婷婷电影| 2018国产大陆天天弄谢| 亚洲自偷自拍三级| 久久久国产一区二区| 观看av在线不卡| 少妇丰满av| 最新中文字幕久久久久| 免费播放大片免费观看视频在线观看| 中文在线观看免费www的网站| 黑人巨大精品欧美一区二区蜜桃 | 亚洲欧美一区二区三区黑人 | 肉色欧美久久久久久久蜜桃| 亚洲精品日韩av片在线观看| 亚洲内射少妇av| 亚洲,一卡二卡三卡| 精品一区二区三区视频在线| 亚洲伊人久久精品综合| 久久综合国产亚洲精品| 一本久久精品| 亚洲精品国产av蜜桃| 日韩在线高清观看一区二区三区| .国产精品久久| 黄片无遮挡物在线观看| 大码成人一级视频| 久久精品久久精品一区二区三区| 国产成人a∨麻豆精品| 国产黄色视频一区二区在线观看| 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看 | 国产国拍精品亚洲av在线观看| 少妇丰满av| 一本色道久久久久久精品综合| 成人无遮挡网站| 欧美高清成人免费视频www| 久久鲁丝午夜福利片| 美女视频免费永久观看网站| 亚洲美女视频黄频| 黑人猛操日本美女一级片| 春色校园在线视频观看| 国产乱人偷精品视频| 夫妻性生交免费视频一级片| 亚洲国产av新网站| 成人漫画全彩无遮挡| 成人漫画全彩无遮挡| 国产老妇伦熟女老妇高清| 性色av一级| 国产精品一二三区在线看| 狂野欧美激情性bbbbbb| 婷婷色av中文字幕| 午夜精品国产一区二区电影| 亚洲精品456在线播放app| 国产精品伦人一区二区| 51国产日韩欧美| 成年人午夜在线观看视频| 少妇熟女欧美另类| 久久久久国产精品人妻一区二区| 国语对白做爰xxxⅹ性视频网站| 免费大片黄手机在线观看| 国产欧美亚洲国产| 久久鲁丝午夜福利片| 伊人亚洲综合成人网| 在线观看免费日韩欧美大片 | 国产精品偷伦视频观看了| 你懂的网址亚洲精品在线观看| 夫妻性生交免费视频一级片| 欧美日韩综合久久久久久| 99re6热这里在线精品视频| 亚洲性久久影院| 午夜免费男女啪啪视频观看| 人人妻人人看人人澡| 一本色道久久久久久精品综合| 亚洲一级一片aⅴ在线观看| 亚洲精品色激情综合| 国产精品人妻久久久影院| av免费在线看不卡| 亚洲精品久久久久久婷婷小说| 免费高清在线观看视频在线观看| 亚洲精品,欧美精品| 国产一区二区在线观看日韩| 3wmmmm亚洲av在线观看| 如日韩欧美国产精品一区二区三区 | 欧美最新免费一区二区三区| 亚洲性久久影院| 国产高清不卡午夜福利| 久久免费观看电影| av福利片在线| 亚洲av.av天堂| 国国产精品蜜臀av免费| 少妇猛男粗大的猛烈进出视频| 热99国产精品久久久久久7| 国产在线一区二区三区精| 18+在线观看网站| 岛国毛片在线播放| 黑丝袜美女国产一区| 男女边摸边吃奶| 成人毛片60女人毛片免费| 精品久久久久久久久亚洲| 夫妻性生交免费视频一级片| 肉色欧美久久久久久久蜜桃| 人人妻人人澡人人看| 亚洲伊人久久精品综合| av国产精品久久久久影院| 如日韩欧美国产精品一区二区三区 | 高清视频免费观看一区二区| 黄色配什么色好看| 最近中文字幕2019免费版| 亚洲欧美精品专区久久| 一区二区av电影网| 自拍偷自拍亚洲精品老妇| 18禁在线播放成人免费| 六月丁香七月| 久久国内精品自在自线图片| 国产黄色视频一区二区在线观看| 香蕉精品网在线| 日韩视频在线欧美| 亚洲久久久国产精品| 欧美日韩一区二区视频在线观看视频在线| 纵有疾风起免费观看全集完整版| 午夜日本视频在线| 久久精品久久久久久噜噜老黄| 又粗又硬又长又爽又黄的视频| 高清欧美精品videossex| 国产无遮挡羞羞视频在线观看| 精品人妻一区二区三区麻豆| 在线精品无人区一区二区三| 肉色欧美久久久久久久蜜桃| 人妻制服诱惑在线中文字幕| 97超视频在线观看视频| 99九九线精品视频在线观看视频| 亚洲,一卡二卡三卡| 国产黄色视频一区二区在线观看| 又爽又黄a免费视频| 欧美变态另类bdsm刘玥| 免费av不卡在线播放| 伊人亚洲综合成人网| 婷婷色综合大香蕉| 久久久久视频综合| 美女xxoo啪啪120秒动态图| 一本大道久久a久久精品| 一级av片app| 国产一区有黄有色的免费视频| 国产毛片在线视频| 欧美日韩综合久久久久久| 亚洲av国产av综合av卡| 亚洲av男天堂| 亚洲av电影在线观看一区二区三区| 毛片一级片免费看久久久久| 亚洲av综合色区一区| 久久这里有精品视频免费| 美女脱内裤让男人舔精品视频| 这个男人来自地球电影免费观看 | 亚洲av中文av极速乱| 久久精品国产自在天天线| 免费观看性生交大片5| 美女cb高潮喷水在线观看| 亚洲精品第二区| 99久国产av精品国产电影| 国语对白做爰xxxⅹ性视频网站| 国产精品99久久99久久久不卡 | 老女人水多毛片| 97超碰精品成人国产| 中国美白少妇内射xxxbb| 黄色日韩在线| 蜜臀久久99精品久久宅男| 亚洲精品乱久久久久久| 久久精品国产a三级三级三级| 午夜激情福利司机影院| 女的被弄到高潮叫床怎么办| 久久精品熟女亚洲av麻豆精品| 久久久国产精品麻豆| av网站免费在线观看视频| 欧美bdsm另类| 免费av中文字幕在线| 能在线免费看毛片的网站| 欧美3d第一页| 欧美97在线视频| 成人免费观看视频高清| 婷婷色麻豆天堂久久| 少妇丰满av| 少妇的逼好多水| 在线观看免费高清a一片| av女优亚洲男人天堂| 亚洲av在线观看美女高潮| 亚洲精品一二三| av.在线天堂| 国产精品无大码| 一边亲一边摸免费视频| 毛片一级片免费看久久久久| 韩国高清视频一区二区三区| 久久久久久久亚洲中文字幕| 十八禁网站网址无遮挡 | 汤姆久久久久久久影院中文字幕| 精品国产一区二区久久| 精品一区在线观看国产| 日韩免费高清中文字幕av| 丰满饥渴人妻一区二区三| 国产亚洲av片在线观看秒播厂| 亚洲欧美日韩东京热| 欧美精品一区二区大全| 日韩制服骚丝袜av| 激情五月婷婷亚洲| 美女福利国产在线| 日韩一区二区三区影片| 日本91视频免费播放| 亚洲精品一二三| 不卡视频在线观看欧美| 纵有疾风起免费观看全集完整版| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡动漫免费视频| 一区二区三区四区激情视频| 一级黄片播放器| 国产日韩欧美亚洲二区| 国产精品99久久久久久久久| a级毛片在线看网站| 日本wwww免费看| 久久综合国产亚洲精品| 精品一区二区三区视频在线| 老熟女久久久| videos熟女内射| 国产欧美亚洲国产| 欧美日韩综合久久久久久| 亚洲一区二区三区欧美精品| 国产精品国产三级专区第一集| 亚洲av男天堂| 校园人妻丝袜中文字幕| 婷婷色综合www| 亚洲国产成人一精品久久久| 国产在线一区二区三区精| 日韩成人伦理影院| 22中文网久久字幕| 亚洲精品乱久久久久久| 久久韩国三级中文字幕| 如何舔出高潮| 国产伦理片在线播放av一区| av一本久久久久| 亚洲欧美日韩卡通动漫| 成年人午夜在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 久久久久国产网址| 亚洲精品,欧美精品| 亚洲精品国产成人久久av| 亚洲国产色片| 一级爰片在线观看| 不卡视频在线观看欧美| 国产黄色免费在线视频| 国产精品免费大片| 久久久久精品性色| 日本欧美视频一区| 亚洲精品久久久久久婷婷小说| 国产精品免费大片| 18禁在线无遮挡免费观看视频| 我要看黄色一级片免费的| 涩涩av久久男人的天堂| .国产精品久久| 国产伦精品一区二区三区视频9| 亚洲色图综合在线观看| 国产一级毛片在线| 五月天丁香电影| 国产探花极品一区二区| 成人国产av品久久久| 精品国产乱码久久久久久小说| 在线观看国产h片| 日韩精品有码人妻一区| 激情五月婷婷亚洲| 97在线人人人人妻| 夫妻性生交免费视频一级片| 欧美另类一区| 亚洲精品成人av观看孕妇| 亚洲精华国产精华液的使用体验| 在线观看av片永久免费下载| 亚洲内射少妇av| 麻豆成人午夜福利视频| 欧美日韩国产mv在线观看视频| 亚洲va在线va天堂va国产| 极品少妇高潮喷水抽搐| 一级毛片黄色毛片免费观看视频| 亚洲综合色惰| 国产精品国产三级国产专区5o| 蜜桃在线观看..| 婷婷色麻豆天堂久久| 狂野欧美激情性xxxx在线观看| 国产中年淑女户外野战色| 国产成人一区二区在线| 国产 一区精品| 久久97久久精品| 熟妇人妻不卡中文字幕| 亚洲成人av在线免费| 国产一区亚洲一区在线观看| 色视频www国产| 欧美人与善性xxx| 美女主播在线视频| 日日爽夜夜爽网站| 国产亚洲一区二区精品| 久久99一区二区三区| 男女国产视频网站| 韩国高清视频一区二区三区| 一个人看视频在线观看www免费| 欧美bdsm另类| 精品国产国语对白av| 午夜影院在线不卡| 22中文网久久字幕| 国精品久久久久久国模美| 极品少妇高潮喷水抽搐| 高清午夜精品一区二区三区| 丰满人妻一区二区三区视频av| 国产精品久久久久成人av| 久久国产精品男人的天堂亚洲 | 九色成人免费人妻av| 亚洲高清免费不卡视频| 新久久久久国产一级毛片| 大陆偷拍与自拍| 啦啦啦在线观看免费高清www| 99久国产av精品国产电影| 欧美 日韩 精品 国产| 在线观看一区二区三区激情| 亚洲欧美精品自产自拍| 少妇猛男粗大的猛烈进出视频| 国产又色又爽无遮挡免| 大码成人一级视频| 老熟女久久久| 观看av在线不卡| 三级国产精品片| 中文天堂在线官网| 成人无遮挡网站| 精品亚洲成a人片在线观看| 最近中文字幕2019免费版| 色婷婷av一区二区三区视频| 久久韩国三级中文字幕| 日日啪夜夜爽| 全区人妻精品视频| 一级毛片 在线播放| av在线播放精品| 极品少妇高潮喷水抽搐| 欧美bdsm另类| 午夜免费鲁丝| 亚洲无线观看免费| 成人特级av手机在线观看| 欧美97在线视频| 十分钟在线观看高清视频www | 少妇熟女欧美另类| 欧美精品一区二区大全| 夜夜骑夜夜射夜夜干| .国产精品久久| 亚洲欧美日韩卡通动漫| 人人妻人人爽人人添夜夜欢视频 | 大片免费播放器 马上看| 人妻制服诱惑在线中文字幕| 9色porny在线观看| 欧美一级a爱片免费观看看| 日本色播在线视频| 在线观看www视频免费| 日韩av不卡免费在线播放| 久久精品久久精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 中文天堂在线官网| 久久久久久久精品精品| 欧美成人午夜免费资源| 午夜福利网站1000一区二区三区| 狂野欧美激情性xxxx在线观看| 三上悠亚av全集在线观看 | 亚洲真实伦在线观看| 在线观看美女被高潮喷水网站| 熟妇人妻不卡中文字幕| 国产白丝娇喘喷水9色精品| 边亲边吃奶的免费视频| 亚洲av电影在线观看一区二区三区| 中文乱码字字幕精品一区二区三区| 国产精品不卡视频一区二区| 交换朋友夫妻互换小说| 综合色丁香网| 丝袜喷水一区| 日韩av在线免费看完整版不卡| 亚洲欧美清纯卡通| 亚洲精品一二三| 人人妻人人澡人人看| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 观看av在线不卡| 女性被躁到高潮视频| 国产又色又爽无遮挡免| 26uuu在线亚洲综合色| 欧美 亚洲 国产 日韩一| 日韩在线高清观看一区二区三区| 亚洲美女黄色视频免费看| 国产高清有码在线观看视频| 久久久a久久爽久久v久久| 伊人亚洲综合成人网| 最近中文字幕高清免费大全6| 乱人伦中国视频| av福利片在线观看| 亚洲欧洲日产国产| 色视频www国产| 欧美区成人在线视频| 亚洲内射少妇av| 极品教师在线视频| 欧美日韩一区二区视频在线观看视频在线| 一级a做视频免费观看|