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

    基于D-P準(zhǔn)則的盾構(gòu)隧道圍巖與襯砌結(jié)構(gòu)相互作用分析1)

    2017-03-20 11:31:55川齊春封肖明清
    力學(xué)學(xué)報(bào) 2017年1期
    關(guān)鍵詞:管片滲流塑性

    何 川齊 春封 坤,3)肖明清

    ?(西南交通大學(xué)交通隧道工程教育部重點(diǎn)實(shí)驗(yàn)室,成都610031)

    ?(中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司,武漢430063)

    基于D-P準(zhǔn)則的盾構(gòu)隧道圍巖與襯砌結(jié)構(gòu)相互作用分析1)

    何 川?,2)齊 春?封 坤?,3)肖明清?,4)

    ?(西南交通大學(xué)交通隧道工程教育部重點(diǎn)實(shí)驗(yàn)室,成都610031)

    ?(中鐵第四勘察設(shè)計(jì)院集團(tuán)有限公司,武漢430063)

    在進(jìn)行盾構(gòu)隧道管片襯砌結(jié)構(gòu)載荷計(jì)算時(shí),常采用全土柱或壓力拱理論計(jì)算圍巖松動(dòng)壓力,但當(dāng)盾構(gòu)隧道面臨深埋條件且需計(jì)入形變壓力時(shí),該方法難以適用.鑒于此,基于Drucker-Prager屈服準(zhǔn)則,推導(dǎo)了考慮滲流效應(yīng)影響下圍巖與襯砌結(jié)構(gòu)相互作用的彈塑性解析解,給出了圍巖彈、塑性區(qū)應(yīng)力與位移、塑性區(qū)半徑等關(guān)鍵參數(shù)與支護(hù)阻力間關(guān)系的解析式.闡述了上述解析結(jié)果在確定襯砌結(jié)構(gòu)載荷中的應(yīng)用,即建立圍巖與襯砌結(jié)構(gòu)靜力平衡狀態(tài)并求二者對(duì)應(yīng)曲線的交點(diǎn).進(jìn)一步地,考慮接頭引起管片襯砌結(jié)構(gòu)整體剛度降低對(duì)圍巖與襯砌結(jié)構(gòu)相互作用的影響,引入剛度折減系數(shù),并在襯砌結(jié)構(gòu)圍巖壓力確定中對(duì)施工期流固耦合效應(yīng)的影響和滲流力對(duì)襯砌結(jié)構(gòu)支護(hù)特性曲線的影響進(jìn)行了簡(jiǎn)化處理.最后,通過(guò)算例將解析解與水下盾構(gòu)隧道載荷實(shí)測(cè)值和數(shù)值計(jì)算值進(jìn)行了比較.結(jié)果表明:用解析方法得到的施工期和穩(wěn)定期的管片襯砌結(jié)構(gòu)圍巖壓力比現(xiàn)場(chǎng)實(shí)測(cè)值分別大28%和12%,穩(wěn)定期圍巖壓力比數(shù)值計(jì)算值大5%,可為類(lèi)似工程的設(shè)計(jì)施工提供一定的參考價(jià)值.

    盾構(gòu)隧道,圍巖,襯砌結(jié)構(gòu),相互作用,圍巖壓力,D-P準(zhǔn)則,收斂約束法

    引言

    管片襯砌作為盾構(gòu)隧道的永久支護(hù)結(jié)構(gòu),在盾構(gòu)機(jī)進(jìn)行洞室開(kāi)挖后通常作為主要承載單元承載.在進(jìn)行盾構(gòu)隧道管片襯砌結(jié)構(gòu)載荷分析時(shí),常采用全土柱法或壓力拱法計(jì)算圍巖松動(dòng)壓力[1].但當(dāng)盾構(gòu)隧道面臨深埋條件且需計(jì)入形變壓力時(shí),該方法難以適用,尤其是在當(dāng)前盾構(gòu)工法廣泛應(yīng)用于我國(guó)城市“排水深隧”、深埋交通隧道、深埋油氣管道及煤礦深長(zhǎng)井巷建設(shè)[2-8]的前提下,深埋盾構(gòu)隧道的載荷理論亟待深入探討.

    深埋條件下盾構(gòu)隧道的載荷,也即圍巖與襯砌結(jié)構(gòu)的相互作用,取決于施工開(kāi)挖對(duì)圍巖的擾動(dòng)及其引起圍巖中原始應(yīng)力場(chǎng)的調(diào)整和重新分布,最終形成次生應(yīng)力場(chǎng)的力學(xué)狀態(tài).然而,由于施工方式的特殊性,其力學(xué)過(guò)程極為復(fù)雜,施工期流固耦合效應(yīng)對(duì)圍巖應(yīng)力場(chǎng)的影響不可忽略[9-11].當(dāng)盾構(gòu)機(jī)通過(guò)時(shí),刀盤(pán)掘削、開(kāi)挖排土,擾動(dòng)地層初始應(yīng)力場(chǎng),在此過(guò)程中施工引起圍巖滲流場(chǎng)的變化,進(jìn)一步影響圍巖應(yīng)力的變化;當(dāng)盾構(gòu)機(jī)通過(guò)后,管片襯砌脫環(huán)后支護(hù)圍巖,圍巖與管片襯砌結(jié)構(gòu)相互作用、協(xié)調(diào)變形,形成圍巖次生應(yīng)力場(chǎng).在此過(guò)程中圍巖滲流場(chǎng)漸趨穩(wěn)定,圍巖的力學(xué)狀態(tài)隨之發(fā)生變化.可見(jiàn),盾構(gòu)隧道施工過(guò)程的特殊性決定了管片襯砌脫環(huán)后圍巖與管片襯砌相互作用力學(xué)過(guò)程的復(fù)雜性,并決定了最終作用在管片襯砌結(jié)構(gòu)上的圍巖壓力.而要探求盾構(gòu)隧道圍巖與管片襯砌相互作用關(guān)系,必須合理考慮施工過(guò)程的特殊性,對(duì)其力學(xué)過(guò)程加以分析,主要包括:(1)開(kāi)挖過(guò)程中圍巖滲流場(chǎng)變化對(duì)應(yīng)力場(chǎng)的影響;(2)施工過(guò)程中開(kāi)挖與襯砌支護(hù)對(duì)于圍巖應(yīng)力場(chǎng)的影響.

    目前,對(duì)于第1個(gè)問(wèn)題的求解,蔡勇平等[12]指出了用Mohr-Coulomb(M-C)破壞條件表達(dá)式計(jì)算水工壓力隧洞彈塑性應(yīng)力的載荷條件,并推導(dǎo)了含水圍巖與襯砌應(yīng)力的計(jì)算式.李宗利等[13]以承受內(nèi)水壓的深埋圓形隧洞為對(duì)象,通過(guò)考慮和不考慮滲流場(chǎng)影響兩種方法對(duì)隧洞開(kāi)展了彈塑性解析,結(jié)果表明隨著洞內(nèi)外水頭差的逐漸增大,滲流場(chǎng)對(duì)應(yīng)力場(chǎng)的影響將顯著增大.張常光[14]基于平面應(yīng)變狀態(tài)下統(tǒng)一強(qiáng)度理論和彈脆塑性軟化模型,推導(dǎo)了水工隧洞施工期和運(yùn)行期不同載荷條件下圍巖和襯砌應(yīng)力、位移的統(tǒng)一解.從上述研究來(lái)看,圍巖與襯砌結(jié)構(gòu)間的相互作用受滲流場(chǎng)變化的影響有明顯差異,可見(jiàn)施工方式、施工過(guò)程導(dǎo)致圍巖滲流場(chǎng)變化對(duì)應(yīng)力場(chǎng)的影響不可忽略.

    在第2個(gè)問(wèn)題的求解中,M-C和Hoek-Brown(HB)準(zhǔn)則被廣泛采用.如吳順川等[15]對(duì)圍巖和襯砌分別采用M-C屈服準(zhǔn)則和雙線性本構(gòu)模型,推導(dǎo)了廣義載荷作用下深埋圓形巷道中圍巖和襯砌相互作用的彈塑性解析解.孫闖等[16]基于M-C應(yīng)變軟化模型采用數(shù)值模擬手段將收斂約束法應(yīng)用于高應(yīng)力軟巖巷道圍巖與支護(hù)結(jié)構(gòu)的相互作用.鄒金鋒等[17]基于M-C屈服準(zhǔn)則和應(yīng)力應(yīng)變軟化模型并考慮軸向應(yīng)力和滲透力的共同作用,采用微小徑向應(yīng)力增量逐步求解法推導(dǎo)出軟化圍巖應(yīng)力應(yīng)變的解,結(jié)果表明滲流力對(duì)隧道圍巖應(yīng)力應(yīng)變分布及塑性半徑和圍巖位移的影響不可忽略.Carranza-Torres等[18-19]基于H-B準(zhǔn)則對(duì)收斂約束法(convergence-confinemen method,CCM)的三大基本要素——圍巖縱向變形曲線(longitudinal deformation profile LDP)、圍巖收斂曲線(ground reaction curve,GRC)和支護(hù)特性曲線(support characteristic curve,SCC)進(jìn)行理論解析,并應(yīng)用該方法進(jìn)行了某圓形斷面隧道的支護(hù)結(jié)構(gòu)設(shè)計(jì).蘇永華等[20]基于H-B失效準(zhǔn)則,根據(jù)開(kāi)挖空間圍巖變形壓力與收斂的關(guān)系導(dǎo)出了圓形隧道圍巖載荷計(jì)算公式.

    然而,M-C準(zhǔn)則和H-B準(zhǔn)則均忽略了中間主應(yīng)力的影響,對(duì)圍巖屈服或破壞的解釋與實(shí)際情況難免存在較大偏差.Drucker-Prager(D-P)準(zhǔn)則引入了中間主應(yīng)力的影響,使得對(duì)此類(lèi)問(wèn)題的求解結(jié)果更符合實(shí)際.侯公羽等[21]基于 Levy-Mises本構(gòu)關(guān)系及D-P屈服準(zhǔn)則對(duì)軸對(duì)稱(chēng)圓巷進(jìn)行了理想彈塑性條件的求解.張小波等[22]考慮不同程度中間主應(yīng)力的影響,采用D-P屈服準(zhǔn)則對(duì)圓形巷道圍巖進(jìn)行彈塑性區(qū)分析,并與M-C準(zhǔn)則解、統(tǒng)一強(qiáng)度準(zhǔn)則解和FLAC3D軟件數(shù)值解的結(jié)果進(jìn)行了對(duì)比,認(rèn)為適當(dāng)?shù)膽?yīng)用D-P屈服準(zhǔn)則更能保證工程實(shí)踐的安全性.

    鑒于此,本文基于Drucker-Prager準(zhǔn)則,推導(dǎo)了盾構(gòu)隧道圍巖的滲流力及滲流影響下圍巖與襯砌結(jié)構(gòu)相互作用的彈塑性解析解,包括圍巖彈、塑性區(qū)應(yīng)力與位移、塑性區(qū)半徑等.在此基礎(chǔ)上,以收斂約束法為載體對(duì)大埋深水下盾構(gòu)隧道施工期和穩(wěn)定期的載荷開(kāi)展了分析,并與現(xiàn)場(chǎng)實(shí)測(cè)和數(shù)值計(jì)算結(jié)果進(jìn)行了對(duì)比驗(yàn)證.

    1 圍巖和襯砌結(jié)構(gòu)的滲透水壓力

    運(yùn)用滲流場(chǎng)基本理論推導(dǎo)隧道從襯砌拼裝完成至水壓穩(wěn)定期間圍巖和襯砌結(jié)構(gòu)的滲透水壓力,作為等效考慮滲流效應(yīng)影響下圍巖與襯砌結(jié)構(gòu)相互作用分析的邊界條件.

    對(duì)于深埋隧道,為簡(jiǎn)化計(jì)算,只考慮開(kāi)挖面為圓形的情況.假設(shè)圍巖在半徑r2的范圍內(nèi)形成的穩(wěn)定滲流場(chǎng),水頭為h2(相應(yīng)水壓為pw0),襯砌的內(nèi)徑為r0,外徑為r1,圍巖和襯砌的滲透系數(shù)分別為kr和kc.分析模型如圖1所示.由于隧道的長(zhǎng)度遠(yuǎn)大于斷面尺寸,可以認(rèn)為是無(wú)限長(zhǎng)的圓柱,兩端受約束,按平面問(wèn)題求解.

    圖1 隧道--圍巖滲流分析模型Fig.1 Model of tunnel-surrounding rock seepage analysis

    襯砌內(nèi)任一半徑r處的滲流量為

    式(1)的邊界條件為

    聯(lián)立式(1)和式(2)可解得

    同理可解得圍巖內(nèi)任一半徑r處的滲流量為

    由于隧道襯砌和圍巖處于穩(wěn)定的滲流場(chǎng)內(nèi),流過(guò)每一斷面的滲流量都相等,即Qc=Qr,聯(lián)立式(3)和式(4)可解得襯砌、圍巖任一半徑處的滲透水壓力

    上式為隧道開(kāi)挖穩(wěn)定后的滲流狀態(tài).對(duì)于采用盾構(gòu)法修建的水下隧道而言,盾構(gòu)施工過(guò)程中流固耦合效應(yīng)引起的水壓滯后導(dǎo)致施工期管片襯砌所受水壓小于穩(wěn)定期[10-11,23-25].該降低程度受地層條件和盾構(gòu)施工參數(shù)等多種因素影響,但具體降低多少,目前的研究成果均未量化.為簡(jiǎn)化起見(jiàn),本文采用將式(5)乘以折減系數(shù)的方法作近似處理,折減系數(shù)取值見(jiàn)后文第5節(jié).

    2 圍巖應(yīng)力與位移分析

    圍巖掘進(jìn)后,由于邊界條件的改變和水壓作用,圍巖應(yīng)力狀態(tài)發(fā)生了變化,如果該應(yīng)力狀態(tài)大于圍巖的抗壓強(qiáng)度,會(huì)在一定范圍內(nèi)形成塑性區(qū).當(dāng)襯砌施作后,襯砌的支護(hù)阻力會(huì)對(duì)圍巖的塑性區(qū)范圍產(chǎn)生影響.

    為方便計(jì)算,作如下假定:(1)巖體為連續(xù)、均質(zhì)、各向同性的理想彈塑性材料;(2)忽略圍巖自重對(duì)塑性區(qū)的影響;(3)圍巖側(cè)壓系數(shù)為1,隧道各向承受等壓;(4)隧道斷面為圓形,縱向水平延伸且長(zhǎng)度無(wú)限大.基于以上假設(shè),問(wèn)題簡(jiǎn)化為幾何結(jié)構(gòu)和載荷都軸對(duì)稱(chēng)的平面應(yīng)變厚壁圓筒問(wèn)題,其力學(xué)分析模型如圖2所示.

    圖2 隧道--圍巖彈塑性力學(xué)分析模型Fig.2 Mechanical model of tunnel-surrounding rock elasto-plastic analysis

    2.1 圍巖塑性區(qū)應(yīng)力

    在半徑r處取一個(gè)微小的單元體,由于為軸對(duì)稱(chēng)問(wèn)題,在單元體上只作用徑向和環(huán)向正應(yīng)力而無(wú)剪應(yīng)力,其力學(xué)模型如圖2所示,其中滲透水壓力以體積力的形式作用在單元體上.利用平衡條件可建立含水圍巖單元體的平衡方程,即

    式中,σr和σθ分別為徑向和環(huán)向有效應(yīng)力,拉應(yīng)力為正,壓應(yīng)力為負(fù).后文中應(yīng)力若無(wú)特別說(shuō)明均為有效應(yīng)力,正負(fù)號(hào)規(guī)定相同.β為滲透水壓力的面積作用系數(shù),與材料的孔隙率有關(guān).對(duì)于混凝土,β一般取2/3~1;對(duì)于接近破壞的巖石,β值接近于1[23].

    當(dāng)圍巖在支承壓力作用下達(dá)到屈服狀態(tài)時(shí),應(yīng)滿足Drucker-Prager準(zhǔn)則,即

    I為第一應(yīng)力不變量,J2為第二偏應(yīng)力不變量;α和k是Drucker-Prager準(zhǔn)則系數(shù).對(duì)于地下工程平面應(yīng)變問(wèn)題,可采用三維應(yīng)力狀態(tài)下的壓縮錐擬合條件確定[21,26],即

    式中,c和φ為圍巖的黏聚力和內(nèi)摩擦角.

    圍巖中間主應(yīng)力方向?yàn)樗淼揽v向,根據(jù)平面應(yīng)變條件,得到在含水條件下圍巖的各主應(yīng)力

    式中 μ為圍巖泊松比,因假設(shè)圍巖應(yīng)力為各向等壓,所以有μ=0.5.將式(10)代入式(8)中,得到

    將式(11)代入式(7),可以求得

    將式(5)第2式和式(12)代入式(6),得到關(guān)于σr的一階線性非齊次微分方程

    代入邊界條件

    2.2 圍巖彈性區(qū)應(yīng)力

    設(shè)圍巖彈性區(qū)和塑性區(qū)交界處的半徑為R0,界面上作用有徑向應(yīng)力σR0;圍巖和襯砌的交界的半徑為r1,在r1處作用有支護(hù)阻力pa.將圍巖彈性區(qū)視為一厚壁圓筒,其無(wú)窮遠(yuǎn)處作用有初始地應(yīng)力q0,內(nèi)半徑R0處承受徑向應(yīng)力σR0,滲透水壓pw作用在整個(gè)區(qū)域內(nèi),則由Kirsch公式可得圍巖彈性區(qū)的總應(yīng)力表達(dá)式為

    2.3 圍巖塑性區(qū)半徑與支護(hù)阻力的關(guān)系

    圍巖在彈、塑性區(qū)交界面處(r=R0)的總應(yīng)力連續(xù),即

    上式左右兩側(cè)分別相加,將式(12)和式(16)代入并化簡(jiǎn),可得

    將式(15)代入式(18)并注意到r=R0,可求得洞周徑向支護(hù)阻力關(guān)于圍巖塑性區(qū)半徑的關(guān)系式

    式(19)表達(dá)了在圍巖巖性特征參數(shù)已知時(shí),徑向支護(hù)阻力pa與塑性區(qū)大小R0之間的關(guān)系.

    取q0=10MPa,pw0=1MPa,c=0.5MPa,φ分別為10°,20°,30°,40°,則pa與R0/r1的關(guān)系如圖3(a)所示.可以看出,在同樣R0/r1的條件下,所需支護(hù)阻力隨著φ的減小而增大;反之,當(dāng)pa一定時(shí),隨著φ的減小,塑性區(qū)范圍將逐漸增大.

    取q0=10MPa,pw0=1MPa,φ=30°,c分別為0.25MPa,0.505MPa,0.75MPa,1.00MPa,pa與R0/r1的關(guān)系如圖3(b)所示.可以看出,R0/r1一定時(shí),所需支護(hù)阻力隨著c的減小而增大;反之,當(dāng)pa一定時(shí),隨著c的減小,塑性區(qū)范圍將逐漸增大.

    圖3說(shuō)明,徑向支護(hù)阻力限制了塑性區(qū)域的發(fā)展,這是支護(hù)阻力一個(gè)很重要的作用.

    圖3 徑向支護(hù)阻力與圍巖塑性區(qū)大小的關(guān)系Fig.3 Relationship between radical support pressure and radius of plastic zone of surrounding rock

    2.4 圍巖的位移分析

    彈性區(qū)內(nèi)的圍巖位移可采用彈性力學(xué)中厚壁圓筒的公式確定,即

    式中,μr為圍巖泊松比,Er為圍巖彈性模量.將圍巖視為在彈、塑性邊界上承受內(nèi)壓力σR0的厚壁圓筒,則彈性區(qū)的位移為

    由圍巖彈、塑性邊界上的位移協(xié)調(diào)條件可知

    假定圍巖在塑性區(qū)內(nèi)產(chǎn)生變形的過(guò)程中體積不變,則圍巖周邊產(chǎn)生的位移近似等于

    聯(lián)立式(15)、式(19)和式(23),可得隧道周邊位移關(guān)于支護(hù)阻力pa的關(guān)系式

    由于式(19)為關(guān)于R0的超越方程,無(wú)法寫(xiě)出R0關(guān)于pa的顯式表達(dá)式,故式(24)也無(wú)法寫(xiě)出與pa的顯式表達(dá)式.對(duì)于此類(lèi)問(wèn)題,一般可由式(19)采用迭代法或試算法先求得某一pa作用時(shí)的R0,再代入式(15)求得σR0,最后代入式(23)求得此pa作用時(shí)的由于采用迭代法或試算法求解較為繁瑣,因此本文擬采用更簡(jiǎn)捷的方法進(jìn)行求解,即將R0作為溝通pa和的“橋梁”,先對(duì)R0賦某一值,通過(guò)式(19)、式(15)求得相應(yīng)的pa和σR0,再通過(guò)式(23)求得當(dāng)R0取值的增量足夠小時(shí),對(duì)應(yīng)的數(shù)值解趨近于精確解,從而使計(jì)算過(guò)程得到簡(jiǎn)化.

    3 圍巖和襯砌結(jié)構(gòu)之間平衡的建立

    3.1 圍巖彈塑性收斂曲線

    圍巖收斂曲線(GRC)是對(duì)不同支護(hù)阻力作用下洞周?chē)鷰r變形的描述,反映了洞周?chē)鷰r位移隨支護(hù)力減小而增大的關(guān)系,其典型形態(tài)如圖4所示[27-28].

    圖4 圍巖收斂曲線和襯砌支護(hù)特性曲線Fig.4 Ground reaction curve(GRC)and support characteristic curve(SCC)

    若定義隧道開(kāi)挖后圍巖周邊不出現(xiàn)塑性破裂的最低支護(hù)阻力為臨界阻力,則當(dāng)支護(hù)結(jié)構(gòu)能夠提供的支護(hù)阻力不小于臨界阻力時(shí),開(kāi)挖后圍巖仍處于彈性狀態(tài),洞周位移與支護(hù)阻力的關(guān)系曲線可用式(21)描述.當(dāng)支護(hù)結(jié)構(gòu)能夠提供的支護(hù)阻力小于臨界阻力時(shí),開(kāi)挖后圍巖出現(xiàn)塑性區(qū),若考慮滲流效應(yīng)并基于D-P屈服準(zhǔn)則,此時(shí)洞周位移與支護(hù)阻力的關(guān)系可用式(24)描述,其對(duì)應(yīng)的曲線可通過(guò)2.4節(jié)介紹的簡(jiǎn)化方法處理.

    3.2 支護(hù)特性曲線

    支護(hù)特性曲線(SCC)是指作用在支護(hù)結(jié)構(gòu)上的載荷與支護(hù)結(jié)構(gòu)變形的關(guān)系曲線,當(dāng)不考慮支護(hù)結(jié)構(gòu)與圍巖的接觸狀態(tài)對(duì)支護(hù)結(jié)構(gòu)剛度的影響時(shí),可認(rèn)為作用在支護(hù)結(jié)構(gòu)上的徑向壓力和它的徑向位移成正比.通常支護(hù)結(jié)構(gòu)都是在隧道圍巖已經(jīng)發(fā)生一定量值的收斂變形后才施設(shè)的,并一直作用到支護(hù)結(jié)構(gòu)強(qiáng)度為止,如上圖4所示.支護(hù)特性曲線可用下式表示

    式中,ur0為支護(hù)結(jié)構(gòu)的徑向位移,u0為支護(hù)結(jié)構(gòu)的初始徑向位移,pa為作用在支護(hù)結(jié)構(gòu)上的徑向壓力,Kc為支護(hù)結(jié)構(gòu)剛度.

    當(dāng)混凝土支護(hù)結(jié)構(gòu)厚度ds較小時(shí) (ds≤0.04r1),可采用薄壁圓筒公式計(jì)算,其支護(hù)剛度和可提供的最大支護(hù)阻力pamax分別為

    式中,pamax為支護(hù)結(jié)構(gòu)的最大支護(hù)阻力;Ec為混凝土彈性模量;μc為混凝土泊松比;fc為混凝土抗壓強(qiáng)度.當(dāng)混凝土支護(hù)結(jié)構(gòu)厚度ds(ds≥0.04r1)較大時(shí),可采用厚壁圓筒公式計(jì)算,其支護(hù)剛度和可提供的最大支護(hù)阻力pamax分別為[5]

    對(duì)于裝配式管片襯砌結(jié)構(gòu),其厚度約為外徑的4%~6%[29-30],屬厚壁圓筒問(wèn)題.由于接頭的存在,管片襯砌結(jié)構(gòu)整體抗彎性能有所降低,可通過(guò)引入結(jié)構(gòu)橫向剛度有效率η(剛度折減系數(shù))[31]來(lái)表征這一特性,于是式(28)可修正為下式

    需要指出,前文對(duì)圍巖進(jìn)行彈塑性分析時(shí),考慮了滲流作用,得到的圍巖收斂曲線方程中含有滲透水壓力項(xiàng).然而在得出管片襯砌支護(hù)特性曲線時(shí)采用的是傳統(tǒng)的分析方法,尚未考慮滲透水壓力的影響,這顯然不盡合理.因此此處采用簡(jiǎn)化方法考慮這一作用,認(rèn)為滲透水壓力使管片結(jié)構(gòu)產(chǎn)生初始位移,即其擴(kuò)大了圍巖與管片結(jié)構(gòu)相互作用前二者的“間隙”.

    文獻(xiàn)[32]指出,由于襯砌厚度相對(duì)較小,可以把以體積力形式作用在襯砌范圍內(nèi)的滲流力以合力的形式表達(dá)為作用在襯砌背后的表面力.對(duì)于本文,其值可通過(guò)對(duì)式(5)第2式積分得到

    由厚壁圓筒僅受外壓作用時(shí)的位移公式可得,滲透水壓作用下管片襯砌的附加徑向位移為

    將上式計(jì)算結(jié)果疊加到式(25)的u0中,即可在管片襯砌支護(hù)特性曲線中考慮滲流作用的影響.

    3.3 圍巖與襯砌結(jié)構(gòu)靜力平衡狀態(tài)的建立

    圍巖開(kāi)挖剛完成時(shí),若要保持圍巖穩(wěn)定,需要很大的支護(hù)阻力才能夠?qū)崿F(xiàn),往往遠(yuǎn)大于支護(hù)結(jié)構(gòu)所能提供的最大支護(hù)力,因此,圍巖會(huì)繼續(xù)發(fā)生收斂變形.與此同時(shí),由于支護(hù)結(jié)構(gòu)徑向變形逐漸增大,所能提供的支護(hù)阻力也隨之增大.最終,圍巖收斂曲線將與管片襯砌的支護(hù)特性曲線相交于一點(diǎn),達(dá)到平衡狀態(tài).交點(diǎn)橫坐標(biāo)對(duì)應(yīng)平衡體系形成時(shí)的協(xié)調(diào)位移,也即洞周發(fā)生的位移;交點(diǎn)縱坐標(biāo)為管片襯砌結(jié)構(gòu)承擔(dān)的載荷,其以上部分為圍巖自身需要承擔(dān)的載荷.

    4 工程實(shí)例分析

    獅子洋隧道位于廣深港客運(yùn)專(zhuān)線東涌站至虎門(mén)站之間,穿越珠江入??诘莫{子洋.隧道工程范圍全長(zhǎng)10800m,其中盾構(gòu)段9340m.隧道主體結(jié)構(gòu)采用單層裝配式管片襯砌結(jié)構(gòu),內(nèi)徑9.8m,外徑10.8m,管片厚度50cm,平均幅寬為2m,采用“7+1”分塊方式的通用楔形環(huán)鋼筋混凝土管片錯(cuò)縫襯砌,混凝土強(qiáng)度等級(jí)為C50,抗?jié)B等級(jí)為P12[33].

    在右線基巖淺埋段設(shè)置試驗(yàn)斷面測(cè)試襯砌承受的水土壓力.試驗(yàn)斷面隧道大范圍通過(guò)弱風(fēng)化泥質(zhì)粉砂巖,隧道埋深約為38m(超過(guò)2倍隧道外直徑),水頭約為38m.實(shí)測(cè)施工期和穩(wěn)定后作用在管片襯砌上的水壓力分別為0.33MPa,0.39MPa,圍巖壓力分別為0.10MPa,0.16MPa[25].因此,施工期流固耦合效應(yīng)引起的水壓折減系數(shù)可暫取0.85.

    監(jiān)測(cè)斷面處圍巖及管片襯砌結(jié)構(gòu)的幾何和力學(xué)參數(shù)根據(jù)文獻(xiàn)[33-35]取值.管片結(jié)構(gòu)橫向剛度有效率η根據(jù)文獻(xiàn)[31]取0.70;管片襯砌混凝土抗?jié)B等級(jí)為P12,可認(rèn)為管片不透水;支護(hù)施設(shè)時(shí)圍巖徑向收斂值u0未知,暫由刀盤(pán)開(kāi)挖直徑和管片外徑確定,取u0=191mm[36].

    綜上,計(jì)算所取參數(shù)分別為

    根據(jù)式(21)、式(24)、式(29)和式(30)求解并在同一坐標(biāo)系中繪制圍巖彈塑性收斂曲線和支護(hù)特性曲線,如圖5所示.從圖中可以看出,對(duì)于本例,由于管片襯砌剛度大,而滲透水壓相對(duì)較小,產(chǎn)生的管片襯砌附加變形很小.分別求出施工期和穩(wěn)定期圍巖收斂曲線與支護(hù)特性曲線的交點(diǎn),其橫坐標(biāo)即為管片結(jié)構(gòu)所承擔(dān)的圍巖載荷.本算例中,計(jì)算結(jié)果分別為0.128MPa,0.179MPa.

    圖5 圍巖與襯砌結(jié)構(gòu)的相互作用Fig.5 Interaction between surrounding rock and lining structure

    為進(jìn)一步驗(yàn)證本文方法的有效性,采用FLAC3D建立數(shù)值分析模型,計(jì)算隧道開(kāi)挖后作用在管片襯砌上的圍巖壓力.計(jì)算所取幾何參數(shù)及材料物理力學(xué)參數(shù)均與理論分析取值相同.對(duì)于理論解析中未考慮的管片壁后注漿材料,數(shù)值分析模型中也采用“彈性等代層”,等代層參數(shù)見(jiàn)表1.

    計(jì)算按平面應(yīng)變問(wèn)題考慮,模型中地層、管片襯砌和注漿材料均采用實(shí)體單元模擬,地層材料采用D-P模型,其他材料采用彈性模型.模型上表面為自由邊界,取實(shí)際埋深38m,其他邊界取4D(D為管片外徑),并施加法向約束.初始地應(yīng)力取自重應(yīng)力,初始水壓取靜水壓力,水位與地表平齊.

    表1 等代層參數(shù)Table 1 Parameters of equivalent layer

    由于施工過(guò)程中漿液逐漸硬化,其材料參數(shù)會(huì)產(chǎn)生較大變化,加之漿液在周?chē)貙又辛鲃?dòng)、擴(kuò)散,施工中某一時(shí)刻的漿液參數(shù)很難確定,因此僅列出穩(wěn)定后作用在管片襯砌的圍巖壓力,計(jì)算結(jié)果為0.170MPa.

    采用本文方法計(jì)算得到的結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)值相比,施工期和穩(wěn)定期分別偏大28%和12%;與數(shù)值計(jì)算結(jié)果相比,穩(wěn)定期偏大5%.初步分析導(dǎo)致差異的因素有4點(diǎn):(1)本文推導(dǎo)時(shí)假定圍巖與襯砌結(jié)構(gòu)緊密接觸,而現(xiàn)場(chǎng)實(shí)際并不能保證二者完全密貼,因此測(cè)得的接觸壓力偏小;(2)現(xiàn)場(chǎng)測(cè)試時(shí)管片壁后注漿漿液容易包裹在土壓力盒周?chē)箛鷰r壓力不能完全傳遞至傳感器受力面,造成實(shí)測(cè)壓力值小于實(shí)際值;(3)推導(dǎo)時(shí)假定圍巖初始應(yīng)力場(chǎng)均勻分布,且圍巖為各項(xiàng)同性的理想彈塑性材料,而實(shí)際很難保證;(4)數(shù)值分析模型所取等代層參數(shù)為經(jīng)驗(yàn)值,條件所限而未能根據(jù)實(shí)測(cè)結(jié)果進(jìn)行參數(shù)反演,與實(shí)際情況可能存在偏差.因此,可認(rèn)為本文解析方法得到的結(jié)果可信,并具有較高的精度.但由于施工期流固耦合效應(yīng)造成的水壓降低程度不易量化,采用本文方法計(jì)算施工期管片襯砌承擔(dān)的圍巖壓力有一定偏差,尚需作進(jìn)一步研究.

    5 結(jié)論

    本文基于 Drucker-Prager準(zhǔn)則,針對(duì)盾構(gòu)隧道施工的力學(xué)過(guò)程,等效考慮滲流效應(yīng)影響推導(dǎo)了圍巖與襯砌結(jié)構(gòu)相互作用的彈塑性解析解,包括圍巖彈、塑性區(qū)應(yīng)力與位移、塑性區(qū)半徑等,并通過(guò)建立圍巖與襯砌結(jié)構(gòu)靜力平衡狀態(tài)得到管片襯砌結(jié)構(gòu)承擔(dān)的圍巖壓力,最后通過(guò)算例與典型大埋深水下盾構(gòu)隧道載荷實(shí)測(cè)及數(shù)值計(jì)算值進(jìn)行比較分析,得到以下結(jié)論:

    (1)等效考慮滲流效應(yīng)影響并基于 Drucker-Prager屈服準(zhǔn)則,推導(dǎo)了無(wú)限均質(zhì)各項(xiàng)同性圍巖受等壓作用平面應(yīng)變問(wèn)題的解析解,得到了圍巖彈、塑性區(qū)應(yīng)力與位移、塑性區(qū)半徑等與支護(hù)壓力間的關(guān)系式,彌補(bǔ)了以往研究成果未能反映滲流力影響的缺憾.

    (2)闡述了圍巖與支護(hù)結(jié)構(gòu)靜力平衡狀態(tài)的建立方法,并作為特例,介紹了盾構(gòu)隧道管片襯砌支護(hù)特性曲線的處理方法,即采用結(jié)構(gòu)橫向剛度有效率η對(duì)整環(huán)剛度進(jìn)行折減,以反映接頭對(duì)結(jié)構(gòu)整體剛度的影響.

    (3)選取廣深港客運(yùn)專(zhuān)線獅子洋穿越基巖段水下盾構(gòu)隧道為對(duì)象,采用本文解析方法計(jì)算了施工期和穩(wěn)定期作用在管片襯砌上的圍巖水土壓力.與現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果相比,施工期和穩(wěn)定期的計(jì)算值分別大28%和12%.

    (4)本文解析方法與現(xiàn)場(chǎng)實(shí)測(cè)及數(shù)值計(jì)算結(jié)果存在差異的原因?yàn)椋含F(xiàn)場(chǎng)圍巖與襯砌間并不能完全滿足計(jì)算假定的緊密接觸條件;現(xiàn)場(chǎng)測(cè)試時(shí)管片壁后注漿的包裹作用使測(cè)試結(jié)果只能反映實(shí)際圍巖壓力的一部分;實(shí)際情況很難保證圍巖為無(wú)限均質(zhì)各項(xiàng)同性介質(zhì),并受等壓作用;數(shù)值計(jì)算模型采用的等代層參數(shù)與實(shí)際情況可能存在偏差.

    1 Ita WG.Guidelines for the design of shield tunnel lining.Tunnelling and Underground Space Technology,2000,15(3):303-331

    2 2020年前廣州將新建多條排水防澇深隧支隧.http://www.zgsz.org. cn/2014/0612/11205.html. 2014-6-12(Guangzhou will build a number of deep buried drainage tunnels and anti-waterlogging branch tunnels until 2020.http://news.xinhuanet.com/energy/2016-01/18/c_1117801885.htm.2014-6-12(in Chinese))

    3 黃鵬程.武漢地下將建首條污水深隧.http://news.xinhuanet.com/ energy/2016-01/18/c1117801885.htm. 2016-1-8(Huang Pengcheng.Wuhan will build its firs deep buried sewage tunnel.http:// news.xinhuanet.com/energy/2016-01/18/c1117801885.htm.2016-1-8(in Chinese))

    4 賴(lài)芳杰.掘地30米成都將建78.8公里污水深隧防內(nèi)澇.http:// www.sc.gov.cn/10462/12771/2016/4/3/10374826.shtml.2016-4-3 (Lai Fangjie.Digging 30 m,Chengdu will build 78.8 km of deep buried sewage tunnel to prevent urban waterlogging.http:// www.sc.gov.cn/10462/12771/2016/4/3/10374826.shtml.2016-4-3 (in Chinese))

    5 何川,封坤,方勇.盾構(gòu)法修建地鐵隧道的技術(shù)現(xiàn)狀與展望.西南交通大學(xué)學(xué)報(bào),2015,50(1):97-109(He Chuan,Feng Kun,Fang Yong.Review and prospect on the constructing technologies of metro tunnel using shield tunneling method.Journal of Southwest Jiaotong University,2015,50(1):97-109(in Chinese))

    6 何川.盾構(gòu)/TBM施工煤礦長(zhǎng)距離斜井的技術(shù)挑戰(zhàn)與展望.隧道建設(shè),2014,34(4):287-297(He Chuan.Challenges and prospectives of construction of long distance inclined shafts of coal mines by shield/TBM.Tunnel Construction,2014,34(4):287-297(in Chinese))

    7 《中國(guó)公路學(xué)報(bào)》編輯部.中國(guó)隧道工程學(xué)術(shù)研究綜述·2015.中國(guó)公路學(xué)報(bào),2015,28(5):1-65(Editorial Office of China Journal of Highway and Transport.Review on China’s tunnel engineering research:2015.China Journal of Highway and Transport,2015, 28(5):1-65(in Chinese))

    8 洪開(kāi)榮.我國(guó)隧道及地下工程發(fā)展現(xiàn)狀與展望.隧道建設(shè),2015, 35(2):95-107(Hong Kairong.State-of-art and prospect of tunnels and underground works in China.Tunnel Construction,2015,35(2):95-107(in Chinese))

    9 何川,封坤.大型水下盾構(gòu)隧道結(jié)構(gòu)研究現(xiàn)狀與展望.西南交通大學(xué)學(xué)報(bào),2011,46(1):1-11(He Chuan,Feng Kun.Review and prospect of structure research of underwater shield tunnel with large cross-section.Journal of Southwest Jiaotong University,2011, 46(1):1-11(in Chinese))

    10 齊春,何川,封坤.考慮流固耦合效應(yīng)的水下盾構(gòu)隧道受力特性.西南交通大學(xué)學(xué)報(bào),2015,50(2):306-311,330(Qi Chun,He Chuan,Feng Kun.Fluid-solid interaction-based mechanical characteristics of underwater shield tunnel.Journal of Southwest Jiaotong University,2015,50(2):306-311,330(in Chinese))

    11 夏煒洋.盾構(gòu)法隧道施工期流固耦合問(wèn)題研究.[博士論文].成都:西南交通大學(xué),2012(Xia Weiyang.Study on coupled solidflui problemofshieldtunnelduringconstructionperiod.[PhDThesis].Chengdu:Southwest Jiaotong University,2012(in Chinese))

    12 蔡勇平,蔡曉鴻.水工壓力隧洞結(jié)構(gòu)應(yīng)力計(jì)算.北京:中國(guó)水利水電出版社,2004(Cai Yongping,Cai Xiaohong.Stress Calculation of Hydraulic Pressure Tunnel Structure.Beijing:China Water &Power Press,2004(in Chinese))

    13 李宗利,任青文,王亞紅.考慮滲流場(chǎng)影響深埋圓形隧洞的彈塑性解.巖石力學(xué)與工程學(xué)報(bào),2004,23(8):1291-1295(Li Zongli, Ren Qingwen,Wang Yahong.Elasto-plastic analytical solution of deep-buried circle tunnel considering flui fl w fieldChinese Journal of Rock Mechanics and Engineering,2004,23(8):1291-1295 (in Chinese))

    14 張常光.圓形壓力隧洞彈塑性應(yīng)力和位移分析.[碩士論文].西安:長(zhǎng)安大學(xué),2008.(Zhang Changguang.Analysis of elasticplastic stress and displacement for circular pressure tunnel.[Master Thesis].Xi’an:Chang’an University,2008(in Chinese))

    15 吳順川,潘旦光,高永濤.深埋圓形巷道圍巖和襯砌相互作用解析解.工程力學(xué),2011,28(3):136-142(Wu Shunchuan,Pan Danguang,Gao Yongtao.Analytic solution for rock-liner interaction of deep circular tunnel.Engineering Mechanics,2011,28(3):136-142 (in Chinese))

    16 孫闖,張向東,李永靖.高應(yīng)力軟巖巷道圍巖與支護(hù)結(jié)構(gòu)相互作用分析.巖土力學(xué),2013,34(9):2601-2607(Sun Chuang,Zhang Xiangdong,Li Yongjing.Analysis of interaction between surrounding rock and support structure in high stressed soft rock roadway.Rock and Soil Mechanics,2013,34(9):2601-2607(in Chinese))

    17 鄒金鋒,李帥帥,張勇等.考慮軸向力和滲透力時(shí)軟化圍巖隧道解析.力學(xué)學(xué)報(bào),2014,46(5):747-755(Zou Jinfeng,Li Shuaishuai, Zhang Yong,et al.Solution and analysis of circular tunnel for the strain-softening rock masses considering the axial in situ stress and seepage force.Chinese Journal of Theoretical and Applied Mechanics,2014,46(5):747-755(in Chinese))

    18 Carranza-Torres C,Fairhurst C.The elasto-plastic response of underground excavations in rock masses that satisfy the Hoek-Brown failure criterion.International Journal of Rock Mechanics and Mining Sciences,1999,36(6):777-809

    19 Carranza-Torres C,Fairhurst C.Application of the convergenceconfinemen method of tunnel design to rock masses that satisfy the hoek-brown failure criterion.Tunnelling and Underground Space Technology,2000,15(2):187-213

    20 蘇永華,劉少峰,王凱旋等.基于收斂--約束原理的地下結(jié)構(gòu)穩(wěn)定性分析.巖土工程學(xué)報(bào),2014,36(11):2002-2009(Su Yonghua, Liu Shaofeng,Wang Kaixuan,et al.Stability analysis of underground structures based on convergence-confinemen method.Chinese Journal of Geotechnical Engineering,2014,36(11):2002-2009(in Chinese))

    21侯公羽,牛曉松.基于Levy-Mises本構(gòu)關(guān)系及D-P屈服準(zhǔn)則的軸對(duì)稱(chēng)圓巷理想彈塑性解.巖土力學(xué),2009,30(6):1555-1562(Hou Gongyu,NiuXiaosong.Perfectelastoplasticsolutionofaxisymmetric circular openings in rock mass based on Levy-Mises constitutive relation and D-P yield criterion.Rock and Soil Mechanics,2009, 30(6):1555-1562(in Chinese))

    22張小波,趙光明,孟祥瑞.基于Drucker-Prager屈服準(zhǔn)則的圓形巷道圍巖彈塑性分析.煤炭學(xué)報(bào),2013,38(S1):30-37(Zhang Xiaobo,Zhao Guangming,Meng Xiangrui.Elastoplastic analysis of surrounding rock on circular roadway based on Drucker-Prager yield criterion.Journal of China Coal Society,2013,38(S1):30-37 (in Chinese))

    23 謝紅強(qiáng).隧道工程熱液固多場(chǎng)耦合效應(yīng)研究.[博士論文].成都:西南交通大學(xué),2006(Xie Hongqiang.Study on multi-fiel coupled ef f ects of heat,liquid and solid in tunnel works.[PhD Thesis]. Chengdu:Southwest Jiaotong University,2006(in Chinese))

    24 何川,謝紅強(qiáng).多場(chǎng)耦合分析在隧道工程中的應(yīng)用.成都:西南交通大學(xué)出版社,2007(He Chuan,Xie Hongqiang.Application of Analysis of Multi-fiel Coupled Ef f ects in Tunnel Works.Chengdu:Southwest Jiaotong University Press,2007(in Chinese))

    25 肖明清.大型水下盾構(gòu)隧道結(jié)構(gòu)設(shè)計(jì)關(guān)鍵問(wèn)題研究.[博士論文].成都:西南交通大學(xué),2014(Xiao Mingqing.Research on key issues of segmental lining structure design for underwater shield tunnel with large cross-section.[PhD Thesis].Chengdu:Southwest Jiaotong University,2014(in Chinese))

    26 鄭穎人,沈珠江,龔曉南.巖土塑性力學(xué)原理——廣義塑性力學(xué).北京:中國(guó)建筑工業(yè)出版社,2002(Zheng Yingren,Shen Zhujiang, Gong Xiaonan.Principle of Plastic Mechanics of Rock and Soil——Generalized Plastic Mechanics.Beijing:China Building Industry Press,2002(in Chinese))

    27 Alejano LR,Alonso E,Rodr′?guez-Dono A,et al.Application of the convergence-confinemen method to tunnels in rock masses exhibiting Hoek-Brown strain-softening behaviour.International Journal of Rock Mechanics and Mining Sciences,2010,47(1):150-160

    28 Gonz′alez-Nicieza C,′alvarez-Vigil AE,Men′endez-D′?az A,et al.Influencof the depth and shape of a tunnel in the application of theconvergence-confinemen method.TunnellingandUnderground Space Technology,2008,23(1):25-37

    29 劉建航,侯學(xué)淵.盾構(gòu)法隧道.北京:中國(guó)鐵道出版社,1991(Liu Jianhang,Hou Xueyuan.Shield Tunnel.Beijing:China Railway Publishing House,1991(in Chinese))

    30 日本土木學(xué)會(huì).盾構(gòu)隧道管片設(shè)計(jì)——從容許應(yīng)力法到極限狀態(tài)法.官林星譯.北京:中國(guó)建筑工業(yè)出版社,2012(JSCE.Design of Shield Tunnel Segment——From Allowable Stresses Method to Limit State Method.Guan LX transl.Beijing:China Building Industry Press,2012(in Chinese))

    31 封坤,何川,夏松林.大斷面盾構(gòu)隧道結(jié)構(gòu)橫向剛度有效率的原型試驗(yàn)研究.巖土工程學(xué)報(bào),2011,33(11):1750-1758(Feng Kun, He Chuan,Xia Songlin.Prototype tests on ef f ective bending rigidity ratios of segmental lining structure for shield tunnel with large cross-section.Chinese Journal of Geotechnical Engineering,2011, 33(11):1750-1758(in Chinese))

    32 王建宇.隧道圍巖滲流和襯砌水壓力載荷.鐵道建筑技術(shù),2008, (2):1-6(Wang Jianyu.Problems on external water pressure on tunnel lining.Railway Construction Technology,2008,(2):1-6(in Chinese))

    33 鐵道第四勘察設(shè)計(jì)院.廣深港客運(yùn)專(zhuān)線廣州至深圳段獅子洋隧道盾構(gòu)隧道結(jié)構(gòu)設(shè)計(jì).武漢:鐵道第四勘察設(shè)計(jì)院,2006(China Railway Siyuan Survey and Design Group Co.,Ltd.Structure design of Shiziyang shield tunnel of Guangzhou-Shenzhen section of Guangzhou-Shenzhen-Hong Kong PDL.Wuhan:China Railway Siyuan Survey and Design Group Co.,Ltd,2006(in Chinese))

    34 鐵道第四勘察設(shè)計(jì)院.廣深港客運(yùn)專(zhuān)線廣州至深圳段獅子洋隧道工程地質(zhì)勘察報(bào)告.武漢:鐵道第四勘察設(shè)計(jì)院,2006(China Railway Siyuan Survey and Design Group Co.,Ltd.Geological survey report of Shiziyang shield tunnel of Guangzhou-Shenzhen section of Guangzhou-Shenzhen-Hong Kong PDL.Wuhan:China Railway Siyuan Survey and Design Group Co.,Ltd,2006(in Chinese))

    35 中華人民共和國(guó)鐵道部.TB10003-2005,J449-2005鐵路隧道設(shè)計(jì)規(guī)范.北京:中國(guó)鐵道出版社,2005(Ministry of Railways of the PRC.TB10003-2005,J449-2005 Code for design of Railway Tunnels.Beijing:China Railway Publishing House,2005(in Chinese))

    36 洪開(kāi)榮.高速鐵路特長(zhǎng)水下盾構(gòu)隧道施工技術(shù).北京:中國(guó)鐵道出版社,2013(Hong Kairong.Construction Technologies of Extra Long Underwater Shield-Bored High Speed Railway Tunnels.Beijing:China Railway Publishing House,2013(in Chinese))

    THEORETICAL ANALYSIS OF INTERACTION BETWEEN SURROUNDING ROCKS AND LINGING STRCTURE OF SHIELD TUNNEL BASED ON DRUCKER-PRAGER YIELD CRITERIA1)

    He Chuan?,2)Qi Chun?Feng Kun?,3)Xiao Mingqing?,4)

    ?(Key Laboratory of Transportation Tunnel Engineering,Ministry of Education,Southwest Jiaotong University,Chengdu610031,China)

    ?(China Railway Siyuan Survey and Design Group Co.,Ltd,Wuhan430063,China)

    In calculating the load for segment lining structure of shield tunnel,the soil column or pressure arch theory are often used to describe surrounding rock loosening pressure.But when it comes to the deep-buried condition and the deformation pressure should be calculated,this method is considered as unreasonable and is difficult to apply.In view ofthis,based on the Drucker-Prager yield criterion,an analytical elasto plastic solution is derived for the interaction between surrounding rock and lining structure considering the influenc of seepage ef f ect,and the formula of the relationship between the support pressure and some key parameters such as the stress and displacement of surrounding rock elastic and plastic zone as well as the radius of the plastic zone is given.The analytical results can be applied in the determination oftheloadoftheliningstructurebyestablishingthestaticequilibriumstateofthesurroundingrockandtheliningstructure and findin the intersection point of the two curves.Further,rigidity reduction factor is induced to equivalent considering the influenc of the overall rigidity decrease of assembled segment liner for the existence of joints to the interaction of surrounding rock and lining structure.In the load determination,the influenc of the fluid-soli coupling ef f ect during construction and seepage force on support characteristic curve of lining structure is also simplifie considered.Finally, an engineering instance of underwater shield tunnel is introduced to compare the calculated load with the measured value and numerical simulated value.The results show that the caculated analytical load values of segment lining structure are greater by 28%and 12%comparing to measured values during the construction period and the stable period,and are greater by 5%comparing to numerical simulated value during the stable period,respectively.The research results of this paper can provide some reference for the design and construction of similar projects.

    shield tunnel,surrounding rock,lining structure,interaction,surrounding rock pressure,Drucker-Prager criterion,convergence-confinemen method

    U25,U451

    A doi:10.6052/0459-1879-16-344

    2016-11-25收稿,2016-12-04錄用,2016-12-04網(wǎng)絡(luò)版發(fā)表.

    1)國(guó)家自然科學(xué)基金(U1361210,51578462)、國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0802200)、中央高校基本科研業(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金(2682015CX077)和中國(guó)鐵路總公司科技研究開(kāi)發(fā)計(jì)劃(2014G004-0)資助項(xiàng)目.

    2)何川,教授,主要研究方向:地鐵及水下盾構(gòu)隧道結(jié)構(gòu)理論、大型及復(fù)雜交通隧道結(jié)構(gòu)安全及長(zhǎng)大交通隧道運(yùn)營(yíng)控制等. E-mail:chuanhe21@163.com

    3)封坤,副教授,主要研究方向:現(xiàn)代盾構(gòu)隧道技術(shù)與結(jié)構(gòu)設(shè)計(jì)理論.E-mail:windfeng813@163.com

    4)肖明清,教授級(jí)高級(jí)工程師,主要研究方向:水下隧道及高速鐵路隧道設(shè)計(jì).E-mail:tsyxmq@163.com

    何川,齊春,封坤,肖明清.基于D-P準(zhǔn)則的盾構(gòu)隧道圍巖與襯砌結(jié)構(gòu)相互作用分析.力學(xué)學(xué)報(bào),2017,49(1):31-40

    He Chuan,Qi Chun,Feng Kun,Xiao Mingqing.Theoretical analysis of interaction between surrounding rocks and linging strcture of shield tunnel based on Drucker-Prager yield criteria.Chinese Journal of Theoretical and Applied Mechanics,2017,49(1):31-40

    猜你喜歡
    管片滲流塑性
    基于應(yīng)變梯度的微尺度金屬塑性行為研究
    硬脆材料的塑性域加工
    鈹材料塑性域加工可行性研究
    管片拼裝技術(shù)研究
    盾構(gòu)管片封頂塊拼裝施工技術(shù)研究
    石英玻璃的熱輔助高效塑性域干磨削
    地鐵盾構(gòu)管片受力分析及管片破損的控制措施研究
    簡(jiǎn)述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計(jì)算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    尾礦壩滲流計(jì)算及排滲設(shè)計(jì)
    金屬礦山(2013年6期)2013-03-11 16:54:05
    亚洲美女视频黄频| 99热这里只有精品一区| 久久久久免费精品人妻一区二区| ponron亚洲| 亚洲av男天堂| 亚洲精品影视一区二区三区av| 免费看av在线观看网站| 久久鲁丝午夜福利片| 欧美+亚洲+日韩+国产| 国产亚洲91精品色在线| 欧美色视频一区免费| 级片在线观看| 国产老妇女一区| 男女下面进入的视频免费午夜| 久久久久久久午夜电影| 精品人妻熟女av久视频| 亚洲国产日韩欧美精品在线观看| 99国产极品粉嫩在线观看| 亚洲内射少妇av| 久久午夜亚洲精品久久| 美女被艹到高潮喷水动态| 欧美三级亚洲精品| 99久久中文字幕三级久久日本| 一个人观看的视频www高清免费观看| 久久精品国产99精品国产亚洲性色| 在线国产一区二区在线| 欧美日本亚洲视频在线播放| 欧美一区二区亚洲| 久久精品国产99精品国产亚洲性色| 亚洲无线观看免费| 日产精品乱码卡一卡2卡三| 国产一区二区亚洲精品在线观看| 一级黄片播放器| 国产伦精品一区二区三区四那| 亚洲成人久久性| 国产一区亚洲一区在线观看| 在线免费十八禁| 亚洲在线自拍视频| 亚洲一区二区三区色噜噜| 免费不卡的大黄色大毛片视频在线观看 | 亚洲不卡免费看| 成人高潮视频无遮挡免费网站| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区三区| 一进一出抽搐动态| 欧美高清成人免费视频www| 中文字幕av成人在线电影| 亚洲内射少妇av| 麻豆成人午夜福利视频| 97超碰精品成人国产| 波多野结衣高清无吗| 成年女人永久免费观看视频| 嫩草影院新地址| 久久国内精品自在自线图片| 国产精品久久视频播放| 校园人妻丝袜中文字幕| 免费人成视频x8x8入口观看| 在线观看午夜福利视频| 亚洲精品456在线播放app| 精品久久久久久成人av| 国产大屁股一区二区在线视频| 看片在线看免费视频| 日韩,欧美,国产一区二区三区 | 国产女主播在线喷水免费视频网站 | 成人性生交大片免费视频hd| 午夜久久久久精精品| 人体艺术视频欧美日本| 九色成人免费人妻av| 国产精品久久久久久精品电影小说 | 尾随美女入室| 国产不卡一卡二| 日韩av在线大香蕉| 老师上课跳d突然被开到最大视频| 小说图片视频综合网站| 国产成人精品婷婷| 99热这里只有是精品在线观看| 欧美日本视频| 国产高清激情床上av| 丝袜美腿在线中文| 久久久久久久久久久丰满| 国产老妇女一区| 美女黄网站色视频| 国产真实乱freesex| 午夜视频国产福利| 少妇人妻一区二区三区视频| 夜夜爽天天搞| 天堂网av新在线| 成年版毛片免费区| 人人妻人人澡欧美一区二区| 99国产精品一区二区蜜桃av| 精品无人区乱码1区二区| 日韩大尺度精品在线看网址| av天堂在线播放| 超碰av人人做人人爽久久| 男人狂女人下面高潮的视频| 小说图片视频综合网站| 国产黄a三级三级三级人| 可以在线观看的亚洲视频| 久久久欧美国产精品| 亚洲真实伦在线观看| 青青草视频在线视频观看| 99热这里只有是精品在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 99热网站在线观看| 国产精品嫩草影院av在线观看| 免费一级毛片在线播放高清视频| 哪个播放器可以免费观看大片| 免费人成视频x8x8入口观看| 国产精品麻豆人妻色哟哟久久 | av福利片在线观看| 国产精品三级大全| 又爽又黄无遮挡网站| 国产成人freesex在线| 日韩精品有码人妻一区| 亚洲欧美清纯卡通| 小蜜桃在线观看免费完整版高清| 欧美日韩在线观看h| 亚洲欧美日韩卡通动漫| 国产色婷婷99| 欧美日本亚洲视频在线播放| 能在线免费看毛片的网站| 啦啦啦观看免费观看视频高清| 日韩一本色道免费dvd| 成人漫画全彩无遮挡| 欧美成人免费av一区二区三区| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 听说在线观看完整版免费高清| 免费人成视频x8x8入口观看| 麻豆成人午夜福利视频| 国产91av在线免费观看| 在线播放无遮挡| 国内少妇人妻偷人精品xxx网站| 美女国产视频在线观看| 日日干狠狠操夜夜爽| 韩国av在线不卡| 国产 一区 欧美 日韩| 国产毛片a区久久久久| 久久99蜜桃精品久久| 色哟哟·www| 99riav亚洲国产免费| 亚洲av电影不卡..在线观看| 看黄色毛片网站| 色哟哟哟哟哟哟| 国产私拍福利视频在线观看| 午夜亚洲福利在线播放| 国产老妇伦熟女老妇高清| 久久人妻av系列| 久久精品国产亚洲av涩爱 | 亚洲成人久久爱视频| 99热这里只有是精品50| 波野结衣二区三区在线| 精华霜和精华液先用哪个| 欧美一区二区亚洲| 麻豆乱淫一区二区| 国产av麻豆久久久久久久| 夜夜看夜夜爽夜夜摸| 久久精品国产亚洲av香蕉五月| 日韩成人av中文字幕在线观看| 欧美色欧美亚洲另类二区| 午夜久久久久精精品| 网址你懂的国产日韩在线| 久久99蜜桃精品久久| 久久久久久国产a免费观看| 一区二区三区免费毛片| 又爽又黄无遮挡网站| 91精品一卡2卡3卡4卡| 12—13女人毛片做爰片一| 亚洲综合色惰| 69av精品久久久久久| 美女高潮的动态| 看片在线看免费视频| av在线蜜桃| 99久久中文字幕三级久久日本| 日韩av在线大香蕉| 一级av片app| kizo精华| av天堂在线播放| 亚洲精品国产成人久久av| 成人午夜高清在线视频| 亚洲av免费高清在线观看| 久久精品国产鲁丝片午夜精品| 欧美精品一区二区大全| 爱豆传媒免费全集在线观看| 亚洲高清免费不卡视频| 亚洲欧美精品综合久久99| 日韩一区二区视频免费看| 亚洲人成网站在线播放欧美日韩| 热99在线观看视频| 欧美zozozo另类| 成人亚洲精品av一区二区| 免费看光身美女| 久久人人精品亚洲av| 九九热线精品视视频播放| 亚洲av不卡在线观看| 夜夜看夜夜爽夜夜摸| 99九九线精品视频在线观看视频| 联通29元200g的流量卡| 久久久成人免费电影| 国产美女午夜福利| 国产亚洲精品久久久久久毛片| 日韩强制内射视频| 国内精品久久久久精免费| 日韩欧美精品v在线| 春色校园在线视频观看| 成人高潮视频无遮挡免费网站| 欧美精品国产亚洲| av黄色大香蕉| 亚洲成人中文字幕在线播放| 国产真实乱freesex| 大香蕉久久网| 欧美性猛交╳xxx乱大交人| 又爽又黄a免费视频| 亚洲真实伦在线观看| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 床上黄色一级片| 免费黄网站久久成人精品| 国产精品综合久久久久久久免费| 性插视频无遮挡在线免费观看| av专区在线播放| 少妇丰满av| 熟女人妻精品中文字幕| 你懂的网址亚洲精品在线观看 | 干丝袜人妻中文字幕| 成年女人看的毛片在线观看| 夜夜爽天天搞| 免费看光身美女| 国产成人福利小说| 国产又黄又爽又无遮挡在线| avwww免费| 欧美xxxx性猛交bbbb| 国产精品1区2区在线观看.| 国产一区二区亚洲精品在线观看| 久久精品国产鲁丝片午夜精品| 2022亚洲国产成人精品| 午夜爱爱视频在线播放| 99久久九九国产精品国产免费| 26uuu在线亚洲综合色| 黄色视频,在线免费观看| 亚洲成a人片在线一区二区| av在线亚洲专区| 日本-黄色视频高清免费观看| 久久久久免费精品人妻一区二区| 人妻制服诱惑在线中文字幕| 国内揄拍国产精品人妻在线| 亚洲欧美日韩卡通动漫| 欧美成人免费av一区二区三区| 国产黄片美女视频| av免费在线看不卡| 六月丁香七月| 欧洲精品卡2卡3卡4卡5卡区| 亚洲色图av天堂| 久久久久久久久大av| 一级毛片我不卡| 中文字幕免费在线视频6| 国产高清不卡午夜福利| 久久久成人免费电影| 99久久精品国产国产毛片| 欧美最黄视频在线播放免费| 此物有八面人人有两片| 美女xxoo啪啪120秒动态图| 99久国产av精品| 人人妻人人澡人人爽人人夜夜 | 国产在线精品亚洲第一网站| 男女做爰动态图高潮gif福利片| 99热6这里只有精品| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 看片在线看免费视频| 丰满人妻一区二区三区视频av| 成熟少妇高潮喷水视频| 村上凉子中文字幕在线| 青春草国产在线视频 | 插阴视频在线观看视频| 欧美最新免费一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲欧美日韩东京热| 国产精品日韩av在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 精品人妻视频免费看| 国产精品一区www在线观看| 精品不卡国产一区二区三区| 一个人观看的视频www高清免费观看| 中国美白少妇内射xxxbb| 变态另类丝袜制服| 啦啦啦观看免费观看视频高清| 免费av毛片视频| 国产一级毛片在线| 中文欧美无线码| 可以在线观看的亚洲视频| 亚洲第一区二区三区不卡| 亚洲精品影视一区二区三区av| 国产免费男女视频| 一级二级三级毛片免费看| 免费一级毛片在线播放高清视频| 精品久久久久久久久av| 亚洲欧美日韩高清在线视频| 亚洲18禁久久av| av黄色大香蕉| 亚洲av.av天堂| 校园春色视频在线观看| 五月伊人婷婷丁香| 国产日本99.免费观看| 成人毛片60女人毛片免费| 成人无遮挡网站| 国产精品.久久久| 日本免费a在线| 中文字幕av在线有码专区| 精品午夜福利在线看| 色尼玛亚洲综合影院| 亚洲精品成人久久久久久| 亚洲成a人片在线一区二区| 久久精品人妻少妇| 好男人在线观看高清免费视频| av福利片在线观看| 能在线免费看毛片的网站| 亚洲色图av天堂| 狠狠狠狠99中文字幕| 国产午夜精品论理片| 中文字幕熟女人妻在线| 18+在线观看网站| 男女边吃奶边做爰视频| 小蜜桃在线观看免费完整版高清| 淫秽高清视频在线观看| 91精品一卡2卡3卡4卡| a级毛片免费高清观看在线播放| 国产大屁股一区二区在线视频| 午夜福利高清视频| 久久久久久国产a免费观看| 又粗又爽又猛毛片免费看| 美女黄网站色视频| 欧美在线一区亚洲| 18禁在线播放成人免费| 韩国av在线不卡| 在现免费观看毛片| 免费观看在线日韩| 久久国内精品自在自线图片| 国内精品美女久久久久久| 亚洲欧美清纯卡通| 亚洲一区二区三区色噜噜| 国产亚洲av嫩草精品影院| 久久久久久久久久久免费av| 嫩草影院入口| 国产精品1区2区在线观看.| 99国产精品一区二区蜜桃av| 免费无遮挡裸体视频| 国产高清有码在线观看视频| 国产毛片a区久久久久| 日本av手机在线免费观看| 国产成人一区二区在线| 亚洲av第一区精品v没综合| 日韩成人伦理影院| 免费观看人在逋| 91狼人影院| 成人漫画全彩无遮挡| 久久亚洲国产成人精品v| 我要看日韩黄色一级片| 内射极品少妇av片p| av在线天堂中文字幕| 人人妻人人看人人澡| 九九热线精品视视频播放| 看黄色毛片网站| 91久久精品国产一区二区三区| 午夜免费男女啪啪视频观看| 久久久精品94久久精品| 中国美白少妇内射xxxbb| 99九九线精品视频在线观看视频| 日韩国内少妇激情av| 精品免费久久久久久久清纯| 婷婷亚洲欧美| 能在线免费看毛片的网站| 欧美区成人在线视频| 午夜福利在线观看吧| 日本成人三级电影网站| 日本一二三区视频观看| 一区二区三区高清视频在线| 久久精品夜色国产| 久久99精品国语久久久| 性欧美人与动物交配| 国产真实伦视频高清在线观看| 欧美丝袜亚洲另类| 久久国产乱子免费精品| 九九爱精品视频在线观看| 草草在线视频免费看| 亚洲不卡免费看| 色综合亚洲欧美另类图片| 国产亚洲av片在线观看秒播厂 | 国产免费男女视频| 99久久无色码亚洲精品果冻| 欧美日韩综合久久久久久| 亚洲精品国产av成人精品| 国产精品久久久久久久电影| 国产真实乱freesex| 国产乱人视频| 干丝袜人妻中文字幕| 哪个播放器可以免费观看大片| 成熟少妇高潮喷水视频| 亚洲av.av天堂| 99热这里只有是精品50| 亚洲精品粉嫩美女一区| 尾随美女入室| 卡戴珊不雅视频在线播放| 在线观看66精品国产| 在线观看午夜福利视频| 国产免费一级a男人的天堂| 亚洲国产精品成人久久小说 | 精品99又大又爽又粗少妇毛片| 久久这里只有精品中国| 亚洲在线观看片| 亚洲欧洲日产国产| 欧美人与善性xxx| 最近中文字幕高清免费大全6| www.av在线官网国产| 欧美成人a在线观看| 亚洲人成网站在线播放欧美日韩| 在线免费十八禁| 午夜福利在线在线| 亚洲av电影不卡..在线观看| 国产精品免费一区二区三区在线| 国产精品电影一区二区三区| av在线天堂中文字幕| 亚洲欧美日韩东京热| 亚洲婷婷狠狠爱综合网| 我的老师免费观看完整版| 可以在线观看的亚洲视频| 国产av在哪里看| 婷婷色av中文字幕| 久久婷婷人人爽人人干人人爱| 天堂影院成人在线观看| 天堂√8在线中文| 麻豆一二三区av精品| 乱人视频在线观看| 看十八女毛片水多多多| 毛片女人毛片| 2022亚洲国产成人精品| 99热这里只有精品一区| 国产精品一区二区三区四区久久| 国产午夜精品一二区理论片| 精品免费久久久久久久清纯| 美女脱内裤让男人舔精品视频 | 亚洲人成网站在线播放欧美日韩| 国产精品一区www在线观看| 蜜臀久久99精品久久宅男| 99久国产av精品| 免费av毛片视频| 青春草视频在线免费观看| 淫秽高清视频在线观看| 国产精品一区二区三区四区免费观看| 国产色爽女视频免费观看| 国产男人的电影天堂91| 国产精品人妻久久久久久| 搡老妇女老女人老熟妇| 99视频精品全部免费 在线| 国产老妇伦熟女老妇高清| 亚洲精品久久久久久婷婷小说 | 国产色婷婷99| 一级毛片久久久久久久久女| 亚洲成人av在线免费| 国产真实乱freesex| 亚洲乱码一区二区免费版| 欧美高清性xxxxhd video| 在线免费观看的www视频| 麻豆一二三区av精品| 国产黄色小视频在线观看| 日本成人三级电影网站| 狂野欧美白嫩少妇大欣赏| 精华霜和精华液先用哪个| 国产综合懂色| 成人特级av手机在线观看| 不卡视频在线观看欧美| 亚洲va在线va天堂va国产| ponron亚洲| 蜜桃亚洲精品一区二区三区| 91麻豆精品激情在线观看国产| 美女被艹到高潮喷水动态| 日本免费a在线| 欧美日韩综合久久久久久| 午夜爱爱视频在线播放| 精品久久久久久久久亚洲| 色尼玛亚洲综合影院| 男女做爰动态图高潮gif福利片| 精品久久久久久成人av| 五月伊人婷婷丁香| 男人狂女人下面高潮的视频| 人妻夜夜爽99麻豆av| 午夜老司机福利剧场| 网址你懂的国产日韩在线| 夜夜爽天天搞| 国产精品一区二区在线观看99 | 国产精品国产高清国产av| 狠狠狠狠99中文字幕| 亚洲国产精品合色在线| 日本爱情动作片www.在线观看| 国产伦在线观看视频一区| 99九九线精品视频在线观看视频| 国产亚洲精品av在线| 亚洲最大成人av| 中国美白少妇内射xxxbb| 欧美一区二区亚洲| 国产伦在线观看视频一区| 亚洲av一区综合| 免费看光身美女| 欧美zozozo另类| 三级经典国产精品| 一区二区三区免费毛片| 在线观看免费视频日本深夜| 一区二区三区免费毛片| av福利片在线观看| 六月丁香七月| .国产精品久久| 欧美在线一区亚洲| 看免费成人av毛片| 69人妻影院| 国产精品国产高清国产av| 女人十人毛片免费观看3o分钟| 97超碰精品成人国产| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添av毛片| 欧美性感艳星| 亚洲性久久影院| 欧美日韩乱码在线| 精品少妇黑人巨大在线播放 | 午夜福利在线观看吧| 在线观看免费视频日本深夜| 午夜福利视频1000在线观看| 久久久久久久久久成人| 我要搜黄色片| 黄色配什么色好看| 欧美成人a在线观看| 欧美区成人在线视频| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 久久精品国产亚洲网站| 神马国产精品三级电影在线观看| 麻豆乱淫一区二区| 国产午夜精品久久久久久一区二区三区| 久久午夜亚洲精品久久| 美女脱内裤让男人舔精品视频 | 高清毛片免费看| 丰满的人妻完整版| 亚洲第一电影网av| 日韩欧美三级三区| 日韩大尺度精品在线看网址| 一边亲一边摸免费视频| 国产精品野战在线观看| 狂野欧美激情性xxxx在线观看| 亚洲精华国产精华液的使用体验 | 色5月婷婷丁香| 国产黄片视频在线免费观看| 少妇高潮的动态图| 午夜福利在线观看吧| 日本-黄色视频高清免费观看| 在线免费十八禁| 日产精品乱码卡一卡2卡三| 嘟嘟电影网在线观看| 亚洲美女视频黄频| 亚洲人成网站高清观看| 国产一区二区亚洲精品在线观看| 亚洲无线观看免费| 免费搜索国产男女视频| av.在线天堂| 国产精品久久久久久久久免| 国产黄a三级三级三级人| 网址你懂的国产日韩在线| 精品无人区乱码1区二区| 99久久中文字幕三级久久日本| 成人永久免费在线观看视频| 最近手机中文字幕大全| av免费观看日本| 一级毛片aaaaaa免费看小| 免费搜索国产男女视频| 99国产极品粉嫩在线观看| av在线老鸭窝| 亚洲欧美日韩卡通动漫| 亚洲精品自拍成人| 久久久a久久爽久久v久久| 国产黄片视频在线免费观看| 精品人妻偷拍中文字幕| 在线播放国产精品三级| 最近视频中文字幕2019在线8| 亚洲第一区二区三区不卡| 中文字幕制服av| 国产一区二区亚洲精品在线观看| 国产老妇女一区| 婷婷色av中文字幕| 欧美精品国产亚洲| 婷婷精品国产亚洲av| 老司机影院成人| 99国产极品粉嫩在线观看| 成人亚洲精品av一区二区| 成人三级黄色视频| 综合色丁香网| 成年版毛片免费区| 久久久国产成人免费| 国产精品福利在线免费观看| 免费在线观看成人毛片| 看非洲黑人一级黄片| av在线老鸭窝| 亚洲av男天堂| 日韩一区二区视频免费看| 免费看美女性在线毛片视频| 国产老妇女一区| 久久久久网色| 中文欧美无线码| 又爽又黄a免费视频| 亚洲无线在线观看| 人人妻人人澡人人爽人人夜夜 | 两性午夜刺激爽爽歪歪视频在线观看| 18禁裸乳无遮挡免费网站照片| 男人舔奶头视频| 国产精品久久电影中文字幕| 亚州av有码| 在线播放无遮挡| 成人无遮挡网站|