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

    汽輪機末級葉片氣流激振力的計算方法

    2021-05-12 04:34:32葛存飛劉慶龍徐志強
    關(guān)鍵詞:動葉激振力末級

    葛存飛,劉慶龍,徐志強,陳 榴

    (1.杭州中能汽輪動力有限公司汽輪機研究所,杭州 310018;2.無錫市厚德自動化儀表有限公司,江蘇江陰 214413;3.上海理工大學能源與動力工程學院,上海 200093)

    葉輪機械是能源轉(zhuǎn)化為電力最重要的動力機械,葉片作為實現(xiàn)動能轉(zhuǎn)化為機械能的重要做功部件,其安全性關(guān)系到整個葉輪機械的順利運行。葉片長期運行在極為復雜的工作條件下,并承受氣流沖擊以及旋轉(zhuǎn)機械振動等一系列力,因此葉片的振動安全性分析一直是葉片設計研發(fā)的難題之一。在蒸汽輪機中,特別是低壓級組中的長扭葉片設計最為復雜。從國內(nèi)數(shù)據(jù)看,末級葉片損壞占葉片故障統(tǒng)計總數(shù)的52.94%,其中60%~80%的葉片損壞的原因是振動疲勞損傷或者斷裂,因此研究長扭葉片振動安全性能尤為重要。

    葉片激振力是指使葉片產(chǎn)生振動的脈沖氣流力等的作用力。研究激振力就是研究其發(fā)生機制,通過分析其大小和頻率,進而得到葉片振動疲勞特性。葉輪機械內(nèi)部流動本質(zhì)上是三維非定常復雜流動,非定常作用將直接在葉片上產(chǎn)生非定常負荷,引起葉片振動。早期的非定常研究主要依靠試驗和理論,很多學者把試驗結(jié)果與理論進行了對比,發(fā)現(xiàn)了試驗得到的結(jié)果與當時的一元理論相去甚遠[1]。隨著氣動理論的完善和計算流體力學(CFD)的發(fā)展,很多學者通過理論分析和數(shù)值模擬的方法對葉片非定常的受力進行分析研究[2-6],并且得出了很多具有工程意義的結(jié)論。

    目前,對于葉片受力的研究著重在各種不同的因素影響下的葉片振動激振力計算方法的研究[7]。但是這些計算方法只能得到非定常激振力的頻率,卻無法準確得到非定常激振力的數(shù)值。在理論推導中一般采用放大系數(shù)等的經(jīng)驗參數(shù)來估算。工程應用中針對具體葉片具體機組,經(jīng)驗參數(shù)波動較大,導致沒法得到精確結(jié)果。特別是末級扭葉片,受到排氣管或者排缸的徑向排汽結(jié)構(gòu)影響較大,末級扭葉片的研發(fā)通常需要攜帶排缸進行三維非定常分析[8]。這項數(shù)值分析任務需要末級動靜整圈葉柵和排缸耦合建模,規(guī)模龐大,研發(fā)工作可能需要幾個月甚至幾年,在實際設計研發(fā)中很少采用。多數(shù)研發(fā)和校核工作是在理論計算結(jié)合定常分析的基礎上進行的,所以確定性較差,對葉片疲勞壽命評估較為吃力。因此,有效地預估汽輪機中的非定常激振力特性,對提高汽輪機葉片的優(yōu)化設計具有重要的意義。本文在結(jié)合 TBR(transient blade row)[9]計算和定常 FFT(fast fourier transformation)變換的基礎上,提出了一種有效的非定常激振力的計算方法。

    1 分析方法和理論

    末級動葉受到的激振力有很多種,末級來說一般不存在部分進汽,故除開由于加工制造和安裝裝配引起的周向結(jié)構(gòu)不均勻以外,影響末級最大的兩個激振力因素是靜葉噴嘴出口的尾跡不均勻和排缸徑向排汽導致的周向不均勻。由于大功率發(fā)電用汽輪機的末級動葉片頻率基本上低于150 Hz。若靜葉噴嘴只數(shù)為60只,則靜葉噴嘴出口尾跡的不均勻激振力頻率至少應該為50 Hz×60=3 000 Hz。該頻率遠遠大于150 Hz,因此靜葉噴嘴出口高頻激振力對末級影響很小。對于末級來說,主要的研究對象著重于排缸結(jié)構(gòu)引起的周向不均勻。

    凝汽式汽輪機的背壓一般由凝汽器來維持穩(wěn)定。排缸一般都設計為徑向排汽,導致排缸在周向的靠近排汽位置側(cè)壓力較低,相應的另外一側(cè)壓力較高。末級動葉片在一個旋轉(zhuǎn)周期內(nèi)將會經(jīng)歷一次壓力震蕩,每個旋轉(zhuǎn)周期所經(jīng)歷的壓力震蕩一致,因此引入傅里葉變換。傅里葉變換的思想可以簡述為,如果一個時域中的周期性函數(shù),則可以轉(zhuǎn)化為一個繞著中心旋轉(zhuǎn)的柱坐標系上的頻域函數(shù)。其中的θ角即該頻域函數(shù)的相位角,其旋轉(zhuǎn)的時間函數(shù)t自變量轉(zhuǎn)變?yōu)棣兀魏我粋€滿足狄里赫萊條件的周期性函數(shù)可以用傅里葉級數(shù)進行展開。

    傅里葉變換公式:

    傅里葉變換實質(zhì)是時域轉(zhuǎn)化為頻域函數(shù)的過程,將原來關(guān)于時間t的時域函數(shù)f(t)轉(zhuǎn)變?yōu)轭l率ω的頻域函數(shù)F(ω)。

    傅里葉逆變換公式:

    傅里葉逆變換的實質(zhì)是頻域轉(zhuǎn)化為時域函數(shù)的過程。

    在末級動靜葉片與排缸組成的系統(tǒng)中,排缸對末級動葉的影響主要體現(xiàn)在壓力周向不均勻性。在末級扭葉片的激振力分析中引入傅里葉逆變換的思想:采用三維定常方法得到排缸參數(shù)(主要是壓力)沿周向分布,通過傅里葉逆變換,將動葉出口周向壓力分布轉(zhuǎn)化為關(guān)于時間的壓力脈動函數(shù),作為單流道出口邊界條件,實現(xiàn)單流道模擬整圈排缸對動葉片的影響。末級動葉出口壓力隨時間的參數(shù)脈動是以汽輪機旋轉(zhuǎn)一周為周期而且實際運行連續(xù),可劃分為有限個單調(diào)區(qū)間(圈數(shù)),壓力函數(shù)為有限值,且絕對可積,滿足狄里赫萊條件。

    將末級動葉出口壓力P,沿周向按傅里葉級數(shù)展開,則有

    式中:ω是旋轉(zhuǎn)角速度;ˉP是末級動葉出口壓力整圈按時間的平均值;k是末級動葉出口壓力脈動階次;PK是第k階末級動葉出口壓力值;φK是排汽壓力值的相位角。

    動葉出口壓力脈動變?yōu)闀r域函數(shù)后,結(jié)合TBR(transient blade row)方法,利用單流道,采用非定常方法計算背壓不均勻性對動葉流場的影響,而不再使用動靜葉只數(shù)的最小公倍數(shù)流道,令網(wǎng)格規(guī)模成倍縮減,計算量大幅降低,計算效率增加。

    2 計算流程和實例

    圖1給出了計算流程圖。動葉出口壓力周向分布,需要通過整圈加排缸的三維定常耦合計算獲得。相對于非定常的耦合計算,其計算量及所需的計算資源都較少。一般的非定常計算需要多個時間步,除去前面參數(shù)波動未趨于穩(wěn)定的時間步以外,至少需要5~6個穩(wěn)定周期,而每個周期都需要多個時間步來計算。楊杰等[10]研究兩級渦輪葉片的定常與非定常比較發(fā)現(xiàn),同樣的網(wǎng)格和計算域,非定常綜合成本大約是定常計算成本的176倍。按照一般的非定常耦合計算,本文計算的綜合成本可能超過200倍。而采用本文提出的在定常FFT變換的基礎上結(jié)合TBR計算,則可以節(jié)約大量的時間和綜合成本。

    圖1 計算流程框圖

    以某汽輪機組N50-8.83的末級葉片激振力計算為例,其末級葉片長度是668 mm,設計排汽壓力為8 kPa。假定定常計算結(jié)果經(jīng)過FFT變換以后,整圈動葉出口壓力隨時間的脈動函數(shù)P只有傅里葉一階級數(shù),高階級數(shù)均為零。此時由于僅有一個Sin函數(shù),其相位只與起始位置有關(guān),而和其他波形沒有關(guān)系,故可令ψ1=0,且可以假定脈動幅值p1為的1%。則有

    圖2展示了作為P的壓力脈動隨時間的變化圖表,實際應用中分離出來末級出口壓力可能存在二階甚至高階級數(shù),可以查看幅值大小確定影響,如果影響較小適當予以忽略。

    圖2 壓力-時間脈動

    將式(4)的壓力值放入單流道的模型中作為出口條件,末級進出口參數(shù)展示在表1中。其單流道模型,僅有一個動靜葉組流道,網(wǎng)格規(guī)模通過無關(guān)性驗證。靜葉片模型網(wǎng)格n數(shù)為20萬,動葉網(wǎng)格數(shù)為25萬。計算采用SST(shear stress transport)湍流模型,二階迎風空間格式計算,瞬態(tài)流場分析采用TBR模塊,瞬態(tài)計算時間步長為1e-5 s。

    表1 末級邊界條件及非定常計算總時

    3 結(jié)果分析

    整個瞬態(tài)計算耗費資源見表2,可以看到對瞬態(tài)來說,由于僅采用單流道,其內(nèi)存和CPU使用量,運算總時間比較少,綜合成本下降較快。以綜合成本計算(計算總和成本為CPU內(nèi)存和時間乘積),單流道TBR計算資源是全環(huán)非定常計算的2.3‰,不計內(nèi)存資源,計算速度提升約8倍。

    表2 末級單流道瞬態(tài)計算資源消耗

    通對瞬態(tài)計算的結(jié)果進行后處理,提取壓力已經(jīng)穩(wěn)定的一個周期進行分析。圖3給出了4個不同時刻的動葉壓力云圖??梢园l(fā)現(xiàn)在t=0.025 s時刻,壓力均勻,出口面型底壓力出現(xiàn)了低壓區(qū)域,可能是瞬時出現(xiàn)了一定的回流。而在t=0.027 5 s時入口面型底出現(xiàn)壓力峰值區(qū)域,呈長條形出現(xiàn)在壓力面的頭緣附近,其出氣邊型底壓力提升到背壓平均值上下。

    圖3 不同時刻的動葉表面壓力云圖

    在t=0.03 s時刻,壓力峰值區(qū)域向葉片吸力面移動,尾緣附近出現(xiàn)同圖3(a)的小壓力情況,而在t=0.032 5 s時刻,型底背壓出現(xiàn)壓力峰值,甚至其大小超過了入口壓力,可能是背壓調(diào)整和越過尾跡區(qū)域后突變引起的。可以從幾個圖中看到:壓力隨時間的變化呈周期性且變化較為明顯。瞬態(tài)葉片的受力和穩(wěn)態(tài)受力結(jié)果很不一樣,因此分析葉片詳細受力的時候必須通過瞬態(tài)流場計算再結(jié)合瞬態(tài)結(jié)構(gòu)分析才能得到最接近真實的結(jié)果。

    為了更好地說明葉片受力隨時間的變化,對葉片表面力做FFT變換。在耦合計算中,直接提取葉片表面壓力隨時間的變化數(shù)值導入瞬態(tài)結(jié)構(gòu)分析中。在本文中由于不再具體進行瞬態(tài)結(jié)構(gòu)耦合分析,直接提取了葉片表面力(Force)作為研究對象。葉片表面力其實分為3個分量,本文提取了其中2個z和y、z向是軸向力,y是切向力(周向力)。圖4展示了2個方向力隨著時間的變化,其中圖4(a)顯示了葉片的軸向受力,僅提取了其中一個時間周期。可以看到:葉片受力的周期性比較明顯,且幅值也相應變化。圖4(c)顯示的結(jié)果幾乎與圖4(a)相同,唯一區(qū)別是幅值和相位有所差別。圖4(b)和圖4(d)分別顯示了軸向和周向受力葉片的FFT變換圖,其縱坐標是振幅(magnitude)。這兩圖中橫坐標0 Hz位置均出現(xiàn)非常高的幅值,可認為這就是葉片受力的傅里葉展開后的第0階級數(shù)的幅值,即葉片受力的時間平均值。Fz的時間平均值約為33 N,F(xiàn)y的時間平均值約為51 N。同時在0~200 Hz之間出現(xiàn)了一個較小的峰值,其頻率在50~100之間??烧J為這個是由于背壓變化引起的傅里葉級數(shù),但由于背壓變化可能影響級前壓力和靜葉的出口壓力分布,因此其產(chǎn)生的對葉片表面受力的激振頻率可能發(fā)生遷移。

    圖4 軸向和周向葉片表面力及FFT變換

    由于本次模擬的背壓一階幅值僅取了1%,其振幅較小。剩余兩個較大的峰值出現(xiàn)在1 000 Hz和2 000 Hz附近,這應該是由于靜葉噴嘴出口的尾跡不均勻?qū)е碌膭尤~壓力脈動頻率,由于背壓脈動導致了高階倍數(shù)頻率的出現(xiàn)。不過這個頻率距離葉片的固有頻率相去甚遠,不能激發(fā)葉片共振。

    背壓頻域變化的加入導致了原有動葉片受力發(fā)生改變,且可能在動靜葉片之間也出現(xiàn)了壓力波動,因此對于該問題還需要進行更加詳細的分析。本文方法已被驗證是可行的,且計算時間和占用的系統(tǒng)資源遠遠小于傳統(tǒng)的非定常計算。

    4 結(jié)論

    1)FFT將空間轉(zhuǎn)換到時間,使定常數(shù)據(jù)可用于非定常計算。

    2)利用FFT變換結(jié)合TBR方法實現(xiàn)瞬態(tài)單流道非定常計算方法模擬實際問題是可行的,TBR方法的計算綜合成本是全環(huán)非定常計算的2.3‰,不計內(nèi)存消耗,計算速度提升約8倍。

    3)通過實例驗證了該方法的可行性和效率,并給出了計算結(jié)果。結(jié)果表明:其背壓壓力波動的存在影響了靜葉噴嘴出口壓力脈動大小和頻率。

    猜你喜歡
    動葉激振力末級
    末級壓出室水力結(jié)構(gòu)對多級離心泵水力性能的影響
    運載火箭末級離軌控制策略優(yōu)化
    M701DA燃機壓氣機動葉鎖鍵的配合方式研究
    超臨界機組鍋爐末級過熱器管爆管分析
    發(fā)電廠汽輪機振動異常增大的原因分析
    DF100A發(fā)射機末級電子管的計算與分析
    電子制作(2018年8期)2018-06-26 06:43:32
    液壓磚機激振力計算方法及影響因素
    軸-艇耦合系統(tǒng)的力傳遞特性分析
    試析引風機動葉故障原因及處理措施
    雙級動調(diào)引風機變頻改造節(jié)能安全運行實踐
    中文欧美无线码| 国产亚洲一区二区精品| 日本一本二区三区精品| 全区人妻精品视频| 亚洲欧洲国产日韩| 在线观看免费高清a一片| 99热6这里只有精品| 夜夜爽夜夜爽视频| 国内精品美女久久久久久| 精品国产乱码久久久久久小说| 国产亚洲最大av| 成年人午夜在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 99re6热这里在线精品视频| freevideosex欧美| 国产国拍精品亚洲av在线观看| 国语对白做爰xxxⅹ性视频网站| 99热这里只有精品一区| 亚洲最大成人中文| 日韩不卡一区二区三区视频在线| 天天躁日日操中文字幕| 高清在线视频一区二区三区| 国产有黄有色有爽视频| av免费在线看不卡| 99久国产av精品国产电影| 高清在线视频一区二区三区| 高清午夜精品一区二区三区| 亚洲av免费在线观看| 久久午夜福利片| 国产欧美日韩一区二区三区在线 | 看十八女毛片水多多多| av天堂中文字幕网| 精品少妇黑人巨大在线播放| 99久久九九国产精品国产免费| .国产精品久久| 制服丝袜香蕉在线| 午夜日本视频在线| 久久99精品国语久久久| 久久精品久久久久久噜噜老黄| 国产伦精品一区二区三区视频9| 亚洲av福利一区| 国产精品爽爽va在线观看网站| 欧美成人a在线观看| 亚洲精品日本国产第一区| 韩国高清视频一区二区三区| 别揉我奶头 嗯啊视频| 另类亚洲欧美激情| 熟女电影av网| 久久久成人免费电影| 亚洲婷婷狠狠爱综合网| 亚洲在线观看片| 日本黄大片高清| 国产精品国产三级国产av玫瑰| videossex国产| 国产乱人偷精品视频| 亚洲经典国产精华液单| 最近最新中文字幕免费大全7| 欧美成人一区二区免费高清观看| 成人美女网站在线观看视频| 中文字幕av成人在线电影| 久久精品综合一区二区三区| 免费看av在线观看网站| 三级男女做爰猛烈吃奶摸视频| 亚洲怡红院男人天堂| 好男人视频免费观看在线| 国产精品久久久久久精品古装| 成人鲁丝片一二三区免费| 日本爱情动作片www.在线观看| 久久99热这里只频精品6学生| 校园人妻丝袜中文字幕| 18禁裸乳无遮挡免费网站照片| 建设人人有责人人尽责人人享有的 | 精品国产一区二区三区久久久樱花 | 春色校园在线视频观看| 少妇人妻精品综合一区二区| 亚洲成人一二三区av| 新久久久久国产一级毛片| kizo精华| 国产男女超爽视频在线观看| 日韩一本色道免费dvd| 少妇人妻 视频| 精品一区在线观看国产| 女人十人毛片免费观看3o分钟| 在线 av 中文字幕| 亚洲av不卡在线观看| 中文字幕久久专区| 精品少妇黑人巨大在线播放| 一级毛片黄色毛片免费观看视频| 欧美精品一区二区大全| 亚洲精品,欧美精品| 免费观看a级毛片全部| 久久久久久久久大av| av在线天堂中文字幕| 天美传媒精品一区二区| 在线免费十八禁| 在线亚洲精品国产二区图片欧美 | 熟女电影av网| 亚洲欧美日韩无卡精品| 精品午夜福利在线看| 色视频在线一区二区三区| videossex国产| av女优亚洲男人天堂| 久久人人爽人人片av| 亚洲国产成人一精品久久久| 新久久久久国产一级毛片| 国产精品人妻久久久久久| 国产精品一区二区在线观看99| 国产乱人视频| 国产成人91sexporn| av国产精品久久久久影院| 熟女人妻精品中文字幕| 51国产日韩欧美| 精品久久久久久电影网| 夫妻性生交免费视频一级片| 日本免费在线观看一区| 夜夜看夜夜爽夜夜摸| 午夜福利视频1000在线观看| 成人二区视频| 老师上课跳d突然被开到最大视频| 国产高清三级在线| 99视频精品全部免费 在线| 免费看a级黄色片| 国产精品成人在线| 国产日韩欧美在线精品| 人妻 亚洲 视频| 欧美激情国产日韩精品一区| 亚洲国产高清在线一区二区三| 又大又黄又爽视频免费| 午夜福利高清视频| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久人人人人人人| 91久久精品电影网| 一级毛片 在线播放| 久久影院123| 嫩草影院入口| 日韩亚洲欧美综合| 精品久久久久久久人妻蜜臀av| 国产av码专区亚洲av| 26uuu在线亚洲综合色| 视频中文字幕在线观看| 美女cb高潮喷水在线观看| 欧美 日韩 精品 国产| 国产综合懂色| 久久精品国产亚洲av涩爱| 天堂中文最新版在线下载 | 欧美日韩国产mv在线观看视频 | 亚洲精品乱码久久久v下载方式| 欧美人与善性xxx| 美女脱内裤让男人舔精品视频| 欧美亚洲 丝袜 人妻 在线| 亚洲内射少妇av| 日韩 亚洲 欧美在线| 免费看av在线观看网站| 可以在线观看毛片的网站| 成年版毛片免费区| 亚洲精品日韩av片在线观看| 2018国产大陆天天弄谢| 亚洲国产最新在线播放| 又爽又黄a免费视频| 国产精品不卡视频一区二区| av.在线天堂| 日韩中字成人| 听说在线观看完整版免费高清| 亚洲国产最新在线播放| 国产在线一区二区三区精| 日韩,欧美,国产一区二区三区| 亚洲av国产av综合av卡| 亚洲精品国产成人久久av| 亚洲综合精品二区| 女的被弄到高潮叫床怎么办| 午夜免费观看性视频| www.av在线官网国产| 男人舔奶头视频| 色吧在线观看| 美女内射精品一级片tv| 丝袜脚勾引网站| 国产日韩欧美亚洲二区| 伦理电影大哥的女人| 日韩一本色道免费dvd| 亚洲国产精品专区欧美| 一二三四中文在线观看免费高清| 午夜激情福利司机影院| 男人添女人高潮全过程视频| 我的老师免费观看完整版| 亚洲精品乱久久久久久| 免费观看无遮挡的男女| 性色av一级| 久久ye,这里只有精品| 男人爽女人下面视频在线观看| 一级毛片aaaaaa免费看小| 天堂中文最新版在线下载 | 深爱激情五月婷婷| 狠狠精品人妻久久久久久综合| 一个人看视频在线观看www免费| 亚洲av中文av极速乱| 成人国产麻豆网| 97热精品久久久久久| 狠狠精品人妻久久久久久综合| 亚洲欧美成人精品一区二区| 蜜臀久久99精品久久宅男| 国产精品三级大全| 午夜爱爱视频在线播放| 高清在线视频一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲一级一片aⅴ在线观看| 性色avwww在线观看| 国产91av在线免费观看| av在线亚洲专区| 亚洲成人中文字幕在线播放| 青青草视频在线视频观看| 麻豆精品久久久久久蜜桃| 黄色欧美视频在线观看| 欧美精品一区二区大全| 99九九线精品视频在线观看视频| 国产人妻一区二区三区在| 色5月婷婷丁香| 亚洲天堂av无毛| 色视频在线一区二区三区| 极品教师在线视频| 久久鲁丝午夜福利片| av播播在线观看一区| 丰满乱子伦码专区| 在线免费十八禁| av在线老鸭窝| 波野结衣二区三区在线| 黄片wwwwww| 一区二区三区四区激情视频| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验| 久久精品人妻少妇| 国产男女内射视频| 热99国产精品久久久久久7| 久久精品国产a三级三级三级| 久久鲁丝午夜福利片| 久久国产乱子免费精品| 观看美女的网站| 蜜桃亚洲精品一区二区三区| 少妇高潮的动态图| 日韩不卡一区二区三区视频在线| 久久韩国三级中文字幕| 国产毛片在线视频| 精品久久久久久久久av| 国产成人精品一,二区| 欧美日韩精品成人综合77777| 国产 精品1| h日本视频在线播放| 久久久精品94久久精品| 精品一区二区三卡| 少妇人妻久久综合中文| 干丝袜人妻中文字幕| 1000部很黄的大片| 亚洲高清免费不卡视频| 久久鲁丝午夜福利片| 男女边吃奶边做爰视频| 欧美97在线视频| 精品国产三级普通话版| 美女视频免费永久观看网站| 国内少妇人妻偷人精品xxx网站| 18禁在线无遮挡免费观看视频| 日日摸夜夜添夜夜添av毛片| 久久久久久久午夜电影| tube8黄色片| 欧美bdsm另类| 啦啦啦啦在线视频资源| 边亲边吃奶的免费视频| 美女cb高潮喷水在线观看| www.av在线官网国产| 久久ye,这里只有精品| 久久久久久久亚洲中文字幕| 免费播放大片免费观看视频在线观看| 久久精品国产亚洲网站| 国内少妇人妻偷人精品xxx网站| 能在线免费看毛片的网站| 不卡视频在线观看欧美| 99九九线精品视频在线观看视频| 国产成人午夜福利电影在线观看| 久久久久久国产a免费观看| 亚洲精品一区蜜桃| 日韩精品有码人妻一区| 91精品国产九色| 久久精品人妻少妇| 国产美女午夜福利| 1000部很黄的大片| 国产老妇伦熟女老妇高清| 国产黄片美女视频| 久久6这里有精品| 在线播放无遮挡| 一边亲一边摸免费视频| 国产综合懂色| 免费电影在线观看免费观看| 日日啪夜夜爽| 亚洲国产精品专区欧美| 久久精品国产a三级三级三级| 久久国产乱子免费精品| 又爽又黄无遮挡网站| 欧美精品人与动牲交sv欧美| 在现免费观看毛片| 寂寞人妻少妇视频99o| 日韩欧美精品免费久久| 99精国产麻豆久久婷婷| videossex国产| 九草在线视频观看| 久久久久久久精品精品| 日本熟妇午夜| av在线播放精品| 97超碰精品成人国产| 亚洲欧洲日产国产| 免费人成在线观看视频色| 中国三级夫妇交换| 国产综合懂色| 亚洲精品乱久久久久久| 日本一本二区三区精品| 亚洲精品色激情综合| 伊人久久国产一区二区| 黄片wwwwww| 亚洲成人一二三区av| 日产精品乱码卡一卡2卡三| 熟女人妻精品中文字幕| 久久精品人妻少妇| 午夜激情福利司机影院| freevideosex欧美| 青春草国产在线视频| 久久精品人妻少妇| 国产中年淑女户外野战色| 日韩欧美 国产精品| 亚洲精品国产av成人精品| 网址你懂的国产日韩在线| 啦啦啦中文免费视频观看日本| 国产成人午夜福利电影在线观看| 日韩强制内射视频| 国产伦理片在线播放av一区| 嫩草影院精品99| 成人午夜精彩视频在线观看| 国产一区二区亚洲精品在线观看| 80岁老熟妇乱子伦牲交| 亚洲精品国产色婷婷电影| 男女无遮挡免费网站观看| 国产精品人妻久久久影院| 在线观看人妻少妇| 国产一区亚洲一区在线观看| 欧美一级a爱片免费观看看| 成人国产av品久久久| 精品酒店卫生间| 男女边吃奶边做爰视频| 国产在线男女| 成人特级av手机在线观看| av专区在线播放| 一级片'在线观看视频| 成人美女网站在线观看视频| 久久久国产一区二区| 天堂网av新在线| 我的女老师完整版在线观看| 亚洲av成人精品一区久久| 中国三级夫妇交换| 黄色欧美视频在线观看| 一级a做视频免费观看| 联通29元200g的流量卡| 欧美另类一区| 久久亚洲国产成人精品v| 搡老乐熟女国产| 国产又色又爽无遮挡免| 久久人人爽av亚洲精品天堂 | 亚洲精品,欧美精品| 亚洲国产色片| 亚洲精品乱码久久久久久按摩| 中文乱码字字幕精品一区二区三区| 午夜福利在线观看免费完整高清在| 婷婷色综合www| 国产熟女欧美一区二区| 2021少妇久久久久久久久久久| 人体艺术视频欧美日本| 男插女下体视频免费在线播放| 欧美日韩一区二区视频在线观看视频在线 | 中文字幕制服av| 免费观看a级毛片全部| 国产精品久久久久久久久免| 少妇人妻一区二区三区视频| 国产爽快片一区二区三区| 成人国产麻豆网| 国产久久久一区二区三区| 日本wwww免费看| 亚洲国产精品国产精品| 免费电影在线观看免费观看| 精品一区二区三卡| 综合色av麻豆| 91狼人影院| 偷拍熟女少妇极品色| av在线亚洲专区| 黄色怎么调成土黄色| 久久鲁丝午夜福利片| 91久久精品国产一区二区成人| 九草在线视频观看| 日韩欧美精品免费久久| 在线看a的网站| 成年人午夜在线观看视频| 男的添女的下面高潮视频| 国产黄色视频一区二区在线观看| 午夜福利在线在线| 香蕉精品网在线| 又大又黄又爽视频免费| 观看免费一级毛片| 大片电影免费在线观看免费| kizo精华| 亚洲欧美日韩卡通动漫| 插逼视频在线观看| 下体分泌物呈黄色| 久久99精品国语久久久| 午夜精品国产一区二区电影 | 美女视频免费永久观看网站| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆| 久久人人爽人人爽人人片va| 亚洲欧美成人精品一区二区| 国产av国产精品国产| 白带黄色成豆腐渣| 禁无遮挡网站| 日韩av在线免费看完整版不卡| 麻豆精品久久久久久蜜桃| 久久国产乱子免费精品| 伊人久久国产一区二区| 国产国拍精品亚洲av在线观看| 成人午夜精彩视频在线观看| 综合色丁香网| av福利片在线观看| 国产老妇伦熟女老妇高清| 日韩av在线免费看完整版不卡| 久久99热6这里只有精品| 婷婷色综合大香蕉| 高清在线视频一区二区三区| 国产国拍精品亚洲av在线观看| 波多野结衣巨乳人妻| 在线观看国产h片| 亚洲最大成人中文| 一级毛片我不卡| 1000部很黄的大片| 久久99热6这里只有精品| 美女主播在线视频| 国产色爽女视频免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | eeuss影院久久| 尤物成人国产欧美一区二区三区| 2018国产大陆天天弄谢| 国产伦在线观看视频一区| av国产精品久久久久影院| av天堂中文字幕网| 国产精品一区www在线观看| 国产av国产精品国产| 亚洲最大成人中文| 最近手机中文字幕大全| 亚洲精品自拍成人| 午夜精品国产一区二区电影 | 亚洲国产成人一精品久久久| 亚洲一级一片aⅴ在线观看| 成人欧美大片| 嫩草影院新地址| 久久久久久久久久人人人人人人| av女优亚洲男人天堂| 国产黄频视频在线观看| 亚洲精品乱码久久久v下载方式| 亚洲精品色激情综合| 欧美变态另类bdsm刘玥| 亚洲欧美日韩卡通动漫| 国产午夜福利久久久久久| 亚洲人与动物交配视频| 看非洲黑人一级黄片| 国产综合懂色| 久久ye,这里只有精品| 九草在线视频观看| 日本欧美国产在线视频| 国产黄频视频在线观看| 又爽又黄a免费视频| 日本猛色少妇xxxxx猛交久久| 久久亚洲国产成人精品v| 在线观看免费高清a一片| 永久网站在线| 欧美3d第一页| 成人国产av品久久久| 午夜亚洲福利在线播放| 九九在线视频观看精品| 亚洲精品成人av观看孕妇| 精品国产乱码久久久久久小说| 国产高清三级在线| 99九九线精品视频在线观看视频| 老司机影院毛片| 丝袜脚勾引网站| 免费观看a级毛片全部| 国产 一区 欧美 日韩| av一本久久久久| 天天一区二区日本电影三级| 国内精品宾馆在线| 一区二区三区精品91| 别揉我奶头 嗯啊视频| 九色成人免费人妻av| 在线精品无人区一区二区三 | 亚洲av成人精品一区久久| 亚洲最大成人av| h日本视频在线播放| 国产精品久久久久久精品电影小说 | 少妇人妻 视频| 大又大粗又爽又黄少妇毛片口| 舔av片在线| .国产精品久久| 国产精品99久久99久久久不卡 | 国产欧美日韩精品一区二区| 国产av国产精品国产| 寂寞人妻少妇视频99o| 欧美xxxx性猛交bbbb| 免费看a级黄色片| 亚洲精品中文字幕在线视频 | 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| 日韩 亚洲 欧美在线| 久久热精品热| 久久精品国产亚洲网站| 中文天堂在线官网| 亚洲精品影视一区二区三区av| 亚洲在久久综合| 黄色怎么调成土黄色| 日本熟妇午夜| av在线播放精品| 国产一区二区三区综合在线观看 | 搞女人的毛片| 精品少妇久久久久久888优播| 伊人久久精品亚洲午夜| 色视频在线一区二区三区| 精品少妇黑人巨大在线播放| 69人妻影院| 国产欧美亚洲国产| 久久韩国三级中文字幕| 久久女婷五月综合色啪小说 | 男女无遮挡免费网站观看| 精品久久久久久久久av| 午夜免费男女啪啪视频观看| 成人漫画全彩无遮挡| 尾随美女入室| 色哟哟·www| 看黄色毛片网站| 偷拍熟女少妇极品色| 国产高清不卡午夜福利| 久久99热6这里只有精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲综合色惰| 亚洲国产精品999| 午夜福利视频1000在线观看| 欧美精品人与动牲交sv欧美| 亚洲欧美日韩另类电影网站 | 黑人高潮一二区| 亚洲久久久久久中文字幕| 最新中文字幕久久久久| 亚洲精品成人av观看孕妇| 边亲边吃奶的免费视频| 人妻少妇偷人精品九色| 国产高清三级在线| 一级毛片黄色毛片免费观看视频| 欧美成人a在线观看| 大陆偷拍与自拍| 六月丁香七月| 九草在线视频观看| 色播亚洲综合网| 国产精品.久久久| 97热精品久久久久久| 国产老妇女一区| 久久午夜福利片| 日本-黄色视频高清免费观看| 成年人午夜在线观看视频| 久久热精品热| 色视频www国产| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 美女xxoo啪啪120秒动态图| 女人久久www免费人成看片| 天天躁日日操中文字幕| 精品熟女少妇av免费看| 国产黄色视频一区二区在线观看| 亚洲欧美精品专区久久| 成年av动漫网址| 各种免费的搞黄视频| 国产在线男女| 嫩草影院新地址| 欧美区成人在线视频| 日本一二三区视频观看| 国产免费福利视频在线观看| 男男h啪啪无遮挡| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99| 王馨瑶露胸无遮挡在线观看| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 夜夜看夜夜爽夜夜摸| 人妻少妇偷人精品九色| 婷婷色麻豆天堂久久| 欧美3d第一页| 国产精品秋霞免费鲁丝片| 亚洲在线观看片| 狂野欧美激情性xxxx在线观看| 人妻制服诱惑在线中文字幕| 亚洲va在线va天堂va国产| 久久热精品热| 极品教师在线视频| 99久久中文字幕三级久久日本| 80岁老熟妇乱子伦牲交| 免费观看性生交大片5| 国产av码专区亚洲av| 国产黄色视频一区二区在线观看| 亚洲精品日韩av片在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产欧美人成| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 蜜臀久久99精品久久宅男| 九九爱精品视频在线观看| 一级片'在线观看视频| 人人妻人人爽人人添夜夜欢视频 |