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

    單管凍結(jié)溫度場演變規(guī)律分析及試驗驗證

    2019-10-20 00:43:21馬蕓夏錦紅桂超
    人民長江 2019年2期

    馬蕓 夏錦紅 桂超

    摘要:為了研究土中水相變(潛熱)熱傳導(dǎo)過程中凍結(jié)溫度場演變規(guī)律,構(gòu)建了穩(wěn)定線性冷源作用下單管凍結(jié)相變溫度場模型,由變量替換法求得凍結(jié)溫度場隨凍結(jié)時間演變規(guī)律近似解析解,同時利用指數(shù)積分函數(shù)解析式表達(dá)了凍結(jié)區(qū)和未凍結(jié)區(qū)的溫度場分布,并將理論解析解與凍結(jié)槽試驗結(jié)果進(jìn)行比對。研究表明:理論與試驗結(jié)果誤差在合理范圍內(nèi),驗證了瞬態(tài)凍結(jié)溫度場分布近似解析解的有效可靠性;單管凍結(jié)鋒面半徑和凍結(jié)時間有著平方根關(guān)系;凍結(jié)溫度場在時空維度上服從對數(shù)分布,且曲線曲率隨凍結(jié)時間的增加而逐漸減小。研究成果對認(rèn)識穩(wěn)定線性冷源作用下凍結(jié)管周圍不良地基土溫度場演變規(guī)律和指導(dǎo)凍結(jié)法施工技術(shù)有重要意義。

    關(guān)鍵詞:凍結(jié)溫度場; 相變; 凍結(jié)鋒面半徑; 理論解析; 凍結(jié)槽試驗

    中圖法分類號: TD265.34文獻(xiàn)標(biāo)志碼: ADOI:10.16232/j.cnki.1001-4179.2019.02.030

    目前,限于對凍融過程、凍土參數(shù)、凍結(jié)區(qū)域的有限認(rèn)識,導(dǎo)致在施工過程中為了所謂“絕對安全”而采取過大的凍結(jié)區(qū)和過低的凍結(jié)溫度,從而引起凍結(jié)區(qū)土體過大的凍脹和融沉,最終導(dǎo)致過長的工期和能源浪費。因此,在保障凍結(jié)法施工安全的前提下,合理確定凍結(jié)區(qū)域和凍結(jié)溫度就顯得非常重要[1]。

    目前,人工凍結(jié)溫度場演變規(guī)律研究主要有理論解析、現(xiàn)場實測(凍結(jié)試驗)和數(shù)值模擬3種方法。其中,在理論解析方面,不少學(xué)者[2-6]針對凍結(jié)溫度場演變規(guī)律進(jìn)行了研究。另外,文獻(xiàn)[7-10]基于現(xiàn)場實測資料(凍結(jié)試驗)方法及文獻(xiàn)[11-14]基于數(shù)值模擬計算方法研究凍結(jié)溫度場變化結(jié)果對于凍結(jié)溫度場演變規(guī)律的理解、分析和解釋也較為重要。然而,單管凍結(jié)溫度場演變規(guī)律的準(zhǔn)確解析表達(dá)是研究多凍結(jié)管溫度場的基礎(chǔ)前提。目前,考慮土中水相變(潛熱)熱傳導(dǎo)過程的凍結(jié)溫度場演變規(guī)律研究還不成熟,且已有研究結(jié)果大多未得到有效可靠驗證。

    由此,將單根凍結(jié)管簡化為軸對稱平面熱傳導(dǎo)問題,結(jié)合凍結(jié)鋒面將受影響的土體溫度范圍劃分為凍結(jié)區(qū)和降溫區(qū),采用變量代換法求解相變熱傳導(dǎo)方程,同時基于指數(shù)積分函數(shù)的級數(shù)展開式得到土體凍結(jié)區(qū)及降溫區(qū)的溫度分布表達(dá)式。結(jié)合天津地鐵2號線機(jī)場延長線1號聯(lián)絡(luò)通道凍結(jié)工程現(xiàn)場資料,進(jìn)行了單根凍結(jié)管凍結(jié)槽試驗。最后將單根凍結(jié)管凍結(jié)溫度場近似解析解與凍結(jié)試驗結(jié)果進(jìn)行比對,驗證了理論分析的有效可靠性。

    1軸對稱熱源溫度場方程的解析解

    1.1基本假定

    本文主要對凍結(jié)區(qū)和降溫區(qū)(凍結(jié)管引起的周圍傳熱介質(zhì)凍結(jié)區(qū)和降溫區(qū)如圖1所示)采用圓管穩(wěn)定導(dǎo)熱方程進(jìn)行計算。為便于理論推導(dǎo)作如下假定:① 將研究問題簡化為單管軸對稱平面熱傳導(dǎo)問題;② 凍結(jié)管周圍土體為均質(zhì)連續(xù)體;③ 凍結(jié)管所吸收熱量不隨時間變化;④ 忽略凍結(jié)管徑影響(為冷源點);⑤ 遠(yuǎn)場土體溫度為恒值。

    1.2相變溫度場控制方程

    凍結(jié)區(qū)溫度場如下

    T1t=a1(?2T1r?2+1r?2T1r?2)(ξ≤r<∞)(1)

    降溫區(qū)溫度場如下

    T2t=a2(?2T2r?2+1r?2T2r?2)(r0≤r<ξ)(2)

    邊界條件

    T1(r,0)=T0(3)

    T1(∞,t)=T0(4)

    1.3相變熱傳導(dǎo)方程

    T1(ξ,t)=T2(ξ,t)=Td(5)

    T2(R0,t)=Tc(6)

    λ2T2(ξ,t)r-λ1T1(ξ,t)r=ψt(7)

    式中,T1(r, t),T2(r, t)分別為凍結(jié)區(qū)、降溫區(qū)的溫度,℃;r是以凍結(jié)管中心為原點的圓柱坐標(biāo),m;t為凍結(jié)時間,d;T0為土體的初始溫度,℃;Td,Tc分別為土體的凍結(jié)溫度和凍結(jié)管外壁溫度,℃;λ1,λ2分別為降溫區(qū)、凍結(jié)區(qū)的導(dǎo)熱系數(shù),kJ/(m·d·℃),其中λn=an·cn;an為導(dǎo)溫系數(shù),m?2/d;cn為容積比熱,kJ/(m?2·℃);R0為凍結(jié)管布置圈半徑,m;R為凍結(jié)影響半徑,m;ξN為凍結(jié)鋒面在N區(qū)域內(nèi)的坐標(biāo),m,當(dāng)N=1時,表示0≤ξ1≤R0,當(dāng)N=2時,表示R0≤ξ1≤R;ξ=ξ(t)為凍結(jié)鋒面半徑,m;ψ為單位體積傳熱介質(zhì)的結(jié)冰潛熱,kJ/m?3,ψ=(W-WH)·Ω·ρs,W為土中總含水量,%;WH為凍土中未凍水含量,%;ρs為土體密度,kg/m?3;Ω為水的相變潛熱,通常取336 kJ/kg。

    引入替換變量x1,x2,即

    x1=r?24a1tx2=r?24a2t(8)

    由式(1)~(6),相變熱傳導(dǎo)方程解可表示為

    T1=T0+Td-T0E1(β)E1(r?24a1t)(ξ≤r<∞)(9)

    T2=Tc+Td-Tcg(aβ)g(r?24a2t)(R0≤r<ξ)(10)

    式中,a=a1a2;β為待求參數(shù),E1(x)=∫?∞xe-xxdx,E1(x)為指數(shù)積分函數(shù),用級數(shù)可表示為

    E1(x)=-0.577216-lnx-∞i=1(-1)?ix?ii·i!(11)

    為更容易求得函數(shù)值,E1(x)還用多項式表示為

    E1(x)=-lnx-5i=0bix?i+ε(x)|ε(x)|<2.0×10-70

    式中,b0=-0.577 215 66,b1=0.999 991 93,b2=-0.249 910 55,b3=0.055 199 68,b4=-0.009 760 04,b5=0.001 078 57。

    xe?xE1(x)=x?2+d1x+d2x?2+d3x+d4+ε(x)|ε(x)|<5.0×10-51≤x≤∞(13)

    式中,d1=2.334 733,d2=0.250 621,d3=3.330 657,d4=1.681 534。

    β=ξ?24a1t(14)

    將式(9)、(10)代入式(7),得

    a1ψβ-λ1(Td-T0)E1(β)e-β-λ2(Td-Tc)e-aβg(aβ)=0(15)

    進(jìn)而有

    ξ=2a1βt(16)

    凍結(jié)能力Q為

    Q=2πR0λ2T2(R0,t)r(17)

    將式(10)代入式(17),則有

    Q=4πλ2Td-Tc(R0,t)g(aβ)e-R?20a2(18)

    由式(18)可知,對Q影響最大的是λ2和(Td-Tc),其呈現(xiàn)出正比例關(guān)系?;谟嬎阏`差不超過

    1%,有

    Q=4πλ2Td-Tcg(aβ)(t>25R?20a2)(19)

    2解適用條件及無量綱表達(dá)

    作變量替換如下

    Tn=Td+(T0-Td)θnn=1,2t=R?20a1τR=R0rξ=R0bq=Qλ2(T0-Td)(20)

    式中,θ1,θ2分別為未凍土和凍結(jié)土無量綱溫度;τ為無量綱時間;r,b分別為無量綱徑向坐標(biāo)和凍結(jié)鋒面半徑;q為凍結(jié)管壁無量綱熱流密度。

    將式(20)分別代入式(9)和式(10)得

    θ(r,τ)=1-f(r?24τ)f(β)(b≤r<∞)(21)

    θ2(r,τ)=θc[1-g(ar?2/4τ)g(aβ)](0

    式中,θc=(Tc-Td)/(T0-Td)。

    參數(shù)β為

    β=b?24τ(23)

    此時,確定參數(shù)β的超越方程式(15)為

    Kβ+e-βf(β)-λθce-aβg(aβ)=0(24)

    式中,K=ψ/[c1(T0-Td)];λ=λ1/λ2。

    由超越方程(24),得參數(shù)β,將其再代入式(23)可得凍結(jié)鋒面的無量綱半徑b為

    b=2βτ(25)

    由此,無量綱熱流密度為

    q=-4πθcg(aβ)e-a4τ(26)

    q=-4πθcg(aβ)(τ>25a)(27)

    3溫度場凍結(jié)槽試驗

    為驗證理論推導(dǎo)的正確性,在凍結(jié)模型槽內(nèi)進(jìn)行穩(wěn)定冷源作用下單根凍結(jié)管溫度場演變試驗。凍結(jié)試驗土樣取自天津地鐵2號線機(jī)場延長線工程項目。聯(lián)絡(luò)通道底板以上位于含黏性土粉砂層中,通過對不同地層土樣進(jìn)行熱物理力學(xué)參數(shù)試驗[15],所得凍土熱物理參數(shù)如表1所列。

    單根凍結(jié)管溫度場試驗槽如圖2(a)所示,槽內(nèi)加墊石棉纖維毯(厚度25 mm),如圖2(b)所示。試驗槽內(nèi)部尺寸為1.1 m×0.25 m×0.8 m。在試驗槽中央插入Φ42 mm×0.8 m的凍結(jié)管,沿試驗槽長度方向0.4 m土層深度處每隔0.1m布置一個測溫點,由凍結(jié)管徑向邊界測點編號分別為C1~C5。凍結(jié)試驗中,室溫為10℃,空氣濕度為40%。使用鹽水作為低溫冷媒劑。凍結(jié)孔單孔流量為5 m?3/h,7d的積極凍結(jié)使鹽水溫度降至-20℃以下,去、回路的鹽水溫差控制在2℃以內(nèi),保證提供穩(wěn)定冷源溫度為-20℃,積極凍結(jié)時間為60 d。

    4結(jié)果分析與討論

    據(jù)凍結(jié)資料及凍土和未凍土的熱物理性質(zhì),計算參數(shù)取值如下:R0=0.08 m,T0=10℃,Td=0℃,Tc=-20℃,λ1=100 kJ/(m·d·℃),λ2=135 kJ/(m·d·℃),c1=2 900 kJ/(m?3·℃),c2=2 100kJ/(m?3·℃),ψ=2.0×10?5 kJ/m?3,ρs=1 600 kg/m?3,W=32%,WH=5.0%。由式(15)或(24)迭代計算求得β=0.2。由式(16)得ξ=0.058t。當(dāng)t>3 d后,凍結(jié)管壁熱流密度Q即處于穩(wěn)定狀態(tài),此時Q=1.12×10?4 kJ/(m·d)。

    圖3為與凍結(jié)管間距為0.3 m的測溫孔C3凍結(jié)溫度隨凍結(jié)時間變化曲線。由圖3可知,測溫孔凍結(jié)溫度與凍結(jié)時間關(guān)系為瞬態(tài)溫度演變,降溫區(qū)和凍結(jié)區(qū)溫度演變均符合對數(shù)變化規(guī)律,且在凍結(jié)初期,理論解析封閉解大于凍結(jié)試驗實測值。

    圖4為凍結(jié)過程中分別為20,40 d和60 d時刻的凍土圓柱面內(nèi)各測溫點凍結(jié)溫度分布情況。凍結(jié)瞬時凍土柱面內(nèi)各測溫點凍結(jié)溫度為穩(wěn)態(tài)的,凍結(jié)初期理論計算值下降較快,凍結(jié)后期其與凍結(jié)試驗實測溫度值接近。

    圖5位凍結(jié)鋒面半徑隨凍結(jié)時間由凍結(jié)管徑向向外邊界擴(kuò)展曲線。凍結(jié)鋒面半徑隨凍結(jié)時間增長呈現(xiàn)出冪指數(shù)關(guān)系,凍結(jié)初期理論計算值大于試驗結(jié)果,而凍結(jié)后期,凍結(jié)試驗得出的凍結(jié)鋒面半徑略大于理論解析計算結(jié)果。

    試驗過程中土樣的初始溫度、凍結(jié)溫度、土水耦合傳熱介質(zhì)的相變潛熱及鹽水溫度等因素的變化對凍結(jié)鋒面均有顯著影響,但理論解析解與試驗結(jié)果總體變化規(guī)律和趨勢基本一致。理論計算與凍結(jié)試驗所得凍結(jié)鋒面半徑結(jié)果誤差為10%以內(nèi),為可接受范圍。由此可知,所推導(dǎo)理論計算結(jié)果是有效可靠的。由于凍結(jié)試驗過程耗時較長,時間成本和經(jīng)濟(jì)成本均較高,可采用經(jīng)由試驗有效可靠性檢驗的理論解析解分析更多凍結(jié)溫度場演變過程。

    5結(jié) 論

    (1) 針對穩(wěn)定線性冷源作用下土體凍結(jié)過程中的瞬態(tài)溫度場進(jìn)行了解析表達(dá),其中,單管凍結(jié)鋒面半徑和凍結(jié)時間有著平方根關(guān)系,并由凍結(jié)槽試驗驗證了理論解析解的有效可靠性。

    (2) 由指數(shù)積分函數(shù)近似級數(shù)形式求得了距離凍結(jié)管不同測點的溫度。凍結(jié)溫度場在時間和空間上呈對數(shù)型函數(shù)分布,且曲率隨凍結(jié)時間的增加而逐漸減小。

    (3) 由凍結(jié)相變過程分析得出了凍結(jié)鋒面半徑與時間之間的理論解析表達(dá)式。基于凍結(jié)鋒面半徑隨凍結(jié)時間變化規(guī)律參數(shù)化的敏感性分析,優(yōu)化凍結(jié)參數(shù),有效地評估凍結(jié)法對周圍環(huán)境的影響。在保證滿足凍結(jié)法施工要求的前提下,施工可做到節(jié)約能源、縮短工期。

    今后研究應(yīng)考慮多根凍結(jié)管不同布置方式(并行、交叉等),研究凍結(jié)溫度場交圈、凍結(jié)管間距、凍結(jié)管排數(shù)等更為復(fù)雜因素引起的傳熱介質(zhì)凍結(jié)溫度場隨凍結(jié)時間的變化規(guī)律。研究成果對認(rèn)識不良地基土體凍結(jié)溫度場分布瞬態(tài)演變規(guī)律和進(jìn)一步指導(dǎo)凍結(jié)法施工設(shè)計等具有重要的工程意義。

    參考文獻(xiàn):

    [1]馬芹永.人工凍結(jié)法的理論與施工技術(shù)[M].北京:人民交通出版社,2007.

    [2]周揚,周國慶.考慮未凍水單管凍結(jié)溫度場解析解[J].煤炭學(xué)報,2012,37(10):1649-1653.

    [3]李方政,夏明萍.基于指數(shù)積分函數(shù)的人工凍土溫度場解析研究[J].東南大學(xué)學(xué)報:自然科學(xué)版,2004,34(4):469-473.

    [4]胡向東,黃峰,白楠.考慮上層凍結(jié)溫度時人工凍結(jié)溫度場模型[J].中國礦業(yè)大學(xué)學(xué)報,2008,37(4):550-555.

    [5]張弛,張濤,韓濤,等.管壁溫度非恒定條件下單管凍結(jié)溫度場解析計算[J].煤炭科學(xué)技術(shù),2012,40(3):20-23.

    [6]蔣斌松,王金鴿,周國慶.單管凍結(jié)溫度場解析計算[J].中國礦業(yè)大學(xué)學(xué)報,2009,38(4):463-466.

    [7]鮑永亮,鄭七振,唐建忠.地鐵隧道旁通道凍結(jié)法施工監(jiān)測分析[J].鐵道工程學(xué)報,2009,(3):93-95.

    [8]王志良,申林方,李明宇.凍結(jié)法施工的地鐵旁通道實測數(shù)據(jù)分析[J].地下空間與工程學(xué)報,2010,6(1):138-143.

    [9]李大勇,呂愛鐘,張慶賀,等.南京地鐵旁通道凍結(jié)實測分析研究[J].巖石力學(xué)與工程學(xué)報,2004,23(2):334-338.

    [10]程學(xué)磊,李順群,孫世娟,等.圍壓和負(fù)溫對凍土強(qiáng)度和微結(jié)構(gòu)的影響研究[J].廣西大學(xué)學(xué)報:自然科學(xué)版,2014,39(1):95-104.

    [11]畢貴權(quán),程形燕,石磊,等.地鐵隧道水平凍結(jié)法施工凍結(jié)壁溫度場影響參數(shù)分析[J].蘭州理工大學(xué)學(xué)報,2009,35(3):121-126.

    [12]尚軍,程學(xué)磊,張來棟,等.引水溝渠凍結(jié)法施工中的溫度場、凍結(jié)壁演變和開挖穩(wěn)定性分析[J].水電能源科學(xué),2015(7):139-143.

    [13]王效賓,楊平,張婷,等.盾構(gòu)出洞水平凍結(jié)解凍溫度場三維有限元分析[J].解放軍理工大學(xué)學(xué)報:自然科學(xué)版,2009,10(6):586-590.

    [14]靳巍巍,陳有亮,李磊,等.隧道聯(lián)絡(luò)通道凍結(jié)法施工三維有限元溫度場分析[J].上海大學(xué)學(xué)報:自然科學(xué)版,2008,14(1):85-90.

    [15]陳之祥,李順群,夏錦紅,等.基于緊密排列土柱模型的凍土熱參數(shù)計算[J].深圳大學(xué)學(xué)報:理工版,2017,34(4):393-399.

    (編輯:鄭 毅)

    引用本文:馬蕓,夏錦紅,桂超.單管凍結(jié)溫度場演變規(guī)律分析及試驗驗證[J].人民長江,2019,50(2):169-173.

    Development law of soil temperature field caused by freezing pipe and its freezing test verification

    MA Yun?, XIA Jinhong1,2, GUI Chao1,3

    (1.School of Civil Engineering and Architecture, Xinxiang University, Xinxiang 453003, China;2.School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China;3.School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

    Abstract: A phase-transition temperature field model for artificial single pipe freezing with a stable linear cold source was established and approximate analytical solution for freezing temperature development with time was obtained by variable substitution method. The temperature distributions in the frozen and non-frozen fields were expressed by an exponential integral function, and the theoretical results were compared with the freezing model test. The research results show that the error between theoretical calculation and test result is within a reasonable range, which verified the effectiveness and reliability of the approximate analytical solution for transient freezing temperature field distribution; the relation between the freezing-front radius and freezing time follows square root equation; the temperature field for a single freezing pipe obeys a logarithmic function in time and space and its curvature decreases with time. The research results have an important engineering significance to understanding the temperature field evolution law of poor foundation soil surrounding a freezing pipe under a stable linear cold source and guiding the freezing construction technology.

    Key words:freezing temperature field; phase change; frozen front radius; theoretical analysis; freezing tank test

    亚洲天堂av无毛| 亚洲欧美清纯卡通| 在现免费观看毛片| 亚洲中文字幕日韩| 又紧又爽又黄一区二区| 午夜福利影视在线免费观看| 日韩中文字幕视频在线看片| 波多野结衣一区麻豆| 久久av网站| 欧美成人午夜精品| 制服诱惑二区| 国产亚洲欧美在线一区二区| 伊人久久大香线蕉亚洲五| 青青草视频在线视频观看| 日韩中文字幕视频在线看片| 别揉我奶头~嗯~啊~动态视频 | 免费日韩欧美在线观看| 大片电影免费在线观看免费| 国产国语露脸激情在线看| 又大又黄又爽视频免费| 五月开心婷婷网| 久久久久精品人妻al黑| 美女脱内裤让男人舔精品视频| 人成视频在线观看免费观看| 欧美xxⅹ黑人| 在线精品无人区一区二区三| 久久亚洲国产成人精品v| 丝袜喷水一区| 国语对白做爰xxxⅹ性视频网站| 日本猛色少妇xxxxx猛交久久| 后天国语完整版免费观看| 国产成人影院久久av| 高清黄色对白视频在线免费看| 欧美在线一区亚洲| 日本a在线网址| h视频一区二区三区| 免费看不卡的av| 成人免费观看视频高清| 欧美在线一区亚洲| 狠狠婷婷综合久久久久久88av| 午夜免费男女啪啪视频观看| 国产免费又黄又爽又色| 一级黄色大片毛片| 免费一级毛片在线播放高清视频 | 永久免费av网站大全| 亚洲精品美女久久久久99蜜臀 | 久久精品亚洲av国产电影网| 免费高清在线观看日韩| 蜜桃在线观看..| 国产日韩欧美亚洲二区| 国产日韩欧美视频二区| 丰满少妇做爰视频| 免费一级毛片在线播放高清视频 | 大片免费播放器 马上看| 多毛熟女@视频| 日本午夜av视频| 天天操日日干夜夜撸| 国产亚洲欧美在线一区二区| 在线观看免费视频网站a站| 自线自在国产av| 午夜两性在线视频| 夜夜骑夜夜射夜夜干| 久久久精品94久久精品| 国产精品成人在线| 看十八女毛片水多多多| 又黄又粗又硬又大视频| av国产精品久久久久影院| 国产女主播在线喷水免费视频网站| 母亲3免费完整高清在线观看| 国产成人精品在线电影| 波多野结衣av一区二区av| 国产精品久久久久久精品古装| 亚洲一区中文字幕在线| 国产精品99久久99久久久不卡| 久久 成人 亚洲| av天堂久久9| 午夜激情久久久久久久| 亚洲国产av影院在线观看| 亚洲精品日韩在线中文字幕| 永久免费av网站大全| 亚洲九九香蕉| 国产真人三级小视频在线观看| 国产又色又爽无遮挡免| 免费女性裸体啪啪无遮挡网站| 久9热在线精品视频| 亚洲国产中文字幕在线视频| 丁香六月欧美| 各种免费的搞黄视频| 天天添夜夜摸| 欧美激情极品国产一区二区三区| 久久久精品免费免费高清| 尾随美女入室| 自线自在国产av| 久久久久久免费高清国产稀缺| 在线 av 中文字幕| 一边摸一边做爽爽视频免费| 国产成人免费无遮挡视频| 亚洲国产精品一区三区| 国产一区亚洲一区在线观看| 亚洲av男天堂| 日本av免费视频播放| 亚洲欧美一区二区三区黑人| 国产精品香港三级国产av潘金莲 | 夫妻性生交免费视频一级片| 少妇被粗大的猛进出69影院| 亚洲av在线观看美女高潮| 日韩 亚洲 欧美在线| 日本欧美国产在线视频| 久久久久久久国产电影| 亚洲 国产 在线| 一级黄片播放器| 秋霞在线观看毛片| 国产视频一区二区在线看| 婷婷色av中文字幕| 五月天丁香电影| 一区二区三区四区激情视频| 看免费av毛片| 免费久久久久久久精品成人欧美视频| 一边亲一边摸免费视频| 人人澡人人妻人| 亚洲欧美中文字幕日韩二区| 午夜视频精品福利| 午夜福利乱码中文字幕| 亚洲精品美女久久久久99蜜臀 | 精品一区二区三区av网在线观看 | 国产精品成人在线| 国产人伦9x9x在线观看| 精品少妇一区二区三区视频日本电影| 久久亚洲精品不卡| 亚洲av日韩在线播放| 精品久久蜜臀av无| 不卡av一区二区三区| 美女午夜性视频免费| 欧美激情高清一区二区三区| 日韩av不卡免费在线播放| 亚洲激情五月婷婷啪啪| 一级片免费观看大全| 波多野结衣一区麻豆| 亚洲欧洲国产日韩| 午夜激情久久久久久久| 亚洲欧美一区二区三区久久| 99久久人妻综合| 老司机影院毛片| 欧美黄色淫秽网站| 女性被躁到高潮视频| 日韩 欧美 亚洲 中文字幕| 亚洲av日韩在线播放| 国产激情久久老熟女| 一级黄色大片毛片| 亚洲男人天堂网一区| 国产一区二区三区av在线| 黄色怎么调成土黄色| 亚洲精品国产av成人精品| 国产精品一二三区在线看| 婷婷色av中文字幕| 波多野结衣av一区二区av| 久久精品aⅴ一区二区三区四区| 99re6热这里在线精品视频| 老鸭窝网址在线观看| 黑人猛操日本美女一级片| 欧美成人午夜精品| 在线精品无人区一区二区三| 婷婷色综合大香蕉| 91精品国产国语对白视频| www.自偷自拍.com| 免费高清在线观看日韩| 久久久久视频综合| 成年人午夜在线观看视频| 欧美黑人精品巨大| 久久精品成人免费网站| 黄色a级毛片大全视频| 免费在线观看影片大全网站 | 久久亚洲国产成人精品v| 老司机影院毛片| 青草久久国产| 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产色婷婷电影| 丰满少妇做爰视频| 久久国产精品男人的天堂亚洲| 精品国产乱码久久久久久男人| 大片电影免费在线观看免费| 99热全是精品| 精品人妻熟女毛片av久久网站| 在线观看免费午夜福利视频| 精品少妇内射三级| 精品人妻一区二区三区麻豆| 性色av乱码一区二区三区2| 日韩av免费高清视频| 亚洲成人免费电影在线观看 | 精品少妇久久久久久888优播| 午夜两性在线视频| 涩涩av久久男人的天堂| 少妇被粗大的猛进出69影院| 国产精品一区二区在线观看99| av网站免费在线观看视频| 视频区图区小说| √禁漫天堂资源中文www| 欧美日韩av久久| kizo精华| 久9热在线精品视频| 91九色精品人成在线观看| 日本欧美视频一区| 麻豆av在线久日| 丁香六月欧美| 在线 av 中文字幕| 欧美激情 高清一区二区三区| 亚洲人成电影免费在线| 人妻人人澡人人爽人人| av在线app专区| 国产欧美亚洲国产| 欧美亚洲 丝袜 人妻 在线| 人人妻人人澡人人看| 免费av中文字幕在线| 一本久久精品| 久久人人97超碰香蕉20202| 亚洲国产精品999| 久久精品国产亚洲av高清一级| 亚洲免费av在线视频| 免费在线观看日本一区| 亚洲国产最新在线播放| 国产日韩欧美视频二区| 亚洲国产精品国产精品| 欧美激情 高清一区二区三区| 午夜福利免费观看在线| 精品久久久久久久毛片微露脸 | 午夜福利视频精品| 欧美性长视频在线观看| 黄色视频在线播放观看不卡| 免费日韩欧美在线观看| 91老司机精品| 亚洲av片天天在线观看| 777久久人妻少妇嫩草av网站| 99热全是精品| 欧美+亚洲+日韩+国产| 天天影视国产精品| 精品人妻在线不人妻| 久久天堂一区二区三区四区| a级片在线免费高清观看视频| 欧美人与善性xxx| 精品久久久久久久毛片微露脸 | 九色亚洲精品在线播放| 国产精品国产三级国产专区5o| 美女视频免费永久观看网站| 天天躁夜夜躁狠狠躁躁| 女人高潮潮喷娇喘18禁视频| 午夜福利乱码中文字幕| 夫妻性生交免费视频一级片| 免费高清在线观看日韩| 最近中文字幕2019免费版| 国产av精品麻豆| 天天躁夜夜躁狠狠久久av| 别揉我奶头~嗯~啊~动态视频 | 美女脱内裤让男人舔精品视频| 欧美激情高清一区二区三区| 波野结衣二区三区在线| 久久国产精品大桥未久av| 精品国产乱码久久久久久男人| 老鸭窝网址在线观看| 国产熟女午夜一区二区三区| 成在线人永久免费视频| 国产免费福利视频在线观看| 日本色播在线视频| 国产熟女午夜一区二区三区| 可以免费在线观看a视频的电影网站| 国产男女内射视频| 9色porny在线观看| 亚洲欧美中文字幕日韩二区| 蜜桃国产av成人99| 久久久久久久国产电影| 欧美国产精品一级二级三级| 七月丁香在线播放| 久久热在线av| xxxhd国产人妻xxx| 国产精品久久久久成人av| 一区二区三区激情视频| 亚洲熟女毛片儿| 国产日韩欧美视频二区| 一二三四在线观看免费中文在| 日韩一区二区三区影片| 成年av动漫网址| 亚洲国产日韩一区二区| 如日韩欧美国产精品一区二区三区| 国产亚洲午夜精品一区二区久久| 精品国产国语对白av| 亚洲精品久久午夜乱码| 这个男人来自地球电影免费观看| 久久热在线av| 韩国精品一区二区三区| 精品免费久久久久久久清纯 | 亚洲欧美色中文字幕在线| 国产一卡二卡三卡精品| 最近最新中文字幕大全免费视频 | 乱人伦中国视频| 国产主播在线观看一区二区 | 国产精品欧美亚洲77777| 99久久99久久久精品蜜桃| 欧美 日韩 精品 国产| 午夜免费观看性视频| 国产有黄有色有爽视频| 中文字幕人妻熟女乱码| svipshipincom国产片| 婷婷色综合www| 成年av动漫网址| 日本黄色日本黄色录像| 超色免费av| 精品一区在线观看国产| 欧美日韩av久久| 国产精品香港三级国产av潘金莲 | 啦啦啦视频在线资源免费观看| 亚洲三区欧美一区| 国产一区二区三区综合在线观看| 日韩人妻精品一区2区三区| 久久久久久亚洲精品国产蜜桃av| 一区二区三区四区激情视频| 天堂俺去俺来也www色官网| 欧美变态另类bdsm刘玥| 亚洲av日韩在线播放| 女人精品久久久久毛片| 中文字幕高清在线视频| 午夜视频精品福利| 91精品三级在线观看| 国产一区二区三区av在线| 午夜福利视频精品| 国产福利在线免费观看视频| 大香蕉久久网| 国产成人一区二区三区免费视频网站 | 国产一卡二卡三卡精品| 这个男人来自地球电影免费观看| 又粗又硬又长又爽又黄的视频| 人妻人人澡人人爽人人| 国产野战对白在线观看| 亚洲人成77777在线视频| 后天国语完整版免费观看| 夫妻午夜视频| 最近手机中文字幕大全| 另类亚洲欧美激情| 一本一本久久a久久精品综合妖精| 国产色视频综合| 成人18禁高潮啪啪吃奶动态图| 少妇 在线观看| 国产午夜精品一二区理论片| 母亲3免费完整高清在线观看| 精品少妇久久久久久888优播| 一个人免费看片子| 精品国产乱码久久久久久男人| 69精品国产乱码久久久| 多毛熟女@视频| 国产亚洲一区二区精品| 欧美日本中文国产一区发布| 啦啦啦在线观看免费高清www| 秋霞在线观看毛片| 久久av网站| 纵有疾风起免费观看全集完整版| 青春草亚洲视频在线观看| 欧美精品亚洲一区二区| 精品视频人人做人人爽| 亚洲熟女精品中文字幕| 久久精品久久精品一区二区三区| 飞空精品影院首页| 看免费成人av毛片| 久久精品久久久久久噜噜老黄| 久久国产亚洲av麻豆专区| 自线自在国产av| 亚洲九九香蕉| 男女高潮啪啪啪动态图| 热99国产精品久久久久久7| 1024香蕉在线观看| 国产成人91sexporn| 国产精品久久久久成人av| 在线观看免费视频网站a站| 亚洲美女黄色视频免费看| 久久国产亚洲av麻豆专区| 欧美97在线视频| 婷婷色麻豆天堂久久| 久久午夜综合久久蜜桃| 久久国产精品影院| 男女床上黄色一级片免费看| 免费在线观看黄色视频的| av不卡在线播放| 亚洲国产看品久久| 国产女主播在线喷水免费视频网站| 国产精品三级大全| 啦啦啦啦在线视频资源| 十八禁网站网址无遮挡| 黄色片一级片一级黄色片| 亚洲欧美中文字幕日韩二区| 亚洲成人国产一区在线观看 | 中文字幕av电影在线播放| 亚洲欧美中文字幕日韩二区| 国产在线免费精品| 成年女人毛片免费观看观看9 | 亚洲欧洲日产国产| 老熟女久久久| 伊人亚洲综合成人网| 亚洲国产中文字幕在线视频| 亚洲国产精品成人久久小说| 91精品三级在线观看| 日韩一卡2卡3卡4卡2021年| av不卡在线播放| 一区二区日韩欧美中文字幕| 日韩大码丰满熟妇| 久久午夜综合久久蜜桃| 老司机靠b影院| 免费观看av网站的网址| 大码成人一级视频| 男人爽女人下面视频在线观看| 一区二区日韩欧美中文字幕| 一区二区三区精品91| 在线观看免费高清a一片| 老熟女久久久| 亚洲国产精品一区二区三区在线| 久久久精品94久久精品| 亚洲 国产 在线| 免费观看人在逋| 午夜久久久在线观看| 另类精品久久| 欧美少妇被猛烈插入视频| 久久精品亚洲熟妇少妇任你| 国产又色又爽无遮挡免| 免费在线观看黄色视频的| 一本大道久久a久久精品| 两个人免费观看高清视频| 欧美日韩视频精品一区| 久久精品国产综合久久久| 久久精品熟女亚洲av麻豆精品| 丝袜人妻中文字幕| 亚洲色图综合在线观看| 亚洲少妇的诱惑av| 交换朋友夫妻互换小说| 看免费av毛片| 亚洲欧美色中文字幕在线| 欧美亚洲 丝袜 人妻 在线| 中文字幕av电影在线播放| 国产精品麻豆人妻色哟哟久久| 中国美女看黄片| 久久国产亚洲av麻豆专区| 精品一区二区三卡| 国产成人影院久久av| 精品久久久精品久久久| 国产又色又爽无遮挡免| xxx大片免费视频| 啦啦啦啦在线视频资源| 国产欧美亚洲国产| 狂野欧美激情性xxxx| 丝袜在线中文字幕| 亚洲黑人精品在线| 亚洲精品久久午夜乱码| 无遮挡黄片免费观看| 两性夫妻黄色片| 午夜激情av网站| 国产在线免费精品| 久久久久精品国产欧美久久久 | 亚洲午夜精品一区,二区,三区| 国产又爽黄色视频| 18禁黄网站禁片午夜丰满| 久久久久国产精品人妻一区二区| 女性生殖器流出的白浆| 久久综合国产亚洲精品| 在线观看免费视频网站a站| 18禁黄网站禁片午夜丰满| 亚洲熟女毛片儿| 交换朋友夫妻互换小说| 99久久99久久久精品蜜桃| 久久精品国产亚洲av涩爱| 成人三级做爰电影| 精品人妻一区二区三区麻豆| 日本色播在线视频| 亚洲精品第二区| 欧美精品一区二区大全| 91精品国产国语对白视频| 久久天堂一区二区三区四区| 色综合欧美亚洲国产小说| 精品人妻一区二区三区麻豆| 欧美日韩福利视频一区二区| 日韩一区二区三区影片| 国产成人精品无人区| 国产成人精品久久二区二区免费| 国产91精品成人一区二区三区 | 国产成人系列免费观看| 日韩伦理黄色片| 妹子高潮喷水视频| 嫁个100分男人电影在线观看 | 亚洲九九香蕉| 国产99久久九九免费精品| 尾随美女入室| 欧美国产精品一级二级三级| 丁香六月天网| 一区在线观看完整版| 亚洲精品一卡2卡三卡4卡5卡 | 欧美成人精品欧美一级黄| 99国产精品一区二区蜜桃av | 无限看片的www在线观看| 人人妻人人爽人人添夜夜欢视频| 99久久精品国产亚洲精品| 国产亚洲一区二区精品| 欧美xxⅹ黑人| 午夜福利视频精品| 欧美激情 高清一区二区三区| 国产成人欧美在线观看 | 女人高潮潮喷娇喘18禁视频| 久久精品亚洲av国产电影网| 在现免费观看毛片| 成人亚洲精品一区在线观看| 国产片内射在线| 国产亚洲av片在线观看秒播厂| 一区在线观看完整版| 国产精品一区二区精品视频观看| 亚洲九九香蕉| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三 | 天天操日日干夜夜撸| 午夜91福利影院| 国产精品二区激情视频| 国产精品三级大全| tube8黄色片| 少妇人妻 视频| 老汉色∧v一级毛片| 亚洲精品自拍成人| 我要看黄色一级片免费的| 成人免费观看视频高清| 久久人人爽av亚洲精品天堂| 一区二区日韩欧美中文字幕| 一本—道久久a久久精品蜜桃钙片| 免费在线观看影片大全网站 | 亚洲精品国产一区二区精华液| 日韩制服丝袜自拍偷拍| 免费人妻精品一区二区三区视频| a 毛片基地| 亚洲七黄色美女视频| 国产精品熟女久久久久浪| 一级毛片黄色毛片免费观看视频| 午夜福利视频精品| 老司机影院毛片| 9191精品国产免费久久| 色94色欧美一区二区| 国产精品一区二区免费欧美 | 精品视频人人做人人爽| 精品久久蜜臀av无| 丁香六月天网| e午夜精品久久久久久久| 新久久久久国产一级毛片| 十分钟在线观看高清视频www| 狠狠精品人妻久久久久久综合| 别揉我奶头~嗯~啊~动态视频 | 成人18禁高潮啪啪吃奶动态图| 精品国产一区二区三区四区第35| 十分钟在线观看高清视频www| 一个人免费看片子| 狂野欧美激情性xxxx| 老熟女久久久| av片东京热男人的天堂| 久久久久久久国产电影| 大香蕉久久成人网| 蜜桃国产av成人99| 国产极品粉嫩免费观看在线| 成人国产一区最新在线观看 | 精品少妇一区二区三区视频日本电影| 久久久久久久久久久久大奶| 欧美国产精品va在线观看不卡| 中文字幕人妻丝袜制服| 国产精品熟女久久久久浪| 大话2 男鬼变身卡| 高清欧美精品videossex| 午夜福利在线免费观看网站| 免费在线观看视频国产中文字幕亚洲 | 久久狼人影院| 国产高清videossex| 国产淫语在线视频| www.精华液| 男女无遮挡免费网站观看| 久久久精品94久久精品| 18禁国产床啪视频网站| 黄片播放在线免费| 日韩,欧美,国产一区二区三区| 国产av国产精品国产| 国产免费现黄频在线看| 亚洲精品久久午夜乱码| 高清视频免费观看一区二区| 欧美另类一区| 亚洲精品久久成人aⅴ小说| 两个人看的免费小视频| 欧美老熟妇乱子伦牲交| 亚洲国产精品一区三区| 91麻豆精品激情在线观看国产 | 欧美日韩综合久久久久久| 99久久综合免费| 99精国产麻豆久久婷婷| 精品国产国语对白av| 国产男女超爽视频在线观看| 一本久久精品| 成年美女黄网站色视频大全免费| 欧美性长视频在线观看| 欧美精品一区二区大全| 亚洲精品第二区| 成年人黄色毛片网站| 午夜精品国产一区二区电影| 国产97色在线日韩免费| 美女午夜性视频免费| 亚洲欧洲日产国产| 成人午夜精彩视频在线观看| 黑人欧美特级aaaaaa片| 十八禁网站网址无遮挡| 亚洲成人手机| 成年美女黄网站色视频大全免费| 国产精品一区二区精品视频观看| 国产黄色免费在线视频| 蜜桃国产av成人99| 精品国产一区二区久久| 亚洲av在线观看美女高潮| 777米奇影视久久| 80岁老熟妇乱子伦牲交| 欧美久久黑人一区二区| 亚洲精品美女久久久久99蜜臀 | 国产91精品成人一区二区三区 | 国产精品一区二区在线观看99|