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

    靜壓孔位置對油膜支承可傾瓦軸承承載性能的影響分析

    2022-05-31 04:27:50華爾天趙天城謝榮盛許永利徐高歡
    中國機械工程 2022年10期
    關(guān)鍵詞:瓦塊偏心率軸瓦

    華爾天 趙天城 謝榮盛 許永利 項 春 徐高歡

    1.浙江工業(yè)大學(xué)先進(jìn)制造研究所,杭州,3100142.先進(jìn)水利裝備浙江省工程研究中心,杭州,3100183.浙江水利水電學(xué)院機械與汽車工程學(xué)院,杭州,310018

    0 引言

    可傾瓦徑向軸承作為轉(zhuǎn)子系統(tǒng)的支承,其潤滑油膜不僅會產(chǎn)生油膜承載力,而且能夠減少摩擦和振動,由此被大量用在汽輪機等高速重載的機械設(shè)備中[1-3]。油膜支承可傾瓦徑向軸承將傳統(tǒng)的可傾瓦徑向軸承的機械支點改進(jìn)為柔性支點,利用流體來支撐軸瓦,具有更好的穩(wěn)定性和抗振性[4-6]。

    HOLLINGSWORTH[7]于20世紀(jì)70年代提出了一種新型結(jié)構(gòu)的滑動軸承,即油膜支承可傾瓦徑向軸承;NELSON等[8]論述了油膜支承可傾瓦徑向軸承的工作原理,并搭建實驗裝置開展實驗,研究結(jié)果表明靜壓孔的存在不會削弱軸承承載能力,空載下該軸承的剛度和阻尼分別是機械支承可傾瓦軸承相應(yīng)值的1.7倍和2.3倍,而功耗僅為其40%~60%。HARANGOZO等[9]通過建立雙質(zhì)量彈簧阻尼動力學(xué)模型討論了油膜支承可傾瓦徑向滑動軸承的動態(tài)特性,發(fā)現(xiàn)在一定頻率范圍內(nèi),相比普通滑動軸承,油膜支承可傾瓦徑向滑動軸承能顯著改善轉(zhuǎn)子動力學(xué)性能。羅建彬[10]在已有實驗臺的基礎(chǔ)上搭建測控系統(tǒng),探索油膜支承可傾瓦徑向軸承動靜態(tài)性能,并對其軸瓦的靜壓腔坡度、靜壓孔直徑等幾何參數(shù)進(jìn)行了9組實驗,發(fā)現(xiàn)當(dāng)靜壓腔坡度和靜壓孔孔徑分別為2°和6 mm時,其底部壓力和最高溫度達(dá)到相對最優(yōu)。李偉光等[11]提出了一種前緣溝槽型(leading edge groove,LEG)三瓦可傾瓦軸承,通過對其動態(tài)特性參數(shù)進(jìn)行理論計算和實驗研究,發(fā)現(xiàn)在承載方向有較高的剛度和阻尼。楊期江等[12-13]、趙偉[14]對常規(guī)滑動軸承與油膜支承可傾瓦徑向軸承進(jìn)行了仿真計算和動平衡實驗,發(fā)現(xiàn)重載低速下油膜支承可傾瓦徑向軸承的瓦塊更易浮起,相比普通滑動軸承有更好的減振效果。LOU等[15]發(fā)現(xiàn)油膜支承可傾瓦徑向軸承在單向浮動狀態(tài)下,其瓦塊擺角隨靜壓孔孔徑的增大而增大,其承載能力也有所提高,而靜壓腔面積比對承載能力的影響不大。

    目前對油膜支承可傾瓦軸承進(jìn)行的承載性能分析僅考慮靜壓孔的孔徑等參數(shù)的影響,往往為了適應(yīng)各種結(jié)構(gòu)工況而將孔直接設(shè)定在軸瓦正中間;孔的位置作為關(guān)鍵影響參數(shù),對軸承承載性能影響的研究較少。此外,以往對油膜支承可傾瓦徑向軸承的研究主要通過求解二維Reynolds方程來分析其特性,這種方法忽略了慣性項、油膜曲率、徑向流場變化等因素的影響[16-17],而油膜支承可傾瓦徑向軸承表現(xiàn)出明顯的三維流動特性,因此,有必要直接求解軸承的三維N-S方程,更加準(zhǔn)確地分析其內(nèi)在性能。

    本文基于三維N-S方程,應(yīng)用計算流體動力學(xué)(computational fluid dynamics,CFD)求解油膜支承可傾瓦軸承流場,通過構(gòu)建油膜動態(tài)平衡模型,在指定載荷下求解出平衡位置,獲得相關(guān)性能參數(shù)。在此基礎(chǔ)上預(yù)測靜壓孔的最優(yōu)位置區(qū)域并進(jìn)行優(yōu)化設(shè)計,以偏心率、油膜承載剛度等作為承載性能評價指標(biāo),得到靜壓孔最優(yōu)位置區(qū)域。

    1 平衡方程及剛度計算方法

    1.1 物理模型

    圖1為油膜支承可傾瓦軸承結(jié)構(gòu)示意圖,可知,該軸承由3片軸瓦以及外殼構(gòu)成。在外載荷的作用下,轉(zhuǎn)子產(chǎn)生偏心,潤滑油進(jìn)入并充滿轉(zhuǎn)子與軸瓦的間隙中。由于靜壓孔的存在,部分潤滑油流入軸瓦背面的靜壓腔。當(dāng)靜壓腔中的靜壓力足夠大時,軸瓦浮起,軸瓦背面形成外層油膜。

    圖1 油膜支承可傾瓦軸承結(jié)構(gòu)示意圖Fig.1 Structural diagram of the fluid pivot tilting pad bearing

    1.2 穩(wěn)態(tài)平衡方程

    轉(zhuǎn)子運動模型可簡化為運轉(zhuǎn)時動壓油膜支承旋轉(zhuǎn)的軸頸[18-19],取Oxy平面,如圖2所示。軸頸分別受到三瓦塊的動壓油膜力Fd1、Fd2、Fd3。圖2為瓦1支承軸頸的示意圖。

    圖2 轉(zhuǎn)子運動模型Fig.2 Rotor kinematic model

    穩(wěn)定運轉(zhuǎn)時軸頸受力保持平衡,滿足以下平衡方程組:

    (1)

    式中,W為載荷;Fxi、Fyi分別為Fdi在x、y方向上的油膜力分量。

    以軸瓦1為例,F(xiàn)d1滿足下式:

    (2)

    (3)

    式中,R為軸頸半徑;z為軸向坐標(biāo);φs、φe分別為軸瓦1的起始和終止角度;pd1為油膜動壓力;B為軸承的寬度。

    油膜支承可傾瓦軸承達(dá)到穩(wěn)定運轉(zhuǎn)狀態(tài)時,除軸平穩(wěn)轉(zhuǎn)動外,還需滿足各瓦塊的受力平衡和扭矩平衡。對單瓦進(jìn)行受力分析,如圖3所示,軸瓦受到動壓油膜的壓力Fd、靜壓油膜的壓力Fs和軸承座對它的支承力Fn。

    圖3 軸瓦受力模型Fig.3 Model of bearing bush stress

    平穩(wěn)運轉(zhuǎn)時,瓦塊上下表面受力滿足平衡方程:

    Fd+Fs+Fn=0

    (4)

    當(dāng)瓦塊為擺動狀態(tài)時,F(xiàn)n大于0,動壓力Fd與靜壓力Fs的壓力差由支撐力Fn補償;當(dāng)瓦塊為浮動狀態(tài)時,F(xiàn)n等于0,軸瓦在Fd與Fs作用下保持平衡。瓦塊需滿足扭矩平衡方程:M=0。M為總力矩,可表示為

    M=Md+Ms

    (5)

    式中,Md為動壓油膜力Fd所產(chǎn)生的合力矩;Ms為靜壓油膜力Fs所產(chǎn)生的合力矩。

    1.3 剛度計算方法

    在小擾動條件下,軸承油膜剛度可采用差分法求解。根據(jù)攝動理論,可將油膜力線性化處理[20-21],剛度的計算公式為

    (6)

    式中,ΔFdij為位移擾動引起油膜力的變化,i、j為x、y中的某一個。

    為了簡化計算,給定軸心x方向一微小位移+Δx和-Δx,迭代計算軸心偏離平衡位置后的油膜壓力分布,求出Fx1、Fx2、Fy1、Fy2??汕蟪鲆唤M剛度系數(shù):

    (7)

    同理,給定軸心y方向一微小位移+Δy和-Δy,迭代計算軸心偏離平衡位置后的油膜壓力分布,求出Fx3、Fx4、Fy3、Fy4??汕蟪鲆唤M剛度系數(shù)[21]:

    (8)

    2 油膜模型的建立

    2.1 油膜模型及網(wǎng)格劃分

    油膜支承可傾瓦軸承的結(jié)構(gòu)參數(shù)及工況見表1[8]。根據(jù)軸承的結(jié)構(gòu)參數(shù),利用三維建模軟件建立軸承油膜模型,如圖4所示。

    表1 軸承的結(jié)構(gòu)參數(shù)及工況Tab.1 Structural parameters and working conditions of bearings

    圖4 油膜支承可傾瓦軸承油膜模型Fig.4 Oil film model of the fluid pivot tilting pad bearing

    如圖5所示,采用ICEM CFD進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分。由于軸承結(jié)構(gòu)相對復(fù)雜以及油膜間隙過小,故仿真計算結(jié)果對網(wǎng)格質(zhì)量非常敏感。本文采用分塊劃分網(wǎng)格的方式,網(wǎng)格劃分主要考慮兩方面的影響:一方面,要保證網(wǎng)格的密度,特別是關(guān)鍵區(qū)域的網(wǎng)格,以確保計算的精度;另一方面,需要保證網(wǎng)格的質(zhì)量要求[22-23]。考慮到油膜厚度方向的網(wǎng)格質(zhì)量對計算結(jié)果影響較大,對該區(qū)域的網(wǎng)格采用等分形式劃分并進(jìn)行加密。在890 N定載荷下進(jìn)行網(wǎng)格無關(guān)性分析:當(dāng)厚度方向節(jié)點數(shù)達(dá)到12以上時,偏心率和油膜最大壓力值基本無變化,因此,油膜厚度方向采用12個節(jié)點為宜。

    (a)油膜模型網(wǎng)格劃分

    (b)油膜厚度方向及連接處網(wǎng)格劃分 (c)靜壓孔區(qū)域網(wǎng)格劃分圖5 油膜網(wǎng)格劃分結(jié)果Fig.5 Results of the oil film meshing

    2.2 流態(tài)分析

    在滑動軸承理論中,通常用雷諾數(shù)來表征潤滑油膜的流動狀態(tài):

    Re=Uρhm/μ≈152

    (9)

    式中,U為軸頸轉(zhuǎn)速,m/s;μ為流體的動力黏度,Pa·s;hm為流體膜厚均值,m;ρ為流體密度,kg/m3。

    徑向滑動軸承的臨界雷諾數(shù)為[16]

    (10)

    式中,r為軸頸半徑;c為軸頸半徑間隙。

    經(jīng)計算Re?Rec,軸承內(nèi)部為層流狀態(tài)。因此,采用層流模型對油膜模型進(jìn)行求解。由上,建立以下流體基本方程:

    (11)

    (12)

    式中,ρ為密度;t為時間;u、v、w為速度矢量在x、y、z方向的分量;p為流體微元體上的壓力;μ為動力黏度。

    2.3 邊界條件及參數(shù)設(shè)置

    進(jìn)口邊界為壓力入口,設(shè)置總壓為0.1 MPa;設(shè)置軸承軸向兩端的開敞式邊界為Opening,設(shè)置相對壓力為0;油膜內(nèi)壁面設(shè)為旋轉(zhuǎn)邊界,采用無滑移邊界條件,其他壁面設(shè)為固定邊界。設(shè)置最大迭代步數(shù)為400,殘差收斂精度設(shè)為10-5。

    2.4 計算流程

    因油膜厚度與油膜半徑比值較大,在三維流場仿真中采用有限元體積法求解時,網(wǎng)格相對變形較大,動態(tài)計算不能保證計算收斂。為了保證網(wǎng)格質(zhì)量和計算收斂,應(yīng)用MATLAB對網(wǎng)格節(jié)點進(jìn)行了變換。因此,本文數(shù)值求解時將軸承工作載荷和轉(zhuǎn)速作為初始條件,并考慮軸承的瓦塊擺角,通過MATLAB調(diào)用外部流體求解軟件CFX,實現(xiàn)數(shù)值計算。

    求解得到初始位置的軸頸受力和瓦塊力矩后,將其導(dǎo)入MATLAB中,利用牛頓迭代法預(yù)估下一軸心位置和瓦塊擺角,并調(diào)用CFX獲得N-S方程的數(shù)值解,如此循環(huán)迭代得到平衡狀態(tài)下的軸心位移、瓦塊擺角及油膜壓力分布等。程序流程圖見圖6。

    圖6 平衡位置計算流程圖Fig.6 Flow chat of the equilibrium position

    3 計算方法驗證及分析

    3.1 方法驗證

    為驗證本文計算方法的可靠性,在與文獻(xiàn)[8]相同的結(jié)構(gòu)參數(shù)條件下,分別在890 N、1335 N的載荷下尋找軸承平衡位置并計算偏心率、偏位角等,并與文獻(xiàn)[8]實驗數(shù)據(jù)進(jìn)行比對,具體見表2、表3。

    由上可知,本文計算結(jié)果與文獻(xiàn)[8]實驗結(jié)果的誤差均在12%以內(nèi),且總體上比文獻(xiàn)[8]的計算誤差更小,因此,本文計算方法的準(zhǔn)確性更高。

    3.2 軸承油膜壓力分析

    在該結(jié)構(gòu)下,計算穩(wěn)態(tài)時軸瓦上的油膜壓力

    表2 文獻(xiàn)[8]與本文數(shù)據(jù)比較(載荷890 N)Tab.2 Comparison between reference[8] and the data inthis paper(load 890 N)

    表3 文獻(xiàn)[8]與本文數(shù)據(jù)比較(載荷1335 N)Tab.3 Comparison between reference[8] and the datain this paper(load 1335 N)

    分布,其壓力云圖見圖7??芍?,軸瓦1最易浮起,且每塊瓦的油膜分布有兩個壓力集中區(qū),一個呈現(xiàn)正壓,另一個呈現(xiàn)負(fù)壓,在最小油膜厚度區(qū)被分隔開來[24]。由于軸的偏心和轉(zhuǎn)動,在瓦1正壓區(qū)產(chǎn)生了明顯的收斂楔形,從而產(chǎn)生了較大的正壓力。由圖7b可知,靜壓孔的存在使部分高壓油流向軸瓦背面的靜壓區(qū),使得瓦1較易浮起。而當(dāng)靜壓孔越接近油膜最大壓力區(qū)時,由于孔兩端壓力差的存在,動壓層潤滑油越容易通過靜壓孔流向軸瓦背面,使得軸瓦擺角增大;軸與軸瓦的相對油膜厚度減小,油膜剛度增大。

    (a)軸

    (b)瓦1圖7 油膜壓力分布云圖Fig.7 Distributions of oil film pressure

    基于以上分析,若能改變靜壓孔位置,使其接近油膜最大壓力區(qū),即能提高軸承油膜剛度。

    4 靜壓孔位置的影響

    為了探究靜油孔位置對軸承承載性能的影響,在油膜模型建立時將靜壓腔與軸瓦分離,以便修改模型。據(jù)此,改變靜油腔的旋轉(zhuǎn)角度(-15°~15°,間隔為5°)來實現(xiàn)靜壓孔相對位置γ的變化,其他結(jié)構(gòu)與原模型相同,如圖8所示。

    圖8 靜壓孔相對位置示意圖Fig.8 Diagram of the relative position of static pressure hole

    4.1 對油膜壓力分布的影響

    對于不同靜壓孔位置下的油膜支承可傾瓦軸承油膜模型,在890 N的載荷下尋找其平衡位置,得到其油膜壓力分布云圖,如圖9所示。

    由圖9a可知,當(dāng)γ為-15°、-10°時,油膜壓力分布差別不大,最大油膜壓力值相近。這是因為孔的位置接近于壓力分布的負(fù)壓區(qū),軸瓦內(nèi)外面無法形成壓力差,靜壓孔的存在不能使動壓油流向軸瓦背面,從而使各瓦塊無法浮起,其性能與

    (a)γ=-15°、-10° (b)γ=-5° (c)γ=5°

    (d)γ=10° (e)γ=15°圖9 不同靜壓孔相對位置下的油膜壓力分布云圖Fig.9 Oil film pressure distribution nephogram of the different position of static pressure hole

    普通滑動軸承類似。

    由圖9b、圖9c可知,當(dāng)γ取-5°~5°時,油膜最大壓力隨著γ的增大而增大。孔的位置逐步靠近瓦1的正壓區(qū),此時靜壓孔兩端形成壓力差,動壓油通過靜壓孔流向瓦1背面。當(dāng)瓦1背面的靜壓油達(dá)到一定壓力時,將瓦塊浮起,且越接近瓦1的正壓區(qū),瓦塊越易浮起。

    由圖9d、圖9e可知,當(dāng)γ取10°、15°時,瓦塊均浮起。圖中瓦1、瓦2的壓力分布顯示,瓦1孔的位置在油膜最大壓力區(qū)附近,此時瓦塊受到來自軸瓦背面的油膜力較大,且越接近油膜最大壓力區(qū),瓦1擺角越大,靜壓油膜支承效果越明顯;瓦2孔的位置位于正壓區(qū),在靜壓油膜力的作用下將瓦塊浮起,且越接近瓦2的正壓區(qū),瓦塊越易浮起。

    4.2 對承載性能的影響

    為了便于分析靜壓孔相對位置對承載性能的影響規(guī)律,將最小油膜厚度和偏心率隨靜壓孔相對位置γ的變化整理成圖10,油膜承載剛度隨靜壓孔相對位置γ的變化整理成圖11。

    圖10 靜壓孔相對位置對最小油膜厚度、偏心率的影響Fig.10 Influence of relative position of static pressure hole on minimum oil film thickness and eccentricity

    圖11 靜壓孔相對位置對剛度的影響Fig.11 Influence of relative position of static pressure hole on stiffness

    由圖10可看出,當(dāng)γ取-15°~-10°時,最小膜厚和偏心率變化不大。由于此時瓦塊均無浮起,軸承內(nèi)只存在動壓油膜層,故均無明顯變化。

    當(dāng)γ取-10°~5°時,最小膜厚先增大后減小,偏心率逐漸減小。當(dāng)靜壓孔相對位置γ逐漸增大時,作為主承載的瓦1逐漸浮起,其擺角的增加使最小膜厚減小,承載能力增加,即在載荷下軸心的偏移量減小,其偏心率減小。

    當(dāng)γ取5°~15°時,最小膜厚繼續(xù)減小,偏心率逐漸增大。當(dāng)γ>5°時,瓦塊2、3均浮起,使得軸的相對位置對比單瓦浮起的狀態(tài)有所不同。而γ越大,靜壓孔越接近油膜最大壓力區(qū),瓦塊更容易浮起;瓦塊擺角的增加使軸與瓦塊間的間隙減小,最小膜厚逐漸減小。

    由圖9、圖11可以看出,隨著靜壓孔相對位置γ的增大,孔位置逐漸接近油膜最大壓力分布區(qū),載荷方向的主剛度kyy與非載荷方向的主剛度kxx均逐漸增大,交叉剛度kyx與kxy呈波動上升/下降的趨勢;當(dāng)γ=-15°時,此時瓦塊均無浮起,kyy、kxx分別為3.66×107N/m和2.18×107N/m;當(dāng)γ=10°或15°時,kyy、kxx都較大,其中當(dāng)γ=15°時,kyy、kxx均為最大,分別為1.23×108N/m和5.58×107N/m;但對比圖10可知,此結(jié)構(gòu)下的最小膜厚比原結(jié)構(gòu)的最小膜厚減少了1/3,偏心率接近甚至大于原結(jié)構(gòu)下的偏心率。

    相比之下,在載荷890 N下,當(dāng)靜壓孔相對位置γ=5°時,與原結(jié)構(gòu)(γ=0°)相比,最小膜厚減小了約9.8%,偏心率減小約48%;而主剛度kyy接近原結(jié)構(gòu)的1.4倍,kxx接近原結(jié)構(gòu)的1.1倍。由此可見,軸承在靜壓孔相對位置γ=5°左右時為最優(yōu)位置區(qū)域,能有效提高承載性能,可靠性較好。

    5 試驗驗證

    相關(guān)企業(yè)根據(jù)本文得到的關(guān)于靜壓孔相對位置對油膜支承可傾瓦軸承的性能影響規(guī)律,設(shè)計制造了新型油膜支承可傾瓦軸承(靜壓孔相對位置γ=5°),如圖12所示。

    圖12 油膜支承可傾瓦軸承實物圖Fig.12 Physical picture of fluid pivot tilting pad bearing

    為驗證該產(chǎn)品軸承的可靠性,在杭州汽輪機公司開展相關(guān)試驗研究。對高轉(zhuǎn)速(5000~7000 r/min)下的普通滑動軸承和油膜支承可傾瓦軸承進(jìn)行了承載性能對比試驗,結(jié)果如圖13所示。可以看出,高轉(zhuǎn)速下油膜支承可傾瓦軸承的相對剛度比普通滑動軸承的相對剛度至少提高25%以上,這說明新型油膜支承可傾瓦軸承具有良好的承載性能。

    圖13 相對剛度試驗結(jié)果Fig.13 Test results of relative stiffness

    圖14 軸振試驗結(jié)果Fig.14 Test results of shaft vibration

    油膜支承可傾瓦軸承的減振特性也是綜合性能評定的重要指標(biāo)之一[13]。為了檢驗新型油膜支承可傾瓦軸承的減振性能,采用T5136機組在32T動平衡室加做普通滑動軸承和油膜支承可傾瓦軸承的不平衡響應(yīng)試驗。轉(zhuǎn)子系統(tǒng)中的軸振試驗結(jié)果如圖14所示。可以看出,在一階臨界轉(zhuǎn)速(3400 r/min)區(qū)域內(nèi),普通滑動軸承的振動峰峰值大于油膜支承可傾瓦軸承的振動峰峰值,但過了一階臨界區(qū)域后,油膜支承可傾瓦軸承的振動幅值總體上要比普通滑動軸承的振動幅值低35%以上,這說明新型油膜支承可傾瓦軸承具有良好的減振性能。

    綜上所述,高轉(zhuǎn)速下新型油膜支承可傾瓦軸承相比普通滑動軸承有更好的綜合性能。

    6 結(jié)論

    (1)軸承靜壓孔的存在使運轉(zhuǎn)中部分高壓油流向軸瓦背面的靜壓區(qū),主承載瓦較易擺起。

    (2)原模型改變了軸承靜壓孔相對位置,得到了不同位置下的有關(guān)性能參數(shù)。載荷為890 N時,當(dāng)靜壓孔相對位置γ=5°左右時,最小膜厚約減小了9.8%,偏心率下降約48%,主剛度kyy、kxx接近原結(jié)構(gòu)的1.4倍和1.1倍,此時孔位置為最優(yōu)位置區(qū)域。

    (3)在油膜支承可傾瓦軸承靜壓孔參數(shù)設(shè)計中,可根據(jù)其油膜壓力分布規(guī)律設(shè)計孔位置,使其接近最大壓力分布區(qū),以提高軸承承載性能。

    (4)試驗驗證發(fā)現(xiàn),優(yōu)化后的新型油膜支承可傾瓦軸承具有良好的綜合性能。

    猜你喜歡
    瓦塊偏心率軸瓦
    Hansen系數(shù)遞推的效率?
    一種高效的頂點偏心率計算方法
    2號發(fā)電機#3軸瓦軸振超標(biāo)原因分析及處理
    某低速機主軸瓦連接螺栓預(yù)緊力計算
    無縫鋼管壁厚偏心率的測量分析及降低方法
    鋼管(2016年1期)2016-05-17 06:12:44
    基于幾何條件的可傾瓦軸承油膜邊界條件判定方法
    測量計算法測定軸瓦間隙
    大偏心率軌道星上快速計算方法
    道奇軸瓦損壞與檢修
    汽輪機可傾瓦顫振引發(fā)低頻振蕩問題研究
    制服人妻中文乱码| 久久这里只有精品中国| 麻豆国产97在线/欧美| 夜夜爽天天搞| 成人亚洲精品av一区二区| 亚洲午夜理论影院| 午夜视频精品福利| 99热只有精品国产| 亚洲欧美日韩卡通动漫| 色av中文字幕| 一个人看视频在线观看www免费 | 制服丝袜大香蕉在线| 黄色片一级片一级黄色片| 国产成人av激情在线播放| 亚洲欧美精品综合一区二区三区| 国产精品一及| 国产亚洲精品久久久com| 特大巨黑吊av在线直播| 欧美日韩瑟瑟在线播放| 啦啦啦免费观看视频1| 久久久久久久午夜电影| 久久午夜亚洲精品久久| 欧美最黄视频在线播放免费| 全区人妻精品视频| 午夜a级毛片| 性色av乱码一区二区三区2| 国产成人一区二区三区免费视频网站| 人人妻人人看人人澡| 一区二区三区激情视频| 一二三四社区在线视频社区8| a在线观看视频网站| 午夜福利成人在线免费观看| 网址你懂的国产日韩在线| 无遮挡黄片免费观看| 亚洲无线在线观看| 精品一区二区三区四区五区乱码| 一本综合久久免费| 国产成人系列免费观看| 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| 在线a可以看的网站| 精品国产三级普通话版| 日本一二三区视频观看| 成人欧美大片| 午夜免费成人在线视频| 亚洲欧美一区二区三区黑人| 午夜成年电影在线免费观看| 全区人妻精品视频| 女生性感内裤真人,穿戴方法视频| 禁无遮挡网站| 99精品在免费线老司机午夜| 一区二区三区高清视频在线| 久久精品综合一区二区三区| 最新在线观看一区二区三区| 国内精品一区二区在线观看| 婷婷六月久久综合丁香| 久久久久国内视频| 啦啦啦韩国在线观看视频| 国产一区二区三区视频了| 久久久久久久久久黄片| 99精品久久久久人妻精品| 在线看三级毛片| 精品国内亚洲2022精品成人| 我的老师免费观看完整版| 亚洲国产欧洲综合997久久,| 99热6这里只有精品| 一级毛片高清免费大全| 亚洲国产精品999在线| 午夜亚洲福利在线播放| 亚洲第一电影网av| 999精品在线视频| 在线免费观看的www视频| 69av精品久久久久久| 精品99又大又爽又粗少妇毛片 | 老鸭窝网址在线观看| 亚洲国产精品久久男人天堂| 亚洲中文字幕一区二区三区有码在线看 | 成年免费大片在线观看| 久久久久免费精品人妻一区二区| 亚洲国产精品sss在线观看| 国产黄色小视频在线观看| 精品久久蜜臀av无| 亚洲成人免费电影在线观看| 国产精品九九99| 日韩有码中文字幕| 特级一级黄色大片| 最新中文字幕久久久久 | 一本综合久久免费| 日韩 欧美 亚洲 中文字幕| 亚洲av熟女| 国产成人啪精品午夜网站| 黄色日韩在线| 国产黄a三级三级三级人| 操出白浆在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 丰满的人妻完整版| 国产麻豆成人av免费视频| 亚洲中文日韩欧美视频| 亚洲美女黄片视频| 欧美zozozo另类| 日本三级黄在线观看| 99视频精品全部免费 在线 | 免费无遮挡裸体视频| 亚洲美女视频黄频| 欧美日韩一级在线毛片| 熟妇人妻久久中文字幕3abv| 18禁国产床啪视频网站| 久9热在线精品视频| 国产欧美日韩精品一区二区| 午夜福利在线在线| 久久亚洲真实| 99热这里只有是精品50| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 一本精品99久久精品77| 婷婷精品国产亚洲av在线| 偷拍熟女少妇极品色| 欧美丝袜亚洲另类 | 午夜免费激情av| 精品国产乱子伦一区二区三区| 欧美一区二区精品小视频在线| 久久久久久大精品| 成人鲁丝片一二三区免费| 在线十欧美十亚洲十日本专区| 国产真人三级小视频在线观看| 欧美黑人巨大hd| 国产高清三级在线| 精品久久久久久,| 夜夜爽天天搞| 夜夜躁狠狠躁天天躁| 日韩欧美 国产精品| 综合色av麻豆| 成人18禁在线播放| 人妻夜夜爽99麻豆av| 国产三级黄色录像| 女人被狂操c到高潮| 国产精品98久久久久久宅男小说| 宅男免费午夜| 午夜精品久久久久久毛片777| 看片在线看免费视频| 亚洲 国产 在线| 在线观看美女被高潮喷水网站 | 国内毛片毛片毛片毛片毛片| 日本在线视频免费播放| 国产爱豆传媒在线观看| 久久久久久久精品吃奶| 脱女人内裤的视频| 一级a爱片免费观看的视频| 免费搜索国产男女视频| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 天天添夜夜摸| 精品久久久久久成人av| 欧美三级亚洲精品| 香蕉丝袜av| 亚洲国产精品久久男人天堂| 又粗又爽又猛毛片免费看| 99国产综合亚洲精品| 成人三级做爰电影| 国产极品精品免费视频能看的| 国产黄a三级三级三级人| 桃色一区二区三区在线观看| 夜夜看夜夜爽夜夜摸| 亚洲熟妇熟女久久| 老司机福利观看| 小说图片视频综合网站| 中文字幕熟女人妻在线| 午夜久久久久精精品| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 国内久久婷婷六月综合欲色啪| 不卡一级毛片| 午夜亚洲福利在线播放| 亚洲欧洲精品一区二区精品久久久| ponron亚洲| 在线观看美女被高潮喷水网站 | 麻豆av在线久日| 人妻丰满熟妇av一区二区三区| 此物有八面人人有两片| 18美女黄网站色大片免费观看| 免费在线观看日本一区| 精华霜和精华液先用哪个| 亚洲五月天丁香| 午夜精品一区二区三区免费看| 亚洲 国产 在线| 一夜夜www| 午夜成年电影在线免费观看| 亚洲 国产 在线| 91在线精品国自产拍蜜月 | 九九在线视频观看精品| 亚洲九九香蕉| 中文字幕最新亚洲高清| 国产三级中文精品| 欧美成人一区二区免费高清观看 | 99国产综合亚洲精品| 亚洲国产精品久久男人天堂| 日本五十路高清| 精品国产乱码久久久久久男人| 国产精品久久电影中文字幕| 亚洲专区中文字幕在线| 国产精品久久视频播放| 一级毛片精品| 巨乳人妻的诱惑在线观看| 成年人黄色毛片网站| 午夜日韩欧美国产| 一进一出好大好爽视频| 狂野欧美白嫩少妇大欣赏| 免费看光身美女| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 国产高潮美女av| 日本a在线网址| 最新美女视频免费是黄的| ponron亚洲| 国产高清有码在线观看视频| 免费在线观看日本一区| 国产免费av片在线观看野外av| 真人一进一出gif抽搐免费| 好看av亚洲va欧美ⅴa在| 亚洲激情在线av| 在线观看免费午夜福利视频| 国产黄片美女视频| 变态另类丝袜制服| 欧美中文综合在线视频| 我的老师免费观看完整版| 国产黄色小视频在线观看| 欧美一级a爱片免费观看看| 2021天堂中文幕一二区在线观| 99久久无色码亚洲精品果冻| 99久久99久久久精品蜜桃| 一a级毛片在线观看| 99精品久久久久人妻精品| 亚洲av五月六月丁香网| 舔av片在线| 级片在线观看| 美女免费视频网站| 桃色一区二区三区在线观看| 久久久久久人人人人人| 国产亚洲av嫩草精品影院| 久久精品国产综合久久久| 精品国产美女av久久久久小说| 99久久国产精品久久久| 日本黄大片高清| 91在线精品国自产拍蜜月 | 国产高潮美女av| 色播亚洲综合网| 国产欧美日韩精品一区二区| 精品国产美女av久久久久小说| 无人区码免费观看不卡| 九九在线视频观看精品| 欧美激情在线99| 色综合婷婷激情| 国产91精品成人一区二区三区| 中文字幕精品亚洲无线码一区| 欧美日韩精品网址| 最近视频中文字幕2019在线8| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 十八禁人妻一区二区| 又黄又粗又硬又大视频| 天天添夜夜摸| 中文字幕av在线有码专区| 国产精品,欧美在线| 国产视频内射| 亚洲中文字幕一区二区三区有码在线看 | 男女那种视频在线观看| 欧美日本亚洲视频在线播放| bbb黄色大片| 亚洲熟女毛片儿| 国产精品精品国产色婷婷| 日韩欧美国产在线观看| 巨乳人妻的诱惑在线观看| 亚洲国产色片| 午夜福利在线观看免费完整高清在 | 国产男靠女视频免费网站| 黄片大片在线免费观看| 色老头精品视频在线观看| 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 国产三级在线视频| 丝袜人妻中文字幕| 国产av麻豆久久久久久久| 免费av毛片视频| 一本一本综合久久| 欧美高清成人免费视频www| 日韩欧美三级三区| 国产 一区 欧美 日韩| 亚洲欧美日韩无卡精品| 久久久色成人| 99视频精品全部免费 在线 | 18禁观看日本| 午夜福利视频1000在线观看| 国产黄片美女视频| 夜夜爽天天搞| 亚洲最大成人中文| 嫩草影院精品99| 国产人伦9x9x在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产日本99.免费观看| 黄色 视频免费看| 国模一区二区三区四区视频 | 日韩高清综合在线| 婷婷亚洲欧美| 狠狠狠狠99中文字幕| 久久久久久久久中文| 岛国视频午夜一区免费看| 精品日产1卡2卡| 亚洲美女黄片视频| 亚洲成a人片在线一区二区| 99久久久亚洲精品蜜臀av| 久久久久九九精品影院| 国产精品一区二区精品视频观看| 成人午夜高清在线视频| 搡老熟女国产l中国老女人| 变态另类丝袜制服| 国产精品电影一区二区三区| 国产蜜桃级精品一区二区三区| 婷婷丁香在线五月| 国产真人三级小视频在线观看| 夜夜夜夜夜久久久久| 国内少妇人妻偷人精品xxx网站 | 啦啦啦免费观看视频1| 欧美色欧美亚洲另类二区| 国产乱人视频| 午夜激情福利司机影院| 日韩人妻高清精品专区| 久99久视频精品免费| 动漫黄色视频在线观看| 久久久久九九精品影院| 亚洲人成网站在线播放欧美日韩| 丝袜人妻中文字幕| 一二三四在线观看免费中文在| 香蕉丝袜av| 真人做人爱边吃奶动态| 99在线人妻在线中文字幕| 久久久久久久久免费视频了| 在线永久观看黄色视频| 亚洲avbb在线观看| 国产一区二区在线观看日韩 | 日本三级黄在线观看| 国产又色又爽无遮挡免费看| 亚洲国产中文字幕在线视频| 久久精品夜夜夜夜夜久久蜜豆| 99久国产av精品| 国产乱人视频| aaaaa片日本免费| 老熟妇仑乱视频hdxx| www.www免费av| 在线a可以看的网站| 欧美黑人巨大hd| 网址你懂的国产日韩在线| 欧美激情在线99| 757午夜福利合集在线观看| 黄频高清免费视频| 天天一区二区日本电影三级| 狠狠狠狠99中文字幕| 国产成+人综合+亚洲专区| 噜噜噜噜噜久久久久久91| 一级毛片女人18水好多| 中文字幕人妻丝袜一区二区| 成人一区二区视频在线观看| 久久香蕉国产精品| 免费在线观看日本一区| 少妇裸体淫交视频免费看高清| 在线国产一区二区在线| 国产高清视频在线观看网站| 中文字幕久久专区| 国产一区二区激情短视频| 国内精品久久久久精免费| 午夜免费观看网址| 欧美大码av| 国产精品一及| 黑人欧美特级aaaaaa片| 老汉色av国产亚洲站长工具| 一级毛片高清免费大全| 天天添夜夜摸| 免费看光身美女| 国产三级黄色录像| 天堂av国产一区二区熟女人妻| 久久久水蜜桃国产精品网| 啦啦啦韩国在线观看视频| 亚洲av免费在线观看| 99久久成人亚洲精品观看| 99re在线观看精品视频| 欧美极品一区二区三区四区| 久久婷婷人人爽人人干人人爱| 久久久久国产一级毛片高清牌| 日本五十路高清| 在线免费观看不下载黄p国产 | 可以在线观看的亚洲视频| 国产精品综合久久久久久久免费| 狂野欧美白嫩少妇大欣赏| 亚洲 国产 在线| 成人av在线播放网站| 亚洲 欧美 日韩 在线 免费| 男人的好看免费观看在线视频| 麻豆成人av在线观看| a级毛片在线看网站| 国产精品1区2区在线观看.| 亚洲欧美日韩卡通动漫| 成年版毛片免费区| 啪啪无遮挡十八禁网站| 手机成人av网站| av视频在线观看入口| 特级一级黄色大片| 两个人视频免费观看高清| 欧美日韩瑟瑟在线播放| 欧美在线一区亚洲| netflix在线观看网站| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 国产黄a三级三级三级人| 国产探花在线观看一区二区| 久久香蕉国产精品| 国产亚洲精品一区二区www| 免费大片18禁| 曰老女人黄片| 嫩草影院入口| 国产精品亚洲一级av第二区| 亚洲成av人片免费观看| 欧美一区二区精品小视频在线| 久9热在线精品视频| 成人av一区二区三区在线看| 两人在一起打扑克的视频| 国产精品久久久久久亚洲av鲁大| 国产成人欧美在线观看| 久久亚洲精品不卡| 在线播放国产精品三级| 国产高清视频在线观看网站| 看片在线看免费视频| 成人特级av手机在线观看| 黑人操中国人逼视频| 免费观看人在逋| 国产亚洲精品一区二区www| 国产伦人伦偷精品视频| 国产乱人伦免费视频| 亚洲精品久久国产高清桃花| 身体一侧抽搐| 一级作爱视频免费观看| 久久久久亚洲av毛片大全| 欧美黄色淫秽网站| 国产精品香港三级国产av潘金莲| 中文在线观看免费www的网站| 日本精品一区二区三区蜜桃| 免费在线观看日本一区| 嫩草影院精品99| 免费在线观看视频国产中文字幕亚洲| 亚洲精品中文字幕一二三四区| 免费无遮挡裸体视频| 一边摸一边抽搐一进一小说| 亚洲国产高清在线一区二区三| 香蕉国产在线看| 成人av一区二区三区在线看| 国产视频内射| 99re在线观看精品视频| 制服人妻中文乱码| 亚洲无线在线观看| 国产伦在线观看视频一区| 一夜夜www| 在线a可以看的网站| 亚洲成av人片免费观看| 成人国产一区最新在线观看| e午夜精品久久久久久久| 美女午夜性视频免费| 久久午夜综合久久蜜桃| 日韩三级视频一区二区三区| 少妇的逼水好多| 久久久精品欧美日韩精品| 免费看光身美女| 久久天躁狠狠躁夜夜2o2o| 全区人妻精品视频| bbb黄色大片| 无遮挡黄片免费观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品,欧美在线| 天堂av国产一区二区熟女人妻| 免费看美女性在线毛片视频| 日韩精品中文字幕看吧| 十八禁人妻一区二区| 国产精品综合久久久久久久免费| 首页视频小说图片口味搜索| 热99在线观看视频| 午夜精品一区二区三区免费看| 成在线人永久免费视频| 一个人看的www免费观看视频| 青草久久国产| 成人国产一区最新在线观看| 老司机福利观看| 宅男免费午夜| 午夜日韩欧美国产| 免费人成视频x8x8入口观看| 日本黄大片高清| 日本撒尿小便嘘嘘汇集6| 国产又色又爽无遮挡免费看| 精品一区二区三区av网在线观看| 最近视频中文字幕2019在线8| 婷婷亚洲欧美| 亚洲欧美日韩高清专用| 国产野战对白在线观看| 亚洲欧洲精品一区二区精品久久久| 床上黄色一级片| 久久精品国产综合久久久| 久久久久久久精品吃奶| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 97人妻精品一区二区三区麻豆| 精品久久久久久成人av| 欧美黑人欧美精品刺激| 好看av亚洲va欧美ⅴa在| 免费看十八禁软件| 午夜视频精品福利| 熟女少妇亚洲综合色aaa.| 91麻豆精品激情在线观看国产| 免费一级毛片在线播放高清视频| 欧美在线一区亚洲| av在线蜜桃| 桃色一区二区三区在线观看| 91麻豆精品激情在线观看国产| 一级毛片精品| 两个人视频免费观看高清| 国产一区二区三区视频了| 成年免费大片在线观看| 久久精品国产清高在天天线| 成人一区二区视频在线观看| 国产高清有码在线观看视频| 夜夜看夜夜爽夜夜摸| 国产亚洲欧美98| 精品国产美女av久久久久小说| 国产伦精品一区二区三区四那| 国产高清有码在线观看视频| 亚洲一区二区三区色噜噜| 成年免费大片在线观看| 国产蜜桃级精品一区二区三区| 国产激情偷乱视频一区二区| 欧美成人性av电影在线观看| 亚洲电影在线观看av| 国产视频内射| 亚洲av免费在线观看| 国产午夜精品论理片| 免费观看的影片在线观看| 波多野结衣高清无吗| 熟女少妇亚洲综合色aaa.| АⅤ资源中文在线天堂| 91九色精品人成在线观看| 夜夜爽天天搞| 久久久久久九九精品二区国产| 国产精品免费一区二区三区在线| 日本成人三级电影网站| 国产综合懂色| 身体一侧抽搐| 国产蜜桃级精品一区二区三区| av女优亚洲男人天堂 | 久久久久久久久中文| 国产高潮美女av| 免费观看人在逋| 国产人伦9x9x在线观看| 久久精品人妻少妇| 国产高清有码在线观看视频| 99精品在免费线老司机午夜| 亚洲国产欧洲综合997久久,| 一卡2卡三卡四卡精品乱码亚洲| 亚洲真实伦在线观看| 亚洲国产精品成人综合色| 中文字幕av在线有码专区| 久久精品国产清高在天天线| 国产视频内射| 成年女人看的毛片在线观看| 日本 欧美在线| 欧美乱色亚洲激情| 琪琪午夜伦伦电影理论片6080| 亚洲一区高清亚洲精品| 18禁观看日本| 亚洲五月天丁香| 看免费av毛片| 精华霜和精华液先用哪个| av女优亚洲男人天堂 | 日本精品一区二区三区蜜桃| 色视频www国产| 在线观看66精品国产| 国产成年人精品一区二区| 欧美xxxx黑人xx丫x性爽| 成人精品一区二区免费| 久久精品夜夜夜夜夜久久蜜豆| 两性夫妻黄色片| 99国产精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 在线播放国产精品三级| 12—13女人毛片做爰片一| 一a级毛片在线观看| 国产蜜桃级精品一区二区三区| 欧美国产日韩亚洲一区| 午夜a级毛片| 91麻豆av在线| 欧美三级亚洲精品| 真人一进一出gif抽搐免费| 亚洲一区二区三区色噜噜| 99久久无色码亚洲精品果冻| 小蜜桃在线观看免费完整版高清| 中文字幕人妻丝袜一区二区| 欧洲精品卡2卡3卡4卡5卡区| 久久国产乱子伦精品免费另类| 久久这里只有精品19| 亚洲国产高清在线一区二区三| 午夜免费激情av| 亚洲 欧美 日韩 在线 免费| 亚洲色图 男人天堂 中文字幕| 亚洲欧美日韩高清专用| 在线观看舔阴道视频| 少妇裸体淫交视频免费看高清| 欧美乱码精品一区二区三区| a级毛片a级免费在线| 欧美日韩国产亚洲二区| 少妇裸体淫交视频免费看高清| 99国产精品一区二区三区| 国产精品久久久久久人妻精品电影|