• <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)水傳輸實驗研究綜述
    迷宮篇
    日本欧美国产在线视频| 黄色一级大片看看| 在线天堂中文资源库| 亚洲欧美一区二区三区黑人 | www日本在线高清视频| 国产精品久久久久久精品古装| 十八禁高潮呻吟视频| 99国产精品免费福利视频| 亚洲国产日韩一区二区| 精品人妻在线不人妻| 亚洲内射少妇av| 久久这里只有精品19| 一级片免费观看大全| 国产成人精品久久二区二区91 | 婷婷色av中文字幕| 激情视频va一区二区三区| 一级a爱视频在线免费观看| 亚洲国产精品成人久久小说| 亚洲综合精品二区| 午夜激情av网站| 另类精品久久| 欧美日韩一区二区视频在线观看视频在线| 男女无遮挡免费网站观看| 婷婷色av中文字幕| 国产午夜精品一二区理论片| 亚洲国产精品一区二区三区在线| xxx大片免费视频| 亚洲国产欧美网| 如何舔出高潮| 日本黄色日本黄色录像| 99香蕉大伊视频| 成人漫画全彩无遮挡| 欧美日韩av久久| 日本欧美视频一区| 男女下面插进去视频免费观看| 色吧在线观看| 午夜日韩欧美国产| 欧美成人午夜免费资源| 欧美日本中文国产一区发布| 两性夫妻黄色片| av一本久久久久| 黄色视频在线播放观看不卡| 天天影视国产精品| 中文天堂在线官网| 精品一区在线观看国产| 日韩一区二区三区影片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩中文字幕视频在线看片| 在现免费观看毛片| 国产成人精品一,二区| 国产高清不卡午夜福利| 99久久中文字幕三级久久日本| 午夜福利视频精品| 高清在线视频一区二区三区| 亚洲国产精品一区二区三区在线| 免费av中文字幕在线| 久久久久网色| 成年动漫av网址| 亚洲欧美一区二区三区黑人 | 秋霞伦理黄片| 伊人亚洲综合成人网| 中国三级夫妇交换| 不卡av一区二区三区| 国产乱来视频区| 国产精品一区二区在线不卡| 丝袜人妻中文字幕| 视频区图区小说| 人人妻人人澡人人看| 国产精品无大码| 日韩中字成人| 国产一区二区 视频在线| 国产精品免费视频内射| 日韩在线高清观看一区二区三区| 久久久亚洲精品成人影院| 欧美bdsm另类| 日本wwww免费看| 久久久久视频综合| 波多野结衣一区麻豆| 欧美日韩成人在线一区二区| av免费观看日本| 18禁裸乳无遮挡动漫免费视频| 色哟哟·www| 久久久久久久精品精品| 热re99久久国产66热| 中文乱码字字幕精品一区二区三区| 少妇人妻精品综合一区二区| 晚上一个人看的免费电影| videos熟女内射| 少妇的逼水好多| 女人久久www免费人成看片| 精品国产乱码久久久久久小说| 黄色毛片三级朝国网站| 日本91视频免费播放| 肉色欧美久久久久久久蜜桃| 一级毛片 在线播放| 咕卡用的链子| 80岁老熟妇乱子伦牲交| 香蕉国产在线看| 9色porny在线观看| 免费黄网站久久成人精品| 成人亚洲精品一区在线观看| 十分钟在线观看高清视频www| 日韩免费高清中文字幕av| 国产亚洲av片在线观看秒播厂| 色94色欧美一区二区| 日本vs欧美在线观看视频| 一级毛片我不卡| 久久久a久久爽久久v久久| 大香蕉久久网| 欧美日本中文国产一区发布| 国产av精品麻豆| 麻豆av在线久日| 99香蕉大伊视频| 国产一区二区 视频在线| 在线 av 中文字幕| 国产成人欧美| 亚洲第一av免费看| 国产精品 国内视频| 中文字幕最新亚洲高清| 欧美 日韩 精品 国产| 免费大片黄手机在线观看| 搡女人真爽免费视频火全软件| 男男h啪啪无遮挡| 男人爽女人下面视频在线观看| 在线天堂最新版资源| 国产成人精品福利久久| 日韩熟女老妇一区二区性免费视频| 天堂中文最新版在线下载| 超色免费av| 国产一区二区在线观看av| 一级毛片电影观看| 欧美人与性动交α欧美精品济南到 | 国产亚洲最大av| 国产av国产精品国产| 久久国内精品自在自线图片| 亚洲三区欧美一区| 在线精品无人区一区二区三| 色婷婷av一区二区三区视频| 国产精品99久久99久久久不卡 | 久久99精品国语久久久| 国产成人av激情在线播放| 久久久久久人人人人人| videosex国产| 午夜激情久久久久久久| 婷婷色av中文字幕| 欧美日韩综合久久久久久| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品成人久久小说| 不卡视频在线观看欧美| 国产精品亚洲av一区麻豆 | 亚洲,欧美精品.| 91aial.com中文字幕在线观看| 91精品国产国语对白视频| 日日啪夜夜爽| 午夜日本视频在线| 嫩草影院入口| 80岁老熟妇乱子伦牲交| 国产一区二区在线观看av| 国产片特级美女逼逼视频| 1024香蕉在线观看| 一级毛片 在线播放| 下体分泌物呈黄色| 五月开心婷婷网| 99久久中文字幕三级久久日本| 丝袜脚勾引网站| 亚洲欧美成人精品一区二区| 国产97色在线日韩免费| 性色av一级| 国产成人精品在线电影| 欧美日韩亚洲国产一区二区在线观看 | 制服诱惑二区| 午夜福利视频精品| 久久毛片免费看一区二区三区| 亚洲欧美一区二区三区黑人 | 黄网站色视频无遮挡免费观看| tube8黄色片| 亚洲激情五月婷婷啪啪| 亚洲成人手机| 日韩在线高清观看一区二区三区| 性少妇av在线| 久久久精品区二区三区| 欧美xxⅹ黑人| 观看美女的网站| 新久久久久国产一级毛片| 久久人妻熟女aⅴ| 精品一区在线观看国产| 午夜精品国产一区二区电影| 麻豆av在线久日| 国精品久久久久久国模美| 国产日韩欧美亚洲二区| 天天躁夜夜躁狠狠久久av| 中文字幕av电影在线播放| 亚洲欧美色中文字幕在线| 1024视频免费在线观看| 午夜福利影视在线免费观看| 精品少妇黑人巨大在线播放| www.精华液| 国产成人欧美| 老鸭窝网址在线观看| xxxhd国产人妻xxx| 久热这里只有精品99| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区在线不卡| 久久久久人妻精品一区果冻| 久久精品久久久久久噜噜老黄| 久久久国产一区二区| 亚洲欧美一区二区三区国产| 在线天堂中文资源库| 在线观看三级黄色| 纯流量卡能插随身wifi吗| 香蕉丝袜av| 国产片内射在线| 自线自在国产av| 高清黄色对白视频在线免费看| 国产精品一区二区在线不卡| 久久久久久久亚洲中文字幕| 99热国产这里只有精品6| 亚洲av免费高清在线观看| 日韩在线高清观看一区二区三区| 国产亚洲欧美精品永久| 自线自在国产av| 电影成人av| 久久97久久精品| 久久狼人影院| 亚洲综合色惰| 亚洲精品中文字幕在线视频| 男的添女的下面高潮视频| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂| 久久精品亚洲av国产电影网| 国产精品嫩草影院av在线观看| 捣出白浆h1v1| 亚洲综合色网址| 女的被弄到高潮叫床怎么办| 伊人亚洲综合成人网| 国产成人免费无遮挡视频| 午夜福利乱码中文字幕| 日韩制服骚丝袜av| 999精品在线视频| 亚洲av日韩在线播放| xxxhd国产人妻xxx| 蜜桃在线观看..| 精品一区二区免费观看| 热re99久久精品国产66热6| 熟女电影av网| 国产免费视频播放在线视频| 18禁国产床啪视频网站| 一级片'在线观看视频| 美女视频免费永久观看网站| 婷婷色麻豆天堂久久| 亚洲国产欧美网| 国产精品不卡视频一区二区| av网站在线播放免费| 香蕉精品网在线| 在线看a的网站| 亚洲精品久久午夜乱码| 免费看不卡的av| 菩萨蛮人人尽说江南好唐韦庄| 麻豆乱淫一区二区| 久久99一区二区三区| 久久 成人 亚洲| 午夜福利视频在线观看免费| 黄频高清免费视频| 寂寞人妻少妇视频99o| 日韩在线高清观看一区二区三区| 亚洲一码二码三码区别大吗| 一本色道久久久久久精品综合| 久久午夜福利片| av在线app专区| 丝袜美腿诱惑在线| 中国三级夫妇交换| 又黄又粗又硬又大视频| 亚洲精品一区蜜桃| 一级毛片黄色毛片免费观看视频| 国产成人91sexporn| 日韩欧美精品免费久久| tube8黄色片| 国产成人精品无人区| 日本猛色少妇xxxxx猛交久久| 久久久久精品性色| 久久久久久人妻| 啦啦啦视频在线资源免费观看| 国产 精品1| 熟妇人妻不卡中文字幕| 欧美激情极品国产一区二区三区| 成人国产麻豆网| av电影中文网址| 久久久久久久亚洲中文字幕| 边亲边吃奶的免费视频| 国产亚洲最大av| 香蕉精品网在线| 国产成人精品在线电影| 各种免费的搞黄视频| 男人爽女人下面视频在线观看| 久久久精品免费免费高清| 在线天堂最新版资源| 最近最新中文字幕免费大全7| 日本欧美国产在线视频| 亚洲成人av在线免费| 人人妻人人添人人爽欧美一区卜| 两个人免费观看高清视频| 中国三级夫妇交换| 极品少妇高潮喷水抽搐| 人人妻人人澡人人看| 春色校园在线视频观看| 免费看不卡的av| 亚洲综合色惰| 亚洲图色成人| 久久这里只有精品19| 999精品在线视频| 亚洲精品视频女| 热99国产精品久久久久久7| 国产成人精品福利久久| 啦啦啦中文免费视频观看日本| 两个人免费观看高清视频| 成年女人毛片免费观看观看9 | 中文精品一卡2卡3卡4更新| 久久久久久久国产电影| 国产一区二区三区av在线| av女优亚洲男人天堂| 最近中文字幕2019免费版| xxxhd国产人妻xxx| 99re6热这里在线精品视频| 男女高潮啪啪啪动态图| 80岁老熟妇乱子伦牲交| 国产亚洲午夜精品一区二区久久| 久久国产精品男人的天堂亚洲| 人人妻人人澡人人看| 老汉色∧v一级毛片| 日韩一本色道免费dvd| 在现免费观看毛片| 久久精品国产亚洲av高清一级| 国产精品久久久久久av不卡| 看十八女毛片水多多多| 国产av国产精品国产| 国产毛片在线视频| 久久热在线av| 亚洲欧美成人精品一区二区| av天堂久久9| 两个人免费观看高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲图色成人| 极品人妻少妇av视频| av女优亚洲男人天堂| 狠狠精品人妻久久久久久综合| 人妻少妇偷人精品九色| 久久精品国产亚洲av天美| 国产黄频视频在线观看| 精品少妇一区二区三区视频日本电影 | 涩涩av久久男人的天堂| 最近最新中文字幕免费大全7| 久久久久久久国产电影| 日韩精品有码人妻一区| 韩国精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一二三区在线看| 天天躁夜夜躁狠狠久久av| 国产精品熟女久久久久浪| 香蕉丝袜av| 99九九在线精品视频| 看免费成人av毛片| 久久青草综合色| 亚洲精品美女久久av网站| 天堂8中文在线网| 国产又色又爽无遮挡免| 91久久精品国产一区二区三区| 亚洲精品久久成人aⅴ小说| 久久久久视频综合| 国产色婷婷99| 一本—道久久a久久精品蜜桃钙片| 三上悠亚av全集在线观看| 欧美人与性动交α欧美精品济南到 | 国产成人精品一,二区| 日韩大片免费观看网站| 一个人免费看片子| 中国三级夫妇交换| 日本欧美视频一区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久av美女十八| 亚洲欧洲日产国产| 18+在线观看网站| 国产亚洲午夜精品一区二区久久| 一区福利在线观看| 亚洲综合色网址| 午夜影院在线不卡| 国产激情久久老熟女| 精品一品国产午夜福利视频| 男女高潮啪啪啪动态图| 18+在线观看网站| 国产一区二区在线观看av| 免费黄色在线免费观看| 久久鲁丝午夜福利片| 中文精品一卡2卡3卡4更新| 免费少妇av软件| 一本色道久久久久久精品综合| 国产野战对白在线观看| 最近最新中文字幕免费大全7| 国产日韩欧美亚洲二区| 成年女人毛片免费观看观看9 | 成年人午夜在线观看视频| 黄色一级大片看看| 久久久久人妻精品一区果冻| 国产不卡av网站在线观看| 日韩av免费高清视频| 免费看av在线观看网站| 久久久亚洲精品成人影院| 一级a爱视频在线免费观看| 午夜日本视频在线| 久久久久久免费高清国产稀缺| 免费人妻精品一区二区三区视频| 国产乱来视频区| 国产精品成人在线| 国产乱来视频区| 街头女战士在线观看网站| 久久久精品区二区三区| 这个男人来自地球电影免费观看 | 色婷婷久久久亚洲欧美| 久久精品熟女亚洲av麻豆精品| 搡女人真爽免费视频火全软件| www.自偷自拍.com| 国产精品亚洲av一区麻豆 | 久久久久久久久久人人人人人人| 免费在线观看黄色视频的| 亚洲人成电影观看| 免费观看无遮挡的男女| 欧美精品高潮呻吟av久久| 18在线观看网站| 免费观看性生交大片5| 久久毛片免费看一区二区三区| 久久ye,这里只有精品| 亚洲成av片中文字幕在线观看 | 18禁观看日本| 中文字幕色久视频| 午夜免费观看性视频| √禁漫天堂资源中文www| 亚洲精品一二三| 免费观看性生交大片5| 久久女婷五月综合色啪小说| 亚洲四区av| 麻豆av在线久日| 婷婷成人精品国产| 精品一区二区三区四区五区乱码 | 亚洲一区中文字幕在线| 婷婷色麻豆天堂久久| 香蕉精品网在线| 日本vs欧美在线观看视频| 国产精品麻豆人妻色哟哟久久| 午夜久久久在线观看| 国产不卡av网站在线观看| 免费观看性生交大片5| 成人亚洲欧美一区二区av| 美国免费a级毛片| 精品福利永久在线观看| 蜜桃国产av成人99| av免费在线看不卡| 国产午夜精品一二区理论片| 啦啦啦在线免费观看视频4| 叶爱在线成人免费视频播放| 丝袜脚勾引网站| 黄片无遮挡物在线观看| 18禁观看日本| 在线观看三级黄色| 国产成人一区二区在线| 国产精品av久久久久免费| 色婷婷久久久亚洲欧美| 国产精品成人在线| 极品少妇高潮喷水抽搐| 丰满乱子伦码专区| 欧美日本中文国产一区发布| 欧美人与性动交α欧美软件| 国产男女超爽视频在线观看| 一区二区三区乱码不卡18| av网站免费在线观看视频| 麻豆乱淫一区二区| 人人妻人人澡人人看| 巨乳人妻的诱惑在线观看| 国产精品一区二区在线不卡| 国产精品一区二区在线观看99| av不卡在线播放| 一级片免费观看大全| 久热久热在线精品观看| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 日本av手机在线免费观看| 国产精品一区二区在线观看99| a级毛片黄视频| 亚洲人成电影观看| 激情视频va一区二区三区| 久久99蜜桃精品久久| 黄色 视频免费看| 亚洲色图 男人天堂 中文字幕| 人妻人人澡人人爽人人| 日本-黄色视频高清免费观看| 久久久久国产精品人妻一区二区| 少妇被粗大猛烈的视频| 男人舔女人的私密视频| 精品久久蜜臀av无| 久久亚洲国产成人精品v| 少妇猛男粗大的猛烈进出视频| 久久久久久久久久人人人人人人| 久久精品夜色国产| 日韩中字成人| 免费看av在线观看网站| 一级毛片 在线播放| 午夜免费男女啪啪视频观看| 国产免费视频播放在线视频| 大陆偷拍与自拍| 免费黄色在线免费观看| 日韩制服骚丝袜av| 久久精品国产a三级三级三级| 欧美中文综合在线视频| 亚洲国产日韩一区二区| 亚洲av免费高清在线观看| 边亲边吃奶的免费视频| 中文字幕制服av| 日韩一卡2卡3卡4卡2021年| 久久久国产欧美日韩av| 在线观看三级黄色| 捣出白浆h1v1| 免费观看性生交大片5| 亚洲人成77777在线视频| 国产精品免费视频内射| 精品一区二区三区四区五区乱码 | 一区在线观看完整版| 国产亚洲av片在线观看秒播厂| 中国国产av一级| 国产av精品麻豆| 老司机亚洲免费影院| 男人操女人黄网站| 国产一区有黄有色的免费视频| 国产视频首页在线观看| www.av在线官网国产| av女优亚洲男人天堂| 中文字幕人妻熟女乱码| 十八禁高潮呻吟视频| 中文字幕最新亚洲高清| 丝袜美腿诱惑在线| 国产精品国产三级国产专区5o| 婷婷色综合大香蕉| 一区福利在线观看| 男男h啪啪无遮挡| 欧美日本中文国产一区发布| 亚洲精品av麻豆狂野| 免费大片黄手机在线观看| 欧美日韩视频高清一区二区三区二| 在线精品无人区一区二区三| 我要看黄色一级片免费的| 美女大奶头黄色视频| 中文字幕制服av| √禁漫天堂资源中文www| 亚洲成色77777| 国产精品嫩草影院av在线观看| 欧美人与性动交α欧美精品济南到 | 美女午夜性视频免费| xxxhd国产人妻xxx| 满18在线观看网站| 亚洲精品在线美女| 777米奇影视久久| 免费大片黄手机在线观看| 好男人视频免费观看在线| 国产精品国产三级专区第一集| 午夜福利视频在线观看免费| 黄色怎么调成土黄色| 老司机影院毛片| 一区二区三区激情视频| 国产一区亚洲一区在线观看| 99re6热这里在线精品视频| 日韩欧美精品免费久久| 欧美精品一区二区免费开放| 美女午夜性视频免费| 黄色毛片三级朝国网站| 免费观看在线日韩| 国产精品蜜桃在线观看| 亚洲精品国产一区二区精华液| av视频免费观看在线观看| 日韩欧美精品免费久久| 亚洲色图 男人天堂 中文字幕| 国产精品香港三级国产av潘金莲 | 日本色播在线视频| 青春草视频在线免费观看| 国产精品久久久久久av不卡| 亚洲人成电影观看| 韩国精品一区二区三区| 久久热在线av| 亚洲一码二码三码区别大吗| 国产爽快片一区二区三区| 中国三级夫妇交换| 看免费成人av毛片| 侵犯人妻中文字幕一二三四区| 色婷婷av一区二区三区视频| 久久精品aⅴ一区二区三区四区 | 91精品三级在线观看| 久久韩国三级中文字幕| 亚洲精品,欧美精品| 精品国产一区二区三区四区第35| av卡一久久| 国产女主播在线喷水免费视频网站| 老熟女久久久| 精品第一国产精品| 不卡视频在线观看欧美| 久久精品夜色国产| 可以免费在线观看a视频的电影网站 | 亚洲,欧美精品.| 夜夜骑夜夜射夜夜干| av网站免费在线观看视频| 精品人妻一区二区三区麻豆| 精品少妇久久久久久888优播| 伊人亚洲综合成人网| 欧美成人精品欧美一级黄| 伦理电影免费视频| 男女免费视频国产| 成年av动漫网址| 亚洲内射少妇av| 精品国产露脸久久av麻豆| 中文字幕制服av|