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

    錨固體系作用下隧道圍巖力學(xué)響應(yīng)的全過(guò)程解析

    2022-07-04 07:39:22孫振宇張頂立劉道平侯艷娟
    工程力學(xué) 2022年7期
    關(guān)鍵詞:施作塑性錨索

    孫振宇,張頂立,2,劉道平,侯艷娟,李 奧,3

    (1. 北京交通大學(xué)城市地下工程教育部重點(diǎn)實(shí)驗(yàn)室,北京 100044;2. 北京交通大學(xué)川藏鐵路研究院,北京 100044;3. 蘇交科集團(tuán)股份有限公司,江蘇,南京 210019)

    隧道臨空面形成后開(kāi)挖面后方圍巖失去支撐而圍巖變形速率較快,因此需要采用及時(shí)施作的初期支護(hù)以控制圍巖變形。錨固體系作為一種及時(shí)有效的支護(hù)手段,因其靈活性強(qiáng)、成本低廉且能主動(dòng)控制圍巖變形而得到設(shè)計(jì)者的青睞,在隧道與地下工程中得到了廣泛應(yīng)用[1?4]。因此,研究錨固體系作用機(jī)理及其變形控制效果,對(duì)于隧道圍巖穩(wěn)定性控制及錨固體系的科學(xué)設(shè)計(jì)具有重要意義。

    目前對(duì)于錨固體系作用機(jī)理的認(rèn)識(shí)主要包括懸吊理論、組合梁理論、壓縮拱理論以及地層加固理論[5?8]。為了研究錨固體系的支護(hù)效果,國(guó)內(nèi)外學(xué)者采用數(shù)值模擬、現(xiàn)場(chǎng)實(shí)測(cè)、室內(nèi)試驗(yàn)和理論解析等手段進(jìn)行了研究。Freeman[9]通過(guò)對(duì)錨桿受力過(guò)程及錨桿應(yīng)力分布進(jìn)行現(xiàn)場(chǎng)試驗(yàn),提出了全長(zhǎng)粘結(jié)錨桿中性點(diǎn)的概念。郭軍等[10]通過(guò)數(shù)值模擬和現(xiàn)場(chǎng)試驗(yàn)方法分析了黃土隧道中系統(tǒng)錨桿的作用效果,指出邊墻錨桿作用效果顯著而拱部錨桿較弱。Zuo等[11]采用室內(nèi)拉拔試驗(yàn),研究了不同直徑匹配下錨固體的破壞模式,提出了保證錨固效果的桿體與鉆孔直徑匹配建議值。蔡躍等[12]基于拉拔試驗(yàn)結(jié)果提出了錨桿與圍巖相互作用分析模型,揭示了桿體與圍巖接觸面間的耦合作用機(jī)理。與數(shù)值分析和物理模擬相比,理論解析過(guò)程透明且使用方便,在進(jìn)行錨固效果分析時(shí)更為直觀,有助于識(shí)別出關(guān)鍵錨固參數(shù),因此在錨固機(jī)理的分析中得到了廣泛應(yīng)用[13?15]。

    針對(duì)隧道錨固機(jī)理的理論解析,主要包括等效支護(hù)力和等效加固圈兩種處理方式。其中等效支護(hù)力將錨桿作用考慮為作用于隧道洞壁的均布力或圍巖內(nèi)部體積力[16],本質(zhì)上是錨固體系的支護(hù)作用;而等效加固圈則將錨桿作用等效為錨固區(qū)內(nèi)圍巖力學(xué)參數(shù)的提高[17],本質(zhì)上是圍巖層面的加固作用。此外,相關(guān)學(xué)者對(duì)于錨桿與錨索的組合支護(hù)形式也進(jìn)行了研究[18?21],對(duì)于認(rèn)識(shí)錨固體系的作用機(jī)理具有重要意義。但上述研究大多未考慮錨固體系作用過(guò)程中與圍巖相對(duì)位置的變化情況,且對(duì)于錨桿與錨索支護(hù)時(shí)機(jī)相對(duì)滯后性的影響機(jī)制未見(jiàn)報(bào)道。

    事實(shí)上,錨固體系伸入圍巖內(nèi)部,兼有支護(hù)與加固的雙重作用,且其與圍巖相互作用過(guò)程復(fù)雜,現(xiàn)有研究?jī)H考慮單一作用而可操作性較差[22]。此外,由于錨固體系施作的滯后性,且考慮到錨桿、錨索作用范圍及其與圍巖結(jié)構(gòu)的位態(tài)關(guān)系,二者作用功能有所區(qū)別[23]。而目前對(duì)于錨固體系在大跨隧道中的協(xié)同作用全過(guò)程機(jī)理仍不十分明確,缺乏科學(xué)合理的設(shè)計(jì)方法。

    為此,本文基于開(kāi)挖面空間效應(yīng),綜合考慮錨固體系的等效支護(hù)力和加固圈效應(yīng),對(duì)錨桿、錨索先后施作時(shí)圍巖力學(xué)狀態(tài)進(jìn)行分析,根據(jù)錨固構(gòu)件與圍巖結(jié)構(gòu)空間位態(tài)關(guān)系的發(fā)展規(guī)律,建立錨固體系的變形控制機(jī)理模型,揭示了錨固作用下隧道圍巖力學(xué)響應(yīng)的全過(guò)程演化機(jī)制,進(jìn)一步分析錨固體系對(duì)于圍巖變形的控制效果及其影響因素,可為隧道支護(hù)體系協(xié)同設(shè)計(jì)提供參考。

    1 分析模型與基本假設(shè)

    1.1 隧道錨固體系作用分析模型

    對(duì)于一深埋隧道,采用等代圓法將其等效為圓形斷面[24]。假定圍巖為各向同性均勻的理想彈塑性體,無(wú)窮遠(yuǎn)處受靜水應(yīng)力p0作用,隧道內(nèi)邊界作用均勻支護(hù)反力pi,隧道半徑為r0,圍巖屈服后在隧道洞周形成半徑為rp的塑性區(qū),考慮塑性圍巖的剪脹效應(yīng)。考慮到錨固體系滯后于隧道開(kāi)挖面施作,且錨桿支護(hù)施作更為及時(shí),故假定錨桿、錨索施作時(shí)距開(kāi)挖面距離分別為x1、x2,如圖1所示。

    圖1 隧道錨固系統(tǒng)的布置Fig. 1 Arrangement of tunnel anchoring system

    錨固區(qū)和原巖區(qū)圍巖均滿足Hoek-Brown強(qiáng)度準(zhǔn)則如下[25]:

    式中:σ1、σ3分別為第一主應(yīng)力和第三主應(yīng)力,在本章假設(shè)中分別為σθ和σr;m和s分別為巖體材料常數(shù);σc為完整巖石單軸抗壓強(qiáng)度。

    1.2 隧道開(kāi)挖面空間效應(yīng)及其表征方式

    隧道圍巖彈塑性變形的釋放并非瞬時(shí)完成,而是隨著開(kāi)挖面推進(jìn)逐漸發(fā)展,因此在錨桿施作時(shí)隧道圍巖已發(fā)生一定程度的變形,此即隧道開(kāi)挖面空間效應(yīng)[26]。此約束效應(yīng)可由虛擬支護(hù)力表征,其求解思路為[27]:首先根據(jù)圍巖特性曲線獲得分析斷面徑向位移ur0與支護(hù)力pi關(guān)系,并得到無(wú)支護(hù)條件下圍巖理論最大位移量ur0(∞),再代入位移釋放率u?x表達(dá)式,進(jìn)而獲得虛擬支護(hù)力與空間位置的關(guān)系。

    關(guān)于隧道位移釋放率的描述,有眾多學(xué)者進(jìn)行了研究,而Hoek基于Chern等[28]的實(shí)測(cè)數(shù)據(jù)采用最佳擬合方法得到的經(jīng)驗(yàn)公式因其表達(dá)式簡(jiǎn)單且將開(kāi)挖面前后統(tǒng)一描述而得到廣泛應(yīng)用。此外,李二兵等[29]通過(guò)現(xiàn)場(chǎng)原位試驗(yàn)表明,Hoek公式對(duì)于圍巖全過(guò)程變形的擬合結(jié)果與實(shí)測(cè)數(shù)據(jù)較為接近,為最優(yōu)求解方程,因此本文采用該公式來(lái)進(jìn)行求解:

    式中,x 為分析斷面距開(kāi)挖面距離,x=0處為開(kāi)挖面,負(fù)值表示開(kāi)挖面前方,正值表示開(kāi)挖面后方。

    1.3 錨固體系全過(guò)程作用工況

    由于錨桿長(zhǎng)度增加時(shí)其柔性特征越顯著,因此實(shí)際工程中常用的錨桿長(zhǎng)度是有限的,而對(duì)于大斷面隧道圍巖穿越不良地質(zhì)體時(shí)塑性區(qū)范圍往往較大,此時(shí)錨桿僅能作用于塑性區(qū),其錨固效果十分有限。此外,對(duì)于高速鐵路隧道往往具有高可靠性需求,則需施作錨索以提高錨固體系的整體剛度。采用鋼絞線組成的錨索一般剛度較大而延伸率較小,長(zhǎng)錨索可伸入穩(wěn)定巖層從而調(diào)動(dòng)圍巖承載。

    鑒于錨桿伸入塑性區(qū)時(shí)有必要施作錨索以保證錨固體系協(xié)同承載,因此本文主要針對(duì)錨索伸入彈性區(qū)且錨桿伸入塑性區(qū)的最終狀態(tài)進(jìn)行解析。根據(jù)錨桿、錨索施作時(shí)隧道圍巖彈塑性狀態(tài)及相互作用過(guò)程中錨桿、錨索與塑性區(qū)空間位置關(guān)系,可將錨固體系與圍巖相互作用全過(guò)程分為6種工況,如圖2所示。

    圖2 隧道錨固體系與圍巖作用全過(guò)程機(jī)理模型Fig. 2 Mechanism model of whole process between tunnel anchorage system and surrounding rock

    工況A,錨桿施作時(shí)圍巖處于彈性狀態(tài)且錨桿施作伸入彈性區(qū),錨索施作時(shí)圍巖彈性;

    工況B,錨桿施作時(shí)圍巖處于彈性狀態(tài),錨索施作時(shí)圍巖塑性且錨桿伸入塑性區(qū);

    工況C,錨桿施作時(shí)圍巖處于彈性狀態(tài),錨索施作時(shí)圍巖塑性且錨桿伸入彈性區(qū);

    工況D,錨桿施作時(shí)圍巖處于塑性且錨桿施作伸入彈性區(qū),錨索施作時(shí)圍巖塑性且錨桿伸入彈性區(qū);

    工況E,錨桿施作時(shí)圍巖處于塑性且錨桿施作伸入彈性區(qū),錨索施作時(shí)圍巖塑性且錨桿伸入塑性區(qū);

    工況F,圍巖塑性階段錨桿施作伸入塑性區(qū),錨索施作時(shí)圍巖塑性。

    2 錨固體系與圍巖作用全過(guò)程解析

    2.1 工況A解析

    無(wú)支護(hù)條件下的隧道彈性位移可表示為:

    式中,E和μ分別為圍巖彈性模量和泊松比。

    將r=r0,pi=0代入式(3)中可得:

    將x=x1代入式(2)可得錨桿施作時(shí)圍巖已釋放位移量:

    結(jié)合式(3)可得此時(shí)虛擬支護(hù)力為:

    因此,對(duì)于x<x1時(shí)圍巖變形采用式(3)計(jì)算。當(dāng)x1≤x<x2時(shí)僅有錨桿作用,根據(jù)Fahimifar等[16],錨桿的等效加固圈效應(yīng)可考慮為圍巖彈性模量的增加,則錨桿單一錨固區(qū)圍巖彈性模量為:

    若錨桿在圍巖變形發(fā)生前即已施作,則圍巖變形可表示為:

    式中,rsb=Lsb+r0,Lsb為錨桿長(zhǎng)度,prsb為錨桿單一錨固區(qū)與原巖區(qū)徑向接觸力,有[17]:

    考慮到彈性區(qū)虛擬支護(hù)力與圍巖參數(shù)無(wú)關(guān),將式(6)代入式(8)并考慮式(5)可得由于錨桿施作的滯后性而產(chǎn)生的增量位移為:

    則錨桿施作后圍巖位移表達(dá)式為:

    對(duì)于某一分析斷面,虛擬支護(hù)力表達(dá)式為[27]:

    鑒于錨索長(zhǎng)度大于錨桿,因此,將錨桿與錨索共同作用區(qū)域稱(chēng)為雙重錨固區(qū),錨索長(zhǎng)于錨桿部分區(qū)域稱(chēng)為錨索單一錨固區(qū)。設(shè)隧道縱向上彈塑性交界臨界距離為xcr,當(dāng)x1<x2<xcr時(shí)錨索施作時(shí)圍巖彈性,則錨索施作后雙重錨固區(qū)彈性模量可表示為:

    式中:下標(biāo)b和c分別表示錨桿和錨索;Es、As、Sc和Sl分別表示錨桿(索)彈性模量、截面積、間距和排距。

    同理可得錨索單一錨固區(qū)圍巖彈性模量為:

    對(duì)于內(nèi)徑r1、外徑r2的厚壁圓筒,當(dāng)內(nèi)外荷載分別為p1和p2時(shí)彈性徑向變形可表示為[30]:

    考慮到錨索單一錨固區(qū)和雙重錨固區(qū)及錨索單一錨固區(qū)與原巖區(qū)位移連續(xù)條件,則有:

    將式(16)代入式(15)可得錨固區(qū)與原巖區(qū)界面上徑向應(yīng)力prsc和錨索單一錨固區(qū)與雙重錨固區(qū)徑向接觸力prsb表達(dá)式為:

    式中,各系數(shù)表達(dá)式如下:

    根據(jù)式(10)易推知由于錨桿和錨索施作的滯后性而產(chǎn)生的增量位移為:

    當(dāng)x>xcr時(shí)圍巖塑性區(qū)產(chǎn)生,無(wú)支護(hù)條件下隧道塑性區(qū)圍巖變形表達(dá)式為[31]:

    將r=r0代入式(22)可得虛擬支護(hù)力表達(dá)式為:

    當(dāng)虛擬支護(hù)力小于臨界值時(shí)將發(fā)生塑性位移,該臨界值為:

    將式(24)與式(25)聯(lián)立并結(jié)合式(5)可得隧道縱向彈塑性分界處距開(kāi)挖面距離為:

    考慮錨桿支護(hù)的等效支護(hù)力效應(yīng),將此作用力以附加體積力作用于圍巖。當(dāng)塑性區(qū)開(kāi)始產(chǎn)生時(shí)錨桿最初伸入彈性區(qū),此時(shí)錨固區(qū)圍巖應(yīng)力為[17]:

    式中,σsb、σsc分別為錨桿、錨索的徑向約束力,有:

    假定錨索與圍巖變形協(xié)調(diào),則有:

    則有:

    式中:

    將式(31)代入平衡方程:

    并積分可得雙重錨固塑性區(qū)圍巖徑向應(yīng)力為:

    式中:f1(r0)為與邊界條件相關(guān)的待定參數(shù);f1(r)需用數(shù)值方法求解。

    隧道圍巖雙重錨固彈性區(qū)應(yīng)力表達(dá)式為:

    式中,σp為彈塑性界面處徑向應(yīng)力??紤]到彈性區(qū)有σr+σθ=2p0,將其代入式(1)可得:

    由式(17)易推得prsc和prsb表達(dá)式為:

    式中:

    考慮到彈塑性界面及隧道洞壁應(yīng)力邊界條件:

    則聯(lián)立式(33)~式(37)可得隧道塑性區(qū)半徑表達(dá)式,此時(shí)隧道洞壁圍巖變形表達(dá)式為:

    當(dāng)x>x?時(shí),錨桿端部位于塑性區(qū)內(nèi),此時(shí)錨索單一錨固塑性區(qū)徑向應(yīng)力滿足:

    代入式(32)并積分可得錨索單一錨固塑性區(qū)圍巖徑向應(yīng)力表達(dá)式為:

    式中:f2(r0)為與邊界條件相關(guān)的待定參數(shù);f2(r)需用數(shù)值方法求解。

    雙重錨固區(qū)與單一錨固區(qū)應(yīng)力連續(xù)性條件為:

    則聯(lián)立式(33)、式(37)、式(35)、式(40)和式(41)可得隧道圍巖塑性區(qū)半徑,再代入式(38)并結(jié)合式(24)即可得到隧道圍巖變形量沿隧道縱向的演化規(guī)律。

    2.2 工況B解析

    當(dāng)x2>x?>xcr時(shí)錨索施作時(shí)圍巖處于塑性狀態(tài)且錨桿伸入圍巖塑性區(qū),因此圍巖變形的縱向發(fā)展規(guī)律分為五段。

    1)對(duì)于x<x1段圍巖變形采用式(3)計(jì)算,而虛擬支護(hù)力表達(dá)式為:

    2)對(duì)于x1<x<xcr段圍巖變形采用式(11)計(jì)算,虛擬支護(hù)力采用式(42)。

    3)對(duì)于xcr≤x<x?段,根據(jù)前文分析可得錨桿單一錨固區(qū)內(nèi)圍巖徑向應(yīng)力滿足:

    錨桿施作時(shí)圍巖徑向應(yīng)變?yōu)椋?/p>

    由于塑性區(qū)全部屬于錨固區(qū),將式(7)代入式(22)可得錨固塑性區(qū)圍巖徑向應(yīng)變表達(dá)式:

    將式(44)和式(45)代入式(43)可得:

    則錨桿單一錨固區(qū)圍巖徑向應(yīng)力表達(dá)式為:

    式中:f(r0)為待定系數(shù);f(r)需用數(shù)值方法求解。

    此時(shí)錨固區(qū)與原巖區(qū)徑向接觸力prsb為:

    聯(lián)立式(34)、式(47)~式(48)并考慮式(37)即可得塑性區(qū)半徑rp,進(jìn)而得到塑性區(qū)應(yīng)力解答。代入式(22)并結(jié)合式(7)可得錨桿施作后塑性圍巖變形表達(dá)式為:

    式中,虛擬支護(hù)力采用式(24)計(jì)算。

    4)對(duì)于x?≤x<x2段,錨桿端部位于塑性區(qū)內(nèi),將式(1)與式(32)聯(lián)立并代入邊界條件r=rp,σr=σp可得塑性原巖區(qū)應(yīng)力表達(dá)式為:

    將r=rs代入式(50)可得錨固區(qū)與原巖區(qū)界面處徑向應(yīng)力為:

    因此,圍巖變形采用式(49)計(jì)算,塑性區(qū)半徑由式(37)、式(47)和式(51)求得,虛擬支護(hù)力采用式(24)。

    5)對(duì)于x≥x2段圍巖變形采用式(38)計(jì)算,塑性區(qū)半徑由式(33)、式(37)、式(40)~式(41)求得,虛擬支護(hù)力采用式(24)。

    2.3 工況C解析

    當(dāng)xcr<x2<x?時(shí)錨索施作時(shí)圍巖處于塑性狀態(tài)且錨桿伸入圍巖彈性區(qū),因此圍巖變形的縱向發(fā)展規(guī)律亦分為五段:

    1)對(duì)于x<x1段圍巖變形采用式(3)計(jì)算,而虛擬支護(hù)力采用式(42);

    2)對(duì)于x1≤x<xcr段圍巖處于彈性狀態(tài),圍巖變形采用式(11)計(jì)算,虛擬支護(hù)力采用式(42);

    3)對(duì)于xcr≤x<x2段圍巖處于塑性狀態(tài)且錨桿伸入彈性區(qū),圍巖變形采用式(49)計(jì)算,塑性區(qū)半徑由式(37)、式(47)~式(48)聯(lián)立求得,虛擬支護(hù)力采用式(24);

    4)對(duì)于x2≤x<x?段圍巖處于塑性狀態(tài)且錨桿、錨索均伸入彈性區(qū),圍巖變形采用式(38)計(jì)算,塑性區(qū)半徑由式(37)聯(lián)立式(33)~式(36)求得,虛擬支護(hù)力采用式(24);

    5)對(duì)于x≥x?段圍巖變形采用式(38)計(jì)算,塑性區(qū)半徑由式(33)、式(37)和式(40)~式(41)求得,虛擬支護(hù)力采用式(24)。

    2.4 工況D解析

    當(dāng)xcr<x1<x2<x?時(shí)圍巖塑性階段錨桿施作且伸入彈性區(qū),錨索施作時(shí)圍巖處于塑性且錨桿伸入彈性區(qū),因此圍巖變形縱向發(fā)展規(guī)律亦分為五段:

    1)對(duì)于x<xcr段圍巖處于彈性階段,圍巖變形由式(3)計(jì)算,而虛擬支護(hù)力采用式(42);

    2)對(duì)于xcr≤x<x1段圍巖處于塑性階段,圍巖變形由式(22)計(jì)算,塑性區(qū)半徑表達(dá)式為:

    式中,虛擬支護(hù)力由式(24)計(jì)算。

    3)對(duì)于x1≤x<x2段圍巖處于塑性階段且錨桿伸入彈性區(qū),圍巖變形由式(49)計(jì)算,塑性區(qū)半徑由由式(37)、式(47)~式(48)聯(lián)立求得,虛擬支護(hù)力采用式(24);

    4)對(duì)于x2≤x<x?段圍巖處于塑性階段且錨桿、錨索均伸入彈性區(qū),圍巖變形采用式(38)計(jì)算,塑性區(qū)半徑由式(37)聯(lián)立式(33)~式(36)求得,虛擬支護(hù)力采用式(24);

    5)對(duì)于x≥x?段圍巖處于塑性階段且錨桿伸入塑性區(qū),錨索伸入彈性區(qū),圍巖變形采用式(38)計(jì)算,塑性區(qū)由式(33)、式(37)和式(40)~式(41)求得,虛擬支護(hù)力采用式(24)。

    2.5 工況E解析

    當(dāng)xcr<x1<x?<x2時(shí)圍巖變形的縱向發(fā)展規(guī)律亦分為五個(gè)階段,其中前兩個(gè)階段與工況D相同,而對(duì)于x1≤x<x?時(shí)圍巖處于塑性且錨桿伸入彈性區(qū),此時(shí)圍巖變形由式(49)計(jì)算,塑性區(qū)半徑由式(37)、式(47)~式(48)聯(lián)立求得,虛擬支護(hù)力采用式(24)。

    當(dāng)x?≤x<x2時(shí)圍巖處于塑性狀態(tài)且錨桿伸入塑性區(qū),圍巖變形采用式(49)計(jì)算,塑性區(qū)半徑由式(37)、式(47)和式(51)求得,虛擬支護(hù)力采用式(24);

    當(dāng)x≥x2時(shí)錨桿伸入塑性區(qū)且錨索伸入彈性區(qū),此時(shí)計(jì)算方法與工況D最后一階段相同。

    2.6 工況F解析

    當(dāng)xcr<x?<x1<x2時(shí)錨桿在圍巖塑性階段施作且伸入塑性區(qū),圍巖變形的縱向發(fā)展分為四個(gè)階段,其中前兩個(gè)階段與工況D相同,而對(duì)于x1≤x<x2時(shí)圍巖處于塑性狀態(tài)且錨桿伸入塑性區(qū),當(dāng)x1>xcr時(shí)錨桿施作處圍巖已進(jìn)入塑性,圍巖已釋放位移為:

    則錨桿施作時(shí)圍巖徑向應(yīng)變?yōu)椋?/p>

    將式(54)和式(45)代入式(43)可得:

    將式(55)代入式(1)并積分可得塑性區(qū)圍巖徑向應(yīng)力表達(dá)式為:

    式中:f′(r0)為待定系數(shù);f′(r)需用數(shù)值方法求解。

    此時(shí)圍巖變形由式(49)計(jì)算,塑性區(qū)半徑由式(37)、式(51)和式(56)聯(lián)立求得,虛擬支護(hù)力采用式(24)。

    當(dāng)x≥x2時(shí)錨桿伸入塑性區(qū)且錨索伸入彈性區(qū),此時(shí)計(jì)算方法與工況D最后一階段相同。

    3 驗(yàn)證與分析

    3.1 分析工況與數(shù)值模型的建立

    為了驗(yàn)證理論模型的正確性,選取三種典型工況與FLAC3D數(shù)值模擬進(jìn)行對(duì)比分析,其中包含圖2中錨固體系施作時(shí)圍巖力學(xué)狀態(tài)及其與圍巖塑性區(qū)所有可能的位置關(guān)系。針對(duì)半徑為7.5 m的隧道,靜水應(yīng)力場(chǎng)為10 MPa,隧道開(kāi)挖進(jìn)尺取1 m。由于模型的中心對(duì)稱(chēng)性,為節(jié)省計(jì)算時(shí)間僅需對(duì)其1/4進(jìn)行建模分析,如圖3所示。

    圖3 數(shù)值計(jì)算模型及邊界條件Fig. 3 Numerical model and boundary conditions

    三種工況所采用的統(tǒng)一圍巖參數(shù)為E=1 GPa,μ=0.3 ,f=1,錨桿參數(shù)為Fsb=50 kN ,Esb=210 GPa ,Ab=5 cm2,Scb=Slb=1 m,錨索參數(shù)為Fsc=200 kN,Esc=200 GPa,Ac=10 cm2,Scc=Slc=2 m ,其他參數(shù)如表1所示。

    表1 典型工況圍巖和錨固體系參數(shù)Table 1 Parameters of surrounding rock and anchorage system of typical cases

    3.2 計(jì)算結(jié)果分析

    各工況下錨固體系與圍巖相互作用的演化路徑判斷過(guò)程如表2所示。以工況2為例,首先判斷錨桿施作時(shí)圍巖彈塑性狀態(tài),由于xcr=1.58 m<x1,因此錨桿施作時(shí)圍巖處于塑性狀態(tài),進(jìn)而計(jì)算錨桿施作時(shí)圍巖塑性區(qū)半徑可得rp|x=x1=7.84 m,由于rp|x=x1<rsb,因此錨桿施作時(shí)伸入彈性區(qū),然后錨索施作時(shí)圍巖塑性區(qū)半徑進(jìn)行計(jì)算,可得rp|x=x2=10.21 m,由于rp|x=x2<rsb,因此錨索施作時(shí)錨桿仍作用與彈性區(qū),因此可確定演化路徑為路徑D,接著對(duì)錨桿前端恰位于彈塑性界面上的臨界距離進(jìn)行計(jì)算,可得x?=14.56 m,則對(duì)于x>x?即為最終階段。需要指出,對(duì)于工況3,圍巖在開(kāi)挖面前方即已進(jìn)入塑性,且錨桿施作時(shí)塑性區(qū)范圍已大于錨桿長(zhǎng)度,因此為路徑F。

    表2 錨固體系與圍巖相互作用演化路徑的判定Table 2 Determination of evolution path of interaction between anchorage system and surrounding rock

    計(jì)算可得各工況下圍巖縱向變形和塑性區(qū)半徑計(jì)算結(jié)果如圖4所示。計(jì)算表明,基于本文計(jì)算模型得到的錨固體系作用下圍巖縱向變形與數(shù)值模擬結(jié)果基本一致,而數(shù)值模擬的塑性區(qū)半徑受網(wǎng)格影響較大無(wú)法比較,但其粗略值和發(fā)展趨勢(shì)與本文計(jì)算結(jié)果是一致的,說(shuō)明本文推導(dǎo)是正確的、可靠的。此外,由于錨固體系與圍巖參數(shù)的差異,相互作用路徑有所不同,由此造成圍巖全過(guò)程變形演化過(guò)程的差異。遺憾的是,在以往的理論推導(dǎo)中并未考慮錨固體系與圍巖相對(duì)位置差異與錨固時(shí)機(jī)的綜合影響,這顯然與實(shí)際情況是不符的。

    圖4 各工況理論解析與數(shù)值模擬結(jié)果對(duì)比Fig. 4 Comparison of theoretical analysis and numerical simulation results of each working condition

    進(jìn)一步以工況1為例與現(xiàn)有研究中考慮錨固體系等效加固圈或等效支護(hù)力單一作用的方法進(jìn)行對(duì)比,可見(jiàn),僅考慮錨固體系某一方面作用時(shí)圍巖變形和塑性區(qū)半徑計(jì)算值較本文結(jié)果偏大,且當(dāng)圍巖變形量越大時(shí)誤差也隨之增大,這表明以往理論低估了錨固體系的作用,而在大變形圍巖條件中則更為明顯,這顯然是不可接受的。因此,采用本文解進(jìn)行錨固體系的設(shè)計(jì)將更為合理,尤其對(duì)于大變形圍巖條件優(yōu)越性則更為顯著。

    4 隧道錨固體系變形控制效果分析

    4.1 錨固體系對(duì)圍巖全過(guò)程變形的控制原理

    施工影響下隧道圍巖縱向變形可依次劃分為初始變形、急劇變形、變形減緩和變形穩(wěn)定等4個(gè)階段,如圖5所示。其中急劇變形階段所占比例最大,一般可占到總變形量的60%左右,根據(jù)控制重點(diǎn)的不同又可劃分為開(kāi)挖面前方和后方兩個(gè)亞階段,即圖中II-A和II-B階段。II-A階段控制的核心為防止開(kāi)挖面失穩(wěn),為隧道開(kāi)挖創(chuàng)造條件,這一功能通常由超前支護(hù)實(shí)現(xiàn);而II-B段則對(duì)應(yīng)隧道圍巖大變形控制,通過(guò)及時(shí)支護(hù)減小圍巖急劇變形的時(shí)間,從而控制圍巖變形量。錨固體系工作核心則是對(duì)II-B段圍巖變形進(jìn)行控制,減小急劇變形時(shí)間。

    圖5 隧道錨固體系變形控制原理示意圖Fig. 5 Schematic diagram of deformation control principle of tunnel anchorage system

    隧道錨固體系作業(yè)靈活性高且安裝方便,可在開(kāi)挖面形成后及時(shí)發(fā)揮作用。由前文分析可知,錨桿一方面通過(guò)對(duì)圍巖加固從而提高錨固區(qū)圍巖力學(xué)性能,即等效加固圈效應(yīng),另一方面施加徑向約束從而控制圍巖變形,即等效支護(hù)力效應(yīng)。由圖5可見(jiàn),對(duì)于一般圍巖條件,當(dāng)僅有錨桿作用時(shí)圍巖變形已得到控制。而當(dāng)圍巖變形控制較為嚴(yán)格時(shí)則需進(jìn)一步采用剛度較大的錨索,通過(guò)錨固體系的協(xié)同作用實(shí)現(xiàn)對(duì)圍巖變形的主動(dòng)控制。

    本質(zhì)上來(lái)講,隧道圍巖變形的產(chǎn)生是圍巖荷載與支護(hù)抗力博弈過(guò)程中不平衡力作用的結(jié)果,根據(jù)牛頓第二定律,此不平衡力越大則變形加速度越大,進(jìn)而影響到圍巖變形。而支護(hù)抗力與支護(hù)剛度相關(guān),顯然支護(hù)與圍巖系統(tǒng)剛度越大則圍巖變形速率越小,變形控制效果也越好。而錨固體系一方面通過(guò)等效支護(hù)力作用施加支護(hù)剛度,另一方面通過(guò)等效加固圈作用提高圍巖剛度和強(qiáng)度,此雙重作用使得支護(hù)—圍巖系統(tǒng)剛度提高,從而有效控制圍巖變形。此外,由于圍巖變形的時(shí)空效應(yīng),錨固體系的施作時(shí)機(jī)和空間分布對(duì)圍巖變形發(fā)展也有重要影響。

    4.2 影響因素分析

    錨固體系各參數(shù)對(duì)圍巖變形的影響已有相應(yīng)研究,而本節(jié)旨在說(shuō)明錨桿與錨索對(duì)于圍巖變形控制效果的相對(duì)影響程度。由于實(shí)際工程中錨桿與錨索密度的設(shè)計(jì)值有時(shí)相差較大,且對(duì)于錨固體系何時(shí)施作也尚無(wú)統(tǒng)一標(biāo)準(zhǔn)??紤]到錨固體系的等效支護(hù)力與加固程度均與密度有關(guān),而施作時(shí)機(jī)則反映了圍巖位移釋放情況,因此選取這兩個(gè)參數(shù)進(jìn)行分析具有足夠的代表性。

    計(jì)算選取的圍巖參數(shù)與節(jié)4.1中工況1相同,由于錨固體系間距和排距對(duì)圍巖變形影響程度相當(dāng),因此設(shè)錨固體系支護(hù)密度為C=1/(ScSl),參考實(shí)際工程中錨固體系的設(shè)計(jì)方案并結(jié)合規(guī)范[32],確定各參數(shù)取值情況如表3所示。

    表3 錨固體系計(jì)算方案Table 3 Calculation scheme of anchorage system

    采用前文公式計(jì)算可得錨固體系典型參數(shù)對(duì)隧道圍巖全過(guò)程變形曲線的影響如圖6所示??梢?jiàn),隨著錨桿支護(hù)時(shí)機(jī)的延遲和錨桿密度的減小,隧道圍巖變形量顯著增大,x1=5 m時(shí)圍巖最大變形量比x1=1 m時(shí)增大約23%,而Cb=0.25根/m2時(shí)圍巖最大變形量比Cb=4根/m2時(shí)增大約54%,這表明錨桿支護(hù)參數(shù)對(duì)圍巖變形控制效果影響十分明顯。隨著錨索密度的減小,圍巖最大變形量也有一定程度的增大,Cb=0.25根/m2時(shí)圍巖最大變形量比Cb=4根/m2時(shí)增大11%。特別地,錨索密度為Cb=4根/m2與Cb=2根/m2時(shí)圍巖變形相差不大,因此實(shí)際工程中錨索無(wú)需布置過(guò)密,而錨索支護(hù)時(shí)機(jī)的延遲對(duì)圍巖變形量無(wú)明顯影響。

    圖6 隧道錨固體系典型參數(shù)對(duì)圍巖全過(guò)程變形影響分析Fig. 6 Influence of typical parameters of tunnel anchorage system on the whole process of surrounding rock deformation

    因此,錨固體系對(duì)圍巖變形的控制效果主要受錨桿參數(shù)的影響,而錨索對(duì)圍巖變形影響不如錨桿顯著,由此可確定錨索的安全儲(chǔ)備作用,應(yīng)在錨桿—圍巖變形發(fā)展至一定程度后施作,其意義主要在于對(duì)于大變形圍巖情況增大錨固系統(tǒng)的整體剛度,并保證運(yùn)營(yíng)期間錨固體系的安全儲(chǔ)備系數(shù)。

    5 結(jié)論

    本文基于隧道開(kāi)挖面空間效應(yīng),考慮錨固體系支護(hù)時(shí)機(jī)的相對(duì)滯后性及錨固區(qū)與塑性區(qū)邊界的空間位態(tài)關(guān)系,建立了錨固體系與圍巖相互作用的全過(guò)程演化機(jī)制模型,進(jìn)一步分析了錨固體系的變形控制效果及其影響因素,對(duì)錨固體系與圍巖相互作用全過(guò)程給出了全新的認(rèn)識(shí),據(jù)此可實(shí)現(xiàn)錨固體系的定量化設(shè)計(jì)。得到如下主要結(jié)論:

    (1)由于錨桿、錨索支護(hù)時(shí)機(jī)及長(zhǎng)度的差異性,二者施作時(shí)圍巖可能處于不同的應(yīng)力狀態(tài),而相互作用過(guò)程中錨固區(qū)與塑性區(qū)的空間位態(tài)也在不斷變化,據(jù)此將錨固體系與隧道圍巖作用全過(guò)程分為6種發(fā)展路徑并進(jìn)行解析,獲得了不同路徑下圍巖變形與塑性區(qū)半徑的演化規(guī)律。

    (2)通過(guò)選取典型的計(jì)算工況與數(shù)值模擬結(jié)果進(jìn)行對(duì)比,驗(yàn)證了本文綜合考慮等效支護(hù)力和加固圈效應(yīng)進(jìn)行推導(dǎo)的正確性,并與考慮錨固體系單一作用的結(jié)果進(jìn)行對(duì)比分析,指出現(xiàn)有理論低估了錨固體系作用,而對(duì)于大變形圍巖條件是不可接受的。

    (3)分析了錨固體系的變形控制原理,錨固體系作為主動(dòng)支護(hù)形式,一方面通過(guò)提供徑向約束改善圍巖受力狀態(tài),另一方面提高錨固區(qū)圍巖力學(xué)性能,本質(zhì)上是增大支護(hù)—圍巖系統(tǒng)剛度,從而實(shí)現(xiàn)對(duì)圍巖急劇變形量的控制。

    (4)由于錨固體系與圍巖相對(duì)位置的變化,在不同的錨固時(shí)機(jī)和錨固長(zhǎng)度影響下隧道圍巖計(jì)算工況具有較大差異,而此差異在以往理論推導(dǎo)中并未考慮,由此造成圍巖變形和塑性區(qū)計(jì)算結(jié)果與實(shí)際情況有所出入。

    (5)計(jì)算分析表明,錨固體系的變形控制效果主要受錨桿影響,而由于錨桿施作時(shí)圍巖已發(fā)生一定程度的變形,為保證支護(hù)效果應(yīng)盡早施作。對(duì)于一般圍巖條件,錨索延伸率低滯后于錨桿施作,此時(shí)圍巖變形已得到較大釋放,錨索參數(shù)對(duì)圍巖變形影響較弱,由此可確定錨索的安全儲(chǔ)備作用。

    猜你喜歡
    施作塑性錨索
    泥石流沖溝隧道洞口適宜工法及支護(hù)施作時(shí)機(jī)研究*
    基于應(yīng)變梯度的微尺度金屬塑性行為研究
    趙莊煤業(yè)13102巷圍巖補(bǔ)強(qiáng)支護(hù)技術(shù)研究與應(yīng)用
    煤(2021年10期)2021-10-11 07:06:12
    隧道濕噴混凝土回彈率影響因素及施作工藝研究
    硬脆材料的塑性域加工
    鈹材料塑性域加工可行性研究
    擠壓性圍巖大跨隧道預(yù)應(yīng)力錨索作用機(jī)制研究
    煤礦深部高地壓巷道錨索支護(hù)技術(shù)對(duì)策研究
    不同錨索參數(shù)對(duì)預(yù)應(yīng)力場(chǎng)影響的數(shù)值模擬研究
    石英玻璃的熱輔助高效塑性域干磨削
    在线观看美女被高潮喷水网站| 亚洲 国产 在线| 国产伦一二天堂av在线观看| 啦啦啦韩国在线观看视频| 午夜精品一区二区三区免费看| 九九爱精品视频在线观看| 亚洲美女搞黄在线观看 | 免费人成视频x8x8入口观看| 天天躁日日操中文字幕| 亚洲va在线va天堂va国产| 国产真实乱freesex| 韩国av在线不卡| 色视频www国产| 色尼玛亚洲综合影院| 尤物成人国产欧美一区二区三区| 国产色婷婷99| 国产伦一二天堂av在线观看| 国产精品久久久久久久电影| 国产av在哪里看| www.色视频.com| 中文字幕熟女人妻在线| 国模一区二区三区四区视频| 97超级碰碰碰精品色视频在线观看| 婷婷丁香在线五月| 内地一区二区视频在线| 国产高清三级在线| 韩国av一区二区三区四区| 国产精品乱码一区二三区的特点| 久久香蕉精品热| 22中文网久久字幕| 国产伦精品一区二区三区视频9| 精品久久国产蜜桃| 欧美性感艳星| 校园人妻丝袜中文字幕| 嫩草影院入口| 一区二区三区免费毛片| 嫩草影院精品99| 亚洲内射少妇av| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一小说| 能在线免费观看的黄片| 久久午夜亚洲精品久久| 欧美色欧美亚洲另类二区| 老司机午夜福利在线观看视频| 特大巨黑吊av在线直播| 午夜精品久久久久久毛片777| 日韩亚洲欧美综合| 在线看三级毛片| 欧美丝袜亚洲另类 | 国产精品亚洲美女久久久| 99久久九九国产精品国产免费| 国产白丝娇喘喷水9色精品| .国产精品久久| av.在线天堂| 999久久久精品免费观看国产| 国内精品美女久久久久久| 一进一出抽搐gif免费好疼| 色在线成人网| 国产男靠女视频免费网站| 国产精品免费一区二区三区在线| 国产精品一区二区三区四区免费观看 | 久久久久久久精品吃奶| 欧美成人一区二区免费高清观看| 国产成人影院久久av| 亚洲欧美激情综合另类| 真实男女啪啪啪动态图| 免费在线观看日本一区| 欧美激情国产日韩精品一区| 伊人久久精品亚洲午夜| 伦精品一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 99久久精品热视频| 国产亚洲精品综合一区在线观看| 亚洲精品在线观看二区| 精品99又大又爽又粗少妇毛片 | 深夜a级毛片| 成人永久免费在线观看视频| 精品人妻一区二区三区麻豆 | 99久久中文字幕三级久久日本| 九九在线视频观看精品| 看黄色毛片网站| 亚洲美女搞黄在线观看 | 亚洲精品色激情综合| 免费无遮挡裸体视频| 亚洲图色成人| 免费在线观看成人毛片| 综合色av麻豆| 亚洲美女黄片视频| 免费无遮挡裸体视频| 啦啦啦观看免费观看视频高清| 国内精品久久久久久久电影| 丝袜美腿在线中文| 色综合亚洲欧美另类图片| 最新中文字幕久久久久| 午夜福利高清视频| 亚洲性久久影院| 亚洲在线观看片| 午夜爱爱视频在线播放| 91久久精品国产一区二区成人| 亚洲精品色激情综合| 亚洲乱码一区二区免费版| 身体一侧抽搐| 天天一区二区日本电影三级| 日韩大尺度精品在线看网址| 免费看日本二区| 国产亚洲精品久久久com| 欧美成人性av电影在线观看| 国产国拍精品亚洲av在线观看| 97人妻精品一区二区三区麻豆| 久久热精品热| av专区在线播放| 久久久久精品国产欧美久久久| 成年免费大片在线观看| 亚洲精品影视一区二区三区av| 国产精品一区二区三区四区久久| 亚洲乱码一区二区免费版| 亚洲专区国产一区二区| 亚洲三级黄色毛片| 啦啦啦啦在线视频资源| 免费搜索国产男女视频| 久久精品综合一区二区三区| 国产精品,欧美在线| 男女做爰动态图高潮gif福利片| 在线观看av片永久免费下载| 国产一区二区三区av在线 | 精品久久久久久,| 欧美三级亚洲精品| 白带黄色成豆腐渣| 中文字幕久久专区| 观看免费一级毛片| 日本撒尿小便嘘嘘汇集6| 毛片女人毛片| 一边摸一边抽搐一进一小说| 国产高清有码在线观看视频| 天天躁日日操中文字幕| 一边摸一边抽搐一进一小说| 一本一本综合久久| 男人和女人高潮做爰伦理| 91久久精品国产一区二区成人| 亚洲无线在线观看| 国产精品一区二区三区四区久久| 在线免费观看的www视频| 色吧在线观看| 成人国产麻豆网| 国产精品福利在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 欧美日本视频| 99久久精品一区二区三区| 又黄又爽又刺激的免费视频.| 亚洲人成网站在线播| 性欧美人与动物交配| 国产精品一区www在线观看 | 黄色一级大片看看| 色5月婷婷丁香| 色综合婷婷激情| 欧美潮喷喷水| 最近最新中文字幕大全电影3| 三级毛片av免费| 精品一区二区三区视频在线| 97超级碰碰碰精品色视频在线观看| 一级av片app| 99热这里只有精品一区| 国产大屁股一区二区在线视频| 免费观看精品视频网站| 亚洲午夜理论影院| 又紧又爽又黄一区二区| 欧美+亚洲+日韩+国产| 亚洲国产精品sss在线观看| 久久婷婷人人爽人人干人人爱| 成年女人永久免费观看视频| 噜噜噜噜噜久久久久久91| 免费不卡的大黄色大毛片视频在线观看 | 99精品在免费线老司机午夜| av在线老鸭窝| 级片在线观看| 校园春色视频在线观看| 一级av片app| 干丝袜人妻中文字幕| 女人被狂操c到高潮| 99久国产av精品| 国产精品久久电影中文字幕| 国产成人a区在线观看| 热99re8久久精品国产| 亚洲av第一区精品v没综合| 日韩欧美一区二区三区在线观看| 他把我摸到了高潮在线观看| 精品人妻1区二区| 久久人人精品亚洲av| 亚洲国产精品sss在线观看| www.色视频.com| 床上黄色一级片| 国产欧美日韩精品亚洲av| 国产一区二区在线av高清观看| 少妇的逼好多水| 黄片wwwwww| 久久亚洲精品不卡| 精品一区二区三区av网在线观看| 99久久成人亚洲精品观看| 成人美女网站在线观看视频| 日韩欧美在线乱码| .国产精品久久| 三级国产精品欧美在线观看| 日韩一本色道免费dvd| 精品国产三级普通话版| 免费电影在线观看免费观看| 热99在线观看视频| 日日撸夜夜添| 亚洲无线在线观看| 少妇裸体淫交视频免费看高清| 日本 欧美在线| 国模一区二区三区四区视频| 欧美成人一区二区免费高清观看| 一区二区三区免费毛片| 国产精品嫩草影院av在线观看 | 成年免费大片在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲真实伦在线观看| 18禁黄网站禁片免费观看直播| 久久亚洲精品不卡| 亚洲av免费在线观看| 成年女人永久免费观看视频| 最近中文字幕高清免费大全6 | 一卡2卡三卡四卡精品乱码亚洲| 内地一区二区视频在线| 成人美女网站在线观看视频| 成人国产综合亚洲| 一进一出抽搐gif免费好疼| 国产午夜精品久久久久久一区二区三区 | 久久久精品欧美日韩精品| 午夜亚洲福利在线播放| 国产精品女同一区二区软件 | 韩国av一区二区三区四区| 日本-黄色视频高清免费观看| 人人妻,人人澡人人爽秒播| 久久久成人免费电影| 国产黄a三级三级三级人| 搡老妇女老女人老熟妇| 欧美zozozo另类| 麻豆精品久久久久久蜜桃| 美女 人体艺术 gogo| 成人二区视频| 亚洲在线观看片| 午夜日韩欧美国产| 国产精品久久久久久久电影| 午夜福利欧美成人| 成年女人永久免费观看视频| 男女视频在线观看网站免费| 97超级碰碰碰精品色视频在线观看| 亚洲av成人精品一区久久| 日韩,欧美,国产一区二区三区 | 一区二区三区激情视频| 亚洲欧美日韩无卡精品| 精品久久久久久久久av| 亚洲aⅴ乱码一区二区在线播放| 久久国内精品自在自线图片| 亚洲一区高清亚洲精品| 真实男女啪啪啪动态图| 日韩国内少妇激情av| 亚洲电影在线观看av| 国产真实伦视频高清在线观看 | 久久国产精品人妻蜜桃| 少妇的逼水好多| 日韩精品中文字幕看吧| 22中文网久久字幕| 91午夜精品亚洲一区二区三区 | 三级毛片av免费| 国产v大片淫在线免费观看| 十八禁国产超污无遮挡网站| 又黄又爽又免费观看的视频| 日本-黄色视频高清免费观看| 91久久精品电影网| 国产一区二区三区av在线 | 午夜福利成人在线免费观看| 亚洲图色成人| 国产精品人妻久久久久久| 麻豆成人av在线观看| 亚洲精品亚洲一区二区| 色噜噜av男人的天堂激情| 亚洲精品粉嫩美女一区| 在线观看美女被高潮喷水网站| 九九热线精品视视频播放| 一级毛片久久久久久久久女| 亚洲七黄色美女视频| 在线观看午夜福利视频| 成年女人毛片免费观看观看9| 国产精品人妻久久久久久| 国产黄片美女视频| 最好的美女福利视频网| 国产精品日韩av在线免费观看| 一边摸一边抽搐一进一小说| 特大巨黑吊av在线直播| 免费在线观看日本一区| 一个人看视频在线观看www免费| 国内少妇人妻偷人精品xxx网站| 婷婷六月久久综合丁香| 18禁在线播放成人免费| 精品一区二区三区视频在线| 一进一出抽搐gif免费好疼| 国产男靠女视频免费网站| 国产午夜精品久久久久久一区二区三区 | 国产欧美日韩一区二区精品| 国产精品久久电影中文字幕| 亚洲最大成人中文| 婷婷精品国产亚洲av在线| 国产男靠女视频免费网站| 亚洲经典国产精华液单| 中文字幕av成人在线电影| 我的女老师完整版在线观看| 超碰av人人做人人爽久久| 国产真实乱freesex| 啦啦啦啦在线视频资源| 欧美潮喷喷水| 久久热精品热| 精品久久国产蜜桃| 三级毛片av免费| 观看免费一级毛片| 在线观看66精品国产| 国内揄拍国产精品人妻在线| 欧美成人性av电影在线观看| 亚洲一区高清亚洲精品| 国产免费男女视频| 国产不卡一卡二| 久久久精品欧美日韩精品| 亚洲经典国产精华液单| 99热这里只有是精品50| 小蜜桃在线观看免费完整版高清| 国产老妇女一区| 精品久久久久久久久亚洲 | 韩国av在线不卡| 高清在线国产一区| 欧美区成人在线视频| 他把我摸到了高潮在线观看| 最近最新免费中文字幕在线| 久久精品国产鲁丝片午夜精品 | 国产精品日韩av在线免费观看| av女优亚洲男人天堂| 99热这里只有是精品50| 在线天堂最新版资源| 国产主播在线观看一区二区| or卡值多少钱| 嫩草影院入口| 婷婷丁香在线五月| 日韩欧美在线乱码| av.在线天堂| 精品久久久久久久久亚洲 | 淫妇啪啪啪对白视频| .国产精品久久| avwww免费| 一夜夜www| 成人鲁丝片一二三区免费| 国产大屁股一区二区在线视频| 亚洲成av人片在线播放无| 亚洲专区国产一区二区| www日本黄色视频网| 黄色视频,在线免费观看| 国产精品久久久久久久久免| 免费黄网站久久成人精品| 亚洲三级黄色毛片| 变态另类成人亚洲欧美熟女| www日本黄色视频网| 免费观看人在逋| 成熟少妇高潮喷水视频| 一个人看的www免费观看视频| 一a级毛片在线观看| 国产精品嫩草影院av在线观看 | 国产精品野战在线观看| 国产三级中文精品| 最近视频中文字幕2019在线8| 亚洲一区高清亚洲精品| 少妇熟女aⅴ在线视频| eeuss影院久久| 久久午夜福利片| 国产色爽女视频免费观看| 九九在线视频观看精品| 国产成人福利小说| 真人一进一出gif抽搐免费| 真实男女啪啪啪动态图| 日本黄色视频三级网站网址| 香蕉av资源在线| a在线观看视频网站| 99久久成人亚洲精品观看| 啦啦啦观看免费观看视频高清| 免费av不卡在线播放| 亚洲av.av天堂| 高清毛片免费观看视频网站| 久久久精品大字幕| 精品人妻偷拍中文字幕| 黄色视频,在线免费观看| 亚洲国产精品久久男人天堂| 亚洲av五月六月丁香网| 精品久久久久久,| 观看美女的网站| 免费黄网站久久成人精品| 老熟妇乱子伦视频在线观看| 成人三级黄色视频| 男女之事视频高清在线观看| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av| 日本精品一区二区三区蜜桃| 久99久视频精品免费| 免费av观看视频| 亚洲av中文av极速乱 | 日本三级黄在线观看| 99视频精品全部免费 在线| 国产精品一区二区三区四区久久| 99在线视频只有这里精品首页| 波多野结衣巨乳人妻| 国内精品美女久久久久久| 深夜精品福利| 在线免费十八禁| 午夜视频国产福利| 中文字幕av成人在线电影| a级毛片a级免费在线| 99久国产av精品| 欧美精品啪啪一区二区三区| 精品一区二区免费观看| 岛国在线免费视频观看| 校园春色视频在线观看| 色哟哟哟哟哟哟| 校园人妻丝袜中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 成人高潮视频无遮挡免费网站| 最新中文字幕久久久久| 日本欧美国产在线视频| 免费观看精品视频网站| 人人妻人人澡欧美一区二区| 亚洲av第一区精品v没综合| 国产在线精品亚洲第一网站| 嫩草影视91久久| 亚洲五月天丁香| 久久久精品大字幕| 午夜福利18| 免费人成视频x8x8入口观看| 人妻丰满熟妇av一区二区三区| 无遮挡黄片免费观看| 校园人妻丝袜中文字幕| 成人美女网站在线观看视频| 日韩精品青青久久久久久| 国产视频内射| 制服丝袜大香蕉在线| 午夜精品久久久久久毛片777| 级片在线观看| 久久精品国产鲁丝片午夜精品 | 成人亚洲精品av一区二区| 日韩在线高清观看一区二区三区 | 亚洲精品一区av在线观看| 国产成年人精品一区二区| 国内精品久久久久久久电影| 国产精品不卡视频一区二区| www.色视频.com| 99riav亚洲国产免费| 亚洲久久久久久中文字幕| 一区二区三区高清视频在线| 亚洲一级一片aⅴ在线观看| 久久99热6这里只有精品| 亚洲欧美日韩东京热| 国产国拍精品亚洲av在线观看| 国产成人aa在线观看| 啦啦啦韩国在线观看视频| 国产一区二区三区av在线 | 欧美高清成人免费视频www| 丝袜美腿在线中文| 午夜爱爱视频在线播放| 日本撒尿小便嘘嘘汇集6| 国产白丝娇喘喷水9色精品| 国产欧美日韩精品亚洲av| 久久精品国产鲁丝片午夜精品 | 精品午夜福利在线看| 蜜桃亚洲精品一区二区三区| 中文字幕av成人在线电影| 欧美色视频一区免费| 91在线观看av| 校园人妻丝袜中文字幕| 在线观看舔阴道视频| 国产一区二区激情短视频| 中文字幕久久专区| 久久久久久久久久黄片| 最近最新中文字幕大全电影3| 日韩中字成人| 亚洲欧美激情综合另类| 深夜a级毛片| 亚洲人成伊人成综合网2020| eeuss影院久久| 18禁黄网站禁片午夜丰满| 日本与韩国留学比较| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 精品99又大又爽又粗少妇毛片 | 日日干狠狠操夜夜爽| 免费无遮挡裸体视频| 淫秽高清视频在线观看| 亚洲国产精品合色在线| 国产av不卡久久| 啦啦啦韩国在线观看视频| 国产精品精品国产色婷婷| av黄色大香蕉| www.色视频.com| 九九爱精品视频在线观看| 午夜福利欧美成人| 亚洲天堂国产精品一区在线| 精品人妻视频免费看| 欧美高清性xxxxhd video| 欧美bdsm另类| 一卡2卡三卡四卡精品乱码亚洲| 无遮挡黄片免费观看| 毛片女人毛片| 日本黄色视频三级网站网址| 国产蜜桃级精品一区二区三区| 人妻久久中文字幕网| 午夜免费成人在线视频| 最新在线观看一区二区三区| 亚洲专区中文字幕在线| 一级黄片播放器| 亚洲 国产 在线| 可以在线观看的亚洲视频| 国产 一区精品| 黄色丝袜av网址大全| 亚洲精品在线观看二区| 啦啦啦啦在线视频资源| 国产白丝娇喘喷水9色精品| 国产精品嫩草影院av在线观看 | 日本-黄色视频高清免费观看| 欧美xxxx性猛交bbbb| 99久久久亚洲精品蜜臀av| 国产真实乱freesex| 国产精品一区二区免费欧美| 亚洲精品粉嫩美女一区| 欧美中文日本在线观看视频| 日本色播在线视频| 老司机福利观看| 日韩中字成人| 免费人成在线观看视频色| 在线观看av片永久免费下载| 无人区码免费观看不卡| 欧美bdsm另类| 日日干狠狠操夜夜爽| 精品一区二区三区视频在线观看免费| 在线免费十八禁| 国产高清视频在线播放一区| 亚洲av一区综合| 亚洲性久久影院| а√天堂www在线а√下载| 在线a可以看的网站| 精品久久久久久久末码| 在线观看66精品国产| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区性色av| 色视频www国产| 麻豆久久精品国产亚洲av| 在现免费观看毛片| 亚洲成人久久性| 免费看美女性在线毛片视频| 中文亚洲av片在线观看爽| 嫁个100分男人电影在线观看| 99视频精品全部免费 在线| 欧美成人一区二区免费高清观看| 色综合色国产| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av涩爱 | 两性午夜刺激爽爽歪歪视频在线观看| 春色校园在线视频观看| 成人永久免费在线观看视频| 国产伦人伦偷精品视频| 亚洲国产高清在线一区二区三| 国产精品爽爽va在线观看网站| 国产伦一二天堂av在线观看| 美女cb高潮喷水在线观看| 成人高潮视频无遮挡免费网站| 好男人在线观看高清免费视频| av.在线天堂| 99精品久久久久人妻精品| 色噜噜av男人的天堂激情| 俺也久久电影网| 日本免费一区二区三区高清不卡| 午夜激情欧美在线| 国产精品综合久久久久久久免费| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 久久久久国内视频| 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 国产真实伦视频高清在线观看 | 麻豆成人av在线观看| 国产成人aa在线观看| 日本五十路高清| 国产精品一区www在线观看 | 国产午夜精品久久久久久一区二区三区 | 欧美成人一区二区免费高清观看| 国产视频一区二区在线看| 国产精品久久久久久av不卡| 91午夜精品亚洲一区二区三区 | 特级一级黄色大片| 亚洲七黄色美女视频| 精品国内亚洲2022精品成人| 久久久久久大精品| av在线老鸭窝| a级毛片a级免费在线| 不卡视频在线观看欧美| 少妇裸体淫交视频免费看高清| 日韩一区二区视频免费看| 亚洲av.av天堂| 亚洲av电影不卡..在线观看| 日本五十路高清| 欧美中文日本在线观看视频| 精品人妻视频免费看| 天堂网av新在线| 少妇丰满av| 大又大粗又爽又黄少妇毛片口| 一级毛片久久久久久久久女| 久久人妻av系列| 欧美性猛交╳xxx乱大交人| 99视频精品全部免费 在线| 日韩强制内射视频| 国产精品1区2区在线观看.| x7x7x7水蜜桃|