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

    磨損工況下可傾瓦軸承性能研究

    2021-11-11 13:52:20唐心昊王玉君許偉偉王振波
    關鍵詞:瓦塊偏心率油膜

    李 強, 唐心昊, 張 碩, 王玉君, 許偉偉, 王振波

    (中國石油大學(華東) 新能源學院,山東青島 266580)

    隨著石油機械行業(yè)的發(fā)展,壓縮機、離心泵等旋轉(zhuǎn)機械逐漸大型化、高速化,對滑動軸承性能要求不斷提高,可傾瓦軸承因其高穩(wěn)定性得到了廣泛應用。滑動軸承性能的研究一直受到學者們的高度重視。Kayar等[1]和Li等[2]基于Reynolds方程的求解研究了偏心率、轉(zhuǎn)速和潤滑介質(zhì)特性對軸承性能的影響。在求解Reynolds方程時耦合求解能量方程進而研究考慮熱效應下固定瓦軸承[3-4]和可傾瓦軸承[5-6]性能。計算機技術的發(fā)展使求解N-S方程成為可能,Keogh等[7]開發(fā)了CFD程序,研究了軸向開槽的軸承溫度分布;李強等[8-11]采用CFD方法,結(jié)合動網(wǎng)格技術研究了軸承動靜特性以及軸承內(nèi)不凝結(jié)氣體;Ding等[12]基于CFD方法研究了考慮空氣夾帶下可傾瓦軸承的摩擦損失;Armentrout等[13]分別采用Reynolds方程和CFD方法分析了湍流及對流慣性項對可傾瓦軸承的影響;Parovay等[14]采用CFD方法研究了可傾瓦軸承性能與載荷、轉(zhuǎn)速和徑向間隙的關系?;赗eynolds方程和CFD方法,逐漸形成了滑動軸承熱流體潤滑(THD)模型。在啟?;蜣D(zhuǎn)子失穩(wěn)的工況下,軸瓦發(fā)生磨損引起其形狀的變化,進而改變軸承性能。磨損工況下軸承性能的研究逐漸展開[15-16]。但目前軸承磨損的研究大多為基于Reynolds方程的固定瓦軸承[17-18],采用CFD方法的可傾瓦軸承的研究較為少見。筆者采用CFD方法、結(jié)合自編動網(wǎng)格程序建立可傾瓦軸承三維熱流體動壓潤滑模型(thermo-hydrodynamic, THD),實現(xiàn)潤滑流場、轉(zhuǎn)子運動以及瓦塊轉(zhuǎn)動三者之間的耦合。在此基礎上研究不同磨損深度對可傾瓦軸承潤滑性能的影響。

    1 控制方程

    1.1 流體控制方程

    基于CFD方法對可傾瓦軸承三維流場進行求解,同時考慮潤滑介質(zhì)黏溫效應及空化的影響。

    質(zhì)量守恒方程為

    (1)

    式中,ρm為混合相密度,kg/m3;vm為混合相速度矢量,m/s。

    動量守恒方程為

    (2)

    其中

    式中,p為壓力,Pa;μm為混合相黏度,Pa·s;g為自由落體加速度矢量,m/s2;F為外體積力,N;ρv為氣相密度,kg/m3。

    轉(zhuǎn)軸的偏心以及各瓦塊的轉(zhuǎn)動,導致軸頸與軸承之間油膜間隙分布不均勻,進而產(chǎn)生收斂楔和發(fā)散楔。潤滑油進入軸承發(fā)散楔,油膜壓力降低。當油膜壓力低于潤滑油的飽和蒸氣壓時,則會發(fā)生空化,產(chǎn)生氣態(tài)潤滑油。為了更好地預測軸承性能,采用組分輸運方程描述潤滑介質(zhì)氣、液兩相組分之間的輸運,表示為

    (3)

    式中,fk為第k相質(zhì)量分數(shù);Re和Rc分別為蒸發(fā)和冷凝速率。

    質(zhì)量平均速度和混合密度可描述為

    (4)

    (5)

    式中,ρk為第k相密度,kg/m3;vk為第k相速度矢量,m/s。

    組分輸運方程中源項由Singhal等[19]提出的“全空化模型”計算。該模型基于多相流框架,考慮了所有一階效應,可以處理與流體相變相關的大密度變化,而無需預先確定空化的位置、程度或類型。空化模型描述為

    (6)

    (7)

    式中,Ce和Cc為經(jīng)驗常數(shù);fg為空氣質(zhì)量分數(shù);pv為空化壓力,Pa。

    1.2 熱傳遞控制方程

    能量守恒方程為

    (8)

    式中,v為運動黏度,m2/s;keff為有效導熱率,W/(m·K);hk為第k相焓,kJ/kg;Shf為體積熱源項,W/m3。

    流體與邊界之間能量的傳遞采用對流換熱邊界條件。其熱通量可表示為

    q=hf(Tw-Tf)+qrad=hext(Text-Tw).

    (9)

    式中,q為熱通量,W/m2;hf為流體側(cè)局部傳熱系數(shù),W/(m2·K);Tw為邊界表面溫度,K;Tf為局部流體溫度,K;qrad為輻射熱通量,W/m2;hext為外部傳熱系數(shù),W/(m2·K);Text為外部散熱溫度,K。

    潤滑油黏度對溫度變化十分敏感,黏度隨溫度升高而迅速降低。黏溫效應可以用Walther[20]方程表征,其形式為

    loglog(v+c)=a-blog(T).

    (10)

    式中,T表示絕對溫度,K;c通常取0.7;a、b由40、100 ℃時潤滑油黏度計算得到,分別為11.88和4.679。

    1.3 邊界條件及離散方法

    軸承進出口邊界條件為壓力邊界條件,進口壓力為0.2 MPa,出口壓力為大氣壓,瓦塊上表面導熱系數(shù)為350 W/(m2·℃),其他為絕熱,固壁為無滑移邊界?;谟邢摅w積法離散控制方程,選擇基于壓力求解器,壓力速度耦合采用SIMPLE算法,壓力差分格式采用PRESTO!。

    在本文所有的計算工況下,可傾瓦軸承內(nèi)流場雷諾數(shù)均小于400,因此流動狀態(tài)可認為層流。為進行區(qū)分,雷諾數(shù)Ren表示為

    Ren=ρvl/μ.

    (11)

    式中,ρ為潤滑介質(zhì)密度,kg/m3;v為特征速度,m/s;l為特征長度,m;μ為動力黏度,Pa·s。

    由式(11)計算最大雷諾數(shù)為238.31。

    2 流域模型及計算方法

    2.1 流域模型

    本文的研究對象為工程上應用最廣泛的五瓦塊可傾瓦軸承,結(jié)構參考文獻[18],軸承及潤滑介質(zhì)主要參數(shù)如下:軸承半徑和長度分別為30和35 mm,軸瓦數(shù)為5,支點偏心為0.5,半徑間隙為0.09 mm, 瓦塊包角為57°, 瓦塊厚度為20 mm, 瓦塊轉(zhuǎn)動慣量為0.000 12 kg·m2, 轉(zhuǎn)子質(zhì)量為55.7 kg, 轉(zhuǎn)速為10 200 r/min。三維流域計算模型如圖1所示,包括進油口、瓦塊間間隙、油膜間隙以及瓦塊背

    圖1 流域模型Fig.1 Watershed model

    部間隙等。軸承加載方式為瓦間加載,供油方式為瓦間徑向供油,潤滑油從軸承兩側(cè)間隙流出。32#潤滑油主要參數(shù)為密度、定壓比熱容、導熱系數(shù)、氣態(tài)密度和空化壓力,其值分別為860 kg/m3、2 025 J/(kg·℃)、0.129 W/(m·℃)、1.2 kg/m3和29 185 Pa。

    加載瓦塊油膜厚度較小,容易發(fā)生磨損。因此假設加載瓦塊(3、4)發(fā)生磨損,磨損位置為瓦塊上表面中心,磨損面積為瓦塊上表面總面積的20%。不同的磨損程度通過磨損深度表征,磨損深度分別為10、20、30、40、50 μm。為滿足網(wǎng)格要求,磨損邊緣采用直邊段過渡。瓦塊磨損示意圖及磨損處網(wǎng)格如圖2所示。

    圖2 瓦塊磨損示意圖及網(wǎng)格劃分Fig.2 Schematic diagram and mesh of pad wear

    2.2 瞬態(tài)計算過程

    在轉(zhuǎn)子由初始位置運動至各工況下靜平衡位置的瞬態(tài)過程中,既包括轉(zhuǎn)子的轉(zhuǎn)動和擾動,也包括各瓦塊的自由擺動。為了準確模擬上述瞬態(tài)過程,通過自編動網(wǎng)格程序?qū)崿F(xiàn)潤滑流場、轉(zhuǎn)子運動以及各瓦塊轉(zhuǎn)動三者之間的耦合。

    假設轉(zhuǎn)子僅受到由轉(zhuǎn)子質(zhì)量產(chǎn)生的重力及軸承產(chǎn)生的油膜力作用。轉(zhuǎn)子所受合力為零時,可認為其處于靜平衡位置。其動力學方程為

    (12)

    式中,M為轉(zhuǎn)子質(zhì)量,kg;Fx和Fy分別為油膜力x、y方向的分力,N。

    由轉(zhuǎn)子表面油膜壓力積分得

    (13)

    式中,R為轉(zhuǎn)子半徑,m;L為軸承長度,m;p為油膜壓力,Pa;θ為壓力與x軸正方向夾角,rad。

    各瓦塊繞支點的轉(zhuǎn)動角速度由剛體繞定軸轉(zhuǎn)動微分方程得

    (14)

    式中,Jop為瓦塊對支點的轉(zhuǎn)動慣量,kg·m2;ωp為瓦塊轉(zhuǎn)動角速度,rad/s;Mop(Fp)為瓦塊受到的合力矩,N·m。

    瞬態(tài)計算過程如圖3所示,通過式(13)積分得轉(zhuǎn)子表面的瞬態(tài)油膜力,求解式(12)得到轉(zhuǎn)子x、y方向的加速度、速度及位移。由瓦塊上下表面單位面積受到的油膜力及與瓦塊支點的距離得到瓦塊受到的力矩,代入式(14)得到瓦塊的角加速度、角速度及轉(zhuǎn)動角度。通過自編動網(wǎng)格程序(UDF)實現(xiàn)可傾瓦軸承網(wǎng)格更新及瓦塊轉(zhuǎn)動。瓦塊轉(zhuǎn)動后網(wǎng)格如圖4所示,為更好顯示網(wǎng)格更新效果,將軸承間隙及瓦塊背部間隙放大并加大了瓦塊轉(zhuǎn)動角度。在轉(zhuǎn)子移動和各個瓦塊轉(zhuǎn)動后,網(wǎng)格節(jié)點實現(xiàn)了精確移動,即使較大的瓦塊轉(zhuǎn)動角度也能保證較高的網(wǎng)格質(zhì)量。

    圖3 瞬態(tài)計算過程Fig.3 Transient calculation process

    圖4 網(wǎng)格更新效果Fig.4 Mesh update effect

    2.3 網(wǎng)格無關性驗證

    采用的動網(wǎng)格程序需要實時讀取網(wǎng)格節(jié)點位置并進行更新。因此需要對計算模型進行網(wǎng)格無關性驗證,以在保證計算準確性的同時提高計算效率。參考文獻[21],表1給出了偏心率為0.4工況下軸承最大油膜壓力隨油膜間隙網(wǎng)格數(shù)變化情況,在5層徑向網(wǎng)格時軸承最大壓力變化低于1%,因此選擇徑向網(wǎng)格層數(shù)為5層,網(wǎng)格數(shù)為336 000。

    表1 網(wǎng)格無關性

    3 結(jié)果與討論

    3.1 模型準確性驗證

    將供油溫度37 ℃工況下軸承瓦塊(1、3、4)最高溫度數(shù)值計算結(jié)果與試驗數(shù)據(jù)[18]進行對比,如圖5所示。由圖5可知,不同磨損深度下瓦塊溫度模擬結(jié)果與試驗數(shù)據(jù)吻合較好,加載瓦塊(3、4)溫度隨磨損深度增加逐漸升高,非加載瓦塊(1)反之,最大溫度誤差出現(xiàn)在瓦塊3,數(shù)值為2.62 ℃,小于文獻[18]中瓦塊1最大溫度誤差5 ℃,證明了采用CFD方法建立的三維THD模型可以更好地預測可傾瓦軸承性能。

    圖5 溫度模擬與試驗結(jié)果對比Fig.5 Comparison of temperature simulation and experimental results

    3.2 軸承壓力、溫度分布

    圖6給出了軸承在供油溫度27 ℃、未發(fā)生磨損工況下瓦塊上表面壓力、溫度分布云圖。由圖6(a)可知,在未發(fā)生磨損的工況下,由于瓦塊的擺動,每個瓦塊均存在高壓區(qū),產(chǎn)生動壓效應,從而提高了可傾瓦軸承的穩(wěn)定性。瓦塊上表面高壓區(qū)沿轉(zhuǎn)子旋轉(zhuǎn)方向(逆時針方向)略有偏移,加載瓦塊(3、4)壓力明顯高于其他瓦塊,起主要承載作用,非加載瓦塊1壓力最低。由圖6(b)可知,在周向方向瓦塊上表面溫度沿轉(zhuǎn)子旋轉(zhuǎn)方向逐漸升高,且在瓦塊后緣處達到最高,瓦塊間間隙供油,潤滑油在間隙內(nèi)混合使得在進入下一油膜間隙時溫度得到了有效降低;在軸向方向,由于供油位置在瓦塊間隙的中間,瓦塊上表面溫度出現(xiàn)了中間低兩邊高的現(xiàn)象;轉(zhuǎn)子偏心導致加載瓦塊(3、4)油膜厚度減小,黏性耗散更加嚴重,因此其溫升明顯高于非加載瓦塊。

    圖6 瓦塊上表面壓力、溫度分布云圖Fig.6 Pressure and temperature contours of pad

    3.3 磨損深度對壓力、溫度影響

    圖7給出了供油溫度27 ℃、不同磨損深度下各瓦塊上表面壓力及溫度變化,數(shù)據(jù)取自軸承中面與瓦塊上表面交線,其中橫坐標表示軸承圓周方向瓦塊位置。

    圖7 不同磨損深度下各瓦塊表面壓力及溫度Fig.7 Surface pressure and temperature of pad for various wear depth

    由圖7可知,在未發(fā)生磨損工況下,瓦塊表面壓力、溫度分布連續(xù)性變化;在磨損工況下,加載瓦塊(3、4)上表面壓力及溫度的連續(xù)性遭到破壞,且磨損深度越大連續(xù)性破壞越嚴重。在瓦塊磨損處壓力及溫度產(chǎn)生局部降低,但與未磨損工況相比,加載瓦塊最高壓力和最高溫度明顯升高,非加載瓦塊反之。隨著磨損深度增加,靜平衡位置逐漸下降,加載瓦塊油膜厚度減小,非加載瓦塊油膜厚度增大,因此加載瓦塊的最高壓力和溫度逐漸升高,非加載瓦塊最高壓力和溫度逐漸降低。

    3.4 磨損深度對軸承承載力影響

    為了研究磨損深度對軸承承載力影響,通過自編動網(wǎng)格程序?qū)⑥D(zhuǎn)子移動至偏心率0.5處,保證潤滑流場、瓦塊轉(zhuǎn)動與轉(zhuǎn)子運動之間的耦合,最終形成穩(wěn)定的油膜間隙和瓦塊轉(zhuǎn)動角,積分轉(zhuǎn)子表面油膜力即為軸承承載力。圖8給出了不同磨損深度下軸承承載力,由圖8可知,隨著磨損深度增加,軸承的承載力不斷降低,但軸承承載力的變化趨勢逐漸減緩。磨損導致瓦塊表面壓力連續(xù)性遭到破壞并引起局部油膜厚度增大,使得加載瓦塊局部壓力降低,因此軸承承載力不斷降低。由于磨損面積僅占瓦塊上表面面積的20%,磨損深度的變化對壓力的影響有限,因此承載力變化趨勢逐漸減緩。在一定的磨損面積(20%)下,軸承承載力對磨損深度的敏感性逐漸降低。

    圖8 不同磨損深度下承載力Fig.8 Loading capacity under different wear depths

    3.5 磨損深度對偏心率和靜平衡位置影響

    圖9給出了供油溫度27 ℃時不同磨損深度下軸心偏心率,由圖9可知:隨著磨損深度增加,軸心的偏心率不斷增大,但變化趨勢逐漸減緩,磨損深度在0~30 μm時,軸心偏心率變化約為10%;而磨損深度在30~50 μm時,軸心偏心率變化率約為2%。隨著磨損深度增加,軸承承載力降低,因此需要更大的偏心以獲得足夠的油膜力支撐轉(zhuǎn)子。在磨損深度為0~30 μm時,一方面,軸承承載力降低且變化較大,另一方面,在較小的偏心率下,由偏心變化引起壓力升高較小,則需要更大的偏心變化獲得足夠的壓力升高,因此在這兩方面因素的共同作用下,偏心率變化較大。磨損深度為30~50 μm時,承載力變化較小,且在高偏心率下,偏心變化引起的壓力升高較大,因此軸心偏心率變化較小。

    圖9 不同磨損深度下軸心偏心率Fig.9 Eccentricity under different wear depths

    圖10給出了不同磨損深度下軸心運動軌跡。

    圖10 不同磨損深度下軸心靜平衡位置及瓦塊轉(zhuǎn)動角Fig.10 Static balance position of shaft and rotation angle of pad under different wear depths

    由圖10可知,軸心由起始位置經(jīng)過擾動后最終到達靜平衡位置。隨著磨損深度增加,軸心運動軌跡不斷擴大,靜平衡位置逐漸降低且在x軸方向先增大后減小。為進一步分析靜平衡位置在x軸方向的變化,圖10(b)給出了不同靜平衡位置下各瓦塊的轉(zhuǎn)動角度。由圖10可知,隨著磨損深度增加,加載方向左側(cè)的瓦塊(2、3)由順時針轉(zhuǎn)動逐漸變?yōu)槟鏁r針轉(zhuǎn)動,加載方向右側(cè)的瓦塊(4、5)始終逆時針轉(zhuǎn)動且轉(zhuǎn)動角度不斷增大,由于不同磨損深度下各瓦塊產(chǎn)生不同的轉(zhuǎn)動,改變了轉(zhuǎn)子表面受力,使得軸心靜平衡位置在x軸方向呈現(xiàn)先增大后減小的變化趨勢。

    4 結(jié) 論

    (1)在所建立的磨損工況下可傾瓦軸承的三維熱流體動壓潤滑模型中,通過自編動網(wǎng)格程序?qū)崿F(xiàn)了潤滑流場、瓦塊轉(zhuǎn)動、轉(zhuǎn)子運動的固體-流體-固體的多向耦合計算,計算結(jié)果與試驗結(jié)果吻合較好,驗證了該計算模型的有效性。

    (2)可傾瓦軸承每個瓦塊均存在高壓區(qū),產(chǎn)生動壓效應,瓦塊上表面溫度沿轉(zhuǎn)子旋轉(zhuǎn)方向逐漸升高,且出現(xiàn)中間低兩邊高的現(xiàn)象。

    (3)瓦塊的磨損破壞了瓦塊表面壓力和溫度的連續(xù)性,出現(xiàn)了局部的壓力、溫度降低,磨損深度越大,連續(xù)性破壞越嚴重。

    (4)隨著磨損深度增加,軸承承載力逐漸降低,偏心率不斷增大,但變化趨勢逐漸減緩,軸心靜平衡位置逐漸降低且在x軸方向先增大后減小。

    猜你喜歡
    瓦塊偏心率油膜
    Hansen系數(shù)遞推的效率?
    天文學報(2021年5期)2021-10-09 06:38:52
    一種高效的頂點偏心率計算方法
    長城油膜軸承油在高速棒材生產(chǎn)線的應用
    石油商技(2021年1期)2021-03-29 02:36:04
    大型數(shù)控立式磨床靜壓轉(zhuǎn)臺油膜熱特性仿真及其實驗分析
    冷軋軋機油膜軸承系統(tǒng)故障分析與對策
    新疆鋼鐵(2016年4期)2016-05-25 00:37:38
    無縫鋼管壁厚偏心率的測量分析及降低方法
    鋼管(2016年1期)2016-05-17 06:12:44
    基于幾何條件的可傾瓦軸承油膜邊界條件判定方法
    測量計算法測定軸瓦間隙
    大偏心率軌道星上快速計算方法
    汽輪機可傾瓦顫振引發(fā)低頻振蕩問題研究
    亚洲人成网站在线播放欧美日韩| 99久久精品一区二区三区| 国产大屁股一区二区在线视频| 国产又黄又爽又无遮挡在线| 亚洲三级黄色毛片| 99久久人妻综合| 国产精品麻豆人妻色哟哟久久 | 人人妻人人看人人澡| 国产精品电影一区二区三区| 成人综合一区亚洲| 中文字幕av在线有码专区| 亚洲精品乱码久久久v下载方式| 99久国产av精品| 日本撒尿小便嘘嘘汇集6| 免费在线观看成人毛片| 精品久久久久久成人av| 亚洲美女搞黄在线观看| 日本三级黄在线观看| 精品久久久噜噜| 国产又黄又爽又无遮挡在线| 国产激情偷乱视频一区二区| 国产精品99久久久久久久久| 国产精品一区二区在线观看99 | 国产黄片视频在线免费观看| 综合色av麻豆| 日韩成人av中文字幕在线观看| 亚洲精品粉嫩美女一区| 欧美一区二区亚洲| 久久韩国三级中文字幕| 校园人妻丝袜中文字幕| 久久韩国三级中文字幕| 亚洲欧美日韩高清专用| 麻豆精品久久久久久蜜桃| 99热这里只有精品一区| 欧美最黄视频在线播放免费| 欧美日韩国产亚洲二区| 日本成人三级电影网站| 国产爱豆传媒在线观看| 国产精品无大码| 国产精品,欧美在线| 色综合站精品国产| 久久韩国三级中文字幕| 国产精品一区www在线观看| 特大巨黑吊av在线直播| 又黄又爽又刺激的免费视频.| 人人妻人人看人人澡| 内射极品少妇av片p| 全区人妻精品视频| 欧美又色又爽又黄视频| 成人国产麻豆网| 精品久久久久久久久av| 搡老妇女老女人老熟妇| 国产人妻一区二区三区在| 日本黄色片子视频| 尤物成人国产欧美一区二区三区| 看片在线看免费视频| 人妻制服诱惑在线中文字幕| 中文亚洲av片在线观看爽| 日韩精品青青久久久久久| 国产一区二区在线观看日韩| 中文精品一卡2卡3卡4更新| 欧美性猛交╳xxx乱大交人| 国模一区二区三区四区视频| 国产黄色视频一区二区在线观看 | 亚洲av电影不卡..在线观看| 99热这里只有是精品50| 伦精品一区二区三区| 高清午夜精品一区二区三区 | 性欧美人与动物交配| 亚洲丝袜综合中文字幕| 亚洲人成网站在线播放欧美日韩| 日日啪夜夜撸| 美女国产视频在线观看| 激情 狠狠 欧美| 亚洲精品久久国产高清桃花| 少妇人妻精品综合一区二区 | 欧美三级亚洲精品| 村上凉子中文字幕在线| 免费av毛片视频| 中国国产av一级| videossex国产| 美女 人体艺术 gogo| 亚洲国产精品久久男人天堂| 人人妻人人看人人澡| 国产老妇女一区| 日本与韩国留学比较| 我的女老师完整版在线观看| 国产精品爽爽va在线观看网站| 18禁在线无遮挡免费观看视频| 中文字幕精品亚洲无线码一区| 九九爱精品视频在线观看| 天堂av国产一区二区熟女人妻| 国产伦一二天堂av在线观看| 免费av不卡在线播放| 亚洲av中文av极速乱| 久久婷婷人人爽人人干人人爱| 国产视频首页在线观看| 欧美成人一区二区免费高清观看| 三级毛片av免费| 级片在线观看| 久久久久久久久久成人| av在线蜜桃| 亚洲欧洲国产日韩| 欧美性猛交黑人性爽| 日韩欧美精品免费久久| 色视频www国产| 男女那种视频在线观看| 国产免费男女视频| 高清午夜精品一区二区三区 | 如何舔出高潮| 久久这里有精品视频免费| 国产黄片美女视频| 人妻夜夜爽99麻豆av| 男女下面进入的视频免费午夜| 在线播放国产精品三级| 亚洲人成网站在线观看播放| 国产精品一区二区在线观看99 | 亚洲无线在线观看| 欧美成人一区二区免费高清观看| 欧美色视频一区免费| 特大巨黑吊av在线直播| 亚洲图色成人| 成熟少妇高潮喷水视频| 啦啦啦韩国在线观看视频| 三级经典国产精品| 亚洲真实伦在线观看| 五月玫瑰六月丁香| 嫩草影院入口| 干丝袜人妻中文字幕| 一级毛片aaaaaa免费看小| 夫妻性生交免费视频一级片| 91aial.com中文字幕在线观看| 长腿黑丝高跟| 欧美xxxx黑人xx丫x性爽| 国产亚洲精品久久久com| 亚洲久久久久久中文字幕| 国产伦理片在线播放av一区 | 成年女人永久免费观看视频| 午夜福利在线观看吧| 久久精品久久久久久久性| 大又大粗又爽又黄少妇毛片口| av天堂中文字幕网| 日韩中字成人| 国产视频首页在线观看| 高清毛片免费观看视频网站| 成人二区视频| 国产在线男女| 日本-黄色视频高清免费观看| 久久婷婷人人爽人人干人人爱| 日韩强制内射视频| 性色avwww在线观看| 国产单亲对白刺激| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久久电影| 国产不卡一卡二| 99九九线精品视频在线观看视频| 亚洲国产精品成人久久小说 | 哪个播放器可以免费观看大片| 国国产精品蜜臀av免费| 男插女下体视频免费在线播放| 午夜福利成人在线免费观看| 久久久久久久久久成人| 亚洲精品自拍成人| 哪里可以看免费的av片| 男人的好看免费观看在线视频| 校园春色视频在线观看| 伦精品一区二区三区| 亚洲精品自拍成人| 少妇的逼好多水| 亚洲欧美清纯卡通| 亚洲精品乱码久久久v下载方式| 成人高潮视频无遮挡免费网站| 99国产极品粉嫩在线观看| 少妇的逼好多水| 18禁在线无遮挡免费观看视频| 久久草成人影院| 欧美日韩国产亚洲二区| 国产精品久久久久久亚洲av鲁大| 色哟哟·www| 亚洲国产精品成人久久小说 | 亚洲熟妇中文字幕五十中出| 国产成人影院久久av| 99热这里只有精品一区| 少妇熟女aⅴ在线视频| 午夜福利在线在线| 中文精品一卡2卡3卡4更新| 99久久精品国产国产毛片| 精品日产1卡2卡| 色5月婷婷丁香| 亚洲av中文av极速乱| 国产精品爽爽va在线观看网站| 国产黄色视频一区二区在线观看 | 国产av一区在线观看免费| 国产精品一区二区三区四区免费观看| 2021天堂中文幕一二区在线观| 最近视频中文字幕2019在线8| 午夜激情福利司机影院| 小蜜桃在线观看免费完整版高清| av在线蜜桃| 国产国拍精品亚洲av在线观看| 高清毛片免费观看视频网站| av黄色大香蕉| 99视频精品全部免费 在线| 精品一区二区免费观看| 午夜福利在线观看吧| 精品久久久久久久久久久久久| 黄色日韩在线| 久久99热这里只有精品18| 亚洲精品影视一区二区三区av| 在线播放国产精品三级| 小说图片视频综合网站| 亚洲真实伦在线观看| 中文字幕熟女人妻在线| 亚洲国产色片| 大又大粗又爽又黄少妇毛片口| 边亲边吃奶的免费视频| 欧美精品国产亚洲| 国内精品宾馆在线| 黄色欧美视频在线观看| 欧美三级亚洲精品| 亚洲av成人av| 欧美色视频一区免费| av专区在线播放| 国产高潮美女av| 亚洲精品粉嫩美女一区| 搡老妇女老女人老熟妇| 亚洲国产精品成人久久小说 | 哪里可以看免费的av片| 亚洲国产色片| 亚洲丝袜综合中文字幕| 男人狂女人下面高潮的视频| 亚洲性久久影院| 亚洲国产精品sss在线观看| av在线播放精品| 欧美高清性xxxxhd video| 两性午夜刺激爽爽歪歪视频在线观看| av在线观看视频网站免费| 深爱激情五月婷婷| 91精品国产九色| 伦理电影大哥的女人| 一区福利在线观看| 亚洲va在线va天堂va国产| av在线老鸭窝| 女同久久另类99精品国产91| 国产精品av视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 观看免费一级毛片| 寂寞人妻少妇视频99o| АⅤ资源中文在线天堂| 午夜激情福利司机影院| 特级一级黄色大片| 国产精品日韩av在线免费观看| 精品人妻一区二区三区麻豆| 一区福利在线观看| 国产v大片淫在线免费观看| 边亲边吃奶的免费视频| 给我免费播放毛片高清在线观看| 成人亚洲欧美一区二区av| 少妇裸体淫交视频免费看高清| 久久久精品大字幕| 色噜噜av男人的天堂激情| 三级毛片av免费| 久久精品91蜜桃| 日日摸夜夜添夜夜爱| 国产精品人妻久久久久久| 亚洲av免费在线观看| 中文字幕免费在线视频6| 男女啪啪激烈高潮av片| 麻豆国产av国片精品| 免费av不卡在线播放| 美女xxoo啪啪120秒动态图| 日韩大尺度精品在线看网址| 美女黄网站色视频| 国内久久婷婷六月综合欲色啪| 亚洲成人精品中文字幕电影| 麻豆乱淫一区二区| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 欧美日韩综合久久久久久| 国产男人的电影天堂91| 美女 人体艺术 gogo| 国产在线男女| 男人舔奶头视频| 偷拍熟女少妇极品色| 亚洲国产精品久久男人天堂| 午夜福利在线观看免费完整高清在 | 久久久久久大精品| 99精品在免费线老司机午夜| 中文精品一卡2卡3卡4更新| 久久久久久久久久久丰满| 深夜a级毛片| 中文字幕精品亚洲无线码一区| 日本黄大片高清| 欧美xxxx黑人xx丫x性爽| 少妇熟女aⅴ在线视频| 中国国产av一级| 日韩一区二区视频免费看| 国产精品,欧美在线| 欧美高清性xxxxhd video| 欧美日韩在线观看h| 亚洲人与动物交配视频| 久久热精品热| 成人av在线播放网站| 人体艺术视频欧美日本| 97热精品久久久久久| 国产免费男女视频| 观看美女的网站| 国产成人午夜福利电影在线观看| 欧美一区二区国产精品久久精品| 日韩一区二区视频免费看| 日韩三级伦理在线观看| 国产色爽女视频免费观看| 亚洲美女视频黄频| 99久久人妻综合| 日韩中字成人| 性欧美人与动物交配| 亚洲精品久久久久久婷婷小说 | 如何舔出高潮| 3wmmmm亚洲av在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产极品天堂在线| 国产视频内射| 91aial.com中文字幕在线观看| 一级黄片播放器| 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 寂寞人妻少妇视频99o| 一本久久中文字幕| 天堂中文最新版在线下载 | 色综合色国产| 国产精品久久久久久av不卡| .国产精品久久| 女人十人毛片免费观看3o分钟| 国产伦一二天堂av在线观看| 听说在线观看完整版免费高清| av专区在线播放| 国产成人精品婷婷| 日韩一区二区视频免费看| 又粗又硬又长又爽又黄的视频 | avwww免费| 中出人妻视频一区二区| 精品久久久久久久久亚洲| 一进一出抽搐动态| 99国产极品粉嫩在线观看| 少妇人妻精品综合一区二区 | av在线天堂中文字幕| 久久久久久久亚洲中文字幕| 久久久久久久久久久丰满| 你懂的网址亚洲精品在线观看 | 亚洲成人精品中文字幕电影| 国产高清有码在线观看视频| 国产在线精品亚洲第一网站| 男女啪啪激烈高潮av片| 日本爱情动作片www.在线观看| 高清在线视频一区二区三区 | 国产v大片淫在线免费观看| 波野结衣二区三区在线| 国产精品三级大全| 欧美丝袜亚洲另类| 在线a可以看的网站| or卡值多少钱| 嘟嘟电影网在线观看| 久久精品久久久久久久性| 国语自产精品视频在线第100页| 国产成人精品婷婷| 女的被弄到高潮叫床怎么办| 精品一区二区三区人妻视频| 精品免费久久久久久久清纯| 中文字幕久久专区| 少妇丰满av| 中文字幕精品亚洲无线码一区| 女同久久另类99精品国产91| 免费看美女性在线毛片视频| 99热这里只有是精品50| 日韩精品青青久久久久久| 亚洲精品影视一区二区三区av| 在线国产一区二区在线| 婷婷六月久久综合丁香| 狂野欧美白嫩少妇大欣赏| 91狼人影院| 黑人高潮一二区| 国产中年淑女户外野战色| 亚洲最大成人中文| 亚洲av.av天堂| 色视频www国产| 亚洲国产欧美在线一区| 国产69精品久久久久777片| 日本黄色视频三级网站网址| 国产黄片视频在线免费观看| 欧美极品一区二区三区四区| 边亲边吃奶的免费视频| 国产黄片视频在线免费观看| 成人欧美大片| 国产精品久久久久久精品电影| 日韩 亚洲 欧美在线| 免费无遮挡裸体视频| 激情 狠狠 欧美| 老师上课跳d突然被开到最大视频| 最好的美女福利视频网| 久久午夜亚洲精品久久| 人妻制服诱惑在线中文字幕| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 精品人妻偷拍中文字幕| 日本av手机在线免费观看| 国产精品电影一区二区三区| 久久久成人免费电影| 久久精品夜色国产| 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频| 麻豆一二三区av精品| 丰满人妻一区二区三区视频av| 国产精品伦人一区二区| 天天躁夜夜躁狠狠久久av| 一本久久中文字幕| 亚洲,欧美,日韩| 国产成人福利小说| 高清毛片免费看| 亚洲精品久久久久久婷婷小说 | 久久精品久久久久久噜噜老黄 | 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩无卡精品| 国产真实乱freesex| 全区人妻精品视频| 夜夜爽天天搞| 免费看光身美女| 亚洲欧洲日产国产| 久久久久久久久久久丰满| av国产免费在线观看| 免费观看的影片在线观看| 欧美日本视频| 免费不卡的大黄色大毛片视频在线观看 | 免费av观看视频| 99热网站在线观看| 男人狂女人下面高潮的视频| 又粗又爽又猛毛片免费看| 一本久久中文字幕| 亚洲色图av天堂| 综合色丁香网| 亚洲人成网站在线观看播放| 五月玫瑰六月丁香| 亚洲国产精品久久男人天堂| 中文欧美无线码| 99热全是精品| 国产午夜福利久久久久久| 一夜夜www| 18+在线观看网站| 久久精品国产亚洲网站| 日韩,欧美,国产一区二区三区 | 99九九线精品视频在线观看视频| 一级毛片aaaaaa免费看小| 国产精品,欧美在线| 亚洲欧美日韩无卡精品| 亚洲欧美中文字幕日韩二区| 国产熟女欧美一区二区| 欧美激情在线99| av在线亚洲专区| 久久精品综合一区二区三区| 天堂网av新在线| 日日撸夜夜添| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 一级黄片播放器| 国产三级中文精品| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频 | 人妻久久中文字幕网| 久久久久久九九精品二区国产| 性色avwww在线观看| 欧美日本视频| 国产黄a三级三级三级人| 精品久久久久久久久久免费视频| 国产一区二区三区在线臀色熟女| 老熟妇乱子伦视频在线观看| 久久久久网色| 久久99蜜桃精品久久| 亚洲美女搞黄在线观看| 国产精品一区二区三区四区久久| 日韩欧美在线乱码| 六月丁香七月| 天天躁夜夜躁狠狠久久av| 嘟嘟电影网在线观看| 国产私拍福利视频在线观看| 国产成人午夜福利电影在线观看| 岛国在线免费视频观看| 老熟妇乱子伦视频在线观看| 日日摸夜夜添夜夜爱| 久久久久国产网址| 超碰av人人做人人爽久久| 亚洲欧美成人精品一区二区| a级毛片免费高清观看在线播放| 成人美女网站在线观看视频| 爱豆传媒免费全集在线观看| 久久精品国产鲁丝片午夜精品| 久久精品久久久久久噜噜老黄 | 真实男女啪啪啪动态图| 男女边吃奶边做爰视频| 九草在线视频观看| 久久久国产成人免费| 国产一区二区三区av在线 | a级毛片免费高清观看在线播放| 亚洲国产精品合色在线| 男人和女人高潮做爰伦理| 最好的美女福利视频网| 久久久久久久久久久丰满| 国产私拍福利视频在线观看| 国产av麻豆久久久久久久| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2| 日本黄色片子视频| 美女黄网站色视频| 97超碰精品成人国产| 六月丁香七月| 日韩精品青青久久久久久| 九九爱精品视频在线观看| 精品午夜福利在线看| 国产精品99久久久久久久久| 一本久久精品| 别揉我奶头 嗯啊视频| 有码 亚洲区| 国产激情偷乱视频一区二区| 久久韩国三级中文字幕| 日韩av不卡免费在线播放| 国产精品久久电影中文字幕| 久久久久久国产a免费观看| 国产精品一区二区三区四区免费观看| 在线免费十八禁| 亚洲国产精品sss在线观看| 99热这里只有是精品50| 免费看美女性在线毛片视频| 日本-黄色视频高清免费观看| 啦啦啦啦在线视频资源| 亚洲av免费在线观看| 在线播放无遮挡| 午夜老司机福利剧场| 2022亚洲国产成人精品| 一本精品99久久精品77| 一级二级三级毛片免费看| 亚洲中文字幕日韩| 国产单亲对白刺激| 又粗又硬又长又爽又黄的视频 | 丰满乱子伦码专区| 波野结衣二区三区在线| 一级黄片播放器| 高清日韩中文字幕在线| 色5月婷婷丁香| 久久国内精品自在自线图片| 日韩一区二区视频免费看| 丝袜美腿在线中文| 男女啪啪激烈高潮av片| 99久久成人亚洲精品观看| 亚洲最大成人av| 欧美区成人在线视频| 色视频www国产| 国产一区亚洲一区在线观看| 亚洲最大成人中文| 熟女人妻精品中文字幕| 色综合亚洲欧美另类图片| 午夜免费男女啪啪视频观看| 99riav亚洲国产免费| 国产精品福利在线免费观看| 国产在线精品亚洲第一网站| 两个人视频免费观看高清| 国产爱豆传媒在线观看| 久久久久久久久久久丰满| 我要看日韩黄色一级片| 成年女人看的毛片在线观看| 国产精品伦人一区二区| 国产成人福利小说| 变态另类成人亚洲欧美熟女| 成人性生交大片免费视频hd| eeuss影院久久| 亚洲欧美日韩高清在线视频| 国产视频首页在线观看| 少妇熟女aⅴ在线视频| 久久午夜福利片| 亚洲电影在线观看av| av女优亚洲男人天堂| 亚洲精品自拍成人| 一级av片app| 九九久久精品国产亚洲av麻豆| 国产色爽女视频免费观看| 国产精品日韩av在线免费观看| 国产午夜精品论理片| 日韩av在线大香蕉| 日本黄色视频三级网站网址| 国产精品一区www在线观看| 2021天堂中文幕一二区在线观| 亚洲美女搞黄在线观看| 国产亚洲91精品色在线| 亚洲欧美成人精品一区二区| 成人午夜高清在线视频| 九九在线视频观看精品| 国产黄色视频一区二区在线观看 | 菩萨蛮人人尽说江南好唐韦庄 | 午夜精品国产一区二区电影 | 国产一区二区三区在线臀色熟女| 伊人久久精品亚洲午夜| 一级毛片我不卡| 欧美一区二区精品小视频在线| a级毛色黄片| 国产高清激情床上av| 国产成人a区在线观看| 成人漫画全彩无遮挡| 三级经典国产精品| 国产老妇伦熟女老妇高清| 亚洲真实伦在线观看| 九九爱精品视频在线观看| 国产亚洲91精品色在线| 色5月婷婷丁香| 好男人视频免费观看在线| 日韩欧美 国产精品| 亚洲一区高清亚洲精品| 美女国产视频在线观看| 国产激情偷乱视频一区二区|