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

    雙圓盾構(gòu)隧道杯型凍結(jié)壁溫度場數(shù)值優(yōu)化分析

    2022-06-15 04:20:18潘名芳董勤喜胡俊曾暉
    森林工程 2022年3期
    關(guān)鍵詞:數(shù)值模擬隧道

    潘名芳 董勤喜 胡俊 曾暉

    摘 要:為得到雙圓盾構(gòu)隧道始發(fā)端頭應(yīng)用水平杯型凍結(jié)法進行加固時的溫度場變化規(guī)律,該文結(jié)合上海市軌道交通6號線某雙圓盾構(gòu)隧道區(qū)間工程,通過使用有限元數(shù)值模擬軟件建立數(shù)值模型,分析初始模型設(shè)計102根凍結(jié)管時雙圓盾構(gòu)隧道端頭杯型凍結(jié)壁溫度場的發(fā)展規(guī)律,并提出2種優(yōu)化方案:方案1為內(nèi)圈、中圈、外圈凍結(jié)管根數(shù)分別減少2、4、6根;方案2為內(nèi)圈、中圈、外圈凍結(jié)管根數(shù)分別減少4、6、10根,分析對比初始模型以及2個優(yōu)化方案的凍結(jié)效果。結(jié)果表明,在凍結(jié)開始后,凍土帷幕圍繞各凍結(jié)管呈圓柱形逐漸向外延伸,距離凍結(jié)管越近,凍結(jié)效果越好;凍土帷幕依次按外圈管到中圈管,再到內(nèi)圈管的次序開始交圈,凍結(jié)第18天形成一個形狀規(guī)則、強度分布均勻的近似雙圓杯型凍結(jié)壁;雙圓杯型凍結(jié)壁杯底長度10.2 m,杯身高度4.4 m,杯壁厚度2 m,杯底厚度2 m;優(yōu)化前初始凍結(jié)方案使用102根凍結(jié)管符合施工要求,凍結(jié)效果良好,在實際施工中具有可行性,但偏于保守,2種優(yōu)化方案在減少凍結(jié)管數(shù)量后依然滿足盾構(gòu)始發(fā)杯型凍結(jié)壁杯體要求。

    關(guān)鍵詞:隧道;雙圓盾構(gòu);端頭加固;數(shù)值模擬;杯型凍結(jié)壁;凍結(jié)法

    中圖分類號:S 773;U 455.43??? 文獻標識碼:A?? 文章編號:1006-8023(2022)03-0122-012

    Numerical Optimization Analysis of Temperature Field of Cup-shaped

    Frozen Wall of Double-Circle Shield Tunnel

    PAN Mingfang1, DONG Qinxi1, HU Jun1*, ZENG Hui2

    (1.College of Civil Engineering and Architecture, Hainan University, Haikou 570228, China;

    2.College of Civil Engineering and Architecture, Wuyi University, Jiangmen 529020, China)

    Abstract:In order to obtain the variation law of temperature field at the starting end of double-circle shield tunnel reinforced by horizontal cup freezing method, combined with a double-circle shield tunnel section project of Shanghai Rail Transit Line 6, this paper established a numerical model by using finite element numerical simulation software, analyzed the development law of temperature field of cup-shaped frozen wall at the end of double-circle shield tunnel when the initial model designed with 102 frozen pipes. Two optimization schemes were put forward: scheme 1 was reducing the number of freezing pipes in inside, center and outside circles by 2, 4 and 6 respectively; scheme 2 was reducing the number of freezing pipes in inside, center and outside circles by 4, 6 and 10 respectively. The freezing effects of the initial model and two optimization schemes were analyzed and compared. The results showed that: at the beginning of freezing, the frozen soil curtain gradually extended outward around each freezing pipe in a cylindrical shape, and the closer it was to the freezing pipe, the better the freezing effect. The frozen soil curtain began to cross in the order of outer circle, center circle and inside circle, and formed an approximate double-round cup-shaped frozen wall with regular shape and uniform strength distribution on the 18th day of freezing. The length of double round cup-shaped frozen wall cup bottom was 10.2 m, the height of cup body was 4.4 m, the thickness of cup wall was 2m, and the thickness of cup bottom was 2 m.The initial freezing scheme which using 102 freezing pipes met the construction requirements, the freezing effect was good, and it was feasible in actual construction, but it was conservative. The two optimization schemes still met the starting cup-shaped freezing wall cup of shield machine requirements after reducing the number of freezing pipes.

    Keywords:Tunnel; double circle shield; end reinforcement; numerical simulation; cup type freezing wall; freezing method

    0 引言

    在地鐵建設(shè)中,盾構(gòu)掘進施工具有極高程度實現(xiàn)自動化、減少人力和快速施工等優(yōu)點,且盾構(gòu)施工受溫度氣候條件影響小,施工過程中可減小對鄰近建筑物影響,易于控制土體沉降。但盾構(gòu)掘進法也有許多施工難點,如始發(fā)時對隧道端頭土體的加固問題 [1-3]。雙圓盾構(gòu)相比于傳統(tǒng)的單圓盾構(gòu),作為一種典型的異形盾構(gòu)(雙圓、多圓盾構(gòu)),具有更高的施工效率,不需要掘削太多的土體,減少了人力物力的投入,更重要的是在掘進完成后,只需要打通海鷗型管片下方的立柱支撐,就可形成聯(lián)絡(luò)通道,省去了很大一部分工程量。由于雙圓盾構(gòu)隧道的形狀特點,其在很大程度上繼承了單圓盾構(gòu)隧道的受力特點,但是由于隧道頂部海鷗型管片在土中具有嚴重的背土效應(yīng),因此其頂部附近土體相比于其他部分會有更大變形,這在實際工程中是不可忽視的問題[4-7]。

    雖然雙圓盾構(gòu)隧道的受力變形規(guī)律與傳統(tǒng)盾構(gòu)隧道相比,會帶來更大的土體位移,更大的土體擾動,但是這仍然是一種未來極具發(fā)展前景的異形盾構(gòu)模式。如今,國內(nèi)外已有許多專家學(xué)者對雙圓盾構(gòu)隧道的工作特點及其給周圍土體帶來的環(huán)境擾動進行分析。孫統(tǒng)立等[8]結(jié)合上海軌道交通8號線的雙圓盾構(gòu)掘進工程,將數(shù)值分析結(jié)果與現(xiàn)場監(jiān)測實際數(shù)值進行比對,發(fā)現(xiàn)雙圓盾構(gòu)的土體擾動位移場與傳統(tǒng)單圓盾構(gòu)相似,但卻具有更大的影響范圍及更高的位移值。陳小亮等[9]依托于實際雙圓盾構(gòu)地鐵隧道區(qū)間工程,應(yīng)用數(shù)值模擬、現(xiàn)場監(jiān)測等方法,建立三維彈塑性有限差分模型計算由雙圓盾構(gòu)施工所引起的土體應(yīng)力分布和地表位移沉降,發(fā)現(xiàn)地表橫向變形影響范圍以及盾首到達和盾尾離開時的沉降特點。劉湘[10]分析對比了雙圓盾構(gòu)隧道與傳統(tǒng)盾構(gòu)隧道在施工過程中對地層、交疊隧道的影響并研究了交疊段施工過程中的影響因素。洪杰[11]研究了盾尾注漿力、正面附加推力,以及側(cè)摩阻力對原本土層上方隧道、鄰近建筑物基礎(chǔ)、地下水電綜合管線的影響。魯漢新[12]提出結(jié)合雙圓盾構(gòu)施工特點的半理論半經(jīng)驗計算模型,并用于地表沉降計算,與實際數(shù)據(jù)對比,證實該計算模型與實測數(shù)據(jù)相對誤差滿足工程要求。

    在土層為軟弱黏土等不利土層地區(qū),無論是單圓盾構(gòu)還是異形盾構(gòu),都不可避免要面對盾構(gòu)端頭加固這一難題,這是盾構(gòu)始發(fā)乃至整個盾構(gòu)法隧道開挖中的重要環(huán)節(jié)。人工凍結(jié)法可以在地面垂直向下打入垂直凍結(jié)管,也可以通過豎井在豎井中水平向土體打入水平凍結(jié)管,通過在凍結(jié)管中循環(huán)通入的冷媒劑的作用,附近含水土層便會形成凍結(jié)區(qū)域以加固土體[13-14]。杯型凍結(jié)壁的形成是通過在工作井中向土體水平打入長度不一的多圈凍結(jié)管,最終形成杯狀的凍土帷幕,目前已經(jīng)在大量工程中得到應(yīng)用[15-21]。如杯型水平凍結(jié)法被創(chuàng)新性地應(yīng)用于上海地鐵2號線某標段,很好地保護了隧道上方的大型污水箱。英旭等[22]系統(tǒng)地描述了凍結(jié)方案設(shè)計、鉆孔施工布置并介紹杯型凍結(jié)壁施工理念。夏江濤等[23]在南京地鐵2號線逸仙橋站溫度場數(shù)值分析實驗中,使用了凍土熱傳導(dǎo)理論來驗證計算方法與數(shù)值模型的正確性。王效賓等[24]將土體材料比熱容及導(dǎo)熱系數(shù)等熱力學(xué)參數(shù)加入到溫度場數(shù)值分析實驗中。王杰[25]在蘇州地鐵2號線火車站東風(fēng)井進出洞數(shù)值分析實驗中,對比分析了數(shù)種盾構(gòu)進出洞端頭土體加固方式,選擇出最優(yōu)解,建立數(shù)值模型,并結(jié)合溫度實測值對比分析。

    目前國內(nèi)將人工凍結(jié)法應(yīng)用于雙圓盾構(gòu)隧道始發(fā)端頭土層加固的工程實例還鮮見報道,本研究以上海市軌道交通6號線某標段雙圓盾構(gòu)隧道為例,將雙圓盾構(gòu)與水平凍結(jié)法結(jié)合起來,以期為今后相關(guān)實際工程和科學(xué)研究提供理論依據(jù)。

    1 凍結(jié)方案設(shè)計

    1.1 工程概況

    本凍結(jié)方案以上海地鐵某雙圓盾構(gòu)隧道掘進工程為例,該掘進工程使用的盾構(gòu)機是土壓平衡盾構(gòu)機,這種盾構(gòu)機有著可以將挖掘出來的土體用于支護結(jié)構(gòu)的特點。該雙圓盾構(gòu)隧道掘進工程由上海市地鐵6號線的濱州路站到成山路站,全長將近2 600 m。穿越粉質(zhì)黏土以及淤泥質(zhì)黏土等多種不利地層,其中隧道經(jīng)過的大部分地層剖面圖如圖1所示。

    根據(jù)地質(zhì)勘探報告,該雙圓盾構(gòu)隧道所穿越土體的物理力學(xué)參數(shù)見表 1。

    1.2 凍結(jié)孔布置

    面向上海地鐵6號線某標段的雙圓盾構(gòu)隧道掘進工程,該雙圓形盾構(gòu)隧道的外直徑為6.2 m、寬度為11.12 m,在盾構(gòu)始發(fā)端頭處,從豎井中向洞門方向打入2根中心凍結(jié)管和3圈外圍凍結(jié)管(內(nèi)圈短管,外圈長管),構(gòu)成“雙圓形”,并循環(huán)通入冷鹽水,最終在盾構(gòu)隧道始發(fā)開挖面附近形成杯型凍土帷幕,對土體進行加固。在有限元數(shù)值模擬軟件中,為了方便布置凍結(jié)管,將最外圈凍結(jié)孔布置成一個外直徑6 m、寬度11 m的“雙圓形”。凍結(jié)孔布置一共102個,從豎井向洞門方向水平布置,其中外圈布置52個凍結(jié)孔(孔徑5 m);中圈布置32個凍結(jié)孔(孔徑2 m);內(nèi)圈布置16個凍結(jié)孔(孔徑2 m);洞門外圈凍結(jié)管構(gòu)成的雙圓形的2個圓心處各布置1個凍結(jié)孔,共2個(孔徑2 m);除中心凍結(jié)孔外共100個外圍凍結(jié)孔都沿著各自所在圓弧均勻布置,各圈之間間距為1 m。所有凍結(jié)管均選用127 mm×5 mm低碳鋼無縫鋼管,凍結(jié)管連接方式選用箍焊連接,全部選用48 mm×3.5 mm鋼管作為供液管;全部102根凍結(jié)管直徑均為0.1 m,如圖2所示。

    1.3 施工工藝流程

    水平凍結(jié)施工工藝難度較低、工作流程簡單、自動化程度高、施工速度快且不影響其他工序。在施工準備完成后,為了節(jié)約時間,可以在土體中開挖凍結(jié)孔的同時,安裝對冷媒介進行降溫的凍結(jié)設(shè)備;凍結(jié)孔施工完成以及凍結(jié)設(shè)備準備就緒之后,可以逐步開始調(diào)試系統(tǒng)設(shè)備,尤其是保溫效果等;然后進行積極凍結(jié),嚴格安排凍結(jié)監(jiān)測;凍結(jié)完成并通過監(jiān)測,確認可以施工后,盾構(gòu)機開始掘進,施工流程如圖3所示。

    2 有限元軟件數(shù)值模型的建立

    2.1 計算基本假定

    為了數(shù)值模型的計算方便以及相對精確性,計算中提出以下基本假定。

    (1)將土體視為均質(zhì)、連續(xù)、各向同性的材料。

    (2)假設(shè)土體原始地溫為18 ℃。

    (3)假設(shè)土體開始凍結(jié)時溫度為-1 ℃,完全凍結(jié)時溫度為-10 ℃。

    (4)忽略水分遷移以及內(nèi)熱阻的影響。

    (5)溫度荷載可以直接作用到凍結(jié)管管壁上。

    2.2 計算模型和參數(shù)選取

    對外徑為6.2 m、寬度為11.12 m的雙圓盾構(gòu)隧道建立三維溫度場數(shù)值模型,為了方便布置凍結(jié)管,將最外圈凍結(jié)孔布置成一個外徑6 m、寬度11 m的“雙圓形”,在考慮凍結(jié)范圍不應(yīng)超過土體三維模型后,將土體幾何尺寸設(shè)置為:以盾構(gòu)隧道開挖面上凍結(jié)孔布置成“雙圓形”的2個圓心連線之間的中點為原點,沿X軸方向為6 m、沿Y軸方向為20 m、沿Z軸方向為10 m的立方體。102根凍結(jié)管插入6 m×20 m×10 m的土體中,如圖4所示。

    根據(jù)地質(zhì)勘探報告顯示,隧道所經(jīng)過土層為粉土、黏土層,但為了模擬最不利情況,選擇傳熱最不利的粉砂、細砂層作為模型土體材料,由于未凍結(jié)狀態(tài)和凍結(jié)狀態(tài)具有不同的密度、導(dǎo)熱系數(shù)、比熱,在參考羅婷等[18]對相似土體進行室內(nèi)試驗與選定的熱物理參數(shù)之后,設(shè)定土體材料參數(shù),見表2。

    假設(shè)凍結(jié)前原始地層溫度為18 ℃,凍結(jié)管壁為熱負荷邊界,邊界負荷是鹽水溫度。設(shè)置全部凍結(jié)時長分為40步,每步時長為一天24 h,一共有960 h,共計40 d;在查閱大量相關(guān)資料以及參考類似實際工程之后,充分考慮實際情況來安排凍結(jié)管降溫速率,最后得出鹽水降溫計劃,見表3[26-27]。凍結(jié)開始時鹽水溫度為原始地溫18 ℃,1 d后降為0 ℃,第5天時降為-15 ℃,從第10天以后一直到凍結(jié)結(jié)束,均降為-28 ℃。

    3 溫度場計算結(jié)果和分析

    3.1 凍土帷幕交圈情況

    根據(jù)圖4的三維模型以及網(wǎng)格劃分,進行數(shù)值計算,計算結(jié)果顯示出凍結(jié)結(jié)束后雙圓盾構(gòu)隧道杯型凍結(jié)壁溫度場的橫截面以及剖面的計算云圖,如圖5所示。在X=0截面上,溫度場呈“雙圓形”,以2根中心凍結(jié)管為圓心,由內(nèi)逐漸向外輻射分布;由最外圈凍結(jié)管形成明顯的分界線,分界線以內(nèi)凍結(jié)溫度平均為-23 ℃,分界線以外隨著距離變遠溫度逐漸升高,在距離最外圈凍結(jié)管0.8 m時溫度升至0 ℃;顯然距離2根中心凍結(jié)管越近,凍結(jié)溫度越低,凍結(jié)效果越好,在最外圈凍結(jié)管內(nèi)已形成完整凍土帷幕。在Z=0截面上,溫度場呈寬口“杯形”,以2根5 m長最外圈凍結(jié)管為杯壁,10根中心及內(nèi)圈2 m長凍結(jié)管為杯底;顯然在杯底范圍內(nèi)以及杯壁附近溫度已經(jīng)降到-23 ℃左右,越遠離凍結(jié)管溫度越高,在距離杯底凍結(jié)管端部0.8 m時溫度已升至0 ℃。

    圖6為杯形凍結(jié)壁X=0剖面不同時間段內(nèi)-1 ℃等溫線圖(紅線表示)。

    由圖6可知,在凍結(jié)初期(第1天到第5天),X=0剖面-1 ℃等溫線基本以凍結(jié)管為圓心呈圓形向外分布,各圈之間沒有接觸,并未形成交圈。在凍結(jié)開始的第6天時,在最外圈凍結(jié)管附近,-1℃等溫線開始交圈,然而中圈和內(nèi)圈以及中心凍結(jié)管均未開始交圈,一直到第8天時,中圈和內(nèi)圈以及中心凍結(jié)管才開始交圈。凍結(jié)開始第9天時,外圈凍結(jié)管附近-1 ℃等溫線已經(jīng)徹底完成交圈。到第12天時,中圈和內(nèi)圈以及中心凍結(jié)管-1 ℃等溫線也完成了交圈,至此,各圈凍結(jié)管均完成交圈,并且隨著時間推移,各圈凍土帷幕開始逐漸與其他凍結(jié)管形成的凍土帷幕相互融合直至形成一道完整的凍土帷幕。到凍結(jié)第18天時,內(nèi)部所有-1 ℃等溫線已經(jīng)完全消失,只有最外圈凍結(jié)管附近形成一圈完整的雙圓形-1 ℃等溫線。

    圖7為杯形凍結(jié)壁X=0剖面不同時間段內(nèi)-10 ℃等溫線圖(藍線表示)。

    由圖7可知,在凍結(jié)初期(以第4天為例),X=0剖面-10 ℃等溫線開始出現(xiàn)并基本以凍結(jié)管為圓心呈圓形向外分布,各圈之間沒有接觸,并未形成交圈。隨著時間推移,-10 ℃等溫線圈徑逐漸增大。在凍結(jié)開始的第10天時,在最外圈凍結(jié)管附近,-10 ℃等溫線開始交圈,持續(xù)到凍結(jié)開始第14天時,外圈凍結(jié)管附近-10 ℃等溫線已經(jīng)徹底完成交圈,中圈和內(nèi)圈以及中心凍結(jié)管此時才開始交圈,到第17天時,中圈和內(nèi)圈以及中心凍結(jié)管-10 ℃等溫線也完成了交圈,至此,各圈凍結(jié)管均完成交圈并且隨著時間推移,各圈凍土帷幕開始逐漸與其他凍結(jié)管形成的凍土帷幕相互融合直至形成一道完整的凍土帷幕。到凍結(jié)第23天時,內(nèi)部所有-10 ℃等溫線已經(jīng)完全消失,只有最外圈凍結(jié)管附近形成一圈完整的雙圓形-10 ℃等溫線。

    此后,杯形凍結(jié)壁X=0剖面-1 ℃等溫線和-10 ℃等溫線交圈所形成的凍土帷幕都隨著時間的推移持續(xù)變大,一直到凍結(jié)終期第40天時,達到峰值。

    3.2 凍結(jié)壁杯底杯身尺寸

    圖8為杯形凍結(jié)壁Z=0剖面-1℃等溫線圖,一般情況下,土體溫度到達-1℃時開始凍結(jié),據(jù)此可以計算出最終在凍結(jié)40 d之后雙圓杯形凍結(jié)壁杯體的杯底、杯身、厚度尺寸:杯底長度為12.4 m,杯身高度為5 m,杯底厚度為2 m,滿足盾構(gòu)隧道始發(fā)條件(盾構(gòu)始發(fā)時凍結(jié)壁杯底厚度條件為大于或等于2 m),杯身厚度為2 m,滿足盾構(gòu)隧道始發(fā)條件(盾構(gòu)始發(fā)時凍結(jié)壁杯身厚度條件為大于或等于1.2 m)。綜上所述,杯形凍結(jié)壁滿足土體加固要求,可以進行盾構(gòu)掘進。

    4 溫度場優(yōu)化分析

    在凍結(jié)進行到第18天時已形成完整凍土帷幕,因此在考慮經(jīng)濟效應(yīng)等因素的前提下,經(jīng)過多次篩選減少凍結(jié)管數(shù)量的優(yōu)化方案,選出2組凍結(jié)方案來進行優(yōu)化,見表4。方案1為內(nèi)圈、中圈、外圈凍結(jié)管根數(shù)分別減少2、4、6根;方案2為內(nèi)圈、中圈、外圈凍結(jié)管根數(shù)分別減少4、6、10根,研究不同方案下溫度場的分布以及發(fā)展規(guī)律。

    4.1 方案1溫度場數(shù)值模擬

    保持優(yōu)化方案1與初始方案的其他尺寸參數(shù)不變,內(nèi)圈、中圈、外圈凍結(jié)管根數(shù)分別減少2、4、6根,觀察優(yōu)化方案1下溫度場分布情況以及發(fā)展規(guī)律,概括分析凍土帷幕的交圈情況等,使用圖形轉(zhuǎn)換法計算出杯型凍結(jié)壁的杯身厚度、高度;杯底厚度、長度等相關(guān)指標參數(shù),并與設(shè)計值進行對比分析來判斷其是否滿足盾構(gòu)始發(fā)必備條件。

    凍結(jié)開始到第10天、第14天以及第20天時,X=0剖面與Z=0剖面上凍土帷幕的交圈情況如圖9所示(紅線表示-1 ℃等溫線,藍線表示-10 ℃等溫線)。在使用凍結(jié)優(yōu)化方案1中,凍土帷幕的變化

    發(fā)展和交圈情況,以及最后形成的杯型凍結(jié)壁溫度場分布及發(fā)展規(guī)律與優(yōu)化前初始模型均基本一致,但是凍土帷幕外圈交圈、中圈交圈、內(nèi)圈交圈、凍土帷幕封閉等步驟所需時間較初始模型均有增加。初始模型開始凍結(jié)第6天外圈交圈,優(yōu)化方案1到第10天外圈交圈;初始模型開始凍結(jié)第8天中圈及內(nèi)圈交圈,優(yōu)化方案1到第14天中圈及內(nèi)圈交圈;初始模型開始凍結(jié)第18天凍土帷幕封閉,優(yōu)化方案1到第20天凍土帷幕封閉。

    一般情況下,土體溫度到達-1 ℃時含水地層中的土體開始凍結(jié),將-1 ℃等溫線內(nèi)的土體視為已凍結(jié)土體。將圖10的凍結(jié)優(yōu)化方案1中凍結(jié)結(jié)束后形成的雙圓杯型凍結(jié)壁進行圖形轉(zhuǎn)換,計算出在凍結(jié)40 d之后雙圓杯型凍結(jié)壁杯體的杯底、杯身、厚度尺寸。經(jīng)計算,雙圓杯型凍結(jié)壁杯底在2個圓心所在連線上的總長度為12.3 m;最外圈凍結(jié)管構(gòu)成的杯身高度5.1 m,杯壁厚度1.8 m;內(nèi)圈中圈及中心凍結(jié)管構(gòu)成的杯底厚度2.2 m,在與盾構(gòu)隧道開始掘進時所需的杯身、杯底厚度最低限值進行對比后,發(fā)現(xiàn)凍結(jié)優(yōu)化方案1模型計算結(jié)果中雙圓杯型凍結(jié)壁滿足土體加固要求,不影響施工安全,可以進行盾構(gòu)掘進。

    4.2 方案2溫度場數(shù)值模擬

    保持優(yōu)化方案2與原始凍結(jié)方案其他尺寸參數(shù)不變,內(nèi)圈、中圈、外圈凍結(jié)管根數(shù)分別減少4、6、10,觀察優(yōu)化方案2下溫度場分布情況以及發(fā)展規(guī)律,分析凍土帷幕的交圈情況等,使用圖形轉(zhuǎn)換法來計算出杯型凍結(jié)壁的杯身厚度和高度、杯底厚度和長度等相關(guān)指標參數(shù),并與設(shè)計值進行對比分析來判斷其是否滿足盾構(gòu)始發(fā)必備條件。

    凍結(jié)開始到第11天、第16天以及第21天時,X=0剖面與Z=0剖面上凍土帷幕的交圈情況如圖11所示,以紅線表示-1 ℃等溫線,以藍線表示-10 ℃等溫線。

    由圖11可知,在使用凍結(jié)優(yōu)化方案2中,凍土帷幕的變化發(fā)展和交圈情況,以及最后形成的杯型凍結(jié)壁溫度場分布、發(fā)展規(guī)律與原模型均基本一致,但是凍土帷幕外圈交圈、中圈交圈、內(nèi)圈交圈和凍土帷幕封閉等步驟所需時間較優(yōu)化前初始模型均有增加。初始模型開始凍結(jié)第6天外圈交圈和優(yōu)化方案2到第11天外圈交圈;初始模型開始凍結(jié)第8天中圈及內(nèi)圈交圈,優(yōu)化方案2到第16天中圈及內(nèi)圈交圈;初始模型開始凍結(jié)第18天凍土帷幕封閉,優(yōu)化方案2到第21天凍土帷幕封閉,各步驟所需時間均略有增加。

    將圖12中凍結(jié)方案2模型凍結(jié)結(jié)束后形成的雙圓杯型凍結(jié)壁進行圖形轉(zhuǎn)換,計算出在凍結(jié)40 d之后雙圓杯型凍結(jié)壁杯體的杯底、杯身和厚度尺寸。經(jīng)計算,雙圓杯型凍結(jié)壁杯底在2個圓心所在連線上的總長度為12.4 m;最外圈凍結(jié)管構(gòu)成的杯身高度5.1 m,杯壁厚度1.2 m;內(nèi)圈中圈及中心凍結(jié)管構(gòu)成的杯底厚度2.3 m,在與盾構(gòu)隧道開始掘進時所需的杯身杯底厚度最低限值進行對比后,發(fā)現(xiàn)凍結(jié)優(yōu)化方案2模型計算結(jié)果中雙圓杯型凍結(jié)壁滿足土體加固要求,不影響施工安全,可以進行盾構(gòu)掘進。

    4.3 不同方案下凍結(jié)壁凍結(jié)狀況比較

    對比初始模型以及優(yōu)化方案1、優(yōu)化方案2,溫度場數(shù)值計算結(jié)果見表5。由表5可知,相比較于初始模型,適當減少凍結(jié)管數(shù)量后,杯型凍結(jié)壁杯底、杯身和杯底厚度等尺寸仍然滿足盾構(gòu)始發(fā)要求(杯底厚度大于等于2.0 m;杯體厚度大于等于1.2 m;

    凍土的平均溫度小于等于-10 ℃;各探孔溫度小于等于-2 ℃;鹽水溫度-25~-28 ℃)。因此,在充分考慮經(jīng)濟效益的前提下,可以適當減少凍結(jié)管布置根數(shù),在初始模型上可以從外圈、中圈、內(nèi)圈分別減少4、6、10根來進行杯型凍結(jié)加固。

    5 結(jié)論

    本文通過ADINA有限元分析軟件溫度場模塊,基于雙圓盾構(gòu)隧道始發(fā)端頭使用杯型凍結(jié)壁水平凍結(jié)法對土體進行加固,對該雙圓形凍結(jié)壁溫度場的發(fā)展規(guī)律進行數(shù)值分析,并對比分析2種優(yōu)化方案。主要得出以下結(jié)論。

    (1)距離雙圓型盾構(gòu)面2個圓心越近凍結(jié)效果越好,距離增大效果變差;凍結(jié)結(jié)束后,形成一個形狀規(guī)則、分布均勻的近似雙圓“杯形”凍結(jié)壁,且強度分布均勻。

    (2)在凍結(jié)初期,凍土帷幕圍繞各凍結(jié)管呈圓柱形逐漸向外延伸。開始交圈的先后次序為:外圈管、中圈管、內(nèi)圈管,在優(yōu)化前初始凍結(jié)方案下,形成閉合的-1 ℃和-10 ℃等溫線分別為凍結(jié)開始第18天和第23天,杯型凍結(jié)壁交圈時間為凍結(jié)第18天。

    (3)凍結(jié)結(jié)束之后,原始模型雙圓杯型凍結(jié)壁的杯底、杯身、厚度尺寸:長度為12.4? m,杯身高度為5 m,厚度為2 m,杯底厚度為2 m。

    (4)相對于始發(fā)基座上的盾構(gòu)端頭推入加固土體,進入地層原狀土區(qū)段的始發(fā)凍結(jié)條件,優(yōu)化前初始模型設(shè)計凍結(jié)方案符合要求,凍結(jié)效果良好,在實際施工中具有可行性,但偏于保守,不利于成本控制。在適量減少凍結(jié)管的情況下杯型凍結(jié)壁的厚度等參數(shù)依然滿足盾構(gòu)始發(fā)條件,因此建議實際工程中適量減少凍結(jié)管根數(shù),以達到更高的經(jīng)濟效益。

    【參 考 文 獻】

    [1]胡俊.高水壓砂性土層地鐵大直徑盾構(gòu)始發(fā)端頭加固方式研究[D].南京:南京林業(yè)大學(xué),2012.

    HU J. Study on the reinforcement methods of subway large-diameter shield launching in the sandy clay with high water pressure[D]. Nanjing: Nanjing Forestry University, 2012.

    [2]胡俊,楊平,董朝文,等.盾構(gòu)始發(fā)端頭化學(xué)加固范圍及加固工藝研究[J].鐵道建筑,2010,50(2):47-51.

    HU J, YANG P, DONG (C /Z) W, et al. Study on chemical reinforcement range and reinforcement technology of shield starting end[J]. Railway Engineering, 2010, 50(2): 47-51.

    [3]HU J, ZENG H, WANG X B. Numerical analysis of temperature field of cup-shaped frozen soil wall reinforcement at shield shaft[J]. Applied Mechanics and Materials, 2013, 341/342: 1467-1471.

    [4]溫瑜琴,蘇棟,鄧碧,等.雙圓盾構(gòu)掘進地層擾動的三維數(shù)值模擬[J].現(xiàn)代隧道技術(shù),2020,57(S1):450-457.

    WEN Y Q, SU D, DENG B, et al. Three dimensional numerical simulation of ground disturbance caused by DOT shield tunnelling[J]. Modern Tunnelling Technology, 2020, 57(S1): 450-457.

    [5]周文波.雙圓(DOT)隧道技術(shù)及應(yīng)用[J].都市快軌交通,2005,18(2):33-38.

    ZHOU W B. DOT technology and its application[J]. Urban Rapid Rail Transit, 2005, 18(2): 33-38.

    [6]管攀峰.上海軌道交通12號線工程區(qū)間隧道設(shè)計優(yōu)化[J].城市軌道交通研究,2017,20(1):137-143.

    GUAN P F. Design optimization of the running tunnel on Shanghai rail transit line 12[J]. Urban Mass Transit, 2017, 20(1): 137-143.

    [7]汪宏,劉湘.過江雙圓盾構(gòu)隧道施工過程的分析和研究[J].中國水運(下半月),2015,15(4):216-218.

    WANG H, LIU X. Analysis and research on construction process of double round shield tunnel crossing the river[J]. China Water Transport, 2015, 15(4): 216-218.

    [8]孫統(tǒng)立,張海濤.雙圓盾構(gòu)隧道施工力學(xué)行為三維數(shù)值模擬[J].城市軌道交通研究,2013,16(5):36-40.

    SUN T L, ZHANG H T. 3D numerical simulation of the mechanical behavior of double-O-tube shield tunneling[J]. Urban Mass Transit, 2013, 16(5): 36-40.

    [9]陳小亮,周松國,朱世煜.軟土地層雙圓盾構(gòu)掘進對周圍土體擾動分析[J].城市道橋與防洪,2019(5):232-236,27.

    CHEN X L, ZHOU S G, ZHU S Y. Analysis of soil disturbance caused by tunneling of double-circle shield tunnel (DOT) in soft soil stratum[J]. Urban Roads Bridges & Flood Control, 2019(5): 232-236, 27.

    [10]劉湘.雙圓盾構(gòu)施工對交疊隧道影響的研究[D].鎮(zhèn)江:江蘇科技大學(xué),2015.

    LIU X. Study on influence of DOT shield excavation on overlapping tunnel[D]. Zhenjiang: Jiangsu University of Science and Technology, 2015.

    [11]洪杰.雙圓盾構(gòu)隧道施工擾動及對周邊構(gòu)筑物影響研究[D].杭州:浙江大學(xué),2013.

    HONG J. Study on DOT shield tunnel construction disturbance and its impact on surrounding structures[D]. Hangzhou: Zhejiang University, 2013.

    [12]魯漢新.雙圓盾構(gòu)法隧道地表沉降的半理論半經(jīng)驗預(yù)測[J].城市道橋與防洪,2013(6):272-274,282.

    LU H X. Semi-empirical and semi-theoretical prediction method of surface settlement caused by DOT tunneling[J]. Urban Roads Bridges & Flood Control, 2013(6): 272-274, 282.

    [13]HU J, WANG X B, JIANG B R. Numerical analysis of temperature field of vertical frozen soil wall reinforcement at shield shaft[J]. Advanced Materials Research, 2014, 918: 218-223.

    [14]HU J, JIANG B S, ZENG H, et al. Temperature field numerical analysis of different freeze pipe spacing of vertical frozen soil wall reinforcement at shield shaft[J]. Applied Mechanics and Materials, 2014, 580-583: 738-741.

    [15]胡俊,楊平.大直徑杯型凍土壁溫度場數(shù)值分析[J].巖土力學(xué),2015,36(2):523-531.

    HU J, YANG P. Numerical analysis of temperature field within large-diameter cup-shaped frozen soil wall[J]. Rock and Soil Mechanics, 2015, 36(2): 523-531.

    [16]曾華.盾構(gòu)進出洞處地基水平杯型凍結(jié)加固的溫度場研究[J].城市軌道交通研究,2016,19(6):112-116.

    ZENG H. On frozen temperature field of horizontal cup-type freezing technology in shield entry engineering[J]. Urban Mass Transit, 2016, 19(6): 112-116.

    [17]胡俊,張皖湘,曾暉.盾構(gòu)隧道端頭杯型凍結(jié)壁加固溫度場數(shù)值分析[J].路基工程,2015(4):20-22.

    HU J, ZHANG W X, ZENG H. Numerical analysis of temperature field at cup-shaped frozen wall reinforcement at the end of shield tunnel[J]. Subgrade Engineering, 2015(4): 20-22.

    [18]羅婷,胡俊,衛(wèi)宏.盾構(gòu)隧道端頭杯型凍結(jié)壁溫度場數(shù)值優(yōu)化分析[J].森林工程,2017,33(4):83-88.

    LUO T, HU J, WEI H. Numerical optimization analysis of temperature field with shield tunnelling cup-shaped freezing wall[J]. Forest Engineering, 2017, 33(4): 83-88.

    [19]周雋,胡俊,衛(wèi)宏,等.盾構(gòu)隧道端頭垂直杯型凍結(jié)壁溫度場數(shù)值分析[J].森林工程,2017,33(3):74-79.

    ZHOU J, HU J, WEI H, et al. Numerical analysis of shield tunnel end’s vertical cup-shaped and frozen wall’s temperature field[J]. Forest Engineering, 2017, 33(3): 74-79.

    [20]付財,崔兵兵,王翔.水泥土杯型凍結(jié)加固盾構(gòu)隧道端頭溫度場發(fā)展規(guī)律研究[J].山西建筑,2020,46(16):125-127.

    FU C, CUI B B, WANG X. Research on the development rule of temperature field at the end of shield tunnel strengthened by cement modified soil cup[J]. Shanxi Architecture, 2020, 46(16): 125-127.

    [21]韓霖磊,張孟喜,吳惠明,等.盾構(gòu)隧道施工誘發(fā)地層變形特征的透明土模型試驗[J].公路工程,2021,46(4):46-51.

    HAN L L, ZHANG M X, WU H M, et al. Model test on subsurface ground settlement caused by shield tunnel construction based on transparent soil[J]. Highway Engineering, 2021, 46(4): 46-51.

    [22]英旭,蔣岳成,李曦.杯型水平凍結(jié)工法在盾構(gòu)進洞施工中的應(yīng)用[J].中國市政工程,2006(2):52-55,94.

    YING X, JIANG Y C, LI X. “Cup-shaped” horizontal freezing method for shield machine driving[J]. China Municipal Engineering, 2006(2): 52-55, 94.

    [23]夏江濤,楊平.盾構(gòu)出洞水平凍結(jié)加固杯型凍土壁溫度場數(shù)值分析[J].河南科技大學(xué)學(xué)報(自然科學(xué)版),2010,31(3):58-62,110.

    XIA J T, YANG P. Numerical analysis of temperature field in cup-shaped frozen soil wall of horizontal freezing reinforcement for shield tunneling setting out[J]. Journal of Henan University of Science & Technology (Natural Science), 2010, 31(3): 58-62, 110.

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

    WANG X B, YANG P, ZHANG T, et al. 3-D finite element analysis of melting temperature field in field tunneling horizontal freezing reinforcing engineering[J]. Journal of PLA University of Science and Technology (Natural Science Edition), 2009, 10(6): 586-590.

    [25]王杰.盾構(gòu)進出洞水泥土加固后水平凍結(jié)溫度場研究[D].南京:南京林業(yè)大學(xué),2012.

    WANG J. Numerical study of soil thermal field in horizontal freezing shield portal reinforced by cement[D]. Nanjing: Nanjing Forestry University, 2012.

    [26]胡俊,曾暉,佳琳.一種大直徑盾構(gòu)隧道端頭杯型水平凍結(jié)孔布設(shè)方法:CN104405400A[P].2015-03-11.

    HU J, ZENG H, JIA L. Method for arranging cup-shaped horizontal freezing holes in end of large-diameter shield tunnel: CN104405400A[P]. 2015-03-11.

    [27]胡俊,唐益群,周潔.一種盾構(gòu)隧道端頭垂直杯型凍結(jié)加固結(jié)構(gòu):CN204609889U[P].2015-09-02.

    HU J, TANG Y Q, ZHOU J. Perpendicular cup type stabilization by freezing structure in shield tunnel end: CN204609889U[P]. 2015-09-02.

    猜你喜歡
    數(shù)值模擬隧道
    與隧道同行
    預(yù)見2019:隧道的微光
    中國外匯(2019年6期)2019-07-13 05:44:02
    神奇的泥巴山隧道
    小讀者(2019年24期)2019-01-10 23:00:37
    黑乎乎的隧道好可怕
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    一種基于液壓緩沖的減震管卡設(shè)計與性能分析
    科技視界(2016年20期)2016-09-29 11:08:27
    男女免费视频国产| 香蕉国产在线看| 99久久中文字幕三级久久日本| 精品一区在线观看国产| 少妇 在线观看| 新久久久久国产一级毛片| 伦精品一区二区三区| 91午夜精品亚洲一区二区三区| 只有这里有精品99| 亚洲av电影在线观看一区二区三区| 亚洲综合精品二区| 少妇精品久久久久久久| 国产成人精品福利久久| 国产极品粉嫩免费观看在线| 亚洲精品久久午夜乱码| 麻豆精品久久久久久蜜桃| av线在线观看网站| 久久久久久久久久人人人人人人| 亚洲美女视频黄频| av卡一久久| 丝袜在线中文字幕| 捣出白浆h1v1| 日本wwww免费看| 精品少妇久久久久久888优播| 亚洲精品成人av观看孕妇| 一区福利在线观看| 伊人亚洲综合成人网| 欧美精品国产亚洲| 老司机影院成人| videos熟女内射| 色94色欧美一区二区| 日韩精品有码人妻一区| 国产精品 欧美亚洲| 国产一区二区 视频在线| 伊人久久国产一区二区| 久久精品夜色国产| 少妇人妻精品综合一区二区| 久久精品夜色国产| 国产日韩欧美视频二区| av线在线观看网站| 国产淫语在线视频| 亚洲视频免费观看视频| 国产精品二区激情视频| 国产av精品麻豆| 老鸭窝网址在线观看| 亚洲人成电影观看| 日本色播在线视频| 最近最新中文字幕免费大全7| 99久久中文字幕三级久久日本| 中文字幕最新亚洲高清| 男男h啪啪无遮挡| 春色校园在线视频观看| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av天美| 老熟女久久久| 天天影视国产精品| 亚洲国产精品999| 国产成人精品在线电影| 久久久久国产一级毛片高清牌| 久久精品国产综合久久久| 老熟女久久久| 日本-黄色视频高清免费观看| 桃花免费在线播放| 热99国产精品久久久久久7| 亚洲第一av免费看| 男女下面插进去视频免费观看| xxx大片免费视频| 亚洲四区av| 国产xxxxx性猛交| 亚洲一码二码三码区别大吗| 欧美中文综合在线视频| 人妻系列 视频| 久久毛片免费看一区二区三区| 老鸭窝网址在线观看| 久久久精品94久久精品| 久久久久久免费高清国产稀缺| 午夜福利,免费看| 香蕉国产在线看| 久久精品久久精品一区二区三区| 九九爱精品视频在线观看| av有码第一页| 成人国语在线视频| 国产精品久久久久久精品电影小说| 久久婷婷青草| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 秋霞在线观看毛片| 欧美激情 高清一区二区三区| 黄片无遮挡物在线观看| 高清不卡的av网站| 精品亚洲成国产av| 高清黄色对白视频在线免费看| 欧美日韩亚洲国产一区二区在线观看 | 秋霞在线观看毛片| 99国产精品免费福利视频| 欧美人与性动交α欧美精品济南到 | 人妻系列 视频| 亚洲成av片中文字幕在线观看 | 午夜免费鲁丝| 久久青草综合色| 9色porny在线观看| 少妇猛男粗大的猛烈进出视频| 最近2019中文字幕mv第一页| 热re99久久精品国产66热6| 国产乱人偷精品视频| 一区福利在线观看| 精品酒店卫生间| 永久网站在线| 亚洲天堂av无毛| 亚洲欧美清纯卡通| 国产精品久久久久久精品电影小说| 国产成人欧美| 久久国产精品男人的天堂亚洲| 岛国毛片在线播放| 如日韩欧美国产精品一区二区三区| 中国国产av一级| 一二三四在线观看免费中文在| 欧美97在线视频| 99re6热这里在线精品视频| 麻豆av在线久日| 七月丁香在线播放| 一二三四在线观看免费中文在| 久久av网站| 欧美人与性动交α欧美精品济南到 | 晚上一个人看的免费电影| 在线精品无人区一区二区三| 精品久久久久久电影网| 日韩 亚洲 欧美在线| 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 亚洲国产最新在线播放| 久久久欧美国产精品| 免费人妻精品一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 日本91视频免费播放| 欧美日韩一区二区视频在线观看视频在线| 久久午夜综合久久蜜桃| 熟妇人妻不卡中文字幕| 国产有黄有色有爽视频| 91久久精品国产一区二区三区| 亚洲精品中文字幕在线视频| 丝袜美腿诱惑在线| 亚洲综合色惰| 少妇熟女欧美另类| 90打野战视频偷拍视频| 久久久精品94久久精品| 黑人巨大精品欧美一区二区蜜桃| 电影成人av| 日本猛色少妇xxxxx猛交久久| 2021少妇久久久久久久久久久| 久久久久久久精品精品| 国语对白做爰xxxⅹ性视频网站| 国产亚洲av片在线观看秒播厂| 国产白丝娇喘喷水9色精品| 国产麻豆69| 日日爽夜夜爽网站| 侵犯人妻中文字幕一二三四区| 亚洲三级黄色毛片| 成人亚洲精品一区在线观看| 麻豆av在线久日| 精品久久蜜臀av无| 校园人妻丝袜中文字幕| 下体分泌物呈黄色| av线在线观看网站| 看免费成人av毛片| 只有这里有精品99| 有码 亚洲区| 不卡av一区二区三区| 国产成人欧美| 中文字幕人妻丝袜一区二区 | freevideosex欧美| 日韩中文字幕视频在线看片| av卡一久久| 国产爽快片一区二区三区| 日韩精品免费视频一区二区三区| 国产成人精品在线电影| 精品国产一区二区久久| 最近中文字幕高清免费大全6| 老汉色∧v一级毛片| 国产免费又黄又爽又色| 啦啦啦啦在线视频资源| 极品人妻少妇av视频| 午夜福利网站1000一区二区三区| 国产精品亚洲av一区麻豆 | 久久人人爽av亚洲精品天堂| 免费黄色在线免费观看| 亚洲欧美成人精品一区二区| 91aial.com中文字幕在线观看| 一区二区三区四区激情视频| 99热全是精品| 国产一区二区三区综合在线观看| 色婷婷av一区二区三区视频| 日本欧美国产在线视频| 久久毛片免费看一区二区三区| 亚洲精品久久久久久婷婷小说| 黄片无遮挡物在线观看| 一本色道久久久久久精品综合| 曰老女人黄片| 亚洲精品中文字幕在线视频| 亚洲国产成人一精品久久久| 久久精品人人爽人人爽视色| 人人妻人人澡人人看| 亚洲精品av麻豆狂野| 熟妇人妻不卡中文字幕| 亚洲av电影在线进入| 亚洲成国产人片在线观看| 亚洲一区中文字幕在线| 国产成人午夜福利电影在线观看| 黄片无遮挡物在线观看| 少妇猛男粗大的猛烈进出视频| 精品亚洲成a人片在线观看| 欧美激情高清一区二区三区 | 中文字幕人妻熟女乱码| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡 | 熟女电影av网| 日韩中文字幕视频在线看片| 午夜激情av网站| 午夜福利一区二区在线看| 美女大奶头黄色视频| 国产激情久久老熟女| 一级黄片播放器| www.av在线官网国产| 少妇人妻久久综合中文| 亚洲欧美一区二区三区久久| 男人舔女人的私密视频| 看十八女毛片水多多多| 精品国产一区二区久久| 午夜老司机福利剧场| 免费少妇av软件| 视频在线观看一区二区三区| 免费av中文字幕在线| 99久久综合免费| 亚洲精品中文字幕在线视频| 亚洲精品一二三| 丝袜人妻中文字幕| 国产精品三级大全| 国产黄色视频一区二区在线观看| 97在线人人人人妻| www.av在线官网国产| 亚洲国产av影院在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一区蜜桃| 久久精品人人爽人人爽视色| 婷婷色综合大香蕉| 丝袜美足系列| 久久久a久久爽久久v久久| 亚洲欧美一区二区三区黑人 | 91成人精品电影| 老鸭窝网址在线观看| 亚洲精品av麻豆狂野| videos熟女内射| 高清黄色对白视频在线免费看| www.精华液| 青春草国产在线视频| 国产在线视频一区二区| 热99国产精品久久久久久7| 爱豆传媒免费全集在线观看| 午夜福利影视在线免费观看| 精品亚洲乱码少妇综合久久| 在线天堂最新版资源| www.精华液| 亚洲欧美一区二区三区国产| 成人免费观看视频高清| av网站在线播放免费| 亚洲精品中文字幕在线视频| 欧美激情高清一区二区三区 | 国产亚洲午夜精品一区二区久久| 老汉色∧v一级毛片| 国产高清不卡午夜福利| 日韩大片免费观看网站| 校园人妻丝袜中文字幕| 九九爱精品视频在线观看| av.在线天堂| 久久久久人妻精品一区果冻| 亚洲av在线观看美女高潮| 欧美xxⅹ黑人| 一本色道久久久久久精品综合| 热99国产精品久久久久久7| 久久精品aⅴ一区二区三区四区 | 成人国语在线视频| 丝袜美足系列| 日本91视频免费播放| 国产欧美日韩综合在线一区二区| 久久午夜综合久久蜜桃| 久久精品夜色国产| 不卡视频在线观看欧美| 人人妻人人澡人人看| videossex国产| 国产又爽黄色视频| 男女免费视频国产| 日韩av免费高清视频| 不卡av一区二区三区| 国产免费又黄又爽又色| 亚洲人成电影观看| 亚洲精品久久午夜乱码| 国产一级毛片在线| 亚洲av国产av综合av卡| 免费黄色在线免费观看| 高清黄色对白视频在线免费看| 久久精品国产综合久久久| 久久鲁丝午夜福利片| 国产成人精品婷婷| 高清在线视频一区二区三区| 少妇人妻 视频| 免费看av在线观看网站| 亚洲av中文av极速乱| 亚洲一区中文字幕在线| 亚洲一码二码三码区别大吗| 男人爽女人下面视频在线观看| 国产男女超爽视频在线观看| 久久99热这里只频精品6学生| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 日韩一区二区视频免费看| 最近最新中文字幕免费大全7| 熟妇人妻不卡中文字幕| 亚洲天堂av无毛| 建设人人有责人人尽责人人享有的| 成人毛片60女人毛片免费| 国产在视频线精品| 肉色欧美久久久久久久蜜桃| 国产在视频线精品| av.在线天堂| 亚洲国产最新在线播放| 叶爱在线成人免费视频播放| 国产成人免费无遮挡视频| 国产探花极品一区二区| 黄片播放在线免费| 免费在线观看完整版高清| av一本久久久久| 精品一品国产午夜福利视频| 国产男女内射视频| 成人午夜精彩视频在线观看| √禁漫天堂资源中文www| 欧美在线黄色| 一边亲一边摸免费视频| 男人添女人高潮全过程视频| 日韩中文字幕视频在线看片| 美女午夜性视频免费| 成年女人毛片免费观看观看9 | 黑人猛操日本美女一级片| 波多野结衣av一区二区av| 亚洲伊人久久精品综合| 久久精品国产亚洲av天美| 国产极品粉嫩免费观看在线| 黑人巨大精品欧美一区二区蜜桃| 极品人妻少妇av视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人久久精品综合| 国产深夜福利视频在线观看| 色婷婷av一区二区三区视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产av新网站| 91成人精品电影| 国产成人91sexporn| 天天躁夜夜躁狠狠躁躁| 搡女人真爽免费视频火全软件| 18禁观看日本| 亚洲欧美成人综合另类久久久| 精品99又大又爽又粗少妇毛片| 亚洲天堂av无毛| 999精品在线视频| 中文欧美无线码| 秋霞在线观看毛片| 制服诱惑二区| 成年av动漫网址| 99香蕉大伊视频| 晚上一个人看的免费电影| 免费观看a级毛片全部| 精品国产超薄肉色丝袜足j| 亚洲精品,欧美精品| 国产精品一二三区在线看| 最新的欧美精品一区二区| 在线免费观看不下载黄p国产| 成人毛片a级毛片在线播放| www日本在线高清视频| 狠狠婷婷综合久久久久久88av| 咕卡用的链子| 免费观看av网站的网址| 男人爽女人下面视频在线观看| 人成视频在线观看免费观看| 97人妻天天添夜夜摸| 中文天堂在线官网| 最近2019中文字幕mv第一页| 一区二区三区精品91| 国产精品国产三级国产专区5o| a级毛片黄视频| 欧美精品人与动牲交sv欧美| 亚洲欧美精品综合一区二区三区 | 日韩免费高清中文字幕av| 黑丝袜美女国产一区| av电影中文网址| 欧美日韩一区二区视频在线观看视频在线| 久久久久网色| 日韩av不卡免费在线播放| 九色亚洲精品在线播放| 国产精品国产三级专区第一集| 天堂俺去俺来也www色官网| 日本av免费视频播放| 女人被躁到高潮嗷嗷叫费观| 少妇人妻久久综合中文| 欧美97在线视频| 国产亚洲午夜精品一区二区久久| 欧美激情 高清一区二区三区| 欧美日韩视频高清一区二区三区二| 亚洲男人天堂网一区| 国产激情久久老熟女| 亚洲国产精品国产精品| 最新的欧美精品一区二区| 看免费av毛片| 午夜福利一区二区在线看| av片东京热男人的天堂| 国产一级毛片在线| 丝袜人妻中文字幕| 国产精品一区二区在线不卡| 人妻人人澡人人爽人人| 欧美日本中文国产一区发布| 啦啦啦在线观看免费高清www| 午夜福利视频在线观看免费| 国产成人免费观看mmmm| 韩国高清视频一区二区三区| 亚洲成人av在线免费| 高清不卡的av网站| 欧美黄色片欧美黄色片| 日本猛色少妇xxxxx猛交久久| 国产精品麻豆人妻色哟哟久久| 欧美人与性动交α欧美精品济南到 | av国产久精品久网站免费入址| 亚洲成色77777| 精品少妇内射三级| 免费少妇av软件| 18禁裸乳无遮挡动漫免费视频| 十八禁高潮呻吟视频| 国产极品粉嫩免费观看在线| 香蕉精品网在线| 大香蕉久久成人网| 亚洲精品成人av观看孕妇| 成人免费观看视频高清| 黑人欧美特级aaaaaa片| 赤兔流量卡办理| 精品人妻在线不人妻| 精品国产乱码久久久久久小说| 777久久人妻少妇嫩草av网站| 老司机影院毛片| 精品第一国产精品| 精品人妻在线不人妻| 大码成人一级视频| 男人舔女人的私密视频| 街头女战士在线观看网站| 天天躁日日躁夜夜躁夜夜| 久久亚洲国产成人精品v| av网站在线播放免费| 建设人人有责人人尽责人人享有的| 国产福利在线免费观看视频| 国产成人a∨麻豆精品| 在线精品无人区一区二区三| 国产亚洲欧美精品永久| 妹子高潮喷水视频| 精品人妻偷拍中文字幕| 啦啦啦中文免费视频观看日本| 新久久久久国产一级毛片| 最新中文字幕久久久久| 最近2019中文字幕mv第一页| 国产在线免费精品| 丰满迷人的少妇在线观看| 如何舔出高潮| 制服丝袜香蕉在线| 久久国产亚洲av麻豆专区| 999久久久国产精品视频| 亚洲激情五月婷婷啪啪| 国产一区二区 视频在线| 久久久国产欧美日韩av| 久久99蜜桃精品久久| 色播在线永久视频| 国产午夜精品一二区理论片| 亚洲av成人精品一二三区| 亚洲一区中文字幕在线| 美女高潮到喷水免费观看| 国产精品免费大片| 亚洲国产成人一精品久久久| 国产av精品麻豆| 天堂8中文在线网| 男人添女人高潮全过程视频| 男人操女人黄网站| 在线观看免费高清a一片| 久热这里只有精品99| 一区二区av电影网| 最近的中文字幕免费完整| 亚洲经典国产精华液单| 制服人妻中文乱码| 国产精品国产av在线观看| 国产精品一区二区在线不卡| 欧美黄色片欧美黄色片| 欧美日韩视频高清一区二区三区二| 亚洲男人天堂网一区| 久久久久国产一级毛片高清牌| 精品亚洲成a人片在线观看| 少妇的丰满在线观看| 欧美97在线视频| 亚洲av.av天堂| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 美国免费a级毛片| 成年美女黄网站色视频大全免费| 女性被躁到高潮视频| 男女免费视频国产| 免费日韩欧美在线观看| 狠狠婷婷综合久久久久久88av| av国产精品久久久久影院| 日韩熟女老妇一区二区性免费视频| 啦啦啦在线免费观看视频4| 免费在线观看完整版高清| 亚洲av福利一区| 伊人久久大香线蕉亚洲五| 亚洲美女黄色视频免费看| 国产精品99久久99久久久不卡 | 美女午夜性视频免费| 大码成人一级视频| 91精品三级在线观看| 一本大道久久a久久精品| 1024视频免费在线观看| 亚洲图色成人| 黑人猛操日本美女一级片| 黄色一级大片看看| 成年人午夜在线观看视频| 五月伊人婷婷丁香| 久久精品熟女亚洲av麻豆精品| 黄色怎么调成土黄色| 成年美女黄网站色视频大全免费| 婷婷成人精品国产| 亚洲综合色惰| 亚洲国产精品999| 亚洲欧美清纯卡通| 香蕉国产在线看| 一级爰片在线观看| 久久青草综合色| 午夜福利网站1000一区二区三区| 中文字幕制服av| 日韩av免费高清视频| 精品酒店卫生间| 欧美精品亚洲一区二区| 国产精品秋霞免费鲁丝片| 国产午夜精品一二区理论片| av国产久精品久网站免费入址| 最新的欧美精品一区二区| 极品少妇高潮喷水抽搐| 亚洲情色 制服丝袜| 黄频高清免费视频| 欧美亚洲日本最大视频资源| 国产男人的电影天堂91| 国产成人精品福利久久| 最近的中文字幕免费完整| 亚洲精品在线美女| 国产一区二区三区av在线| 久久国产精品男人的天堂亚洲| 亚洲色图综合在线观看| 中文天堂在线官网| 国产一级毛片在线| 国产男女内射视频| 波野结衣二区三区在线| 韩国av在线不卡| 亚洲成国产人片在线观看| 十八禁网站网址无遮挡| 亚洲欧美中文字幕日韩二区| 久久久精品国产亚洲av高清涩受| 1024香蕉在线观看| 国产精品久久久久成人av| 超色免费av| 国产日韩欧美在线精品| 啦啦啦视频在线资源免费观看| 一级,二级,三级黄色视频| av线在线观看网站| 久久热在线av| 国产成人精品久久久久久| 亚洲精品国产av蜜桃| 可以免费在线观看a视频的电影网站 | 极品少妇高潮喷水抽搐| 久久综合国产亚洲精品| 午夜福利网站1000一区二区三区| 韩国精品一区二区三区| 亚洲国产精品一区三区| 久久99一区二区三区| 纯流量卡能插随身wifi吗| 一级毛片电影观看| 你懂的网址亚洲精品在线观看| 亚洲国产欧美日韩在线播放| 国产精品欧美亚洲77777| 色网站视频免费| 9热在线视频观看99| 日本vs欧美在线观看视频| 在线精品无人区一区二区三| 9热在线视频观看99| 色婷婷av一区二区三区视频| 亚洲在久久综合| 777久久人妻少妇嫩草av网站| 亚洲综合色惰| 久久久久精品人妻al黑| 一二三四在线观看免费中文在| 最近最新中文字幕免费大全7| 日韩,欧美,国产一区二区三区| 午夜日韩欧美国产| 日本-黄色视频高清免费观看| 亚洲欧美中文字幕日韩二区| 午夜福利乱码中文字幕| 亚洲久久久国产精品| 国精品久久久久久国模美| 亚洲国产精品一区二区三区在线| 国产日韩欧美视频二区| 亚洲国产精品一区二区三区在线| 99热全是精品| 国产精品秋霞免费鲁丝片| 国产淫语在线视频| 夫妻午夜视频| 国产伦理片在线播放av一区|