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

    基于三維點云螺紋中軸線的擬合方法研究

    2020-09-08 08:48:50謝張寧姜志華朱惠臣李智瑋雷李華王道檔
    計量學(xué)報 2020年8期
    關(guān)鍵詞:單位向量外螺紋中軸線

    謝張寧, 姜志華, 朱惠臣, 李智瑋, 雷李華,,孔 明, 王道檔, 張 波

    (1. 中國計量大學(xué)計量測試工程學(xué)院,浙江杭州310018;2. 上海計量測試技術(shù)研究院,上海201203)

    1 引 言

    高精度螺紋規(guī)是精密機械、航天航空、石油化工、核電等領(lǐng)域中極為重要的基礎(chǔ)零部件,其幾何精度對裝配質(zhì)量具有決定性影響[1]?,F(xiàn)有的高精度測量方法眾多,但均是基于一維、二維的螺紋測量技術(shù),無法獲取三維綜合參數(shù),檢測結(jié)果難以反映實際精度,導(dǎo)致裝配質(zhì)量不能準(zhǔn)確控制。要改變“粗糙螺紋”生產(chǎn)的現(xiàn)狀,需要對三維螺紋測量儀器設(shè)備進(jìn)行研究開發(fā)。

    當(dāng)前的三維測量方法在實際測量前,首先要實現(xiàn)機器坐標(biāo)系MCS與工件坐標(biāo)系WCS轉(zhuǎn)換[2],在螺紋面的測量中,螺紋中軸線的準(zhǔn)確測定對螺紋的后續(xù)測量至關(guān)重要。然而,在對螺紋規(guī)測量時,仍舊是依靠測量人員的經(jīng)驗,手動調(diào)整螺紋規(guī)位置,然后進(jìn)行后續(xù)測量。這種測量方法不僅增加了測量難度,也增大了測量結(jié)果的不確定度。

    本文研究了一種能夠快速精準(zhǔn)擬合螺紋中軸線的方法。在測量時,螺紋規(guī)放置在工作臺上進(jìn)行調(diào)整,但是工件坐標(biāo)系不可能與機器坐標(biāo)系完全重合,所以要使用螺紋規(guī)的基準(zhǔn)來建立工件坐標(biāo)系。將測量儀器對準(zhǔn)螺紋規(guī)上下底面的頂針孔實現(xiàn)定位,在機器坐標(biāo)系MCS下,非接觸掃描得到的螺紋規(guī)外表面點云;經(jīng)過隨機增量算法計算點云凸包,制成螺紋凸包點集;基于最小二乘的概念建立螺紋中軸線的擬合數(shù)學(xué)模型,得出中軸線的單位向量v(a,b,c);再基于單位向量得到螺紋底面中心A,坐標(biāo)旋轉(zhuǎn)矩陣R,坐標(biāo)偏移矩陣S;螺紋測量儀器根據(jù)上述參數(shù)得出工件坐標(biāo)系。

    由于螺紋復(fù)雜的三維形體結(jié)構(gòu),在建立數(shù)學(xué)模型對螺紋面進(jìn)行直接計算時,牙槽表面、螺紋表面缺陷點云數(shù)據(jù)均會影響螺紋中軸線的擬合,導(dǎo)致結(jié)果不準(zhǔn)確。為防止干擾,本文用隨機增量算法計算螺紋表面三維點云凸包(convex hull)[3];將螺紋牙槽表面、螺紋表面缺陷點云數(shù)據(jù)篩去,對于空間離散點集,求解最小的凸多面體,其中把空間中點集所有的點包含在里面的凸多面體就是凸包[4]?;谌S點云凸包建立中軸線擬合數(shù)學(xué)模型,能有效濾除掉螺紋牙槽、硬件缺陷、掃描誤差等形成的干擾點云,減少了后續(xù)算法的運算量,加快了運行時間,提升了運算結(jié)果。

    經(jīng)過計算凸包的螺紋表面點云實質(zhì)上是以公稱直徑為上下底直徑的圓柱面。要擬合該圓柱面的軸線就要以面上的點與軸線的位置關(guān)系建立數(shù)學(xué)模型,并采用優(yōu)化算法求解。常用的優(yōu)化算法有最小二乘法、爬山算法、模擬退火算法、遺傳算法等[5~8]。用爬山算法、模擬退火算法、遺傳算法等隨機優(yōu)化算法確定唯一解的方法雖然可行,但隨機優(yōu)化算法大多處理高階次線性問題;目前通用、認(rèn)可度高的算法仍為最小二乘法,在工程測量上,用最小二乘法解決超定方程[9]、線性擬合[10~12]等數(shù)據(jù)量大的低階次線性問題速度快、結(jié)果準(zhǔn)確;且隨機算法優(yōu)化過程具有隨機性,容易出現(xiàn)早熟或過擬合等問題。因此,擬合螺紋中軸線采用最小二次算法具有較強的實用性。

    2 螺紋中軸線數(shù)學(xué)模型

    在空間機器坐標(biāo)系中,螺紋中軸線l可以表示為經(jīng)過點A(x0,y0,z0),單位向量為v(a,b,c)的直線l,其點斜式為[13]:

    根據(jù)最小二乘原理,使得經(jīng)過三維凸包計算的螺紋表面點云P(xi,yi,zi),i=1,…,N到擬合直線距離d的平方和最?。?/p>

    (1)

    (2)

    假設(shè):

    (3)

    將式(2)、式(3)代入式(1),可得最小二乘意義上的螺紋中軸線數(shù)學(xué)模型:

    (bδz-cδy)2]+f(a,b,c)

    (4)

    由式(4)可以看出,對于任意的a,b,c均有:

    (5)

    當(dāng)δx=δy=δz=0時,有:

    min{f(a,b,c)}
    a2+b2+c2=1

    其中:

    f(a,b,c)=(1-a2)B11+(1-b2)B22+

    (1-c2)B33-2abB12-2acB13-2bcB23

    (6)

    3 基于三維點云的最小二乘擬合

    3.1 制成三維凸包點集

    由于螺紋表面三維結(jié)構(gòu)復(fù)雜以及可能存在缺陷的原因,要對點云數(shù)據(jù)計算凸包。計算凸包是根據(jù)算法在點云中篩選出能構(gòu)造包含所有點的最小多面體的點集。

    螺紋表面點云是三維數(shù)據(jù)點,且數(shù)據(jù)量龐大,用隨機增量算法計算凸包效果好、速度最快[15]。計算流程為:

    (1)首先用篩選函數(shù)找到x,y,z軸上的邊界點,如x,y,z值最小的點,設(shè)其為P1。

    (2)隨機找到與P1不共線的2個點P2,P3,再找到不共面的第4個點P4,構(gòu)成四面體。

    (3)若新點Pi在凸包內(nèi)部或邊界,則忽略。若新點Pi在凸包外部,見圖1(a),就要計算原凸包相對于Pi的面,見圖1(b)所示;將Pi點能看到的所有面刪除,同時將邊上的點與Pi相連組成新的凸包,見圖1(c)。

    (4)依次遍歷余下的點,制成三維凸包點集。

    圖1 隨機增量算法計算凸包Fig.1 Random incremental algorithm to calculate convex hull

    3.2 單位向量的求解

    將式(6)引入矩陣B=[Bij],整理成矩陣式,則可得:

    由于B是對稱方陣,不難看出,當(dāng)[a,b,c]T是矩陣B的最大特征值λmax(B)對應(yīng)的單位特征向量時,上式可以取得最小[16]。此時,[a,b,c]T為螺紋中軸線單位方向向量:

    3.3 底面中心點的求解

    在轉(zhuǎn)換后的坐標(biāo)系o-x′y′z′中,取圓環(huán)最外層,設(shè)其擬合曲線為:

    式中M′為螺紋外徑擬合結(jié)果。將上式展開,令:

    整理得:

    最小二乘法求解上述超定方程解得[18]:

    其中:

    3.4 測量坐標(biāo)系到工件坐標(biāo)系的轉(zhuǎn)換

    測量坐標(biāo)系到工件坐標(biāo)系的轉(zhuǎn)換分兩步:

    (1)先將測量坐標(biāo)系進(jìn)行原點偏移。即將A(x0,y0,z0)作為工件坐標(biāo)系的新原點o′,只需在原坐標(biāo)系各軸減去相應(yīng)的偏移量,可得偏移矩陣S:

    (2)將坐標(biāo)軸旋轉(zhuǎn),使z軸與對中軸線單位向量重合。坐標(biāo)旋轉(zhuǎn)矩陣R在求解中軸線底面中心時已經(jīng)求出。

    4 仿真實驗

    4.1 有效性驗證

    用MATLAB軟件編寫確定螺紋中軸線的程序。用MATLAB自帶的[v,λ]=eig(X)函數(shù)求解特征向量λ對應(yīng)矩陣的特征值,v中3個列向量為λ對應(yīng)的特征向量。求解單位向量程序框圖見圖2。

    圖2 求解單位向量程序框圖Fig.2 Block diagram for solving unit vector

    用MATLAB編程求解底面中心,將螺紋點云對應(yīng)的三個軸數(shù)據(jù)向量經(jīng)過矩陣R后得到處于o-x′y′z′坐標(biāo)系的點云數(shù)據(jù),再由MATLAB自帶的X-1=inv(X)函數(shù)進(jìn)行坐標(biāo)反變換,得到原坐標(biāo)系的底面中心數(shù)據(jù)。

    圖3 求解底面中心坐標(biāo)程序框圖Fig.3 Block diagram of the bottom center coordinates

    用圓柱點云圖檢驗本文中軸線擬合算法的有效性。已知該圓柱點云中軸線參數(shù)為:

    點云數(shù)據(jù)由scatter3函數(shù)繪制,在計算完成的中軸線上取若干點用MATLAB自帶函數(shù)plot3繪制,如圖4所示。

    圖4 基于三維凸包的最小二乘擬合結(jié)果Fig.4 Results of least square fitting based on 3D convex hull

    中軸線擬合算法解得單位向量vt:

    圓柱底面中心點At:

    坐標(biāo)旋轉(zhuǎn)矩陣R:

    坐標(biāo)偏移矩陣S:

    在有限值域區(qū)間內(nèi),以一條直線上的點到另一條直線的距離的最大值作為指針,定量評定不同直線與已知直線之間的符合程度。本文以兩直線在螺紋旋合長度內(nèi),某一直線兩端點到另一條直線的距離的最大者為評定參數(shù),評價兩條直線的符合程度。其算法過程為:

    (1)首先以任意一條直線l的單位向量為法向量,取點云邊界點Pα生成上截面α;下截面法向量相同。

    (2)將上截面方程轉(zhuǎn)換為參數(shù)方程,將參數(shù)方程代入兩直線lt、lz的點斜式,可解到擬合直線lt在上截面α的端點坐標(biāo)(xlt,α,ylt,α,zlt,α),已知直線lz在上截面α的端點坐標(biāo)(xlz,α,ylz,α,zlz,α)。

    (3)同理按步驟(1)生成下截面β,擬合直線lt在下界面β的端點坐標(biāo)(xlt,β,ylt,β,zlt,β),已知直線lz在下截面β的端點坐標(biāo)(xlz,β,ylz,β,zlz,β)。

    (4)分別計算截面α,β上,兩直線端點距離Eα、Eβ,用MATLAB自帶的Max函數(shù)取兩者的較大值即為評定參數(shù)E。

    運算所得E=0.21 μm,所得中軸線參數(shù)與lz直線基本重合,算法可行。

    4.2 中軸線擬合仿真實驗

    為了比較研究中軸線擬合算法的可行性,本文用MATLAB軟件編譯基于三維凸包的最小二乘算法擬合螺紋中軸線,并與投影法比較驗證其可行性。根據(jù)國家標(biāo)準(zhǔn)《普通螺紋基本尺寸》表1[19]里的螺紋小徑、中徑數(shù)據(jù)建立了公稱直徑為80 mm,螺距為4 mm的普通外螺紋面的三維模型,并在模型表面隨機成塊刪除點云形成螺紋自身表面缺陷、掃描失真等造成的干擾因素,如圖5所示。

    圖5 帶有干擾因素的外螺紋三維模型面Fig.5 Three-dimensional model surface of external thread with interference factors

    將三維模型面保存為obj格式,在Visual Studio 2017環(huán)境下用PCL18.1編程庫將其生成外螺紋面三維點云。外螺紋面底面圓心坐標(biāo)為(-40,150,-40),頂端圓心坐標(biāo)為(-40,0,-40),如圖6所示。

    圖6 帶有干擾的外螺紋三維點云Fig.6 3D point cloud image of external thread with interference

    4.2.1 投影法擬合中軸線

    投影法是軸體測量中常見的方法,在測量時,其將螺紋器具投影至xOy面、yOz面和zOx面,基于平面測量[20]。

    本文將三維點云投影至xOy面、yOz面經(jīng)過計算邊界坐標(biāo),將xOy面、yOz面的邊界點y值取中值,分別擬合直線。求得投影面上的中線,過所得中線作面,根據(jù)2個不平行的平面相交確定1條直線的原則[13]確定空間螺紋中軸線,用MATLAB編程得中軸線一般方程:

    4.2.2 最小二乘擬合螺紋中軸線

    用隨機增量算法對帶有干擾因素的外螺紋三維點云作convex hull計算,計算結(jié)果如圖7所示。

    圖7 經(jīng)過convex hull處理后外螺紋三維點云Fig.7 Thread 3D point cloud after convex hull processed

    經(jīng)過convex hull處理后的外螺紋三維點云導(dǎo)入MATLAB,將參數(shù)區(qū)間限制為點云數(shù)值區(qū)域內(nèi),用投影法和基于三維凸包的最小二乘算法擬合中軸線。圖8、圖9分別為螺紋中軸線三維擬合結(jié)果以及其在xOy面、yOz面的投影圖。

    圖8 基于三維凸包的最小二乘算法和投影法擬合結(jié)果Fig.8 Fitting results of least squares algorithm and projection method based on three-dimensional convex hull

    圖9 擬合結(jié)果的投影Fig.9 Projection of the fitting result

    用基于三維凸包的最小二乘算法擬合的螺紋中軸線所確定的中軸線lls參數(shù)為:

    單位向量vls

    螺紋底面中心Als

    坐標(biāo)旋轉(zhuǎn)矩陣R

    坐標(biāo)偏移矩陣S

    4.3 結(jié)果驗證

    由圖8、圖9可知,投影法擬合的直線與基于三維凸包最小二乘法擬合的直線基本重合。為了驗證基于三維凸包的最小二乘算法擬合中軸線的精度和符合程度,本文采用計算lty,lls與外螺紋點云距離的方差[21]來表示算法擬合精度:

    符合程度評價標(biāo)準(zhǔn)與上節(jié)相同,在有限值域區(qū)間內(nèi),將1條線段的2個端點到另1條直線的距離的最大者作為評定標(biāo)準(zhǔn)E。以max函數(shù)、min函數(shù)對點云的x,y,z值取最大或最小,找到上下邊界點Pb1,Pb2,再按照第4.1節(jié)E算法過程求解。Xdata,Ydata,Zdata分別為點云的x,y,z數(shù)據(jù)集,上、下邊界點的生成代碼如下:

    %生成上邊界點

    max_X=max(Xdata); max_Y=max(Xdata);

    max_Z=max(Xdata);

    P_a=[max_X,max_Y,max_Z];

    %生成下邊界點

    min_X=min(Xdata); min_Y=min(Xdata);

    min_Z=min(Xdata);

    P_b=[min_X,min_Y,max_Z];

    為比較投影法和最小二乘算法的計算速度,在處理器為Intel(R)Core(TM)i5-3230M的計算機上運行,結(jié)果見表1。

    表1 擬合確定的中軸線比較表Tab.1 Table of fitted central axis variances

    從擬合結(jié)果來看,用二維投影法和最小二乘法擬合的中軸線與外螺紋點云距離的方差結(jié)果都為0.34,且兩直線在旋合長度范圍內(nèi)端面點距離的最大者,值為0.15 μm,達(dá)到了三維高精度測量標(biāo)準(zhǔn),說明最小二乘法能夠獲得投影法一致的效果;在計算機執(zhí)行速度上來看,最小二乘法擬合中軸線程序在MATLAB上運行時間只有15.3 ms,在同樣的硬件條件下,逐次求解的投影法計算速度為1 196.7 ms。由此可以得出基于三維點云的最小二乘擬合算法是一種快速、準(zhǔn)確的螺紋中軸線的擬合方法。

    綜上,最小二乘擬合中軸線的方法是從三維的角度上擬合中軸線,其囊括的維度更高更具有說服性,且計算時間短,在高精度螺紋規(guī)的測量上更具有實用性和推廣價值。

    5 結(jié) 論

    本文從螺紋中軸線特征入手,建立了螺紋中軸線數(shù)學(xué)模型,提出了基于三維凸包的最小二乘算法擬合螺紋中軸線的方法,通過MATLAB軟件編程驗證算法的可行性,最后結(jié)合三維模擬仿真實驗,并經(jīng)二維投影法驗證,得出以下結(jié)論:

    (1)經(jīng)過三維凸包計算的最小二乘算法在擬合具有三維復(fù)雜結(jié)構(gòu)的螺紋中軸線具有與投影法一致的結(jié)果。

    (2)經(jīng)過三維凸包處理的最小二乘法從三維的角度解決問題,計算速度快,在高精度螺紋測量上,具有實用性和推廣價值。

    猜你喜歡
    單位向量外螺紋中軸線
    漫畫北京中軸線(三)
    《穿越北京中軸線》簡介
    行走中軸線 尋找城市靈魂
    中國收藏(2023年6期)2023-06-08 21:13:31
    聚焦單位向量的常見題型
    單位向量用途大
    不容忽視的基本概念—單位向量
    全新的圓柱螺紋聯(lián)接防松原理及防松技術(shù)
    外螺紋冷滾壓精密成形工藝研究進(jìn)展
    平分集與球面的交集的連通性及其應(yīng)用
    更正啟事
    毛片女人毛片| 午夜激情福利司机影院| 国内揄拍国产精品人妻在线| av.在线天堂| 高清午夜精品一区二区三区 | 国产精品亚洲美女久久久| 欧美日韩国产亚洲二区| 亚洲中文日韩欧美视频| 午夜视频国产福利| a级一级毛片免费在线观看| 国产精品野战在线观看| 国产色婷婷99| 国产精品电影一区二区三区| 中文字幕av在线有码专区| av在线亚洲专区| 看黄色毛片网站| 日韩成人伦理影院| 精品久久久久久久人妻蜜臀av| 亚洲欧美成人精品一区二区| 伦精品一区二区三区| 久久精品国产亚洲av涩爱 | 内射极品少妇av片p| 亚洲欧美日韩东京热| 美女 人体艺术 gogo| 校园人妻丝袜中文字幕| 女人十人毛片免费观看3o分钟| 亚洲av熟女| 六月丁香七月| 亚洲精品国产av成人精品 | .国产精品久久| 青春草视频在线免费观看| 日韩一区二区视频免费看| 精品久久久久久久人妻蜜臀av| 中文字幕免费在线视频6| 免费看日本二区| 亚洲av免费高清在线观看| 久久九九热精品免费| 亚洲国产精品成人久久小说 | 国产69精品久久久久777片| 极品教师在线视频| 别揉我奶头 嗯啊视频| 男女那种视频在线观看| av专区在线播放| 看片在线看免费视频| 欧美最黄视频在线播放免费| 波多野结衣巨乳人妻| 伦理电影大哥的女人| 老熟妇乱子伦视频在线观看| 色综合色国产| 不卡一级毛片| 男女啪啪激烈高潮av片| 少妇人妻一区二区三区视频| 亚洲美女视频黄频| 变态另类丝袜制服| 噜噜噜噜噜久久久久久91| 日本a在线网址| 午夜福利在线观看吧| 99久国产av精品国产电影| 一进一出抽搐动态| 久久精品国产99精品国产亚洲性色| 亚洲av电影不卡..在线观看| 九九爱精品视频在线观看| 国产又黄又爽又无遮挡在线| 少妇的逼水好多| 搞女人的毛片| av福利片在线观看| 免费在线观看影片大全网站| 日日啪夜夜撸| 国产不卡一卡二| 亚州av有码| 久久精品综合一区二区三区| 日韩欧美免费精品| 亚洲五月天丁香| 老司机影院成人| 99热这里只有是精品在线观看| 欧美成人免费av一区二区三区| 成人综合一区亚洲| 伦精品一区二区三区| 国产视频一区二区在线看| 欧美+亚洲+日韩+国产| 一a级毛片在线观看| 亚洲av美国av| 欧美性感艳星| 久久午夜福利片| 可以在线观看的亚洲视频| 中出人妻视频一区二区| 又黄又爽又刺激的免费视频.| 国产亚洲精品av在线| h日本视频在线播放| 在线国产一区二区在线| 亚洲国产精品久久男人天堂| 在线免费观看不下载黄p国产| 久久天躁狠狠躁夜夜2o2o| 国产一区二区三区在线臀色熟女| 男人舔女人下体高潮全视频| 国产爱豆传媒在线观看| 看十八女毛片水多多多| 啦啦啦观看免费观看视频高清| 夜夜夜夜夜久久久久| 插阴视频在线观看视频| 美女xxoo啪啪120秒动态图| 久久久久久伊人网av| 亚洲av中文av极速乱| h日本视频在线播放| 久久久久久久久中文| 九九爱精品视频在线观看| 午夜爱爱视频在线播放| av在线播放精品| 亚洲18禁久久av| 日本色播在线视频| 99国产精品一区二区蜜桃av| 亚洲国产精品成人综合色| 在线免费十八禁| 久久精品国产99精品国产亚洲性色| 国产 一区 欧美 日韩| 麻豆一二三区av精品| 国产白丝娇喘喷水9色精品| 日韩欧美 国产精品| 亚洲成人中文字幕在线播放| 亚洲一级一片aⅴ在线观看| 国产一区二区三区在线臀色熟女| 欧美成人a在线观看| 可以在线观看的亚洲视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲在线观看片| 久久久久久久久久成人| 日本黄色视频三级网站网址| 久久人人精品亚洲av| 香蕉av资源在线| 激情 狠狠 欧美| 欧美成人一区二区免费高清观看| 国产精品美女特级片免费视频播放器| 午夜精品国产一区二区电影 | 亚洲精品粉嫩美女一区| 久久韩国三级中文字幕| 一进一出抽搐动态| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| 成人亚洲欧美一区二区av| av天堂在线播放| 亚洲一级一片aⅴ在线观看| 一级毛片久久久久久久久女| 国产亚洲av嫩草精品影院| 色视频www国产| 亚洲人成网站高清观看| 一个人看视频在线观看www免费| 1000部很黄的大片| 菩萨蛮人人尽说江南好唐韦庄 | 国产乱人偷精品视频| 国产伦在线观看视频一区| 精品熟女少妇av免费看| 精品久久久久久久久亚洲| 91久久精品国产一区二区三区| 天天躁日日操中文字幕| 国产日本99.免费观看| 国产精品一区二区免费欧美| 免费看光身美女| 在线观看av片永久免费下载| 久久99热这里只有精品18| 国产精品爽爽va在线观看网站| 18禁裸乳无遮挡免费网站照片| 成人鲁丝片一二三区免费| 欧美一级a爱片免费观看看| 日韩精品中文字幕看吧| 此物有八面人人有两片| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av香蕉五月| 午夜福利18| 久久热精品热| 91久久精品电影网| 夜夜看夜夜爽夜夜摸| 波野结衣二区三区在线| 久久久久久久久久成人| 亚洲性久久影院| 亚洲中文日韩欧美视频| av专区在线播放| 国产精品电影一区二区三区| 精品久久久噜噜| 国产大屁股一区二区在线视频| 精品日产1卡2卡| 亚洲欧美精品综合久久99| 国产精品久久久久久久久免| 干丝袜人妻中文字幕| 亚洲在线观看片| 国产精品亚洲一级av第二区| 色噜噜av男人的天堂激情| 国产精品爽爽va在线观看网站| 老司机影院成人| 老女人水多毛片| 黄色日韩在线| 国产精品无大码| 国产一区二区亚洲精品在线观看| 最近的中文字幕免费完整| 最近最新中文字幕大全电影3| 久久鲁丝午夜福利片| 91av网一区二区| 国产高清有码在线观看视频| 国产淫片久久久久久久久| 日本在线视频免费播放| 中国美白少妇内射xxxbb| 中文字幕av成人在线电影| 日韩三级伦理在线观看| 卡戴珊不雅视频在线播放| 欧美色欧美亚洲另类二区| 免费人成在线观看视频色| 国产精品久久久久久久电影| 五月伊人婷婷丁香| 我要看日韩黄色一级片| 最近2019中文字幕mv第一页| 精品久久久噜噜| 国产熟女欧美一区二区| 搡女人真爽免费视频火全软件 | 国产精品伦人一区二区| 亚洲人成网站在线观看播放| 免费观看在线日韩| 99久久精品热视频| av卡一久久| 99热这里只有精品一区| 三级经典国产精品| 国产精品综合久久久久久久免费| 国产高清视频在线观看网站| 国产真实乱freesex| 亚洲中文日韩欧美视频| 婷婷精品国产亚洲av在线| 观看美女的网站| 久久精品国产亚洲av天美| 国内精品美女久久久久久| 97在线视频观看| 国产精品伦人一区二区| .国产精品久久| 国产男人的电影天堂91| 能在线免费观看的黄片| 91狼人影院| 人人妻人人澡欧美一区二区| 女生性感内裤真人,穿戴方法视频| 给我免费播放毛片高清在线观看| 久久久久精品国产欧美久久久| 亚洲av二区三区四区| 一a级毛片在线观看| 欧美zozozo另类| 精品午夜福利视频在线观看一区| 日日摸夜夜添夜夜添av毛片| 又粗又爽又猛毛片免费看| 麻豆国产97在线/欧美| 国产淫片久久久久久久久| 欧美一级a爱片免费观看看| 亚洲电影在线观看av| 丰满人妻一区二区三区视频av| 中国美女看黄片| 日韩欧美国产在线观看| 国产高清不卡午夜福利| 久久久成人免费电影| 成人无遮挡网站| 中文在线观看免费www的网站| 一个人看的www免费观看视频| 亚洲欧美成人综合另类久久久 | 12—13女人毛片做爰片一| 国产黄a三级三级三级人| 国产真实乱freesex| 成人精品一区二区免费| 97碰自拍视频| 中文字幕熟女人妻在线| aaaaa片日本免费| 欧美日本视频| 午夜影院日韩av| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 国内精品一区二区在线观看| 少妇丰满av| 国产精品久久视频播放| 在线播放国产精品三级| 在线观看免费视频日本深夜| 国产精品,欧美在线| 久久九九热精品免费| 看片在线看免费视频| 波多野结衣高清无吗| av在线天堂中文字幕| 干丝袜人妻中文字幕| 伊人久久精品亚洲午夜| 一a级毛片在线观看| videossex国产| av在线老鸭窝| 三级国产精品欧美在线观看| 国产国拍精品亚洲av在线观看| 午夜激情福利司机影院| 日日摸夜夜添夜夜添av毛片| 国产精品无大码| 国产成人a∨麻豆精品| 午夜久久久久精精品| 日本黄色视频三级网站网址| 国产精品三级大全| 亚洲成a人片在线一区二区| 成人高潮视频无遮挡免费网站| 色噜噜av男人的天堂激情| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人精品一区久久| 国产高清视频在线播放一区| 亚洲精品国产成人久久av| 成人鲁丝片一二三区免费| 黄色一级大片看看| 日韩欧美三级三区| 亚洲国产色片| 亚洲国产日韩欧美精品在线观看| 亚洲欧美日韩无卡精品| 一个人观看的视频www高清免费观看| 亚洲第一电影网av| 在线观看免费视频日本深夜| 精品午夜福利在线看| 久久精品夜色国产| 特级一级黄色大片| 热99在线观看视频| 美女cb高潮喷水在线观看| 日韩一区二区视频免费看| 亚洲自拍偷在线| 性插视频无遮挡在线免费观看| av专区在线播放| 免费看a级黄色片| 一区二区三区四区激情视频 | 国产在线精品亚洲第一网站| 99久久久亚洲精品蜜臀av| 亚洲最大成人av| 国产精品伦人一区二区| 女同久久另类99精品国产91| 哪里可以看免费的av片| 日本黄色片子视频| 看十八女毛片水多多多| 深夜精品福利| 热99re8久久精品国产| 欧美日韩国产亚洲二区| 午夜a级毛片| 亚洲四区av| 久久亚洲国产成人精品v| 最新中文字幕久久久久| 性色avwww在线观看| 欧美潮喷喷水| 99国产极品粉嫩在线观看| 一区二区三区四区激情视频 | 亚洲三级黄色毛片| 亚洲精品在线观看二区| 免费人成视频x8x8入口观看| 精品久久久久久久久久免费视频| 久久精品人妻少妇| 日韩成人av中文字幕在线观看 | 九九久久精品国产亚洲av麻豆| 最近在线观看免费完整版| 久久久久免费精品人妻一区二区| 成人国产麻豆网| 精品一区二区三区视频在线| 欧美人与善性xxx| 久久精品国产清高在天天线| 久久人人爽人人片av| 婷婷亚洲欧美| 国产国拍精品亚洲av在线观看| 日本黄色片子视频| 在线观看66精品国产| 欧美一区二区国产精品久久精品| 男人狂女人下面高潮的视频| 国产高清不卡午夜福利| 在线观看免费视频日本深夜| 麻豆国产97在线/欧美| 一级av片app| 亚洲第一区二区三区不卡| 俄罗斯特黄特色一大片| 亚洲国产精品久久男人天堂| 日本免费一区二区三区高清不卡| 一区二区三区四区激情视频 | 亚洲丝袜综合中文字幕| eeuss影院久久| 久久午夜亚洲精品久久| 亚洲无线观看免费| 菩萨蛮人人尽说江南好唐韦庄 | 白带黄色成豆腐渣| 哪里可以看免费的av片| 91精品国产九色| 不卡一级毛片| 日韩成人av中文字幕在线观看 | 国产精品国产三级国产av玫瑰| aaaaa片日本免费| 一级黄色大片毛片| 亚洲最大成人av| 国产综合懂色| 午夜影院日韩av| 国产亚洲精品久久久com| 欧美人与善性xxx| 国产男人的电影天堂91| 国产精品伦人一区二区| 亚洲一区二区三区色噜噜| 乱人视频在线观看| 国产aⅴ精品一区二区三区波| 人人妻人人澡欧美一区二区| 成人无遮挡网站| 91狼人影院| 亚洲av美国av| 亚洲精品456在线播放app| 午夜福利在线观看免费完整高清在 | 一个人看视频在线观看www免费| 国产精品,欧美在线| 最新在线观看一区二区三区| 亚洲第一区二区三区不卡| 一个人免费在线观看电影| 国产成人a区在线观看| 97人妻精品一区二区三区麻豆| 最近在线观看免费完整版| 婷婷亚洲欧美| 亚洲性夜色夜夜综合| 国产大屁股一区二区在线视频| 亚洲乱码一区二区免费版| 午夜精品在线福利| 婷婷亚洲欧美| 综合色av麻豆| 乱码一卡2卡4卡精品| 蜜臀久久99精品久久宅男| 亚洲欧美清纯卡通| 国产成人一区二区在线| 亚洲国产欧洲综合997久久,| 一卡2卡三卡四卡精品乱码亚洲| 老熟妇仑乱视频hdxx| 少妇猛男粗大的猛烈进出视频 | 中文字幕精品亚洲无线码一区| 在线观看66精品国产| 俄罗斯特黄特色一大片| 亚洲欧美日韩高清在线视频| 久久精品久久久久久噜噜老黄 | 久久亚洲精品不卡| 99热全是精品| 亚洲中文字幕一区二区三区有码在线看| 久久综合国产亚洲精品| 午夜精品一区二区三区免费看| eeuss影院久久| 美女 人体艺术 gogo| 国产极品精品免费视频能看的| av视频在线观看入口| 午夜日韩欧美国产| 日本精品一区二区三区蜜桃| 欧美成人精品欧美一级黄| 搞女人的毛片| 99久久中文字幕三级久久日本| 久久99热6这里只有精品| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 中文在线观看免费www的网站| 夜夜看夜夜爽夜夜摸| 日韩av在线大香蕉| 久久精品国产鲁丝片午夜精品| .国产精品久久| 我要搜黄色片| 午夜视频国产福利| 亚州av有码| 欧美国产日韩亚洲一区| 亚洲欧美日韩无卡精品| 国产人妻一区二区三区在| 在线a可以看的网站| 最近中文字幕高清免费大全6| 日日摸夜夜添夜夜添小说| 99久久中文字幕三级久久日本| 精品午夜福利视频在线观看一区| 亚洲精品亚洲一区二区| 亚洲欧美精品综合久久99| 日韩欧美 国产精品| 亚洲精品久久国产高清桃花| 日本熟妇午夜| 丰满乱子伦码专区| 国产精品一区二区免费欧美| avwww免费| 国产视频内射| 欧美一区二区精品小视频在线| 99riav亚洲国产免费| 看免费成人av毛片| 日韩强制内射视频| 啦啦啦观看免费观看视频高清| 99热网站在线观看| 午夜福利成人在线免费观看| 人人妻,人人澡人人爽秒播| 国产人妻一区二区三区在| 成人特级黄色片久久久久久久| 日韩av在线大香蕉| 国产伦一二天堂av在线观看| 成年女人永久免费观看视频| 成人美女网站在线观看视频| av国产免费在线观看| 成年免费大片在线观看| 亚洲激情五月婷婷啪啪| 好男人在线观看高清免费视频| 欧美日韩综合久久久久久| а√天堂www在线а√下载| 国产亚洲精品av在线| 91久久精品国产一区二区三区| 亚洲三级黄色毛片| 欧美精品国产亚洲| 日本一本二区三区精品| 九九久久精品国产亚洲av麻豆| 91精品国产九色| 国产av在哪里看| 久久精品夜色国产| 国产爱豆传媒在线观看| 久久九九热精品免费| 一级黄色大片毛片| 久久精品综合一区二区三区| 男女那种视频在线观看| 最新中文字幕久久久久| 久久久国产成人精品二区| 亚洲真实伦在线观看| 久99久视频精品免费| 身体一侧抽搐| 亚洲最大成人手机在线| 欧美最新免费一区二区三区| 免费人成视频x8x8入口观看| 国产探花极品一区二区| 免费大片18禁| 国产精品不卡视频一区二区| 男女边吃奶边做爰视频| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 又爽又黄a免费视频| 久久精品国产亚洲av香蕉五月| 欧美成人免费av一区二区三区| 国产一区二区三区av在线 | 成人二区视频| 亚洲乱码一区二区免费版| 日本撒尿小便嘘嘘汇集6| 午夜福利在线在线| 一级黄片播放器| 人人妻,人人澡人人爽秒播| 国产精品美女特级片免费视频播放器| 日韩欧美 国产精品| 国产成人aa在线观看| 91久久精品电影网| 人妻少妇偷人精品九色| 国产女主播在线喷水免费视频网站 | 亚洲无线在线观看| 看片在线看免费视频| 国产精品三级大全| 国产激情偷乱视频一区二区| 亚洲内射少妇av| 国产精品电影一区二区三区| 哪里可以看免费的av片| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 亚洲av成人精品一区久久| 国产成年人精品一区二区| 精品一区二区三区av网在线观看| 免费不卡的大黄色大毛片视频在线观看 | 午夜激情欧美在线| 男人舔奶头视频| 韩国av在线不卡| 蜜臀久久99精品久久宅男| 国产精品一二三区在线看| 最好的美女福利视频网| 又爽又黄无遮挡网站| 91av网一区二区| 久久国产乱子免费精品| 欧美色欧美亚洲另类二区| 亚洲乱码一区二区免费版| 女人被狂操c到高潮| 日本五十路高清| 干丝袜人妻中文字幕| 日本色播在线视频| 校园春色视频在线观看| 久久亚洲精品不卡| 国产av一区在线观看免费| av免费在线看不卡| 国产色爽女视频免费观看| 日韩三级伦理在线观看| 欧美一级a爱片免费观看看| 精品久久久久久久久av| 久久99热这里只有精品18| 成年免费大片在线观看| 国产 一区精品| 啦啦啦啦在线视频资源| 欧美日韩在线观看h| 亚洲欧美精品综合久久99| 国内揄拍国产精品人妻在线| 在线播放国产精品三级| 亚洲av免费在线观看| 啦啦啦韩国在线观看视频| 久久久精品94久久精品| 国产亚洲精品av在线| 男女之事视频高清在线观看| 日韩,欧美,国产一区二区三区 | 欧美三级亚洲精品| 在线免费十八禁| 久久久午夜欧美精品| 99热这里只有是精品在线观看| 国产午夜精品久久久久久一区二区三区 | 女人被狂操c到高潮| 亚洲精品成人久久久久久| 免费黄网站久久成人精品| 一进一出好大好爽视频| 真实男女啪啪啪动态图| 有码 亚洲区| 搡老熟女国产l中国老女人| 一个人免费在线观看电影| 国产精品一区二区三区四区免费观看 | 晚上一个人看的免费电影| 国产一区二区在线观看日韩| 日韩欧美国产在线观看| 天堂√8在线中文| 欧美日本视频| 免费人成视频x8x8入口观看| 九九在线视频观看精品| 日韩亚洲欧美综合| 97碰自拍视频| 亚洲精品456在线播放app| 国产成人91sexporn| 亚洲三级黄色毛片| 中文字幕熟女人妻在线| 成人高潮视频无遮挡免费网站| 久久午夜福利片| 欧美一区二区精品小视频在线| a级毛色黄片| 久久精品久久久久久噜噜老黄 | 露出奶头的视频| 久久久久国产精品人妻aⅴ院| 性色avwww在线观看| 一进一出抽搐gif免费好疼|