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

    水泵水輪機(jī)泵工況停機(jī)過程流動(dòng)特性分析

    2022-12-20 03:45:04姜?jiǎng)傥?/span>歐傳奇李騰飛楚士冀
    小水電 2022年6期
    關(guān)鍵詞:導(dǎo)葉轉(zhuǎn)輪水輪機(jī)

    姜?jiǎng)傥?,歐傳奇,趙 越,李騰飛,楚士冀

    (國際小水電中心, 浙江 杭州 310002)

    0 引 言

    目前,平穩(wěn)過渡綠色發(fā)展,穩(wěn)妥推動(dòng)可再生能源參與市場競爭成為能源利用的主題[1]。小水電作為清潔能源,是實(shí)現(xiàn)碳達(dá)峰、碳中和的重要選擇[2]。我國小水電資源豐富,全國小水電現(xiàn)有裝機(jī)容量約為世界各國小水電裝機(jī)容量的總和。但是大部分的小水電由于建設(shè)時(shí)間較為久遠(yuǎn),自身的安全問題日益凸顯,時(shí)下許多小水電都進(jìn)行了增容改造且對原有機(jī)組進(jìn)行完善。當(dāng)前,由于大中型抽蓄技術(shù)較為成熟,為將小水電的作用進(jìn)一步擴(kuò)大,國家鼓勵(lì)小水電機(jī)組向抽蓄機(jī)組改造。

    國內(nèi)對于大中型抽蓄的研究較為豐富。朱迪等[3]發(fā)現(xiàn)水泵水輪機(jī)在小開度大流量工況和大開度小流量工況下受許多低頻影響,導(dǎo)致壓力脈動(dòng)幅值較大,運(yùn)行穩(wěn)定性較差。游光華等[4]發(fā)現(xiàn)在泵工況停機(jī)過程中活動(dòng)導(dǎo)葉間的主流形態(tài)會發(fā)生改變。吳亞軍等[5]發(fā)現(xiàn)隨著導(dǎo)葉開度的增大,駝峰現(xiàn)象顯著,在駝峰區(qū)壓力脈動(dòng)較大,泵工況起動(dòng)過程中應(yīng)避免進(jìn)入駝峰區(qū)。李琪飛等[6]對不同開度下水泵水輪機(jī)泵工況的流場進(jìn)行研究,發(fā)現(xiàn)導(dǎo)葉開度越小,導(dǎo)葉水力矩和靜壓值的周期性越差,幅值越大。張春澤等[7]分析了水輪機(jī)內(nèi)的壓力脈動(dòng)和流態(tài)隨流量變化的特性,并得到了旋轉(zhuǎn)失速的演變規(guī)律。

    綜上所述,對于水泵水輪機(jī)過渡過程的研究對象主要是大中型機(jī)組,對于小容量機(jī)組的過渡過程研究很少。本文對小水電改造為抽蓄的小容量機(jī)組下的泵工況停機(jī)過程進(jìn)行研究,以期為小水電的穩(wěn)定運(yùn)行提供參考。

    1 數(shù)值計(jì)算方法

    1.1 基本方程

    水泵水輪機(jī)泵工況停機(jī)過程采用CFD軟件進(jìn)行模擬,其基于有限體積法求解N_S方程,N_S方程見公式(1)、(2)。在計(jì)算過程中,一般忽略水的密度變化所產(chǎn)生的影響[8]。

    連續(xù)性方程為:

    ·u=0

    (1)

    動(dòng)量方程為:

    (2)

    式中,u為流體速度(m/s);ρ為流體密度(kg/m3);p為靜壓(Pa);μ為動(dòng)力黏度系數(shù)(kg/(m·s) );為哈密頓算子,2為拉普拉斯算符,t為時(shí)間(s),g為重力加速度(m/s2)。

    1.2 邊界條件及初場設(shè)置

    湍流方程采用SST k_ω模型[9,10],SST k_ω模型主要應(yīng)用于內(nèi)部流動(dòng)、射流、大曲率流、分離流,對于泵工況停機(jī)過程的內(nèi)部流動(dòng)模擬效果較好,該方程的具體形式見公式(3)、(4)。蝸殼進(jìn)口及尾水管出口分別采用壓力進(jìn)口和壓力出口,數(shù)值根據(jù)進(jìn)出口的水頭差設(shè)置。對擴(kuò)散項(xiàng)和對流項(xiàng)采用二階迎風(fēng)格式進(jìn)行離散化處理,時(shí)間項(xiàng)采用一階迎風(fēng)格式進(jìn)行離散化處理,速度項(xiàng)和壓力項(xiàng)采用SIMPIC算法解耦。不同結(jié)構(gòu)之間的數(shù)據(jù)傳遞采用interface,參考壓力設(shè)置為1個(gè)大氣壓。時(shí)間步長設(shè)置為0.04 s,相當(dāng)于以轉(zhuǎn)輪的額定轉(zhuǎn)速旋轉(zhuǎn)6°的時(shí)間。

    為保證計(jì)算的穩(wěn)定性,在水泵水輪機(jī)作泵工況停機(jī)前,先進(jìn)行2.4 s(10個(gè)轉(zhuǎn)動(dòng)周期)的非定常計(jì)算設(shè)定初場。非定常計(jì)算中轉(zhuǎn)輪區(qū)域采用無滑移網(wǎng)格方法:

    (3)

    (4)

    lk-ω=k1/2βkω

    (5)

    式中,ρ為密度,P為生成項(xiàng),μ為動(dòng)力粘度,μt、μl為渦粘性系數(shù),σk、σω、σω2、βk為模型常數(shù),Cω為生成項(xiàng)系數(shù),F(xiàn)l為混合系數(shù),lk-ω為湍流尺度,k為湍流動(dòng)能強(qiáng)度,t為時(shí)間,ui為速度,xi為單位坐標(biāo),ω為湍動(dòng)能耗散項(xiàng)。

    1.3 算法建立

    水泵水輪機(jī)泵工況停機(jī)過程中,轉(zhuǎn)輪在t=9 s時(shí)與電動(dòng)機(jī)失去聯(lián)系,在此之前保持額定轉(zhuǎn)速不變。在t=9 s后轉(zhuǎn)輪的轉(zhuǎn)速變化根據(jù)力矩平衡方程得出:

    (6)

    其中,J表示轉(zhuǎn)輪的轉(zhuǎn)動(dòng)慣量,ω為轉(zhuǎn)輪的角速度。MH表示水力矩,Mf表示軸承摩擦等阻力矩,由于在水泵水輪機(jī)做泵工況過程中,摩擦產(chǎn)生的阻力矩較小,因此本次忽略阻力矩對轉(zhuǎn)輪轉(zhuǎn)動(dòng)的影響,轉(zhuǎn)速的變化快慢僅取決于水力矩的大小[11_12]。

    轉(zhuǎn)輪轉(zhuǎn)速和導(dǎo)葉開度角的變化均采用Fluent軟件的二次開發(fā)UDF實(shí)現(xiàn),本次泵工況停機(jī)的模擬總時(shí)長為25 s。轉(zhuǎn)輪轉(zhuǎn)速的變化情況見公式(7),泵工況停機(jī)過程中導(dǎo)葉開度角變化情況見公式(8),導(dǎo)葉開度角采用兩段關(guān)閉規(guī)律,第一段為線性關(guān)閉,第二段為曲線關(guān)閉。

    (7)

    式中,ωi為某一時(shí)刻的角速度,nr為額定轉(zhuǎn)速,ωi-1為上一時(shí)刻的角速度,Mi-1為上一時(shí)刻的力矩,J為機(jī)組的轉(zhuǎn)動(dòng)慣量,Δt為時(shí)間步長。

    (8)

    式中,Δγi為導(dǎo)葉單位時(shí)間的關(guān)閉角度,Δt為時(shí)間步長。

    2 計(jì)算模型及模擬可行性分析

    2.1 計(jì)算模型及網(wǎng)格劃分

    本次模擬的對象為某小水電站改造為抽蓄電站后的水泵水輪機(jī)(見圖1),包括蝸殼、活動(dòng)導(dǎo)葉和固定導(dǎo)葉、轉(zhuǎn)輪和尾水管。為保證計(jì)算的穩(wěn)定性,在蝸殼進(jìn)口和尾水管出口分別增加1段延伸段(見表1)。

    圖1 水泵水輪機(jī)示意圖

    表1 水泵水輪機(jī)基本參數(shù)

    網(wǎng)格劃分采用ICEM CFD軟件實(shí)現(xiàn)。根據(jù)各部分的復(fù)雜程度,對蝸殼采用非結(jié)構(gòu)網(wǎng)格進(jìn)行劃分,對導(dǎo)葉采用楔形網(wǎng)格進(jìn)行劃分,對轉(zhuǎn)輪和尾水管及延伸段采用結(jié)構(gòu)網(wǎng)格進(jìn)行劃分,并對靠近壁面處的網(wǎng)格進(jìn)行加密處理,使靠近壁面的網(wǎng)格y+平均值小于20(見圖2)。

    (a)網(wǎng)格劃分示意圖

    2.2 網(wǎng)格無關(guān)性驗(yàn)證

    網(wǎng)格數(shù)量的多少影響計(jì)算的準(zhǔn)確性和計(jì)算時(shí)長。采用6套不同網(wǎng)格數(shù)量的網(wǎng)格,利用水泵效率進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證(見圖2(b))。可以看出,當(dāng)網(wǎng)格數(shù)量超過543萬后,隨網(wǎng)格數(shù)量的增大,水泵水輪機(jī)泵工況的效率變化較??;因此模擬采用的網(wǎng)格數(shù)量最終選取為543萬。

    3 結(jié)果分析

    3.1 模擬結(jié)果對比

    為驗(yàn)證三維數(shù)值模擬的準(zhǔn)確性,對該水泵水輪機(jī)在導(dǎo)葉開度為18°下的模擬參數(shù)與試驗(yàn)數(shù)據(jù)進(jìn)行對比(見圖3)。可以看出,三維數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的誤差較小,證明了三維數(shù)值模結(jié)果的可信度。

    圖3 模擬結(jié)果與試驗(yàn)對比圖

    3.2 外特性變化情況分析

    圖4為水泵水輪機(jī)泵工況停機(jī)下的外特性變化曲線(見圖4),外特性參數(shù)包括轉(zhuǎn)速、流量、轉(zhuǎn)矩和導(dǎo)葉開度角;為便于定性分析,圖中曲線參數(shù)采用與初始時(shí)刻參數(shù)的百分比表示。圖中泵工況停機(jī)前的初始轉(zhuǎn)輪轉(zhuǎn)速為250 r/min,初始流量為44.964 m3/s,初始力矩為1.864×106N·m,初始導(dǎo)葉開度角為22.44°。

    圖4 水泵水輪機(jī)泵工況停機(jī)過程外特性變化

    初始時(shí)刻導(dǎo)葉收到關(guān)機(jī)命令發(fā)生動(dòng)作,導(dǎo)葉開度角逐漸變小。由于此時(shí)轉(zhuǎn)輪與電網(wǎng)相連,因此轉(zhuǎn)輪的轉(zhuǎn)速不發(fā)生變化,保持額定轉(zhuǎn)速不變,流量和轉(zhuǎn)矩隨導(dǎo)葉開度角變小的同時(shí)也逐漸變小。t=9 s時(shí)轉(zhuǎn)輪與電動(dòng)機(jī)失去聯(lián)系,轉(zhuǎn)輪轉(zhuǎn)速在水力矩的作用下逐漸變小,力矩變化幅度變大。在導(dǎo)葉進(jìn)入第二段曲線關(guān)閉后,力矩變化幅度變小,流量繼續(xù)減小;在t=22 s附近流量減小為0;在t=25 s時(shí),力矩繼續(xù)減小,轉(zhuǎn)輪轉(zhuǎn)速下降至初始轉(zhuǎn)速的65.2%。

    3.3 內(nèi)流特性分析

    水泵水輪機(jī)外特性的變化與內(nèi)部流態(tài)的變化有關(guān),圖5、圖6、圖7分別為轉(zhuǎn)輪及導(dǎo)葉區(qū)域流態(tài)和尾水管流態(tài)變化示意圖(見圖5~圖7)。

    圖5 轉(zhuǎn)輪及導(dǎo)葉區(qū)域流線圖和導(dǎo)葉位置壓力云圖

    初始時(shí)刻,轉(zhuǎn)輪區(qū)域及導(dǎo)葉區(qū)域的流動(dòng)情況較好,導(dǎo)葉位置兩側(cè)的壓力分布較為均勻,所受到的水力矩很小。隨著導(dǎo)葉開度的減小,在t=8 s時(shí)轉(zhuǎn)輪區(qū)域流線仍然保持良好的流動(dòng)狀態(tài),但出口水流導(dǎo)葉區(qū)域流線發(fā)生改變,這是由于此時(shí)導(dǎo)葉開度角較小,導(dǎo)致部分水流在流過導(dǎo)葉后發(fā)生脫流,使得在活動(dòng)導(dǎo)葉與固定導(dǎo)葉之間形成漩渦。導(dǎo)葉區(qū)域尾端的壓力變大,這是由于從轉(zhuǎn)輪流出的水流與導(dǎo)葉尾端發(fā)生碰撞,使得該區(qū)域的壓力得到一定的提高。t=15 s時(shí),轉(zhuǎn)輪葉道內(nèi)形成了明顯的渦結(jié)構(gòu),不同葉道內(nèi)渦結(jié)構(gòu)的數(shù)量不同,且流道內(nèi)流態(tài)不對稱,因此可以推斷出該工況的流動(dòng)穩(wěn)定性較差,轉(zhuǎn)輪受力不均。t=25 s時(shí),由于此時(shí)導(dǎo)葉完全關(guān)閉,導(dǎo)致由導(dǎo)葉間隙流出的水流在無葉區(qū)內(nèi)高速旋轉(zhuǎn),葉道內(nèi)渦結(jié)構(gòu)的分布更加混亂,且此時(shí)導(dǎo)葉所受到的水力矩最大。

    圖6 轉(zhuǎn)輪區(qū)域壓力云圖

    圖6為轉(zhuǎn)輪區(qū)域壓力云圖。在初始時(shí)刻,葉片壓力沿流動(dòng)方向均勻增大,在葉片出口邊位置壓力發(fā)生了一定的突增。當(dāng)導(dǎo)葉開度減小后,由于轉(zhuǎn)輪區(qū)域近乎于密閉空間,葉片流出的水流在無葉區(qū)受到嚴(yán)重的堵塞,使得無葉區(qū)的壓力產(chǎn)生一定的升高。而從導(dǎo)葉間隙流出的水流由于流速較大,導(dǎo)葉固定導(dǎo)葉和活動(dòng)導(dǎo)葉之間的區(qū)域壓力變小,在兩個(gè)活動(dòng)導(dǎo)葉之間水流的壓力最低,在流過活動(dòng)導(dǎo)葉后壓力慢慢得到恢復(fù),中心低壓區(qū)面積的大小與停機(jī)前近乎相同。t=15 s時(shí),導(dǎo)葉開度進(jìn)一步減小,此時(shí)葉片上的壓力分布不均勻,轉(zhuǎn)輪中心低壓區(qū)面積明顯增大。在t=25 s時(shí),轉(zhuǎn)輪區(qū)域壓力明顯減小,中心低壓區(qū)面積進(jìn)一步增大。

    圖7為尾水管區(qū)域流線圖。從圖中可以看出,在初始時(shí)刻水流從尾水管流入,在整個(gè)尾水管區(qū)域水流流態(tài)較好,能夠順暢地流入轉(zhuǎn)輪。由于在流進(jìn)方向上尾水管的截面面積逐漸減小,因此水流流速逐漸增大。t=8 s時(shí),在尾水管直錐段流態(tài)發(fā)生改變,部分水流緊貼壁面流向轉(zhuǎn)輪,在尾水管直錐段形成了渦流,而尾水管擴(kuò)散段的流態(tài)受到的影響較小。t=15 s時(shí),由于流量很小,直錐段渦流的影響效果逐漸擴(kuò)散至下游,在尾水管擴(kuò)散段水流流態(tài)異?;靵y,大部分水流緊貼壁面,形成渦流,僅有部分水流流向轉(zhuǎn)輪。由于流入到轉(zhuǎn)輪的水流減少,導(dǎo)致轉(zhuǎn)輪的低壓區(qū)域面積增大。t=25 s時(shí),轉(zhuǎn)輪的轉(zhuǎn)速逐漸減小,尾水管擴(kuò)散段渦流區(qū)域范圍減小,尾水管流態(tài)得到一定的改善。

    圖7 尾水管區(qū)域流線圖

    3.4 壓力脈動(dòng)分析

    水泵水輪機(jī)由于運(yùn)行過程中流態(tài)變化復(fù)雜,因此不同位置產(chǎn)生的壓力脈動(dòng)也各不相同。為分析水泵水輪機(jī)在泵工況停機(jī)過程中的壓力脈動(dòng)分布情況,在蝸殼位置、無葉區(qū)位置及尾水管位置處分別建立監(jiān)測點(diǎn)(見圖8)。

    圖8 監(jiān)測點(diǎn)位置選取

    圖9為不同位置監(jiān)測點(diǎn)在整個(gè)時(shí)間段內(nèi)的壓力變化情況(見圖9)。從圖中可以看出,在泵工況停機(jī)的初始階段,由于導(dǎo)葉關(guān)閉的影響,尾水管處的壓力無法得到釋放,形成了水錘效應(yīng),使得尾水管位置的壓力得到一定程度的升高;而在蝸殼位置和無葉區(qū)位置,壓力有一定程度的降低。t=9 s后,由于轉(zhuǎn)輪與電動(dòng)機(jī)失去聯(lián)系,轉(zhuǎn)輪轉(zhuǎn)速逐漸降低,水錘效應(yīng)程度逐漸變小,蝸殼位置的壓力得到恢復(fù),尾水管位置的壓力有所降低;而無葉區(qū)位置的壓力在升高一段時(shí)間后呈減低趨勢。這是由于這一階段的初始由于水錘效應(yīng)減弱,使得壓力得到一定的提升,后期由于轉(zhuǎn)輪轉(zhuǎn)速逐漸降低,導(dǎo)葉開度逐漸減小,使得無葉區(qū)的壓力逐漸降低。從圖中可以看出,在t=9 s時(shí),無葉區(qū)的壓力波動(dòng)程度很大,這是由于此處的壓力受到動(dòng)靜干涉的影響。在t=16.5 s后,蝸殼、無葉區(qū)和尾水管位置的壓力逐漸平穩(wěn)并在一定范圍內(nèi)波動(dòng)。

    圖9 泵工況停機(jī)過程中監(jiān)測點(diǎn)壓力變化

    4 結(jié) 論

    (1)水泵水輪機(jī)在泵工況停機(jī)過程中隨著導(dǎo)葉開度的減小,流量逐漸減小,力矩逐漸減小,轉(zhuǎn)速逐漸減小,導(dǎo)葉的關(guān)閉規(guī)律會影響外特性參數(shù)的變化速率。

    (2)在泵工況停機(jī)過程中,隨著導(dǎo)葉開度角的減小,轉(zhuǎn)輪內(nèi)部流態(tài)逐漸變差。開始時(shí),僅在導(dǎo)葉區(qū)域有漩渦產(chǎn)生,隨著轉(zhuǎn)輪轉(zhuǎn)速的降低,在轉(zhuǎn)輪區(qū)域也產(chǎn)生了不同情況的漩渦;漩渦的增多使得流量下降。隨著導(dǎo)葉的關(guān)閉,活動(dòng)導(dǎo)葉兩側(cè)的壓力差逐漸增大,活動(dòng)導(dǎo)葉所受到的水力矩逐漸增大。

    (3)初始時(shí)刻葉片壓力分布均勻,隨后由于導(dǎo)葉開度的減小,無葉區(qū)壓力變大。隨著停機(jī)過程的進(jìn)行,中心低壓區(qū)域范圍逐漸變大。

    (4)隨著導(dǎo)葉開度角的減小,首先在尾水管直錐段水流流態(tài)發(fā)生變化,形成渦流。隨著停機(jī)過程的進(jìn)行,渦流向下游發(fā)展,當(dāng)轉(zhuǎn)速逐漸變小后,尾水管流態(tài)得到一定的改善。

    (5)蝸殼位置、無葉區(qū)位置和尾水管位置的壓力變化規(guī)律在一定程度上受水錘的影響,無葉區(qū)壓力脈動(dòng)由于受動(dòng)靜干涉的影響,變化區(qū)間范圍較大。

    猜你喜歡
    導(dǎo)葉轉(zhuǎn)輪水輪機(jī)
    水輪機(jī)過流面非金屬材料的修復(fù)及防護(hù)
    大中型水斗式水輪機(jī)的關(guān)鍵技術(shù)
    詞語大轉(zhuǎn)輪
    ——“AABC”和“無X無X”式詞語
    寺廟里有座大書架——神奇的轉(zhuǎn)輪藏
    水輪機(jī)虛擬仿真動(dòng)畫制作的研究
    壓氣機(jī)導(dǎo)葉用GTD222合金鑄造技術(shù)研究
    模具制造(2019年4期)2019-06-24 03:36:50
    我國第一臺分半鑄造的30萬千瓦水輪發(fā)電機(jī)轉(zhuǎn)輪實(shí)制成功
    蘭臺世界(2017年12期)2017-06-22 14:17:25
    水輪機(jī)過流部件改造與節(jié)能增效
    混流式模型水輪機(jī)導(dǎo)葉水力矩試驗(yàn)
    一種新型的電子式導(dǎo)葉主令控制器的設(shè)計(jì)及應(yīng)用
    可以免费在线观看a视频的电影网站| 亚洲国产欧洲综合997久久, | 日韩欧美一区视频在线观看| 高清在线国产一区| 日本 欧美在线| 波多野结衣高清无吗| 熟女少妇亚洲综合色aaa.| 成在线人永久免费视频| 日本免费一区二区三区高清不卡| 亚洲 欧美一区二区三区| 18禁裸乳无遮挡免费网站照片 | 午夜福利在线在线| 黄色丝袜av网址大全| 久久久精品欧美日韩精品| 欧美+亚洲+日韩+国产| 中文字幕最新亚洲高清| 国产一卡二卡三卡精品| 欧美三级亚洲精品| 亚洲人成网站高清观看| 最近最新中文字幕大全免费视频| 久久人人精品亚洲av| 亚洲中文字幕日韩| 精品国内亚洲2022精品成人| 18禁黄网站禁片午夜丰满| 嫩草影视91久久| 人人妻人人看人人澡| 亚洲第一青青草原| www.精华液| 男人舔女人的私密视频| 观看免费一级毛片| 国产av不卡久久| 国产亚洲精品久久久久久毛片| 国产在线精品亚洲第一网站| 亚洲熟妇中文字幕五十中出| 亚洲精品av麻豆狂野| 两个人看的免费小视频| 高清毛片免费观看视频网站| 亚洲自偷自拍图片 自拍| 精品久久久久久久毛片微露脸| 美女免费视频网站| 亚洲国产高清在线一区二区三 | 欧美成人免费av一区二区三区| 1024手机看黄色片| 淫秽高清视频在线观看| 宅男免费午夜| 欧美黑人精品巨大| 2021天堂中文幕一二区在线观 | 哪里可以看免费的av片| 人人妻,人人澡人人爽秒播| 亚洲国产精品久久男人天堂| 日本一区二区免费在线视频| 变态另类丝袜制服| www.自偷自拍.com| 白带黄色成豆腐渣| 中文字幕最新亚洲高清| tocl精华| 亚洲精品粉嫩美女一区| 男人舔奶头视频| 欧美精品啪啪一区二区三区| 久久久久久久久免费视频了| 亚洲欧美激情综合另类| 中亚洲国语对白在线视频| 国产亚洲精品av在线| 亚洲第一电影网av| 夜夜夜夜夜久久久久| 香蕉丝袜av| 午夜福利在线在线| 精品不卡国产一区二区三区| 又黄又粗又硬又大视频| 亚洲国产精品久久男人天堂| 这个男人来自地球电影免费观看| 视频区欧美日本亚洲| 国产免费av片在线观看野外av| 嫩草影视91久久| 欧美不卡视频在线免费观看 | 成人av一区二区三区在线看| 国产免费男女视频| 在线av久久热| 久久久精品欧美日韩精品| 黄色视频不卡| 国产免费av片在线观看野外av| 亚洲成人免费电影在线观看| netflix在线观看网站| 欧美成人免费av一区二区三区| 真人一进一出gif抽搐免费| 中文字幕人妻丝袜一区二区| 久久精品人妻少妇| 亚洲 欧美一区二区三区| 久久伊人香网站| 天堂影院成人在线观看| 在线观看免费视频日本深夜| 久久 成人 亚洲| 成人三级做爰电影| 久久中文字幕一级| av欧美777| 真人做人爱边吃奶动态| 成熟少妇高潮喷水视频| 一级作爱视频免费观看| 性色av乱码一区二区三区2| 在线永久观看黄色视频| 好男人电影高清在线观看| 中文字幕精品免费在线观看视频| 午夜成年电影在线免费观看| 国产av又大| 欧美性猛交╳xxx乱大交人| www日本黄色视频网| 身体一侧抽搐| 99国产精品一区二区三区| 亚洲人成网站高清观看| 欧美乱妇无乱码| 亚洲欧美一区二区三区黑人| 欧美黄色片欧美黄色片| 国产又色又爽无遮挡免费看| 午夜视频精品福利| 天天一区二区日本电影三级| 国产成人av激情在线播放| 国产单亲对白刺激| 中文字幕人妻熟女乱码| 成人午夜高清在线视频 | 精品免费久久久久久久清纯| 欧美日本亚洲视频在线播放| 亚洲 欧美 日韩 在线 免费| 欧美乱色亚洲激情| 久久热在线av| 免费在线观看黄色视频的| 国产精品av久久久久免费| 久久精品国产亚洲av香蕉五月| 日韩欧美国产一区二区入口| 亚洲欧美精品综合一区二区三区| 99精品在免费线老司机午夜| 国产单亲对白刺激| 国产极品粉嫩免费观看在线| 中文字幕另类日韩欧美亚洲嫩草| 国产精品,欧美在线| 美女大奶头视频| 波多野结衣高清无吗| 不卡一级毛片| 成人国产综合亚洲| 91在线观看av| 亚洲人成网站高清观看| 亚洲欧美精品综合久久99| 欧美色欧美亚洲另类二区| 久久香蕉激情| 一个人观看的视频www高清免费观看 | 最近最新免费中文字幕在线| 欧美成人性av电影在线观看| 高清毛片免费观看视频网站| 久久中文字幕一级| 波多野结衣高清作品| 久久精品影院6| 岛国视频午夜一区免费看| 亚洲片人在线观看| 人人妻,人人澡人人爽秒播| 黄色 视频免费看| 日韩中文字幕欧美一区二区| 欧美乱码精品一区二区三区| 美国免费a级毛片| 制服人妻中文乱码| 色av中文字幕| 亚洲第一欧美日韩一区二区三区| 婷婷丁香在线五月| 免费在线观看黄色视频的| 1024手机看黄色片| 成人午夜高清在线视频 | 久久精品成人免费网站| 国产精品av久久久久免费| 亚洲九九香蕉| 黄频高清免费视频| 琪琪午夜伦伦电影理论片6080| 欧美黄色片欧美黄色片| 日韩欧美一区二区三区在线观看| 成人三级黄色视频| 自线自在国产av| 亚洲精品一卡2卡三卡4卡5卡| 日本熟妇午夜| 精品久久久久久久人妻蜜臀av| 国内久久婷婷六月综合欲色啪| 人妻久久中文字幕网| 久久久国产成人免费| 色播亚洲综合网| 免费高清在线观看日韩| 国产私拍福利视频在线观看| 手机成人av网站| 精品日产1卡2卡| 中文字幕另类日韩欧美亚洲嫩草| 无遮挡黄片免费观看| 人人澡人人妻人| 欧美乱妇无乱码| 天堂动漫精品| 在线观看免费视频日本深夜| www国产在线视频色| 天天添夜夜摸| 成人国产一区最新在线观看| 男男h啪啪无遮挡| 亚洲avbb在线观看| 大型av网站在线播放| 伊人久久大香线蕉亚洲五| 日本一区二区免费在线视频| 男女下面进入的视频免费午夜 | 免费人成视频x8x8入口观看| 亚洲精品久久成人aⅴ小说| 国产男靠女视频免费网站| 午夜福利18| 日韩欧美三级三区| 国产精品免费一区二区三区在线| 波多野结衣av一区二区av| 人人妻,人人澡人人爽秒播| 好男人电影高清在线观看| 亚洲专区中文字幕在线| 午夜福利高清视频| 欧美日韩黄片免| 亚洲五月色婷婷综合| 看黄色毛片网站| 啦啦啦韩国在线观看视频| 日韩欧美在线二视频| 亚洲avbb在线观看| 成熟少妇高潮喷水视频| 老司机午夜福利在线观看视频| 两人在一起打扑克的视频| 久久久国产成人免费| 午夜激情福利司机影院| 久9热在线精品视频| 久99久视频精品免费| 国产精品美女特级片免费视频播放器 | 一级黄色大片毛片| 天天一区二区日本电影三级| 不卡av一区二区三区| 国产精品乱码一区二三区的特点| 搡老熟女国产l中国老女人| 国内精品久久久久精免费| 丁香六月欧美| 精品第一国产精品| 亚洲人成网站高清观看| 久久国产精品人妻蜜桃| 我的亚洲天堂| svipshipincom国产片| svipshipincom国产片| 亚洲欧美日韩高清在线视频| 久久伊人香网站| 久久伊人香网站| 国产高清videossex| 欧美中文综合在线视频| 丰满人妻熟妇乱又伦精品不卡| 国产伦一二天堂av在线观看| 久久精品国产综合久久久| 看片在线看免费视频| 老熟妇乱子伦视频在线观看| 国产亚洲av高清不卡| 黄色 视频免费看| 国内少妇人妻偷人精品xxx网站 | 搡老岳熟女国产| 婷婷精品国产亚洲av在线| 韩国精品一区二区三区| www.自偷自拍.com| 久热这里只有精品99| 淫秽高清视频在线观看| 亚洲av熟女| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看成人毛片| 国产精品久久电影中文字幕| 窝窝影院91人妻| 成人精品一区二区免费| 黄色丝袜av网址大全| 国产极品粉嫩免费观看在线| 欧美性猛交黑人性爽| 亚洲国产看品久久| 精品国产亚洲在线| 久久久久久久久中文| 首页视频小说图片口味搜索| 免费人成视频x8x8入口观看| 久久中文看片网| 欧美成人午夜精品| 窝窝影院91人妻| 国产色视频综合| 亚洲人成77777在线视频| 日韩精品免费视频一区二区三区| 成人特级黄色片久久久久久久| 亚洲国产日韩欧美精品在线观看 | 听说在线观看完整版免费高清| 精品久久久久久久毛片微露脸| 亚洲精品中文字幕一二三四区| 久久天躁狠狠躁夜夜2o2o| 国产真实乱freesex| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品久久男人天堂| 麻豆一二三区av精品| 欧美日韩精品网址| www.自偷自拍.com| 国产精品久久久人人做人人爽| 精品久久久久久久毛片微露脸| 中国美女看黄片| 国产色视频综合| xxxwww97欧美| 91大片在线观看| 成人欧美大片| 欧美乱色亚洲激情| 欧美日韩黄片免| 亚洲精品一区av在线观看| 午夜免费观看网址| 亚洲熟妇熟女久久| 国产黄色小视频在线观看| 国产av又大| 成年免费大片在线观看| 黄片播放在线免费| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩精品网址| 亚洲美女黄片视频| 他把我摸到了高潮在线观看| av在线天堂中文字幕| 久久天躁狠狠躁夜夜2o2o| 国产av又大| 中文在线观看免费www的网站 | 美女高潮喷水抽搐中文字幕| 日日干狠狠操夜夜爽| 亚洲国产欧美一区二区综合| 黄色片一级片一级黄色片| 可以免费在线观看a视频的电影网站| 久久草成人影院| 黄片播放在线免费| 精品国产国语对白av| 亚洲五月色婷婷综合| 午夜免费成人在线视频| 男女做爰动态图高潮gif福利片| 国产精品日韩av在线免费观看| 日韩大码丰满熟妇| 国内精品久久久久久久电影| 午夜久久久久精精品| 真人做人爱边吃奶动态| 黄色成人免费大全| 国产蜜桃级精品一区二区三区| 国产一区二区三区在线臀色熟女| 久久久国产精品麻豆| 国产精品 国内视频| 美女 人体艺术 gogo| √禁漫天堂资源中文www| 国产伦一二天堂av在线观看| 亚洲五月色婷婷综合| 成人三级黄色视频| 国产一区二区在线av高清观看| 99精品久久久久人妻精品| 国产区一区二久久| 国产成人av激情在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 99在线人妻在线中文字幕| 夜夜爽天天搞| 亚洲第一欧美日韩一区二区三区| 热re99久久国产66热| 亚洲成人免费电影在线观看| 日韩欧美一区二区三区在线观看| 久久精品国产亚洲av高清一级| 日韩一卡2卡3卡4卡2021年| 国产高清videossex| 亚洲一区二区三区不卡视频| 国产午夜福利久久久久久| 亚洲第一av免费看| 国产精品久久久久久精品电影 | 免费看十八禁软件| 欧美日韩一级在线毛片| 久久伊人香网站| 一级毛片精品| 亚洲午夜理论影院| 亚洲人成网站在线播放欧美日韩| 最好的美女福利视频网| 亚洲,欧美精品.| av视频在线观看入口| 精品午夜福利视频在线观看一区| 热re99久久国产66热| 亚洲成人免费电影在线观看| 国产成人欧美| 亚洲中文字幕日韩| 成人国产一区最新在线观看| 成人av一区二区三区在线看| 99精品欧美一区二区三区四区| 嫩草影院精品99| 好男人在线观看高清免费视频 | 中文在线观看免费www的网站 | 最近最新免费中文字幕在线| 日韩有码中文字幕| 美女扒开内裤让男人捅视频| 9191精品国产免费久久| 最好的美女福利视频网| 精品久久久久久久毛片微露脸| 99久久精品国产亚洲精品| 午夜免费鲁丝| 久久久久久免费高清国产稀缺| 久久婷婷人人爽人人干人人爱| 一个人免费在线观看的高清视频| 欧美成狂野欧美在线观看| 最新美女视频免费是黄的| 欧美日韩瑟瑟在线播放| 日本免费一区二区三区高清不卡| 一区二区日韩欧美中文字幕| 精品国产乱码久久久久久男人| 18禁裸乳无遮挡免费网站照片 | 国产av又大| 久久午夜综合久久蜜桃| 精品国产亚洲在线| 久久久久久大精品| 国产在线精品亚洲第一网站| 亚洲av成人不卡在线观看播放网| 午夜福利18| 日韩欧美国产在线观看| 亚洲一区二区三区色噜噜| 亚洲精品国产精品久久久不卡| 精品久久蜜臀av无| 老汉色∧v一级毛片| 91av网站免费观看| 日本三级黄在线观看| 国产99白浆流出| 一个人免费在线观看的高清视频| 在线av久久热| 国产精品99久久99久久久不卡| 欧美黄色片欧美黄色片| av福利片在线| 中文字幕高清在线视频| 亚洲国产精品sss在线观看| 制服丝袜大香蕉在线| 精品久久久久久久久久久久久 | 黄色 视频免费看| 日韩欧美免费精品| 久久久久国产一级毛片高清牌| 熟女电影av网| 精品久久久久久久毛片微露脸| 欧美国产日韩亚洲一区| 日韩免费av在线播放| 免费人成视频x8x8入口观看| 高清在线国产一区| 精品国产超薄肉色丝袜足j| 视频区欧美日本亚洲| 国产精品永久免费网站| 精品电影一区二区在线| 国产精品 国内视频| 免费一级毛片在线播放高清视频| 国产精品久久久久久亚洲av鲁大| 欧美乱妇无乱码| 手机成人av网站| 国产精品九九99| 99精品久久久久人妻精品| 757午夜福利合集在线观看| 99热只有精品国产| 国产成人一区二区三区免费视频网站| 色综合婷婷激情| 美女午夜性视频免费| 中国美女看黄片| 精品一区二区三区av网在线观看| 国产男靠女视频免费网站| 高清毛片免费观看视频网站| 在线观看舔阴道视频| 久久中文看片网| av片东京热男人的天堂| 国产av在哪里看| 久久久国产精品麻豆| 特大巨黑吊av在线直播 | АⅤ资源中文在线天堂| 国产精品亚洲av一区麻豆| 成人三级黄色视频| 国产黄色小视频在线观看| 亚洲 欧美一区二区三区| 国产主播在线观看一区二区| 麻豆成人午夜福利视频| 欧美精品啪啪一区二区三区| 亚洲国产欧洲综合997久久, | 欧美三级亚洲精品| 91av网站免费观看| 91老司机精品| 一a级毛片在线观看| 别揉我奶头~嗯~啊~动态视频| 老司机在亚洲福利影院| 欧美激情极品国产一区二区三区| 色av中文字幕| 成人av一区二区三区在线看| 亚洲av第一区精品v没综合| 人人澡人人妻人| 国产成人系列免费观看| 免费在线观看黄色视频的| 日本成人三级电影网站| 男女做爰动态图高潮gif福利片| aaaaa片日本免费| 动漫黄色视频在线观看| 亚洲一码二码三码区别大吗| 亚洲国产欧美网| 久久婷婷人人爽人人干人人爱| 99久久国产精品久久久| 一级黄色大片毛片| 制服丝袜大香蕉在线| 国内精品久久久久精免费| 18禁黄网站禁片免费观看直播| 中文字幕人妻熟女乱码| 午夜激情av网站| 精品国产乱码久久久久久男人| 黑人欧美特级aaaaaa片| 啦啦啦免费观看视频1| xxxwww97欧美| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 99国产精品99久久久久| 91麻豆精品激情在线观看国产| 欧美丝袜亚洲另类 | 中文字幕高清在线视频| 免费av毛片视频| 久久久国产成人精品二区| 午夜福利高清视频| 这个男人来自地球电影免费观看| 国产亚洲精品综合一区在线观看 | www日本黄色视频网| 国产真实乱freesex| 久久香蕉精品热| 最近最新中文字幕大全免费视频| 久久精品国产综合久久久| 18禁国产床啪视频网站| 少妇的丰满在线观看| 亚洲成人久久性| 最近最新中文字幕大全免费视频| 久久青草综合色| 草草在线视频免费看| 亚洲成人久久爱视频| 国产一级毛片七仙女欲春2 | 国产成人av激情在线播放| 村上凉子中文字幕在线| 国产国语露脸激情在线看| 黄网站色视频无遮挡免费观看| 叶爱在线成人免费视频播放| 国产99久久九九免费精品| 在线永久观看黄色视频| 精品午夜福利视频在线观看一区| 日韩欧美国产一区二区入口| 99在线视频只有这里精品首页| 国产黄片美女视频| 国产精品av久久久久免费| 淫妇啪啪啪对白视频| 国产免费av片在线观看野外av| 首页视频小说图片口味搜索| 亚洲 国产 在线| 男人操女人黄网站| 黑人巨大精品欧美一区二区mp4| 亚洲欧美一区二区三区黑人| 精品久久久久久久人妻蜜臀av| 亚洲精品色激情综合| 国产精品永久免费网站| 18禁黄网站禁片免费观看直播| 非洲黑人性xxxx精品又粗又长| 欧美亚洲日本最大视频资源| 免费看十八禁软件| 美国免费a级毛片| 成人免费观看视频高清| 国内精品久久久久精免费| 在线永久观看黄色视频| 99热这里只有精品一区 | 国产免费av片在线观看野外av| 99热6这里只有精品| 久久久久久人人人人人| 在线观看免费日韩欧美大片| 又黄又爽又免费观看的视频| 嫩草影院精品99| 中文字幕av电影在线播放| 18禁国产床啪视频网站| 伦理电影免费视频| 久热这里只有精品99| 丰满人妻熟妇乱又伦精品不卡| 国产又爽黄色视频| 国产97色在线日韩免费| 久久久久亚洲av毛片大全| 可以在线观看的亚洲视频| 国产99白浆流出| 午夜福利视频1000在线观看| 国产乱人伦免费视频| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲精品综合一区在线观看 | 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 精品国产美女av久久久久小说| 三级毛片av免费| 人人妻,人人澡人人爽秒播| 亚洲精华国产精华精| 啦啦啦韩国在线观看视频| 久久99热这里只有精品18| 欧美色视频一区免费| 一本大道久久a久久精品| 18禁国产床啪视频网站| 高清在线国产一区| 黄色女人牲交| 国产激情久久老熟女| 97人妻精品一区二区三区麻豆 | 可以免费在线观看a视频的电影网站| 97碰自拍视频| 99在线视频只有这里精品首页| 欧美日韩黄片免| 美女高潮到喷水免费观看| 91大片在线观看| 中文字幕人妻丝袜一区二区| 我的亚洲天堂| 露出奶头的视频| 岛国在线观看网站| 俺也久久电影网| 每晚都被弄得嗷嗷叫到高潮| xxxwww97欧美| 一区福利在线观看| 两个人看的免费小视频| 国产精品日韩av在线免费观看| 亚洲 欧美一区二区三区| 国产91精品成人一区二区三区| 午夜久久久久精精品| 亚洲精品色激情综合| 最新美女视频免费是黄的| 岛国视频午夜一区免费看| 夜夜躁狠狠躁天天躁| 女性生殖器流出的白浆| 成人精品一区二区免费| 天天添夜夜摸| 日韩国内少妇激情av| 满18在线观看网站| 老司机在亚洲福利影院| 十八禁人妻一区二区| 悠悠久久av| 免费在线观看成人毛片|