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

    基于圓形斷面的隧道溫度場有限差分計算模型

    2016-12-13 06:34:35周小涵曾艷華魏英杰
    隧道建設(中英文) 2016年11期
    關鍵詞:寒區(qū)風流圓形

    周小涵, 曾艷華, 白 赟, 魏英杰, 何 省

    (1. 西南交通大學交通隧道工程教育部重點實驗室, 四川 成都 610031; 2. 中鐵大橋勘測設計院集團有限公司, 湖北 武漢 430050; 3. 成都建筑工程集團總公司, 四川 成都 610014;4. 中國中鐵二院工程集團有限責任公司, 四川 成都 610031)

    ?

    基于圓形斷面的隧道溫度場有限差分計算模型

    周小涵1,2, 曾艷華1,*, 白 赟1, 魏英杰3, 何 省4

    (1. 西南交通大學交通隧道工程教育部重點實驗室, 四川 成都 610031; 2. 中鐵大橋勘測設計院集團有限公司, 湖北 武漢 430050; 3. 成都建筑工程集團總公司, 四川 成都 610014;4. 中國中鐵二院工程集團有限責任公司, 四川 成都 610031)

    除了少數(shù)圓形隧道外,大部分山嶺交通隧道斷面采用馬蹄形或端墻式等形狀。在隧道溫度場的預測中,圓形斷面模型能否代替馬蹄形等實際隧道模型,其適應性值得研究。應用基于空氣-襯砌-圍巖的對流-導熱耦合作用控制方程的有限差分方法,建立圓形斷面模型對東北寒區(qū)馬蹄形隧道溫度場進行計算,并與現(xiàn)場實測溫度場進行對比。結(jié)果表明: 1) 圓形斷面隧道模型有限差分計算方法克服了通用有限元軟件建模復雜、對硬件要求高的弊端,考慮了隧道內(nèi)風流速度和入口風流溫度的影響,在隧道溫度場的預測計算中能夠滿足工程使用要求。2)隧道內(nèi)風流速度和入口風流溫度對隧道溫度場影響較大。本文算例中,入口風流溫度每升高10 ℃,二次襯砌表面溫度升高約7.2 ℃,增幅均勻;從1~5 m/s,洞內(nèi)風流速度每增大1 m/s,二次襯砌表面溫度降低的幅度為6.6、2.7、1.5、0.9 ℃,降幅越來越小。

    馬蹄形隧道; 圓形隧道; 溫度場; 有限差分; 現(xiàn)場實測; 風流溫度; 風流速度

    0 引言

    準確預測隧道溫度場在隧道貫通后的變化規(guī)律是寒區(qū)和高地溫隧道設計的前提,國內(nèi)外學者對隧道溫度場展開過多方面的研究。

    在隧道結(jié)構(gòu)的熱傳導方面: R.N.Pande等[1]推導了兩相介質(zhì)的整合熱傳導率取值方法;J.A.Krumhansl[2]討論了有序和無序固體介質(zhì)的熱傳導問題;張耀等[3]考慮正弦曲線規(guī)律變化的對流換熱邊界條件,建立了復合結(jié)構(gòu)的圓形隧道熱傳導方程。

    在理論求解方面: 賴遠明等[4]應用無量綱量和攝動技術,給出了圓形隧道凍結(jié)過程的近似解析解;S.Singh等[5-6]利用疊加原理和分離變量方法得到了圓形斷面隧道在考慮對流換熱邊界條件下的瞬態(tài)溫度場解析解;張國柱等[7]利用疊加原理及貝塞爾特征函數(shù)的正交及展開定理,得到了寒區(qū)隧道圍巖徑向溫度的理論解;J.L.Battaglia等[8]進行了傳熱模型的拉普拉斯變換反演研究;邵珠山等[9]利用無量綱化和微分方程級數(shù)求解法,得到了溫度場、位移場及應力場的熱彈性理論解;鄭陽等[10]通過Fluent與Ansys的假耦合計算,建立了馬蹄形隧道圍巖、襯砌及洞內(nèi)氣體的流固對流換熱模型。

    在保溫隔熱層計算方法方面: K.Okada[11]提出了寒區(qū)隧道隔熱層的一種設計方法;陳建勛[12]推導了圓形隧道模型防凍隔溫層厚度的計算方法;夏才初等[13]采用理論分析和數(shù)值模擬計算等手段,研究了圓形隧道不同類型凍土段隔熱(保溫)層鋪設厚度計算方法。

    在溫度實測方面: 陳建勛等[14]對隧道拱頂、拱腰邊墻和路面4個部位的溫度進行1.5年的測試,總結(jié)了其溫度變化規(guī)律;郝飛等[15]利用Ansys軟件對馬蹄形寒區(qū)隧道進行了模擬分析,將計算的理論值與實測值進行了對比。

    除了少數(shù)圓形隧道外,大部分公路和鐵路隧道斷面采用的是馬蹄形或端墻式等形狀,圓形斷面模型能否代替馬蹄形等實際隧道模型,其適應性值得研究。在當前的實際應用中,若使用通用計算軟件(Fluent等),可采用隧道實際斷面模型;但在實際建模中,由于該類軟件無法同時考慮對流換熱、圍巖導熱和空氣傳熱等問題,隧道圍巖原始溫度等邊界條件設置困難,模型網(wǎng)格量大,對電腦硬件要求高,工程適應性較差。

    文獻[16]建立了基于空氣-襯砌-圍巖的對流-導熱耦合控制方程的有限差分計算模型,本文旨在通過研究以圓形斷面代替馬蹄形斷面計算隧道溫度場的適應性,以驗證基于空氣-襯砌-圍巖的對流-導熱耦合控制方程的有限差分計算模型在隧道溫度場預測中的可行性,同時重點分析隧道內(nèi)風流速度和隧道入口風流溫度2大因素對隧道溫度場的影響。

    1 有限差分控制方程

    文獻[16-17]在已有研究中基于空氣-襯砌-圍巖的對流-導熱耦合控制方程,用有限差分法實現(xiàn)了計算,同時考慮了隧道入口風流溫度和洞內(nèi)風流速度對隧道溫度場的影響。本文列出圍巖和氣流溫度場主要控制方程,其他方程不再贅述。

    1.1 圓形有限差分模型控制方程

    在流體力學中,把水力半徑相等的圓管直徑定義為非圓管的當量直徑

    R=A/L=0.25d。

    (1)

    式中:A=(π×d2)/4,為截面積;L=π×d,為圓周長度。與圓形斷面相比,非圓形斷面的當量直徑也可以用4倍的水力半徑表示。

    隧道縱向瞬態(tài)溫度場圍巖的控制方程(柱坐標系下)可以表示為

    (2)

    邊界條件

    將隧道內(nèi)氣流視為無黏性、不可壓縮的穩(wěn)定流體,根據(jù)隧道內(nèi)縱向氣流的質(zhì)量守恒、動量守恒和能量守恒定律,得到隧道縱向瞬態(tài)溫度場氣流部分的控制方程

    (5)

    邊界條件為

    tf(x,τ)=tb(τ=0)。

    (6)

    式(2)—(6)中: t是襯砌或圍巖內(nèi)部任一點的溫度,℃; r是襯砌或圍巖中任一點到隧道中心的距離,m; τ是時間參數(shù),s; ρ是空氣的密度,kg/m3; a為圍巖的導溫系數(shù); cp是空氣的比熱容,J/(kg·℃); tf是風流溫度,℃; tb是隧道壁面溫度,℃; tw是圍巖溫度,℃; h是對流換熱系數(shù),W/(m2·℃); U是隧道斷面周長,m; qs是隧道內(nèi)熱源的發(fā)熱量,W; v是風流速度,m/s。

    1.2 有限差分相關參數(shù)

    為了在實際隧道溫度場計算中應用,必須根據(jù)隧道的實際結(jié)構(gòu)對其空間域進行離散處理。

    沿徑向,將二次襯砌、初期支護分成5等分,并將圍巖劃分為若干個子區(qū)域,每個子區(qū)域間距為1m,節(jié)點1為對流邊界,節(jié)點6和節(jié)點11為復合材料傳熱邊界,其余點均為內(nèi)部點。計算時沿縱向每30m設置一個節(jié)點。橫斷面內(nèi)差分節(jié)點劃分如圖1所示。縱斷面內(nèi)差分節(jié)點劃分如圖2所示。

    圖1 橫斷面內(nèi)差分節(jié)點劃分

    Fig. 1 Partition of primary support, secondary lining and surrounding rock

    圖2 縱斷面內(nèi)差分節(jié)點劃分

    用偏微分方程的替代法建立隧道凈空與襯砌之間、襯砌內(nèi)部、襯砌與圍巖之間、圍巖內(nèi)部的傳熱及隧道氣流傳熱的中心有限差分方程,結(jié)合能量守恒定律建立三維瞬態(tài)傳熱的異步長顯示差分格式。利用差分格式穩(wěn)定性判據(jù)式,確定時間步長Δτ與隧道縱向空間步長Δx,編制程序進行計算。

    2 寒區(qū)隧道溫度實測

    2015年1月,對東北牡綏線鐵路改造工程在建長度為6 170 m的綏陽隧道(于2013年9月貫通)進行溫度實測。根據(jù)隧道長度在綏陽隧道設置了13個測試斷面,如圖3所示。隧道二次襯砌表面溫度測試點由拱頂?shù)杰壍烂嬉来尾荚O,如圖4所示。使用紅外線測溫儀GM1150對隧道斷面二次襯砌表面的4個測點進行溫度實測和記錄分析。使用機械風表對隧道內(nèi)風流速度和風向進行了測試,測試結(jié)果表明,綏陽隧道內(nèi)的自然風流速度為3.8~4.0 m/s,風向都是自小里程端吹向大里程端。

    3 溫度場計算值與實測值對比

    綏陽隧道橫斷面凈空面積約為94 m2,橫斷面周長約為37.5 m,取其等效水力半徑為2.5 m。對綏陽隧道所在地區(qū)從隧道貫通到溫度實測期間的溫度進行統(tǒng)計分析,得出隧道入口風流溫度tf(℃)隨時間i(d)的關系

    tf=2.5+18.5×sin (2×π×i /365)。

    (7)

    圖3 綏陽隧道溫度測試斷面隨縱向分布示意圖

    圖4 隧道二次襯砌表面溫度測試點橫斷面分布圖

    使用tf作為隧道入口風流溫度邊界條件進行計算,計算風流速度取隧道內(nèi)實測平均風流速度3.8 m/s,方向由小里程端吹向大里程端。隧道原始巖溫按照地勘報告中埋深和地溫梯度3.0 ℃/100 m選取。根據(jù)風流速度與對流換熱系數(shù)的關系[18],本次計算對流換熱系數(shù)取為21 W/(m2·℃),其他介質(zhì)材料計算參數(shù)取值見表1。由于缺少實驗參數(shù),本文計算不考慮相變問題。

    應用式(2)—(6)的有限差分控制方程,編程計算得出綏陽隧道由2013年9月貫通至2015年1月隧道襯砌表面溫度場分布情況。將綏陽隧道2015年1月7—13日各襯砌表面測點的實測溫度平均值和數(shù)值計算結(jié)果中相應時段的襯砌表面平均溫度進行比較分析,如圖5所示。

    表1 介質(zhì)材料參數(shù)

    圖5 綏陽隧道襯砌表面溫度實測值和圓形斷面模型計算值

    Fig. 5 Comparison between measured lining surface temperature of Suiyang Tunnel and circular cross-section model calculated results

    通過現(xiàn)場溫度測試數(shù)據(jù)可以看出,隧道拱頂二次襯砌表面溫度最高而軌道面測點溫度最低。從拱頂?shù)杰壍烂妫?個測點溫度逐漸降低。各測點溫度平均值與測點2和測點3的溫度接近,測點平均溫度比拱頂測點溫度低約0.7 ℃,比軌道面測點溫度高約0.4 ℃。

    對比有限差分計算結(jié)果和實測溫度結(jié)果可知,兩者隨隧道縱向分布趨勢一致且差值較小。隧道出入口段計算溫度和實測溫度有較大差異,主要是因為洞口處圍巖受外界大氣影響顯著,而該因素在模型計算中考慮不足以及計算中沒有考慮相變所致。在洞身段,數(shù)值計算襯砌平均溫度高于實際測點平均溫度約1.9 ℃,滿足工程使用要求。

    4 風流溫度和風流速度對溫度場的影響

    基于本文有限差分算法,利用綏陽隧道實際模型,重點研究隧道內(nèi)風流溫度和風流速度對溫度場的影響。為了便于研究,分2種情況: 1)保持隧道內(nèi)風流速度3.8 m/s 不變,隧道入口風流溫度采用-20、-10、0、10、20 ℃ 5個恒定值;2)保持入口風流溫度-20 ℃不變,隧道內(nèi)風流速度分別采用1、2、3、4、5 m/s 5個恒定值。選擇隧道橫斷面DK493+529的計算結(jié)果進行分析。不同風流溫度和風流速度條件下隧道二次襯砌表面溫度如圖6和圖7所示。

    圖6 不同風流溫度時綏陽隧道襯砌表面溫度(DK493+529)

    Fig. 6 Lining surface temperature of Suiyang Tunnel under different wind flow temperatures at cross-section DK493+529

    圖7 不同洞內(nèi)風流速度條件下綏陽隧道襯砌表面溫度(DK493+529)

    Fig. 7 Lining surface temperature of Suiyang Tunnel under different wind flow velocities at cross-section DK493+529

    從圖6可知,隧道入口風流溫度對隧道襯砌溫度影響較大。本文算例中,計算時間為30 d時,風流溫度每升高10 ℃,二次襯砌表面溫度升高約7.2 ℃,增幅均勻。

    從圖7可知,隧道洞內(nèi)風流速度對隧道襯砌溫度影響較大。當有低溫風流對隧道進行降溫作用時,隨著風流速度的增大,二次襯砌表面溫度越來越低,但是其降低的幅度越來越小。計算時間為30 d時,從1~5 m/s,洞內(nèi)風流速度每增大1 m/s,二次襯砌表面溫度降低的幅度為6.6、2.7、1.5、0.9 ℃。

    5 結(jié)論與討論

    應用基于空氣-襯砌-圍巖的對流-導熱耦合作用控制方程的有限差分方法,建立圓形斷面模型對東北寒區(qū)馬蹄形隧道溫度場進行計算,并與實測溫度場進行對比,得出如下結(jié)論。

    1)根據(jù)現(xiàn)場溫度測試數(shù)據(jù)可知,隧道拱頂二次襯砌表面的溫度最高而軌道面測點溫度最低。各測點溫度平均值與拱腰的溫度接近。對比寒區(qū)隧道溫度場有限差分計算結(jié)果和實測結(jié)果可知,兩者隨隧道縱向分布趨勢一致且差值較小。圓形斷面模型有限差分計算方法克服了通用有限元軟件無法同時考慮對流換熱和傳熱問題、建模復雜和對硬件要求高等弊端,在隧道溫度場長期分布規(guī)律的預測計算中易于使用。

    2)隧道入口風流溫度和洞內(nèi)風流速度對隧道溫度場有較大影響。本文算例中,入口風流溫度每升高10 ℃,二次襯砌表面溫度升高約7.2 ℃,增幅均勻; 從1~5 m/s,洞內(nèi)風流速度每增大1 m/s,二次襯砌表面溫度降低的幅度為6.6、2.7、1.5、0.9 ℃,降幅越來越小。

    3)在收集隧道當?shù)貧夂蛸Y料而得到預測隧道入口風流溫度的條件下,本文計算模型可以用來計算不同隧道內(nèi)風流速度條件下的溫度場分布,為具體設計提供幫助。本方法對凍土(季節(jié)性凍土和永久性凍土)隧道的溫度場及隔熱層厚度預測的適應性值得進一步研究。

    [1] Pande R N, Kumar V, Chaudhary D R. Thermal conduction in a homogeneous two-phase system[J]. Pramana, 1984, 22(1): 63-70.

    [2] Krumhansl J A. On the general physics of heat conduction in ordered and disordered solids[J]. Springer Berlin Heidelberg, 1993, 112: 3-8.

    [3] 張耀,何樹生,李靖波. 寒區(qū)有隔熱層的圓形隧道溫度場解析解[J]. 冰川凍土, 2009, 31(1): 113-118. (ZHANG Yao,HE Shusheng,LI Jingbo. Analytic solutions for the temperature fields of a circular tunnel with insulation layer in cold region[J]. Journal of Glaciology and Geocryology,2009,31(1): 113-118. (in Chinese))

    [4] 賴遠明,喻文兵,吳紫汪,等. 寒區(qū)圓形截面隧道溫度場的解析解[J]. 冰川凍土, 2001, 23(2): 126-130. (LAI Yuanming, YU Wenbing, WU Ziwang, et al. Approximate analytical solution for the temperature fields of a circular tunnel in cold regions[J]. Journal of Glaciology and Geocryology, 2001, 23(2): 126-130. (in Chinese))

    [5] Singh S,Jain P K,Rizwan-uddin. Analytical solution to transient heat conduction in polar coordinates with multiple layers in radial direction[J]. International Journal of Thermal Sciences,2008,47(3): 261-273.

    [6] Jain P K,Singh S,Rizwan-uddin. Analytical solution to transient asymmetric heat conduction in a multilayer annulus[J]. Journal of Heat Transfer,2009,131(1): 73-92.

    [7] 張國柱,夏才初,殷卓. 寒區(qū)隧道軸向及徑向溫度分布理論解[J]. 同濟大學學報(自然科學版), 2010, 38(8): 1117-1122. (ZHANG Guozhu, XIA Caichu, YIN Zhuo. Analytical solution to axial and radial temperature of tunnel in cold region[J]. Journal of Tongji University (Natural Science), 2010, 38(8): 1117-1122. (in Chinese))

    [8] Battaglia J L, Toutain J, Pradere C, et al. Numerical inversion of Laplace transform for time resolved thermal characterization experiment[J]. Journal of Heat Transfer, 2011, 133 (4): 626-638.

    [9] 邵珠山,喬汝佳,王新宇. 高地溫隧道溫度與熱應力場的彈性理論解[J]. 巖土力學, 2013(增刊1): 1-8. (SHAO Zhushan,QIAO Rujia,WANG Xinyu. Elasticity solution for temperature and stress fields of tunnels with high geothermal temperature[J]. Rock and Soil Mechanics, 2013(S1): 1-8. (in Chinese))

    [10] 鄭陽,王道遠,朱永全. 寒區(qū)運營隧道溫度場數(shù)值模擬[J]. 公路工程, 2014, 39(1): 40-43. (ZHENG Yang, WANG Daoyuan, ZHU Yongquan. Numerical simulation of temperature field in cold region operational tunnel[J]. Highway Engineering, 2014, 39(1): 40-43. (in Chinese))

    [11] Okada K. Prevention by adiatic treatment of tunnel lining[J]. Japanese Railway Engineering, 1985, 26 (3): 75-80.

    [12] 陳建勛. 寒冷地區(qū)隧道防凍隔溫層設計計算方法及應用[J]. 土木工程學報, 2004, 37(11): 85-88. (CHEN Jianxun. Design method and application for anti-freezing thermal-protective layers in cold-area tunnels [J]. China Civil Engineering Journal, 2004, 37(11): 85-88. (in Chinese))

    [13] 夏才初,范東方,韓常領. 寒區(qū)隧道不同類型凍土段隔熱(保溫)層鋪設厚度計算方法[J]. 中國公路學報, 2013, 26(5): 131-139. (XIA Caichu, FAN Dongfang, HAN Changling. Piecewise calculation method for insulation layer thickness in cold region tunnels[J]. China Journal of Highway and Transport, 2013, 26(5): 131-139. (in Chinese))

    [14] 陳建勛,羅彥斌. 寒冷地區(qū)隧道溫度場的變化規(guī)律[J]. 交通運輸工程學報, 2008, 8(2): 44-48. (CHEN Jianxun, LUO Yanbin. Changing rules of temperature field for tunnel in cold area [J]. Journal of Traffic and Transportation Engineering,2008,8(2):44-48. (in Chinese))

    [15] 郝飛,孫全勝. 寒區(qū)公路隧道溫度場及保溫層的研究[J]. 現(xiàn)代隧道技術, 2012, 49(1): 39-43. (HAO Fei, SUN Quansheng. Research on the temperature fields and insulating layers of highway tunnels in a cold region [J]. Modern Tunnelling Technology, 2012,49(1): 39-43. (in Chinese))

    [16] 周小涵. 通風對寒區(qū)隧道抗防凍的影響及合理設防范圍研究[D]. 成都: 西南交通大學, 2012: 29-31. (ZHOU Xiaohan. Study of the influence of ventilation on frost resistance and reasonable range of resistance for tunnel in cold region[D]. Chengdu: Southwest Jiaotong University, 2012: 29-31. (in Chinese))

    [17] ZHOU Xiaohan, ZENG Yanhua, FAN Lei. Temperature field analysis of a cold-region railway tunnel considering mechanical and train-induced ventilation effects [J]. Applied Thermal Engineering, 2016, 100: 114-124.

    [18] 張建榮,劉照球. 混凝土對流換熱系數(shù)的風洞實驗研究[J]. 土木工程學報, 2006, 39(9): 39-42. (ZHANG Jianrong, LIU Zhaoqiu. A study of the convective heat transfer coefficient of concrete in wind tunnel experiment [J]. China Civil Engineering Journal, 2006,39(9): 39-42. (in Chinese))

    Finite Difference Calculation Model for Tunnel Temperature Field Based on Circular Cross-section

    ZHOU Xiaohan1, 2, ZENG Yanhua1,*, BAI Yun1, WEI Yingjie3, HE Xing4

    (1.KeyLaboratoryofTransportationTunnelEngineeringofMinistryofEducation,SouthwestJiaotongUniversity,Chengdu610031,Sichuan,China; 2.ChinaRailwayMajorBridgeReconnaissance&DesignInstituteCo.,Ltd.,Wuhan430050,Hubei,China; 3.ChengduConstructionEngineeringGroupCorporation,Chengdu610014,Sichuan,China; 4.ChinaRaiwayEryuanEngineeringGroupCo.,Ltd.,Chengdu610031,Sichuan,China)

    The horseshoe-shaped tunnels and end wall type tunnels are common in mountain areas. It’s necessary to study the suitability of tunnel temperature field calculation model using circular cross-section. The temperature fields of cold region horseshoe-shaped tunnels in Northeast China are calculated by circular cross-section model using finite difference method of convection-conduction equations based on air-lining-surrounding rock. The calculation results are compared with measured data. The results show that: 1) The finite difference calculation method by using circular cross-section model and considering effects of wind flow velocity and wind flow temperature, is easy and convenient; it can meet the construction requirements. 2) The wind flow velocity in tunnel and wind flow field at tunnel portal have a significant influence on tunnel temperature field. The temperature of secondary lining surface increases by 7.2 ℃ when wind flow temperature at tunnel portal increases by 10 ℃ . The temperature of secondary lining surface reduces by 6.6 ℃, 2.7 ℃, 1.5 ℃ and 0.9 ℃ when the wind flow velocity is 2 m/s, 3 m/s, 4 m/s and 5 m/s respectively.

    horseshoe-shaped tunnel; circular tunnel; temperature field; finite difference; field monitoring; wind flow temperature; wind flow velocity

    2016-02-03;

    2016-10-31

    國家自然科學基金資助項目(51278426)

    周小涵(1988—),男,重慶忠縣人,西南交通大學橋梁與隧道工程專業(yè)在讀博士,從事大型交通隧道的防災救援技術和地鐵結(jié)構(gòu)研究。E-mail: zhouxh2008@126.com。 *通訊作者: 曾艷華, E-mail: zengyhua@163.com。

    10.3973/j.issn.1672-741X.2016.11.007

    U 45

    A

    1672-741X(2016)11-1332-05

    猜你喜歡
    寒區(qū)風流圓形
    數(shù)你最風流
    黃河之聲(2021年14期)2021-10-28 03:53:04
    不同施肥處理對寒區(qū)苜蓿草地產(chǎn)量的影響
    留白——不著點墨,盡顯風流
    兼職者的風流史
    風流總被雨打風吹去
    海峽姐妹(2019年7期)2019-07-26 00:50:50
    為什么窨井蓋大多都是圓形的
    肥皂泡為什么是圓形?
    圓形題
    圓形變身喵星人
    寒區(qū)某水電站大管徑PCCP管施工過程
    精品国产乱码久久久久久男人| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 欧美一区二区精品小视频在线| 久久久久国产精品人妻aⅴ院| 午夜精品久久久久久毛片777| 久久香蕉精品热| www.精华液| 99热6这里只有精品| 88av欧美| 亚洲成人精品中文字幕电影| 妹子高潮喷水视频| 久久精品国产综合久久久| 成人亚洲精品av一区二区| 1024手机看黄色片| 日韩高清综合在线| 午夜亚洲福利在线播放| 黄色女人牲交| 国产又色又爽无遮挡免费看| 午夜福利免费观看在线| 麻豆成人午夜福利视频| 国产人伦9x9x在线观看| 国产亚洲精品一区二区www| 小说图片视频综合网站| 国产成人精品久久二区二区91| 亚洲人成伊人成综合网2020| 91麻豆av在线| 黑人欧美特级aaaaaa片| 亚洲av成人一区二区三| 国产成人aa在线观看| 动漫黄色视频在线观看| 精品一区二区三区av网在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品 国内视频| 好男人在线观看高清免费视频| 成熟少妇高潮喷水视频| 国产熟女xx| 少妇熟女aⅴ在线视频| 精品国产超薄肉色丝袜足j| 日本熟妇午夜| 久久久久精品国产欧美久久久| 精品久久久久久久末码| 国产精品av久久久久免费| 毛片女人毛片| 两个人看的免费小视频| 国产亚洲精品久久久久久毛片| 97人妻精品一区二区三区麻豆| 19禁男女啪啪无遮挡网站| 啪啪无遮挡十八禁网站| 人妻久久中文字幕网| 色综合婷婷激情| 99久久久亚洲精品蜜臀av| 免费av毛片视频| 69av精品久久久久久| 亚洲欧美精品综合久久99| 欧美 亚洲 国产 日韩一| 神马国产精品三级电影在线观看 | 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 久久婷婷成人综合色麻豆| 深夜精品福利| 午夜福利视频1000在线观看| 日本在线视频免费播放| 欧美激情久久久久久爽电影| 免费看a级黄色片| 亚洲av成人av| 麻豆成人午夜福利视频| 国产99久久九九免费精品| 亚洲av五月六月丁香网| 国产男靠女视频免费网站| 一二三四社区在线视频社区8| 久久精品91无色码中文字幕| www.精华液| 老司机深夜福利视频在线观看| 好男人电影高清在线观看| 黄片小视频在线播放| 亚洲aⅴ乱码一区二区在线播放 | 99re在线观看精品视频| 日本a在线网址| 国产精品久久久久久人妻精品电影| 别揉我奶头~嗯~啊~动态视频| svipshipincom国产片| 久久精品亚洲精品国产色婷小说| 日韩精品青青久久久久久| 五月玫瑰六月丁香| 成人手机av| 级片在线观看| 久热爱精品视频在线9| 久久久国产欧美日韩av| 国内毛片毛片毛片毛片毛片| www国产在线视频色| 国产熟女xx| 久久欧美精品欧美久久欧美| 黄色视频不卡| 日本一二三区视频观看| 夜夜躁狠狠躁天天躁| 亚洲专区国产一区二区| 免费av毛片视频| 日韩高清综合在线| 黄色a级毛片大全视频| 又紧又爽又黄一区二区| 国产成人系列免费观看| 国产欧美日韩一区二区三| 99久久精品热视频| 99热只有精品国产| 亚洲精品色激情综合| 中文字幕熟女人妻在线| 国产熟女午夜一区二区三区| 色综合欧美亚洲国产小说| 男女那种视频在线观看| 亚洲男人的天堂狠狠| 国产三级中文精品| 一区二区三区激情视频| 亚洲国产精品合色在线| 中文在线观看免费www的网站 | 成人欧美大片| 正在播放国产对白刺激| 精品久久久久久久久久久久久| 最近最新中文字幕大全电影3| 最新在线观看一区二区三区| 欧美一级毛片孕妇| 精品无人区乱码1区二区| 99久久国产精品久久久| 久久久久久亚洲精品国产蜜桃av| 母亲3免费完整高清在线观看| www国产在线视频色| 又黄又粗又硬又大视频| 2021天堂中文幕一二区在线观| 女生性感内裤真人,穿戴方法视频| 男插女下体视频免费在线播放| 亚洲男人的天堂狠狠| 色老头精品视频在线观看| 亚洲真实伦在线观看| 琪琪午夜伦伦电影理论片6080| 大型av网站在线播放| 久久久久久亚洲精品国产蜜桃av| 美女 人体艺术 gogo| 一级毛片女人18水好多| 午夜福利高清视频| 亚洲精品国产精品久久久不卡| 日本一二三区视频观看| 成人欧美大片| av天堂在线播放| 日本 av在线| 一进一出好大好爽视频| 亚洲国产欧洲综合997久久,| av有码第一页| 极品教师在线免费播放| 欧美成人性av电影在线观看| 在线观看日韩欧美| 国产亚洲欧美98| 色尼玛亚洲综合影院| 国产久久久一区二区三区| 久久婷婷人人爽人人干人人爱| 在线永久观看黄色视频| 免费电影在线观看免费观看| 91成年电影在线观看| 黄色视频不卡| av有码第一页| 国产三级黄色录像| 国产精品九九99| 国产一区二区在线av高清观看| 亚洲狠狠婷婷综合久久图片| 久久精品国产亚洲av高清一级| 岛国在线免费视频观看| 麻豆av在线久日| 欧美在线一区亚洲| 制服诱惑二区| 淫妇啪啪啪对白视频| 18美女黄网站色大片免费观看| 国产av不卡久久| 亚洲人与动物交配视频| 国产1区2区3区精品| 国产av麻豆久久久久久久| 狂野欧美白嫩少妇大欣赏| 欧美最黄视频在线播放免费| 色老头精品视频在线观看| 亚洲成人中文字幕在线播放| а√天堂www在线а√下载| 蜜桃久久精品国产亚洲av| 欧美绝顶高潮抽搐喷水| 日韩 欧美 亚洲 中文字幕| 成人午夜高清在线视频| 欧美黄色淫秽网站| 可以在线观看毛片的网站| 亚洲全国av大片| 欧美性猛交黑人性爽| 三级毛片av免费| 亚洲精品av麻豆狂野| 国产午夜福利久久久久久| 久久这里只有精品中国| 国产精品日韩av在线免费观看| 国产黄a三级三级三级人| 两个人的视频大全免费| 国产午夜精品论理片| 一区二区三区激情视频| 一级a爱片免费观看的视频| 琪琪午夜伦伦电影理论片6080| 一个人观看的视频www高清免费观看 | 性欧美人与动物交配| 亚洲中文日韩欧美视频| 国产精品一区二区精品视频观看| xxx96com| 欧美日本视频| 国产99白浆流出| 美女大奶头视频| 麻豆国产97在线/欧美 | 久久这里只有精品19| 精品国产乱子伦一区二区三区| 男女床上黄色一级片免费看| 女生性感内裤真人,穿戴方法视频| 99在线视频只有这里精品首页| 国产伦一二天堂av在线观看| 黄色毛片三级朝国网站| 国产成人系列免费观看| 三级毛片av免费| av福利片在线观看| 性色av乱码一区二区三区2| 91麻豆av在线| 亚洲av五月六月丁香网| 国产三级中文精品| 久久久久久久久久黄片| 久久精品成人免费网站| www国产在线视频色| 特级一级黄色大片| 国产在线观看jvid| 久久精品综合一区二区三区| 日本一二三区视频观看| 性欧美人与动物交配| 麻豆一二三区av精品| 小说图片视频综合网站| 亚洲激情在线av| 色播亚洲综合网| 老司机福利观看| 精品高清国产在线一区| 国产一级毛片七仙女欲春2| 久久99热这里只有精品18| 俄罗斯特黄特色一大片| 91国产中文字幕| 成年女人毛片免费观看观看9| 亚洲午夜理论影院| 男人舔女人的私密视频| 极品教师在线免费播放| 国内精品久久久久久久电影| 欧美高清成人免费视频www| 久久久久久九九精品二区国产 | 国产视频内射| 天堂影院成人在线观看| 99久久99久久久精品蜜桃| 亚洲精品中文字幕在线视频| 国产野战对白在线观看| videosex国产| 别揉我奶头~嗯~啊~动态视频| 色综合婷婷激情| 欧美日本亚洲视频在线播放| 国产精品九九99| 老司机在亚洲福利影院| 亚洲中文日韩欧美视频| 日韩精品免费视频一区二区三区| a级毛片在线看网站| 久久人妻av系列| 日本成人三级电影网站| 日韩成人在线观看一区二区三区| 精品久久久久久久毛片微露脸| 真人一进一出gif抽搐免费| 久久久久久久久免费视频了| 12—13女人毛片做爰片一| 人成视频在线观看免费观看| 亚洲av成人一区二区三| 久久精品综合一区二区三区| 脱女人内裤的视频| 亚洲18禁久久av| 美女扒开内裤让男人捅视频| 亚洲自拍偷在线| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| av视频在线观看入口| 在线国产一区二区在线| 久久精品综合一区二区三区| 最新美女视频免费是黄的| 久久久国产欧美日韩av| 久久久久性生活片| 国产精品亚洲一级av第二区| 十八禁网站免费在线| 妹子高潮喷水视频| 一进一出好大好爽视频| 亚洲五月天丁香| 最新在线观看一区二区三区| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久久电影 | 正在播放国产对白刺激| 日韩欧美精品v在线| 国产真实乱freesex| 亚洲欧美精品综合一区二区三区| 亚洲av第一区精品v没综合| 一级毛片精品| 亚洲精品粉嫩美女一区| 亚洲 欧美一区二区三区| 婷婷六月久久综合丁香| 亚洲精品国产一区二区精华液| 午夜福利在线在线| 国产麻豆成人av免费视频| 精品国产亚洲在线| 久久性视频一级片| 国产精品亚洲一级av第二区| 久久精品成人免费网站| 亚洲欧美精品综合一区二区三区| 欧美黑人巨大hd| 欧美性长视频在线观看| 丝袜美腿诱惑在线| 在线观看舔阴道视频| 国产黄a三级三级三级人| 午夜久久久久精精品| 日韩 欧美 亚洲 中文字幕| xxxwww97欧美| 97碰自拍视频| 日本一区二区免费在线视频| 男人舔女人的私密视频| 日本一区二区免费在线视频| 一进一出抽搐gif免费好疼| 99riav亚洲国产免费| 亚洲色图av天堂| 在线观看美女被高潮喷水网站 | 亚洲成人中文字幕在线播放| 精品久久蜜臀av无| 亚洲欧美日韩无卡精品| www国产在线视频色| 国产av在哪里看| 日本一区二区免费在线视频| 又粗又爽又猛毛片免费看| 午夜精品久久久久久毛片777| 亚洲中文字幕日韩| 亚洲五月婷婷丁香| 999久久久精品免费观看国产| 丝袜人妻中文字幕| 精品少妇一区二区三区视频日本电影| 亚洲精华国产精华精| 曰老女人黄片| 日韩av在线大香蕉| 亚洲午夜精品一区,二区,三区| 性色av乱码一区二区三区2| 国产亚洲精品一区二区www| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩高清在线视频| 变态另类丝袜制服| 天堂av国产一区二区熟女人妻 | 成人18禁高潮啪啪吃奶动态图| 国产蜜桃级精品一区二区三区| 欧美中文综合在线视频| 国产精品久久久久久精品电影| 欧美成人免费av一区二区三区| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| 亚洲 国产 在线| 欧美人与性动交α欧美精品济南到| 亚洲 欧美 日韩 在线 免费| 欧美成狂野欧美在线观看| 久热爱精品视频在线9| 色精品久久人妻99蜜桃| 中文亚洲av片在线观看爽| 一边摸一边抽搐一进一小说| 在线观看免费视频日本深夜| 久久中文字幕人妻熟女| 91av网站免费观看| 欧美国产日韩亚洲一区| 在线观看66精品国产| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| av福利片在线观看| 90打野战视频偷拍视频| 天堂影院成人在线观看| 在线视频色国产色| 欧美色欧美亚洲另类二区| 嫩草影视91久久| 国产成人av教育| 好男人电影高清在线观看| 国产成人影院久久av| 香蕉久久夜色| 欧美黑人精品巨大| tocl精华| 国产精品久久久久久久电影 | 99国产精品一区二区蜜桃av| 精品福利观看| 一边摸一边抽搐一进一小说| 熟女电影av网| 夜夜爽天天搞| 看片在线看免费视频| 亚洲精品久久成人aⅴ小说| aaaaa片日本免费| 国产精品综合久久久久久久免费| 欧美黑人巨大hd| 久久久久免费精品人妻一区二区| 精品国产美女av久久久久小说| 又爽又黄无遮挡网站| 黄色 视频免费看| 岛国在线免费视频观看| 亚洲人成网站高清观看| 国产一区在线观看成人免费| 男插女下体视频免费在线播放| 国产精品亚洲美女久久久| 国产91精品成人一区二区三区| 很黄的视频免费| 国产精品乱码一区二三区的特点| 日韩三级视频一区二区三区| 国产精品 欧美亚洲| 久久欧美精品欧美久久欧美| 青草久久国产| 激情在线观看视频在线高清| av在线天堂中文字幕| 91麻豆av在线| 波多野结衣高清无吗| 成人手机av| 欧美日韩精品网址| 国产精品av久久久久免费| 可以在线观看毛片的网站| 精品乱码久久久久久99久播| 欧美日韩亚洲国产一区二区在线观看| 一级黄色大片毛片| 欧美日韩瑟瑟在线播放| 少妇被粗大的猛进出69影院| 99热6这里只有精品| 亚洲国产高清在线一区二区三| 丝袜美腿诱惑在线| 成人欧美大片| 丰满的人妻完整版| 欧美一级毛片孕妇| 亚洲国产精品999在线| 成人av在线播放网站| 国产精品一区二区三区四区免费观看 | 女人爽到高潮嗷嗷叫在线视频| 天天添夜夜摸| 老熟妇乱子伦视频在线观看| 18美女黄网站色大片免费观看| 亚洲精品在线美女| 一进一出抽搐动态| 最好的美女福利视频网| 亚洲男人的天堂狠狠| 国产91精品成人一区二区三区| 一级毛片女人18水好多| 亚洲九九香蕉| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频| av天堂在线播放| 91麻豆av在线| 精品乱码久久久久久99久播| 国产成人精品久久二区二区免费| 18禁国产床啪视频网站| 欧美中文综合在线视频| 日韩国内少妇激情av| 老熟妇乱子伦视频在线观看| 九色成人免费人妻av| 久久久久久大精品| √禁漫天堂资源中文www| 日韩欧美三级三区| 亚洲国产欧美网| 国产高清有码在线观看视频 | 美女高潮喷水抽搐中文字幕| 一级毛片高清免费大全| 岛国在线免费视频观看| 在线观看www视频免费| 精品熟女少妇八av免费久了| 久久人人精品亚洲av| www.精华液| 五月玫瑰六月丁香| 搡老熟女国产l中国老女人| 国产精品久久久久久人妻精品电影| 亚洲成人免费电影在线观看| 美女 人体艺术 gogo| 欧美在线黄色| 久久人妻av系列| 欧美性猛交黑人性爽| 少妇熟女aⅴ在线视频| a在线观看视频网站| 亚洲 国产 在线| 国产v大片淫在线免费观看| 久久99热这里只有精品18| 亚洲 欧美一区二区三区| 久久人妻福利社区极品人妻图片| 在线观看66精品国产| 成人av在线播放网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产又黄又爽又无遮挡在线| 精品久久久久久久毛片微露脸| 婷婷亚洲欧美| 制服诱惑二区| 亚洲国产精品999在线| 日韩欧美国产在线观看| a级毛片a级免费在线| 成人高潮视频无遮挡免费网站| 黄片大片在线免费观看| 亚洲真实伦在线观看| svipshipincom国产片| 19禁男女啪啪无遮挡网站| 一本大道久久a久久精品| 丝袜人妻中文字幕| 久久久久亚洲av毛片大全| 法律面前人人平等表现在哪些方面| 波多野结衣高清作品| 国产伦在线观看视频一区| 国模一区二区三区四区视频 | 欧美一区二区国产精品久久精品 | av在线天堂中文字幕| 国产精品免费视频内射| 无人区码免费观看不卡| 亚洲国产看品久久| 精品一区二区三区av网在线观看| 欧美日韩亚洲综合一区二区三区_| 三级男女做爰猛烈吃奶摸视频| 国产真实乱freesex| 18禁黄网站禁片免费观看直播| 久久热在线av| 日韩欧美免费精品| 动漫黄色视频在线观看| 女人被狂操c到高潮| 丰满人妻一区二区三区视频av | 性欧美人与动物交配| 宅男免费午夜| 女人爽到高潮嗷嗷叫在线视频| 又大又爽又粗| 日韩精品青青久久久久久| 国产久久久一区二区三区| 亚洲中文av在线| 亚洲午夜精品一区,二区,三区| 欧美色视频一区免费| 五月伊人婷婷丁香| 啦啦啦韩国在线观看视频| 制服丝袜大香蕉在线| 精品国产亚洲在线| 成年女人毛片免费观看观看9| 岛国视频午夜一区免费看| 91在线观看av| 黄片大片在线免费观看| 村上凉子中文字幕在线| 日本在线视频免费播放| 黄片小视频在线播放| 51午夜福利影视在线观看| 天天一区二区日本电影三级| 麻豆国产av国片精品| www日本黄色视频网| 脱女人内裤的视频| 亚洲国产精品成人综合色| 国产高清视频在线播放一区| 亚洲男人的天堂狠狠| 国产亚洲欧美在线一区二区| 老熟妇乱子伦视频在线观看| 亚洲人成电影免费在线| 久久这里只有精品中国| 亚洲av片天天在线观看| a在线观看视频网站| 午夜福利18| 日本一本二区三区精品| 国内毛片毛片毛片毛片毛片| 欧美日韩精品网址| 久久久久九九精品影院| 少妇粗大呻吟视频| 国内精品久久久久久久电影| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利成人在线免费观看| 男女下面进入的视频免费午夜| 色av中文字幕| 一进一出好大好爽视频| 91老司机精品| 99在线人妻在线中文字幕| 国产成人精品久久二区二区免费| 国内毛片毛片毛片毛片毛片| 亚洲午夜精品一区,二区,三区| 男人的好看免费观看在线视频 | 丰满人妻熟妇乱又伦精品不卡| 中文字幕熟女人妻在线| 亚洲电影在线观看av| 久久久国产成人精品二区| 国产成年人精品一区二区| 99国产精品一区二区三区| 亚洲人成电影免费在线| 色综合亚洲欧美另类图片| 国产午夜精品论理片| 真人一进一出gif抽搐免费| 色噜噜av男人的天堂激情| 男男h啪啪无遮挡| 十八禁人妻一区二区| 亚洲美女黄片视频| 国产探花在线观看一区二区| 69av精品久久久久久| 此物有八面人人有两片| 天堂影院成人在线观看| 搡老妇女老女人老熟妇| 成人av一区二区三区在线看| 丰满人妻熟妇乱又伦精品不卡| 国内毛片毛片毛片毛片毛片| 免费在线观看影片大全网站| 日本在线视频免费播放| 日韩欧美免费精品| ponron亚洲| 一本精品99久久精品77| 日日干狠狠操夜夜爽| 日韩精品青青久久久久久| 欧美久久黑人一区二区| 大型av网站在线播放| 人人妻,人人澡人人爽秒播| 淫妇啪啪啪对白视频| 一个人观看的视频www高清免费观看 | 曰老女人黄片| 老司机午夜十八禁免费视频| 亚洲国产欧洲综合997久久,| 国语自产精品视频在线第100页| 精品无人区乱码1区二区| 亚洲18禁久久av| 国产精品98久久久久久宅男小说| 亚洲免费av在线视频| 丰满人妻一区二区三区视频av | 法律面前人人平等表现在哪些方面| 日韩精品免费视频一区二区三区| 亚洲天堂国产精品一区在线| 国产亚洲欧美在线一区二区| 日本熟妇午夜|