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

    圍護結(jié)構(gòu)熱濕耦合傳遞模型及簡便求解方法

    2016-12-07 06:58:51劉向偉陳友明陳國杰郭興國羅娜
    土木與環(huán)境工程學報 2016年4期
    關(guān)鍵詞:含濕量液態(tài)水水蒸氣

    劉向偉,陳友明,陳國杰,3,郭興國,羅娜

    (1. 南昌大學 建筑工程學院,南昌 330031;2. 湖南大學 土木工程學院,長沙 410082;3. 南華大學 城市建設學院,湖南 衡陽 421001)

    ?

    圍護結(jié)構(gòu)熱濕耦合傳遞模型及簡便求解方法

    劉向偉1,陳友明2,陳國杰2,3,郭興國1,羅娜1

    (1. 南昌大學 建筑工程學院,南昌 330031;2. 湖南大學 土木工程學院,長沙 410082;3. 南華大學 城市建設學院,湖南 衡陽 421001)

    為了預測圍護結(jié)構(gòu)內(nèi)的溫度和濕度分布,以連續(xù)變量,相對濕度和溫度為驅(qū)動勢,考慮熱傳遞與濕傳遞之間的耦合作用,建立了圍護結(jié)構(gòu)熱濕耦合傳遞非穩(wěn)態(tài)模型,并提出了基于多物理場耦合仿真模擬軟件COMSOL的熱濕耦合傳遞模型簡便求解方法。通過對比新建模型模擬結(jié)果與HAMSTAD標準驗證實例,驗證了模型及求解方法的準確性。

    圍護結(jié)構(gòu);熱濕耦合傳遞;相對濕度;含濕量

    圍護結(jié)構(gòu)多由多孔介質(zhì)材料構(gòu)建而成,其內(nèi)熱傳遞與濕傳遞屬于典型的多孔介質(zhì)熱質(zhì)傳遞過程。1957年P(guān)hilip等[1]首次提出以溫度和含濕量為驅(qū)動勢,考慮多孔介質(zhì)材料內(nèi)熱傳遞、濕遷移及其耦合作用,建立了多孔介質(zhì)材料熱濕耦合傳遞模型。在Philip研究的基礎上,Pedersen[2],Liesen等[3],Lu[4],Belarbi等[5],Zhong[6],郭興國等[7],Chu等[8],孔凡紅等[9-10],Qin等[11],Leskovsek等[12],Liu等[13]和Vasilyev等[14]分別建立了多孔介質(zhì)材料的熱濕耦合傳遞模型。盡管對圍護結(jié)構(gòu)內(nèi)的熱濕傳遞現(xiàn)象進行了大量研究,但由于圍護結(jié)構(gòu)熱濕耦合傳遞模型為變系數(shù)偏微分方程組,高度非線性且相互耦合,如何簡便求解圍護結(jié)構(gòu)熱濕耦合傳遞模型仍是圍護結(jié)構(gòu)熱濕耦合研究領域的一大難題。

    目前,求解圍護結(jié)構(gòu)熱濕耦合傳遞模型常用的方法是先將控制方程和邊界條件用有限差分法[15]、有限容積法[16]或有限元法[17]進行離散,將偏微分方程組轉(zhuǎn)化為代數(shù)方程組,然后設計算法,編程求解離散后得到的代數(shù)方程組,從而得到溫度和濕度在離散點處的值。由于熱濕耦合傳遞模型的離散過程繁瑣、復雜,需深厚的數(shù)學基礎,而程序算法設計需熟練掌握某種程序設計語言,且模型修改后需重新對其進行離散,這在很大程度上使研究僅停留在理論層面,難以應用于工程實際。

    首先建立圍護結(jié)構(gòu)熱濕耦合傳遞數(shù)學模型,然后提出圍護結(jié)構(gòu)熱濕耦合傳遞模型簡便求解方法,最后通過對比模型模擬結(jié)果與HAMSTAD驗證實例,驗證模型及簡便求解方法的準確性。

    1 圍護結(jié)構(gòu)熱濕耦合傳遞模型

    1.1 控制方程

    為了避免驅(qū)動勢在交界面處不連續(xù),采用連續(xù)變量,溫度和相對濕度,作為驅(qū)動勢,根據(jù)單元體質(zhì)量和能量守恒建立建筑圍護結(jié)構(gòu)熱濕耦合傳遞非穩(wěn)態(tài)模型。

    1.1.1 濕控制方程 雖然圍護結(jié)構(gòu)內(nèi)氣液兩相濕流動不能嚴格的區(qū)分開來,但計算圍護結(jié)構(gòu)內(nèi)濕流動的一個有效方法是將濕流量分為水蒸氣擴散和液態(tài)水遷移兩部分來計算。根據(jù)單元體質(zhì)量守恒,得

    (1)

    式中:ω為含濕量,kg/m3;t為時間,s;jv為水蒸氣擴散速率,kg/(m2·s);jl為液態(tài)水傳遞速率,kg/(m2·s)。

    根據(jù)菲克定律,得

    jv=-δp▽Pv

    (2)

    式中:δp為水蒸氣滲透率,kg/(m·s·Pa);Pv為水蒸氣分壓力,Pa。

    根據(jù)達西定律,得

    jl=Kl▽Pk

    (3)

    式中:Kl為液態(tài)水滲透率,kg/(m·s·Pa);Pk為毛細水壓力,Pa。

    將式(2)、(3)代入式(1),得

    (4)

    體積含濕量是溫度和相對濕度的函數(shù),即

    ω=f(φ,T)

    (5)

    式中:T為溫度,K;φ為相對濕度。

    將方程(5)兩邊同時對時間t求偏導,得

    (6)

    忽略溫度對材料平衡含濕量的影響,則方程(6)可以簡化為

    (7)

    式中:ζ為濕容(等溫吸放曲線的斜率),kg/m3。

    多孔介質(zhì)材料孔隙內(nèi)混合氣體(濕空氣)按理想氣體處理,則

    Pv=φPs

    (8)

    式中:Ps為飽和蒸汽壓力,Pa。

    式(8)兩邊同時求導,得

    ▽Pv=Ps▽φ+φ▽Ps

    (9)

    由于飽和水蒸汽壓力Ps為溫度的單值函數(shù)

    (10)

    將式(10)代入式(9)得

    (11)

    根據(jù)開爾文關(guān)系式,得

    Pk=-ρlRDTln(φ)

    (12)

    式中:ρl為液態(tài)水的密度,kg/m3;RD為水蒸氣氣體常數(shù),J(kg·K)。

    方程(12)兩邊同時求導,得

    (13)

    將式(7)、(11)和(13)代入式(4),得

    (14)

    1.1.2 熱控制方程 根據(jù)單元體能量守恒,控制單元內(nèi)焓的變化等于流入控制單元的凈能量。能量守恒方程可表示為

    (15)

    式中:ρm為干材料密度,kg/m3;cp,m為干材料的比熱,J/(kg·K);ωv為水蒸氣形式的含濕量,kg/m3;ωl為液態(tài)水形式的含濕量,kg/m3;hv為水蒸氣的比焓,J/kg;hl為液態(tài)水的比焓,J/kg;qc為熱流密度,W/m2。

    水蒸氣的比焓可以表示為液態(tài)水的比焓與汽化潛熱之和。

    hv=hl+hlv

    (16)

    式中:hlv為水蒸氣的汽化潛熱,J/kg。

    假設水蒸氣汽化潛熱、干材料和液態(tài)水的比熱為常數(shù),將式(16)代入式(15),整理得

    (17)

    式中:cp,l為液態(tài)水的比熱,J/(kg·K)。

    與水蒸氣汽化潛熱相比,水蒸氣和液態(tài)水的顯熱可以忽略不計,故式(17)右邊第二項可以忽略。盡管水蒸氣的汽化潛熱很大,但由于水蒸汽傳遞速率小,水蒸氣形式的含濕量變化率非常小,式(17)右邊第三項也可以忽略不計,式(17)可簡化為

    (18)

    根據(jù)傅里葉定律,得

    qc=-k▽T

    (19)

    式中:k為導熱系數(shù),W/(m·K)。

    將式(19)代入式(18),得

    (20)

    1.2 邊界條件

    通過圍護結(jié)構(gòu)外表面的濕流量gn,e可表示為:

    gn,e=βp,e(φePs,e-φsurfePs,surfe)

    (21)

    式中:βp,e為圍護結(jié)構(gòu)外表面對流傳質(zhì)系數(shù),kg/(m2·s·Pa);φe為室外相對濕度;φsurfe為圍護結(jié)構(gòu)外表面相對濕度;ps,e為室外飽和水蒸汽壓力,Pa;ps,surfe為圍護結(jié)構(gòu)外表面飽和水蒸汽壓力,Pa。

    通過圍護結(jié)構(gòu)外表面的熱流量qn,e包括對流換熱,太陽輻射及水蒸氣潛熱。

    qn,e=he(Te-Tsurfe)+hlvgn,e+αqsolar

    (22)

    式中:he為外表面對流換熱系數(shù),W/(m2·K);Te為室外空氣溫度,K;Tsurfe為圍護結(jié)構(gòu)外表面溫度,K;α為圍護結(jié)構(gòu)外表面太陽輻射吸收率;qsolar為垂直照射在圍護結(jié)構(gòu)外表面上的太陽輻射,W/m2。

    通過圍護結(jié)構(gòu)內(nèi)表面的濕流量gn,i可表示為

    gn,i=βp,i(φiPs,i-φsurfiPs,surfi)

    (23)

    式中:βp,i為圍護結(jié)構(gòu)內(nèi)表面對流傳質(zhì)系數(shù),kg/(m2·s·Pa);φi為室內(nèi)相對濕度;φsurfi為圍護結(jié)構(gòu)內(nèi)表面相對濕度;ps,i為室內(nèi)飽和水蒸汽壓力,Pa;ps,surfi為圍護結(jié)構(gòu)內(nèi)表面飽和水蒸汽壓力,Pa。

    通過圍護結(jié)構(gòu)內(nèi)表面的熱流量qn,i可表示為

    qn,i=hi(Ti-Tsurfi)+hlvgn,i

    (24)

    式中:hi為內(nèi)表面對流換熱系數(shù),W/(m2·K);Ti為室內(nèi)空氣溫度,K;Tsurfi為圍護結(jié)構(gòu)內(nèi)表面溫度,K。

    2 圍護結(jié)構(gòu)熱濕模型簡便求解方法

    圍護結(jié)構(gòu)熱濕耦合傳遞模型為變系數(shù)偏微分方程組,高度非線性且相互耦合,為了獲得圍護結(jié)構(gòu)內(nèi)的溫濕度分布,控制方程需同時求解。為了避免傳統(tǒng)求解方法中繁瑣、復雜的手動離散過程以及求解算法設計,本文采用多物理場耦合仿真模擬軟件COMSOL來求解圍護結(jié)構(gòu)熱濕耦合傳遞模型。

    COMSOL為工程或數(shù)學問題提供了偏微分方程模型(PDEs),其控制方程和邊界條件為

    (25)

    式中:u為因變量;ea為質(zhì)量系數(shù);da為衰減系數(shù);c為擴散系數(shù);α為守恒通量對流系數(shù);γ為守恒通量源;β為對流系數(shù);a為吸收系數(shù);f為源項;n為朝外的單位向量;g為邊界通量/源;q為邊界吸收/阻抗項。

    圍護結(jié)構(gòu)熱濕耦合傳遞模型中因變量T和φ寫成矩陣形式u=(T,φ),則控制方程在COMSOL的偏微分方程模塊中可表示為

    (26)

    其中:C=ρmcp,m+ω cp,l

    外表面邊界條件可表示為

    (27)

    內(nèi)表面邊界條件可表示為

    (28)

    COMSOL用有限元方法自動對模型進行離散,并用數(shù)值求解器求解離散后的代數(shù)方程組,避免了繁瑣、復雜的手動離散過程以及求解算法設計,研究者和用戶能方便地對已有的模型進行修改或二次開發(fā)。在網(wǎng)格劃分方面COMSOL自帶預定義的三角形網(wǎng)格單元劃分方式,用戶只需選擇合理的網(wǎng)格尺寸(極端粗化、特別粗化、較粗化、粗化、正常、細化、較細化、特別細化、極端細化)對研究對象進行網(wǎng)格劃分。網(wǎng)格尺寸越密,計算越精確,但所需計算時間越長。對于圍護結(jié)構(gòu)熱濕耦合傳遞研究而言,正常的網(wǎng)格尺寸是合理的網(wǎng)格尺寸,但分析雨水吸收問題時,為了計算收斂,需要選取特別細化的網(wǎng)格尺寸。COMSOL求解流程圖如圖1所示。

    圖1 COMSOL求解流程圖

    3 熱濕耦合模型及求解方法驗證

    圍護結(jié)構(gòu)熱濕耦合模型是對實際物理現(xiàn)象和過程的抽象與簡化,模型高度非線性且相互耦合,數(shù)值求解過程中難免會有誤差,因而在模型應用前需對其進行驗證。本文通過對比新模型模擬結(jié)果與HAMSTAD驗證實例來驗證模型及求解方法。

    3.1 HAMSTAD驗證實例2

    HAMSTAD驗證實例2分析了200 mm厚的單層各向同性墻體的等溫干燥過程。初始條件為20 ℃,95%。在開始時刻,周圍環(huán)境的相對濕度突然改變,室外相對濕度變?yōu)?5%,室內(nèi)相對濕度變?yōu)?5%,室內(nèi)外對流換熱系數(shù)均為25 W/(m2·K),室內(nèi)外對流質(zhì)傳遞系數(shù)均為1×10-3s/m 。材料熱濕參數(shù)見表1。驗證實例的詳細描述見文獻[18]。100、300、1 000 h時,墻體內(nèi)的含濕量分布如圖2所示。從對比結(jié)果可以看出,新建模型模擬結(jié)果與HAMSTAD驗證實例2中的分析解吻合良好。

    圖2 100、300、1 000 h時墻體內(nèi)的含濕量分布

    表1 各向同性墻體的材料參數(shù)[18]

    3.2 HAMSTAD驗證實例5

    HAMSTAD驗證實例5分析了3層復合墻體內(nèi)含濕量的變化。墻體構(gòu)造為365 mm磚墻,15 mm砂漿層,40 mm保溫層(由外至內(nèi)),初始條件為25 ℃、60%,室外溫濕度為0 ℃、80%,室內(nèi)溫濕度為20 ℃、60%,室內(nèi)外對流換熱系數(shù)分別為8 W/(m2·K)和25 W/(m2·K),室內(nèi)外對流質(zhì)傳遞系數(shù)分別為5.882 3×10-8s/m和1.838 2×10-7s/m。材料的密度、比熱和導熱系數(shù)見表2。材料的其他熱物性參數(shù)及驗證實例的詳細描述見文獻[18]。從圖3可以看出,新建模型模擬結(jié)果與HAMSTAD驗證實例5中的結(jié)果吻合良好。

    表2 墻體部件的密度、比熱容和導熱系數(shù)[18]

    圖3 60 d后墻體內(nèi)的含濕量分布(TUD, Technion, KUL, TUE, CTH和NRC表示參與HAMSTAD項目的其他研究機構(gòu)的模擬結(jié)果)

    4 結(jié) 論

    本文以多孔介質(zhì)材料傳熱傳質(zhì)理論為基礎,根據(jù)單元體質(zhì)量和能量守恒,建立了圍護結(jié)構(gòu)熱濕耦合傳遞非穩(wěn)態(tài)模型,并提出了基于多物理場耦合仿真模擬軟件COMSOL的熱濕耦合傳遞模型簡便求解方法。通過對比新建模型模擬結(jié)果與HAMSTAD驗證實例2和5,驗證了模型及求解方法的準確性。該模型可為優(yōu)化圍護結(jié)構(gòu)熱工性能,預測圍護結(jié)構(gòu)濕損壞提供理論依據(jù)和技術(shù)指導。

    [1] PHILIP J R,DEVRIES D A.Moisture movement in porous materials under temperature gradients [J].Transactions American Geophysical Union,1957,38(2):222-232.

    [2] PEDERSEN C R.Combined heat and moisture transfer in building constructions [D].Copenhagen:Technical University of Denmark,1990:24-56.

    [3] LIESEN R,PEDERSEN C.Modelling the energy effects of combined heat and mass transfer in building elements: Part I—Theory [J].ASHRAE Transactions,1999,105(1):941-953.

    [4] LU X.Modelling of heat and moisture transfer in buildings I. model program [J].Energy and Buildings,2002,34(10):1033-1043.

    [5] BELARBI R,QIN M H,ABDELKARIM A M,et al.Experimental and theoretical investigation of non-isothermal transfer in hygroscopic building materials [J].Building and Environment,2008,43(12):2154-2162.

    [6] ZHONG Z P.Combined heat and moisture transport modeling for residential buildings [D].Indiana:Purdue University,2008:18-24.

    [7] 郭興國,陳友明,鄧永強.Budaiwi模型的修正及實驗驗證[J].土木建筑與環(huán)境工程,2010,32(2);91-95.

    GUO X G,CHEN Y M,DENG Y Q.A modified buidaiwi model and its experimental evaluation[J].Journal of Civil,Architectural and Environmental Engineering,2010,32(2):91-95.(in Chinese)

    [8] CHU S S,F(xiàn)ANG T H,CHANG W J.Modelling of coupled heat and moisture transfer in porous construction materials [J].Mathematical and Computer Modelling, 2009,50(7/8):1195-1204.

    [9] 孔凡紅.多層多孔介質(zhì)層間界面處的熱質(zhì)耦合傳遞[J].太陽能學報,2010,31(10):1281-1286.

    KONG F H.Coupled transfer of heat and mass at interface of multilayer porous media[J],2010,31(10):1281-1286.(in Chinese)

    [10] KONG F H,ZHANG Q L.Effect of heat and mass coupled transfer combined with freezing process on building exterior envelope [J].Energy and Buildings,2013,62:486-495.

    [11] QIN M H,ABDELKARIM A M,BELARBI R.Two-dimensional hygrothermal transfer in porous building materials [J].Applied Thermal Engineering,2010,30(16):2555-2562.

    [12] LESKOVSEK U,MEDVED S.Heat and moisture transfer in fibrous thermal insulation with tight boundaries and a dynamical boundary temperature [J].International Journal of Heat and Mass Transfer, 2011,54(19/20):4333-4340.

    [13] LIU Y F,WANG Y Y,WANG D J,et al.Effect of moisture transfer on internal surface temperature [J].Energy and Buildings,2013,60:83-91.

    [14] VASILYEV G P,LICHMAN V A,PESKOV N V,et al.Simulation of heat and moisture transfer in a multiplex structure [J].Energy and Buildings,2015,86:803-807.

    [15] KALIS H.Efficient finite-difference scheme for solving some heat transfer problems with convection in multilayer media [J].International Journal of Heat and Mass Transfer,2000,43(24):4467-4474.

    [16] 張景欣,郭興國,陳友明,等.墻體內(nèi)熱濕耦合過程的時域遞歸展開算法[J].土木建筑與環(huán)境工程,2015,37(6);147-152.

    ZHANG J X,GUO X G,CHEN Y M,et al.A time domain recursive algorithm for solving the model of coupled heat and moisture transfer model in building wall[J].Journal of Civil,Architectural and Environmental Engineering,2015,37(6):147-152. (in Chinese)

    [17] KHOSHBAKHT M,LIN M,F(xiàn)EICKERT C.A finite element model for hygrothermal analysis of masonry walls with FRP reinforcement [J].Finite Elements in Analysis and Design,2009,45(8/9):511-518.

    [18] HAGENTOFT C E.HAMSTAD-Final report:Methodology of HAM-modeling,Report R-02:8 [R].Gothenburg:Department of Building Physics,Chalmers University of Technology,2002:5-51.

    (編輯 胡玲)

    Coupled heat and moisture transfer model and simple solution method for building envelopes

    Liu Xiangwei1, Chen Youming2, Chen Guojie2,3, Guo Xingguo1, Luo Na1

    (1. School of Civil Engineering and Architecture, Nanchang University, Nanchang 330031, P.R. China; 2. College of Civil Engineering, Hunan University, Changsha 410082, P.R. China; 3. College of City Construction, University of South China, Hengyang 421001,Hunan, P.R. China)

    A coupled heat and moisture transfer model which takes the coupled effect between heat transfer and moisture transfer into consideration is developed to model the temperature and humidity distribution within building envelopes. The continuous variables, relative humidity and temperature are chosen as the driving potentials. And a simple solution method based on a multi-physics simulation soft, COMSOL, is proposed to solve the coupled heat and moisture transfer model simply. The numerical results of this model agree well with HAMSTAD benchmarks which are international accepted.

    building envelope; coupled heat and moisture transfer; relative humidity; moisture content

    2016-04-12

    國家自然科學基金 (51208247、51408294)

    劉向偉(1987-),男,博士,主要從事建筑節(jié)能及圍護結(jié)構(gòu)熱濕耦合傳遞研究,(E-mail)xiangwei.liu@ncu.edu.cn。

    郭興國(通信作者),男,副教授,(Email)guoxingguo@ncu.edu.cn。

    Foundation item:Natural Science Foundation of China (No.51208247, No.51408294)

    TU111.4

    A

    1674-4764(2016)04-0007-06

    10.11835/j.issn.1674-4764.2016.04.002

    Received:2016-04-12

    Author brief:Liu Xiangwei (1987-), PhD, main research interests: building efficiency and coupled heat and moisture transfer in building envelopes, (E-mail)xiangwei.liu@ncu.edu.cn.

    Guo Xingguo(corresponding author), associate professor, (E-mail)guoxingguo@ncu.edu.cn.

    猜你喜歡
    含濕量液態(tài)水水蒸氣
    藏在空氣里的水蒸氣
    科學大眾(2022年23期)2023-01-30 07:04:00
    基于微波輻射計的張掖地區(qū)水汽、液態(tài)水變化特征分析
    Ka/Ku雙波段毫米波雷達功率譜數(shù)據(jù)反演液態(tài)水含量方法研究
    非分散紅外法測定煙氣中二氧化硫常見問題的探討
    濕空氣含濕量的測量方法及其表達式的研究與應用
    新型一體式除濕熱泵空調(diào)循環(huán)夏季工況實驗研究
    制冷學報(2018年5期)2018-10-16 11:06:56
    零下溫度的液態(tài)水
    水蒸氣變戲法
    PEMFC氣體擴散層中液態(tài)水傳輸實驗研究綜述
    迷宮篇
    成人av一区二区三区在线看| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩无卡精品| 亚洲 欧美 日韩 在线 免费| 国产精品国产高清国产av| 啦啦啦观看免费观看视频高清| 最近最新中文字幕大全电影3 | 男人的好看免费观看在线视频 | 中文字幕人妻丝袜一区二区| 亚洲精品色激情综合| 国产精品二区激情视频| 90打野战视频偷拍视频| 国产成人精品无人区| 一本久久中文字幕| 啦啦啦 在线观看视频| 夜夜看夜夜爽夜夜摸| 国产乱人伦免费视频| 99久久99久久久精品蜜桃| 亚洲在线自拍视频| 亚洲一区二区三区不卡视频| 麻豆成人av在线观看| 自线自在国产av| 久久99热这里只有精品18| 欧美激情极品国产一区二区三区| or卡值多少钱| 国产精品久久视频播放| 成年免费大片在线观看| 成人18禁在线播放| 一区二区三区国产精品乱码| 欧美不卡视频在线免费观看 | 午夜日韩欧美国产| 午夜日韩欧美国产| 精品国产乱子伦一区二区三区| 中文在线观看免费www的网站 | 国产精品精品国产色婷婷| 国产精品爽爽va在线观看网站 | 欧美黄色淫秽网站| 国产一区二区激情短视频| 久久狼人影院| 淫妇啪啪啪对白视频| 免费一级毛片在线播放高清视频| 亚洲av片天天在线观看| 人妻丰满熟妇av一区二区三区| 国产精品野战在线观看| 无限看片的www在线观看| 国产成人欧美| 一区二区日韩欧美中文字幕| 久久久久久久久久黄片| 国产精品自产拍在线观看55亚洲| 久久久久久久午夜电影| 51午夜福利影视在线观看| 两性夫妻黄色片| 久久久久九九精品影院| 人成视频在线观看免费观看| 琪琪午夜伦伦电影理论片6080| 免费av毛片视频| 日韩三级视频一区二区三区| 香蕉丝袜av| 国产精品久久视频播放| tocl精华| 国产极品粉嫩免费观看在线| 中亚洲国语对白在线视频| 哪里可以看免费的av片| 欧美成人午夜精品| 亚洲精品中文字幕一二三四区| 国产精品免费一区二区三区在线| av欧美777| 午夜免费成人在线视频| 午夜两性在线视频| 欧美性长视频在线观看| 观看免费一级毛片| 在线播放国产精品三级| 国产精品 国内视频| 国产精品亚洲美女久久久| 一级毛片高清免费大全| 国产精品99久久99久久久不卡| 欧美国产精品va在线观看不卡| 亚洲av五月六月丁香网| 日本黄色视频三级网站网址| av中文乱码字幕在线| 亚洲精品一区av在线观看| 18禁黄网站禁片免费观看直播| 日韩 欧美 亚洲 中文字幕| 91麻豆精品激情在线观看国产| 午夜福利18| 久久亚洲精品不卡| 伊人久久大香线蕉亚洲五| 精品一区二区三区视频在线观看免费| 一本大道久久a久久精品| 欧美精品啪啪一区二区三区| 妹子高潮喷水视频| 午夜视频精品福利| 黄网站色视频无遮挡免费观看| 麻豆成人午夜福利视频| 精品乱码久久久久久99久播| 久久精品国产亚洲av香蕉五月| 午夜福利高清视频| 99国产精品一区二区三区| 亚洲av五月六月丁香网| 亚洲精品中文字幕在线视频| 亚洲一区二区三区色噜噜| 18禁观看日本| 2021天堂中文幕一二区在线观 | 成人18禁在线播放| 亚洲男人的天堂狠狠| 狂野欧美激情性xxxx| 18禁黄网站禁片午夜丰满| av在线天堂中文字幕| 国产精品亚洲一级av第二区| 国产日本99.免费观看| 最近最新免费中文字幕在线| 人成视频在线观看免费观看| 国产黄色小视频在线观看| a级毛片a级免费在线| 草草在线视频免费看| 国产成人av教育| 99国产精品一区二区蜜桃av| 91成人精品电影| 亚洲成人精品中文字幕电影| 亚洲国产精品sss在线观看| 身体一侧抽搐| 亚洲中文字幕日韩| 成人特级黄色片久久久久久久| 中文字幕高清在线视频| 在线观看日韩欧美| 色播亚洲综合网| 精品国产乱码久久久久久男人| 琪琪午夜伦伦电影理论片6080| 欧美在线一区亚洲| 美国免费a级毛片| 国产成人精品久久二区二区免费| 精品福利观看| 韩国精品一区二区三区| 男女午夜视频在线观看| 午夜福利在线在线| 国产蜜桃级精品一区二区三区| 久久香蕉国产精品| a级毛片a级免费在线| 夜夜爽天天搞| 久久九九热精品免费| 国产精品免费视频内射| 桃色一区二区三区在线观看| 啦啦啦韩国在线观看视频| 欧美日韩瑟瑟在线播放| 久热爱精品视频在线9| 成人av一区二区三区在线看| 亚洲国产看品久久| 深夜精品福利| 国产亚洲欧美在线一区二区| 免费高清视频大片| 国产三级在线视频| 久久99热这里只有精品18| 欧美性猛交黑人性爽| 国产高清有码在线观看视频 | av中文乱码字幕在线| 熟妇人妻久久中文字幕3abv| 精品不卡国产一区二区三区| 亚洲av成人不卡在线观看播放网| 久久国产乱子伦精品免费另类| 精品电影一区二区在线| 一区福利在线观看| 亚洲男人天堂网一区| 99久久无色码亚洲精品果冻| 欧美色欧美亚洲另类二区| 99在线视频只有这里精品首页| 欧美另类亚洲清纯唯美| 精品高清国产在线一区| 亚洲欧美精品综合一区二区三区| 精品国产超薄肉色丝袜足j| 亚洲美女黄片视频| 中文字幕人成人乱码亚洲影| 亚洲熟妇熟女久久| 日韩av在线大香蕉| 1024视频免费在线观看| 国产亚洲欧美98| 国产片内射在线| 久久精品亚洲精品国产色婷小说| 久久久久久免费高清国产稀缺| 日韩成人在线观看一区二区三区| 亚洲国产看品久久| 日韩欧美一区二区三区在线观看| 久久国产精品影院| 国产区一区二久久| 亚洲avbb在线观看| 国产精品亚洲美女久久久| 国产一区二区三区视频了| 12—13女人毛片做爰片一| 久久亚洲精品不卡| 久久久久久久久久黄片| 视频在线观看一区二区三区| 免费在线观看黄色视频的| 色av中文字幕| www国产在线视频色| 亚洲成av片中文字幕在线观看| 波多野结衣高清作品| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 亚洲av成人一区二区三| 亚洲专区国产一区二区| 人人妻人人澡欧美一区二区| 两个人免费观看高清视频| 久久精品影院6| 成在线人永久免费视频| 国产国语露脸激情在线看| 级片在线观看| 国产成+人综合+亚洲专区| 女性被躁到高潮视频| 亚洲精品国产区一区二| 怎么达到女性高潮| 一区二区三区精品91| 在线观看免费午夜福利视频| 一进一出抽搐动态| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩高清在线视频| 啦啦啦韩国在线观看视频| 啦啦啦观看免费观看视频高清| 久久精品国产亚洲av高清一级| 又大又爽又粗| 三级毛片av免费| 国产片内射在线| 亚洲欧美一区二区三区黑人| 国产91精品成人一区二区三区| 免费高清视频大片| 国产欧美日韩精品亚洲av| 麻豆久久精品国产亚洲av| 91字幕亚洲| 日本撒尿小便嘘嘘汇集6| 十分钟在线观看高清视频www| 我的亚洲天堂| 亚洲五月色婷婷综合| 亚洲 国产 在线| 日韩精品免费视频一区二区三区| 老汉色av国产亚洲站长工具| 亚洲黑人精品在线| 日韩有码中文字幕| 成人一区二区视频在线观看| 国产精品1区2区在线观看.| 窝窝影院91人妻| 久久久久免费精品人妻一区二区 | 老司机午夜十八禁免费视频| 国产单亲对白刺激| 国产成年人精品一区二区| 欧美激情高清一区二区三区| 午夜老司机福利片| 麻豆一二三区av精品| 成人特级黄色片久久久久久久| 啦啦啦免费观看视频1| 听说在线观看完整版免费高清| 一区二区三区精品91| 国产激情久久老熟女| 又黄又爽又免费观看的视频| 波多野结衣av一区二区av| 波多野结衣高清作品| 国产成人精品无人区| aaaaa片日本免费| 50天的宝宝边吃奶边哭怎么回事| 不卡av一区二区三区| 免费看a级黄色片| 国产精品精品国产色婷婷| 99精品久久久久人妻精品| 在线看三级毛片| 亚洲全国av大片| 日韩免费av在线播放| 欧美又色又爽又黄视频| 日韩欧美一区视频在线观看| 国产亚洲精品久久久久5区| 久久久精品欧美日韩精品| 99riav亚洲国产免费| 精品无人区乱码1区二区| 免费看日本二区| 亚洲精品中文字幕一二三四区| 国产视频内射| 亚洲午夜理论影院| 麻豆成人av在线观看| 亚洲精品在线美女| 日韩欧美一区视频在线观看| 一夜夜www| 亚洲国产中文字幕在线视频| 熟女电影av网| 亚洲精品国产精品久久久不卡| 亚洲成国产人片在线观看| 天天一区二区日本电影三级| 久久亚洲精品不卡| 久久久国产精品麻豆| 欧美色欧美亚洲另类二区| 欧美最黄视频在线播放免费| 最好的美女福利视频网| 国产av不卡久久| 亚洲人成电影免费在线| 黑人欧美特级aaaaaa片| 日本黄色视频三级网站网址| 欧美日本亚洲视频在线播放| 变态另类丝袜制服| 别揉我奶头~嗯~啊~动态视频| 国产97色在线日韩免费| 国产视频内射| 欧美日韩福利视频一区二区| 色尼玛亚洲综合影院| 亚洲国产日韩欧美精品在线观看 | 淫妇啪啪啪对白视频| 国产亚洲av高清不卡| 国产一区在线观看成人免费| 国产高清有码在线观看视频 | 免费在线观看视频国产中文字幕亚洲| www.精华液| 精品国产国语对白av| 女人高潮潮喷娇喘18禁视频| 草草在线视频免费看| 制服丝袜大香蕉在线| 老司机午夜福利在线观看视频| 在线国产一区二区在线| 中文亚洲av片在线观看爽| 国产亚洲欧美在线一区二区| 亚洲成国产人片在线观看| 欧美不卡视频在线免费观看 | 国产精品久久久久久亚洲av鲁大| 90打野战视频偷拍视频| 两人在一起打扑克的视频| 麻豆一二三区av精品| 国产爱豆传媒在线观看 | 一本久久中文字幕| 观看免费一级毛片| 黑人巨大精品欧美一区二区mp4| 无遮挡黄片免费观看| 国产一区二区在线av高清观看| 色精品久久人妻99蜜桃| 人人澡人人妻人| 在线观看免费视频日本深夜| 美女免费视频网站| 99热只有精品国产| 亚洲国产高清在线一区二区三 | 少妇 在线观看| 一进一出好大好爽视频| 精品久久久久久久末码| 亚洲五月婷婷丁香| 又大又爽又粗| 成人18禁高潮啪啪吃奶动态图| 自线自在国产av| 日本 av在线| 日本a在线网址| 欧洲精品卡2卡3卡4卡5卡区| 成人欧美大片| 亚洲中文字幕一区二区三区有码在线看 | 日日爽夜夜爽网站| 国产精品综合久久久久久久免费| 国产一区二区三区在线臀色熟女| 色精品久久人妻99蜜桃| 最近最新中文字幕大全电影3 | 久久精品91蜜桃| 亚洲精品美女久久久久99蜜臀| xxxwww97欧美| 黄片播放在线免费| 久久伊人香网站| 一本一本综合久久| 精品高清国产在线一区| 伊人久久大香线蕉亚洲五| 黄色视频,在线免费观看| 三级毛片av免费| 高清在线国产一区| 不卡一级毛片| 国产欧美日韩精品亚洲av| 夜夜夜夜夜久久久久| 国产精品98久久久久久宅男小说| 国产精品日韩av在线免费观看| 亚洲成国产人片在线观看| 熟女电影av网| 久久香蕉国产精品| 亚洲一区高清亚洲精品| 亚洲午夜精品一区,二区,三区| 欧美在线黄色| 亚洲aⅴ乱码一区二区在线播放 | 人成视频在线观看免费观看| 丁香六月欧美| 亚洲激情在线av| 两性夫妻黄色片| 黑人操中国人逼视频| 国产精品久久久av美女十八| 亚洲国产精品成人综合色| 在线播放国产精品三级| 在线观看免费视频日本深夜| 哪里可以看免费的av片| 两人在一起打扑克的视频| 中文字幕精品免费在线观看视频| 男女做爰动态图高潮gif福利片| 一本一本综合久久| 欧美激情久久久久久爽电影| 亚洲精品美女久久久久99蜜臀| 国产高清有码在线观看视频 | 1024视频免费在线观看| 国产高清videossex| 午夜成年电影在线免费观看| 国产伦一二天堂av在线观看| 日韩精品青青久久久久久| 欧美av亚洲av综合av国产av| 亚洲av电影在线进入| 久久久精品国产亚洲av高清涩受| 一二三四在线观看免费中文在| 正在播放国产对白刺激| 久久草成人影院| 国产精品久久视频播放| 91成年电影在线观看| 国产一区二区激情短视频| 精品免费久久久久久久清纯| 婷婷精品国产亚洲av| 极品教师在线免费播放| 日韩av在线大香蕉| 国产精品久久久久久人妻精品电影| 给我免费播放毛片高清在线观看| 人成视频在线观看免费观看| 久99久视频精品免费| 午夜精品久久久久久毛片777| 后天国语完整版免费观看| 成年免费大片在线观看| 中文字幕av电影在线播放| 欧美绝顶高潮抽搐喷水| 9191精品国产免费久久| 欧美性长视频在线观看| 18美女黄网站色大片免费观看| 久久性视频一级片| 亚洲天堂国产精品一区在线| 麻豆成人av在线观看| 亚洲激情在线av| 国产精华一区二区三区| 亚洲第一欧美日韩一区二区三区| 欧美日韩福利视频一区二区| 法律面前人人平等表现在哪些方面| 看免费av毛片| 亚洲色图av天堂| 丝袜在线中文字幕| 日韩欧美一区视频在线观看| 黑丝袜美女国产一区| 欧美乱妇无乱码| 18禁裸乳无遮挡免费网站照片 | 最近最新免费中文字幕在线| 精品欧美国产一区二区三| 国产成人啪精品午夜网站| 无限看片的www在线观看| 18禁黄网站禁片免费观看直播| 久久香蕉国产精品| 国产精品av久久久久免费| 又大又爽又粗| 久久草成人影院| 欧美绝顶高潮抽搐喷水| 精品久久久久久,| 亚洲激情在线av| 男人的好看免费观看在线视频 | 在线天堂中文资源库| 两人在一起打扑克的视频| 久久久久国产精品人妻aⅴ院| 久久狼人影院| 国产在线观看jvid| videosex国产| 麻豆国产av国片精品| svipshipincom国产片| 91老司机精品| 亚洲成人免费电影在线观看| 亚洲国产看品久久| 欧美成人性av电影在线观看| 非洲黑人性xxxx精品又粗又长| 国产1区2区3区精品| 一级毛片女人18水好多| 不卡一级毛片| 精品乱码久久久久久99久播| 亚洲五月天丁香| 99久久综合精品五月天人人| 黄色毛片三级朝国网站| 国产亚洲av高清不卡| av在线播放免费不卡| 狠狠狠狠99中文字幕| 十分钟在线观看高清视频www| 看片在线看免费视频| 久久久久久免费高清国产稀缺| 亚洲国产看品久久| а√天堂www在线а√下载| 国产国语露脸激情在线看| 中文字幕久久专区| 午夜视频精品福利| 国产精品久久久人人做人人爽| 日本免费a在线| 老司机午夜十八禁免费视频| 午夜久久久久精精品| 欧美性长视频在线观看| 在线视频色国产色| 久久中文字幕人妻熟女| 日本撒尿小便嘘嘘汇集6| 精品国产乱子伦一区二区三区| 不卡一级毛片| 啦啦啦 在线观看视频| 后天国语完整版免费观看| 久久热在线av| 777久久人妻少妇嫩草av网站| 国产高清激情床上av| 亚洲精华国产精华精| 亚洲成人久久爱视频| 日日爽夜夜爽网站| 黄色a级毛片大全视频| 免费无遮挡裸体视频| 啦啦啦免费观看视频1| 久久久久久久精品吃奶| 校园春色视频在线观看| 婷婷亚洲欧美| 精华霜和精华液先用哪个| 精品国产乱码久久久久久男人| 好男人在线观看高清免费视频 | 精品无人区乱码1区二区| 人妻久久中文字幕网| 欧美精品亚洲一区二区| 一进一出好大好爽视频| 国产精品 国内视频| 91大片在线观看| 久久中文看片网| 老司机福利观看| 日本免费一区二区三区高清不卡| 中亚洲国语对白在线视频| 麻豆国产av国片精品| 中文字幕高清在线视频| 国产不卡一卡二| bbb黄色大片| 久久香蕉精品热| а√天堂www在线а√下载| 亚洲av片天天在线观看| 久久久久久久久久黄片| 欧美日韩亚洲国产一区二区在线观看| 久久欧美精品欧美久久欧美| 精品久久久久久久末码| 色老头精品视频在线观看| 嫩草影院精品99| 我的亚洲天堂| 色综合亚洲欧美另类图片| www国产在线视频色| 国产精品久久久久久亚洲av鲁大| 久久久久九九精品影院| 美女国产高潮福利片在线看| 国产精品日韩av在线免费观看| 亚洲第一电影网av| 人人妻,人人澡人人爽秒播| 午夜福利视频1000在线观看| 亚洲最大成人中文| 日韩视频一区二区在线观看| 亚洲 欧美一区二区三区| 少妇被粗大的猛进出69影院| 狂野欧美激情性xxxx| 午夜福利在线在线| 高清毛片免费观看视频网站| 国产精品一区二区精品视频观看| 日韩欧美免费精品| 真人做人爱边吃奶动态| 免费电影在线观看免费观看| 丰满人妻熟妇乱又伦精品不卡| av在线播放免费不卡| 狂野欧美激情性xxxx| 91成人精品电影| 免费在线观看成人毛片| 亚洲专区中文字幕在线| 欧美成狂野欧美在线观看| 婷婷亚洲欧美| 亚洲av第一区精品v没综合| 久久久久久久久中文| 国产午夜精品久久久久久| 国产aⅴ精品一区二区三区波| 一a级毛片在线观看| 亚洲第一欧美日韩一区二区三区| 2021天堂中文幕一二区在线观 | 亚洲av片天天在线观看| 国产精品日韩av在线免费观看| 制服人妻中文乱码| 国产免费av片在线观看野外av| 久久精品国产99精品国产亚洲性色| 欧美一区二区精品小视频在线| 午夜福利视频1000在线观看| 麻豆成人av在线观看| 麻豆av在线久日| 亚洲专区国产一区二区| 久久精品91无色码中文字幕| 老熟妇乱子伦视频在线观看| 叶爱在线成人免费视频播放| 国产1区2区3区精品| 变态另类丝袜制服| 午夜免费鲁丝| 欧美激情 高清一区二区三区| √禁漫天堂资源中文www| 国产不卡一卡二| 亚洲avbb在线观看| 精品国产一区二区三区四区第35| 狂野欧美激情性xxxx| 国产成人系列免费观看| 欧美丝袜亚洲另类 | 午夜福利高清视频| 真人一进一出gif抽搐免费| 黄网站色视频无遮挡免费观看| 亚洲五月天丁香| 又紧又爽又黄一区二区| 久久久久国内视频| 精品少妇一区二区三区视频日本电影| 国产精品久久久久久人妻精品电影| 又黄又爽又免费观看的视频| 亚洲av熟女| 亚洲va日本ⅴa欧美va伊人久久| 国产伦人伦偷精品视频| 久久久久免费精品人妻一区二区 | 亚洲成av人片免费观看| 妹子高潮喷水视频| 国产亚洲av嫩草精品影院| 亚洲熟女毛片儿| 俺也久久电影网| 欧美日韩黄片免| 老鸭窝网址在线观看| 又黄又爽又免费观看的视频| 男女视频在线观看网站免费 | 欧美激情久久久久久爽电影| 色播亚洲综合网| 中文字幕精品免费在线观看视频| 亚洲精品中文字幕一二三四区| 日本熟妇午夜| av视频在线观看入口|