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

    換熱器內(nèi)彈性管束流體組合誘導振動響應的數(shù)值分析

    2015-03-07 00:34:02季家東葛培琪畢文波
    西安交通大學學報 2015年9期
    關鍵詞:管程殼程管束

    季家東,葛培琪,2,畢文波

    ?

    換熱器內(nèi)彈性管束流體組合誘導振動響應的數(shù)值分析

    季家東1,葛培琪1,2,畢文波1

    彈性管束換熱器通過殼、管程兩場流體的組合誘導引起內(nèi)部彈性管束的振動,從而實現(xiàn)強化傳熱[1-5]。然而,這種流體誘導而引起的彈性管束振動會導致內(nèi)部管束的疲勞破壞,影響換熱器的使用壽命。彈性管束換熱器的設計應遵循這樣一條主線:在滿足強化傳熱的參數(shù)范圍內(nèi),對流體誘導彈性管束的振動進行合理誘發(fā)和適當控制,同時兼顧管束的疲勞壽命,確保在強化傳熱的同時不發(fā)生疲勞破壞。因此,對實際運行條件下管束振動的合理誘發(fā)成為彈性管束換熱器設計的關鍵[6]。由于管束結構、殼程管程流場和換熱器工作條件的復雜性,現(xiàn)階段關于彈性管束流體誘導振動的研究大多是實驗研究,具有較大的局限性,因此關于殼、管程兩場流體組合誘導下彈性管束振動響應的分析,對研究強化傳熱機理、管束結構優(yōu)化以及實現(xiàn)振動的有效控制等都具有重要意義。

    為了誘發(fā)彈性管束的振動,姜波在文獻[1-2]提出的彈性管束基礎上進行了適當改進,設計了一種電機驅(qū)動的管程流體誘導脈動裝置,建立了含管程脈動裝置的恒熱流傳熱試驗臺,并對改進型彈性管束的管外傳熱特性進行了實驗研究[7]。研究表明,管內(nèi)流體介質(zhì)對彈性管束的振動特性影響較大,振動特性增強使換熱特性增強,在所討論的參數(shù)范圍內(nèi)彈性管束的管外平均表面?zhèn)鳠嵯禂?shù)可以達到固定管束的3倍以上。閆柯等依據(jù)彈性管束的工作原理,提出了一種新型空間錐螺旋彈性管束,并實驗測試了該管束在一定殼程流速范圍內(nèi)的縱向振動頻率[8-9]。研究表明,管束的工作模態(tài)頻率隨流速的增加而降低,在一定流速范圍內(nèi),管束振幅增長較為緩慢,超過一定流速時,管束的振幅隨著流速的增加迅速增大。此外,宿艷彩對單排彈性管束在殼程流體誘導下的振動響應進行了實驗測試,得到了不同殼程流速條件下單排彈性管束不銹鋼連接體的振動頻率[10]。研究表明,殼程流速對不銹鋼連接體振動頻率的影響較大,且低速流體誘導下彈性管束的振動存在諧頻。

    本文基于流固耦合問題的弱耦合法,采用通用CFD分析軟件CFX和ANSYS Workbench軟件的瞬態(tài)動力學分析模塊,對彈性管束在不同流速的殼、管程兩場流體組合誘導下的振動響應進行了數(shù)值研究。

    1 計算方法

    1.1 幾何模型

    圖1所示為彈性管束的幾何結構。彈性管束由4根純銅彎管(彎管半徑分別為R1、R2、R3、R4;截面外徑為d;壁厚為δ)和2塊不銹鋼連接體(III、IV)組成。I、II兩處為固定端,不銹鋼連接體III、IV為自由端,在換熱器實際工作時,管程流體由I處端口流入,II處端口流出。

    圖1 彈性管束結構示意圖

    圖2所示為彈性管束的殼程流體計算域、管程流體計算域和結構計算域網(wǎng)格。圖2a中,x方向為殼程流體的流動方向。圖2c中,A、B兩點是設立在不銹鋼連接體III、IV上的監(jiān)測點,用來檢測殼、管程流體組合誘導彈性管束的振動情況。殼、管程流體域網(wǎng)格采用網(wǎng)格劃分軟件ICEM進行劃分,均為結構化網(wǎng)格;結構部分網(wǎng)格采用Workbench軟件的網(wǎng)格劃分模塊Mechanical平臺進行劃分,包含結構化網(wǎng)格(銅彎管)及非結構化網(wǎng)格(不銹鋼連接體)。其中,殼程流體計算域網(wǎng)格包含329 708個單元,309 305個節(jié)點;管程流體計算域網(wǎng)格包含52 940個單元,47 880個節(jié)點;結構計算域網(wǎng)格包含9 366個單元,54 777個節(jié)點。經(jīng)過試算,進一步增加網(wǎng)格數(shù)量或管束周圍網(wǎng)格密度對計算結果的影響不大,網(wǎng)格獨立性得到驗證。

    (a)殼程流體計算域網(wǎng)格

    (b)管程流體計算域網(wǎng)格 (c)結構計算域網(wǎng)格圖2 流體和結構計算域網(wǎng)格

    1.2 材料屬性和結構參數(shù)

    彈性管束彎管部分的材料為純銅,連接體部分的材料為不銹鋼,具體材料屬性如表1所示。

    表1 彈性管束的材料屬性

    為了研究彈性管束在不同流速的殼、管程兩場流體組合誘導下的振動特性,分析過程中所用彈性管束的具體尺寸如表2所示。

    表2 彈性管束的具體尺寸

    1.3 數(shù)值分析方法

    基于本文的研究問題,將整個求解域分成流體域(包含殼程流體域和管程流體域)和結構域,如圖2所示。流固耦合求解采用弱耦合方法在每個時間步內(nèi)按順序進行獨立求解。流固耦合計算時,分別將管程流體域的外表面和殼程流體域的內(nèi)表面設置為流固耦合面。結構部分將彈性管束內(nèi)、外表面設置為流固耦合面,分別與管程流體域的外表面和殼程流體域的內(nèi)表面相對應。流體域與結構域間的數(shù)據(jù)傳遞通過流固耦合面完成。

    殼、管程流體域選擇通用CFD分析軟件CFX進行計算。分析過程中,殼、管程流速變化范圍均設為0.6~1.2 m/s(Re>2 300),均屬湍流流動,選用標準的k-ε模型對湍流進行模擬。邊界條件:入口邊界類型為Inlet,并給定入口速度;出口邊界類型為Outlet,設置出口相對靜壓為0 Pa,且計算環(huán)境為標準大氣壓(101.325 kPa)。結構部分采用ANSYS Workbench軟件的瞬態(tài)動力學分析模塊進行計算,兩固定端I、II處截面設為固定約束。

    考慮流體與結構交界面上的能量損耗,得到流固耦合系統(tǒng)的有限元波動方程如下

    (1)

    利用哈密爾頓(Hamilton)原理建立結構的整體運動方程??紤]流體的作用,在液體中離散后彈性體的結構力學方程為

    (2)

    式中:M為質(zhì)量矩陣;C為阻尼矩陣;K為剛度矩陣;下標f、p分別表示流體和結構;p為流體動壓力;u為速度矢量;RT為流固耦合交界面上的耦合矩陣;Fp為結構內(nèi)部作用力;Ff為流體對結構的作用力。

    由式(1)和式(2)可得到完整描述流固耦合問題的有限元離散方程

    (3)

    求解式(3)可得到完全流固耦合條件下流體誘導彈性管束振動的頻率和振型。

    流固耦合計算時,首先采用CFX軟件進行殼程和管程流體域的計算,得到兩耦合面(管程流體域的外表面和殼程流體域的內(nèi)表面)的壓力分布,然后將此壓力分布作為初始條件,采用ANSYS Workbench軟件對結構進行瞬態(tài)動力學分析,得到兩耦合面的位移,以此作為下一個時間步內(nèi)流體域計算的邊界條件,如此交替迭代,直至完成。

    2 管束振動響應分析

    為了對數(shù)值分析的正確性和準確性進行驗證,基于文獻[10]實驗用彈性管束的結構參數(shù),建立與之相應的數(shù)學模型,對不同殼程流速vshell(0.2、0.4 m/s)條件下不銹鋼連接體III、IV上兩監(jiān)測點在x方向(沿殼程流速方向)的振動頻率f進行了求解,結果如表3所示。其中,殼程流體介質(zhì)為水,彈性管束內(nèi)無流體介質(zhì)。

    表3 數(shù)值解與實驗值的對比情況

    從表3可以看出,不同殼程流速條件下的數(shù)值解與實驗數(shù)據(jù)高度一致,求解參數(shù)范圍內(nèi)的最大相對誤差僅為2.14%,數(shù)值求解方法的正確性和準確性得到了驗證。

    2.1 流體速度對振動響應的影響

    為了研究殼、管程兩場流體速度對彈性管束振動響應的影響,研究了不同流速條件下,兩場流體組合誘導下彈性管束的振動特性。其中,殼、管程流體具有相等的流速,即vshell=vtube,且流體介質(zhì)為水。

    圖3所示為不同流速(0.8、1.2 m/s)時,不銹鋼連接體III、IV上監(jiān)測點A、B在x、y、z方向的振動位移S隨計算時間t的變化情況。從圖3可以看出:在殼、管程兩場流體的組合誘導下,彈性管束能夠?qū)崿F(xiàn)均勻的振動;當流速較高(1.2 m/s)時,不銹鋼連接體IV(監(jiān)測點B)在z方向的振動異常劇烈;由于流體的沖擊和管束重力的影響,當流速較低(0.8 m/s)時,管束x方向的振動平衡位置較低,隨著流速的增加,管束x方向的振動平衡位置升高;當流速較低(0.8 m/s)時,監(jiān)測點A在x方向的振動位移曲線存在明顯的“雙峰”現(xiàn)象。因為低速流體對彈性管束的沖擊力較小,重力對彈性管束在x方向振動的影響較大,此時彈性管束振動控制方程中的二次方非線性項所起的作用較大,這是振動位移曲線存在“雙峰”現(xiàn)象的原因。當流速較高(1.2 m/s)時,流體的沖擊力增強,弱化了重力對振動的影響,此時振動控制方程中的二次方非線性項所起的作用較小,所以在振動位移波動圖中未觀察到明顯“雙峰”現(xiàn)象。

    圖4所示為不同流速時,不銹鋼連接體IV上監(jiān)測點B的振動位移頻譜圖。從圖4可以看出:由于殼、管程兩場流體的組合誘導,監(jiān)測點各方向的振動主頻、諧頻的大小一致,當流速較低(0.8 m/s)時,相較于監(jiān)測點x方向的主頻而言,二倍諧頻所起的作用明顯;隨著流速的增加,監(jiān)測點在x、y、z方向的振動頻率均逐漸增加,在z方向的振動幅值A也逐漸增加,而在x、y方向的振動幅值A則先增加后減小;在殼、管程流體組合誘導下,監(jiān)測點的振動以z方向的振動為主,此時彈性管束的陣型為面內(nèi)振動。

    (a)vshell=vtube=0.8 m/s

    (b)vshell=vtube=1.2 m/s

    (a)x方向的頻譜圖 (b)y方向的頻譜圖 (c)z方向的頻譜圖圖4 不同流速時監(jiān)測點B的振動位移頻譜圖

    2.2 殼程流速對振動響應的影響

    為了研究殼程流體速度對彈性管束振動響應的影響,對不同殼程流速條件下彈性管束的流體誘導振動特性進行了研究,計算時管程流速固定為0.8 m/s。

    圖5所示為不同殼程流速(0.6、1.2 m/s)時,不銹鋼連接體III、IV上監(jiān)測點A、B在x、y、z方向的振動位移隨計算時間的變化情況。從圖5可以看出:殼程流速變化對管束振動響應的影響較大,在圖示速度條件下,監(jiān)測點各方向的振動幅值和頻率均隨著殼程流速的增加而增加,當流速為1.2 m/s時,不銹鋼連接體IV(監(jiān)測點B)在z方向的振動異常劇烈;圖5b與圖3b相比,監(jiān)測點各方向振動位移曲線的差距不大,說明彈性管束的振動主要由殼程流體引起;殼程流體速度的變化對彈性管束x方向振動平衡位置的影響較大,隨著殼程流速的增加,彈性管束x方向的振動平衡位置升高。

    圖6所示為不同殼程流速時,不銹鋼連接體IV上監(jiān)測點B的振動位移頻譜圖。從圖6可以看出:不同殼程流速時,不銹鋼連接體IV上監(jiān)測點B各方向的振動位移頻譜圖與圖4所示殼、管程流速相同時的頻譜圖類似,再次表明管程流體誘導對彈性管束振動響應的影響較小,彈性管束的振動主要由殼程流體引起。

    2.3 管程流速對振動響應的影響

    為了研究管程流體速度對彈性管束振動響應的影響,并進一步驗證殼程流體是引起彈性管束振動的主要誘因,對不同管程流速條件下彈性管束的流體誘導振動特性進行了研究,計算時殼程流速固定為0.8 m/s。

    (a)vshell=0.6 m/s

    (b)vshell=1.2 m/s

    (a)x方向的頻譜圖 (b)y方向的頻譜圖 (c)z方向的頻譜圖圖6 不同殼程流速時監(jiān)測點B的振動位移頻譜圖

    (a)vtube=0.6 m/s

    (b)vtube=1.2 m/s

    圖7所示為不同管程流速(0.6、1.2 m/s)時,不銹鋼連接體III、IV上監(jiān)測點A、B在x、y、z方向的振動位移隨計算時間的變化情況。將圖7與圖3a對比可以看出,幾種管程流速下同監(jiān)測點、同方向的位移曲線基本一致,說明管程流體誘導對彈性管束振動響應的影響較小,進一步驗證了殼程流體是引起彈性管束振動的主要誘因。

    3 結 論

    本文基于流固耦合問題的弱耦合法,對彈性管束在不同流速的殼、管程兩場流體組合誘導下的振動響應進行了研究。通過與實驗數(shù)據(jù)對比,驗證了數(shù)值求解方法的正確性和準確性,主要結論如下。

    (1)隨著殼、管程流速的增加,監(jiān)測點在各方向的振動頻率均逐漸增加,在z方向的振動幅值逐漸增加,在x、y方向振幅先增加后減小。綜合看來,監(jiān)測點以z方向振動為主,說明在殼、管程流體的組合誘導下,管束的振動主要體現(xiàn)為面內(nèi)振動。

    (2)由于殼、管程流體的組合誘導作用,監(jiān)測點各方向的振動主頻、諧頻的大小一致。殼程流速較低(0.8 m/s)時,監(jiān)測點沿殼程流速方向(x方向)二倍諧頻所起的作用明顯。

    (3)由于流體的沖擊和管束重力的影響,當殼程流速較低時,管束沿殼程流速方向(x方向)的振動平衡位置較低,隨著流速的增加,管束沿殼程流速方向的振動平衡位置逐漸升高。

    (4)殼程流速變化對管束振動響應的影響較大,管程流速變化對管束振動響應的影響較小,說明彈性管束的振動主要由殼程流體引起。

    [1] 程林. 彈性管束換熱器原理與應用 [M]. 北京: 科學出版社, 2001: 34-67.

    [2] CHENG Lin, LUAN Ting, DU Wenjing, et al. Heat transfer enhancement by flow-induced vibration in heat exchangers [J]. International Journal of Heat and Mass Transfer, 2009, 52(3): 1053-1057.

    [3] YAN Ke, GE Peiqi, HONG Jun. Experimental study of shell side flow-induced vibration of conical spiral tube bundle [J]. Journal of Hydrodynamics, 2013, 25(5): 695-701.

    [4] JIANG Bo, TIAN Maocheng, CHENG Lin. Research on vibration and heat transfer characteristics of an optimized structure of elastic tube bundle [J]. Advanced Materials Research, 2013, 816/817: 260-265.

    [5] 閆柯, 葛夢然, 高軍, 等. 空間錐螺旋管束流體誘導振動換熱器及性能分析 [J]. 西安交通大學學報, 2011, 45(11): 22-26. YAN Ke, GE Mengran, GAO Jun, et al. Performance of flow-induced vibration heat exchanger with conical spiral tube bundle [J]. Journal of Xi’an Jiaotong University, 2011, 45(11): 22-26.

    [6] 閆柯, 葛培琪, 宿艷彩. 基于復模態(tài)的空間錐螺旋管管內(nèi)流固耦合振動特性的有限元分析 [J]. 振動與沖擊, 2011, 31(1): 204-207. YAN Ke, GE Peiqi, SU Yancai. FEM analysis for fluid-structure interaction in a dimensional helical tube base on complex mode theory [J]. Journal of Vibration and Shock, 2011, 31(1): 204-207.

    [7] 姜波, 郝衛(wèi)東, 劉福國, 等. 流體脈動對新型彈性管束傳熱特性影響的實驗研究 [J]. 振動與沖擊, 2012, 31(10): 59-63. JIANG Bo, HAO Weidong, LIU Fuguo, et al. Experiments on effects of fluid pulsation on heat transfer characteristics of the new type elastic tube bundle [J]. Journal of Vibration and Shock, 2012, 31(10): 59-63.

    [8] 閆柯. 換熱器內(nèi)錐螺旋彈性管束振動與傳熱特性研究 [D]. 濟南: 山東大學, 2012.

    [9] YAN Ke, GE Peiqi, HU Ruirong, et al. Heat transfer and resistance characteristics of conical spiral tube bundle based on field synergy principle [J]. Chinese Journal of Mechanical Engineering, 2012, 25(2): 370-376.

    [10]宿艷彩. 彈性管束流體誘導振動及換熱特性研究 [D]. 濟南: 山東大學, 2012.

    [本刊相關文獻鏈接]

    張凡,李兆輝,李曉宇,等.不同材料翅片管換熱器特性的實驗研究.2015,49(5):62-67.[doi:10.7652/xjtuxb201505010]

    文鍵,楊輝著,王斯民,等.旋梯式螺旋折流板換熱器優(yōu)化結構的數(shù)值模擬.2014,48(11):8-14.[doi:10.7652/xjtuxb 201411002]

    楊衛(wèi)衛(wèi),曹興起,何雅玲,等.一種帶中間換熱器的雙熱源高溫熱泵系統(tǒng).2014,48(11):70-74.[doi:10.7652/xjtuxb 201411012]

    文鍵,楊輝著,杜冬冬,等.螺旋折流板換熱器換熱強化的數(shù)值研究.2014,48(9):43-48.[doi:10.7652/xjtuxb201409008]

    韓輝,何雅玲,唐桂華.硫酸蒸氣在換熱器表面沉積特性的數(shù)值研究.2013,47(11):7-13.[doi:10.7652/xjtuxb201311002]

    (編輯 荊樹蓉)

    (1.山東大學機械工程學院,250061,濟南;2.山東大學高效潔凈機械制造教育部重點實驗室,250061,濟南)

    為了探討換熱器內(nèi)彈性管束在流體誘導下的振動特性,采用流固耦合的弱耦合法,研究了彈性管束在不同流速的殼程和管程兩場流體組合誘導下的振動響應,并對比分析了殼程流體和管程流體對彈性管束振動響應的影響。研究表明:當殼程和管程流速一定時,彈性管束兩監(jiān)測點各方向的振動主頻、諧頻大小一致,且振動主要體現(xiàn)為面內(nèi)振動。殼程流速較低時,監(jiān)測點沿殼程流速方向的位移曲線存在明顯的二倍諧頻。殼程流速變化對管束振動響應的影響較大,管程流速變化對管束振動響應的影響較小,說明彈性管束的振動主要由殼程流體引起。由于流體的沖擊和管束重力的影響,當殼程流速較低時,管束沿殼程流速方向的振動平衡位置較低,隨著流速的增加,管束沿殼程流速方向的振動平衡位置逐漸上移。

    彈性管束;換熱器;流固耦合;流體誘導振動

    Numerical Analysis on the Combined Flow Induced Vibration Response of Elastic Tube Bundle in Heat Exchanger

    JI Jiadong1,GE Peiqi1,2,BI Wenbo1

    (1. School of Mechanical Engineering, Shandong University, Jinan 250061, China; 2. Key Laboratory of High-efficiency and Clean Mechanical Manufacture at Ministry of Education, Shandong University, Jinan 250061, China)

    The vibration characteristics of the elastic tube bundle in heat exchanger, under the combined induction of shell-side and tube-side cross flows at different inlet water velocities are studied using weak coupling method of fluid-structure interaction. And the influences of shell-side and tube-side cross flows on the vibration responses of the elastic tube bundle are also analyzed comparatively. Numerical results show that the flow-induced dominant and harmonic frequency vibration frequencies, at the two monitor points of the elastic tube bundle and under the same water velocity, are uniform in all directions; and the flow-induced vibration is mainly an in-plane vibration. There is an obvious second harmonic frequency in the shell-side water flow direction at the two monitor points when the water velocity is relatively low. Compared with the change of the shell-side water velocity, the variation of the tube-side inlet flow velocity has less effect on the vibration response. It shows that the vibration is mainly induced by the shell-side cross flow. Furthermore, the vibration balance position of the elastic tube bundle in shell-side water flow direction is also low when the velocity is low, due to the influences of the gravity of tube bundle and the impact force of fluid; and the vibration balance position gradually moves up with the increase of the shell-side water velocity.

    elastic tube bundle; heat exchanger; fluid-structure interaction; flow-induced vibration

    2015-03-16。 作者簡介:季家東(1982—),男,博士生;葛培琪(通信作者),男,教授,博士生導師。 基金項目:國家自然科學基金資助項目(51475268);國家重點基礎研究發(fā)展計劃資助項目(2007CB206903)。

    時間:2015-06-04

    http:∥www.cnki.net/kcms/detail/61.1069.T.20150604.1711.002.html

    10.7652/xjtuxb201509005

    TH123;TK172

    A

    0253-987X(2015)09-0024-06

    猜你喜歡
    管程殼程管束
    管間距對橫掠管束換熱影響及非線性現(xiàn)象分析
    能源工程(2020年6期)2021-01-26 00:55:16
    管殼式換熱器殼程的傳熱強化探討
    三葉膨脹管換熱器殼程強化傳熱的數(shù)值研究*
    新能源進展(2020年1期)2020-03-09 03:50:16
    放開管束孩子的手
    海峽姐妹(2019年6期)2019-06-26 00:52:38
    管殼式換熱器管束拆卸問題與建議
    加氫精制U形管式換熱器管束泄漏分析與對策
    傾斜折流柵式換熱器殼程熱力特性研究
    換熱管束排列對換熱器殼程流動與傳熱影響的研究
    多管程布置微通道分液冷凝器的熱力性能
    制冷學報(2015年5期)2015-09-01 07:24:22
    利用管程概念求解哲學家進餐問題
    欧美成人一区二区免费高清观看| 婷婷精品国产亚洲av在线| 小蜜桃在线观看免费完整版高清| 91狼人影院| 免费av不卡在线播放| 我的老师免费观看完整版| 国产男靠女视频免费网站| www.色视频.com| 亚洲内射少妇av| 六月丁香七月| 一个人观看的视频www高清免费观看| 禁无遮挡网站| 麻豆乱淫一区二区| 韩国av在线不卡| 亚洲欧美中文字幕日韩二区| 日本撒尿小便嘘嘘汇集6| 麻豆国产av国片精品| 春色校园在线视频观看| 极品教师在线视频| 老熟妇乱子伦视频在线观看| 色综合亚洲欧美另类图片| 亚洲av熟女| 成人鲁丝片一二三区免费| 午夜日韩欧美国产| 亚洲美女黄片视频| 深爱激情五月婷婷| 成年免费大片在线观看| 久久精品国产亚洲av涩爱 | 99热网站在线观看| 午夜福利视频1000在线观看| 精品久久久噜噜| 中文亚洲av片在线观看爽| 久久这里只有精品中国| 中国美女看黄片| 国产精品久久电影中文字幕| 伊人久久精品亚洲午夜| 国产一区二区三区av在线 | 美女高潮的动态| 免费av观看视频| 日韩精品青青久久久久久| 欧美一区二区国产精品久久精品| 国产视频内射| 日韩国内少妇激情av| 亚洲自拍偷在线| 欧美性猛交╳xxx乱大交人| 最近2019中文字幕mv第一页| 精品久久久久久久人妻蜜臀av| 国产精品综合久久久久久久免费| 97人妻精品一区二区三区麻豆| 成人av一区二区三区在线看| 亚洲人与动物交配视频| 热99re8久久精品国产| 1024手机看黄色片| 免费看av在线观看网站| a级毛色黄片| 欧美3d第一页| 性欧美人与动物交配| 免费看日本二区| av黄色大香蕉| 欧美最黄视频在线播放免费| 中文字幕av在线有码专区| 国产伦一二天堂av在线观看| 黄色日韩在线| 日本成人三级电影网站| 国产免费一级a男人的天堂| 97碰自拍视频| 欧美日韩精品成人综合77777| 国产国拍精品亚洲av在线观看| 亚洲精品国产av成人精品 | 午夜a级毛片| 亚洲国产精品成人综合色| 精品日产1卡2卡| 国产一区二区三区av在线 | 一个人观看的视频www高清免费观看| 亚洲av电影不卡..在线观看| 一级黄色大片毛片| 国产男人的电影天堂91| 久久精品国产亚洲网站| 黄色配什么色好看| 丰满的人妻完整版| 麻豆成人午夜福利视频| 精品国内亚洲2022精品成人| 波多野结衣高清无吗| 91久久精品国产一区二区三区| 三级毛片av免费| 久久久久久久久中文| 日日啪夜夜撸| 亚洲av一区综合| 人妻久久中文字幕网| 亚洲经典国产精华液单| 国产精品,欧美在线| 久久久久九九精品影院| 久久久久久久亚洲中文字幕| av天堂在线播放| 俄罗斯特黄特色一大片| 欧美中文日本在线观看视频| 看十八女毛片水多多多| 国产男人的电影天堂91| 久久精品国产99精品国产亚洲性色| 精品少妇黑人巨大在线播放 | 日日啪夜夜撸| 性插视频无遮挡在线免费观看| 午夜久久久久精精品| 精品熟女少妇av免费看| 淫秽高清视频在线观看| 午夜福利成人在线免费观看| 欧美最新免费一区二区三区| 亚洲熟妇熟女久久| 久久久国产成人精品二区| 日本熟妇午夜| 搡女人真爽免费视频火全软件 | 亚洲性久久影院| 成人三级黄色视频| 国产成人影院久久av| 99热精品在线国产| 日韩中字成人| 久久九九热精品免费| 大又大粗又爽又黄少妇毛片口| 一进一出抽搐gif免费好疼| 亚洲熟妇中文字幕五十中出| 99九九线精品视频在线观看视频| 亚洲18禁久久av| 91av网一区二区| 久久这里只有精品中国| 色综合亚洲欧美另类图片| 国产精品女同一区二区软件| 少妇裸体淫交视频免费看高清| 91久久精品国产一区二区成人| 午夜福利视频1000在线观看| 成熟少妇高潮喷水视频| 国内精品久久久久精免费| .国产精品久久| 亚洲av中文av极速乱| 免费在线观看影片大全网站| 18禁黄网站禁片免费观看直播| 免费av毛片视频| 久久久久九九精品影院| 精品久久久噜噜| 两个人的视频大全免费| 看非洲黑人一级黄片| 久久婷婷人人爽人人干人人爱| 波多野结衣高清作品| 欧美zozozo另类| 日本撒尿小便嘘嘘汇集6| 一级毛片aaaaaa免费看小| 熟女电影av网| 久久精品国产自在天天线| 欧美性感艳星| 国产片特级美女逼逼视频| 日本三级黄在线观看| 欧美人与善性xxx| 男女啪啪激烈高潮av片| АⅤ资源中文在线天堂| 热99re8久久精品国产| 国产精品乱码一区二三区的特点| 高清日韩中文字幕在线| 国产视频内射| 国产一级毛片七仙女欲春2| 国产精品日韩av在线免费观看| 麻豆av噜噜一区二区三区| 又黄又爽又免费观看的视频| 亚洲av.av天堂| 国产成人aa在线观看| 亚洲,欧美,日韩| 九九久久精品国产亚洲av麻豆| 亚洲四区av| 日日摸夜夜添夜夜添av毛片| 国产精品一区二区三区四区久久| 极品教师在线视频| 日韩人妻高清精品专区| 亚洲av第一区精品v没综合| 国产在视频线在精品| 91久久精品电影网| 嫩草影院入口| 在线播放无遮挡| 九九久久精品国产亚洲av麻豆| av女优亚洲男人天堂| 97超碰精品成人国产| 国产午夜福利久久久久久| 国产亚洲精品久久久久久毛片| 一个人看的www免费观看视频| 亚洲国产高清在线一区二区三| 色在线成人网| 国产精品爽爽va在线观看网站| 国产精品嫩草影院av在线观看| 国产在视频线在精品| 变态另类成人亚洲欧美熟女| 久久精品久久久久久噜噜老黄 | 成年版毛片免费区| 少妇被粗大猛烈的视频| 日日摸夜夜添夜夜爱| 国产精品美女特级片免费视频播放器| 精品久久久久久成人av| 中文在线观看免费www的网站| 一个人看视频在线观看www免费| 亚洲精品在线观看二区| 久久久久久久午夜电影| 一本一本综合久久| 精品人妻一区二区三区麻豆 | 亚洲图色成人| 特级一级黄色大片| 亚洲自拍偷在线| 69人妻影院| 在线免费观看不下载黄p国产| 欧美日本视频| 长腿黑丝高跟| 成年av动漫网址| 两个人的视频大全免费| 亚洲精品久久国产高清桃花| 亚洲最大成人中文| 色综合站精品国产| 成人美女网站在线观看视频| 久久精品国产亚洲av香蕉五月| 成人性生交大片免费视频hd| 免费无遮挡裸体视频| 变态另类丝袜制服| 观看免费一级毛片| 欧美bdsm另类| 欧美色欧美亚洲另类二区| 在线观看免费视频日本深夜| 观看美女的网站| 国产色爽女视频免费观看| 深夜a级毛片| 亚洲av免费高清在线观看| 日本熟妇午夜| 国产精品精品国产色婷婷| 精品久久久久久久久久久久久| 国产三级中文精品| 久久久国产成人免费| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 国产又黄又爽又无遮挡在线| 黄片wwwwww| 日韩 亚洲 欧美在线| 国产精品乱码一区二三区的特点| 嫩草影院精品99| 日韩欧美国产在线观看| 日产精品乱码卡一卡2卡三| 一本久久中文字幕| 网址你懂的国产日韩在线| 色综合站精品国产| 联通29元200g的流量卡| 香蕉av资源在线| 国产精品一二三区在线看| 99视频精品全部免费 在线| 成人特级av手机在线观看| 夜夜夜夜夜久久久久| 中文字幕久久专区| 国语自产精品视频在线第100页| av在线观看视频网站免费| 狠狠狠狠99中文字幕| 国产久久久一区二区三区| 日韩国内少妇激情av| 午夜影院日韩av| 少妇猛男粗大的猛烈进出视频 | 性色avwww在线观看| 亚洲三级黄色毛片| 日韩精品青青久久久久久| 两个人视频免费观看高清| av黄色大香蕉| 精品熟女少妇av免费看| 国产一区二区激情短视频| 麻豆成人午夜福利视频| 久久人妻av系列| 亚洲av二区三区四区| 99久久中文字幕三级久久日本| 午夜免费激情av| 啦啦啦韩国在线观看视频| 婷婷精品国产亚洲av在线| 毛片女人毛片| h日本视频在线播放| 噜噜噜噜噜久久久久久91| 中文字幕免费在线视频6| 赤兔流量卡办理| 看非洲黑人一级黄片| 免费搜索国产男女视频| 日韩欧美国产在线观看| 久久久久久九九精品二区国产| 亚洲欧美日韩卡通动漫| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线观看免费| 日韩欧美一区二区三区在线观看| 国产精品久久久久久av不卡| 蜜桃亚洲精品一区二区三区| 亚洲丝袜综合中文字幕| 国产精品乱码一区二三区的特点| 国产伦精品一区二区三区视频9| 国产精品日韩av在线免费观看| 少妇人妻精品综合一区二区 | 国产高清视频在线播放一区| 久久国内精品自在自线图片| 日本黄色视频三级网站网址| 一本久久中文字幕| 日本成人三级电影网站| 亚洲aⅴ乱码一区二区在线播放| 午夜日韩欧美国产| 淫秽高清视频在线观看| 久久久久久久久中文| 久久精品影院6| 老女人水多毛片| avwww免费| 欧美中文日本在线观看视频| 国产真实伦视频高清在线观看| 国产黄色视频一区二区在线观看 | 在线观看66精品国产| 国产黄色小视频在线观看| 久久久欧美国产精品| 国产精品一区二区三区四区久久| 精品一区二区三区av网在线观看| 国产私拍福利视频在线观看| 久久久久精品国产欧美久久久| 久久精品久久久久久噜噜老黄 | 非洲黑人性xxxx精品又粗又长| 两个人的视频大全免费| 亚洲欧美中文字幕日韩二区| 国产精品综合久久久久久久免费| 国产精品一区www在线观看| 国产熟女欧美一区二区| 亚洲精品久久国产高清桃花| 国产老妇女一区| av卡一久久| 不卡视频在线观看欧美| av在线播放精品| 亚洲精华国产精华液的使用体验 | 中文字幕精品亚洲无线码一区| 亚洲一区二区三区色噜噜| 高清毛片免费观看视频网站| 欧美区成人在线视频| 亚洲精品影视一区二区三区av| 色5月婷婷丁香| a级毛片免费高清观看在线播放| 亚洲欧美日韩高清在线视频| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 亚洲av美国av| 青春草视频在线免费观看| 国产精品免费一区二区三区在线| 中文字幕人妻熟人妻熟丝袜美| 99久久精品热视频| 又爽又黄a免费视频| 一a级毛片在线观看| 国产单亲对白刺激| 日本色播在线视频| 如何舔出高潮| 久久精品国产自在天天线| 如何舔出高潮| 热99在线观看视频| 别揉我奶头 嗯啊视频| 国产一区二区在线观看日韩| 99热这里只有精品一区| 菩萨蛮人人尽说江南好唐韦庄 | 最后的刺客免费高清国语| 老熟妇乱子伦视频在线观看| 黄片wwwwww| av女优亚洲男人天堂| 麻豆久久精品国产亚洲av| 午夜精品在线福利| 国产亚洲精品久久久久久毛片| 国模一区二区三区四区视频| 中文亚洲av片在线观看爽| 日韩高清综合在线| 亚洲欧美精品自产自拍| 亚洲欧美中文字幕日韩二区| 久久欧美精品欧美久久欧美| 国产69精品久久久久777片| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲网站| 日产精品乱码卡一卡2卡三| 久久久久性生活片| 国产精品福利在线免费观看| 成人鲁丝片一二三区免费| 国产一区二区三区av在线 | 日本熟妇午夜| 国产精品久久久久久av不卡| 在线天堂最新版资源| 男人的好看免费观看在线视频| 九九在线视频观看精品| 中国美女看黄片| 国产视频一区二区在线看| 18禁裸乳无遮挡免费网站照片| 一级黄色大片毛片| 午夜精品在线福利| 大型黄色视频在线免费观看| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 精品日产1卡2卡| 午夜a级毛片| 国内久久婷婷六月综合欲色啪| 男女做爰动态图高潮gif福利片| 内地一区二区视频在线| 美女黄网站色视频| 国产av麻豆久久久久久久| 丰满乱子伦码专区| 久久久久久久久久成人| 99国产精品一区二区蜜桃av| 国产精品1区2区在线观看.| 亚洲精品久久国产高清桃花| 国产一级毛片七仙女欲春2| 久久久国产成人免费| 22中文网久久字幕| 在线看三级毛片| 亚洲自拍偷在线| 女的被弄到高潮叫床怎么办| 亚洲美女视频黄频| 精品一区二区三区av网在线观看| 一级毛片电影观看 | 麻豆av噜噜一区二区三区| 亚洲国产精品sss在线观看| 老司机午夜福利在线观看视频| 亚洲精品成人久久久久久| 精品日产1卡2卡| 男女做爰动态图高潮gif福利片| 日韩欧美国产在线观看| 国产成人福利小说| 亚洲性久久影院| 成人午夜高清在线视频| 高清午夜精品一区二区三区 | 青春草视频在线免费观看| 久久久久久久午夜电影| 国产精品99久久久久久久久| 亚洲婷婷狠狠爱综合网| av.在线天堂| 国产精品一区www在线观看| 欧美+亚洲+日韩+国产| 最新在线观看一区二区三区| 国产高清激情床上av| 最近最新中文字幕大全电影3| 国产一区亚洲一区在线观看| 人人妻人人看人人澡| 我要搜黄色片| 看十八女毛片水多多多| 亚洲av免费高清在线观看| 久久久久久九九精品二区国产| 秋霞在线观看毛片| 97人妻精品一区二区三区麻豆| 内地一区二区视频在线| 亚洲国产欧洲综合997久久,| 久久久久久久久中文| 观看美女的网站| 一个人看的www免费观看视频| 日本一二三区视频观看| 日韩精品青青久久久久久| 久久久国产成人精品二区| 亚洲欧美成人综合另类久久久 | 一卡2卡三卡四卡精品乱码亚洲| www.色视频.com| 成人综合一区亚洲| 欧美一区二区亚洲| 欧美激情久久久久久爽电影| 日韩高清综合在线| 亚洲激情五月婷婷啪啪| 成熟少妇高潮喷水视频| 国产一区亚洲一区在线观看| 国产精品久久电影中文字幕| 日日啪夜夜撸| 秋霞在线观看毛片| 日本黄大片高清| 国产精华一区二区三区| 日本成人三级电影网站| 免费搜索国产男女视频| 亚洲性久久影院| 丰满乱子伦码专区| 国产探花在线观看一区二区| 99久久中文字幕三级久久日本| 精品日产1卡2卡| 成人三级黄色视频| 国产视频内射| 男女啪啪激烈高潮av片| 亚洲国产色片| 在线观看免费视频日本深夜| 国产精品久久久久久av不卡| 成人高潮视频无遮挡免费网站| 亚洲无线在线观看| 啦啦啦观看免费观看视频高清| 久久中文看片网| 成年av动漫网址| 亚洲精品成人久久久久久| 啦啦啦啦在线视频资源| 日韩一本色道免费dvd| 欧美色视频一区免费| 人妻夜夜爽99麻豆av| 麻豆国产av国片精品| 成年女人看的毛片在线观看| 一本久久中文字幕| 国产高清不卡午夜福利| 亚洲精品久久国产高清桃花| 中文亚洲av片在线观看爽| 俺也久久电影网| 插逼视频在线观看| 成人无遮挡网站| 在线观看66精品国产| 欧美在线一区亚洲| 日韩国内少妇激情av| 高清毛片免费观看视频网站| 亚洲精品在线观看二区| 国内久久婷婷六月综合欲色啪| 久久精品影院6| 搞女人的毛片| 国产成人影院久久av| 97碰自拍视频| 国产久久久一区二区三区| 亚洲熟妇熟女久久| 内射极品少妇av片p| 午夜福利在线观看吧| 日韩人妻高清精品专区| 亚洲久久久久久中文字幕| 久久精品影院6| 不卡视频在线观看欧美| 在现免费观看毛片| 亚洲精品在线观看二区| av黄色大香蕉| 我要看日韩黄色一级片| 国产精品亚洲一级av第二区| 热99re8久久精品国产| 久久亚洲精品不卡| 麻豆一二三区av精品| 91久久精品国产一区二区成人| 1000部很黄的大片| 亚洲av不卡在线观看| 国产欧美日韩精品亚洲av| 两性午夜刺激爽爽歪歪视频在线观看| 最近手机中文字幕大全| 色在线成人网| 亚洲在线自拍视频| 一夜夜www| 日本黄色视频三级网站网址| 干丝袜人妻中文字幕| 久久久久久久久久黄片| 欧美绝顶高潮抽搐喷水| 亚洲欧美日韩高清专用| 成人三级黄色视频| 人妻丰满熟妇av一区二区三区| 成人漫画全彩无遮挡| 国产亚洲精品久久久久久毛片| 国产日本99.免费观看| 九九久久精品国产亚洲av麻豆| 少妇的逼水好多| 卡戴珊不雅视频在线播放| 免费在线观看影片大全网站| 午夜福利成人在线免费观看| 成年女人永久免费观看视频| 黄色日韩在线| 日本撒尿小便嘘嘘汇集6| av免费在线看不卡| 日韩成人伦理影院| 日本一二三区视频观看| 美女大奶头视频| 青春草视频在线免费观看| 久久6这里有精品| 精品一区二区三区人妻视频| 老司机午夜福利在线观看视频| 国产亚洲91精品色在线| 色综合色国产| 中文字幕精品亚洲无线码一区| 国产成人91sexporn| 亚洲无线观看免费| 日本一二三区视频观看| 男女那种视频在线观看| 在线观看美女被高潮喷水网站| 亚洲中文字幕一区二区三区有码在线看| 成年女人毛片免费观看观看9| 久久人人爽人人爽人人片va| 老熟妇乱子伦视频在线观看| 午夜激情福利司机影院| 日韩制服骚丝袜av| 91狼人影院| 男人和女人高潮做爰伦理| 黄片wwwwww| 少妇的逼水好多| 俺也久久电影网| 日本色播在线视频| 日日摸夜夜添夜夜添小说| 亚洲高清免费不卡视频| 2021天堂中文幕一二区在线观| 老女人水多毛片| 亚洲av第一区精品v没综合| 老熟妇仑乱视频hdxx| 色尼玛亚洲综合影院| 深爱激情五月婷婷| 麻豆久久精品国产亚洲av| 内地一区二区视频在线| 国产精品精品国产色婷婷| 亚洲欧美精品自产自拍| 国产精品久久电影中文字幕| 嫩草影院精品99| 永久网站在线| 欧美高清性xxxxhd video| 亚洲无线在线观看| 亚洲电影在线观看av| 人妻久久中文字幕网| 久久久久久大精品| 亚洲第一区二区三区不卡| 一卡2卡三卡四卡精品乱码亚洲| a级毛片a级免费在线| 亚洲电影在线观看av| 欧美极品一区二区三区四区| 国产精品日韩av在线免费观看| 亚洲人成网站在线播放欧美日韩| 久久草成人影院| 久久久精品94久久精品| 亚洲国产高清在线一区二区三| 色av中文字幕| 乱人视频在线观看| 1000部很黄的大片| 久久久久久久久久久丰满| 少妇裸体淫交视频免费看高清| 久久久久久九九精品二区国产| 天堂影院成人在线观看| 搡老妇女老女人老熟妇| 欧美+亚洲+日韩+国产| 日本一本二区三区精品| 久久韩国三级中文字幕| 午夜福利18| 天堂av国产一区二区熟女人妻| 亚洲欧美清纯卡通|