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

    半開式離心葉輪氣動(dòng)阻尼的影響規(guī)律研究*

    2021-03-24 03:39:08陳旭東楊樹華王曉放劉海濤
    風(fēng)機(jī)技術(shù) 2021年1期
    關(guān)鍵詞:阻尼比葉輪氣動(dòng)

    陳旭東 楊樹華 王曉放 劉海濤

    (1.大連理工大學(xué)能源與動(dòng)力學(xué)院;2.沈陽鼓風(fēng)機(jī)集團(tuán)股份有限公司)

    0 引言

    高壓比、大流量、高效率、寬運(yùn)行工況是離心壓縮機(jī)設(shè)計(jì)的發(fā)展趨勢[1]。復(fù)雜的運(yùn)行工況帶來的各種激振力可導(dǎo)致葉片動(dòng)應(yīng)力過大,從而引起葉片的高周疲勞破壞,疲勞破壞已是壓縮機(jī)主要的破壞形式之一[2]。準(zhǔn)確計(jì)算葉片的動(dòng)應(yīng)力是葉輪強(qiáng)度設(shè)計(jì)中的重要問題[3]。由于離心壓縮機(jī)葉輪結(jié)構(gòu)的振動(dòng)模態(tài)密集,葉輪轉(zhuǎn)速變化等因素影響,使得氣流激振力的頻率范圍寬泛,葉輪共振校核難以完全滿足要求。對于不滿足共振要求的振動(dòng)模態(tài),需要進(jìn)一步進(jìn)行詳細(xì)的葉輪動(dòng)應(yīng)力計(jì)算。葉輪結(jié)構(gòu)阻尼的計(jì)算是否準(zhǔn)確直接決定了葉輪動(dòng)應(yīng)力的計(jì)算精度。

    對于半開式離心壓縮機(jī),葉輪總阻尼包括材料阻尼和氣動(dòng)阻尼。國外學(xué)者在葉輪阻尼的實(shí)驗(yàn)測量方面做了相關(guān)研究,Kammerer[4]通過實(shí)驗(yàn)測試得到了某離心壓縮機(jī)葉輪的阻尼,以及材料阻尼和氣動(dòng)阻尼在總阻尼中所占的比例。Zemp[5]等利用類似的實(shí)驗(yàn)方法,研究了某離心壓縮機(jī)的氣動(dòng)阻尼以及材料阻尼。相關(guān)實(shí)驗(yàn)測試的主要方法是:首先在真空環(huán)境中測量葉輪的阻尼,此時(shí)得到的是材料阻尼。假設(shè)材料阻尼是常數(shù),再測量實(shí)際運(yùn)行工況下的阻尼,從中減去材料阻尼部分即可得到氣動(dòng)阻尼。這種方法得到的總阻尼雖然比較精確,但實(shí)驗(yàn)測量很難考慮葉輪的所有模態(tài)以及機(jī)組的各種工況,并且實(shí)驗(yàn)成本也比較高。此外,假設(shè)材料阻尼比在各種工況下是常數(shù)也并不合適,Rao[6]得到了材料阻尼比與結(jié)構(gòu)應(yīng)力幅值的非線性關(guān)系。因此,采用實(shí)驗(yàn)測量得到的材料阻尼和氣動(dòng)阻尼的比例可能與實(shí)際情況差別較大。

    利用數(shù)值計(jì)算方法準(zhǔn)確估計(jì)不同工況下葉輪的氣動(dòng)阻尼成為關(guān)鍵問題。Parthasarathy[7]利用數(shù)值計(jì)算,得到了某超音速風(fēng)機(jī)葉片的氣動(dòng)阻尼,并研究了葉間相角對氣動(dòng)阻尼的影響。池志強(qiáng)[8]通過數(shù)值方法分析了風(fēng)力機(jī)柔性葉片在多種隱態(tài)風(fēng)速下一階模態(tài)的氣動(dòng)阻尼。相關(guān)研究[9-10]主要針對軸流透平,對離心透平的研究較少。由于轉(zhuǎn)子葉盤的結(jié)構(gòu)特點(diǎn),軸流透平的葉片往往呈現(xiàn)行波型振動(dòng),特別對于軸流葉片一類高而短的結(jié)構(gòu),可假設(shè)所有葉片均固支于其根部振動(dòng),各葉片以同樣的頻率和振幅做簡諧振動(dòng),相鄰葉片的振動(dòng)只相差一個(gè)葉間相角,該方法被廣泛應(yīng)用于軸流葉片氣動(dòng)阻尼的計(jì)算中。但對于離心葉輪結(jié)構(gòu),因其葉片展向短而弦向長,葉-盤振動(dòng)必須整體分析,不同節(jié)徑的振動(dòng)不僅葉間相角不同,頻率亦發(fā)生變化,故離心葉輪氣動(dòng)阻尼的分析方法與軸流葉片有所不同[11]。近年來,工業(yè)離心壓縮機(jī)葉輪的疲勞破壞時(shí)有發(fā)生,葉輪結(jié)構(gòu)的疲勞強(qiáng)度校核是技術(shù)關(guān)鍵[12],對離心葉輪氣動(dòng)阻尼特性的深入研究有助于提高葉輪動(dòng)應(yīng)力計(jì)算的精度,進(jìn)而提升葉輪可靠性分析的技術(shù)水平。

    本文首先介紹了氣動(dòng)阻尼的數(shù)值計(jì)算方法。通過模態(tài)分析,得到葉輪的固有頻率、振型等數(shù)據(jù)。然后通過線性插值的方式將振型數(shù)據(jù)傳遞到流體仿真模型中,進(jìn)行氣動(dòng)阻尼計(jì)算及分析。在此基礎(chǔ)上,研究了葉片振幅、壓縮機(jī)進(jìn)口的流量、介質(zhì)溫度等對葉輪氣動(dòng)阻尼的影響。

    1 氣動(dòng)阻尼計(jì)算方法

    葉片表面的非定常氣動(dòng)載荷阻礙葉片振動(dòng)時(shí),氣動(dòng)載荷做負(fù)功,消耗結(jié)構(gòu)振動(dòng)能量,即氣體介質(zhì)具有正阻尼效應(yīng);相反,氣動(dòng)載荷對振動(dòng)中的葉片做正功時(shí),將使葉片振動(dòng)能量增加,引發(fā)葉片的顫振。由于小幅振動(dòng)引起的流場擾動(dòng)具有線性特點(diǎn),所以在葉片一個(gè)振動(dòng)周期內(nèi),葉片表面對流體所做的功(即阻尼耗功)為

    式中,t0是一個(gè)振動(dòng)周期的開始時(shí)刻;T是振動(dòng)周期;p是葉片表面壓力;v是振動(dòng)速度;n是葉片表面的單位法向量;A是葉片表面面積。葉片振動(dòng)的周期、頻率、振型等振動(dòng)數(shù)據(jù),由葉片的模態(tài)分析得到。根據(jù)等效黏性阻尼的定義[13],葉片氣動(dòng)阻尼比為

    式中,U0為葉片一個(gè)振動(dòng)周期內(nèi)的最大變形能,定義為

    式中,v是葉片體積。

    本文采用流固弱耦合方法來求解氣動(dòng)阻尼。流固弱耦合方法的主要思想是分別求解結(jié)構(gòu)方程及流體方程,利用線性插值方法將葉片的運(yùn)動(dòng)施加到流固耦合邊界上,實(shí)現(xiàn)葉片在流場中的振動(dòng)。由于Moffatt 和He[14]假設(shè)葉片的模態(tài)和固有頻率不受氣動(dòng)載荷的影響,因此在結(jié)構(gòu)動(dòng)力學(xué)分析中只需考慮無阻尼自由振動(dòng)。主要過程分為以下3步:

    1)計(jì)算葉片的固有模態(tài),得到振型及頻率;

    2)利用線性插值方法,設(shè)置最大振幅,將葉片按某階固有振型振動(dòng)的運(yùn)動(dòng)規(guī)律施加到耦合邊界上,對流場進(jìn)行擾動(dòng);

    3)計(jì)算流場,得到非定常氣動(dòng)力,并結(jié)合已知的葉片振動(dòng)位移,得到一個(gè)葉片振動(dòng)周期內(nèi)的非定常氣動(dòng)功。

    氣動(dòng)阻尼計(jì)算利用的是傅里葉變換法(Fourier Transformation Method)。傅里葉變換法的主要思想是,在不同時(shí)刻的情況下,相鄰節(jié)距邊界是具有周期性的。這種方法不需要儲(chǔ)存整個(gè)周期所有節(jié)距邊界上的信號(hào),只需儲(chǔ)存傅里葉系數(shù)Am來得到任意時(shí)刻的解[15]。時(shí)間傅里葉級數(shù)展開如下

    這種算法需要利用至少兩個(gè)葉片通道,如圖1 所示。計(jì)算中,采集的是相鄰兩個(gè)葉片通道交界面上(邊界2)的傅里葉系數(shù)。交界面離周期邊界遠(yuǎn),信號(hào)質(zhì)量高,計(jì)算比較準(zhǔn)確。

    圖1 葉片通道簡化模型Fig.1 The simplified model for passages

    圖2 相鄰振動(dòng)葉片間的相位差Fig.2 The phase difference between adjacent vibrating blades

    對于典型的葉輪機(jī)械氣體動(dòng)力學(xué)問題(即葉片顫振和強(qiáng)迫響應(yīng)),葉片-輪盤結(jié)構(gòu)通常為周向勻速旋轉(zhuǎn)的循環(huán)對稱模式。對于這類葉片振動(dòng)分析,相鄰振動(dòng)葉片間的相位差如圖2所示,它可以表示為

    2 葉輪模態(tài)分析

    本文研究對象為一離心壓縮機(jī)半開式徑向葉輪,其出口直徑1 179mm,出口寬度90.39mm,葉片數(shù)為19。由于葉輪是循環(huán)對稱結(jié)構(gòu),可以采用葉輪扇區(qū)模型進(jìn)行有限元分析。圖3是葉輪有限元分析模型,葉輪葉片部分采用結(jié)構(gòu)化網(wǎng)格,輪盤部分采用非結(jié)構(gòu)化網(wǎng)格。考慮葉輪離心力的影響,采用有限元分析軟件Ansys進(jìn)行模態(tài)分析。葉輪一階模態(tài)頻率為384.6Hz,一階模態(tài)振型如圖4所示。

    圖3 葉輪有限元模型Fig.3 FE model for impeller

    圖4 葉輪一階模態(tài)振型Fig.4 Deformation for the 1st impeller mode

    3 半開式離心葉輪氣動(dòng)阻尼分析

    3.1 氣動(dòng)阻尼計(jì)算參數(shù)和設(shè)置

    葉輪流體域和固體域交界面網(wǎng)格節(jié)點(diǎn)不一一對應(yīng),需要通過線性插值的方式,將葉片振型的節(jié)點(diǎn)位移數(shù)據(jù)傳遞到交界面上的流體域網(wǎng)格節(jié)點(diǎn)。額定工況下,葉輪的進(jìn)口總壓為99 000Pa,進(jìn)口總溫為298K,質(zhì)量流量為54.88kg/s,轉(zhuǎn)速5 556r/min。Ansys CFX中葉片通道的網(wǎng)格劃分如圖5所示,單個(gè)通道網(wǎng)格單元總數(shù)約為14萬。計(jì)算工質(zhì)為理想空氣。由于離心葉輪中存在較強(qiáng)的逆壓梯度,其可能導(dǎo)致流動(dòng)分離。而SST模型對于強(qiáng)逆壓梯度的邊界層流動(dòng)和流動(dòng)分離預(yù)測有很好的表現(xiàn),故湍流模型采用SST模型[15]。CFD計(jì)算模型如圖6所示。

    圖5 葉片流道網(wǎng)格Fig.5 Mesh of passage for CFD calculation

    圖6 CFD計(jì)算模型Fig.6 CFD calculation model

    瞬態(tài)CFD計(jì)算以定常計(jì)算結(jié)果為初值,計(jì)算中,一個(gè)振動(dòng)周期設(shè)置70 個(gè)時(shí)間步。額定工況下,葉輪一階模態(tài)的氣動(dòng)阻尼計(jì)算結(jié)果如圖7 所示,由式(1)計(jì)算得到氣動(dòng)阻尼在一個(gè)振動(dòng)周期內(nèi)的耗能為0.000 284J。根據(jù)式(3)可得葉片的應(yīng)變能為U0=0.018 2J,再由式(2)得氣動(dòng)阻尼比為0.001 24。

    圖7 氣動(dòng)阻尼耗功收斂曲線Fig.7 Convergence for aerodynamic damping work

    為了盡量消除網(wǎng)格數(shù)量引起的離散誤差,本文選擇4 種網(wǎng)格數(shù)進(jìn)行網(wǎng)格無關(guān)性分析,結(jié)果如表1 所示。由表1可以發(fā)現(xiàn),網(wǎng)格數(shù)14萬時(shí)計(jì)算結(jié)果基本滿足計(jì)算精度要求。因此,綜合考慮計(jì)算時(shí)間和計(jì)算資源,本文決定采用14萬網(wǎng)格數(shù)。

    表1 網(wǎng)格無關(guān)性分析Tab.1 Grid independence analysis

    3.2 葉片振幅對氣動(dòng)阻尼的影響

    在葉輪阻尼及動(dòng)應(yīng)力未知的情況下,葉片的振幅是無法準(zhǔn)確估計(jì)的。這里研究葉片振幅大小對氣動(dòng)阻尼的影響,考察了額定工況半開式離心葉輪前三階模態(tài)的氣動(dòng)阻尼比,計(jì)算結(jié)果如圖8 所示。對于同一模態(tài),在不同振幅下,葉片的氣動(dòng)阻尼比基本保持不變,相關(guān)文獻(xiàn)也得到了類似結(jié)果[16-18]。Bidkar[18]引進(jìn)了一個(gè)無量綱量KC:

    其中,A是振幅;c是結(jié)構(gòu)特征長度。

    當(dāng)KC 值遠(yuǎn)小于1 時(shí),氣動(dòng)阻尼比幾乎與振幅無關(guān)。對于離心壓縮機(jī)葉輪而言,振幅A 的范圍在10-5~10-3m,特征長度c的范圍在10-1~1m,KC值遠(yuǎn)小于1。因此,在同一工況某模態(tài)下葉片氣動(dòng)阻尼比可近似為常數(shù)。

    圖8 前三階模態(tài)不同振幅下的葉片氣動(dòng)阻尼比Fig.8 Aerodynamic damping ratio for different modes

    3.3 進(jìn)口氣體流量及溫度對氣動(dòng)阻尼的影響

    在機(jī)組實(shí)際運(yùn)行時(shí),根據(jù)生產(chǎn)工藝需求,需要改變進(jìn)口氣體流量。為了研究流量對氣動(dòng)阻尼的影響,本文計(jì)算了不同流量下葉片一階模態(tài)的氣動(dòng)阻尼比,結(jié)果如圖9所示。

    圖9 不同流量下的葉片氣動(dòng)阻尼比Fig.9 Aerodynamic damping ratio for different flow rate

    可以發(fā)現(xiàn),同一模態(tài)同一工況下,不同葉片振幅下的葉片氣動(dòng)阻尼比基本不變。表2是不同流量下,取了各振幅的氣動(dòng)阻尼比平均值。可以發(fā)現(xiàn),隨著進(jìn)口流量的增加,葉片氣動(dòng)阻尼比顯著增大。

    表2 不同流量下的平均氣動(dòng)阻尼比Tab.2 Average aerodynamic damping ratio for different flow rate

    從壓縮機(jī)運(yùn)行特性方面分析流量對氣動(dòng)阻尼的影響規(guī)律。離心壓縮機(jī)軸功率隨進(jìn)口氣體流量的增大而增大,功率可表示為

    其中,T0為作用在葉片上的扭矩;ω0為角頻率。轉(zhuǎn)速不變的情況下,葉片上的扭矩隨功率的增大而增大。而作用在葉片上的扭矩由葉片表面的壓力差決定。因此,氣體流量的增大會(huì)使得葉片表面壓差(壓力載荷)增大。

    如圖11 和圖12 所示是不同流量下,0.5 和0.9 倍葉高處(見圖10)的壓力面與吸力面氣動(dòng)壓力差沿流線的分布(葉片一個(gè)振動(dòng)周期內(nèi)的平均值)。可以發(fā)現(xiàn),流量增大時(shí),葉片表面壓差增大。因此,由式(1)可知,氣動(dòng)阻尼耗功增大,即氣動(dòng)阻尼比隨流量的增大而增大。對于本文研究的半開式離心葉輪,進(jìn)口質(zhì)量流量由85%額定流量變化到104%額定流量,氣動(dòng)阻尼比由0.000 50變?yōu)?.001 61,增大了2.22倍。

    圖10 葉輪葉片子午流道Fig.10 Meridional view of blade

    圖11 0.5倍葉高處葉片表面的壓差Fig.11 Pressure difference of blade surface at 0.5 blade height

    圖12 0.9倍葉高處葉片表面的壓差Fig.12 Pressure difference of blade surface at 0.9 blade height

    最后,考察氣體進(jìn)口溫度對氣動(dòng)阻尼的影響。由于壓縮機(jī)進(jìn)口氣體參數(shù)為當(dāng)?shù)卮髿鈪?shù),進(jìn)口氣體溫度隨環(huán)境的變化而變化,不是恒定值。因此,本文研究了兩種工況下的氣動(dòng)阻尼結(jié)果。這兩種工況包括夏季工況,進(jìn)口溫度為25℃;冬季工況,進(jìn)口溫度為-10℃。由于進(jìn)口壓力的變化相對較小,本文忽略壓力對計(jì)算結(jié)果的影響,即兩種工況下的進(jìn)口總壓取值都為99 000Pa。

    壓縮機(jī)運(yùn)行工況的相似條件為:

    其中,Qa,Qa' 為兩種工況下的進(jìn)口氣體體積流量;ml為葉輪結(jié)構(gòu)幾何比例;R' ,R 為氣體常數(shù);Ta' ,Ta為進(jìn)口氣體溫度。對于同一臺(tái)機(jī)組,ml=1;對于同一種氣體,R' =R。因此,冬季工況下,進(jìn)口體積流量為

    圖13 是95%~104%額定流量下的氣動(dòng)阻尼計(jì)算結(jié)果。冬季工況與夏季工況計(jì)算結(jié)果類似,氣動(dòng)阻尼比隨著進(jìn)口流量的增大而增大。流量由95%額定流量增大到104%額定流量,氣動(dòng)阻尼比增大了86.1%。

    表3 不同進(jìn)口氣體溫度下葉片的氣動(dòng)阻尼比Tab.3 The aerodynamic damping ratio for different inlet air temperature

    4 結(jié)論

    圖13 冬季工況下不同流量的氣動(dòng)阻尼比Fig.13 Aerodynamic damping ratio for different flow rate in winter

    對于本文研究的半開式離心葉輪,冬季工況下,進(jìn)口流量為額定流量的90%時(shí),數(shù)值計(jì)算不能收斂。這是由于在冬季工況下,進(jìn)口氣體溫度較低,密度較大,質(zhì)量流量一定的情況下,體積流量減小。機(jī)組在冬季變工況運(yùn)行,如果進(jìn)口體積流量進(jìn)一步減小,機(jī)組很可能進(jìn)入非穩(wěn)定流動(dòng)狀態(tài),甚至發(fā)生喘振。此時(shí),數(shù)值計(jì)算得到的流場也是非穩(wěn)定的,無法得到一個(gè)恒定的氣動(dòng)阻尼值。

    最后,由表3可知,額定流量下,對于葉輪的一階模態(tài),夏季工況(進(jìn)口氣體溫度為25℃)和冬季工況(進(jìn)口氣體溫度為-10℃)葉片的氣動(dòng)阻尼比分別為0.001 24和0.003 03,冬季工況葉片的氣動(dòng)阻尼比夏季工況增大了1.46倍。

    本文研究了某半開式離心葉輪的氣動(dòng)阻尼特性,并且對氣動(dòng)阻尼的影響因素進(jìn)行了理論分析。流固弱耦合計(jì)算可以使固體域與流體域分別建模,通過線性插值完成位移載荷的傳遞,從而使計(jì)算大大簡化,實(shí)現(xiàn)任意固有頻率及不同流場狀態(tài)下的氣動(dòng)阻尼的計(jì)算。

    分析發(fā)現(xiàn)葉輪在同一工況同一模態(tài)下的氣動(dòng)阻尼比與葉片振幅基本無關(guān)。數(shù)值計(jì)算結(jié)果表明,進(jìn)口氣體溫度一定時(shí),葉輪氣動(dòng)阻尼比隨流量的增大而增大;而質(zhì)量流量一定時(shí),葉輪氣動(dòng)阻尼比隨進(jìn)口氣體溫度的減小而增大。因此,在設(shè)計(jì)中需要考慮變工況對氣動(dòng)阻尼的影響。

    半開式葉輪的總阻尼包括氣動(dòng)阻尼和材料阻尼兩部分,本文只對氣動(dòng)阻尼特性進(jìn)行了初步研究,需要進(jìn)一步綜合考慮氣動(dòng)阻尼和材料阻尼的共同作用,才能獲得準(zhǔn)確的動(dòng)應(yīng)力計(jì)算結(jié)果。

    猜你喜歡
    阻尼比葉輪氣動(dòng)
    中寰氣動(dòng)執(zhí)行機(jī)構(gòu)
    1.4317 QT2鋼在高能泵葉輪上的應(yīng)用
    基于NACA0030的波紋狀翼型氣動(dòng)特性探索
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    黏滯阻尼器在時(shí)程分析下的附加有效阻尼比研究
    基于反饋線性化的RLV氣動(dòng)控制一體化設(shè)計(jì)
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    應(yīng)用石膏型快速精密鑄造技術(shù)制造葉輪
    結(jié)構(gòu)構(gòu)件阻尼比對大跨度懸索橋地震響應(yīng)的影響
    離心泵葉輪切割方法
    日本午夜av视频| 男人狂女人下面高潮的视频| 成年人午夜在线观看视频| 国产美女午夜福利| 国产无遮挡羞羞视频在线观看| 欧美3d第一页| 国产精品久久久久久久久免| 国产欧美日韩一区二区三区在线 | 日本猛色少妇xxxxx猛交久久| 亚洲综合色惰| 色吧在线观看| 欧美精品国产亚洲| a级毛片免费高清观看在线播放| 人人妻人人添人人爽欧美一区卜| 在线观看免费高清a一片| 一级毛片我不卡| 99久久精品国产国产毛片| 丰满少妇做爰视频| 欧美日韩av久久| 婷婷色av中文字幕| 老熟女久久久| 久久国产精品男人的天堂亚洲 | 我的女老师完整版在线观看| 最黄视频免费看| 亚洲精品,欧美精品| 岛国毛片在线播放| 久久毛片免费看一区二区三区| 久热久热在线精品观看| 老司机影院成人| 美女xxoo啪啪120秒动态图| √禁漫天堂资源中文www| 91aial.com中文字幕在线观看| 国产一区二区三区综合在线观看 | 日韩一本色道免费dvd| 成人亚洲欧美一区二区av| 一级,二级,三级黄色视频| 久久热精品热| 国产熟女午夜一区二区三区 | 国产高清三级在线| 国产精品一二三区在线看| 在线观看www视频免费| 热re99久久精品国产66热6| 中文在线观看免费www的网站| av专区在线播放| 亚洲国产精品999| 免费大片黄手机在线观看| 欧美三级亚洲精品| 美女大奶头黄色视频| 婷婷色麻豆天堂久久| 只有这里有精品99| 亚洲av男天堂| 亚洲精品国产色婷婷电影| 久久97久久精品| 亚洲av免费高清在线观看| 99热6这里只有精品| 内地一区二区视频在线| 夜夜爽夜夜爽视频| 日日摸夜夜添夜夜添av毛片| 免费看不卡的av| 国产色婷婷99| 亚洲久久久国产精品| 全区人妻精品视频| av播播在线观看一区| 丰满人妻一区二区三区视频av| 国产亚洲最大av| 91在线精品国自产拍蜜月| 美女cb高潮喷水在线观看| a级毛色黄片| 十八禁高潮呻吟视频 | 一级a做视频免费观看| 亚洲,一卡二卡三卡| 国产片特级美女逼逼视频| 欧美成人精品欧美一级黄| 高清在线视频一区二区三区| 午夜视频国产福利| h日本视频在线播放| 亚洲欧美精品自产自拍| 三级国产精品欧美在线观看| 婷婷色麻豆天堂久久| a级毛片在线看网站| 大话2 男鬼变身卡| 国产男女超爽视频在线观看| 国产黄色免费在线视频| 色婷婷av一区二区三区视频| a级毛片免费高清观看在线播放| 日日爽夜夜爽网站| 国产精品99久久99久久久不卡 | 99九九线精品视频在线观看视频| 最近中文字幕高清免费大全6| 亚洲欧美中文字幕日韩二区| 亚洲人成网站在线播| 青春草视频在线免费观看| 日韩,欧美,国产一区二区三区| 亚洲电影在线观看av| 国产黄频视频在线观看| 亚洲av男天堂| 国产高清有码在线观看视频| 五月天丁香电影| 精品99又大又爽又粗少妇毛片| 精品一品国产午夜福利视频| 国产精品.久久久| 亚洲av中文av极速乱| 国产伦精品一区二区三区视频9| √禁漫天堂资源中文www| 99九九在线精品视频 | 寂寞人妻少妇视频99o| 午夜av观看不卡| 有码 亚洲区| 春色校园在线视频观看| 国产一级毛片在线| h视频一区二区三区| 大香蕉久久网| 日日啪夜夜撸| 国产精品人妻久久久影院| 国产高清有码在线观看视频| 一级黄片播放器| 国产乱人偷精品视频| 亚洲精华国产精华液的使用体验| 你懂的网址亚洲精品在线观看| 久久99精品国语久久久| 一本—道久久a久久精品蜜桃钙片| 黄色怎么调成土黄色| 少妇 在线观看| 99久久人妻综合| 黄色一级大片看看| 男人舔奶头视频| 99久久中文字幕三级久久日本| 青青草视频在线视频观看| 日韩 亚洲 欧美在线| 国产高清三级在线| 国产精品一区二区三区四区免费观看| 三级国产精品欧美在线观看| 一区二区三区精品91| 99久久精品热视频| 亚洲av男天堂| 亚洲内射少妇av| 99九九线精品视频在线观看视频| 久久青草综合色| 日韩免费高清中文字幕av| 久久精品夜色国产| 人妻系列 视频| 久久精品国产a三级三级三级| 日日啪夜夜爽| 亚洲人成网站在线观看播放| 久久这里有精品视频免费| 午夜视频国产福利| 人人妻人人澡人人爽人人夜夜| 欧美日韩精品成人综合77777| 狂野欧美激情性xxxx在线观看| 男人和女人高潮做爰伦理| 最近的中文字幕免费完整| 久热久热在线精品观看| 中文字幕人妻熟人妻熟丝袜美| 免费看av在线观看网站| av线在线观看网站| 亚洲精华国产精华液的使用体验| 高清午夜精品一区二区三区| 精华霜和精华液先用哪个| 亚洲欧洲精品一区二区精品久久久 | 少妇裸体淫交视频免费看高清| 狂野欧美白嫩少妇大欣赏| 一个人免费看片子| 欧美区成人在线视频| 成人国产麻豆网| 免费观看无遮挡的男女| 激情五月婷婷亚洲| 精品一区在线观看国产| 高清不卡的av网站| 另类精品久久| 久久久精品免费免费高清| 寂寞人妻少妇视频99o| 美女国产视频在线观看| 国产免费视频播放在线视频| 欧美日韩亚洲高清精品| 日日摸夜夜添夜夜添av毛片| 一级,二级,三级黄色视频| 欧美国产精品一级二级三级 | av.在线天堂| 九九久久精品国产亚洲av麻豆| 中文字幕av电影在线播放| 香蕉精品网在线| 两个人免费观看高清视频 | 国产片特级美女逼逼视频| 日本欧美视频一区| 日本免费在线观看一区| 一级毛片 在线播放| 国产成人a∨麻豆精品| 91精品国产九色| 日本黄色日本黄色录像| 日本午夜av视频| 色视频在线一区二区三区| 亚洲成人一二三区av| 欧美高清成人免费视频www| 汤姆久久久久久久影院中文字幕| 亚洲精品国产色婷婷电影| 久久精品夜色国产| 久久精品熟女亚洲av麻豆精品| av有码第一页| 精品久久久久久久久亚洲| 精品久久国产蜜桃| 一本久久精品| 特大巨黑吊av在线直播| 人妻 亚洲 视频| 免费av不卡在线播放| 国产一区二区三区av在线| 美女xxoo啪啪120秒动态图| 免费黄色在线免费观看| 欧美97在线视频| 我的老师免费观看完整版| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美激情极品国产一区二区三区 | 国产免费又黄又爽又色| 日本免费在线观看一区| www.色视频.com| 夜夜骑夜夜射夜夜干| 涩涩av久久男人的天堂| 国产精品一区二区在线观看99| 日韩制服骚丝袜av| 美女内射精品一级片tv| 十八禁高潮呻吟视频 | 一区二区三区乱码不卡18| 极品少妇高潮喷水抽搐| 插逼视频在线观看| 亚洲av男天堂| 亚洲精品国产av成人精品| 一级爰片在线观看| 人妻夜夜爽99麻豆av| 国产成人精品婷婷| 日韩亚洲欧美综合| 久久婷婷青草| 日韩欧美 国产精品| 精品久久久久久久久av| 狂野欧美激情性bbbbbb| 亚洲欧美成人精品一区二区| 这个男人来自地球电影免费观看 | 亚洲精品自拍成人| 99九九线精品视频在线观看视频| 国产69精品久久久久777片| 97超碰精品成人国产| 日产精品乱码卡一卡2卡三| 久久久久人妻精品一区果冻| 久久精品久久久久久久性| 免费大片黄手机在线观看| 一区二区av电影网| 国产在视频线精品| 国产色婷婷99| 一级,二级,三级黄色视频| 成人影院久久| 欧美丝袜亚洲另类| 精品人妻熟女毛片av久久网站| 欧美日韩精品成人综合77777| 久久韩国三级中文字幕| 中文在线观看免费www的网站| 亚洲国产精品一区二区三区在线| √禁漫天堂资源中文www| 一级毛片我不卡| 亚洲精品成人av观看孕妇| 两个人免费观看高清视频 | 六月丁香七月| 成人二区视频| 91精品一卡2卡3卡4卡| 永久免费av网站大全| 久久久国产一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲国产欧美日韩在线播放 | 日日爽夜夜爽网站| videos熟女内射| 亚洲久久久国产精品| 深夜a级毛片| 久久影院123| 麻豆成人午夜福利视频| av专区在线播放| 少妇人妻一区二区三区视频| 日韩成人av中文字幕在线观看| 丝瓜视频免费看黄片| 午夜福利在线观看免费完整高清在| 少妇人妻久久综合中文| 18禁在线无遮挡免费观看视频| 欧美一级a爱片免费观看看| 男女免费视频国产| 少妇的逼好多水| 中文字幕精品免费在线观看视频 | 一区二区av电影网| 亚洲四区av| 久久久精品免费免费高清| 久久99精品国语久久久| 极品教师在线视频| 国产一区有黄有色的免费视频| 久久99一区二区三区| 激情五月婷婷亚洲| 久热久热在线精品观看| av福利片在线观看| 2018国产大陆天天弄谢| 丝袜脚勾引网站| 人妻一区二区av| 搡女人真爽免费视频火全软件| 一本一本综合久久| 91精品一卡2卡3卡4卡| 国产极品粉嫩免费观看在线 | 一级a做视频免费观看| 建设人人有责人人尽责人人享有的| 性高湖久久久久久久久免费观看| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 精品人妻熟女毛片av久久网站| 亚洲综合色惰| 国产精品一区www在线观看| 国产精品99久久久久久久久| 日本vs欧美在线观看视频 | 久久人妻熟女aⅴ| 91久久精品国产一区二区三区| 中文字幕制服av| 美女脱内裤让男人舔精品视频| 国产亚洲欧美精品永久| 韩国av在线不卡| 亚洲成色77777| 国产一区亚洲一区在线观看| 国产中年淑女户外野战色| 不卡视频在线观看欧美| 一区二区av电影网| 丝袜脚勾引网站| 亚洲欧美日韩东京热| 97超视频在线观看视频| 国产视频内射| 欧美日本中文国产一区发布| 婷婷色综合大香蕉| videossex国产| 亚洲人成网站在线观看播放| 免费观看a级毛片全部| 91在线精品国自产拍蜜月| 亚洲国产成人一精品久久久| 波野结衣二区三区在线| 人人妻人人澡人人看| 久久国内精品自在自线图片| 久久av网站| 日本vs欧美在线观看视频 | 男女免费视频国产| 青春草亚洲视频在线观看| 中文字幕制服av| 精品人妻偷拍中文字幕| 黄色怎么调成土黄色| 午夜影院在线不卡| 免费人妻精品一区二区三区视频| 91在线精品国自产拍蜜月| av在线老鸭窝| 黄片无遮挡物在线观看| 午夜91福利影院| 午夜免费观看性视频| 国产伦精品一区二区三区视频9| 国产精品一区www在线观看| 大香蕉久久网| √禁漫天堂资源中文www| 高清毛片免费看| 2022亚洲国产成人精品| 欧美日韩国产mv在线观看视频| 久久这里有精品视频免费| 一区二区av电影网| 简卡轻食公司| 伦理电影大哥的女人| 亚洲内射少妇av| 自线自在国产av| 亚洲av中文av极速乱| 精品少妇内射三级| 精品一区二区三卡| 黄片无遮挡物在线观看| 日本欧美国产在线视频| 亚洲性久久影院| 五月玫瑰六月丁香| 亚洲国产欧美日韩在线播放 | 女性被躁到高潮视频| 精品人妻熟女毛片av久久网站| 啦啦啦在线观看免费高清www| 晚上一个人看的免费电影| 午夜福利在线观看免费完整高清在| 日韩成人av中文字幕在线观看| 男人添女人高潮全过程视频| 美女主播在线视频| 老司机亚洲免费影院| 内射极品少妇av片p| 大片电影免费在线观看免费| 国产精品伦人一区二区| 大话2 男鬼变身卡| 久久6这里有精品| 赤兔流量卡办理| 我的女老师完整版在线观看| 超碰97精品在线观看| 亚洲成人av在线免费| 国产精品女同一区二区软件| 涩涩av久久男人的天堂| 午夜av观看不卡| 国产黄色视频一区二区在线观看| 日韩熟女老妇一区二区性免费视频| 97在线视频观看| 老女人水多毛片| 久久国产乱子免费精品| 亚洲国产色片| 成人亚洲欧美一区二区av| 色视频在线一区二区三区| h视频一区二区三区| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| av.在线天堂| 美女福利国产在线| 亚洲无线观看免费| 国产精品人妻久久久影院| 成年av动漫网址| 国产成人精品无人区| 国产精品女同一区二区软件| 成年av动漫网址| 自拍偷自拍亚洲精品老妇| 国产淫片久久久久久久久| 99re6热这里在线精品视频| 亚洲一级一片aⅴ在线观看| 国产 精品1| 成人亚洲精品一区在线观看| 极品人妻少妇av视频| 国产亚洲精品久久久com| 亚洲第一av免费看| 日本午夜av视频| 女性被躁到高潮视频| 欧美人与善性xxx| 啦啦啦在线观看免费高清www| 国产日韩欧美视频二区| 免费av不卡在线播放| 久久6这里有精品| 伦精品一区二区三区| 亚洲内射少妇av| 亚洲精品一区蜜桃| 91午夜精品亚洲一区二区三区| 天堂8中文在线网| 99久久人妻综合| 91午夜精品亚洲一区二区三区| 大香蕉97超碰在线| 伦理电影大哥的女人| 亚洲av福利一区| 中文字幕亚洲精品专区| 这个男人来自地球电影免费观看 | 成人亚洲精品一区在线观看| 亚洲精品国产成人久久av| 亚洲国产毛片av蜜桃av| 亚洲av成人精品一二三区| 欧美3d第一页| 成人毛片a级毛片在线播放| 亚洲高清免费不卡视频| 欧美+日韩+精品| 久久午夜综合久久蜜桃| 国产伦精品一区二区三区四那| 国产一区二区在线观看av| 中文字幕亚洲精品专区| 亚洲精品,欧美精品| 久久久久国产精品人妻一区二区| 熟女电影av网| 午夜91福利影院| 三级经典国产精品| 国产日韩一区二区三区精品不卡 | 妹子高潮喷水视频| 91久久精品电影网| av国产久精品久网站免费入址| 一个人看视频在线观看www免费| 噜噜噜噜噜久久久久久91| 秋霞伦理黄片| 777米奇影视久久| 极品少妇高潮喷水抽搐| 精品人妻熟女av久视频| 欧美激情极品国产一区二区三区 | 天堂中文最新版在线下载| 91精品伊人久久大香线蕉| 亚洲,欧美,日韩| 成人漫画全彩无遮挡| 免费少妇av软件| 日韩三级伦理在线观看| av在线老鸭窝| 十分钟在线观看高清视频www | 国产熟女欧美一区二区| 亚洲欧美日韩卡通动漫| 免费看光身美女| 观看美女的网站| 中文资源天堂在线| av天堂久久9| 国产成人午夜福利电影在线观看| 精品久久久久久久久av| 亚洲三级黄色毛片| 各种免费的搞黄视频| 在线亚洲精品国产二区图片欧美 | 欧美日韩在线观看h| 久久精品国产亚洲网站| 91久久精品电影网| 成人无遮挡网站| 亚洲欧美日韩东京热| 丰满乱子伦码专区| 观看免费一级毛片| 3wmmmm亚洲av在线观看| 狂野欧美白嫩少妇大欣赏| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 日韩一区二区视频免费看| 色婷婷av一区二区三区视频| 欧美丝袜亚洲另类| 天天操日日干夜夜撸| 我要看黄色一级片免费的| 亚洲av二区三区四区| 久久精品久久精品一区二区三区| 亚洲欧美成人精品一区二区| 97在线视频观看| 欧美高清成人免费视频www| 中国国产av一级| 人人妻人人爽人人添夜夜欢视频 | 欧美3d第一页| 久久久午夜欧美精品| 亚洲精品第二区| 婷婷色麻豆天堂久久| 日韩一区二区三区影片| 日产精品乱码卡一卡2卡三| 免费在线观看成人毛片| 亚洲av中文av极速乱| 亚洲欧美日韩另类电影网站| 蜜臀久久99精品久久宅男| 久久99热6这里只有精品| 久久99蜜桃精品久久| 亚洲成人一二三区av| 五月伊人婷婷丁香| 午夜精品国产一区二区电影| 又大又黄又爽视频免费| 亚洲精品国产av蜜桃| 亚洲天堂av无毛| 99热这里只有精品一区| 一区二区av电影网| 在线观看人妻少妇| 亚洲av福利一区| av专区在线播放| 在线观看国产h片| 国产一区二区三区综合在线观看 | 在线看a的网站| 国产精品嫩草影院av在线观看| 欧美精品一区二区大全| 狂野欧美激情性bbbbbb| 亚洲真实伦在线观看| 妹子高潮喷水视频| av免费观看日本| 九九爱精品视频在线观看| 亚洲综合色惰| 成人美女网站在线观看视频| 有码 亚洲区| 91精品一卡2卡3卡4卡| 视频中文字幕在线观看| www.av在线官网国产| 18禁在线无遮挡免费观看视频| 免费观看av网站的网址| 国产淫片久久久久久久久| 天天躁夜夜躁狠狠久久av| 久久精品国产亚洲网站| 日本av免费视频播放| 亚洲人成网站在线播| 一本—道久久a久久精品蜜桃钙片| 大又大粗又爽又黄少妇毛片口| 久久久精品94久久精品| 尾随美女入室| 高清不卡的av网站| 免费av不卡在线播放| 99久久精品国产国产毛片| 永久免费av网站大全| 日韩一区二区三区影片| 69精品国产乱码久久久| av线在线观看网站| 色婷婷av一区二区三区视频| 久久国产精品男人的天堂亚洲 | 精品酒店卫生间| 六月丁香七月| 日本黄大片高清| 久久久久国产精品人妻一区二区| 中文欧美无线码| 久久久久久人妻| 一本色道久久久久久精品综合| h视频一区二区三区| 中国国产av一级| 欧美日本中文国产一区发布| 国模一区二区三区四区视频| 国产有黄有色有爽视频| 国内精品宾馆在线| av网站免费在线观看视频| 免费黄色在线免费观看| 午夜av观看不卡| 国产深夜福利视频在线观看| av免费在线看不卡| 日韩一区二区视频免费看| 偷拍熟女少妇极品色| 久久久久久久久久久免费av| 大香蕉97超碰在线| 国产黄色免费在线视频| 亚洲国产成人一精品久久久| 一级毛片 在线播放| a级一级毛片免费在线观看| 久久久久人妻精品一区果冻| 少妇高潮的动态图| 春色校园在线视频观看| 国语对白做爰xxxⅹ性视频网站| 国产男女超爽视频在线观看| 国内精品宾馆在线| 久久久久久久久久久久大奶| 国产精品国产三级国产专区5o| 国产极品天堂在线| 在线亚洲精品国产二区图片欧美 | 我的女老师完整版在线观看| 大又大粗又爽又黄少妇毛片口| 一边亲一边摸免费视频| 久久婷婷青草| 美女xxoo啪啪120秒动态图| 色94色欧美一区二区| 亚洲精品日本国产第一区| 日本与韩国留学比较| 亚洲精品乱码久久久v下载方式| 在线观看三级黄色| 欧美老熟妇乱子伦牲交| 热re99久久国产66热| 免费观看在线日韩|