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

    基于蒙特卡洛法和分形理論的接觸導(dǎo)熱計(jì)算

    2019-08-17 07:39:36馮妍卉王智學(xué)
    關(guān)鍵詞:熱導(dǎo)率熱阻維數(shù)

    張 琳,馮妍卉,王智學(xué)

    (1.中國(guó)航發(fā)四川燃?xì)鉁u輪研究院, 成都 610500;2.北京科技大學(xué) 能源與環(huán)境工程學(xué)院, 北京 100083)

    隨著傳熱技術(shù)的不斷發(fā)展和廣泛應(yīng)用,接觸導(dǎo)熱逐漸成為一個(gè)不可忽略的影響傳熱的因素,成為航空航天、機(jī)械制造、電子、低溫超導(dǎo)等學(xué)科和動(dòng)力工程領(lǐng)域中的研究熱點(diǎn)。接觸導(dǎo)熱是一個(gè)受材料熱物性、機(jī)械特性、表面形貌及接觸環(huán)境(壓力、溫度和間隙填充介質(zhì)等)等眾多因素影響的非線性問(wèn)題。在工程應(yīng)用中,通常需要降低或提高接觸熱阻以控制接觸傳熱,在設(shè)計(jì)過(guò)程中需要采用理論、實(shí)驗(yàn)或數(shù)值分析的方法對(duì)接觸熱導(dǎo)率或接觸熱阻進(jìn)行預(yù)測(cè)。但在不同的應(yīng)用中,材料種類繁多,應(yīng)用環(huán)境也不同,用實(shí)驗(yàn)手段來(lái)獲取接觸熱導(dǎo)率通常需要設(shè)計(jì)復(fù)雜的實(shí)驗(yàn)裝置,具有較高的成本。近年來(lái),由于計(jì)算機(jī)技術(shù)的發(fā)展,數(shù)值分析的方法越來(lái)越多地應(yīng)用于接觸導(dǎo)熱的研究中[1-3],但此類方法通常需要進(jìn)行大量的前、后數(shù)據(jù)處理工作,時(shí)間成本較高。因此,工程領(lǐng)域通常采用理論模型和相應(yīng)的計(jì)算公式對(duì)接觸熱導(dǎo)率進(jìn)行預(yù)測(cè)。國(guó)內(nèi)外學(xué)者對(duì)接觸導(dǎo)熱進(jìn)行了大量的研究,產(chǎn)生了許多具有重要參考價(jià)值和工程應(yīng)用價(jià)值的理論和方法。接觸導(dǎo)熱的理論模型和計(jì)算方法的研究包括對(duì)粗糙表面的形貌分析、接觸形變分析和接觸導(dǎo)熱分析,研究者將不同的粗糙表面形貌模型、形變模型和導(dǎo)熱模型進(jìn)行組合,形成了不同的接觸導(dǎo)熱計(jì)算模型。王安良等[4]對(duì)接觸導(dǎo)熱的理論預(yù)測(cè)方法進(jìn)行了較為全面的綜述。

    早期的接觸導(dǎo)熱計(jì)算模型主要以統(tǒng)計(jì)學(xué)參數(shù)描述固體粗糙表面形貌,認(rèn)為固體粗糙表面是由無(wú)數(shù)個(gè)均勻分布的微凸體堆積而成,這些微凸體的高度符合高斯分布。以此為基礎(chǔ),許多學(xué)者建立了接觸導(dǎo)熱的計(jì)算公式,如Mikic[5]、Yovanovich等[6]、陳劍楠等[7]。此外,Kumar等[8]建立了接觸熱阻的蒙特卡洛模型,該模型根據(jù)粗糙表面粗糙峰分布符合高斯分布的特性,在粗糙表面上隨機(jī)模擬N個(gè)符合高斯分布的粗糙峰,并在不同的接觸界面間距下計(jì)算出發(fā)生接觸的粗糙峰數(shù)目,每個(gè)發(fā)生接觸的粗糙峰都形成一個(gè)單點(diǎn)接觸熱阻,界面間總的接觸熱阻即為這些單點(diǎn)接觸熱阻的并聯(lián)。

    需要指出的是,表面粗糙高度的分布并不完全符合高斯分布,以上所述接觸導(dǎo)熱計(jì)算模型都建立在以統(tǒng)計(jì)學(xué)參數(shù)描述的表面形貌模型的基礎(chǔ)上,而研究表明,表面形貌的統(tǒng)計(jì)學(xué)參數(shù)很明顯地受到測(cè)量?jī)x器的分辨率和取樣長(zhǎng)度的影響[9],不能唯一確定地表征一個(gè)粗糙表面,在此基礎(chǔ)上建立的接觸導(dǎo)熱計(jì)算模型通常具有較大的不確定性,計(jì)算結(jié)果的準(zhǔn)確性與工程應(yīng)用需求差距較大。因此,研究者們開(kāi)始尋求與儀器分辨率和取樣長(zhǎng)度無(wú)關(guān)的粗糙表面表征參數(shù),而分形理論使得粗糙表面的唯一表征成為可能。分形理論采用反映表面不規(guī)則程度的分形維數(shù)D和反映表面輪廓幅值的比例參數(shù)G來(lái)確定一個(gè)表面輪廓,分形參數(shù)D和G都是與儀器分辨率和取樣長(zhǎng)度無(wú)關(guān)的參數(shù)。Majumdar等[9]采用修正的Weierstrass-Mandelbrot函數(shù)(W-M函數(shù))描述粗糙表面輪廓,建立了固體接觸界面的彈塑性分形接觸模型和接觸導(dǎo)熱的M-T分形模型。Warren等[10]采用Cantor集分形形貌模型模擬粗糙表面,建立了彈性-塑性接觸和完全塑性接觸的Cantor集接觸模型和接觸導(dǎo)熱模型。此外,趙蘭萍[11]、徐瑞萍[12]、Zou等[13]、Ji等[14]、馬麗娜[15]和李小彭等[16]也進(jìn)行了基于表面分形理論的接觸導(dǎo)熱計(jì)算模型研究。

    一系列研究表明,基于表面分形理論的接觸導(dǎo)熱計(jì)算模型,不受儀器分辨率和取樣長(zhǎng)度的影響,對(duì)接觸導(dǎo)熱的預(yù)測(cè)具有確定性,能夠較為準(zhǔn)確地預(yù)測(cè)接觸導(dǎo)熱。近年來(lái)將分形理論應(yīng)用于接觸導(dǎo)熱預(yù)測(cè)的研究較少,現(xiàn)有的分形理論和模型較為復(fù)雜,與實(shí)際工程應(yīng)用存在很大差距[17]。因此,本文提出了一種適合于工程應(yīng)用的接觸導(dǎo)熱分形計(jì)算模型。模型基于固體粗糙表面形貌的隨機(jī)特性和分形特性,采用蒙特卡洛法和W-M分形函數(shù)模擬表面形貌和接觸,構(gòu)建了固體界面間接觸熱阻的并聯(lián)網(wǎng)絡(luò)模型,簡(jiǎn)化了接觸熱阻的計(jì)算過(guò)程,模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果基本一致,能夠較為準(zhǔn)確地預(yù)測(cè)固體界面間的接觸導(dǎo)熱,具有較高的工程應(yīng)用價(jià)值。

    1 固體粗糙表面形貌的分形表征

    固體粗糙表面形貌是一個(gè)非穩(wěn)定隨機(jī)過(guò)程,如圖1所示[9]。表面輪廓曲線具有隨機(jī)、多尺度和無(wú)序的特性以及統(tǒng)計(jì)自仿射和自相似的數(shù)學(xué)特性,分形幾何學(xué)中的W-M函數(shù)可以滿足表面輪廓曲線的所有數(shù)學(xué)特征,其表達(dá)式如下:

    (1)

    式中:D是分形維數(shù),它反映的是輪廓z(x)在所有空間尺度上的不規(guī)則性;G是決定輪廓高度幅值的比例系數(shù),決定了z(x)的具體尺寸;γn為空間頻率,曲線的最低頻率取決于取樣長(zhǎng)度L,并有γn1=1/L。對(duì)于任意取樣長(zhǎng)度,決定輪廓高度z(x)的是分形參數(shù)D和G。

    圖1 表面形貌的統(tǒng)計(jì)自仿射和自相似特性

    由式(1)可以得出

    z(γx)=γ(2-D)z(x)

    (2)

    若位置坐標(biāo)x放大γ倍,則輪廓高度放大為原來(lái)的γ(2-D)倍。粗糙表面的分形參數(shù)D和G是尺度獨(dú)立的,可以提供存在于分形表面上的所有尺度范圍內(nèi)的全部粗糙度信息[10]。

    圖2為某粗糙表面實(shí)測(cè)形貌與采用分形分析后得到分形參數(shù)并采用W-M分形函數(shù)模擬的形貌對(duì)比,可以看出二者表面形貌曲線波動(dòng)幅度和頻率具有較高的符合性,W-M函數(shù)對(duì)固體粗糙表面的形貌具有較好的模擬效果。

    圖2 W-M函數(shù)模擬形貌與真實(shí)表面形貌對(duì)比

    2 基于蒙特卡洛法和分形理論的接觸熱阻計(jì)算

    2.1 基本假設(shè)

    本文研究在一定壓力下相接觸的固體表面間的接觸熱阻,基本假設(shè)如下:

    1) 粗糙表面輪廓具有隨機(jī)性和分形特性,固體表面之間的接觸可等效為一個(gè)當(dāng)量粗糙表面和剛性光滑平面之間的接觸。

    2) 固體表面間的接觸和傳熱僅僅發(fā)生在離散分布的一系列大小不同的接觸點(diǎn)上,接觸點(diǎn)發(fā)生塑性變形,忽略間隙介質(zhì)的導(dǎo)熱和熱對(duì)流以及非接觸部分的輻射換熱。

    3) 表面間各接觸點(diǎn)的傳熱互不影響,每個(gè)接觸點(diǎn)相當(dāng)于一個(gè)熱流通道,并形成一個(gè)單點(diǎn)接觸熱阻,總的接觸熱阻是所有單點(diǎn)接觸熱阻的并聯(lián)。

    2.2 接觸界面形貌的模擬

    接觸界面形貌模擬的目的是確定粗糙表面接觸點(diǎn)的數(shù)量、尺寸及分布等參數(shù)。粗糙表面輪廓本質(zhì)上是一個(gè)隨機(jī)過(guò)程,因此本文將分形理論與蒙特卡洛法相結(jié)合,對(duì)表面粗糙峰的分布進(jìn)行模擬。

    首先將2個(gè)粗糙表面等效成一個(gè)當(dāng)量粗糙表面,將粗糙表面的接觸等效為當(dāng)量粗糙表面與剛性光滑平面的接觸,從而求出當(dāng)量表面分形參數(shù)。一條數(shù)字化輪廓曲線可以視為一個(gè)時(shí)間序列,在z(x)曲線上取時(shí)間間隔為Δτ的N個(gè)采樣點(diǎn),令τ=nΔτ,則該曲線的結(jié)構(gòu)函數(shù)為

    E(τ)=〈[z(x+nΔτ)-z(x)]2〉=

    (3)

    等價(jià)于

    E(τ)=CG2(D-1)τ4-2D

    (4)

    其中

    (5)

    式中函數(shù)Γ為Gamma函數(shù),當(dāng)1

    在雙對(duì)數(shù)坐標(biāo)上,結(jié)構(gòu)函數(shù)E(τ)與τ呈線性關(guān)系,即

    lgE(τ)=B+klgτ

    (6)

    式中:

    k=4-2D

    (7)

    B=lgCG2(D-1)

    (8)

    對(duì)于當(dāng)量粗糙表面,其結(jié)構(gòu)函數(shù)為2個(gè)粗糙表面結(jié)構(gòu)函數(shù)之和,即

    E(τ)=E1(τ)+E2(τ)

    (9)

    當(dāng)量硬度H、當(dāng)量彈性模量E和當(dāng)量熱導(dǎo)率k通過(guò)以下各式求得:

    H=min(H1,H2)

    (10)

    (11)

    (12)

    式中ν為材料的泊松比,下標(biāo)1和2分別代表相接觸的2個(gè)表面。

    研究表明,表面粗糙峰均方根高度σ與分形參數(shù)具有如下關(guān)系[9]:

    (13)

    表面粗糙峰高度具有隨機(jī)性,并符合W-M函數(shù),因此對(duì)于一個(gè)邊長(zhǎng)為L(zhǎng)的方形接觸表面,可隨機(jī)選取N個(gè)x值(0

    N=L2η

    (14)

    式中η為表面粗糙峰密度。根據(jù)Hsieh[20]的結(jié)果,表面粗糙峰密度

    (15)

    式中m為表面粗糙峰平均斜率。

    2.3 表面接觸參數(shù)的計(jì)算

    本文采用的接觸點(diǎn)形變模型為塑性圓錐模型,如圖3所示。假設(shè)當(dāng)量粗糙表面與剛性光滑平面的距離為d,則高度為zi的粗糙峰的峰頂與剛性光滑平面的距離δi為

    δi=zi-d

    (16)

    當(dāng)δi>0時(shí),粗糙峰與剛性光滑平面發(fā)生接觸,其接觸半徑ai、接觸面積Aci和接觸載荷Fi分別為:

    ai=δi/m

    (17)

    (18)

    (19)

    通過(guò)對(duì)所有粗糙峰的計(jì)算,可以得到發(fā)生接觸的粗糙峰的個(gè)數(shù)n,進(jìn)而計(jì)算出表面距離為d時(shí)的實(shí)際接觸面積Ar和總接觸壓力Fc:

    (20)

    (21)

    上述計(jì)算過(guò)程是在已知d的情況下求界面間接觸載荷Fc等參數(shù),然而在進(jìn)行實(shí)驗(yàn)研究或?qū)嶋H工程應(yīng)用中很難知道接觸界面的實(shí)際距離,但很容易知道接觸載荷Fc,因此在接觸熱阻的理論計(jì)算中可迭代求得d,進(jìn)而求出其他參數(shù)。

    圖3 塑性圓錐形變模型

    2.4 單點(diǎn)接觸熱阻和總熱阻的計(jì)算

    研究表明,半無(wú)限長(zhǎng)流管換熱模型更適用于接觸熱阻的計(jì)算,Yovanovich[6]采用的單點(diǎn)接觸熱阻計(jì)算公式為

    (22)

    式中

    ψ=(1-ε)1.5

    (23)

    并且有

    (24)

    式(24)中Aa為名義接觸面積,

    Aa=L2

    (25)

    根據(jù)本文的基本假設(shè)可知,總熱阻為各單點(diǎn)接觸熱阻的并聯(lián),因此可以得出總接觸熱阻Rc和接觸熱導(dǎo)率hc:

    (26)

    (27)

    3 結(jié)果與分析

    3.1 計(jì)算方法的驗(yàn)證

    本文采用文獻(xiàn)[11]中的2組實(shí)驗(yàn)數(shù)據(jù)對(duì)接觸熱阻的計(jì)算模型進(jìn)行了驗(yàn)證,并與Mikic模型、Yovanovich模型和MT分形模型進(jìn)行了對(duì)比。實(shí)驗(yàn)材料為鋁合金Al5052,實(shí)驗(yàn)溫度為155 K,表1和表2分別為材料Al5052的機(jī)械性能和熱力學(xué)性能參數(shù)以及4個(gè)Al5052試樣的表面分形參數(shù)。

    圖4和圖5分別為理論計(jì)算結(jié)果與文獻(xiàn)[11]中的實(shí)驗(yàn)結(jié)果的對(duì)比,可以看出:Yovanovich模型的接觸熱導(dǎo)率計(jì)算值高于實(shí)驗(yàn)值,并有很大偏差,而Mikic模型的接觸熱導(dǎo)率計(jì)算值小于實(shí)驗(yàn)值。本文的分形模型和MT分形模型的接觸熱導(dǎo)率計(jì)算值與實(shí)驗(yàn)值比較接近,這表明分形模型能夠較為準(zhǔn)確地預(yù)測(cè)接觸熱阻。在Yovanovich模型中,認(rèn)為粗糙表面微凸體的形變是塑形形變,但并沒(méi)有考慮變形后體積的損失,所以在計(jì)算過(guò)程中,每對(duì)微凸體接觸點(diǎn)的接觸半徑都要小于實(shí)際變形后的接觸半徑,接觸點(diǎn)的對(duì)數(shù)就會(huì)相應(yīng)地增加,因此得出的接觸熱導(dǎo)率計(jì)算值偏大。Mikic模型認(rèn)為微凸體的變形為彈性形變,而在實(shí)際接觸過(guò)程中,隨著壓力的增大,必然會(huì)引起塑形形變,所以接觸熱導(dǎo)率求解的結(jié)果會(huì)偏低。

    表1 鋁試樣Al5052的熱力學(xué)和機(jī)械性能

    表2 鋁試樣Al5052的表面分形參數(shù)

    圖4 接觸熱導(dǎo)率理論計(jì)算值與實(shí)驗(yàn)值對(duì)比(試樣1/試樣2)

    圖5 接觸熱導(dǎo)率理論計(jì)算值與實(shí)驗(yàn)值對(duì)比(試樣3/試樣4)

    表面形貌是接觸導(dǎo)熱的主要影響因素,傳統(tǒng)的粗糙度參數(shù)受到儀器分辨率和取樣長(zhǎng)度的影響,使用尺度獨(dú)立的分形參數(shù)模擬粗糙表面形貌建立的接觸熱阻計(jì)算模型更有利于對(duì)接觸傳熱的預(yù)測(cè)。

    3.2 分形參數(shù)對(duì)接觸熱阻的影響分析

    計(jì)算模擬了2個(gè)具有相同性能的尺寸為2.5 mm×2.5 mm的材料表面的接觸(材料參數(shù)見(jiàn)表3),計(jì)算中設(shè)定的分形參數(shù)為當(dāng)量分形參數(shù)。

    表3 模擬材料的熱力學(xué)和機(jī)械性能

    圖6為G=1.0×10-5m時(shí),不同壓力下的接觸熱導(dǎo)率與分形維數(shù)D的關(guān)系??梢钥闯觯簤毫υ酱蠼缑骈g的接觸熱導(dǎo)率越大,這反映了壓力增大使界面間的接觸面積增大。在定壓下看,接觸熱導(dǎo)率隨著分形維數(shù)D的增大而增大,這是因?yàn)楫?dāng)D增大時(shí),表面輪廓的精細(xì)結(jié)構(gòu)增多,精細(xì)結(jié)構(gòu)增加導(dǎo)致整個(gè)輪廓的整體高度降低,減低了輪廓粗糙度,導(dǎo)致接觸導(dǎo)熱的增強(qiáng),壓力越大接觸導(dǎo)熱的增強(qiáng)越明顯。

    圖7為在接觸壓力p=2 MPa時(shí)、不同分形參數(shù)G下,接觸熱導(dǎo)率與分形維數(shù)D的關(guān)系??梢钥闯觯航佑|熱導(dǎo)率隨分形維數(shù)D的增大而增大,分形比例參數(shù)G越小,增大越明顯。圖8為接觸壓力p=2 MPa時(shí)、不同分形維數(shù)D下接觸熱導(dǎo)率與分形參數(shù)G的關(guān)系,可以看出:分形維數(shù)D一定時(shí),分形參數(shù)G越大,接觸熱導(dǎo)率越小。這是因?yàn)榉中螀?shù)G越大,表面粗糙度越大,增加了界面間的傳熱阻力。分形維數(shù)D越大,G的變化對(duì)接觸熱導(dǎo)率的影響越明顯。

    圖6 不同接觸壓力下分形維數(shù)D對(duì)接觸熱導(dǎo)率的影響(G=1.0×10-5 m)

    圖7 不同分形比例參數(shù)G條件下分形維數(shù)D對(duì)接觸熱導(dǎo)率的影響(P=2 MPa)

    圖8 不同分形維數(shù)D的條件下分形比例參數(shù)G對(duì)接觸熱導(dǎo)率的影響(P=2 MPa)

    綜上所述,在一定的接觸壓力下,接觸界面間的熱導(dǎo)率隨著分形維數(shù)D的增大而增大,隨著分形參數(shù)G的增大而減小。這是因?yàn)樵谙嗤姆中螀?shù)G下,D越大表面粗糙度越小,在相同的分形參數(shù)D下,G越小表面粗糙度越小。表面粗糙度越小,材料表面越光滑,表面之間的接觸越充分,接觸熱導(dǎo)率越大。

    4 結(jié)束語(yǔ)

    本文采用蒙特卡洛法和W-M分形函數(shù)模擬表面形貌,建立了一種能夠反映固體粗糙表面形貌的隨機(jī)特性和分形特性的固體接觸界面間接觸導(dǎo)熱的計(jì)算模型,其計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果基本一致,能夠較為準(zhǔn)確地預(yù)測(cè)固體界面間的接觸導(dǎo)熱。分析了分形參數(shù)對(duì)固體界面間接觸導(dǎo)熱的影響。分析結(jié)果表明:在一定的接觸壓力下,分形參數(shù)G一定時(shí),接觸熱導(dǎo)率隨著分形維數(shù)D的增大而增大,G越小,增大越明顯;分形參數(shù)D一定時(shí),接觸熱導(dǎo)率隨著分形參數(shù)G的增大而減小,分形維數(shù)D越大,G的變化對(duì)接觸熱導(dǎo)率的影響越明顯。由于分形參數(shù)能夠唯一確定地表征固體粗糙表面的形貌,本文提出的模型對(duì)接觸熱阻的預(yù)測(cè)具有確定性,為接觸傳熱的工程設(shè)計(jì)提供了一種有效的方法。

    猜你喜歡
    熱導(dǎo)率熱阻維數(shù)
    豎直單U型地埋管換熱器埋管間負(fù)熱阻現(xiàn)象的參數(shù)化研究與分析
    β-變換中一致丟番圖逼近問(wèn)題的維數(shù)理論
    空位缺陷對(duì)單層石墨烯導(dǎo)熱特性影響的分子動(dòng)力學(xué)
    連續(xù)碳纖維鋁基復(fù)合材料橫向等效熱導(dǎo)率的模擬分析
    Si3N4/BN復(fù)合陶瓷熱導(dǎo)率及其有限元分析
    一類齊次Moran集的上盒維數(shù)
    界面熱阻對(duì)L型鎂合金鑄件凝固過(guò)程溫度場(chǎng)的影響
    關(guān)于齊次Moran集的packing維數(shù)結(jié)果
    換熱設(shè)備污垢熱阻和腐蝕監(jiān)測(cè)技術(shù)綜述
    涉及相變問(wèn)題Julia集的Hausdorff維數(shù)
    日日干狠狠操夜夜爽| 一进一出好大好爽视频| 久久中文看片网| 日韩欧美免费精品| 亚洲一区高清亚洲精品| 欧美乱色亚洲激情| 亚洲全国av大片| 日韩欧美三级三区| 久久久久九九精品影院| 国产在线精品亚洲第一网站| 亚洲av五月六月丁香网| 欧美在线黄色| 亚洲欧美一区二区三区黑人| 精品日产1卡2卡| 亚洲av成人一区二区三| 999久久久国产精品视频| 精品乱码久久久久久99久播| 久久午夜亚洲精品久久| 久久热在线av| 身体一侧抽搐| 国产精品二区激情视频| 给我免费播放毛片高清在线观看| 91大片在线观看| 97人妻天天添夜夜摸| 精品久久久久久,| 露出奶头的视频| 在线观看日韩欧美| 窝窝影院91人妻| 国产99白浆流出| 欧美午夜高清在线| 国产精华一区二区三区| 国产精品一区二区免费欧美| 亚洲精品中文字幕在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 一区二区三区精品91| 午夜久久久久精精品| 777久久人妻少妇嫩草av网站| 丁香欧美五月| 91麻豆精品激情在线观看国产| 久久久久久久午夜电影| 嫁个100分男人电影在线观看| 操出白浆在线播放| 亚洲成国产人片在线观看| 欧美中文综合在线视频| 热re99久久国产66热| 岛国在线观看网站| 久久精品亚洲熟妇少妇任你| 一边摸一边抽搐一进一出视频| 亚洲欧美激情在线| 久久久久久人人人人人| 在线视频色国产色| 一区二区日韩欧美中文字幕| 亚洲精华国产精华精| 久久精品国产99精品国产亚洲性色 | 校园春色视频在线观看| svipshipincom国产片| 成人三级做爰电影| 一级,二级,三级黄色视频| 久久婷婷人人爽人人干人人爱 | 成在线人永久免费视频| 欧美激情高清一区二区三区| 欧美日韩乱码在线| 91国产中文字幕| 日本撒尿小便嘘嘘汇集6| 99精品久久久久人妻精品| 女警被强在线播放| 欧美成狂野欧美在线观看| 亚洲av熟女| 亚洲国产精品成人综合色| 午夜精品在线福利| 极品人妻少妇av视频| 熟女少妇亚洲综合色aaa.| 夜夜夜夜夜久久久久| 美女 人体艺术 gogo| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久蜜臀av无| 久久午夜综合久久蜜桃| 久久 成人 亚洲| 国产精品1区2区在线观看.| 午夜福利高清视频| 9色porny在线观看| 黄片大片在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久精品吃奶| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线观看免费| www.999成人在线观看| tocl精华| 女人爽到高潮嗷嗷叫在线视频| 国产男靠女视频免费网站| 国产av精品麻豆| 亚洲免费av在线视频| 香蕉久久夜色| 国产精品久久久av美女十八| 日韩高清综合在线| 欧美最黄视频在线播放免费| 色在线成人网| 51午夜福利影视在线观看| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| 日本 av在线| 久久精品91无色码中文字幕| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 亚洲精品中文字幕一二三四区| 国产亚洲欧美精品永久| 国产一级毛片七仙女欲春2 | 99国产精品一区二区蜜桃av| 精品一区二区三区av网在线观看| 乱人伦中国视频| 大陆偷拍与自拍| 女人被狂操c到高潮| 成人国产综合亚洲| 91精品三级在线观看| av福利片在线| 国产激情欧美一区二区| 欧美不卡视频在线免费观看 | 9热在线视频观看99| 女同久久另类99精品国产91| 在线视频色国产色| 久热这里只有精品99| 国产欧美日韩精品亚洲av| 欧美成人一区二区免费高清观看 | a在线观看视频网站| av片东京热男人的天堂| 国产熟女午夜一区二区三区| 多毛熟女@视频| 一级毛片高清免费大全| 777久久人妻少妇嫩草av网站| 欧美日本亚洲视频在线播放| 欧美在线黄色| 国产成人精品在线电影| 51午夜福利影视在线观看| 黄频高清免费视频| 日韩大码丰满熟妇| 国产成人一区二区三区免费视频网站| 亚洲欧洲精品一区二区精品久久久| 精品国产国语对白av| 黑人巨大精品欧美一区二区mp4| 男女午夜视频在线观看| 搡老妇女老女人老熟妇| 淫妇啪啪啪对白视频| 色哟哟哟哟哟哟| 视频区欧美日本亚洲| 国产97色在线日韩免费| 黄网站色视频无遮挡免费观看| 久久香蕉国产精品| 性欧美人与动物交配| 手机成人av网站| 国产91精品成人一区二区三区| 日韩有码中文字幕| 国产熟女xx| 日韩 欧美 亚洲 中文字幕| 欧美成人一区二区免费高清观看 | 午夜福利,免费看| 久久香蕉激情| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品美女久久久久99蜜臀| 国产国语露脸激情在线看| 欧美性长视频在线观看| 一进一出好大好爽视频| 亚洲一区高清亚洲精品| 国产成人欧美| 免费人成视频x8x8入口观看| 精品久久久久久久人妻蜜臀av | 午夜老司机福利片| 国产亚洲av高清不卡| 在线天堂中文资源库| e午夜精品久久久久久久| 国产精品亚洲美女久久久| 国产精品1区2区在线观看.| 日本一区二区免费在线视频| 99久久精品国产亚洲精品| 露出奶头的视频| 亚洲黑人精品在线| 久久九九热精品免费| 亚洲精品国产色婷婷电影| 中文字幕精品免费在线观看视频| 在线观看www视频免费| 国产一卡二卡三卡精品| 美女高潮到喷水免费观看| 亚洲一区高清亚洲精品| 国产视频一区二区在线看| 成人亚洲精品一区在线观看| 九色国产91popny在线| 欧美+亚洲+日韩+国产| 亚洲欧美一区二区三区黑人| 精品久久久久久久人妻蜜臀av | 国产精品自产拍在线观看55亚洲| 国产aⅴ精品一区二区三区波| 亚洲全国av大片| 久久久久久久精品吃奶| 9191精品国产免费久久| 日韩欧美三级三区| x7x7x7水蜜桃| 亚洲中文日韩欧美视频| 久久中文字幕人妻熟女| 亚洲av日韩精品久久久久久密| 国产成人精品在线电影| 国产色视频综合| 国产精品久久电影中文字幕| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区mp4| 色综合亚洲欧美另类图片| 久久九九热精品免费| 国产精品亚洲美女久久久| 亚洲片人在线观看| 在线观看免费日韩欧美大片| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线| 欧美中文日本在线观看视频| 亚洲三区欧美一区| 一夜夜www| 久久久久久久久免费视频了| 久久久国产欧美日韩av| 99香蕉大伊视频| 亚洲最大成人中文| 天堂动漫精品| 久久久久国内视频| 搡老妇女老女人老熟妇| 亚洲男人天堂网一区| 91国产中文字幕| 国产成人啪精品午夜网站| 亚洲av片天天在线观看| 精品电影一区二区在线| 久久久精品欧美日韩精品| 午夜福利视频1000在线观看 | 欧美成人午夜精品| 动漫黄色视频在线观看| 久久精品国产亚洲av香蕉五月| 搡老妇女老女人老熟妇| 少妇的丰满在线观看| 欧美乱码精品一区二区三区| 神马国产精品三级电影在线观看 | 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 在线国产一区二区在线| 一二三四社区在线视频社区8| 午夜福利成人在线免费观看| 国产av在哪里看| 亚洲人成电影观看| 亚洲一区中文字幕在线| 欧美日韩乱码在线| 精品国产乱子伦一区二区三区| 丁香欧美五月| 非洲黑人性xxxx精品又粗又长| 亚洲专区国产一区二区| cao死你这个sao货| 99国产精品一区二区三区| 国产午夜精品久久久久久| 人妻久久中文字幕网| 成人国产一区最新在线观看| www.www免费av| 一本综合久久免费| 天堂动漫精品| 九色国产91popny在线| 国产真人三级小视频在线观看| 日本一区二区免费在线视频| 亚洲成人久久性| 麻豆av在线久日| 又黄又粗又硬又大视频| 国产av精品麻豆| 精品国产乱子伦一区二区三区| 国产真人三级小视频在线观看| 麻豆av在线久日| 国产主播在线观看一区二区| 搡老岳熟女国产| 亚洲一码二码三码区别大吗| 高清在线国产一区| 亚洲精品中文字幕在线视频| 国产成人精品在线电影| 国产激情久久老熟女| 亚洲色图av天堂| 18禁国产床啪视频网站| av网站免费在线观看视频| 亚洲欧美精品综合久久99| 国产精品免费视频内射| 十八禁网站免费在线| 黄色片一级片一级黄色片| www日本在线高清视频| 午夜福利一区二区在线看| 18美女黄网站色大片免费观看| 九色国产91popny在线| 97碰自拍视频| av视频在线观看入口| 国内精品久久久久精免费| 搡老岳熟女国产| tocl精华| 精品卡一卡二卡四卡免费| 国产男靠女视频免费网站| 日韩欧美免费精品| 精品日产1卡2卡| 女同久久另类99精品国产91| 中文字幕人妻丝袜一区二区| 黄片大片在线免费观看| 真人做人爱边吃奶动态| 久久人人97超碰香蕉20202| 大陆偷拍与自拍| 黄色丝袜av网址大全| 变态另类成人亚洲欧美熟女 | 久久久国产成人免费| 好看av亚洲va欧美ⅴa在| 一a级毛片在线观看| 久久久久九九精品影院| 国产99久久九九免费精品| 日韩欧美国产一区二区入口| 亚洲精华国产精华精| 成人精品一区二区免费| 色在线成人网| 大香蕉久久成人网| 69精品国产乱码久久久| 成人精品一区二区免费| 99re在线观看精品视频| 日韩大码丰满熟妇| 国产成人免费无遮挡视频| 亚洲欧美日韩无卡精品| 麻豆av在线久日| 热re99久久国产66热| 两个人免费观看高清视频| 日本一区二区免费在线视频| 搞女人的毛片| 久久亚洲精品不卡| 久久午夜亚洲精品久久| 久久 成人 亚洲| 男人舔女人下体高潮全视频| 母亲3免费完整高清在线观看| 久久人人爽av亚洲精品天堂| 国产99白浆流出| 久久人人精品亚洲av| 免费高清在线观看日韩| 国产成人一区二区三区免费视频网站| 国产精品野战在线观看| 最好的美女福利视频网| 搡老妇女老女人老熟妇| 精品久久久久久成人av| 国产一级毛片七仙女欲春2 | 90打野战视频偷拍视频| 免费久久久久久久精品成人欧美视频| 一本大道久久a久久精品| 亚洲第一欧美日韩一区二区三区| 久久国产亚洲av麻豆专区| 午夜福利欧美成人| 亚洲国产精品成人综合色| 国产人伦9x9x在线观看| 国产亚洲精品久久久久久毛片| 欧美激情高清一区二区三区| 久久影院123| www.自偷自拍.com| 美女 人体艺术 gogo| 久久国产精品影院| 制服丝袜大香蕉在线| 欧美最黄视频在线播放免费| 成人免费观看视频高清| 黑人操中国人逼视频| 露出奶头的视频| 亚洲中文日韩欧美视频| 99国产综合亚洲精品| 久久午夜综合久久蜜桃| 欧美av亚洲av综合av国产av| 国产亚洲欧美在线一区二区| 亚洲av美国av| 母亲3免费完整高清在线观看| 色综合欧美亚洲国产小说| 动漫黄色视频在线观看| 国产精品久久视频播放| 中文字幕精品免费在线观看视频| 日本精品一区二区三区蜜桃| 亚洲av熟女| 亚洲狠狠婷婷综合久久图片| 操美女的视频在线观看| 无限看片的www在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 一本综合久久免费| www.999成人在线观看| 精品国产一区二区久久| 亚洲中文字幕日韩| 国产精品爽爽va在线观看网站 | 88av欧美| 极品人妻少妇av视频| 欧美国产精品va在线观看不卡| 久久精品影院6| 国产片内射在线| 亚洲第一av免费看| 一进一出抽搐动态| 黄色成人免费大全| 欧美日韩乱码在线| 久久人妻av系列| 日韩欧美一区视频在线观看| 国产极品粉嫩免费观看在线| 国产精品国产高清国产av| 老汉色av国产亚洲站长工具| 国产亚洲欧美98| 亚洲专区国产一区二区| 啦啦啦 在线观看视频| 亚洲欧美精品综合久久99| 久久香蕉精品热| АⅤ资源中文在线天堂| 成人三级做爰电影| 999精品在线视频| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频| 国产熟女xx| 欧美老熟妇乱子伦牲交| 美女 人体艺术 gogo| 久久久久久免费高清国产稀缺| 日韩欧美免费精品| 黑人欧美特级aaaaaa片| 亚洲国产精品999在线| 99久久国产精品久久久| 国产精品1区2区在线观看.| 日本a在线网址| 精品国产国语对白av| 88av欧美| 在线观看日韩欧美| 亚洲欧美激情在线| 亚洲精华国产精华精| 欧美绝顶高潮抽搐喷水| 淫秽高清视频在线观看| 最近最新中文字幕大全电影3 | 高潮久久久久久久久久久不卡| 男人舔女人的私密视频| 亚洲人成电影观看| 国产精品亚洲美女久久久| avwww免费| 日韩一卡2卡3卡4卡2021年| 日韩中文字幕欧美一区二区| 亚洲精品中文字幕一二三四区| 精品一区二区三区av网在线观看| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区精品| 亚洲 欧美 日韩 在线 免费| 亚洲一码二码三码区别大吗| 老司机深夜福利视频在线观看| 精品久久久久久久毛片微露脸| 多毛熟女@视频| 国产又色又爽无遮挡免费看| 精品欧美国产一区二区三| 九色亚洲精品在线播放| 免费高清在线观看日韩| 两个人看的免费小视频| 男人舔女人下体高潮全视频| 两人在一起打扑克的视频| 精品国产超薄肉色丝袜足j| 丝袜美腿诱惑在线| 色播在线永久视频| 天堂√8在线中文| 国产精品 欧美亚洲| 亚洲国产看品久久| 国产麻豆69| 国产精品一区二区精品视频观看| 国产成人精品无人区| 国产成人精品久久二区二区免费| 亚洲,欧美精品.| 人妻久久中文字幕网| 欧美日韩黄片免| 久久精品国产亚洲av高清一级| 国产亚洲精品综合一区在线观看 | 久久人妻熟女aⅴ| 日韩大尺度精品在线看网址 | 91老司机精品| 国产精品一区二区精品视频观看| 欧美激情 高清一区二区三区| 亚洲国产精品久久男人天堂| 在线观看一区二区三区| 国产单亲对白刺激| 女人高潮潮喷娇喘18禁视频| 日本vs欧美在线观看视频| 亚洲 欧美 日韩 在线 免费| 久久久久久久精品吃奶| 成人永久免费在线观看视频| 91成年电影在线观看| 精品一区二区三区av网在线观看| 手机成人av网站| 男女下面插进去视频免费观看| 色老头精品视频在线观看| 美女高潮到喷水免费观看| 午夜免费鲁丝| 88av欧美| 黑人欧美特级aaaaaa片| 午夜精品国产一区二区电影| 国产黄a三级三级三级人| 国产av一区在线观看免费| 日韩 欧美 亚洲 中文字幕| 首页视频小说图片口味搜索| 日韩国内少妇激情av| 国产精品影院久久| 真人一进一出gif抽搐免费| 激情视频va一区二区三区| 夜夜爽天天搞| 给我免费播放毛片高清在线观看| 手机成人av网站| 满18在线观看网站| 好看av亚洲va欧美ⅴa在| 久久久精品欧美日韩精品| 国产激情欧美一区二区| 国产欧美日韩精品亚洲av| xxx96com| 黑人巨大精品欧美一区二区蜜桃| 久久天躁狠狠躁夜夜2o2o| 一级作爱视频免费观看| 在线观看舔阴道视频| 久久中文看片网| 国产精品美女特级片免费视频播放器 | 成人手机av| 亚洲专区中文字幕在线| 色综合欧美亚洲国产小说| 99精品欧美一区二区三区四区| 欧美乱码精品一区二区三区| 日韩大码丰满熟妇| 国产av在哪里看| 亚洲第一电影网av| 亚洲人成77777在线视频| 色老头精品视频在线观看| 好男人在线观看高清免费视频 | 亚洲七黄色美女视频| 欧美成人午夜精品| 欧美色欧美亚洲另类二区 | 窝窝影院91人妻| 最近最新中文字幕大全免费视频| 老司机靠b影院| 91麻豆av在线| 变态另类丝袜制服| 欧美国产精品va在线观看不卡| 久久午夜亚洲精品久久| 亚洲国产毛片av蜜桃av| 9色porny在线观看| 黄频高清免费视频| 老汉色∧v一级毛片| 首页视频小说图片口味搜索| 99精品欧美一区二区三区四区| 久久精品亚洲精品国产色婷小说| 国产色视频综合| 在线视频色国产色| 怎么达到女性高潮| 亚洲一码二码三码区别大吗| 免费在线观看影片大全网站| 精品欧美一区二区三区在线| 久久久水蜜桃国产精品网| 欧美绝顶高潮抽搐喷水| 波多野结衣一区麻豆| 亚洲第一av免费看| 99精品在免费线老司机午夜| 最近最新中文字幕大全免费视频| 久久久久久人人人人人| 国产成人啪精品午夜网站| 国产精品日韩av在线免费观看 | 亚洲成av人片免费观看| 亚洲第一av免费看| 窝窝影院91人妻| 一个人观看的视频www高清免费观看 | 国产精品二区激情视频| 欧美国产精品va在线观看不卡| 亚洲五月婷婷丁香| 久久人人97超碰香蕉20202| 女人精品久久久久毛片| 亚洲国产看品久久| 国产日韩一区二区三区精品不卡| 手机成人av网站| 久久精品国产亚洲av高清一级| 少妇熟女aⅴ在线视频| 18禁裸乳无遮挡免费网站照片 | 脱女人内裤的视频| 欧美激情高清一区二区三区| 在线观看66精品国产| 99在线视频只有这里精品首页| 欧美+亚洲+日韩+国产| 操美女的视频在线观看| av在线播放免费不卡| 中文字幕最新亚洲高清| 国产欧美日韩精品亚洲av| 久久精品国产综合久久久| 国产成人精品无人区| 动漫黄色视频在线观看| 三级毛片av免费| 久久国产精品男人的天堂亚洲| 午夜福利成人在线免费观看| 在线免费观看的www视频| 久久 成人 亚洲| 51午夜福利影视在线观看| 欧美成人午夜精品| av福利片在线| 亚洲国产精品sss在线观看| 亚洲第一电影网av| 精品乱码久久久久久99久播| 看片在线看免费视频| 国产成人精品无人区| 国产一区二区三区视频了| 亚洲成a人片在线一区二区| av欧美777| xxx96com| 久久精品人人爽人人爽视色| 亚洲专区字幕在线| 日本 欧美在线| 日本黄色视频三级网站网址| 给我免费播放毛片高清在线观看| 中文字幕人成人乱码亚洲影| 久久影院123| 琪琪午夜伦伦电影理论片6080| 在线播放国产精品三级| 精品卡一卡二卡四卡免费| 淫秽高清视频在线观看| 香蕉久久夜色| 在线观看午夜福利视频| 大型黄色视频在线免费观看| 色综合欧美亚洲国产小说| 岛国视频午夜一区免费看| 国产熟女xx| 免费看十八禁软件| 欧美中文综合在线视频| 熟妇人妻久久中文字幕3abv| 亚洲熟妇熟女久久| 99久久综合精品五月天人人| 婷婷六月久久综合丁香| 午夜免费鲁丝|