• <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)信息
    99热这里只有是精品50| 动漫黄色视频在线观看| 日本黄大片高清| 亚洲一区二区三区色噜噜| 国产aⅴ精品一区二区三区波| 精品一区二区三区四区五区乱码| 中出人妻视频一区二区| 成年免费大片在线观看| 国产人伦9x9x在线观看| 精品国内亚洲2022精品成人| 欧美一级毛片孕妇| 国产精品久久电影中文字幕| 午夜精品一区二区三区免费看| 俺也久久电影网| 好男人在线观看高清免费视频| 欧美成狂野欧美在线观看| 日本 欧美在线| 欧美色视频一区免费| 日本撒尿小便嘘嘘汇集6| 色av中文字幕| 亚洲专区字幕在线| 妹子高潮喷水视频| 天天躁夜夜躁狠狠躁躁| 国产男靠女视频免费网站| 美女午夜性视频免费| 国产精品久久久久久精品电影| 久久这里只有精品中国| 久久香蕉精品热| 国产三级在线视频| 午夜福利18| 麻豆国产97在线/欧美 | 免费看a级黄色片| 国产三级在线视频| avwww免费| 五月玫瑰六月丁香| 亚洲aⅴ乱码一区二区在线播放 | 国产精品乱码一区二三区的特点| 国产亚洲av高清不卡| 精品国产亚洲在线| 中文字幕人成人乱码亚洲影| 天堂动漫精品| 欧美成人性av电影在线观看| 亚洲成人久久爱视频| 特大巨黑吊av在线直播| 精品国内亚洲2022精品成人| 黄色丝袜av网址大全| 看片在线看免费视频| 精品福利观看| 男人舔女人的私密视频| 久久中文看片网| 97人妻精品一区二区三区麻豆| 在线播放国产精品三级| 国产99久久九九免费精品| 成人18禁高潮啪啪吃奶动态图| 在线国产一区二区在线| 动漫黄色视频在线观看| 美女免费视频网站| 亚洲熟女毛片儿| 久久精品aⅴ一区二区三区四区| 琪琪午夜伦伦电影理论片6080| 国产精品av视频在线免费观看| 久久精品综合一区二区三区| 在线观看午夜福利视频| 18禁裸乳无遮挡免费网站照片| 叶爱在线成人免费视频播放| 丝袜美腿诱惑在线| 亚洲五月天丁香| 99国产综合亚洲精品| av有码第一页| 日本一二三区视频观看| 91麻豆av在线| 成人欧美大片| 老司机午夜十八禁免费视频| 中文在线观看免费www的网站 | 三级男女做爰猛烈吃奶摸视频| 亚洲中文字幕日韩| 午夜福利视频1000在线观看| 最近在线观看免费完整版| 女人被狂操c到高潮| 两个人看的免费小视频| 成人永久免费在线观看视频| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美精品综合一区二区三区| 国产av在哪里看| 精品乱码久久久久久99久播| 18禁美女被吸乳视频| 欧美日韩一级在线毛片| 欧美久久黑人一区二区| 精品国内亚洲2022精品成人| 免费无遮挡裸体视频| 动漫黄色视频在线观看| 国产高清激情床上av| 在线免费观看的www视频| 国产私拍福利视频在线观看| a级毛片在线看网站| 久久人妻福利社区极品人妻图片| 好男人在线观看高清免费视频| 99久久综合精品五月天人人| 亚洲国产日韩欧美精品在线观看 | 我的老师免费观看完整版| 亚洲国产看品久久| 亚洲美女视频黄频| 国产在线观看jvid| 夜夜爽天天搞| 精品午夜福利视频在线观看一区| 国产亚洲精品av在线| 精品国产亚洲在线| 亚洲av美国av| 国产真实乱freesex| 法律面前人人平等表现在哪些方面| 国产日本99.免费观看| 国产成人av激情在线播放| 欧美激情久久久久久爽电影| 欧美av亚洲av综合av国产av| 亚洲中文av在线| 日日干狠狠操夜夜爽| 99精品久久久久人妻精品| 在线播放国产精品三级| 免费无遮挡裸体视频| 在线观看一区二区三区| 在线观看美女被高潮喷水网站 | 一进一出抽搐gif免费好疼| 国产三级在线视频| 日本 欧美在线| 又爽又黄无遮挡网站| 美女扒开内裤让男人捅视频| 搡老熟女国产l中国老女人| 久久久久国产精品人妻aⅴ院| 亚洲中文字幕一区二区三区有码在线看 | 免费在线观看完整版高清| 免费观看人在逋| 日韩有码中文字幕| 窝窝影院91人妻| 50天的宝宝边吃奶边哭怎么回事| 久久久久久免费高清国产稀缺| 久热爱精品视频在线9| 久久婷婷人人爽人人干人人爱| 中文字幕熟女人妻在线| 欧美黑人巨大hd| 99久久综合精品五月天人人| 免费在线观看视频国产中文字幕亚洲| 国产高清视频在线播放一区| 色综合站精品国产| 99精品久久久久人妻精品| 亚洲天堂国产精品一区在线| 一二三四在线观看免费中文在| 神马国产精品三级电影在线观看 | 久久精品成人免费网站| 婷婷丁香在线五月| 欧美日韩精品网址| 亚洲成人中文字幕在线播放| АⅤ资源中文在线天堂| 在线观看免费日韩欧美大片| 美女扒开内裤让男人捅视频| 欧美一级a爱片免费观看看 | 高潮久久久久久久久久久不卡| 国产亚洲av高清不卡| 午夜福利在线在线| 国产精品爽爽va在线观看网站| 日本一二三区视频观看| 精品国内亚洲2022精品成人| 亚洲va日本ⅴa欧美va伊人久久| 精品国产亚洲在线| 99精品欧美一区二区三区四区| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 天堂av国产一区二区熟女人妻 | 在线观看一区二区三区| 黄色视频不卡| avwww免费| 免费电影在线观看免费观看| 给我免费播放毛片高清在线观看| 欧美乱码精品一区二区三区| 村上凉子中文字幕在线| 夜夜看夜夜爽夜夜摸| 人人妻,人人澡人人爽秒播| 婷婷亚洲欧美| 亚洲专区中文字幕在线| 国产一区二区激情短视频| 亚洲熟妇熟女久久| 久久香蕉激情| 国产精品 国内视频| 特大巨黑吊av在线直播| 亚洲精品一卡2卡三卡4卡5卡| 成人永久免费在线观看视频| 777久久人妻少妇嫩草av网站| 无遮挡黄片免费观看| 亚洲美女视频黄频| 日韩高清综合在线| 色综合欧美亚洲国产小说| 欧美绝顶高潮抽搐喷水| 最好的美女福利视频网| 黄片大片在线免费观看| 欧美久久黑人一区二区| 夜夜躁狠狠躁天天躁| 动漫黄色视频在线观看| 国产黄a三级三级三级人| 搡老岳熟女国产| 黄色视频,在线免费观看| 久99久视频精品免费| 人妻丰满熟妇av一区二区三区| 又黄又粗又硬又大视频| 日本黄色视频三级网站网址| 男女做爰动态图高潮gif福利片| 一个人观看的视频www高清免费观看 | 精品一区二区三区av网在线观看| 日韩中文字幕欧美一区二区| 国产午夜精品久久久久久| 国模一区二区三区四区视频 | a在线观看视频网站| 欧美国产日韩亚洲一区| 丝袜人妻中文字幕| 麻豆一二三区av精品| 老司机靠b影院| 国产91精品成人一区二区三区| 午夜福利在线观看吧| 青草久久国产| 啦啦啦韩国在线观看视频| 久久 成人 亚洲| 欧美乱码精品一区二区三区| 久久国产精品影院| 亚洲欧美日韩高清在线视频| 中文资源天堂在线| 婷婷亚洲欧美| 婷婷亚洲欧美| 麻豆久久精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 日韩 欧美 亚洲 中文字幕| 国产亚洲av高清不卡| 伊人久久大香线蕉亚洲五| 美女免费视频网站| 十八禁网站免费在线| 久久香蕉国产精品| 变态另类丝袜制服| 亚洲激情在线av| 欧美一区二区国产精品久久精品 | 丝袜人妻中文字幕| 在线观看午夜福利视频| 欧美中文日本在线观看视频| av片东京热男人的天堂| 精品高清国产在线一区| 两个人的视频大全免费| 精品电影一区二区在线| 国产精品久久视频播放| 国内久久婷婷六月综合欲色啪| 国产精品99久久99久久久不卡| 久久精品夜夜夜夜夜久久蜜豆 | av片东京热男人的天堂| 欧美绝顶高潮抽搐喷水| 国产又色又爽无遮挡免费看| 日本免费a在线| 色综合欧美亚洲国产小说| 午夜激情av网站| 欧美精品啪啪一区二区三区| 亚洲精品久久成人aⅴ小说| 女生性感内裤真人,穿戴方法视频| 久久午夜亚洲精品久久| 久久草成人影院| 国产精品一区二区免费欧美| 国产成人av教育| 国产在线精品亚洲第一网站| a级毛片在线看网站| 最近最新免费中文字幕在线| 男插女下体视频免费在线播放| 免费在线观看黄色视频的| 观看免费一级毛片| 午夜日韩欧美国产| 亚洲午夜理论影院| 在线观看免费视频日本深夜| 18禁观看日本| 国产精品免费一区二区三区在线| 麻豆久久精品国产亚洲av| 欧美 亚洲 国产 日韩一| 777久久人妻少妇嫩草av网站| 999久久久精品免费观看国产| 亚洲中文字幕一区二区三区有码在线看 | 国产av在哪里看| 亚洲,欧美精品.| 亚洲人成伊人成综合网2020| 久久久久精品国产欧美久久久| 国产精品一区二区三区四区免费观看 | 国产野战对白在线观看| 中文字幕av在线有码专区| 首页视频小说图片口味搜索| 制服人妻中文乱码| 人人妻人人看人人澡| 国产在线精品亚洲第一网站| 三级男女做爰猛烈吃奶摸视频| 免费看a级黄色片| 久久99热这里只有精品18| 日韩欧美精品v在线| 国产一区二区三区视频了| 欧美激情久久久久久爽电影| 亚洲人成网站高清观看| 十八禁网站免费在线| 欧美日本亚洲视频在线播放| 大型黄色视频在线免费观看| 黑人巨大精品欧美一区二区mp4| 老司机深夜福利视频在线观看| 欧美高清成人免费视频www| 午夜日韩欧美国产| 精品久久久久久久人妻蜜臀av| 欧美绝顶高潮抽搐喷水| 成人av一区二区三区在线看| 99国产精品一区二区三区| 午夜福利在线在线| 男女床上黄色一级片免费看| 中文字幕熟女人妻在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品中文字幕在线视频| 欧美日韩乱码在线| 日日干狠狠操夜夜爽| 久久久久性生活片| 精品无人区乱码1区二区| 亚洲国产精品久久男人天堂| 亚洲午夜理论影院| 亚洲精品av麻豆狂野| 成人永久免费在线观看视频| 欧美性长视频在线观看| 亚洲成人久久性| 国产精品久久久av美女十八| 一进一出好大好爽视频| 黄频高清免费视频| 国产一区二区在线av高清观看| 一级片免费观看大全| 国产高清激情床上av| 午夜福利在线在线| 在线免费观看的www视频| 嫁个100分男人电影在线观看| av片东京热男人的天堂| 免费无遮挡裸体视频| 波多野结衣高清作品| 97超级碰碰碰精品色视频在线观看| 亚洲欧美精品综合一区二区三区| 日本 av在线| 毛片女人毛片| 91老司机精品| 国产亚洲精品第一综合不卡| 俺也久久电影网| 窝窝影院91人妻| 国产一区二区在线观看日韩 | 国产精华一区二区三区| 黄色毛片三级朝国网站| 亚洲美女视频黄频| 操出白浆在线播放| 免费在线观看影片大全网站| 亚洲一码二码三码区别大吗| 亚洲人成网站高清观看| 18禁黄网站禁片免费观看直播| 午夜精品在线福利| 亚洲国产精品合色在线| 午夜福利免费观看在线| 级片在线观看| 久久久精品大字幕| 美女免费视频网站| 99久久国产精品久久久| 欧美乱码精品一区二区三区| 欧美黑人欧美精品刺激| 18禁黄网站禁片午夜丰满| 亚洲一区中文字幕在线| av中文乱码字幕在线| 国产成人精品久久二区二区91| 亚洲在线自拍视频| 久久精品aⅴ一区二区三区四区| 一级作爱视频免费观看| 又黄又爽又免费观看的视频| 国产在线观看jvid| 好男人在线观看高清免费视频| 人妻久久中文字幕网| 精品国产亚洲在线| 久久人妻福利社区极品人妻图片| 亚洲av成人精品一区久久| av天堂在线播放| 免费在线观看成人毛片| 国产午夜福利久久久久久| 日韩 欧美 亚洲 中文字幕| 久久久国产成人免费| 免费在线观看黄色视频的| 丰满人妻熟妇乱又伦精品不卡| 午夜福利高清视频| 久久精品影院6| 蜜桃久久精品国产亚洲av| 亚洲中文字幕一区二区三区有码在线看 | 色综合婷婷激情| 色综合亚洲欧美另类图片| xxxwww97欧美| 国产av一区在线观看免费| 亚洲美女黄片视频| 两性夫妻黄色片| 无限看片的www在线观看| www日本黄色视频网| 成年女人毛片免费观看观看9| 我的老师免费观看完整版| 国产精品98久久久久久宅男小说| 床上黄色一级片| 欧美中文综合在线视频| 久久精品91蜜桃| 亚洲成av人片免费观看| 少妇熟女aⅴ在线视频| 丁香欧美五月| 欧美黑人巨大hd| 精品久久久久久久久久久久久| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 香蕉国产在线看| 欧美绝顶高潮抽搐喷水| АⅤ资源中文在线天堂| 看片在线看免费视频| 国产91精品成人一区二区三区| aaaaa片日本免费| 神马国产精品三级电影在线观看 | 不卡av一区二区三区| 免费看美女性在线毛片视频| 一区二区三区国产精品乱码| 欧美午夜高清在线| 亚洲一区高清亚洲精品| 亚洲国产看品久久| 97超级碰碰碰精品色视频在线观看| 日本成人三级电影网站| 老鸭窝网址在线观看| 国产成+人综合+亚洲专区| 变态另类丝袜制服| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利18| 后天国语完整版免费观看| 国产精品日韩av在线免费观看| 色av中文字幕| 国产成人系列免费观看| 欧美一级毛片孕妇| 日韩成人在线观看一区二区三区| 老司机福利观看| 在线观看日韩欧美| 欧美+亚洲+日韩+国产| 亚洲av片天天在线观看| 久久中文字幕一级| 一本大道久久a久久精品| 黄频高清免费视频| 三级男女做爰猛烈吃奶摸视频| 嫁个100分男人电影在线观看| 日本黄大片高清| aaaaa片日本免费| 亚洲精品av麻豆狂野| 日韩成人在线观看一区二区三区| 999精品在线视频| 在线观看www视频免费| 欧美黑人巨大hd| 婷婷丁香在线五月| 老熟妇乱子伦视频在线观看| 九九热线精品视视频播放| 日本三级黄在线观看| 两个人的视频大全免费| a在线观看视频网站| 特级一级黄色大片| 欧美在线一区亚洲| netflix在线观看网站| 一二三四社区在线视频社区8| 久久精品人妻少妇| 色av中文字幕| 国产av在哪里看| 亚洲欧美激情综合另类| 亚洲精品国产精品久久久不卡| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 国产精华一区二区三区| 国产精品一区二区精品视频观看| 成人永久免费在线观看视频| 久久精品综合一区二区三区| 亚洲天堂国产精品一区在线| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久精品国产欧美久久久| www日本黄色视频网| 亚洲成a人片在线一区二区| 欧美成人午夜精品| 免费观看人在逋| 亚洲午夜理论影院| 首页视频小说图片口味搜索| 正在播放国产对白刺激| 国产一区在线观看成人免费| 欧美成狂野欧美在线观看| 亚洲最大成人中文| netflix在线观看网站| 婷婷精品国产亚洲av在线| 国产精品美女特级片免费视频播放器 | 国产一级毛片七仙女欲春2| 1024手机看黄色片| 日韩欧美精品v在线| 久99久视频精品免费| 日韩欧美精品v在线| 免费看a级黄色片| x7x7x7水蜜桃| 嫩草影视91久久| 午夜精品久久久久久毛片777| 可以在线观看毛片的网站| svipshipincom国产片| 啦啦啦免费观看视频1| 中出人妻视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 91大片在线观看| 国产91精品成人一区二区三区| 成人欧美大片| 色综合站精品国产| 好男人在线观看高清免费视频| 国产午夜福利久久久久久| 丰满人妻一区二区三区视频av | 欧美极品一区二区三区四区| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 亚洲国产精品久久男人天堂| 中文字幕人成人乱码亚洲影| 老司机午夜十八禁免费视频| 亚洲 欧美 日韩 在线 免费| 最近最新免费中文字幕在线| 国语自产精品视频在线第100页| 搡老妇女老女人老熟妇| www日本在线高清视频| 狠狠狠狠99中文字幕| 久久久国产成人精品二区| 亚洲av成人不卡在线观看播放网| 欧美日韩福利视频一区二区| 午夜福利18| 国产黄色小视频在线观看| 一级黄色大片毛片| 日本a在线网址| 中文在线观看免费www的网站 | 亚洲aⅴ乱码一区二区在线播放 | 在线国产一区二区在线| 嫩草影院精品99| 国产成人精品久久二区二区免费| 国产一区二区三区在线臀色熟女| 亚洲中文日韩欧美视频| 嫩草影院精品99| 日本a在线网址| 91九色精品人成在线观看| 天天添夜夜摸| 国产成+人综合+亚洲专区| 欧美日韩乱码在线| 欧美 亚洲 国产 日韩一| 99国产综合亚洲精品| 亚洲av成人一区二区三| 国产精品爽爽va在线观看网站| 色综合婷婷激情| 少妇被粗大的猛进出69影院| 国内精品久久久久久久电影| 看免费av毛片| 他把我摸到了高潮在线观看| 精品一区二区三区四区五区乱码| 久久久久国产精品人妻aⅴ院| 久久亚洲真实| 成年版毛片免费区| 久久久久国内视频| 波多野结衣巨乳人妻| 国产一区二区在线av高清观看| 精品国内亚洲2022精品成人| 天堂动漫精品| 欧美日韩亚洲国产一区二区在线观看| 美女高潮喷水抽搐中文字幕| 首页视频小说图片口味搜索| 国产不卡一卡二| av片东京热男人的天堂| 欧美三级亚洲精品| 色精品久久人妻99蜜桃| 观看免费一级毛片| 变态另类丝袜制服| 国产人伦9x9x在线观看| 波多野结衣高清作品| 国产片内射在线| 男女床上黄色一级片免费看| 精品欧美一区二区三区在线| 国产精品香港三级国产av潘金莲| 免费在线观看成人毛片| av国产免费在线观看| av片东京热男人的天堂| 久久国产精品影院| 日韩高清综合在线| 欧美绝顶高潮抽搐喷水| 淫秽高清视频在线观看| 很黄的视频免费| 身体一侧抽搐| 午夜亚洲福利在线播放| 精品福利观看| а√天堂www在线а√下载| 最近在线观看免费完整版| 亚洲人成网站在线播放欧美日韩| 69av精品久久久久久| 成人精品一区二区免费| 国产一区二区三区视频了| 黄色 视频免费看| 免费在线观看影片大全网站| 国产午夜精品久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 男人舔奶头视频| 欧美zozozo另类| 中文字幕av在线有码专区| 国产精品av视频在线免费观看| 无人区码免费观看不卡| 国产精品久久久久久久电影 | 桃红色精品国产亚洲av| 亚洲第一欧美日韩一区二区三区| 黄色 视频免费看| 国产精品久久视频播放| 国产熟女午夜一区二区三区| 国产又黄又爽又无遮挡在线| 国产三级在线视频| 欧美乱妇无乱码| 亚洲欧美日韩东京热| 中国美女看黄片| 妹子高潮喷水视频| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女| 男女视频在线观看网站免费 | 黑人欧美特级aaaaaa片| 午夜免费成人在线视频| 欧美3d第一页| 一a级毛片在线观看| 黄色a级毛片大全视频|