• <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片中文字幕在线观看| 久久人人精品亚洲av| 校园春色视频在线观看| 国产精品久久久久久人妻精品电影| 18禁国产床啪视频网站| 国产精品久久久av美女十八| 午夜精品国产一区二区电影| 亚洲专区中文字幕在线| 身体一侧抽搐| 麻豆成人av在线观看| 一边摸一边抽搐一进一出视频| 一本综合久久免费| 咕卡用的链子| xxx96com| 久久青草综合色| 香蕉久久夜色| 精品久久久久久,| 国产成人精品在线电影| 无人区码免费观看不卡| 午夜日韩欧美国产| 国产野战对白在线观看| 性色av乱码一区二区三区2| x7x7x7水蜜桃| 美国免费a级毛片| 自线自在国产av| 国产成年人精品一区二区| 久久亚洲真实| 国产精品一区二区在线不卡| 在线观看www视频免费| 国产欧美日韩一区二区三| 免费搜索国产男女视频| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 日韩欧美三级三区| 嫁个100分男人电影在线观看| 日本三级黄在线观看| 757午夜福利合集在线观看| 久久国产乱子伦精品免费另类| 变态另类丝袜制服| 俄罗斯特黄特色一大片| 欧美成人性av电影在线观看| 国产一区二区三区在线臀色熟女| 成在线人永久免费视频| 69av精品久久久久久| 两个人视频免费观看高清| 国产成人精品久久二区二区91| 国内精品久久久久久久电影| 久久精品国产综合久久久| 好男人电影高清在线观看| 欧美日本中文国产一区发布| 国产成人一区二区三区免费视频网站| 欧美乱码精品一区二区三区| 国产在线观看jvid| 亚洲黑人精品在线| 美女扒开内裤让男人捅视频| 一边摸一边做爽爽视频免费| 久久香蕉国产精品| 在线观看免费午夜福利视频| 国产亚洲精品av在线| 国产精品久久久久久亚洲av鲁大| 欧美日韩乱码在线| 欧美另类亚洲清纯唯美| 亚洲人成网站在线播放欧美日韩| 亚洲精品美女久久久久99蜜臀| 国产亚洲av嫩草精品影院| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 亚洲无线在线观看| 老司机午夜十八禁免费视频| 欧美 亚洲 国产 日韩一| 国产在线观看jvid| 色精品久久人妻99蜜桃| 日韩大尺度精品在线看网址 | 国产伦人伦偷精品视频| 桃色一区二区三区在线观看| 久久影院123| 久热爱精品视频在线9| 国产午夜福利久久久久久| 一区二区日韩欧美中文字幕| 男人的好看免费观看在线视频 | 男人操女人黄网站| 国产成人免费无遮挡视频| 精品电影一区二区在线| 99国产综合亚洲精品| 日本撒尿小便嘘嘘汇集6| 亚洲av成人一区二区三| 村上凉子中文字幕在线| 亚洲激情在线av| 91麻豆av在线| 老司机福利观看| 一级毛片女人18水好多| 午夜两性在线视频| 亚洲少妇的诱惑av| 午夜影院日韩av| 99国产极品粉嫩在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲中文字幕一区二区三区有码在线看 | 色尼玛亚洲综合影院| 欧美日本视频| 欧美亚洲日本最大视频资源| avwww免费| 制服丝袜大香蕉在线| 99香蕉大伊视频| 成熟少妇高潮喷水视频| 亚洲欧洲精品一区二区精品久久久| 久久国产精品人妻蜜桃| 熟女少妇亚洲综合色aaa.| 岛国在线观看网站| 亚洲精品中文字幕在线视频| 国产私拍福利视频在线观看| 黄色视频不卡| 91av网站免费观看| 高清毛片免费观看视频网站| 日本一区二区免费在线视频| 国产高清videossex| 国产亚洲精品一区二区www| 757午夜福利合集在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲熟妇熟女久久| 亚洲av美国av| 国产成人精品在线电影| 18美女黄网站色大片免费观看| 搡老妇女老女人老熟妇| 久99久视频精品免费| 久久久国产成人精品二区| svipshipincom国产片| 日日摸夜夜添夜夜添小说| 一区二区三区国产精品乱码| 97人妻天天添夜夜摸| 亚洲精品中文字幕在线视频| 欧美av亚洲av综合av国产av| 久久久久久国产a免费观看| 亚洲aⅴ乱码一区二区在线播放 | 久久久久精品国产欧美久久久| 午夜日韩欧美国产| 精品国产超薄肉色丝袜足j| 两个人免费观看高清视频| 久久久久久大精品| 桃红色精品国产亚洲av| 久久精品成人免费网站| av超薄肉色丝袜交足视频| 别揉我奶头~嗯~啊~动态视频| 黄色a级毛片大全视频| 精品日产1卡2卡| 黄色a级毛片大全视频| 精品久久久久久久久久免费视频| 亚洲人成网站在线播放欧美日韩| 亚洲自拍偷在线| 国产午夜精品久久久久久| 宅男免费午夜| 亚洲中文av在线| 久久中文字幕一级| 在线观看免费午夜福利视频| 亚洲在线自拍视频| 不卡av一区二区三区| 亚洲 国产 在线| 深夜精品福利| 露出奶头的视频| 欧美成人一区二区免费高清观看 | 久久中文看片网| 久久精品影院6| a在线观看视频网站| 欧美国产日韩亚洲一区| 在线观看舔阴道视频| 午夜福利高清视频| 久久狼人影院| 久久久水蜜桃国产精品网| 亚洲性夜色夜夜综合| 在线免费观看的www视频| 十八禁网站免费在线| 精品国产乱码久久久久久男人| 精品欧美国产一区二区三| 久久精品成人免费网站| 亚洲一区高清亚洲精品| 国产亚洲精品av在线| 国产99白浆流出| 国产精品久久久av美女十八| 久久久久久大精品| 欧美日韩一级在线毛片| 亚洲精华国产精华精| 国产成人系列免费观看| 国产主播在线观看一区二区| 国产一区二区三区在线臀色熟女| 亚洲人成电影观看| 成人国语在线视频| 在线观看舔阴道视频| 黄色 视频免费看| 自线自在国产av| 在线观看一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产色视频综合| 成人三级做爰电影| 91老司机精品| 色老头精品视频在线观看| 久久精品91蜜桃| 欧美在线黄色| 天天添夜夜摸| 制服诱惑二区| 久久热在线av| 两性午夜刺激爽爽歪歪视频在线观看 | 啦啦啦 在线观看视频| 天天添夜夜摸| 午夜免费激情av| 极品人妻少妇av视频| 欧美色欧美亚洲另类二区 | 亚洲成人免费电影在线观看| 色哟哟哟哟哟哟| 女人高潮潮喷娇喘18禁视频| 午夜福利一区二区在线看| 亚洲国产毛片av蜜桃av| 天堂动漫精品| 两个人视频免费观看高清| 久久午夜亚洲精品久久| 可以免费在线观看a视频的电影网站| 级片在线观看| 免费高清视频大片| 琪琪午夜伦伦电影理论片6080| 亚洲精品国产精品久久久不卡| 嫩草影院精品99| 黑丝袜美女国产一区| 后天国语完整版免费观看| 日韩一卡2卡3卡4卡2021年| 脱女人内裤的视频| 亚洲三区欧美一区| 精品人妻在线不人妻| 天堂√8在线中文| 日日干狠狠操夜夜爽| 亚洲av电影不卡..在线观看| 男男h啪啪无遮挡| 久久香蕉精品热| 日本精品一区二区三区蜜桃| 一个人观看的视频www高清免费观看 | 欧美成人免费av一区二区三区| 日本三级黄在线观看| 亚洲欧洲精品一区二区精品久久久| 丁香六月欧美| 97人妻天天添夜夜摸| 免费在线观看影片大全网站| 一区二区日韩欧美中文字幕| 麻豆国产av国片精品| 国产1区2区3区精品| 精品少妇一区二区三区视频日本电影| 日本免费一区二区三区高清不卡 | 欧美乱妇无乱码| 高清在线国产一区| 亚洲一区中文字幕在线| 伊人久久大香线蕉亚洲五| 久久精品影院6| 制服诱惑二区| 精品久久久久久成人av| 日日爽夜夜爽网站| 老熟妇乱子伦视频在线观看| 午夜影院日韩av| 91av网站免费观看| 亚洲第一电影网av| 一边摸一边抽搐一进一出视频| 美女免费视频网站| 欧美日韩乱码在线| 国产欧美日韩一区二区精品| www.熟女人妻精品国产| 亚洲专区国产一区二区| 黑人巨大精品欧美一区二区蜜桃| 91成年电影在线观看| 精品国产一区二区久久| 国产激情久久老熟女| 亚洲av成人一区二区三| 亚洲一区二区三区色噜噜| 国产精品久久久久久精品电影 | 宅男免费午夜| a级毛片在线看网站| 欧美另类亚洲清纯唯美| 曰老女人黄片| 乱人伦中国视频| 国产成人精品在线电影| 91麻豆av在线| 国产精品免费视频内射| 18禁国产床啪视频网站| 欧美日韩一级在线毛片| 搞女人的毛片| 亚洲视频免费观看视频| 欧美色视频一区免费| av电影中文网址| avwww免费| 波多野结衣一区麻豆| 久久久久久久久中文| 99re在线观看精品视频| 国产色视频综合| 搡老熟女国产l中国老女人| 日本五十路高清| 色精品久久人妻99蜜桃| 啦啦啦韩国在线观看视频| 在线视频色国产色| 久久天躁狠狠躁夜夜2o2o| 色尼玛亚洲综合影院| 免费在线观看影片大全网站| 国产一区二区在线av高清观看| www国产在线视频色| 一本综合久久免费| 涩涩av久久男人的天堂| 亚洲熟妇熟女久久| 国产亚洲欧美精品永久| 精品国产超薄肉色丝袜足j| 免费看美女性在线毛片视频| 校园春色视频在线观看| 午夜两性在线视频| 美女高潮喷水抽搐中文字幕| 国产精华一区二区三区| 看片在线看免费视频| 亚洲国产毛片av蜜桃av| 亚洲七黄色美女视频| 一区在线观看完整版| 午夜久久久久精精品| 亚洲欧美日韩另类电影网站| 久久精品国产亚洲av香蕉五月| 后天国语完整版免费观看| 黑人操中国人逼视频| 国产精品香港三级国产av潘金莲| 日本免费a在线| 国产亚洲欧美98| 精品国产国语对白av| 国产精品综合久久久久久久免费 | 国产亚洲精品久久久久5区| 色综合欧美亚洲国产小说| 欧美国产精品va在线观看不卡| 精品久久久久久,| 老鸭窝网址在线观看| 高清在线国产一区| 色老头精品视频在线观看| 免费看a级黄色片| videosex国产| 操美女的视频在线观看| 激情在线观看视频在线高清| 99热只有精品国产| 国产欧美日韩综合在线一区二区| 女性被躁到高潮视频| 日韩欧美三级三区| 色播在线永久视频| av在线播放免费不卡| 好男人在线观看高清免费视频 | 9191精品国产免费久久| 精品人妻1区二区| 一个人观看的视频www高清免费观看 | 国产一级毛片七仙女欲春2 | 两个人看的免费小视频| 免费在线观看完整版高清| 国产精品一区二区在线不卡| 99热只有精品国产| 精品国产乱子伦一区二区三区| 在线免费观看的www视频| 亚洲天堂国产精品一区在线| 深夜精品福利| 脱女人内裤的视频| 19禁男女啪啪无遮挡网站| 人妻丰满熟妇av一区二区三区| 国产91精品成人一区二区三区| 正在播放国产对白刺激| 在线免费观看的www视频| 最近最新中文字幕大全电影3 | 男女床上黄色一级片免费看| 国产在线精品亚洲第一网站| 在线永久观看黄色视频| 久久这里只有精品19| 亚洲国产精品合色在线| 午夜影院日韩av| 国产成人精品在线电影| 女人被狂操c到高潮| 超碰成人久久| 午夜福利一区二区在线看| 免费久久久久久久精品成人欧美视频| av视频免费观看在线观看| 国内久久婷婷六月综合欲色啪| 男人舔女人下体高潮全视频| 在线播放国产精品三级| 1024香蕉在线观看| 精品国产乱子伦一区二区三区| 亚洲精品美女久久久久99蜜臀| 精品国产乱子伦一区二区三区| 国产欧美日韩精品亚洲av| 国产99久久九九免费精品| 黄色视频不卡| 国产亚洲精品av在线| 俄罗斯特黄特色一大片| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 久久香蕉精品热| 一区二区三区精品91| 亚洲欧美激情在线| 制服诱惑二区| 日韩精品青青久久久久久| 老司机午夜十八禁免费视频| 国产在线观看jvid| 国内久久婷婷六月综合欲色啪| 久久天堂一区二区三区四区| 国产精品一区二区免费欧美| 欧美日本中文国产一区发布| 亚洲午夜理论影院| 亚洲七黄色美女视频| 激情在线观看视频在线高清| 日韩大码丰满熟妇| 亚洲成人久久性| 中文字幕人妻熟女乱码| 国产男靠女视频免费网站| 桃色一区二区三区在线观看| 一a级毛片在线观看| 亚洲avbb在线观看| 亚洲,欧美精品.| 久久人人爽av亚洲精品天堂| 在线天堂中文资源库| 一级毛片高清免费大全| 叶爱在线成人免费视频播放| 老熟妇乱子伦视频在线观看| 午夜免费激情av| 老司机靠b影院| 久99久视频精品免费| 亚洲精品在线观看二区| 精品日产1卡2卡| 两个人免费观看高清视频| 国产精品一区二区在线不卡| 亚洲视频免费观看视频| 最新在线观看一区二区三区| 老熟妇乱子伦视频在线观看| 99久久国产精品久久久| 国产午夜福利久久久久久| 日本免费a在线| 脱女人内裤的视频| 欧美黑人精品巨大| 嫁个100分男人电影在线观看| 亚洲专区字幕在线| 国产av一区二区精品久久| 男人操女人黄网站| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 亚洲午夜精品一区,二区,三区| 一二三四社区在线视频社区8| 精品第一国产精品| 精品一品国产午夜福利视频| 可以免费在线观看a视频的电影网站| 久久久久亚洲av毛片大全| 国产av又大| 天天一区二区日本电影三级 | 亚洲成a人片在线一区二区| 成人三级黄色视频| 午夜久久久久精精品| 变态另类丝袜制服| 又黄又粗又硬又大视频| 在线观看午夜福利视频| 午夜福利在线观看吧| 亚洲成人免费电影在线观看| 国产麻豆成人av免费视频| 在线视频色国产色| av天堂久久9| 亚洲熟妇中文字幕五十中出| 色综合婷婷激情| 日本a在线网址| 国产精品一区二区精品视频观看| 99在线视频只有这里精品首页| 精品人妻1区二区| 校园春色视频在线观看| 99国产精品免费福利视频| 亚洲av片天天在线观看| 国产精品久久久久久人妻精品电影| 一区二区三区激情视频| 日韩一卡2卡3卡4卡2021年| 日韩 欧美 亚洲 中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产一区二区在线av高清观看| 国产伦一二天堂av在线观看| 亚洲aⅴ乱码一区二区在线播放 | 搞女人的毛片| 在线观看午夜福利视频| 成人国产一区最新在线观看| 久久久久国产一级毛片高清牌| 日韩 欧美 亚洲 中文字幕| 久久国产精品人妻蜜桃| 嫩草影院精品99| 亚洲色图av天堂| 久久性视频一级片| 韩国精品一区二区三区| 国产精品久久电影中文字幕| 99在线人妻在线中文字幕| 又黄又粗又硬又大视频| 男人的好看免费观看在线视频 | 精品不卡国产一区二区三区| 99在线视频只有这里精品首页| 村上凉子中文字幕在线| 制服人妻中文乱码| 成人亚洲精品一区在线观看| 国产乱人伦免费视频| 久9热在线精品视频| 村上凉子中文字幕在线| 最新美女视频免费是黄的| 91成年电影在线观看| 欧美不卡视频在线免费观看 | 满18在线观看网站| 午夜两性在线视频| 99久久精品国产亚洲精品| 久久精品国产亚洲av香蕉五月| 亚洲色图综合在线观看| 国产成人精品久久二区二区91| 精品一区二区三区视频在线观看免费| 精品一区二区三区四区五区乱码| 两个人免费观看高清视频| 国产精品电影一区二区三区| 亚洲av熟女| 999精品在线视频| 欧美中文日本在线观看视频| 午夜福利,免费看| 欧美乱码精品一区二区三区| 99riav亚洲国产免费| 12—13女人毛片做爰片一| 精品久久蜜臀av无| www.999成人在线观看| 多毛熟女@视频| 长腿黑丝高跟| 欧美日韩瑟瑟在线播放| 国产熟女午夜一区二区三区| 美女午夜性视频免费| 50天的宝宝边吃奶边哭怎么回事| 午夜亚洲福利在线播放| 国产99白浆流出| 久久久国产成人精品二区| 在线免费观看的www视频| 咕卡用的链子| 一本久久中文字幕| 婷婷精品国产亚洲av在线| 国产av一区在线观看免费| 国产成人精品在线电影| 亚洲五月色婷婷综合| 很黄的视频免费| 国产精品久久久久久精品电影 | 丁香六月欧美| 国产视频一区二区在线看| 国产精品一区二区三区四区久久 | 成人亚洲精品一区在线观看| 99国产精品99久久久久| 99精品欧美一区二区三区四区| 最新美女视频免费是黄的| 99精品欧美一区二区三区四区| 黄片播放在线免费| 长腿黑丝高跟| 中文字幕高清在线视频| 国产精品九九99| 我的亚洲天堂| 亚洲精品粉嫩美女一区| 男男h啪啪无遮挡| 国内毛片毛片毛片毛片毛片| 亚洲最大成人中文| 老鸭窝网址在线观看| 亚洲欧美激情在线| 波多野结衣一区麻豆| 每晚都被弄得嗷嗷叫到高潮| 国产精品美女特级片免费视频播放器 | 国内久久婷婷六月综合欲色啪| 亚洲视频免费观看视频| 亚洲情色 制服丝袜| 丝袜美腿诱惑在线| 国产一区二区三区视频了| 成人三级做爰电影| 色综合站精品国产| 亚洲人成网站在线播放欧美日韩| 欧美黄色淫秽网站| 老汉色av国产亚洲站长工具| 岛国在线观看网站| 国产成人av教育| 午夜免费观看网址| 欧美绝顶高潮抽搐喷水| 中文字幕人妻熟女乱码| av视频免费观看在线观看| 在线观看舔阴道视频| 青草久久国产| 国产真人三级小视频在线观看| 日日夜夜操网爽| 色播亚洲综合网| 亚洲专区中文字幕在线| 日韩欧美国产在线观看| 叶爱在线成人免费视频播放| 999精品在线视频| 欧美中文综合在线视频| 大陆偷拍与自拍| 亚洲精品国产区一区二| www.精华液| 黄色视频,在线免费观看| 999久久久精品免费观看国产| 久久久久国产一级毛片高清牌| 国产一区二区在线av高清观看| 黄色成人免费大全| 亚洲一区高清亚洲精品| 99在线人妻在线中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 日本vs欧美在线观看视频| 99精品欧美一区二区三区四区| 一边摸一边抽搐一进一出视频| x7x7x7水蜜桃| 国产精品二区激情视频| 女性被躁到高潮视频| 精品无人区乱码1区二区| 人人妻人人爽人人添夜夜欢视频| 亚洲午夜精品一区,二区,三区| 很黄的视频免费| 亚洲久久久国产精品| 国产精品国产高清国产av| 国内精品久久久久久久电影| 中文字幕人妻丝袜一区二区| 国产乱人伦免费视频| 国产区一区二久久| 国产一区二区在线av高清观看| 欧美激情久久久久久爽电影 | 成人三级黄色视频| 日日爽夜夜爽网站| 高清毛片免费观看视频网站| 中文字幕人成人乱码亚洲影| 99久久99久久久精品蜜桃| 久久性视频一级片| 国产高清激情床上av|