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

    小流量工況泵閥耦合系統(tǒng)失穩(wěn)機理研究

    2023-02-27 13:14:32率志君李玩幽姜晨醒
    振動與沖擊 2023年4期

    于 濤,率志君,王 曦,簡 潔,李玩幽,姜晨醒,肖 頎

    (1. 哈爾濱工程大學 動力與能源工程學院,哈爾濱 150000;2. 廈門大學 能源學院,福建 廈門 361000;3. 武漢第二船舶設計研究所 熱能動力技術重點實驗室,武漢 430205)

    泵閥耦合系統(tǒng)的穩(wěn)定性對于核電站給水系統(tǒng)的安全運行至關重要。離心泵、彈簧閥以及管路等部件是核電站給水系統(tǒng)的重要組成部分,各元件協(xié)同作用,保證系統(tǒng)正常運行。彈簧止回閥安裝在離心泵出口處,閥門開度根據(jù)不同工況下的流量和壓力等參數(shù)進行自我調節(jié),使得泵閥耦合系統(tǒng)能夠適應多種運行工況,并表現(xiàn)出良好的特性[1]。核電站系統(tǒng)的安全設計至關重要,其中穩(wěn)定可靠的回路系統(tǒng)是保證核電系統(tǒng)平穩(wěn)運行的重要組成部分,因此對低振動的泵閥耦合子系統(tǒng)提出了更高要求。

    根據(jù)系統(tǒng)運行需要,泵閥耦合系統(tǒng)往往需要在小流量工況條件下運行,容易引起系統(tǒng)的流動失穩(wěn)。一方面,離心泵流場內(nèi)部會發(fā)生失速現(xiàn)象,從而誘發(fā)幅值較高的低頻壓力脈動,使得離心泵振動加劇[2-9]。周佩劍等[10]研究了小流量工況下離心泵隔舌區(qū)域的流場特性,發(fā)現(xiàn)失速團常出現(xiàn)在逆壓梯度較大的區(qū)域,造成流道堵塞,引發(fā)高幅值壓力脈動。瞿麗霞等[11]研究了雙吸離心泵的流場特性,發(fā)現(xiàn)偏離額定流量程度越大,葉輪內(nèi)部的壓力脈動幅值越高。Sano等[12]采用油膜法顯示技術捕捉到了離心泵導葉中的失速現(xiàn)象,并研究了導葉失速所引起的低頻壓力脈動。另一方面,當離心泵下游使用受系統(tǒng)參數(shù)實時調節(jié)的被動式彈簧閥時,閥內(nèi)流動特性和閥瓣瞬態(tài)動力學特性受上游壓力波動影響很大,可能導致閥瓣顫振,甚至出現(xiàn)敲擊現(xiàn)象,嚴重影響泵閥耦合系統(tǒng)的可靠性[13-17]。張偉濤[18]和劉華坪等[19]研究了閥門在開啟和關閉過程中的閥瓣受力特性以及逆流沖擊的產(chǎn)生原因。Botros[20]研究了止回閥的動態(tài)過程,發(fā)現(xiàn)當水動力小于彈簧力時,閥門不能全開,閥體振動加劇。Qian等[21]研究了不同進口壓力下閥瓣運動規(guī)律以及閥瓣受力情況,發(fā)現(xiàn)進口波動壓力對閥瓣受力影響很大,并且研究了閥內(nèi)彈簧剛度對開閥過程的影響。

    目前對于泵閥耦合系統(tǒng)流致激勵特性的研究,更多關注于固定閥門開度或閥門開度變化規(guī)律已知的情況[22-26]。張玉良等[27-28]建立了離心泵的循環(huán)管路系統(tǒng),研究了啟動過程中的外特性和內(nèi)部流動演化特性。楊帥[29]采用MOC-CFD耦合模擬,研究了快速開閥和關閥過程中泵的動態(tài)特性和瞬態(tài)流動的演化過程。劉巧玲[30]采用一維與三維耦合計算,研究了不同啟動方式下并聯(lián)泵送系統(tǒng)瞬態(tài)特性。對于使用受離心泵出口壓力實時調節(jié)的被動式彈簧閥的研究仍是較少。因此研究小流量工況泵閥耦合系統(tǒng)的流場特性以及結構動力學特性具有重要的科研價值。本文以某泵閥耦合系統(tǒng)為對象,重點研究小流量工況非穩(wěn)態(tài)流致激勵特性,并闡明泵閥耦合系統(tǒng)顫振現(xiàn)象的產(chǎn)生機理。

    1 計算區(qū)域及方法

    1.1 彈簧止回閥工作原理

    本文所研究的彈簧止回閥是由一個進口和兩個出口所組成的三通結構。彈簧止回閥的閥門開度受到流體介質的壓力,閥瓣的重力和彈簧的阻力等共同作用。當入口壓力小于閥瓣的觸發(fā)壓力時,流體通過閥桿底部進入閥桿內(nèi)部,隨后經(jīng)由空排口流出閥門;當入口壓力達到閥瓣觸發(fā)壓力,閥門完全打開,閥桿隨著閥瓣向上運動,流體介質無法通過閥桿流向空排口,而只能由主出口排出。由圖1閥瓣受力分析可知,閥瓣運動過程是在壓差引起的流體力、運動部件重力和彈簧力的平衡過程中實現(xiàn)的,根據(jù)牛頓第二定律總結出閥瓣運動方程如式(1)所示。

    (1)

    式中:m為運動部件質量;k為彈簧剛度;F為閥瓣所受流體力;c為阻尼;x0為彈簧的初始位移。

    圖1 閥瓣受力分析Fig.1 Valve force analysis

    1.2 計算區(qū)域

    本文的研究的泵閥耦合系統(tǒng)是一個復雜的多部件耦合系統(tǒng),離心泵和彈簧止回閥是系統(tǒng)中兩個主要的運動設備。泵出口流體介質推動彈簧止回閥閥瓣運動,是閥瓣所受流體激勵的主要來源,離心泵出口壓力變化和流量波動情況將直接影響閥門內(nèi)部流場特性。泵閥耦合系統(tǒng)的幾何和性能參數(shù)由表1給出。針對離心泵計算,旋轉區(qū)域葉輪要單獨進行設置,因而將離心泵的流體域分成四部分:入口導葉、葉輪、出口導葉以及蝸殼。對于彈簧止回閥計算,需要將閥瓣壁面和閥桿設置為剛體運動,并賦予一定的運動規(guī)律,進而模擬閥瓣的運動過程,因此將彈簧止回閥的流體域分成三部分:閥體、閥桿以及空排口。離心泵和彈簧止回閥的流體區(qū)域如圖2所示。

    表1 主要參數(shù)Tab.1 Main parameters

    圖2 主要部件流體模型Fig.2 Fluid model of main components

    試驗中系統(tǒng)通過主出口的閘閥控制回路流量,仿真中將其簡化為具有縮擴結構的管路,無需施加流量邊界條件而自耦合求解,流量參數(shù)變化通過內(nèi)部計算獲得。通過減小閥門喉部直徑來增大局部水力損失,進而減弱過流能力,通過調節(jié)閥門喉部直徑即可獲得相應的穩(wěn)定流量和系統(tǒng)阻力。通過試驗數(shù)據(jù)分析發(fā)現(xiàn),當總流量小于0.14倍額定流量時,閥門處于關閉狀態(tài),流體只能通過空排口排出;當總流量大于0.4倍額定流量時,閥門完全打開,流體只能通過主出口排出;當泵閥耦合系統(tǒng)在0.14~0.4倍額定流量之內(nèi)運行時,閥門處于半開半閉狀態(tài),流體既能通過主出口排出又能通過空排口排出,且開度越小,閥瓣處的流動越復雜,受離心泵出口壓力脈動的影響較大。根據(jù)系統(tǒng)運行需要,0.2倍額定流量為典型工況點,更容易引起系統(tǒng)流動失穩(wěn),所以本文選取該流量點進行研究。為了更真實的模擬流場內(nèi)部流動,對泵閥耦合系統(tǒng)進出口流域做相應延長,以便其在邊界處的流動能夠充分發(fā)展。本文采用四面體單元對流體區(qū)域進行網(wǎng)格劃分,計算區(qū)域網(wǎng)格如圖3所示。

    圖3 計算區(qū)域網(wǎng)格Fig.3 Computational domain mesh

    1.3 計算方法

    本文使用ANSYS FLUENT軟件對計算域的控制方程進行求解,并采用RNGk~ε湍流模型封閉方程組[31],其主要特點是考慮了高波數(shù)、小尺度渦對低波數(shù)、大尺度渦的影響,提高了漩渦流動的精度,在求解離心泵和閥門流場方面具有優(yōu)勢[32-35]。進口采取總壓邊界條件,出口采取靜壓邊界條件,均設置為0.1 MPa。葉輪流域設置為旋轉域,轉速為2 980 r/min,同時導葉流域與葉輪前后蓋板外側相接觸的部分,設置為旋轉壁面,其他流體域設置為靜止區(qū)域。彈簧止回閥的閥瓣運動過程受壓差大小以及彈簧位移量實時控制,因此采用UDF宏命令對閥瓣運動方程進行編譯,其具體計算流程如圖4所示。根據(jù)牛頓第二定律求取閥瓣的運動加速度和速度,并將運動速度賦值給閥瓣,進而指定剛體的運動方式。網(wǎng)格更新方法采用Smoothing和Remeshing結合的方式,對網(wǎng)格質量不好的區(qū)域進行網(wǎng)格重新劃分。

    圖4 計算流程圖Fig.4 Calculation process diagram

    在瞬態(tài)求解中,采用滑移網(wǎng)格模型來模擬動靜干涉效應,葉輪每旋轉2度計算一個時間步。在流場計算分析之前,分別對離心泵和彈簧止回閥進行網(wǎng)格無關性以及Yplus值檢驗。由圖5可知,當離心泵和彈簧止回閥的網(wǎng)格數(shù)量分別達到160萬和190萬之后,網(wǎng)格數(shù)目繼續(xù)增長對揚程和閥瓣升力的收斂性貢獻很小。因此綜合考慮計算精度和時間成本,本文采取350萬網(wǎng)格總量的計算模型,各區(qū)域網(wǎng)格數(shù)為:閥體132.1萬、空排口26.7萬、閥桿15.8萬、入口導葉14.9萬、葉輪36.4萬、蝸殼導葉99.5萬和進出口管路24.6萬。所有流域的Yplus值均小于150,滿足RNGk~ε湍流模型的適用條件。

    圖5 模型檢驗Fig.5 Model checking

    1.4 水力特性試驗驗證

    本文試驗回路由離心泵、彈簧止回閥、閘閥和管道回路組成,試驗裝置臺如圖6所示。離心泵由電機驅動,變頻電機轉速為2 980 r/min,滿足試驗要求。

    圖7(a)和圖7(b)為離心泵和彈簧止回閥的外特性曲線,其中流量比為主出口流量和空排口流量的比值。通過對比試驗與仿真數(shù)據(jù)發(fā)現(xiàn),二者基本吻合。圖7(c)為離心泵出口壓力頻域曲線,對比發(fā)現(xiàn)離心泵出口壓力在葉頻248.3 Hz及其倍頻和低頻3.8 Hz處均存在明顯的峰值,仿真結果與試驗結果基本一致,僅存在幅值上的區(qū)別,主要原因是計算過程中對邊界條件的簡化。綜上所述,數(shù)值計算結果可信。

    注:1.水箱;2.壓力表;3.離心泵;4.壓力脈動傳感器;5.彈簧止回閥;6.閘閥;7.流量計。圖6 試驗裝置臺結構圖Fig.6 Test system structure

    圖7 試驗值和仿真值對比Fig.7 Comparison of experimental and simulation results

    2 計算結果及分析

    2.1 內(nèi)部流動分析

    由圖8泵閥耦合系統(tǒng)壓力和速度分布圖可知,由于離心泵對經(jīng)其內(nèi)部流體的增壓作用,離心泵流場內(nèi)部壓力沿流動方向壓力逐漸增加,并且受動靜干涉作用的影響,導葉隔舌處是整個流場內(nèi)壓力變化最劇烈的位置,同時葉輪內(nèi)壓力等值線與葉片表面不再垂直,說明葉片表面處存在著較大壓力梯度。當流體流經(jīng)閘閥和閥桿進出口等過流面積較小的區(qū)域時,受節(jié)流作用的影響,主出口和空排出口處壓力明顯降低。流體在經(jīng)過葉輪做功之后,流速明顯升高,在導葉入口處達到最大,并且葉輪和導葉流道存在明顯的流動分離、漩渦和回流等不穩(wěn)定現(xiàn)象。當流體流經(jīng)狹小的閥瓣間隙和閥桿窗口時,由于過流面積較小,流體被擠壓為一股高速射流,這樣多股水流相遇對撞以及流動死角的存在所引起的銳緣現(xiàn)象,使得閥門流場內(nèi)部出現(xiàn)了尺度各異的漩渦結構。

    2.2 瞬態(tài)流動特性分析

    小流量工況下泵閥耦合系統(tǒng)流場內(nèi)部存在較多的漩渦結構,會對其壓力分布特性產(chǎn)生影響。為了探究泵閥耦合系統(tǒng)瞬態(tài)流動特性,在流場內(nèi)部布置壓力脈動監(jiān)測點,如圖9所示。

    圖8 壓力與速度分布Fig.8 Distribution of pressure and velocity

    圖9 監(jiān)測點分布Fig.9 Distribution of monitoring points

    圖10為導葉和蝸殼區(qū)域監(jiān)測點壓力脈動頻域曲線,其主要頻率為葉頻及其倍頻和低頻3.8 Hz。離心泵出口壓力脈動表現(xiàn)出與蝸殼監(jiān)測點一樣的頻率分布特性,主要集中在低頻3.8 Hz,其產(chǎn)生原因將在后續(xù)進行詳細說明。與導葉區(qū)域相比,蝸殼區(qū)域的壓力脈動幅值明顯下降,主要是因為流體流入蝸殼后,過流面積突然增大,流體發(fā)生強擴散現(xiàn)象,導致壓力脈動產(chǎn)生衰減。但是不同頻率的壓力脈動幅值衰減程度也不一樣,葉頻及其倍頻衰減較大,而低頻3.8 Hz衰減較小,這與其產(chǎn)生機理有關。葉頻的產(chǎn)生源于動靜干涉效應。低頻的產(chǎn)生則與導葉區(qū)域的周期性堵塞密切相關。導葉失速會引起較大的流量脈動,其引起低頻壓力脈動相較于葉頻不易耗散,因此壓力幅值下降幅度較小。

    通過前部分的研究可知,導葉流域可能存在產(chǎn)生低頻壓力脈動的周期性堵塞,進而對離心泵導葉流量進行監(jiān)測分析。由圖11可知,各個導葉流道不斷的發(fā)生阻塞和退出阻塞并且沿著順時針方向不斷傳遞,在39個葉輪旋轉周期內(nèi)存在3個明顯的失速周期,其中1個葉輪旋轉周期為葉輪旋轉1周所需要的時間。由于導葉失速團產(chǎn)生時,會阻塞流道進口的大部分面積,導致各個流道內(nèi)的流量進行重新分配,更多的水流進入到相鄰流道,流量也就發(fā)生周期性的變化。受導葉區(qū)域的周期性堵塞的影響,壓力脈動也呈現(xiàn)出周期性變化,1個導葉失速周期所需的時間為0.26 s,體現(xiàn)在頻率上為3.8 Hz。

    通過不同時刻流線圖,可以更加直觀的觀察到整個失速過程。圖12為不同時刻葉輪截面流線圖,T為一個葉輪旋轉周期,由圖可知,隨著葉輪旋轉,沿著葉輪旋轉方向導葉區(qū)域會發(fā)生周期性的堵塞,并且以13個葉輪旋轉周期為一個導葉失速周期。

    圖10 監(jiān)測點壓力脈動頻譜Fig.10 Pressure pulsation spectrum at the monitoring point

    圖11 導葉流量分布Fig.11 Flow distribution of guide vane

    流體對過流結構表面的作用力包括壓力和黏性切應力兩部分[36]。不過,由于水的黏性比較小,葉輪徑向力和閥瓣受力主要以流體壓力為主,采用直接積分法,可以直接獲取葉輪表面和閥瓣表面上所受到的作用力大小及方向,其計算公式為

    (2)

    式中:p為作用于表面A的壓力;e為表面A的法向量。

    由圖13葉輪徑向力頻譜可知,其主要頻率為低頻3.8 Hz、軸頻49.7 Hz和3倍葉頻745 Hz,并且在低頻3.8 Hz處峰值最高,表明導葉的周期性堵塞也會對其上游部件產(chǎn)生影響。

    圖12 葉輪截面流線圖Fig.12 Streamline diagram of impeller section

    圖13 葉輪徑向力頻譜Fig.13 Spectrum of the impeller radial force

    2.3 閥瓣瞬態(tài)動力學特性

    閥瓣運動特性取決于閥瓣所受的流體力、彈簧力和自身慣性力,其中彈簧力和慣性力都是結構屬性,而流體力則完全取決于閥內(nèi)流動特性。閥瓣運動過程是閥瓣所受流體力、閥瓣重力和彈簧反作用力的平衡過程中實現(xiàn)的,并且離心泵的出口壓力將直接影響閥瓣所受流體力。由圖14(a)可知,在離心泵的運轉過程中閥瓣受力始終在正負波動,且存在較為明顯的周期性,這將直接決定閥瓣的運動過程,進而引發(fā)彈簧止回閥的顫振現(xiàn)象。由圖14(b)可知,閥瓣受力在低頻3.8 Hz和10.2 Hz存在較高峰值,其中低頻3.8 Hz與離心泵出口壓力脈動的低頻成分相一致,而10.2 Hz則與彈簧止回閥內(nèi)部的彈簧質量的固有頻率相接近。

    圖14 閥瓣受力Fig.14 Force analysis during valve movement

    由圖15可知,閥瓣運動過程是一個阻尼震蕩過程。當閥瓣受到流體介質的壓力大于運動部件的重力和彈簧力之和時,此時閥瓣受力為正,閥門加速開啟。隨著閥瓣上移,其彈簧力作為反作用力逐漸增大。當閥瓣受力為負時,閥瓣速度減小直至反向,閥門漸漸關閉,彈簧力也減小,壓差引起的流體力卻快速增大。當閥瓣受力為正時,閥瓣下移速度降低并最終正向,實現(xiàn)1個震蕩周期。

    圖15 閥瓣運動規(guī)律Fig.15 The law of valve movement

    閥瓣運動過程會改變泵閥耦合系統(tǒng)的流量分配,因此需要進一步分析閥瓣位移波動引起的流量脈動。由圖16可知,空排口流量與總流量脈動幅值相當,趨勢大致相同,但總流量略微滯后于空排口流量,主出口流量脈動不明顯。由此可見,在閥瓣運動過程中,閥桿處的節(jié)流作用是引起泵閥耦合系統(tǒng)流量脈動的主要原因。

    圖16 流量分布特性Fig.16 Flow distribution characteristics

    3 結 論

    本文研究了小流量工況下泵閥耦合系統(tǒng)的非穩(wěn)態(tài)流致激勵特性,主要結論如下:

    (1) 小流量工況離心泵導葉區(qū)域沿著葉輪旋轉方向發(fā)生周期性堵塞,并以13個葉輪旋轉周期為1個導葉失速周期,受其影響,流場內(nèi)部呈現(xiàn)明顯的低頻壓力脈動。

    (2) 閥瓣運動過程是在壓差引起的流體力、運動部件重力和彈簧力的平衡過程中實現(xiàn)的,受上游離心泵出口壓力波動影響,閥瓣受力始終在正負之間波動,進而引起彈簧止回閥發(fā)生顫振現(xiàn)象。

    (3) 閥瓣運動過程中,閥桿處的節(jié)流作用是引起泵閥耦合系統(tǒng)流量脈動的主要原因。

    老师上课跳d突然被开到最大视频| 国产三级在线视频| 97在线视频观看| 国产成人影院久久av| 精品午夜福利视频在线观看一区| 久久久久免费精品人妻一区二区| 精品人妻偷拍中文字幕| 国产精品国产三级国产av玫瑰| 久久韩国三级中文字幕| 亚洲av不卡在线观看| 日韩成人伦理影院| 亚洲电影在线观看av| .国产精品久久| 国产精华一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 久久久久国内视频| 精品人妻熟女av久视频| 啦啦啦啦在线视频资源| 亚洲不卡免费看| 九九爱精品视频在线观看| 成人美女网站在线观看视频| 在线国产一区二区在线| 在线播放无遮挡| 黄色配什么色好看| 高清午夜精品一区二区三区 | 国产精品不卡视频一区二区| 黄色配什么色好看| 少妇丰满av| 日本a在线网址| 美女高潮的动态| 97超碰精品成人国产| 天天躁夜夜躁狠狠久久av| 日本精品一区二区三区蜜桃| 性色avwww在线观看| 欧美另类亚洲清纯唯美| 91午夜精品亚洲一区二区三区| 日日撸夜夜添| 1024手机看黄色片| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 国产男人的电影天堂91| 午夜精品在线福利| 天堂av国产一区二区熟女人妻| 身体一侧抽搐| 久久久久久大精品| 少妇熟女欧美另类| 国产免费男女视频| 亚洲国产精品sss在线观看| 99九九线精品视频在线观看视频| 亚洲欧美精品综合久久99| 国产探花极品一区二区| 精品久久久久久久久亚洲| 99九九线精品视频在线观看视频| 女的被弄到高潮叫床怎么办| 国产高清不卡午夜福利| 哪里可以看免费的av片| 老司机影院成人| 国产色婷婷99| 久久久久久九九精品二区国产| 亚洲熟妇中文字幕五十中出| 中国美女看黄片| 久久久国产成人精品二区| 国产伦一二天堂av在线观看| 搡女人真爽免费视频火全软件 | 日韩在线高清观看一区二区三区| 一级a爱片免费观看的视频| 日本爱情动作片www.在线观看 | 日韩欧美三级三区| 少妇人妻精品综合一区二区 | 99热精品在线国产| 亚洲精品一区av在线观看| 舔av片在线| 乱系列少妇在线播放| 国产精品野战在线观看| 日韩欧美精品v在线| 少妇裸体淫交视频免费看高清| 色视频www国产| 日韩国内少妇激情av| 亚洲精品成人久久久久久| 欧美在线一区亚洲| 国产精品一区二区免费欧美| 日韩一本色道免费dvd| 久久久久久国产a免费观看| 国产精品美女特级片免费视频播放器| 国产在线男女| 国产三级在线视频| 欧美激情在线99| 亚洲欧美精品综合久久99| 亚洲色图av天堂| 3wmmmm亚洲av在线观看| 中国美女看黄片| 亚洲天堂国产精品一区在线| 国产爱豆传媒在线观看| 校园人妻丝袜中文字幕| 人妻少妇偷人精品九色| 一个人观看的视频www高清免费观看| 国产精品免费一区二区三区在线| 丰满人妻一区二区三区视频av| 国产一区二区激情短视频| 久久久精品94久久精品| 男女之事视频高清在线观看| 春色校园在线视频观看| 成人二区视频| 青春草视频在线免费观看| 色综合色国产| 久久精品人妻少妇| 我的女老师完整版在线观看| 欧美一区二区亚洲| 久久国内精品自在自线图片| 极品教师在线视频| 少妇的逼水好多| 少妇高潮的动态图| 国产男靠女视频免费网站| 国产成年人精品一区二区| 亚洲高清免费不卡视频| 中文在线观看免费www的网站| 亚洲国产高清在线一区二区三| 久久精品影院6| 中文字幕精品亚洲无线码一区| 久久久久免费精品人妻一区二区| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 欧美高清性xxxxhd video| 国产熟女欧美一区二区| 国产高清视频在线观看网站| 热99re8久久精品国产| 国产精品亚洲美女久久久| 亚洲精品成人久久久久久| 亚洲欧美日韩卡通动漫| 国产v大片淫在线免费观看| h日本视频在线播放| 国产精品国产三级国产av玫瑰| 亚洲国产精品成人久久小说 | 成人美女网站在线观看视频| 久久久久久大精品| 日韩欧美免费精品| 日韩亚洲欧美综合| 又爽又黄a免费视频| 永久网站在线| 亚洲最大成人av| 亚洲一区高清亚洲精品| 黄色配什么色好看| 日本 av在线| 美女cb高潮喷水在线观看| 欧美潮喷喷水| 国产男靠女视频免费网站| 久久6这里有精品| 日韩欧美国产在线观看| 国内少妇人妻偷人精品xxx网站| 国产国拍精品亚洲av在线观看| 秋霞在线观看毛片| 久久久久国产网址| 免费不卡的大黄色大毛片视频在线观看 | 国产成人aa在线观看| 成人永久免费在线观看视频| 久久精品国产亚洲av香蕉五月| 亚洲av免费高清在线观看| a级一级毛片免费在线观看| 亚洲成人精品中文字幕电影| 69人妻影院| 一本精品99久久精品77| 日本a在线网址| 久久草成人影院| 亚洲人成网站在线播放欧美日韩| 哪里可以看免费的av片| 女人被狂操c到高潮| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| 色视频www国产| .国产精品久久| 波多野结衣高清作品| 午夜爱爱视频在线播放| 日韩国内少妇激情av| 99热网站在线观看| 可以在线观看的亚洲视频| 国产伦在线观看视频一区| 久久久成人免费电影| 能在线免费观看的黄片| 亚洲欧美日韩东京热| 久久久久久久久久黄片| 亚洲国产精品成人久久小说 | 精品久久久久久久久亚洲| 亚洲高清免费不卡视频| a级一级毛片免费在线观看| 人人妻人人澡欧美一区二区| 国产精品久久久久久av不卡| 欧美一级a爱片免费观看看| 亚洲精品在线观看二区| 亚洲av一区综合| 国产在视频线在精品| 夜夜看夜夜爽夜夜摸| 激情 狠狠 欧美| 直男gayav资源| 女同久久另类99精品国产91| a级毛片免费高清观看在线播放| 欧美成人精品欧美一级黄| 国产成人精品久久久久久| 超碰av人人做人人爽久久| 99久久成人亚洲精品观看| 亚洲精品亚洲一区二区| 亚洲高清免费不卡视频| 久久久久久久久中文| 午夜a级毛片| 99在线视频只有这里精品首页| 国产成人一区二区在线| 色综合亚洲欧美另类图片| 一进一出抽搐动态| 国语自产精品视频在线第100页| 国产一区二区在线av高清观看| 亚洲欧美精品综合久久99| 欧美三级亚洲精品| 18禁在线播放成人免费| 国产高清三级在线| 深夜a级毛片| 男人的好看免费观看在线视频| 国产精品人妻久久久久久| 色综合亚洲欧美另类图片| 日本撒尿小便嘘嘘汇集6| av天堂中文字幕网| 久久久精品94久久精品| 中文字幕免费在线视频6| 国产高潮美女av| 国产成人a∨麻豆精品| 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 精品无人区乱码1区二区| 成人亚洲欧美一区二区av| 国产精品爽爽va在线观看网站| 天堂影院成人在线观看| 亚洲精品日韩在线中文字幕 | 男女那种视频在线观看| 午夜福利18| 亚洲va在线va天堂va国产| 国产午夜精品论理片| 国产高清视频在线观看网站| 午夜福利成人在线免费观看| 日韩一本色道免费dvd| 欧美成人一区二区免费高清观看| 在线国产一区二区在线| 美女黄网站色视频| 91精品国产九色| 两个人的视频大全免费| 欧美zozozo另类| 小蜜桃在线观看免费完整版高清| 色综合亚洲欧美另类图片| 美女高潮的动态| 三级国产精品欧美在线观看| 69av精品久久久久久| 亚洲av中文字字幕乱码综合| 99久久精品热视频| 色视频www国产| 国产精品一区二区性色av| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站在线观看播放| 免费av毛片视频| 麻豆av噜噜一区二区三区| 成年女人看的毛片在线观看| av在线播放精品| 精品久久久久久久人妻蜜臀av| 99热网站在线观看| 成年女人毛片免费观看观看9| 日本免费一区二区三区高清不卡| 亚洲婷婷狠狠爱综合网| 亚洲成av人片在线播放无| 国产成年人精品一区二区| 精品99又大又爽又粗少妇毛片| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 国产午夜精品论理片| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 国产一区二区三区av在线 | 亚洲中文日韩欧美视频| 亚洲av美国av| 国产久久久一区二区三区| 午夜福利18| 熟女人妻精品中文字幕| 日本黄色视频三级网站网址| 99国产精品一区二区蜜桃av| 亚洲精华国产精华液的使用体验 | 最近手机中文字幕大全| 亚洲人成网站高清观看| av在线蜜桃| 色综合亚洲欧美另类图片| 99riav亚洲国产免费| 欧美最黄视频在线播放免费| 此物有八面人人有两片| 男人的好看免费观看在线视频| 黄色一级大片看看| 中国国产av一级| 免费人成在线观看视频色| 久久久久久久久大av| 国产精品伦人一区二区| 日韩一本色道免费dvd| 毛片女人毛片| 黄色一级大片看看| 淫秽高清视频在线观看| 黄色欧美视频在线观看| 香蕉av资源在线| 精品99又大又爽又粗少妇毛片| 国产精品福利在线免费观看| 国产精品美女特级片免费视频播放器| 日本-黄色视频高清免费观看| 毛片一级片免费看久久久久| 不卡视频在线观看欧美| 中文字幕熟女人妻在线| 一级毛片aaaaaa免费看小| 在线免费观看的www视频| 精品久久久久久成人av| 国产 一区精品| aaaaa片日本免费| 在线观看av片永久免费下载| 免费观看的影片在线观看| h日本视频在线播放| 日本爱情动作片www.在线观看 | 精品久久久久久久末码| 天天一区二区日本电影三级| 国产精品精品国产色婷婷| 可以在线观看毛片的网站| 亚洲一级一片aⅴ在线观看| 婷婷精品国产亚洲av| 欧美日韩乱码在线| 国产精品不卡视频一区二区| 久久久久久国产a免费观看| 国产色婷婷99| 97超视频在线观看视频| 久久精品国产鲁丝片午夜精品| 少妇被粗大猛烈的视频| 国产亚洲欧美98| 国产白丝娇喘喷水9色精品| 床上黄色一级片| 18禁黄网站禁片免费观看直播| 村上凉子中文字幕在线| 久久韩国三级中文字幕| 99久久精品一区二区三区| 国产真实伦视频高清在线观看| 久久午夜亚洲精品久久| 日日摸夜夜添夜夜添av毛片| 黑人高潮一二区| 麻豆乱淫一区二区| 搞女人的毛片| 欧美又色又爽又黄视频| 欧美xxxx黑人xx丫x性爽| 欧洲精品卡2卡3卡4卡5卡区| 亚洲自拍偷在线| 久久精品夜色国产| 永久网站在线| 亚洲精品国产成人久久av| 晚上一个人看的免费电影| 一进一出抽搐gif免费好疼| 久久久久久久午夜电影| 欧美色欧美亚洲另类二区| 人妻夜夜爽99麻豆av| 最近中文字幕高清免费大全6| 色尼玛亚洲综合影院| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品456在线播放app| 欧美一级a爱片免费观看看| 蜜臀久久99精品久久宅男| 欧美成人a在线观看| 在线播放无遮挡| av免费在线看不卡| 99热这里只有是精品在线观看| 插逼视频在线观看| 变态另类丝袜制服| 最新在线观看一区二区三区| 国产精品亚洲美女久久久| 国产精品人妻久久久影院| 老师上课跳d突然被开到最大视频| 国产白丝娇喘喷水9色精品| 99热只有精品国产| 国产中年淑女户外野战色| 看十八女毛片水多多多| 日本五十路高清| 国产精品久久视频播放| 日本一本二区三区精品| 老熟妇乱子伦视频在线观看| 欧美日韩精品成人综合77777| 久久午夜亚洲精品久久| 午夜免费男女啪啪视频观看 | 一本精品99久久精品77| 最好的美女福利视频网| 国内精品宾馆在线| 免费av观看视频| 精品一区二区免费观看| 99久久精品热视频| 中文字幕精品亚洲无线码一区| 亚洲第一区二区三区不卡| 99久国产av精品国产电影| 欧美区成人在线视频| 亚洲欧美清纯卡通| 亚洲国产高清在线一区二区三| 亚洲欧美清纯卡通| 91久久精品国产一区二区成人| 精品久久久久久久久久免费视频| 三级国产精品欧美在线观看| 久久午夜福利片| 两性午夜刺激爽爽歪歪视频在线观看| 女生性感内裤真人,穿戴方法视频| 国产欧美日韩精品一区二区| 亚洲精品国产成人久久av| 91在线观看av| 在线观看美女被高潮喷水网站| 欧美人与善性xxx| 亚洲av第一区精品v没综合| 日韩精品中文字幕看吧| 国产精品久久视频播放| 日本黄色视频三级网站网址| 麻豆国产av国片精品| 两个人的视频大全免费| 最近的中文字幕免费完整| 亚洲精品在线观看二区| 久久精品国产99精品国产亚洲性色| 欧美日韩精品成人综合77777| 免费黄网站久久成人精品| 青春草视频在线免费观看| 亚洲内射少妇av| 99久久精品一区二区三区| 91久久精品国产一区二区三区| 日韩成人伦理影院| 嫩草影院新地址| 欧美不卡视频在线免费观看| 性色avwww在线观看| 不卡一级毛片| 乱系列少妇在线播放| 午夜福利在线在线| 亚洲精品456在线播放app| 久久亚洲精品不卡| av在线老鸭窝| 国产精品三级大全| 日日摸夜夜添夜夜添av毛片| 亚洲七黄色美女视频| 看十八女毛片水多多多| 男插女下体视频免费在线播放| 两个人的视频大全免费| 麻豆乱淫一区二区| 精品久久久久久久久亚洲| 色尼玛亚洲综合影院| 五月玫瑰六月丁香| 日日摸夜夜添夜夜添小说| 最近视频中文字幕2019在线8| 免费看光身美女| 六月丁香七月| 亚洲电影在线观看av| 亚洲成人久久爱视频| 村上凉子中文字幕在线| 性插视频无遮挡在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看| videossex国产| 免费观看在线日韩| 午夜激情福利司机影院| eeuss影院久久| 又爽又黄a免费视频| 长腿黑丝高跟| 看黄色毛片网站| 日本免费a在线| 欧美最新免费一区二区三区| 亚洲欧美中文字幕日韩二区| 深夜a级毛片| 久久精品国产鲁丝片午夜精品| 亚洲国产欧洲综合997久久,| 中国美白少妇内射xxxbb| 久久久久久大精品| 国产欧美日韩精品一区二区| 精品国产三级普通话版| 亚洲无线在线观看| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久精品电影| 国产乱人偷精品视频| av免费在线看不卡| 国产在视频线在精品| 国产伦精品一区二区三区视频9| 国产精品久久久久久久久免| 婷婷精品国产亚洲av在线| 最近视频中文字幕2019在线8| 午夜免费男女啪啪视频观看 | 一级黄片播放器| 91午夜精品亚洲一区二区三区| 好男人在线观看高清免费视频| 大香蕉久久网| 亚洲av第一区精品v没综合| 99久久九九国产精品国产免费| 免费观看的影片在线观看| 亚洲av中文av极速乱| 如何舔出高潮| 亚州av有码| 日韩欧美在线乱码| 欧美zozozo另类| 国产成人精品久久久久久| 久久99热6这里只有精品| 99久久精品热视频| 精品久久久久久久人妻蜜臀av| 狠狠狠狠99中文字幕| 亚洲中文字幕日韩| 国产视频内射| 日韩欧美一区二区三区在线观看| 少妇的逼水好多| 卡戴珊不雅视频在线播放| 人妻夜夜爽99麻豆av| 男人狂女人下面高潮的视频| 亚洲av免费在线观看| 国产亚洲精品综合一区在线观看| 久久天躁狠狠躁夜夜2o2o| 日韩人妻高清精品专区| 内地一区二区视频在线| 久久久久久国产a免费观看| 久久久久国产精品人妻aⅴ院| 久久午夜亚洲精品久久| 久久国产乱子免费精品| 丝袜喷水一区| 亚洲av一区综合| 国产视频一区二区在线看| 久久人人爽人人片av| 久久这里只有精品中国| 可以在线观看毛片的网站| 岛国在线免费视频观看| 国产高潮美女av| 老司机福利观看| 久久精品影院6| 成人高潮视频无遮挡免费网站| 精品欧美国产一区二区三| 深夜a级毛片| 国产在线精品亚洲第一网站| 欧美日韩乱码在线| 亚洲一区二区三区色噜噜| 91精品国产九色| 亚洲av成人av| 大香蕉久久网| 国产麻豆成人av免费视频| 久久这里只有精品中国| 久久久久久久久中文| 久久天躁狠狠躁夜夜2o2o| 亚洲精品日韩在线中文字幕 | 在线观看免费视频日本深夜| 欧美色欧美亚洲另类二区| 久久精品国产清高在天天线| 国产午夜精品论理片| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 一本精品99久久精品77| 亚洲av一区综合| 午夜福利在线观看免费完整高清在 | 久久精品国产99精品国产亚洲性色| 97超视频在线观看视频| 免费看光身美女| 久久久国产成人免费| 3wmmmm亚洲av在线观看| 国产精品野战在线观看| 久久天躁狠狠躁夜夜2o2o| 波多野结衣高清无吗| 久久久久精品国产欧美久久久| 激情 狠狠 欧美| 你懂的网址亚洲精品在线观看 | 非洲黑人性xxxx精品又粗又长| 国产高清激情床上av| 日韩欧美在线乱码| 男插女下体视频免费在线播放| 日韩欧美精品v在线| 日韩大尺度精品在线看网址| 91久久精品国产一区二区成人| 天天躁夜夜躁狠狠久久av| 国内久久婷婷六月综合欲色啪| 舔av片在线| 午夜福利在线观看免费完整高清在 | 别揉我奶头~嗯~啊~动态视频| 97碰自拍视频| 国产精品日韩av在线免费观看| 精品午夜福利视频在线观看一区| 男女下面进入的视频免费午夜| 亚洲国产精品国产精品| 一进一出好大好爽视频| 国产又黄又爽又无遮挡在线| 五月玫瑰六月丁香| 尾随美女入室| 色播亚洲综合网| 日本黄色视频三级网站网址| 亚洲成人久久爱视频| 国产真实伦视频高清在线观看| 女的被弄到高潮叫床怎么办| 99久久精品热视频| 国产色婷婷99| 国产精品精品国产色婷婷| 大香蕉久久网| 99久久精品国产国产毛片| 日本免费一区二区三区高清不卡| 美女高潮的动态| 久久精品国产亚洲av香蕉五月| 啦啦啦观看免费观看视频高清| 国产精品精品国产色婷婷| 色综合色国产| 给我免费播放毛片高清在线观看| 欧美一区二区精品小视频在线| 国产亚洲欧美98| 婷婷六月久久综合丁香| 大又大粗又爽又黄少妇毛片口| 国产精品女同一区二区软件| 又黄又爽又免费观看的视频| 亚洲精品国产av成人精品 | 久久草成人影院| 97在线视频观看| 国产视频一区二区在线看| 91精品国产九色| a级毛片免费高清观看在线播放| 你懂的网址亚洲精品在线观看 | 麻豆av噜噜一区二区三区| 哪里可以看免费的av片| 婷婷精品国产亚洲av| 最新在线观看一区二区三区| 精品久久久久久久久av| 九九在线视频观看精品| 老司机福利观看| 我要看日韩黄色一级片| 国产一级毛片七仙女欲春2| 九九爱精品视频在线观看| 免费大片18禁|