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

    內(nèi)嵌管式輻射地板的頻域熱特性分析

    2015-03-08 06:34:45李安邦徐新華
    關(guān)鍵詞:絕熱層相角管式

    李安邦,徐新華

    (華中科技大學(xué) 建筑環(huán)境與能源工程系,湖北 武漢 430074)

    內(nèi)嵌管式輻射地板的頻域熱特性分析

    李安邦,徐新華?

    (華中科技大學(xué) 建筑環(huán)境與能源工程系,湖北 武漢 430074)

    基于經(jīng)典有限差分原理建立了內(nèi)嵌管式輻射地板的頻域有限差分(FDFD)模型,同時(shí)采用Fluent軟件建立了內(nèi)嵌管式輻射地板的CFD模型作為參考模型,將內(nèi)嵌管式輻射地板FDFD模型在典型頻域點(diǎn)的計(jì)算結(jié)果轉(zhuǎn)換成時(shí)域內(nèi)的幅值與相角,并與CFD模型的計(jì)算結(jié)果進(jìn)行對(duì)比.FDFD模型能準(zhǔn)確預(yù)測(cè)內(nèi)嵌管式輻射地板的熱特性.采用FDFD模型進(jìn)一步計(jì)算了不同厚度絕熱層的輻射地板的頻域熱特性,分析了絕熱層厚度對(duì)內(nèi)嵌管式輻射地板熱特性的影響.結(jié)果表明,在高頻區(qū)域,絕熱層厚度對(duì)地板傳熱的影響較小,而在低頻區(qū)域內(nèi)影響較為明顯,盡管絕熱層厚度取到40 mm,地板下表面仍存在較大的熱流損失,約占管道熱流的16%.

    頻域有限差分法;內(nèi)嵌管式輻射地板;頻域熱特性;CFD

    由于內(nèi)嵌管式地板輻射系統(tǒng)相比于空氣空調(diào)系統(tǒng)具有更好的熱舒適性以及更高的能效,該系統(tǒng)在建筑里越來(lái)越多地被使用[1-3].在內(nèi)嵌管式地板輻射系統(tǒng)的設(shè)計(jì)過(guò)程中,準(zhǔn)確的能耗模擬模型可以幫助設(shè)計(jì)人員確定合適的系統(tǒng)規(guī)模以及選擇更優(yōu)化的系統(tǒng)運(yùn)行策略[4].關(guān)于內(nèi)嵌管式地板傳熱模型的動(dòng)態(tài)解析求解很難實(shí)現(xiàn),目前在內(nèi)嵌管式地板傳熱模擬計(jì)算中主要采用的是半解析模型[5]、簡(jiǎn)化計(jì)算模型[6]以及數(shù)值模型[7-8].半解析模型和簡(jiǎn)化模型雖然計(jì)算比較方便,但由于對(duì)原物理模型進(jìn)行了較大程度的簡(jiǎn)化,計(jì)算結(jié)果偏差較大.時(shí)域數(shù)值模型(比如有限元模型、有限差分模型等)都可以獲得具有較高準(zhǔn)確度的計(jì)算結(jié)果,但是需要耗費(fèi)大量的計(jì)算時(shí)間,且精度越高所需的計(jì)算時(shí)間代價(jià)越高.頻域差分法(frequency-domain finite-difference method, FDFDM) 在電磁領(lǐng)域有著較為廣泛的應(yīng)用[9-10],而在求解傳熱問(wèn)題方面的應(yīng)用較少.相比于時(shí)域的數(shù)值模型,F(xiàn)DFDM主要有兩大特點(diǎn):其一,計(jì)算過(guò)程中沒(méi)有時(shí)間項(xiàng),在求解周期性傳熱問(wèn)題時(shí),無(wú)論周期的長(zhǎng)短,F(xiàn)DFDM計(jì)算熱流以及溫度的衰減和延遲所需要的時(shí)間遠(yuǎn)小于時(shí)域數(shù)值模型;其二,F(xiàn)DFDM計(jì)算出頻域熱響應(yīng)可以作為簡(jiǎn)化模型參數(shù)辨識(shí)的數(shù)據(jù)基礎(chǔ),例如:多項(xiàng)式傳遞函數(shù)模型的辨識(shí)[11],RC網(wǎng)絡(luò)模型的辨識(shí)[12].本文將采用FDFD方法計(jì)算帶有不同絕熱層厚度的內(nèi)嵌管式輻射地板的頻域熱特性,并分析不同絕熱層厚度對(duì)內(nèi)嵌管式輻射地板的頻域熱特性的影響.

    1 內(nèi)嵌管式輻射地板的FDFD模型

    內(nèi)嵌管式輻射地板的結(jié)構(gòu)如圖1所示,其中絕熱層的厚度δc將作為變量來(lái)討論絕熱層厚度對(duì)內(nèi)嵌管式輻射地板的頻域熱特性的影響.以圖1中的abcd結(jié)構(gòu)為研究對(duì)象,并將該結(jié)構(gòu)劃分成有限個(gè)微小的矩形單元.內(nèi)嵌管式輻射地板中各材料的熱物性如表1所示.

    將所有溫度表示成復(fù)數(shù)形式, 如式(1)所示:

    (1)

    (2)

    (3)

    圖2以單元0為例展示了該單元與周?chē)噜弳卧g的二維傳熱機(jī)理(沿管道長(zhǎng)度的溫降較小,不考慮沿管道長(zhǎng)度方向的傳熱),據(jù)此可以建立單元0的熱流平衡方程如式(4)所示.

    圖1 內(nèi)嵌管式輻射地板的離散原理

    表1 混凝土和絕熱材料的熱物性參數(shù)

    s1(θ1-θ0)+s2(θ2-θ0)+s3(θ3-θ0)+

    (4)

    式中:ρ為單元0的密度,kg/m3;c為單元0的熱容,J/kg·K;sj(j=1,2,3,4)為相鄰兩單元間的傳熱系數(shù),W/m2·K,可以表示為:

    (5)

    (6)

    圖2 單元j與周?chē)噜弳卧g的二維傳熱機(jī)理

    其中λj表示單元j的導(dǎo)熱系數(shù),W/m·K.將方程(1), (5), (6)代入方程(4)中可得關(guān)于單元0 熱流平衡的2×2矩陣方程為:

    (7)

    聯(lián)立所有單元的矩陣方程得到關(guān)于整個(gè)計(jì)算區(qū)域熱流平衡的大型稀疏矩陣方程,通過(guò)對(duì)這一大型稀疏矩陣方程進(jìn)行求解就可以得到任意單元頻域熱響應(yīng).

    2 內(nèi)嵌管式輻射地板的CFD模型

    隨著計(jì)算機(jī)科學(xué)和數(shù)值計(jì)算等學(xué)科的發(fā)展,計(jì)算機(jī)數(shù)值模擬日趨成熟,在一些領(lǐng)域內(nèi)實(shí)驗(yàn)成本很高甚至有的實(shí)驗(yàn)難以實(shí)現(xiàn),許多研究[13-14]常采用數(shù)值模擬替代真實(shí)實(shí)驗(yàn)作為檢驗(yàn)?zāi)P褪欠駵?zhǔn)確的參考依據(jù).本文采用該結(jié)構(gòu)CFD模型預(yù)測(cè)結(jié)果作為參考值.建立圖1中abcd結(jié)構(gòu)的CFD模型,通過(guò)CFD模型和FDFD模型計(jì)算數(shù)據(jù)的對(duì)比來(lái)說(shuō)明FDFD模型的準(zhǔn)確性.對(duì)于任意頻率點(diǎn)ω,CFD模型的外擾為以2π/ω為周期的單位諧波,該單位諧波將以用戶(hù)自定義函數(shù)(UDF)的形式輸入到CFD模擬軟件中.在FDFD模型和CFD模型中,abcd結(jié)構(gòu)的上表面ab定義為第3類(lèi)邊界條件,ab表面與室內(nèi)的綜合換熱系數(shù)為7W/m2·K,下表面cd和管道表面(忽略管道壁厚)定義為第1類(lèi)溫度邊界條件(給定溫度),表面ac和bd均看作絕熱面.

    3 FDFD模型與CFD模型的計(jì)算結(jié)果對(duì)比

    對(duì)于線性導(dǎo)熱問(wèn)題, 頻域熱響應(yīng)(板壁對(duì)溫度諧波熱擾的衰減和時(shí)間延遲或者相角差)與溫度諧波熱擾的幅值和相角均無(wú)關(guān)[15].內(nèi)嵌管式輻射地板的導(dǎo)熱微分方程為線性方程,即同樣為線性導(dǎo)熱問(wèn)題, 地板各表面的頻域熱流響應(yīng)的幅值與溫度諧波熱擾的幅值呈線性比例關(guān)系(該比例只與頻率和該結(jié)構(gòu)的熱物性有關(guān)), 兩者的相角差與溫度諧波熱擾的相角無(wú)關(guān).因此關(guān)于溫度諧波熱擾的幅值和相角的選取對(duì)內(nèi)嵌管式輻射地板頻域熱特性的分析并無(wú)影響.

    以不帶絕熱層的內(nèi)嵌管式輻射地板為例,且考慮單位溫度諧波設(shè)置在地板上表面的情況(其余兩表面的溫度均設(shè)為0),由FDFD模型以及CFD模型計(jì)算出的地板上表面和管道表面的周期分別為12h和48h的熱流曲線分別如圖3及圖4所示.通過(guò)比較可以發(fā)現(xiàn),F(xiàn)DFD模型和CFD模型的計(jì)算結(jié)果吻合很好,這說(shuō)明FDFD模型能很好地計(jì)算內(nèi)嵌管式輻射地板的熱特性.計(jì)算時(shí)間表明FDFD模型的計(jì)算時(shí)間遠(yuǎn)遠(yuǎn)比CFD模型少,計(jì)算快速高效.

    時(shí)間/h

    時(shí)間/h

    4 內(nèi)嵌管式輻射地板的頻域熱響應(yīng)

    頻域熱特性包括幅頻特性和相頻特性,幅頻特性是指溫度諧波熱擾傳遞的衰減倍數(shù)(輸入溫度諧波的幅值與輸出溫度諧波的幅值之間的比值)與頻率的關(guān)系,相頻特性是指輸出溫度諧波和輸入的溫度諧波之間的相位差或者相角差與頻率的關(guān)系.衰減倍數(shù)大和相位差大均可以反映系統(tǒng)具有較好的隔熱效果.為了研究絕熱層厚度對(duì)內(nèi)嵌管式輻射地板的頻域熱特性的影響,本文建立了不同厚度絕熱層的內(nèi)嵌管式地板的FDFD模型,并計(jì)算出10-10~10-3rad/s之間的一系列離散頻率點(diǎn)上的熱響應(yīng)(各表面熱流的幅值和相角).在工程實(shí)際中常用的絕熱層厚度分別為20,30,40mm, 為了進(jìn)行對(duì)比研究,本文還計(jì)算分析了絕熱層厚度分別為100,200mm及理想絕熱情況下內(nèi)嵌管式輻射地板的頻域熱特性.圖5~圖7為管道邊界存在幅值為1 ℃,相角為0的溫度諧波作用而其余兩表面(地板上表面和地板下表面)的溫度均為0時(shí),絕熱層厚度對(duì)地板各表面的熱流響應(yīng)的影響.

    頻率/(rad·s-1) (a)幅值

    頻率/(rad·s-1) (b)相角

    頻率/(rad·s-1) (a)幅值

    頻率/(rad·s-1) (b)相角

    頻率/(rad·s-1) (a)幅值

    頻率/(rad·s-1) (b)相角

    由圖5(a)可知,在低頻區(qū)域,隨著絕熱層厚度由20mm增加到200mm,地板上表面熱流熱流幅值平均增加了4%,在高頻區(qū)域熱流幅值幾乎不受絕熱層厚度的影響.圖5(b)顯示絕熱層厚度對(duì)地板上表面熱流相角幾乎沒(méi)有的影響.圖5顯示的結(jié)果說(shuō)明,內(nèi)嵌管式地板絕熱層厚度對(duì)地板上表面與房間之間的換熱(內(nèi)嵌管式地板輻射系統(tǒng)的制冷或者制熱量)沒(méi)有明顯的影響.

    由圖6(a)可以發(fā)現(xiàn),在低頻區(qū)域地板管道表面的熱流幅值隨著絕熱層厚度的增加會(huì)有明顯的減少(絕熱層厚度由20mm增加到200mm,管道表面熱流幅值減少了20%),而在高頻區(qū)域這一變化并不明顯.這主要是由于在低頻區(qū)域,地板主要表現(xiàn)的是穩(wěn)態(tài)特性(此時(shí)地板熱質(zhì)的蓄熱較小),隨著絕熱層厚度的增加,從管道到地板下表面的傳熱熱阻也增加了,導(dǎo)致地板下表面的冷量或熱量損失越來(lái)越小,而在高頻區(qū)域地板表現(xiàn)的是動(dòng)態(tài)特性,但絕熱層的熱質(zhì)較小,對(duì)該動(dòng)態(tài)特性影響很小.因此在高頻區(qū)域絕熱層厚度對(duì)管道表面熱流幅值的影響很小.由圖5(b)可知,絕熱層厚度對(duì)管道表面熱流的相角幾乎沒(méi)有影響,這是由于輕質(zhì)絕熱層材料的熱質(zhì)很小,對(duì)管道表面熱流幾乎沒(méi)有延遲作用.

    由圖7(a)可知,在低頻區(qū)域絕熱層厚度由20mm增加到200mm,地板下表面的熱流幅值由1.68W/m2·K減少為0.18W/m2·K,變化十分明顯,高頻區(qū)域的變化則并不明顯.另外,絕熱層厚度分別為20,30,40,100和200mm時(shí),地板下表面熱流占管道表面熱流的百分比分別為26%,19%,16%,7%和4%.說(shuō)明絕熱層厚度分別為30mm和40mm時(shí),地板下表面仍然存在明顯的冷量或熱量的損失,進(jìn)一步說(shuō)明絕熱層厚度對(duì)內(nèi)嵌管式地板輻射系統(tǒng)熱損失或冷損失不可忽視.

    5 總 結(jié)

    本文建立了內(nèi)嵌管式輻射地板的FDFD模型,并計(jì)算了該輻射地板的頻域熱響應(yīng).建立了內(nèi)嵌管式輻射地板的CFD模型并作為參考模型,兩個(gè)模型的計(jì)算結(jié)果對(duì)比表明FDFD模型能準(zhǔn)確計(jì)算頻域熱特性.通過(guò)比較兩個(gè)模型的計(jì)算時(shí)間表明,在求解周期性傳熱問(wèn)題時(shí),F(xiàn)DFD模型的計(jì)算時(shí)間遠(yuǎn)比CFD模型少、且更加方便.采用FDFD模型計(jì)算了不同絕熱層厚度的內(nèi)嵌管式輻射地板的頻域熱響應(yīng),并分析了絕熱層厚度對(duì)地板各表面的熱流幅值及相角的影響.雖然在實(shí)際工程中絕熱層的厚度取到了20,30及40mm, 地板輻射系統(tǒng)的冷損失或熱損失依然很大.

    [1]LEHAMANNB,DORERV,KOSCHENZM.Applicationrangeofthermallyactivatedbuildingsystemstabs[J].EnergyandBuilding, 2007, 39 (5): 593-598.

    [2]SIMMONDSP,HOLSTS,REUSSS,etal.Using radiant cooled floors to condition large spaces and maintain comfort conditions[J]. ASHRAE Transaction, 2000,106 (1) : 695-701.

    [3] FEUSTEL H E, STETIU C. Hydronic radiant cooling: preliminary assessment[J]. Energy and Buildings, 1995, 22(3): 193-205.

    [4] LAOUADI A. Development of a radiant heating and cooling model for building energy simulation software[J]. Energy and Buildings, 2004, 39 (4): 421-431.

    [5] ZHANG Z,PATE M B.A semi-analytical formulation for heat transfer from structures with embedded tubes[C]//The 24th National Heat Transfer Conference.Pittsburgh, Pennsylvania:ASME-HTD,1987,78:17-25.

    [6] KOSCHENZ M,DORER V.Interaction of an air system with concrete core conditioning[J]. Building and Environment, 1999, 30(2): 139-145.

    [7] ANTONOPOULOS K A,TZIVANIDIS C.Numerical solution of unsteady three-dimensional heat transfer during space cooling using ceiling embedded piping[J]. Energy,1997, 22(1): 59-67.

    [8] HO S Y,HAYES R E,WOOD R K.Simulation of the dynamic behaviour of a hydronic floor heating system[J]. Heat Recovery Systems and CHP,1995,15(6): 505-519.

    [9] ALBANI M,BERNARDI P.A numerical method based on the discretization of maxwell equations in integral form[J]. IEEE Transactions on Microwave Theory and Techniques, 1974, 22(4): 446-450.

    [10]LAVRANOS C S,KYRIACOU G A,SAHALOS J N.A 2-D finite difference frequency domain (FDFD) eigenvalue method for orthogonal curvilinear coordinates [C]//Progress in Electromagnetic Research Symposium.Pisa, Italy:The Electromagnetics Academy, 2004:398-400.

    [11]WANG Sheng-wei, CHEN You-ming. Transient heat flow calculation for multilayer constructions using a frequency-domain regression method[J]. Building and Environment, 2003, 38(1): 45-61.

    [12]WANG Sheng-wei, XU Xin-hua. Simplified building model for transient thermal performance estimation using GA-based parameter identification[J]. International Journal of Thermal Sciences, 2006,45:419-432.

    [13]李偉光,李安邦,徐新華,等.復(fù)雜船舶圍壁傳熱系數(shù)取值探討[J],中國(guó)艦船研究,2014, 9(2): 78-83.

    LI Wei-guang,LI An-bang,XU Xin-hua,etal. Heat transfer coefficient of complex ship casing walls[J]. Chinese Journal of Ship Research, 2014, 9(2): 78-83. (In Chinese)

    [14]LIU Kui-xing, TIAN Zhe, ZHANG Cheng,etal. Establishment and validation of modified star-type RC-network model for concrete core cooling slab[J]. Energy and Buildings, 2011,43( 9): 2378-2384.

    [15]彥啟森, 趙慶珠. 建筑熱過(guò)程[M]. 北京: 中國(guó)建筑工業(yè)出版社, 1996: 57-63.

    YAN Qi-seng,ZHAO Qing-zhu.Thernal process in buildings[M]. Beijing: China Architecture & Building Press, 1996: 57-63. (In Chinese)

    Study of the Frequency Thermal Characteristics of Pipe-embedded Radiant Floors

    LI An-bang, XU Xin-hua?

    (Dept of Building Environment and Energy Engineering, Huazhong Univ of Science and Technology, Wuhan,Hubei 430074, China)

    A two-dimensional frequency-domain finite-difference (FDFD) model was developed on the basis of the classical finite-difference principle for a pipe-embedded concrete radiant floor to predict its frequency thermal responses. In the meantime, a CFD for this floor was also developed for this floor by using a reference model. The calculated frequency thermal responses of the floor by using FDFD model at typical frequency points were expressed in the form of amplitude and phase lag in time domain, and were then compared with the results by using CFD model. The results indicate that the FDFD model can predict the thermal characteristics accurately. The FDFD model was used to predict the frequency thermal characteristics of the pipe-embedded concrete radiant floor with insulation of different thicknesses, and the effects of the insulation thickness on the thermal responses of the pipe-embedded concrete radiant floor were analyzed and presented. The results indicate that the effect of the insulation thickness on the heat transfer of the floor is obvious in low frequency range while not obvious in high frequency range. Even if the thickness of the insulation is 40mm, the heat loss from the bottom surface of the floor constitutes about 16% of the heat transfer on the pipe surface and cannot be neglected.

    FDFD (frequency-domain finite-difference) model; pipe-embedded concrete radiant floor; frequency thermal characteristics; Computational Fluid Dynamics

    1674-2974(2015)01-0115-05

    2014-04-15

    國(guó)家自然科學(xué)基金資助項(xiàng)目 (51178201), National Natural Science Foundation of China(51178201);新世紀(jì)優(yōu)秀人才支持計(jì)劃資助項(xiàng)目(NCET110189); 高等學(xué)校博士學(xué)科點(diǎn)專(zhuān)項(xiàng)科研基金資助項(xiàng)目(20120142110078)

    李安邦(1989-),男,湖北武穴人,華中科技大學(xué)博士研究生?通訊聯(lián)系人,E-mail:bexhxu@mail.hust.edu.cn

    TU832

    A

    猜你喜歡
    絕熱層相角管式
    管式太陽(yáng)能集熱器的顆粒換熱模擬
    固體火箭發(fā)動(dòng)機(jī)EPDM絕熱層產(chǎn)品一體化工作實(shí)踐
    基于Dynamo的熱力管道絕熱層經(jīng)濟(jì)厚度計(jì)算
    甘肅科技(2022年23期)2022-02-07 12:17:48
    基于實(shí)時(shí)服役參數(shù)的交流接觸器電壽命最大化控制策略
    固體發(fā)動(dòng)機(jī)飛行橫向過(guò)載下絕熱層燒蝕探究①
    配電網(wǎng)30°相角差線路不停電轉(zhuǎn)供方案探究
    電子制作(2018年1期)2018-04-04 01:48:40
    基于相角差的絕緣子表面污穢受潮檢測(cè)系統(tǒng)設(shè)計(jì)與仿真
    管式空氣空預(yù)器泄漏分析及改進(jìn)
    小電流接地故障初相角模擬試驗(yàn)控制系統(tǒng)的設(shè)計(jì)
    三元乙丙橡膠絕熱層的燒蝕特性研究
    航空兵器(2015年2期)2015-11-15 05:12:40
    26uuu在线亚洲综合色| 欧美丝袜亚洲另类| 少妇的逼水好多| 国产av码专区亚洲av| 男人舔女人下体高潮全视频| 搡老妇女老女人老熟妇| 91精品国产九色| 18+在线观看网站| 欧美日韩一区二区视频在线观看视频在线 | 伦精品一区二区三区| 久久精品人妻少妇| 亚洲国产欧美人成| 97精品久久久久久久久久精品| 日本三级黄在线观看| 亚洲欧美日韩东京热| 2021天堂中文幕一二区在线观| 亚洲国产欧美在线一区| 成人性生交大片免费视频hd| 美女被艹到高潮喷水动态| 天堂中文最新版在线下载 | 日韩av免费高清视频| 高清欧美精品videossex| 久久精品国产鲁丝片午夜精品| 天天一区二区日本电影三级| 蜜臀久久99精品久久宅男| 亚洲内射少妇av| 2021少妇久久久久久久久久久| 熟妇人妻不卡中文字幕| 国产男女超爽视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 日本爱情动作片www.在线观看| 非洲黑人性xxxx精品又粗又长| 国产欧美日韩精品一区二区| 亚洲av中文av极速乱| 亚洲欧洲日产国产| 国产黄色视频一区二区在线观看| 国产又色又爽无遮挡免| 天天躁夜夜躁狠狠久久av| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 中文字幕久久专区| 成年女人看的毛片在线观看| 亚洲国产欧美在线一区| 狂野欧美白嫩少妇大欣赏| 欧美潮喷喷水| 免费高清在线观看视频在线观看| av在线天堂中文字幕| 亚洲精品第二区| 国产午夜精品一二区理论片| 免费看a级黄色片| 国产高清国产精品国产三级 | 男人和女人高潮做爰伦理| 真实男女啪啪啪动态图| 国内精品美女久久久久久| 97超视频在线观看视频| 亚洲国产色片| 成人一区二区视频在线观看| 22中文网久久字幕| 免费大片18禁| 亚洲久久久久久中文字幕| 精品酒店卫生间| 成人国产麻豆网| 美女内射精品一级片tv| 亚洲av在线观看美女高潮| 久久久久久久久久久免费av| 最近最新中文字幕免费大全7| 97热精品久久久久久| 尤物成人国产欧美一区二区三区| 大话2 男鬼变身卡| 亚洲精品国产成人久久av| 国产黄色视频一区二区在线观看| 麻豆成人av视频| 日韩欧美精品v在线| 国产黄片视频在线免费观看| 亚洲久久久久久中文字幕| 波野结衣二区三区在线| 国产成年人精品一区二区| 久久国内精品自在自线图片| 久久精品熟女亚洲av麻豆精品 | 好男人视频免费观看在线| 日本三级黄在线观看| 国产av码专区亚洲av| 97超视频在线观看视频| 亚洲人与动物交配视频| 久久鲁丝午夜福利片| 国产不卡一卡二| 男女视频在线观看网站免费| 天堂网av新在线| 亚洲av日韩在线播放| 麻豆国产97在线/欧美| 国产人妻一区二区三区在| 国产欧美日韩精品一区二区| av免费观看日本| 一个人免费在线观看电影| 国产精品蜜桃在线观看| 欧美日韩一区二区视频在线观看视频在线 | 哪个播放器可以免费观看大片| 一个人免费在线观看电影| 久久久久久久大尺度免费视频| 国产成人精品福利久久| 国产欧美另类精品又又久久亚洲欧美| 非洲黑人性xxxx精品又粗又长| 亚洲无线观看免费| 国内揄拍国产精品人妻在线| 男人和女人高潮做爰伦理| 秋霞伦理黄片| 久久久亚洲精品成人影院| 26uuu在线亚洲综合色| 免费观看无遮挡的男女| 韩国高清视频一区二区三区| 日韩一区二区视频免费看| 亚洲av中文av极速乱| 午夜亚洲福利在线播放| 天堂影院成人在线观看| 亚洲精品自拍成人| 男女视频在线观看网站免费| 日本一本二区三区精品| 自拍偷自拍亚洲精品老妇| 日本av手机在线免费观看| 女的被弄到高潮叫床怎么办| .国产精品久久| 日韩三级伦理在线观看| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜 | 91精品国产九色| 国产精品1区2区在线观看.| 日韩欧美三级三区| 午夜激情欧美在线| 国产精品一区二区在线观看99 | 狠狠精品人妻久久久久久综合| 欧美最新免费一区二区三区| 色综合站精品国产| 欧美日韩亚洲高清精品| 成人毛片60女人毛片免费| 狂野欧美激情性xxxx在线观看| 国产亚洲一区二区精品| 国内精品一区二区在线观看| 777米奇影视久久| 精品人妻熟女av久视频| 国产一区二区三区综合在线观看 | 国产成人freesex在线| 国内揄拍国产精品人妻在线| 在线免费十八禁| 看黄色毛片网站| 91av网一区二区| 一级毛片我不卡| 免费观看精品视频网站| 精品久久久久久久久亚洲| 97热精品久久久久久| 国产淫语在线视频| 3wmmmm亚洲av在线观看| 婷婷色综合www| 亚洲人成网站在线观看播放| 黄片wwwwww| 午夜免费激情av| 亚洲av中文字字幕乱码综合| 九九久久精品国产亚洲av麻豆| 亚洲四区av| 免费看不卡的av| 97人妻精品一区二区三区麻豆| 18禁在线无遮挡免费观看视频| 国产日韩欧美在线精品| 免费观看a级毛片全部| 亚洲国产最新在线播放| 啦啦啦韩国在线观看视频| av在线观看视频网站免费| 国产免费视频播放在线视频 | 欧美性感艳星| 国产日韩欧美在线精品| 看非洲黑人一级黄片| 亚洲va在线va天堂va国产| 秋霞伦理黄片| 久久这里只有精品中国| 韩国av在线不卡| 国产黄片美女视频| 听说在线观看完整版免费高清| 亚洲久久久久久中文字幕| 亚洲av日韩在线播放| 青青草视频在线视频观看| 日本三级黄在线观看| 国内少妇人妻偷人精品xxx网站| 视频中文字幕在线观看| 精品久久久噜噜| 2021天堂中文幕一二区在线观| 大片免费播放器 马上看| 青春草视频在线免费观看| 亚洲国产精品sss在线观看| 一区二区三区乱码不卡18| 久久这里有精品视频免费| 国产 一区精品| 久久国内精品自在自线图片| 大香蕉久久网| 男女国产视频网站| 精品人妻一区二区三区麻豆| 91aial.com中文字幕在线观看| 久久午夜福利片| 久久这里有精品视频免费| 久久精品国产自在天天线| 大香蕉97超碰在线| 国产精品.久久久| 高清av免费在线| 欧美日韩亚洲高清精品| 国产免费一级a男人的天堂| 国产黄色小视频在线观看| 国产高清三级在线| 日韩伦理黄色片| 在线 av 中文字幕| 一个人看的www免费观看视频| 亚洲国产欧美人成| 国产高清三级在线| 热99在线观看视频| 国产麻豆成人av免费视频| h日本视频在线播放| 麻豆国产97在线/欧美| 国产91av在线免费观看| 2021天堂中文幕一二区在线观| 一本久久精品| 亚洲色图av天堂| 一级片'在线观看视频| 国产一级毛片在线| 久久久午夜欧美精品| 欧美xxxx黑人xx丫x性爽| 人人妻人人看人人澡| 婷婷六月久久综合丁香| 人妻系列 视频| 夫妻性生交免费视频一级片| 中国国产av一级| 国产综合懂色| 又大又黄又爽视频免费| 狂野欧美白嫩少妇大欣赏| 国产男女超爽视频在线观看| 一区二区三区四区激情视频| 免费看不卡的av| 久久精品久久久久久噜噜老黄| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看 | av免费在线看不卡| 午夜精品国产一区二区电影 | 国产三级在线视频| 亚洲国产精品成人久久小说| 国产精品久久久久久久电影| 国产黄片视频在线免费观看| 日日摸夜夜添夜夜添av毛片| 床上黄色一级片| 精品久久久久久成人av| 少妇人妻一区二区三区视频| 午夜激情福利司机影院| 1000部很黄的大片| 99久久精品国产国产毛片| 一个人看视频在线观看www免费| 亚洲av福利一区| 中文字幕制服av| 国产永久视频网站| 成人无遮挡网站| 亚洲性久久影院| 2018国产大陆天天弄谢| 亚洲不卡免费看| 国产极品天堂在线| 中文乱码字字幕精品一区二区三区 | 97在线视频观看| 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕 | 看十八女毛片水多多多| 免费看美女性在线毛片视频| 99久久精品一区二区三区| 精品一区二区免费观看| 特大巨黑吊av在线直播| 亚洲自拍偷在线| av女优亚洲男人天堂| 人人妻人人澡人人爽人人夜夜 | 成人午夜精彩视频在线观看| 国产男人的电影天堂91| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 女的被弄到高潮叫床怎么办| 亚洲熟女精品中文字幕| 久久精品国产亚洲网站| 男女边摸边吃奶| 国语对白做爰xxxⅹ性视频网站| 少妇熟女aⅴ在线视频| 国产一区二区亚洲精品在线观看| 久久精品国产鲁丝片午夜精品| 少妇人妻一区二区三区视频| 色视频www国产| 永久免费av网站大全| 大片免费播放器 马上看| 国产精品99久久久久久久久| 欧美性猛交╳xxx乱大交人| 国产免费又黄又爽又色| 国产精品国产三级国产av玫瑰| 国产v大片淫在线免费观看| 边亲边吃奶的免费视频| 春色校园在线视频观看| 人体艺术视频欧美日本| 国产高清有码在线观看视频| 免费看光身美女| 国产精品三级大全| 亚洲三级黄色毛片| 又粗又硬又长又爽又黄的视频| 99热这里只有精品一区| 黄片wwwwww| 日本黄色片子视频| 久久99精品国语久久久| 久久久久性生活片| 神马国产精品三级电影在线观看| 亚洲精品日本国产第一区| 99视频精品全部免费 在线| 91在线精品国自产拍蜜月| 国产亚洲精品av在线| 国产精品三级大全| 精品久久久噜噜| 亚洲av男天堂| 精品久久久久久久久久久久久| 日本色播在线视频| 成人二区视频| 精华霜和精华液先用哪个| 国产视频内射| 国产午夜精品久久久久久一区二区三区| 水蜜桃什么品种好| 男女啪啪激烈高潮av片| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 欧美日韩国产mv在线观看视频 | 亚洲精品一二三| 国产一区二区亚洲精品在线观看| 中文字幕久久专区| 九色成人免费人妻av| 最近最新中文字幕大全电影3| 婷婷六月久久综合丁香| 又爽又黄a免费视频| 亚洲欧美成人综合另类久久久| 婷婷色av中文字幕| 中文欧美无线码| 麻豆乱淫一区二区| 亚洲真实伦在线观看| 午夜激情福利司机影院| 国产av在哪里看| 亚洲精品第二区| 久久精品人妻少妇| 久久久久性生活片| 亚洲自拍偷在线| 国产亚洲一区二区精品| 亚洲成人精品中文字幕电影| 啦啦啦啦在线视频资源| 97人妻精品一区二区三区麻豆| 免费看av在线观看网站| av国产免费在线观看| 超碰97精品在线观看| 婷婷色麻豆天堂久久| 久久久久精品性色| 亚洲激情五月婷婷啪啪| 亚洲精品第二区| 日韩av不卡免费在线播放| 人妻少妇偷人精品九色| 亚洲精品自拍成人| 久久99热6这里只有精品| 亚洲精品久久久久久婷婷小说| 最近的中文字幕免费完整| 国内少妇人妻偷人精品xxx网站| 麻豆成人午夜福利视频| 男女边摸边吃奶| 午夜激情福利司机影院| 久久99蜜桃精品久久| 久久99热6这里只有精品| 亚洲真实伦在线观看| 免费在线观看成人毛片| 国产三级在线视频| 欧美成人精品欧美一级黄| 少妇人妻一区二区三区视频| 亚洲av电影不卡..在线观看| 亚洲国产av新网站| 极品少妇高潮喷水抽搐| 免费不卡的大黄色大毛片视频在线观看 | 免费高清在线观看视频在线观看| 日本黄色片子视频| 国产一区二区在线观看日韩| av线在线观看网站| 亚洲综合色惰| 国产高清三级在线| 精品久久久久久久久av| 在线播放无遮挡| 中国国产av一级| av在线亚洲专区| 观看美女的网站| 特级一级黄色大片| 一二三四中文在线观看免费高清| 亚洲成色77777| 亚洲精品影视一区二区三区av| 免费无遮挡裸体视频| 日日摸夜夜添夜夜爱| 欧美不卡视频在线免费观看| av卡一久久| 三级男女做爰猛烈吃奶摸视频| av国产久精品久网站免费入址| 国产高清不卡午夜福利| 久久久色成人| 亚洲精品自拍成人| 国产单亲对白刺激| 男女下面进入的视频免费午夜| 狠狠精品人妻久久久久久综合| 热99在线观看视频| 99热这里只有是精品50| 床上黄色一级片| 国产午夜精品一二区理论片| 午夜亚洲福利在线播放| 岛国毛片在线播放| videossex国产| 欧美日韩一区二区视频在线观看视频在线 | 嫩草影院精品99| 国产伦理片在线播放av一区| 欧美激情在线99| 三级经典国产精品| 一本久久精品| 久久久久久九九精品二区国产| 国产美女午夜福利| 国产三级在线视频| 亚洲美女搞黄在线观看| 我的女老师完整版在线观看| 日韩电影二区| 久久久国产一区二区| 国产在线男女| 综合色丁香网| 国产白丝娇喘喷水9色精品| 欧美成人a在线观看| 熟女电影av网| 最近2019中文字幕mv第一页| 日本与韩国留学比较| 久热久热在线精品观看| 久久久久久久久大av| 国产高清三级在线| 男女国产视频网站| 国内精品一区二区在线观看| 一本久久精品| 中文字幕av成人在线电影| 国产亚洲一区二区精品| 亚洲av成人精品一二三区| 老女人水多毛片| 高清av免费在线| 97超碰精品成人国产| 观看美女的网站| 国产白丝娇喘喷水9色精品| 婷婷色综合www| 中文天堂在线官网| 成人二区视频| 国产精品一及| 国产欧美日韩精品一区二区| 午夜亚洲福利在线播放| 大陆偷拍与自拍| 国产老妇伦熟女老妇高清| 亚洲av国产av综合av卡| 亚洲精品国产av蜜桃| 美女黄网站色视频| 波多野结衣巨乳人妻| 欧美一区二区亚洲| 免费av不卡在线播放| 亚洲av免费高清在线观看| 精品人妻偷拍中文字幕| a级毛色黄片| 深夜a级毛片| 久久精品国产亚洲av天美| 亚洲国产欧美在线一区| 尾随美女入室| 午夜福利视频精品| 国产午夜福利久久久久久| 最近2019中文字幕mv第一页| av国产免费在线观看| 亚洲天堂国产精品一区在线| 久久久久久国产a免费观看| 亚洲在久久综合| videossex国产| 中文天堂在线官网| 极品教师在线视频| 爱豆传媒免费全集在线观看| 波多野结衣巨乳人妻| 国产黄片视频在线免费观看| 肉色欧美久久久久久久蜜桃 | 国产精品人妻久久久久久| av在线播放精品| 国产精品久久久久久久电影| 久久精品国产亚洲网站| 亚洲怡红院男人天堂| 国产淫语在线视频| 久久精品久久精品一区二区三区| 欧美区成人在线视频| 国产精品久久久久久精品电影小说 | 欧美高清成人免费视频www| 久久精品久久精品一区二区三区| 免费观看在线日韩| 国产视频内射| av免费在线看不卡| 亚洲精品一二三| 天天躁日日操中文字幕| 91久久精品电影网| 如何舔出高潮| 中文精品一卡2卡3卡4更新| 亚洲精华国产精华液的使用体验| 大片免费播放器 马上看| 欧美精品一区二区大全| 欧美成人精品欧美一级黄| 亚洲精品色激情综合| 中文字幕av在线有码专区| 欧美zozozo另类| 寂寞人妻少妇视频99o| 国产69精品久久久久777片| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片aaaaaa免费看小| 在线观看美女被高潮喷水网站| av天堂中文字幕网| 久久精品夜色国产| 亚洲国产av新网站| 少妇的逼水好多| 毛片一级片免费看久久久久| 日韩人妻高清精品专区| 内地一区二区视频在线| 天堂√8在线中文| 高清视频免费观看一区二区 | 久久久午夜欧美精品| 日本三级黄在线观看| 国国产精品蜜臀av免费| 免费无遮挡裸体视频| 成人性生交大片免费视频hd| 最新中文字幕久久久久| 99久久精品一区二区三区| 一二三四中文在线观看免费高清| 中文天堂在线官网| 亚洲久久久久久中文字幕| 一级a做视频免费观看| 大香蕉久久网| 熟女人妻精品中文字幕| 国内揄拍国产精品人妻在线| 成人午夜精彩视频在线观看| 18禁在线无遮挡免费观看视频| 日韩一本色道免费dvd| 亚洲不卡免费看| 成人特级av手机在线观看| 黄色日韩在线| 51国产日韩欧美| 久久这里只有精品中国| 18+在线观看网站| 成人无遮挡网站| 久久99精品国语久久久| 国产高清三级在线| 日本av手机在线免费观看| 寂寞人妻少妇视频99o| 嘟嘟电影网在线观看| 五月玫瑰六月丁香| 久久鲁丝午夜福利片| 亚洲欧美中文字幕日韩二区| 男女下面进入的视频免费午夜| 秋霞伦理黄片| 精品久久久久久久末码| 亚洲真实伦在线观看| 观看免费一级毛片| 校园人妻丝袜中文字幕| av在线蜜桃| 国产91av在线免费观看| 激情五月婷婷亚洲| 最近中文字幕2019免费版| 韩国av在线不卡| 国产 亚洲一区二区三区 | 久久久欧美国产精品| 欧美日韩国产mv在线观看视频 | ponron亚洲| 精品99又大又爽又粗少妇毛片| 国产精品日韩av在线免费观看| 三级经典国产精品| 久久久久性生活片| 国产欧美另类精品又又久久亚洲欧美| 天天一区二区日本电影三级| 久99久视频精品免费| 精品久久久久久电影网| 亚洲精品一二三| 高清欧美精品videossex| 91av网一区二区| 久久这里有精品视频免费| av国产久精品久网站免费入址| 国产三级在线视频| 成人综合一区亚洲| 性插视频无遮挡在线免费观看| 午夜日本视频在线| av在线老鸭窝| 亚洲欧美日韩卡通动漫| 亚洲精品自拍成人| 国产色爽女视频免费观看| 日日干狠狠操夜夜爽| 中文字幕av成人在线电影| 国产在视频线在精品| 日本免费a在线| 我的老师免费观看完整版| 免费黄网站久久成人精品| 嫩草影院精品99| 国产一级毛片七仙女欲春2| 亚洲人成网站在线播| 乱人视频在线观看| 国产一级毛片七仙女欲春2| 狠狠精品人妻久久久久久综合| 国产视频首页在线观看| 亚洲电影在线观看av| 成人亚洲精品av一区二区| 亚洲最大成人av| 国产午夜福利久久久久久| 亚洲人成网站在线播| 青春草亚洲视频在线观看| 寂寞人妻少妇视频99o| 欧美高清性xxxxhd video| 国产高清国产精品国产三级 | 色网站视频免费| 国产精品伦人一区二区| 久久精品综合一区二区三区| 成人亚洲欧美一区二区av| 久久草成人影院| 国产亚洲最大av| 国产成人精品久久久久久| 中文字幕免费在线视频6| 欧美性感艳星|