摘 要:【目的】軟弱圍巖淺埋大斷面隧道施工會(huì)增大周?chē)h(huán)境的變形風(fēng)險(xiǎn),尋求有效的軟弱圍巖大斷面隧道施工開(kāi)挖技術(shù)具有重要的工程實(shí)踐意義?!痉椒ā坎捎脭?shù)值模擬和現(xiàn)場(chǎng)監(jiān)測(cè)相結(jié)合的手段,分析軟弱圍巖淺埋大斷面隧道施工工法與變形控制技術(shù)?!窘Y(jié)果】結(jié)果表明,三臺(tái)臨時(shí)仰拱法能較好地控制軟弱圍巖淺埋大斷面隧道的變形,采用該工法時(shí)隧道中心處的地表沉降值、拱頂沉降及開(kāi)挖底部最大隆起值均能滿足施工規(guī)范要求,且相應(yīng)的數(shù)值模擬與現(xiàn)場(chǎng)監(jiān)測(cè)具有相似的變化規(guī)律?!窘Y(jié)論】研究成果可為今后類(lèi)似隧道工程設(shè)計(jì)與施工提供技術(shù)支撐。
關(guān)鍵詞:軟弱圍巖;淺埋大斷面隧道;三臺(tái)階臨時(shí)仰拱法;變形控制技術(shù);數(shù)值模擬
中圖分類(lèi)號(hào):U455.4 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1003-5168(2024)19-0048-06
DOI:10.19968/j.cnki.hnkj.1003-5168.2024.19.010
Analysis on Construction Method and Deformation Control Technology of Shallow Buried Large Section Tunnel with Weak Surrounding Rock
LIU Jiacheng
(School of Civil Engineering, Changsha University of Science & Technology, Changsha 410114, China)
Abstract: [Purposes] It is of great practical significance to find effective construction and excavation technology of large section tunnel with weak surrounding rock due to the increasing of the deformation risk of the surrounding environment caused by the construction of that. [Methods] The construction method and deformation control technology of shallow buried large section tunnel with weak surrounding rock are analyzed by means of numerical simulation and field monitoring. [Findings] The results show that the three-step temporary invert methobe3c0c24e2b87c4f47bd2d874728da60d14d0c76604f8e0ddf693d6a6e4dedabd can control the deformation of shallow buried large section tunnel in weak surrounding rock. The surface settlement value, vault settlement and maximum uplift value of excavation bottom at the center of tunnel can meet the requirements of construction specifications when using this method, and the corresponding numerical simulation and field detection have similar variation rules.[Conclusions] The research results can provide technical support for design and construction of similar tunnel engineering in the future.
Keywords: weak surrounding rock; shallow buried large section tunnel; three-step temporary invert method; deformation control technology; numerical simulation
0 引言
大斷面隧道作為現(xiàn)代交通基礎(chǔ)設(shè)施的重要組成部分,其建設(shè)規(guī)模和數(shù)量都在不斷增加[1]。至今,已形成了“復(fù)雜地質(zhì)、大直徑、高水壓、長(zhǎng)距離”的隧道建設(shè)理念,但在處理巖層漏水、破碎等問(wèn)題時(shí),仍會(huì)遇到不少難點(diǎn),尤其是軟弱圍巖大斷面隧道[2-3]。因此,尋求合理有效的軟弱圍巖大斷面隧道施工開(kāi)挖技術(shù)是目前隧道工程亟須解決的問(wèn)題。
大斷面隧道施工方法和小斷面隧道有所區(qū)別,特別是在復(fù)雜條件下的施工方法,不僅影響隧道安全施工,還涉及如何保證施工速度的問(wèn)題[4]。當(dāng)前大斷面隧道的施工方法,主要有全斷面法、臺(tái)階法、雙側(cè)壁導(dǎo)坑法、CD法和CRD法等[5],而針對(duì)軟弱破碎圍巖(Ⅵ、Ⅴ、Ⅳ級(jí))地段的超大斷面隧道施工設(shè)計(jì),常采用雙側(cè)壁導(dǎo)坑、CRD、CD法及三臺(tái)階法作為初期和臨時(shí)支護(hù)的施工方法[6]。例如,崔小鵬等[7]綜合分析了CRD法和三臺(tái)階七步法的優(yōu)缺點(diǎn),并對(duì)相應(yīng)的施工工法進(jìn)行優(yōu)化改造,提出了核心土加臨時(shí)仰拱臺(tái)階法;李波等[8]以拱北隧道為例,制定了適合超大斷面的五臺(tái)階15分部、五臺(tái)階14分部和四臺(tái)階8分部等3種開(kāi)挖工序,提出了與之相匹配的設(shè)計(jì)參數(shù),并最終確定了五臺(tái)階14分部的開(kāi)挖方案為最佳方案。此外,不少學(xué)者借助數(shù)值分析手段評(píng)估了隧道施工過(guò)程中的力學(xué)行為、穩(wěn)定性及施工風(fēng)險(xiǎn)。例如,石熊等[9]模擬了CRD法的施工順序,將數(shù)值結(jié)果與現(xiàn)場(chǎng)監(jiān)測(cè)進(jìn)行了對(duì)比分析,提出了“優(yōu)先施工山體外側(cè),后施工山體內(nèi)側(cè)”的方案,確保了隧道施工過(guò)程的安全穩(wěn)定性;蔣坤等[10]利用離散元數(shù)值軟件分析了CRD法、CD法和雙側(cè)壁導(dǎo)坑法施工的結(jié)構(gòu)力學(xué)響應(yīng),發(fā)現(xiàn)CRD法可以縮短工期,降低成本,增強(qiáng)隧道安全穩(wěn)定性。盡管對(duì)大斷面隧道施工過(guò)程方面的研究取得了豐富的成果,但對(duì)于軟弱破碎圍巖超大斷面隧道施工的研究仍不充分,考慮地質(zhì)條件的不確定性以及有效預(yù)測(cè)和防控隧道施工過(guò)程中的變形仍是隧道工程領(lǐng)域研究的熱點(diǎn)問(wèn)題。
基于此,本研究通過(guò)數(shù)值計(jì)算和現(xiàn)場(chǎng)工程驗(yàn)證相結(jié)合的手段,探討復(fù)雜環(huán)境大斷面隧道三臺(tái)階臨時(shí)仰拱法與變形控制技術(shù),研究成果可為地質(zhì)條件復(fù)雜、施工技術(shù)難度大、環(huán)境保護(hù)要求高及安全管理風(fēng)險(xiǎn)大的大斷面隧道工程施工提供技術(shù)支撐。
1 工程概況
武廣高速金沙洲隧道進(jìn)口至斜井段,地形地貌主要是剝蝕殘丘、沖積平原地貌單元及剝蝕殘丘陵地貌單元。山頂最大標(biāo)高67.626 m,相對(duì)高程67.626~1.303 m。場(chǎng)區(qū)處于南亞熱帶海洋性季風(fēng)氣候區(qū)。武廣高速金沙洲隧道DK2193+340~DK2193+510段共穿越兩個(gè)水塘,平面位置如圖1所示,其中DK2193+337~DK2193+398段魚(yú)塘將水排盡, DK2193+400~DK2193+448魚(yú)塘基本未排水,其中DK2193+440附近埋深僅為3.56 m,淤泥層厚度為1 m,隧道在此處屬于淺埋段。隨后,經(jīng)歷豬圈位置,基礎(chǔ)為回填土,此處拱頂埋深為10 m。DK2193+500以后,隧道穿越一座小山。
金沙洲隧道DK2194+560~DK2194+835區(qū)段地質(zhì)條件處在古生界及新生界巖性地層的區(qū)域內(nèi),場(chǎng)區(qū)的巖土層按其成因主要有第四系人工填土層(Q4ml)、全新統(tǒng)沖積相沉積層(Q4al )、殘積層(Q4el)、石炭系下統(tǒng)(C1)及泥盆系上統(tǒng)(D3)。隧道穿越巖層主要為灰?guī)r、炭質(zhì)頁(yè)巖及強(qiáng)風(fēng)化巖層,巖體破碎程度高,遇水會(huì)軟化崩解,屬于V級(jí)圍巖,其斷面形式如圖2所示,隧道寬度超過(guò)10 m,初期支護(hù)28 cm,邊墻和仰拱二次襯砌分別為50 cm和60 cm,以0.8 m×1.0 m(環(huán)×縱)形式布置長(zhǎng)度4.0 m的錨桿,全環(huán)架設(shè)Ⅰ20型鋼臨時(shí)仰拱鋼架。由于隧道在此處穿越水塘,屬于富水軟弱圍巖地層,隧道開(kāi)挖過(guò)程中圍巖極易坍塌,施工風(fēng)險(xiǎn)較大。
2 施工方案及工法
2.1 施工方案
2.1.1 圍堰抽水。DK2193+340~DK2193+510經(jīng)過(guò)兩個(gè)魚(yú)塘,經(jīng)過(guò)排水后回填土平整,以避免隧道頂?shù)孛娲嬖谶^(guò)多積水。原計(jì)劃在隧道掘進(jìn)魚(yú)塘范圍前,修筑土壩,抽干堰內(nèi)塘水,局部回填黏土,再在地表施工加固措施,但考慮現(xiàn)場(chǎng)工程條件復(fù)雜,僅進(jìn)行了部分圍堰抽水。
2.1.2 地面加固。地面采用袖閥管進(jìn)行注漿加固,注漿孔呈梅花形布置,排距1.5 m,注漿平面和立面位置如圖3和圖4所示。注漿底部進(jìn)入風(fēng)化巖(D3)或隔離層,采用純水泥漿液。含水量較大的地層采用水泥+水玻璃漿液,注漿壓力為0.5~1.0 MPa。注漿完成后,鉆孔取樣檢測(cè),發(fā)現(xiàn)其滲透系數(shù)和土的力學(xué)指標(biāo)有所改善。
2.2 施工工法
2.2.1 超前預(yù)報(bào)。施工采取掌子面地質(zhì)素描法、物探法(TSP、地質(zhì)雷達(dá))及鉆探法(超前深孔鉆探、炮眼孔加深鉆探)等綜合方法,進(jìn)行超前地質(zhì)預(yù)測(cè)預(yù)報(bào)。通過(guò)地質(zhì)描述法分析前方地質(zhì)構(gòu)造的變化,并通過(guò)物探法預(yù)報(bào)前方不良地質(zhì),以不斷提高預(yù)測(cè)預(yù)報(bào)的水平。
2.2.2 開(kāi)挖工法。隧道開(kāi)挖掘進(jìn)采用三臺(tái)階法臨時(shí)仰拱封閉法施工,其橫斷面如圖5所示。通過(guò)土層軟弱地段時(shí),上臺(tái)開(kāi)挖預(yù)留核心土,上臺(tái)開(kāi)挖支護(hù)加密超前小導(dǎo)管間距及搭接。小導(dǎo)管間距為35 cm,每進(jìn)尺1.5 m(即兩榀格柵鋼架距離)施工3.5 m小導(dǎo)管,搭接2 m。小導(dǎo)管注漿壓力為0.8~1.0 MPa。
2.2.3 監(jiān)控量測(cè)。隨時(shí)關(guān)注隧道穿越魚(yú)塘施工的監(jiān)測(cè)數(shù)據(jù),加強(qiáng)洞內(nèi)外觀察,并進(jìn)行掌子面地質(zhì)描述,對(duì)圍巖變形及支護(hù)應(yīng)力進(jìn)行監(jiān)測(cè),及時(shí)反饋檢測(cè)信息,確保工程安全。
3 有限元驗(yàn)證分析
3.1 數(shù)值模型
選用Midas GTS/NX軟件,基于巖土材料中的摩爾-庫(kù)倫本構(gòu)模型建立軟弱圍巖淺埋大斷面隧道開(kāi)挖下穿魚(yú)塘的數(shù)值計(jì)算模型。以DK2193+440斷面作為穿越水庫(kù)的數(shù)值分析斷面,以消除邊界效應(yīng)為原則,上表面取至地表面,下邊界及左右邊界均取5倍以上洞徑。地質(zhì)條件方面僅考慮粉質(zhì)黏土和粉砂巖兩種地層分布,粉質(zhì)黏土覆蓋層5.0 m左右,向下全部為粉砂巖,二維有限元模型如圖6所示。
數(shù)值分析計(jì)算的邊界條件按如下考慮:①應(yīng)力邊界,模型兩側(cè)和底部為法向約束,頂部為自由邊界,初始應(yīng)力場(chǎng)取自重應(yīng)力場(chǎng);②由于采用注漿加固方案,該模型不考慮滲流。
本研究計(jì)算時(shí)針對(duì)粉質(zhì)黏土和粉砂巖采用摩爾-庫(kù)倫模型進(jìn)行模擬,臨時(shí)仰拱、初期支護(hù)噴射混凝土和二襯模注混凝土按彈性計(jì)算,圍巖和結(jié)構(gòu)的物理力學(xué)參數(shù)見(jiàn)表1。
三臺(tái)階臨時(shí)仰拱封閉法施工工序?yàn)槌跏紙?chǎng)地應(yīng)力分析(工況1)—地表注漿(工況2)—上臺(tái)階開(kāi)挖(工況3)—上臺(tái)階初期支護(hù)(工況4)—中臺(tái)階開(kāi)挖(工況5)—中臺(tái)階初期支護(hù)(工況6)—下臺(tái)階開(kāi)挖(工況7)—下臺(tái)階支護(hù)(清底、施做仰拱和回填)(工況8)—澆筑二襯(工況9)。
3.2 計(jì)算結(jié)果
3.2.1 地層豎向位移。經(jīng)過(guò)計(jì)算,獲得施工各階段地層豎向位移值,如圖7所示,各施工階段下隧道中心地表沉降變化曲線如圖8所示。由圖可知,地層豎向位移變化趨勢(shì)如下。
①雖然DK2193+440斷面緊鄰魚(yú)塘,但在開(kāi)挖各階段,開(kāi)挖斷面附近圍巖地層豎向位移基本呈對(duì)稱(chēng)分布。由于拱頂覆土層較淺,從上臺(tái)階開(kāi)挖至下臺(tái)階過(guò)程中,隧道頂部位移下沉,底部位移隆起,地表沉降最大值發(fā)生在隧道拱頂偏右位置,這是因?yàn)樗淼理敳坑覀?cè)地形偏高,土體自重較大。地表隆起最大值發(fā)生在隧道仰拱中心位置,當(dāng)施工完二襯后,地表隆起最大值逐步向仰拱右側(cè)發(fā)展,這是因?yàn)槎r施工后,隧道中間仰拱回填導(dǎo)致隧道中心位置發(fā)生下沉位移,引起最大隆起位移右移。
②上臺(tái)階開(kāi)挖后,隧道中心處地表沉降值為6.029 mm,拱頂沉降值為8.622 mm,開(kāi)挖底部最大隆起值為5.046 mm;中臺(tái)階開(kāi)挖后,隧道中心處地表沉降值為7.957 mm,拱頂沉降值為11.364 mm,開(kāi)挖底部最大隆起值為12.698 mm;下臺(tái)階開(kāi)挖后,隧道中心處地表沉降值為10.185 mm,拱頂沉降值為14.278 mm,開(kāi)挖底部最大隆起值為23.546 mm;最終(二襯施工完畢后),隧道中心處地表沉降值為12.659 mm,拱頂沉降值為16.846 mm,開(kāi)挖底部最大隆起值為19.292 mm。由此可見(jiàn),隨著掌子面的掘進(jìn),地層不斷持續(xù)變形, 豎向位移值也逐漸變大,但隨著仰拱施工、回填及二襯澆筑完畢,開(kāi)挖邊界附近的圍巖變形基本穩(wěn)定。
3.2.2 隧道豎向位移。經(jīng)過(guò)計(jì)算,提取隧道拱頂位置和仰拱中心位置豎向變形曲線,如圖9所示。
從圖中可以看出:①由于拱頂覆土層較淺,從上臺(tái)階開(kāi)挖至下臺(tái)階過(guò)程中,隧道頂部位移下沉,底部位移隆起,隧道頂部沉降最大值發(fā)生在隧道拱頂偏右位置,這是因?yàn)樗淼理敳坑覀?cè)地勢(shì)偏高,土體自重較大。在仰拱施工前,隧道底部隆起最大值發(fā)生在隧道底部中心位置,當(dāng)仰拱施工、回填和二襯后,隧道底部隆起最大值逐步向仰拱右側(cè)發(fā)展,這是因?yàn)槎r施工時(shí),隧道中間仰拱回填荷載增大,導(dǎo)致隧道中心位置發(fā)生下沉位移,引起最大隆起位移右移。②上臺(tái)階開(kāi)挖完成初期支護(hù)后,隧道拱頂最大沉降值為9.305 mm,開(kāi)挖底部最大隆起值為10.988 mm;中臺(tái)階開(kāi)挖完成初期支護(hù)后,隧道拱頂最大沉降值為11.597 mm,開(kāi)挖底部最大隆起值為17.875 mm;下臺(tái)階開(kāi)挖完成仰拱及回填后,隧道拱頂最大沉降值為16.917 mm,開(kāi)挖底部最大隆起值為20.108 mm;最終(二襯施工完畢后),隧道拱頂最大沉降值為16.855 mm,開(kāi)挖底部最大隆起值為19.360 mm。由此可見(jiàn),隨著掌子面的掘進(jìn),隧道拱頂和底部地層不斷持續(xù)變形, 豎向位移值也逐漸變大。其中,拱頂?shù)貙幼畲蟪两底冃螢?4.299 mm,發(fā)生于下臺(tái)階支護(hù)階段(工況8),此時(shí)開(kāi)始施工隧道結(jié)構(gòu),土體開(kāi)挖變形穩(wěn)定,仰拱中心土體最大隆起值為23.546 mm,發(fā)生于下臺(tái)階開(kāi)挖(工況7),此時(shí)土體開(kāi)挖卸荷最嚴(yán)重,但是在仰拱施工與回填完畢后,開(kāi)挖邊界附近的圍巖變形基本穩(wěn)定。
4 現(xiàn)場(chǎng)監(jiān)測(cè)分析
4.1 監(jiān)測(cè)結(jié)果
在實(shí)際施工各階段中,對(duì)DK2193+440緊鄰魚(yú)塘斷面的拱頂沉降與隧道中心地表環(huán)境收斂進(jìn)行監(jiān)控量測(cè),其現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)如圖10所示。由圖10可知,①在各開(kāi)挖施工工序階段,斷面附近拱頂沉降和圍巖的收斂情況變化均較正常,該斷面測(cè)點(diǎn)的布置是在上臺(tái)階開(kāi)挖初期支護(hù)施工以后,此時(shí)拱頂沉降和地表周邊收斂值變化速率較小,開(kāi)挖中臺(tái)階時(shí),拱頂沉降速率及收斂速率變化顯著。結(jié)構(gòu)初期支護(hù)完成以后,圍巖與初期支護(hù)共同變形完成并逐漸趨于平穩(wěn),下臺(tái)階的開(kāi)挖對(duì)拱頂沉降和周邊收斂的影響不大,仰拱施作以后圍巖基本穩(wěn)定。②該斷面的拱頂沉降和地表周邊收斂最大值分別為10.895 mm與9.658 mm,小于安全限定值,意味著采用三臺(tái)階臨時(shí)仰拱封閉法淺埋穿越軟弱圍巖地層時(shí),隧道及周?chē)h(huán)境較穩(wěn)定。
4.2 對(duì)比分析
將實(shí)際施工各階段中DK2193+440緊鄰魚(yú)塘斷面的拱頂沉降與洞內(nèi)收斂的監(jiān)控量測(cè)數(shù)據(jù)與數(shù)值計(jì)算數(shù)據(jù)進(jìn)行對(duì)比,如圖11所示。由圖11可知,①對(duì)比結(jié)果反映的現(xiàn)場(chǎng)監(jiān)測(cè)變形規(guī)律與數(shù)值計(jì)算的變化趨勢(shì)基本相同,實(shí)際施工過(guò)程中,圍巖在開(kāi)挖以后一天左右時(shí)間內(nèi)出現(xiàn)較大變形,初期支護(hù)完成以后,圍巖與初期支護(hù)結(jié)構(gòu)共同變形完成并逐漸趨于平穩(wěn),這與數(shù)dfa9bfd137ecd0d0ea701fc7ab376028cafc10e790017a38feb8884f8fe8ab1a值計(jì)算結(jié)果一致,意味著數(shù)值計(jì)算結(jié)果合理可靠。②對(duì)于拱頂沉降,最大計(jì)算值為14.299 mm,最大監(jiān)測(cè)值為10.626 mm,誤差為3.673 mm;對(duì)于地表沉降,最大計(jì)算值為13.056 mm,最大監(jiān)測(cè)值為9.260 mm,誤差為3.769 mm。這說(shuō)明拱頂沉降和地表沉降均呈現(xiàn)計(jì)算值大于監(jiān)測(cè)值的趨勢(shì),這是因?yàn)橛?jì)算模型較理想化,且在模擬開(kāi)挖與初期支護(hù)的過(guò)程中,未對(duì)隧道結(jié)構(gòu)進(jìn)行約束,導(dǎo)致隧道開(kāi)挖面處于無(wú)約束狀態(tài),而實(shí)際現(xiàn)場(chǎng)隧道開(kāi)挖后會(huì)及時(shí)封閉。
5 結(jié)論
本研究分析了下穿魚(yú)塘軟弱圍巖大斷面隧道的施工技術(shù),詳細(xì)闡述了施工技術(shù)方案可行性,采用數(shù)值模擬的方法驗(yàn)證了該施工技術(shù)的合理性,主要結(jié)論如下。
①隧道DK2193+350~DK2193500段為淺埋段,巖體破碎,節(jié)理裂隙發(fā)育,拱頂易坍頂,施工時(shí)采用三臺(tái)階法與地表加固結(jié)合措施,以保障施工過(guò)程安全穩(wěn)定。
②隨著掌子面的掘進(jìn),地層不斷持續(xù)變形, 豎向位移值也逐漸變大,但在仰拱施工、回填及二襯澆筑完成后,開(kāi)挖邊界附近的圍巖變形基本穩(wěn)定。最終,隧道中心處地表沉降為12.659 mm,拱頂沉降為16.846 mm,開(kāi)挖底部最大隆起值為19.292 mm,滿足安全施工規(guī)范要求。
③隧道變形與地表沉降的數(shù)值計(jì)算與現(xiàn)場(chǎng)監(jiān)測(cè)的變化趨勢(shì)基本相同,拱頂沉降和地表周邊收斂最大監(jiān)測(cè)值分別為10.895 mm與9.658 mm,略小于數(shù)值計(jì)算值,這是由于實(shí)際施工會(huì)及時(shí)封閉隧道開(kāi)挖面,而數(shù)值計(jì)算則是偏理想化模型,但研究成果表明三臺(tái)階法臨時(shí)仰拱封閉法可用于軟弱圍巖大斷面隧道的開(kāi)挖。
參考文獻(xiàn):
[1]王小鎖. Ⅳ級(jí)圍巖超大斷面公路隧道開(kāi)挖方法與施工參數(shù)研究[D].北京: 北京交通大學(xué),2019.
[2]劉聰,李術(shù)才,周宗青,等.復(fù)雜地層超大斷面隧道施工圍巖力學(xué)特征模型試驗(yàn)[J].巖土力學(xué),2018,39(9):3495-3504.
[3]HISATAKE M,OHNO S.Effects of pipe roof supports and the excavation method on the displacements above a tunnel face[J].Tunnelling and Underground Space Technology,2008,23(2):120-127.
[4]HU X Y, CHENG J, JV J W. Influence of the cutterhead configuration and operation parameters on the face stability of epb shield tunnels in dry granular soils[J].International Journal of Geomechanics,2021,21(5): 04021050.
[5]程選生,王建華. 基于圍巖位移控制的超大斷面黃土隧道施工方法研究[J].巖土工程學(xué)報(bào),2013,35(S1):82-89.
[6]于天賜.軟巖隧道大變形控制技術(shù)研究[J].土木工程學(xué)報(bào),2017,50(S2):112-117.
[7]崔小鵬,孫韶峰,王廣宏,等. CRD工法及三臺(tái)階七步開(kāi)挖工法的對(duì)比和改進(jìn)[J].隧道建設(shè),2010,30(2):131-133,178.
[8]李波,魏龍海,舒恒.濱海軟弱地層超大斷面暗挖隧道施工工序研究[J].現(xiàn)代隧道技術(shù),2018,55(S2):283-294.
[9]石熊,張家生,劉寶琛.大斷面淺埋偏壓隧道CRD法施工工序研究[J].現(xiàn)代隧道技術(shù),2015,52(3):193-199.
[10]蔣坤,夏才初,卞躍威.節(jié)理巖體中雙向八車(chē)道小凈距隧道施工方案優(yōu)化分析[J].巖土力學(xué),2012,33(3):841-847.