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

    基于NCGM的粗糙表面數(shù)值模擬與實驗對比

    2014-12-05 06:54:46唐進元廖東日
    中國機械工程 2014年14期
    關(guān)鍵詞:偏態(tài)高斯分布長度

    唐進元 廖東日 周 煒

    中南大學(xué)高性能復(fù)雜制造國家重點實驗室,長沙,410083

    0 引言

    很多數(shù)值方法已用于粗糙表面的模擬。Patir[1]采用線性矩陣變換方法模擬具有指定相關(guān)函數(shù)的表面,計算自相關(guān)矩陣時采用牛頓迭代法求解相關(guān)系數(shù)矩陣,該方法的主要缺點是所需的儲存空間很大,計算相當(dāng)費時,而且其收斂性也一直備受爭議。文獻[2-3]采用時間序列模型來模擬粗糙表面,采用的是自回歸滑動平均模型(ARMA)。Whitehouse[4]采用自回歸模型(AR)來模擬粗糙表面,但是ARMA與MA模型只考慮了低階系統(tǒng),而且只描述了原點附近的自相關(guān)函數(shù)。針對儲存空間和計算速度的限制,快速傅里葉變換(FFT)方法也被用于粗糙表面的模擬。文獻[5-7]采 用 FFT 方 法 模 擬 粗 糙 表 面。Wu[6-7]指出,采用FFT方法模擬粗糙表面不能同時滿足正確的自相關(guān)函數(shù)與循環(huán)自相關(guān)函數(shù),但是可以保證平均自相關(guān)函數(shù)的正確性,F(xiàn)FT方法只能應(yīng)用于自相關(guān)長度較小的粗糙表面的模擬。基于上述文獻的工作與Partir方法,為解決儲存空間與計算時間的限制問題,本文采用非線性共軛梯度法[8](NCGM)計算自相關(guān)矩陣,模擬非高斯型粗糙表面時采用Johnson轉(zhuǎn)換系統(tǒng)[9-11]產(chǎn)生所需的非高斯序列。

    1 粗糙表面統(tǒng)計規(guī)律

    自相關(guān)函數(shù)可用于描述不同位置之間數(shù)據(jù)的相關(guān)程度,該函數(shù)提供了數(shù)據(jù)最基本的空間信息和數(shù)據(jù)相互之間的關(guān)系。因此,所有空間上的表面粗糙度參數(shù)都是從自相關(guān)函數(shù)獲取的。

    三維表面的自相關(guān)函數(shù)形式如下:

    其中,E是數(shù)學(xué)期望值,τx、τy分別是x 方向、y方向的相關(guān)距離,R(0,0)即為σ2(σ為均方根)。

    其離散表達式為

    其中,n、m分別為自相關(guān)長度在x、y方向上的最大值;N、M分別為離散點沿x、y方向的總數(shù)。

    相關(guān)長度τ*定義為在輪廓方向上自相關(guān)函數(shù)衰減到臨界值時的長度,Stout等[12]定義自相關(guān)函數(shù)從原點衰減到其值的10%時的長度為相關(guān)長度,此時相關(guān)函數(shù)的表達式如下:

    2 粗糙表面數(shù)值模擬

    按高度分布函數(shù)的不同,粗糙表面可分為高斯型粗糙表面與非高斯型粗糙表面,兩種粗糙表面生成的理論模型均為時間序列模型中的滑動平均模型(MA)。

    2.1 高斯型粗糙表面模擬

    基于Patir[1]的理論,生成(M,N)的粗糙表面需先產(chǎn)生(M+m,N+n)的獨立分布的隨機數(shù)和(m,n)維的自相關(guān)矩陣,生成粗糙表面的線性變換方程如下:

    其中,ak,l是欲生成自相關(guān)函數(shù)的自相關(guān)系數(shù)矩陣第k行第l列元素,ηi,j是獨立的單位方差隨機數(shù),滿足以下關(guān)系式:

    自相關(guān)函數(shù)的定義如下:

    式(6)是求解維數(shù)為m×n的自相關(guān)系數(shù)矩陣的元素ak,l的非線性方程組,Patir[1]采用牛頓法求解該方程組,給定迭代的近似初值為

    當(dāng)用牛頓法求解式(6)時,Bakolas[13]發(fā)現(xiàn),當(dāng)矩陣維數(shù)一定時,該方程組不收斂。除收斂性問題外,求解式(6)的計算量相當(dāng)大,還需要很大的空間來存儲雅可比矩陣。

    本文采用NCGM來求解該非線性方程組,把式(6)改寫成:

    該方法把求解原函數(shù)替換成求解原函數(shù)的一階導(dǎo)數(shù)。按照NCGM,現(xiàn)有兩種搜索方向選取方法,一種為Fletcher-Reeves方法,另一種為Polak-Ribiere方法,其表達式分別為

    當(dāng)初值選取為非常接近目標(biāo)函數(shù)的解時,F(xiàn)letcher-Reeves方法會收斂,而Polak-Ribiere方法在極少數(shù)情況下會不收斂,但是Polak-Ribiere方法的收斂速度比前者快很多。在迭代過程中,為避免求函數(shù)導(dǎo)數(shù),采用割線法來進行非精確線搜索,其中割線法的步長需不斷地調(diào)整來達到收斂,本文采用Polak-Ribiere方法求解方程組,求解步驟如下:

    (2)d(0)=r(0)-f′(x(0))(d為搜索方向)。

    (3)x(i+1)= x(i)+α(i)d(i),找 到 令 f(x(i)+α(i)d(i))最小的α(i)。

    (4)r(i+1)=-f′(x(i+1))。

    (6)d(i+1)=r(i+1)+β(i+1)d(i)。

    指定自相關(guān)函數(shù)的粗糙表面模擬步驟如下:

    (1)利用計算機產(chǎn)生一高斯分布的隨機序列μi,j。

    (3)根據(jù)式(10),通過NCGM 求解得出自相關(guān)系數(shù)矩陣。

    (4)根據(jù)式(4)計算得出模擬的粗糙表面z(x,y)。

    基于以上理論和方法,圖1和圖2分別給出了相關(guān)長度βx=βy=3μm和βx=βy=10μm時的粗糙表面,粗糙表面的均值均為0,標(biāo)準(zhǔn)差為1,模擬粗糙表面點數(shù)M=N=60,x與y方向采樣間隔均為1μm??煽闯霎?dāng)相關(guān)長度變大時表面形貌更加緩和。當(dāng)βx和βy相等時,為各向同性表面。

    圖1 相關(guān)長度βx=βy=3μm時的高斯分布表面

    當(dāng)βx和βy不相等時為各向異性表面,圖3給出了βx=30μm,βy=3μm時的各向異性粗糙表面。當(dāng)相關(guān)長度不同且數(shù)值上差別較大時,可明顯看出紋理的走向。

    圖2 相關(guān)長度βx=βy=10μm時的高斯分布表面

    圖3 相關(guān)長度βx=30μm、βy=3μm時的高斯分布表面

    2.2 非高斯型粗糙表面模擬

    實際工程中,工程表面大多為非高斯分布表面,表面上的峰相對谷來說比較容易被除去,所以很多粗糙表面呈負(fù)偏態(tài),故相對而言負(fù)偏態(tài)表面的接觸性能較好。因此,為了更加真實地模擬粗糙表面,非高斯型表面的模擬尤為重要。粗糙表面分布的偏態(tài)為

    粗糙表面的峰態(tài)為

    非高斯型表面的模擬相對高斯型表面來說,其輸入隨機序列μi,j與高斯型表面模擬時不同,模擬高斯型表面時,μi,j是高斯分布隨機序列;而模擬非高斯型表面時,μi,j為非高斯分布隨機序列,所以必須產(chǎn)生一個μi,j來滿足需模擬表面的偏態(tài)與峰態(tài)。式(4)本質(zhì)上是一個MA模型,對于一個純MA模型,有

    當(dāng)偏態(tài)與峰態(tài)通過一個數(shù)字濾波器時,其數(shù)值會改變,所以必須對輸入的偏態(tài)與峰態(tài)進行修正,其輸入與輸出存在以下對應(yīng)關(guān)系:

    其中Sη、Kη分別為輸入序列的偏態(tài)與峰態(tài)。Sz、Kz分別為輸出序列的偏態(tài)與峰態(tài),即為所求參數(shù)。根據(jù)式(4),得

    把高斯分布的隨機序列轉(zhuǎn)換為指定偏態(tài)與峰態(tài)非高斯隨機序列需使用Johnson轉(zhuǎn)換系統(tǒng),按照偏態(tài)與峰態(tài)取值范圍的不同,Johnson轉(zhuǎn)換系統(tǒng)可分為SU、SL、SB三種不同的形式,在不同的偏態(tài)與峰態(tài)下有以下公式對應(yīng):

    其中,SU為無界系統(tǒng),SL為對數(shù)正態(tài)系統(tǒng),SB為有界系統(tǒng)。x是標(biāo)準(zhǔn)正態(tài)分布的隨機序列,η是所得到的序列。γ、δ、ξ、λ是系統(tǒng)得到所需偏態(tài)與峰態(tài)的4個參數(shù)。

    非高斯型粗糙表面按以下步驟模擬:

    (1)求解非線性方程組(式(10))。

    (2)產(chǎn)生高斯分布序列。

    (3)通過式(16)和式(17)得出指定生成的偏態(tài)與峰態(tài)序列。

    (4)利用Johnson轉(zhuǎn)換系統(tǒng)和對應(yīng)系統(tǒng)的公式產(chǎn)生非高斯序列。

    (5)根據(jù)式(4)產(chǎn)生非高斯型粗糙表面。

    模擬非高斯型表面時,為了更加接近真實表面的形貌分布,采用實際加工的表面進行模擬,故自相關(guān)矩陣由真實表面計算得出。該工件采用超聲磨削,磨削進給量為0.01mm,超聲振動振幅為4μm,頻率為20kHz,進給速度為69mm/s,砂輪轉(zhuǎn)速為1440r/min,砂輪型號為 WA40L,工件材料為45鋼,低溫回火處理。采用Veeco Wyko NT9100光學(xué)輪廓測量儀進行表面微觀形貌測量,采樣間隔為0.5μm,采樣點數(shù)為640×480。為減少計算量,截取其測量表面點數(shù)為60×60進行分析重構(gòu),經(jīng)測量計算,x與y方向采樣間隔均為4μm。

    按非高斯表面模擬方法,實驗掃描完整表面,實際測量區(qū)域表面與模擬的表面如圖4~圖6所示??煽闯龀暷ハ骷庸け砻嬗忻黠@的加工輪廓。經(jīng)計算,相關(guān)長度在y方向非常大,而x方向上的相關(guān)長度很小,完全符合加工方向的紋理分布。

    圖4 實驗測量表面

    圖5 實際測量區(qū)域表面

    圖6 模擬表面

    3 粗糙表面模擬誤差分析

    粗糙表面模擬與理論模型要求的誤差可從高度分布函數(shù)與自相關(guān)函數(shù)的對比上得出。其中高度分布函數(shù)統(tǒng)計主要進行均值、方差、偏態(tài)和峰態(tài)四個參數(shù)的對比,而相關(guān)函數(shù)則采用實測表面相關(guān)函數(shù)與數(shù)值模擬表面相關(guān)函數(shù)進行對比。

    3.1 表面高度分布統(tǒng)計參數(shù)誤差分析

    實際加工表面的算術(shù)平均偏差Ra為0,均方根偏差Rq為6.12μm,偏態(tài)Rsk為-0.39μm,峰態(tài)Rku為3.30μm。實際加工表面與重構(gòu)表面統(tǒng)計參數(shù)與概率密度分布曲線對比如表1、圖7和圖8所示??煽闯鰞煞N表面在統(tǒng)計參數(shù)上可以很好地符合理論模型的要求。

    表1 實際加工表面與重構(gòu)表面統(tǒng)計參數(shù)對比 μm

    圖7 實際加工表面概率密度分布

    圖8 重構(gòu)表面概率密度分布

    3.2 表面相關(guān)函數(shù)誤差分析

    相關(guān)函數(shù)反映的是表面各點之間的相關(guān)程度,相關(guān)函數(shù)和高度分布統(tǒng)計參數(shù)能完整地表征表面形貌。由于y方向相關(guān)長度很大,是主要紋理方向,故只需對比y方向自相關(guān)函數(shù)即可。實際加工表面與重構(gòu)表面的自相關(guān)函數(shù)對比如圖9所示,可以看出重構(gòu)的表面與實際加工表面的自相關(guān)函數(shù)能較好地吻合。

    圖9 實際加工表面與重構(gòu)表面的相關(guān)函數(shù)

    4 結(jié)語

    表面高度的均值、標(biāo)準(zhǔn)差、偏態(tài)、峰態(tài)是描述高度分布函數(shù)的統(tǒng)計參數(shù),而表面的相關(guān)性則用自相關(guān)函數(shù)來描述。本文采用MA模型和NCGM,能很好地克服牛頓法所帶來的儲存數(shù)據(jù)大、無法收斂的困難,可模擬出高度概率分布函數(shù)和自相關(guān)函數(shù)指定要求的粗糙表面。綜合使用NCGM和Johnson轉(zhuǎn)換系統(tǒng)能有效地模擬出工程表面,從而建立高精度的粗糙表面模型。

    [1]Patir N.A Numerical Procedure for Random Generation of Rough Surfaces[J].Wear,1978,47(2):263-277.

    [2]Watson W,King T,Spedding T,et al.The Machined Surface-time Series Modelling[J].Wear,1979,57(1):195-205.

    [3]Watson W,Spedding T.The Time Series Modelling of Non-Gaussian Engineering Processes[J].Wear,1982,83(2):215-231.

    [4]Whitehouse D.The Generation of Two Dimensional Random Surfaces Having a Specified Function[J].CIRP Annals-Manufacturing Technology,1983,32(1):495-498.

    [5]Hu Y,Tonder K.Simulation of 3-D Random Rough Surface by 2-D Digital Filter and Fourier Analysis[J].International Journal of Machine Tools and Manufacture,1992,32(1/2):83-90.

    [6]Wu J J.Simulation of Non-Gaussian Surfaces with FFT[J].Tribology international,2004,37(4):339-346.

    [7]Wu J J.Simulation of Rough Surfaces with FFT[J].Tribology International,2000,33(1):47-58.

    [8]Manesh K,Ramamoorthy B,Singaperumal M.Numerical Generation of Anisotropic 3DNon-Gaussian Engineering Surfaces with Specified 3DSurface Roughness Parameters[J].Wear,2010,268(11/12):1371-1379.

    [9]Hill I,Hill R,Holder R.Algorithm AS 99:Fitting Johnson Curves by Moments[J].Journal of the Royal Statistical Society Series C(Applied Statistics),1976,25(2):180-189.

    [10]Hill I.Algorithm AS 100:Normal-Johnson and Johnson-normal Transformations[J].Journal of the Royal Statistical Society Series C(Applied Statistics),1976,25(2):190-192.

    [11]Johnson N L.Systems of Frequency Curves Generated by Methods of Translation[J].Biometrika,1949,36(1/2):149-176.

    [12]Stout K J,Sullivan P J,Dong W P,et al.The Development of Methods for the Characterization of Roughness on Three Dimensions:EUR,0015178 EN[P].1994-01-01.

    [13]Bakolas V.Numerical Generation of Arbitrarily Oriented Non-Gaussian Three-dimensional Rough Surfaces[J].Wear,2003,254(5/6):546-554.

    猜你喜歡
    偏態(tài)高斯分布長度
    利用Box-Cox變換對移動通信中小區(qū)級業(yè)務(wù)流量分布的研究
    考試成績轉(zhuǎn)換成偏態(tài)分布量化成績的算法
    1米的長度
    2種非對稱廣義高斯分布模型的構(gòu)造
    隧穿量子點分子的Wigner-Yanase偏態(tài)信息
    愛的長度
    怎樣比較簡單的長度
    一種基于改進混合高斯模型的前景檢測
    不同長度
    讀寫算(上)(2015年6期)2015-11-07 07:17:55
    雙模壓縮真空態(tài)光場作用下耦合雙原子的Wigner-Yanase偏態(tài)信息
    成年人午夜在线观看视频| 中文字幕人妻熟女乱码| 亚洲人成网站在线观看播放| 亚洲欧美中文字幕日韩二区| 国产精品 欧美亚洲| 亚洲激情五月婷婷啪啪| 蜜桃在线观看..| 国产亚洲一区二区精品| 91在线精品国自产拍蜜月| 国产一区二区三区av在线| 亚洲久久久国产精品| 亚洲成人一二三区av| 啦啦啦在线免费观看视频4| 美女主播在线视频| 波野结衣二区三区在线| 天天影视国产精品| 国产精品免费视频内射| 天堂8中文在线网| 免费播放大片免费观看视频在线观看| 美女脱内裤让男人舔精品视频| 日韩视频在线欧美| 人妻少妇偷人精品九色| 国产日韩一区二区三区精品不卡| 波多野结衣一区麻豆| 亚洲人成77777在线视频| 老女人水多毛片| 不卡av一区二区三区| 亚洲欧洲日产国产| 中文字幕色久视频| 亚洲综合色惰| 精品人妻熟女毛片av久久网站| 黄片播放在线免费| 亚洲四区av| 日本vs欧美在线观看视频| 亚洲av.av天堂| 欧美少妇被猛烈插入视频| 免费观看无遮挡的男女| 国产老妇伦熟女老妇高清| 老汉色∧v一级毛片| 国产精品一二三区在线看| 久久精品久久久久久久性| 中文乱码字字幕精品一区二区三区| 国产成人免费观看mmmm| 免费人妻精品一区二区三区视频| 丝袜美腿诱惑在线| a级片在线免费高清观看视频| 老熟女久久久| 秋霞在线观看毛片| 久久午夜福利片| 成人18禁高潮啪啪吃奶动态图| xxx大片免费视频| 中文字幕最新亚洲高清| 黄片播放在线免费| 日本av手机在线免费观看| 久久久久久久久免费视频了| 大片免费播放器 马上看| 国产1区2区3区精品| 免费在线观看视频国产中文字幕亚洲 | 男人爽女人下面视频在线观看| 中文字幕精品免费在线观看视频| 丰满饥渴人妻一区二区三| 久久久久久人妻| 一区福利在线观看| 国产一区二区激情短视频 | 多毛熟女@视频| 久久 成人 亚洲| 最近的中文字幕免费完整| 晚上一个人看的免费电影| 欧美av亚洲av综合av国产av | 老汉色∧v一级毛片| 婷婷色av中文字幕| 欧美日韩精品网址| 老司机亚洲免费影院| 久久这里只有精品19| 亚洲天堂av无毛| 久久免费观看电影| 性高湖久久久久久久久免费观看| 日韩欧美精品免费久久| 婷婷色综合www| 性色av一级| 丝袜在线中文字幕| 日韩免费高清中文字幕av| 日韩大片免费观看网站| 1024视频免费在线观看| 亚洲人成77777在线视频| 两个人免费观看高清视频| 国产深夜福利视频在线观看| 一二三四在线观看免费中文在| 亚洲国产欧美网| 国产精品三级大全| 搡女人真爽免费视频火全软件| 亚洲欧美一区二区三区黑人 | 欧美日韩视频精品一区| 免费观看在线日韩| 亚洲,欧美,日韩| 国产在线一区二区三区精| 亚洲五月色婷婷综合| 精品卡一卡二卡四卡免费| 久久影院123| 欧美国产精品一级二级三级| 亚洲天堂av无毛| 国产av精品麻豆| 久久亚洲国产成人精品v| 国产精品麻豆人妻色哟哟久久| 亚洲欧美日韩另类电影网站| 一区二区三区乱码不卡18| 亚洲欧洲精品一区二区精品久久久 | 亚洲av男天堂| 母亲3免费完整高清在线观看 | 久久久亚洲精品成人影院| 色94色欧美一区二区| 1024香蕉在线观看| 黑人巨大精品欧美一区二区蜜桃| 狂野欧美激情性bbbbbb| 卡戴珊不雅视频在线播放| 咕卡用的链子| 久久精品夜色国产| 欧美av亚洲av综合av国产av | 免费看不卡的av| 午夜福利乱码中文字幕| 久久精品国产a三级三级三级| 麻豆乱淫一区二区| 亚洲成人手机| 国产人伦9x9x在线观看 | 国产亚洲精品第一综合不卡| 国产一级毛片在线| 青春草国产在线视频| 欧美精品亚洲一区二区| 欧美国产精品va在线观看不卡| 亚洲婷婷狠狠爱综合网| 久久韩国三级中文字幕| 最近中文字幕高清免费大全6| 亚洲av电影在线观看一区二区三区| av天堂久久9| av一本久久久久| 国产精品人妻久久久影院| 少妇被粗大猛烈的视频| 国产极品天堂在线| 久久免费观看电影| 狠狠婷婷综合久久久久久88av| 国产免费现黄频在线看| 女人久久www免费人成看片| 国产在线视频一区二区| 18+在线观看网站| 午夜免费男女啪啪视频观看| 日韩欧美一区视频在线观看| 国产色婷婷99| 青春草视频在线免费观看| 国产xxxxx性猛交| 欧美97在线视频| 欧美av亚洲av综合av国产av | 一级黄片播放器| 亚洲国产精品成人久久小说| 国产精品熟女久久久久浪| 美国免费a级毛片| 精品久久久精品久久久| 亚洲国产欧美网| 亚洲国产最新在线播放| av电影中文网址| 亚洲精品国产av成人精品| 夜夜骑夜夜射夜夜干| 亚洲综合精品二区| 亚洲综合精品二区| 中文字幕精品免费在线观看视频| 我要看黄色一级片免费的| 高清欧美精品videossex| 亚洲av福利一区| 一本—道久久a久久精品蜜桃钙片| 色吧在线观看| 精品午夜福利在线看| 久久久国产一区二区| av网站在线播放免费| 日韩av免费高清视频| 欧美另类一区| 成人18禁高潮啪啪吃奶动态图| 亚洲av电影在线观看一区二区三区| 精品第一国产精品| 黑人欧美特级aaaaaa片| 久久久久久久精品精品| 久久99热这里只频精品6学生| 校园人妻丝袜中文字幕| 日日摸夜夜添夜夜爱| 在线观看免费高清a一片| 久久韩国三级中文字幕| av国产久精品久网站免费入址| 欧美日韩国产mv在线观看视频| 亚洲国产最新在线播放| 欧美 亚洲 国产 日韩一| 老司机亚洲免费影院| 男人爽女人下面视频在线观看| 中文乱码字字幕精品一区二区三区| 在现免费观看毛片| 曰老女人黄片| 日韩中文字幕欧美一区二区 | 久久免费观看电影| 精品国产一区二区三区久久久樱花| 性高湖久久久久久久久免费观看| 在线天堂中文资源库| 岛国毛片在线播放| 看免费成人av毛片| 久久久久精品人妻al黑| 老汉色av国产亚洲站长工具| 日本黄色日本黄色录像| 免费大片黄手机在线观看| 伦精品一区二区三区| 国产一区二区三区综合在线观看| 亚洲久久久国产精品| 久久鲁丝午夜福利片| 日韩,欧美,国产一区二区三区| 看免费成人av毛片| 看十八女毛片水多多多| 熟女电影av网| 纵有疾风起免费观看全集完整版| 日韩av在线免费看完整版不卡| 中文欧美无线码| 免费在线观看黄色视频的| 亚洲精品乱久久久久久| 国产亚洲午夜精品一区二区久久| 国产亚洲欧美精品永久| 中文精品一卡2卡3卡4更新| 午夜久久久在线观看| 在线观看国产h片| 精品少妇久久久久久888优播| 色播在线永久视频| 男人爽女人下面视频在线观看| 日韩中文字幕欧美一区二区 | 国产成人精品无人区| 汤姆久久久久久久影院中文字幕| av免费观看日本| 极品人妻少妇av视频| 成年女人毛片免费观看观看9 | 国产高清不卡午夜福利| 你懂的网址亚洲精品在线观看| 男女免费视频国产| 男的添女的下面高潮视频| 这个男人来自地球电影免费观看 | 最近中文字幕2019免费版| 肉色欧美久久久久久久蜜桃| 亚洲av欧美aⅴ国产| 免费高清在线观看视频在线观看| 街头女战士在线观看网站| 精品国产乱码久久久久久小说| 国产成人91sexporn| 国产一区亚洲一区在线观看| 2018国产大陆天天弄谢| 男女边吃奶边做爰视频| 成人影院久久| 日本av手机在线免费观看| 91精品国产国语对白视频| 精品卡一卡二卡四卡免费| 欧美人与性动交α欧美软件| 美女国产高潮福利片在线看| 亚洲精品成人av观看孕妇| 蜜桃国产av成人99| 精品人妻一区二区三区麻豆| 一级爰片在线观看| 蜜桃在线观看..| 精品酒店卫生间| 纵有疾风起免费观看全集完整版| 两性夫妻黄色片| av免费在线看不卡| 有码 亚洲区| 成人影院久久| 观看美女的网站| 亚洲精品国产av成人精品| 欧美精品亚洲一区二区| 黄色 视频免费看| 在线观看人妻少妇| 热99久久久久精品小说推荐| 九草在线视频观看| 制服丝袜香蕉在线| 男人操女人黄网站| 日韩一卡2卡3卡4卡2021年| av一本久久久久| 欧美日韩国产mv在线观看视频| 一二三四在线观看免费中文在| 自线自在国产av| 黄片小视频在线播放| 国产精品av久久久久免费| 国产极品天堂在线| 成人国产av品久久久| 好男人视频免费观看在线| 亚洲一区中文字幕在线| 亚洲欧美日韩另类电影网站| 中文字幕亚洲精品专区| 亚洲第一av免费看| a级片在线免费高清观看视频| 日韩熟女老妇一区二区性免费视频| 午夜福利在线观看免费完整高清在| 1024视频免费在线观看| 日韩不卡一区二区三区视频在线| 免费久久久久久久精品成人欧美视频| 日本爱情动作片www.在线观看| 国产免费一区二区三区四区乱码| 亚洲一区中文字幕在线| 国产成人欧美| 在线观看一区二区三区激情| 99re6热这里在线精品视频| 叶爱在线成人免费视频播放| 亚洲av中文av极速乱| 久久久精品国产亚洲av高清涩受| 一边亲一边摸免费视频| 我的亚洲天堂| av不卡在线播放| 黄色 视频免费看| 亚洲国产看品久久| 美女高潮到喷水免费观看| 亚洲av日韩在线播放| 免费观看a级毛片全部| 中文精品一卡2卡3卡4更新| 亚洲 欧美一区二区三区| 在线 av 中文字幕| 狠狠精品人妻久久久久久综合| 久久精品熟女亚洲av麻豆精品| 午夜福利,免费看| 成年美女黄网站色视频大全免费| 中文欧美无线码| av网站在线播放免费| 卡戴珊不雅视频在线播放| 国产精品免费视频内射| 亚洲国产精品国产精品| 欧美黄色片欧美黄色片| 久久精品国产亚洲av高清一级| 最黄视频免费看| 久久久久人妻精品一区果冻| 中文字幕精品免费在线观看视频| 女的被弄到高潮叫床怎么办| 成人午夜精彩视频在线观看| 国产成人精品在线电影| 久久精品亚洲av国产电影网| 午夜激情久久久久久久| 丝袜美足系列| 亚洲国产欧美日韩在线播放| 赤兔流量卡办理| 亚洲av日韩在线播放| 日日撸夜夜添| 美女福利国产在线| 黄网站色视频无遮挡免费观看| h视频一区二区三区| 国产免费福利视频在线观看| 久久久久久伊人网av| 亚洲成人手机| 999久久久国产精品视频| 尾随美女入室| 亚洲精品日韩在线中文字幕| 可以免费在线观看a视频的电影网站 | 天堂俺去俺来也www色官网| 妹子高潮喷水视频| 久久久精品免费免费高清| 中文字幕人妻丝袜制服| 一本—道久久a久久精品蜜桃钙片| 日韩中文字幕欧美一区二区 | 久久久久久久亚洲中文字幕| 国产不卡av网站在线观看| 久久久久久久久久人人人人人人| 成年女人毛片免费观看观看9 | h视频一区二区三区| 成人毛片a级毛片在线播放| 黄片播放在线免费| 免费黄频网站在线观看国产| 欧美人与性动交α欧美精品济南到 | 丰满少妇做爰视频| 99国产精品免费福利视频| 成年女人毛片免费观看观看9 | 男女午夜视频在线观看| 免费人妻精品一区二区三区视频| 欧美国产精品va在线观看不卡| 亚洲久久久国产精品| 观看av在线不卡| 久久精品国产综合久久久| 午夜免费观看性视频| 一区在线观看完整版| 欧美人与性动交α欧美软件| 国产精品不卡视频一区二区| 最新中文字幕久久久久| av女优亚洲男人天堂| 黄色一级大片看看| 99热网站在线观看| 日韩一区二区三区影片| 香蕉丝袜av| 国产亚洲最大av| 国产色婷婷99| 欧美中文综合在线视频| 久久久欧美国产精品| 成人18禁高潮啪啪吃奶动态图| 久久毛片免费看一区二区三区| 国产精品 欧美亚洲| 日韩一区二区视频免费看| 国产精品免费大片| 久久精品aⅴ一区二区三区四区 | 精品亚洲成国产av| 十分钟在线观看高清视频www| 国产野战对白在线观看| 大码成人一级视频| 国产精品香港三级国产av潘金莲 | 搡女人真爽免费视频火全软件| 熟妇人妻不卡中文字幕| 一本—道久久a久久精品蜜桃钙片| 老女人水多毛片| 亚洲成国产人片在线观看| 国产亚洲av片在线观看秒播厂| 十分钟在线观看高清视频www| 老女人水多毛片| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 丰满迷人的少妇在线观看| 亚洲精品一二三| 国产亚洲午夜精品一区二区久久| 少妇被粗大的猛进出69影院| 大码成人一级视频| 久久狼人影院| 伦理电影免费视频| 欧美精品高潮呻吟av久久| 亚洲成国产人片在线观看| 秋霞在线观看毛片| 永久免费av网站大全| 9191精品国产免费久久| 国产有黄有色有爽视频| 亚洲一区中文字幕在线| 国产男女超爽视频在线观看| 黄色毛片三级朝国网站| 久久久久国产网址| 午夜91福利影院| 国产男女超爽视频在线观看| 亚洲欧美清纯卡通| 交换朋友夫妻互换小说| 最近2019中文字幕mv第一页| 街头女战士在线观看网站| 国产无遮挡羞羞视频在线观看| 中文欧美无线码| 蜜桃国产av成人99| 国产精品秋霞免费鲁丝片| 久久精品国产鲁丝片午夜精品| 99久久中文字幕三级久久日本| 国产精品麻豆人妻色哟哟久久| 中文字幕色久视频| 久热这里只有精品99| 中国三级夫妇交换| 精品国产一区二区三区四区第35| 爱豆传媒免费全集在线观看| 精品人妻在线不人妻| 久久综合国产亚洲精品| 在线精品无人区一区二区三| 亚洲成人手机| 久久久久国产网址| 国产免费一区二区三区四区乱码| 久久精品人人爽人人爽视色| 99久久综合免费| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 久久免费观看电影| 建设人人有责人人尽责人人享有的| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩另类电影网站| 久久99热这里只频精品6学生| 国产又色又爽无遮挡免| 日日摸夜夜添夜夜爱| 免费高清在线观看日韩| 亚洲欧美中文字幕日韩二区| 视频在线观看一区二区三区| 如何舔出高潮| 99久久中文字幕三级久久日本| 精品一区二区三卡| 国产成人av激情在线播放| 美女高潮到喷水免费观看| 一本久久精品| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 纵有疾风起免费观看全集完整版| 激情视频va一区二区三区| 看非洲黑人一级黄片| 女人高潮潮喷娇喘18禁视频| 2018国产大陆天天弄谢| 99re6热这里在线精品视频| 男男h啪啪无遮挡| 男女免费视频国产| 一区二区三区激情视频| 亚洲欧美成人综合另类久久久| 一级爰片在线观看| 国产视频首页在线观看| 中文字幕人妻丝袜一区二区 | 韩国精品一区二区三区| 少妇的丰满在线观看| 蜜桃在线观看..| 亚洲精品av麻豆狂野| 午夜福利在线观看免费完整高清在| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 国产亚洲一区二区精品| 午夜av观看不卡| 中文字幕制服av| 十分钟在线观看高清视频www| 考比视频在线观看| 亚洲av电影在线观看一区二区三区| 欧美+日韩+精品| 亚洲婷婷狠狠爱综合网| 欧美精品人与动牲交sv欧美| 久久久久久久大尺度免费视频| 国产97色在线日韩免费| 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 一区二区日韩欧美中文字幕| av天堂久久9| 高清视频免费观看一区二区| 日韩欧美精品免费久久| 一级,二级,三级黄色视频| 狠狠精品人妻久久久久久综合| 日日撸夜夜添| av在线老鸭窝| 两个人看的免费小视频| 日韩熟女老妇一区二区性免费视频| 亚洲四区av| 香蕉国产在线看| 久久精品国产a三级三级三级| 精品午夜福利在线看| 免费看不卡的av| 亚洲精品久久午夜乱码| 久久人妻熟女aⅴ| 国产乱来视频区| 大香蕉久久成人网| 久久精品国产鲁丝片午夜精品| 黑人猛操日本美女一级片| 国产欧美亚洲国产| 夜夜骑夜夜射夜夜干| 这个男人来自地球电影免费观看 | 夫妻性生交免费视频一级片| 国产 精品1| 纵有疾风起免费观看全集完整版| 一边摸一边做爽爽视频免费| 日产精品乱码卡一卡2卡三| 久久久久久久国产电影| 国产麻豆69| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区在线观看av| 午夜影院在线不卡| 天堂8中文在线网| 日本-黄色视频高清免费观看| 老司机影院成人| 中文欧美无线码| 9热在线视频观看99| 美女主播在线视频| 欧美在线黄色| 久久 成人 亚洲| 中文字幕人妻丝袜制服| 99九九在线精品视频| 日本-黄色视频高清免费观看| 国产免费一区二区三区四区乱码| 中文字幕人妻丝袜制服| 韩国av在线不卡| 亚洲精品视频女| 免费看av在线观看网站| 一本色道久久久久久精品综合| 18禁国产床啪视频网站| 午夜91福利影院| 在线观看免费日韩欧美大片| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 考比视频在线观看| 99久久综合免费| 视频在线观看一区二区三区| 成年av动漫网址| 国产精品 欧美亚洲| 91精品三级在线观看| 亚洲 欧美一区二区三区| 欧美xxⅹ黑人| 欧美日韩av久久| 欧美日韩一级在线毛片| 午夜免费观看性视频| 王馨瑶露胸无遮挡在线观看| 一级,二级,三级黄色视频| 黑人巨大精品欧美一区二区蜜桃| 男女边吃奶边做爰视频| 日韩在线高清观看一区二区三区| 好男人视频免费观看在线| 母亲3免费完整高清在线观看 | 亚洲av日韩在线播放| 老汉色∧v一级毛片| 亚洲av福利一区| 老汉色av国产亚洲站长工具| 成年人午夜在线观看视频| 看非洲黑人一级黄片| 一区二区日韩欧美中文字幕| 久久韩国三级中文字幕| 1024视频免费在线观看| 新久久久久国产一级毛片| 久久午夜综合久久蜜桃| av.在线天堂| 久久精品夜色国产| 久久久精品国产亚洲av高清涩受| 亚洲视频免费观看视频| 2021少妇久久久久久久久久久| 免费人妻精品一区二区三区视频| 一个人免费看片子| 中文字幕最新亚洲高清| 国产精品av久久久久免费| 亚洲国产欧美日韩在线播放| 大片电影免费在线观看免费| 日韩,欧美,国产一区二区三区| 中文乱码字字幕精品一区二区三区| 欧美日韩视频精品一区| 欧美在线黄色| 男男h啪啪无遮挡| 青草久久国产| 丁香六月天网| 久久影院123| 大话2 男鬼变身卡| 在线亚洲精品国产二区图片欧美| 国产人伦9x9x在线观看 | 日韩中文字幕欧美一区二区 | 欧美激情极品国产一区二区三区| 99久久精品国产国产毛片| 免费人妻精品一区二区三区视频| 老女人水多毛片| 国产av国产精品国产| 亚洲av国产av综合av卡|