• <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ā)低頻振蕩問題研究
    欧美一区二区国产精品久久精品| 中文字幕熟女人妻在线| av在线蜜桃| 欧美日韩国产亚洲二区| 蜜桃久久精品国产亚洲av| a在线观看视频网站| 亚洲国产欧洲综合997久久,| 国产男靠女视频免费网站| 国产精品99久久99久久久不卡| 黄色丝袜av网址大全| 国产中年淑女户外野战色| 女人十人毛片免费观看3o分钟| 国内揄拍国产精品人妻在线| 91久久精品电影网| 亚洲无线观看免费| 成人亚洲精品av一区二区| 久久九九热精品免费| 久久久国产精品麻豆| 九九久久精品国产亚洲av麻豆| 色哟哟哟哟哟哟| 久久精品国产亚洲av涩爱 | 国产av麻豆久久久久久久| av片东京热男人的天堂| 夜夜爽天天搞| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人不卡在线观看播放网| 欧美日韩亚洲国产一区二区在线观看| 美女被艹到高潮喷水动态| 精品一区二区三区视频在线 | 两个人看的免费小视频| 国产97色在线日韩免费| 欧美成人一区二区免费高清观看| 高清日韩中文字幕在线| 久久久久久九九精品二区国产| 日韩国内少妇激情av| 欧美一级a爱片免费观看看| 啦啦啦观看免费观看视频高清| 亚洲aⅴ乱码一区二区在线播放| 国产真人三级小视频在线观看| 欧美日韩国产亚洲二区| 精品久久久久久成人av| 性色avwww在线观看| 99久久精品热视频| 国产亚洲精品久久久久久毛片| 婷婷丁香在线五月| 免费观看人在逋| 精品人妻偷拍中文字幕| 少妇的逼水好多| 18禁裸乳无遮挡免费网站照片| 久久久久久久久大av| 国产一区二区三区视频了| 日本在线视频免费播放| 亚洲五月婷婷丁香| 香蕉av资源在线| 夜夜看夜夜爽夜夜摸| 日本黄色视频三级网站网址| 亚洲欧美日韩东京热| 中文在线观看免费www的网站| 亚洲欧美精品综合久久99| 国产精品美女特级片免费视频播放器| 久久久久国内视频| 51午夜福利影视在线观看| 亚洲国产欧美人成| 亚洲国产欧洲综合997久久,| 色吧在线观看| 国产亚洲av嫩草精品影院| 嫩草影院精品99| 国产国拍精品亚洲av在线观看 | 欧美黄色片欧美黄色片| 搡女人真爽免费视频火全软件 | 国产 一区 欧美 日韩| 国产久久久一区二区三区| 免费看日本二区| 午夜激情福利司机影院| 日本一二三区视频观看| 日韩高清综合在线| 欧美xxxx黑人xx丫x性爽| 久久天躁狠狠躁夜夜2o2o| 国产欧美日韩精品亚洲av| 日韩欧美一区二区三区在线观看| 国内精品美女久久久久久| h日本视频在线播放| 中文亚洲av片在线观看爽| 午夜两性在线视频| 中文亚洲av片在线观看爽| 老熟妇仑乱视频hdxx| 亚洲欧美日韩高清在线视频| 久久久久精品国产欧美久久久| 美女 人体艺术 gogo| 精品欧美国产一区二区三| 久久久久国内视频| 51午夜福利影视在线观看| 国产成人福利小说| 国产精品,欧美在线| 九九久久精品国产亚洲av麻豆| 精品人妻1区二区| 麻豆久久精品国产亚洲av| 一个人看的www免费观看视频| 少妇的逼水好多| 欧美另类亚洲清纯唯美| av女优亚洲男人天堂| 一本久久中文字幕| 欧美一区二区国产精品久久精品| 一区福利在线观看| 一区二区三区高清视频在线| 国产午夜福利久久久久久| 日本成人三级电影网站| 激情在线观看视频在线高清| 特级一级黄色大片| 91九色精品人成在线观看| 看片在线看免费视频| 国产欧美日韩精品一区二区| 天堂网av新在线| 国产 一区 欧美 日韩| 色尼玛亚洲综合影院| 国产精品永久免费网站| 亚洲不卡免费看| 亚洲不卡免费看| svipshipincom国产片| 狂野欧美激情性xxxx| 三级毛片av免费| 国产欧美日韩一区二区精品| 男人的好看免费观看在线视频| 最近最新中文字幕大全免费视频| 精品久久久久久久毛片微露脸| 19禁男女啪啪无遮挡网站| 日韩 欧美 亚洲 中文字幕| 色哟哟哟哟哟哟| 一区二区三区高清视频在线| 91久久精品电影网| 日韩免费av在线播放| 成人国产一区最新在线观看| 悠悠久久av| 精品一区二区三区视频在线观看免费| 欧美日韩综合久久久久久 | 国产免费一级a男人的天堂| 国产极品精品免费视频能看的| 三级男女做爰猛烈吃奶摸视频| 午夜精品一区二区三区免费看| 欧美日本亚洲视频在线播放| 性欧美人与动物交配| 国产一区二区激情短视频| 国产精品亚洲一级av第二区| 女人被狂操c到高潮| 国产精品三级大全| 精品福利观看| tocl精华| 亚洲精品成人久久久久久| 少妇的逼水好多| 18禁裸乳无遮挡免费网站照片| 日韩高清综合在线| 悠悠久久av| 亚洲人成网站在线播| 国产精品日韩av在线免费观看| 天天一区二区日本电影三级| 在线观看免费午夜福利视频| 男女之事视频高清在线观看| 超碰av人人做人人爽久久 | 久久6这里有精品| 十八禁人妻一区二区| 久久九九热精品免费| 一本久久中文字幕| 日韩亚洲欧美综合| x7x7x7水蜜桃| 蜜桃亚洲精品一区二区三区| svipshipincom国产片| 中文在线观看免费www的网站| 久久国产精品人妻蜜桃| 国产成人a区在线观看| 一区二区三区国产精品乱码| 伊人久久精品亚洲午夜| 国产伦一二天堂av在线观看| 少妇高潮的动态图| 成人特级黄色片久久久久久久| 午夜福利在线观看免费完整高清在 | 欧美极品一区二区三区四区| 亚洲avbb在线观看| 成熟少妇高潮喷水视频| 日韩欧美在线二视频| 真实男女啪啪啪动态图| 日韩大尺度精品在线看网址| netflix在线观看网站| 午夜免费男女啪啪视频观看 | 热99re8久久精品国产| 狂野欧美激情性xxxx| 级片在线观看| 天堂av国产一区二区熟女人妻| 啦啦啦观看免费观看视频高清| 午夜精品一区二区三区免费看| 高潮久久久久久久久久久不卡| 国产精品一区二区免费欧美| 丁香六月欧美| 一级毛片高清免费大全| 韩国av一区二区三区四区| 欧美激情久久久久久爽电影| 欧美日韩精品网址| 三级男女做爰猛烈吃奶摸视频| 女同久久另类99精品国产91| 亚洲成人久久性| 免费看a级黄色片| 色噜噜av男人的天堂激情| 美女免费视频网站| 免费av毛片视频| 国产av不卡久久| 精品久久久久久久毛片微露脸| 淫秽高清视频在线观看| 欧美日韩乱码在线| 国产精品久久视频播放| 毛片女人毛片| 日本撒尿小便嘘嘘汇集6| 美女高潮的动态| 色尼玛亚洲综合影院| 亚洲第一电影网av| 日日夜夜操网爽| 精品久久久久久久人妻蜜臀av| 欧美色欧美亚洲另类二区| 亚洲aⅴ乱码一区二区在线播放| 丰满的人妻完整版| 嫩草影视91久久| 中文资源天堂在线| 国产精品永久免费网站| bbb黄色大片| 亚洲性夜色夜夜综合| eeuss影院久久| 12—13女人毛片做爰片一| 少妇的逼水好多| 精品久久久久久久久久久久久| 99精品久久久久人妻精品| 深爱激情五月婷婷| 亚洲av免费在线观看| 叶爱在线成人免费视频播放| 免费人成视频x8x8入口观看| 日韩亚洲欧美综合| 亚洲欧美精品综合久久99| 国产精品一区二区免费欧美| 国产精品亚洲美女久久久| 亚洲电影在线观看av| 国产男靠女视频免费网站| 丰满人妻一区二区三区视频av | 国产成人影院久久av| 性色avwww在线观看| 欧美黄色片欧美黄色片| 少妇的逼好多水| 国产aⅴ精品一区二区三区波| 中文在线观看免费www的网站| a级一级毛片免费在线观看| 网址你懂的国产日韩在线| 日韩精品青青久久久久久| 国产国拍精品亚洲av在线观看 | 老鸭窝网址在线观看| 中国美女看黄片| 成年人黄色毛片网站| 亚洲欧美日韩卡通动漫| 日本一二三区视频观看| 操出白浆在线播放| 美女 人体艺术 gogo| 午夜福利成人在线免费观看| 99久久成人亚洲精品观看| 啦啦啦韩国在线观看视频| 成人永久免费在线观看视频| av天堂中文字幕网| 1000部很黄的大片| 色在线成人网| 免费av毛片视频| 国产亚洲欧美在线一区二区| 亚洲无线在线观看| 日本a在线网址| 老鸭窝网址在线观看| 国产麻豆成人av免费视频| 中文在线观看免费www的网站| 伊人久久精品亚洲午夜| 亚洲中文字幕一区二区三区有码在线看| 久久国产乱子伦精品免费另类| 日韩欧美精品v在线| 国产国拍精品亚洲av在线观看 | 欧美日韩福利视频一区二区| 欧美中文综合在线视频| 午夜视频国产福利| 国内精品美女久久久久久| eeuss影院久久| 国产一区二区三区在线臀色熟女| 最后的刺客免费高清国语| 成熟少妇高潮喷水视频| 欧美乱码精品一区二区三区| 午夜两性在线视频| 国产男靠女视频免费网站| 19禁男女啪啪无遮挡网站| 亚洲精品久久国产高清桃花| 真实男女啪啪啪动态图| 女同久久另类99精品国产91| 又粗又爽又猛毛片免费看| 国产精品免费一区二区三区在线| 国产精品精品国产色婷婷| 一本一本综合久久| 制服人妻中文乱码| 波野结衣二区三区在线 | 男人舔女人下体高潮全视频| 免费看美女性在线毛片视频| 三级毛片av免费| 男女那种视频在线观看| 成人特级黄色片久久久久久久| 亚洲自拍偷在线| or卡值多少钱| 很黄的视频免费| 国产激情偷乱视频一区二区| 亚洲 国产 在线| 免费高清视频大片| 男女视频在线观看网站免费| 夜夜夜夜夜久久久久| 怎么达到女性高潮| 成人特级黄色片久久久久久久| av国产免费在线观看| 欧美一区二区亚洲| 国产成人av教育| 最新美女视频免费是黄的| 黄色日韩在线| 日韩 欧美 亚洲 中文字幕| 国产一区二区三区在线臀色熟女| 给我免费播放毛片高清在线观看| 亚洲男人的天堂狠狠| 少妇熟女aⅴ在线视频| 日本与韩国留学比较| 亚洲av电影在线进入| 少妇人妻精品综合一区二区 | 全区人妻精品视频| 他把我摸到了高潮在线观看| 男女视频在线观看网站免费| 校园春色视频在线观看| 亚洲欧美日韩东京热| 午夜福利18| 国内毛片毛片毛片毛片毛片| 可以在线观看的亚洲视频| 窝窝影院91人妻| 国语自产精品视频在线第100页| 午夜福利18| 欧美一区二区亚洲| 精品午夜福利视频在线观看一区| 日韩欧美在线乱码| 男人舔奶头视频| 欧美性猛交黑人性爽| 午夜亚洲福利在线播放| aaaaa片日本免费| 99久久精品一区二区三区| 12—13女人毛片做爰片一| 久久中文看片网| 男女视频在线观看网站免费| 蜜桃久久精品国产亚洲av| 老熟妇乱子伦视频在线观看| 午夜福利欧美成人| 亚洲av熟女| 熟女少妇亚洲综合色aaa.| 国产精品98久久久久久宅男小说| 亚洲成人久久性| 黄片小视频在线播放| 欧美成人一区二区免费高清观看| 黑人欧美特级aaaaaa片| 神马国产精品三级电影在线观看| 日日干狠狠操夜夜爽| 亚洲人成网站高清观看| 人人妻人人澡欧美一区二区| 老汉色av国产亚洲站长工具| 偷拍熟女少妇极品色| 精品一区二区三区视频在线 | 亚洲 欧美 日韩 在线 免费| 午夜a级毛片| 岛国在线免费视频观看| 国产69精品久久久久777片| 亚洲国产欧美网| 手机成人av网站| 日本与韩国留学比较| 欧美性猛交黑人性爽| 欧美色欧美亚洲另类二区| 欧美不卡视频在线免费观看| 日韩欧美一区二区三区在线观看| 亚洲精华国产精华精| 欧美大码av| 国产精品亚洲一级av第二区| 亚洲av成人av| 午夜精品久久久久久毛片777| 国产又黄又爽又无遮挡在线| 国内精品美女久久久久久| 久久久久久国产a免费观看| 99久久精品国产亚洲精品| 小蜜桃在线观看免费完整版高清| 国产精品综合久久久久久久免费| 一二三四社区在线视频社区8| 一本一本综合久久| 看黄色毛片网站| 欧美av亚洲av综合av国产av| 精品免费久久久久久久清纯| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人免费电影在线观看| 岛国视频午夜一区免费看| 女警被强在线播放| 国产激情欧美一区二区| av天堂在线播放| 免费大片18禁| 内地一区二区视频在线| 精品国产超薄肉色丝袜足j| 精品乱码久久久久久99久播| 禁无遮挡网站| 午夜免费男女啪啪视频观看 | av天堂中文字幕网| 99热只有精品国产| 可以在线观看的亚洲视频| 免费看光身美女| 亚洲激情在线av| 国产精品美女特级片免费视频播放器| 亚洲中文日韩欧美视频| av中文乱码字幕在线| 香蕉丝袜av| 久久久久精品国产欧美久久久| 午夜福利在线观看免费完整高清在 | 五月玫瑰六月丁香| 国产亚洲欧美98| 国产视频一区二区在线看| 哪里可以看免费的av片| 美女高潮的动态| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区三| 精品熟女少妇八av免费久了| 9191精品国产免费久久| 搡老岳熟女国产| av视频在线观看入口| 两个人看的免费小视频| 日韩精品青青久久久久久| 成人午夜高清在线视频| 中文字幕人妻熟人妻熟丝袜美 | 69人妻影院| 久久人妻av系列| 午夜福利在线在线| 天天躁日日操中文字幕| 精品国内亚洲2022精品成人| 午夜福利视频1000在线观看| 久久久久精品国产欧美久久久| 国产精品 国内视频| 美女 人体艺术 gogo| 国产伦在线观看视频一区| 国产欧美日韩一区二区三| 18禁黄网站禁片午夜丰满| 香蕉久久夜色| 欧美午夜高清在线| 最近在线观看免费完整版| 国产精品99久久久久久久久| 亚洲,欧美精品.| 三级毛片av免费| 在线十欧美十亚洲十日本专区| 久久婷婷人人爽人人干人人爱| 久久九九热精品免费| 色综合婷婷激情| www.色视频.com| 亚洲人成伊人成综合网2020| 成年版毛片免费区| 国产黄a三级三级三级人| 婷婷丁香在线五月| 国产亚洲精品久久久久久毛片| 久久精品人妻少妇| 乱人视频在线观看| 国产又黄又爽又无遮挡在线| 午夜福利在线在线| 国产成人影院久久av| 久久久精品大字幕| 亚洲精品日韩av片在线观看 | 亚洲av二区三区四区| 免费看a级黄色片| 桃红色精品国产亚洲av| 国产欧美日韩精品一区二区| avwww免费| 国产亚洲精品久久久com| 老司机深夜福利视频在线观看| 91麻豆精品激情在线观看国产| 女警被强在线播放| 国产精品一区二区三区四区免费观看 | 精品久久久久久久人妻蜜臀av| 少妇人妻一区二区三区视频| 午夜精品在线福利| 欧美日韩瑟瑟在线播放| 亚洲精品粉嫩美女一区| 久久久久性生活片| 成人特级黄色片久久久久久久| or卡值多少钱| 国产又黄又爽又无遮挡在线| 国产午夜精品久久久久久一区二区三区 | 中文资源天堂在线| 国产精品久久久人人做人人爽| 久久九九热精品免费| 免费人成在线观看视频色| 97人妻精品一区二区三区麻豆| 最近最新免费中文字幕在线| 亚洲在线观看片| 亚洲成人中文字幕在线播放| 国产亚洲精品久久久com| 国产精品亚洲美女久久久| 午夜免费观看网址| 久久久久性生活片| 国产亚洲精品久久久久久毛片| 亚洲熟妇熟女久久| 色综合欧美亚洲国产小说| 久久香蕉精品热| 欧美性感艳星| 真实男女啪啪啪动态图| 亚洲自拍偷在线| 嫁个100分男人电影在线观看| 日日夜夜操网爽| 美女大奶头视频| 精品一区二区三区视频在线 | 中文字幕人妻丝袜一区二区| 免费在线观看亚洲国产| 国产精品亚洲一级av第二区| 听说在线观看完整版免费高清| 欧美另类亚洲清纯唯美| 欧美+亚洲+日韩+国产| 最近在线观看免费完整版| 亚洲人成网站在线播放欧美日韩| 性色av乱码一区二区三区2| 日本三级黄在线观看| 久久久久亚洲av毛片大全| 欧美成人一区二区免费高清观看| 成人一区二区视频在线观看| 精品国产三级普通话版| 在线免费观看不下载黄p国产 | 深爱激情五月婷婷| 欧美高清成人免费视频www| 国产精品三级大全| 啦啦啦观看免费观看视频高清| eeuss影院久久| 国产视频内射| 亚洲专区中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 亚洲精品国产精品久久久不卡| 国产久久久一区二区三区| 亚洲自拍偷在线| 非洲黑人性xxxx精品又粗又长| 黄色日韩在线| 欧美激情在线99| 91九色精品人成在线观看| 亚洲av免费在线观看| 成人18禁在线播放| 五月伊人婷婷丁香| 久久久国产成人精品二区| 国产三级中文精品| 午夜精品久久久久久毛片777| 欧美日韩精品网址| a级毛片a级免费在线| 免费观看精品视频网站| 宅男免费午夜| 黄色丝袜av网址大全| 久久亚洲精品不卡| e午夜精品久久久久久久| 亚洲成人久久性| 免费观看的影片在线观看| 国产老妇女一区| 最近最新免费中文字幕在线| 一本综合久久免费| 国产三级黄色录像| 尤物成人国产欧美一区二区三区| 亚洲欧美日韩无卡精品| 免费在线观看亚洲国产| 精品电影一区二区在线| 麻豆成人午夜福利视频| 夜夜看夜夜爽夜夜摸| 免费看十八禁软件| www日本在线高清视频| 亚洲av一区综合| 成人永久免费在线观看视频| 一个人免费在线观看的高清视频| 精品一区二区三区av网在线观看| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 国产野战对白在线观看| 51国产日韩欧美| ponron亚洲| 国产成人aa在线观看| 免费在线观看影片大全网站| 一本久久中文字幕| 国产精品98久久久久久宅男小说| www.色视频.com| 又爽又黄无遮挡网站| 精品人妻1区二区| 亚洲人成电影免费在线| 人妻夜夜爽99麻豆av| 桃红色精品国产亚洲av| a级一级毛片免费在线观看| 久久中文看片网| 欧洲精品卡2卡3卡4卡5卡区| 日本免费一区二区三区高清不卡| 狠狠狠狠99中文字幕| 国产伦一二天堂av在线观看| 制服丝袜大香蕉在线| 亚洲最大成人中文| 免费人成视频x8x8入口观看| 老司机在亚洲福利影院| 岛国在线免费视频观看| 手机成人av网站| 久久6这里有精品| 老司机在亚洲福利影院| av视频在线观看入口| 成年版毛片免费区| 欧美xxxx黑人xx丫x性爽| 全区人妻精品视频| 午夜精品在线福利| 精品一区二区三区人妻视频| 别揉我奶头~嗯~啊~动态视频| 精品不卡国产一区二区三区| 亚洲精品久久国产高清桃花| 一区二区三区激情视频| 最近最新免费中文字幕在线| 国产精品一区二区三区四区免费观看 | 两性午夜刺激爽爽歪歪视频在线观看| 黄色视频,在线免费观看| 欧美高清成人免费视频www| 国内毛片毛片毛片毛片毛片| 香蕉久久夜色| 女人十人毛片免费观看3o分钟|