陳鵬飛 張景科 諶文武 楊善龍 和法國(guó) 陸愷

內(nèi)容摘要：潮濕環(huán)境中考古探方的失穩(wěn)破壞不僅對(duì)遺址文物的保護(hù)不利，而且威脅著現(xiàn)場(chǎng)工作人員的人身安全。目前對(duì)潮濕環(huán)境下考古探方穩(wěn)定性預(yù)判的研究較少。本文基于杭州潮濕地區(qū)的巖土參數(shù)采用FLAC3D數(shù)值模擬軟件進(jìn)行探方開(kāi)挖前穩(wěn)定性預(yù)判模擬計(jì)算。模擬探方為直立式，尺寸（長(zhǎng)×寬×深）為5m×5m×2.4m。模擬發(fā)掘過(guò)程采取降排水措施。數(shù)值計(jì)算結(jié)果表明：（1）隨開(kāi)挖不斷進(jìn)行，坑壁水平位移逐漸增大，曲線形狀呈“弓”型，水平位移由頂部至底部先增大到極值后減小，極值點(diǎn)位于距頂部1.75m的坑壁處。（2）最終安全系數(shù)為2.37。（3）得出了潛在滑移帶所在位置，坑壁滑移帶剪出口位于坑壁2.0m—2.4m處。計(jì)算結(jié)果可以為現(xiàn)場(chǎng)考古模擬發(fā)掘試驗(yàn)提供參考，為直立式考古探方穩(wěn)定性預(yù)加固提供依據(jù)。

關(guān)鍵詞：潮濕環(huán)境；史前考古遺址；直立探方；FLAC3D；穩(wěn)定預(yù)判

中圖分類號(hào)：P901 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-4106（2016）04-0135-06

Abstract： The instability of an archaeology excavation unit in moist environments is not only bad for the sites and protection of relics， but also threatens the personal safety of archaeologists. Little research on the stability of this equipment when used in such environments has been conducted at present， however. This study presents research on how to prejudge the stability of excavation units prior to excavation by applying the numerical simulation software FLAC3D to the geological conditions and moist environment of Hangzhou when using a vertical excavation unit that measures 5m×5m×2.4m. The numerical calculation， combined with drainage measures also recorded during the simulation， yielded several results. First， during excavation the deformation of the pit wall increases gradually in the shape of a bow being drawn， and the horizontal displacement first increases to extremum and then decreases from the top to the bottom. The point where the extremum is reached is 1.75 m from the top of the pit wall. Second， the stability coefficient is 2.37， and； finally， the location of the potential slip band is clearly observable， the opening of the slip band being 2.0m-2.4m from the top of the pit wall. These calculation results can provide reference for field archaeological excavation and evidence for preparatory reinforcement of vertical excavation units.

Keywords： moist environment； prehistoric archaeology site； vertical excavation unit； FLAC3D； pre-judging stability；

1 引 言

我國(guó)東南方潮濕環(huán)境地區(qū)的史前考古遺址主要集中分布在長(zhǎng)江流域：如長(zhǎng)江上游地區(qū)的巫山大溪文化遺址、新繁水觀音遺址等，長(zhǎng)江中游地區(qū)的盤(pán)龍城遺址、屈家?guī)X文化遺址等，長(zhǎng)江下游地區(qū)的河姆渡文化遺址、良渚文化遺址等。為確保文化層信息的完整連續(xù)性，直立探方是常用的考古發(fā)掘型式。鑒于探方四壁均保存著重要文化信息，不宜過(guò)多采取工程加固措施；而文化層開(kāi)挖與信息的提取歷時(shí)長(zhǎng)，我國(guó)東南方潮濕環(huán)境中的直立探方極易失穩(wěn)破壞。探方失穩(wěn)不僅對(duì)發(fā)掘的古遺址保護(hù)有威脅，同時(shí)也會(huì)影響到考古人員的安全。因此，在潮濕環(huán)境下史前土遺址的考古探方穩(wěn)定性的預(yù)判尤為重要。

考古發(fā)掘探方與工程基坑有諸多類似之處。目前對(duì)于在潮濕環(huán)境中基坑工程的研究已經(jīng)取得了豐富的成果。隨著計(jì)算機(jī)技術(shù)的發(fā)展，基坑工程穩(wěn)定性預(yù)判多利用數(shù)值計(jì)算方法，有限元及有限差分?jǐn)?shù)值計(jì)算軟件在基坑模擬計(jì)算中的應(yīng)用越來(lái)越廣泛，而且取得了良好的效果。劉杰等通過(guò)ANSYS數(shù)值模擬計(jì)算對(duì)地鐵車站基坑的圍護(hù)結(jié)構(gòu)進(jìn)行了模擬，結(jié)合現(xiàn)場(chǎng)變形監(jiān)測(cè)得到了很好的驗(yàn)證[1]。李磊等則是通過(guò)ABAQUS有限元軟件對(duì)基坑進(jìn)行了數(shù)值模擬[2]。周念清等采用三維有限差分法對(duì)基坑降水進(jìn)行了模擬計(jì)算，模擬結(jié)果與地下水位實(shí)際監(jiān)測(cè)數(shù)據(jù)十分吻合[3]?？椎律韧ㄟ^(guò)建立二維有限元彈塑性分析模型采用強(qiáng)度折減法對(duì)軟土深基坑坑底抗隆起情況進(jìn)行了穩(wěn)定性分析[4]。吳忠誠(chéng)等通過(guò)深基坑復(fù)合土釘墻支護(hù)大型現(xiàn)場(chǎng)原位測(cè)試對(duì)比FLAC3D模擬，得出通過(guò)建立的三維數(shù)值模擬分析模型能夠反映復(fù)合土釘支護(hù)的基坑在開(kāi)挖過(guò)程中的受力和變形特征[5]。蔡海波等通過(guò)對(duì)基坑在原有設(shè)計(jì)基礎(chǔ)上沿深開(kāi)挖深度的工況下進(jìn)行模擬計(jì)算，優(yōu)化了原有設(shè)計(jì)的支護(hù)方案[6]。吳意謙等對(duì)濕陷性黃土地區(qū)深基坑進(jìn)行了模擬計(jì)算，與現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果較為一致，為基坑的合理設(shè)計(jì)與安全施工提供了科學(xué)依據(jù)[7]。丁勇春等通過(guò)計(jì)算模擬出臨江基坑在考慮江內(nèi)潮位變化情況下基坑支護(hù)結(jié)構(gòu)的穩(wěn)定性[8]?？梢?jiàn)數(shù)值模擬方法在基坑工程優(yōu)化設(shè)計(jì)、穩(wěn)定性預(yù)判中應(yīng)用廣泛，通過(guò)現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)驗(yàn)證了計(jì)算結(jié)果的合理性。結(jié)合考古發(fā)掘坑與基坑的類同性，選擇性地引入基坑穩(wěn)定性預(yù)判的理論方法應(yīng)用到考古發(fā)掘探方的穩(wěn)定性預(yù)判之中，對(duì)于考古發(fā)掘具有重要意義。

FLAC3D軟件能較好地模擬地質(zhì)材料在達(dá)到強(qiáng)度極限或屈服極限時(shí)發(fā)生的破壞或塑性流動(dòng)的力學(xué)行為，特別適用于分析漸進(jìn)破壞和失穩(wěn)以及模擬大變形，具有強(qiáng)大的前后處理功能，在巖土工程問(wèn)題中應(yīng)用廣泛[9-12]。本文采用FLAC3D軟件。選取浙江大學(xué)紫金港校區(qū)巖土條件為基本工況，模擬了開(kāi)挖深2.4m、長(zhǎng)寬均為5.0m的模擬考古直立探方，開(kāi)展了FLAC3D計(jì)算軟件在考古探方穩(wěn)定性預(yù)判中的探索。

2 試驗(yàn)場(chǎng)地概況

2.1 試驗(yàn)地點(diǎn)

現(xiàn)場(chǎng)試驗(yàn)選擇在浙江省科技考古與文物保護(hù)技術(shù)研究試驗(yàn)基地。該場(chǎng)地位于浙江省杭州市浙江大學(xué)紫金港校區(qū)。

2.2 試驗(yàn)場(chǎng)地與史前考古遺址的相似性

通過(guò)現(xiàn)場(chǎng)調(diào)查，場(chǎng)地土體為物理力學(xué)性質(zhì)差的黏土、淤泥質(zhì)粉質(zhì)黏土；場(chǎng)地地下水位為0.8m，且試驗(yàn)基地周邊有多處池塘、水田。該場(chǎng)地的土體條件、地下水情況與良渚遺址的考古環(huán)境相似。

2.3 巖土條件及地下水

試驗(yàn)基地現(xiàn)場(chǎng)勘察表明，試驗(yàn)場(chǎng)地地層自上而下依次是：

素填土：黃褐色內(nèi)部有較多裂隙，裂隙寬1—2cm，土體厚度80cm。

雜填土：主要由磚塊、水泥塊及夾雜其中的碎土組成，力學(xué)性質(zhì)較差、透水性好。土體厚度35cm。

素填土：深褐色，土體內(nèi)部有水平發(fā)育空隙，空隙內(nèi)有細(xì)沙充填。土體厚度35cm。

淤泥質(zhì)黏土：灰黑色，土質(zhì)均勻，流塑性，土體性質(zhì)較差。土體厚度20cm。

黏土：黃色，土質(zhì)均勻，軟塑性，小孔隙發(fā)育，土體疏松，性質(zhì)較差。土體厚度70cm。

黏土：青灰色，土質(zhì)均勻，可塑性，無(wú)裂隙發(fā)育，土體性質(zhì)較差。

各土層的物理力學(xué)參數(shù)如表1所示。

場(chǎng)地的地下水豐富，周邊有多處池塘、水田，現(xiàn)場(chǎng)地下水位一般位置在0.8m附近，隨天氣變化而波動(dòng)，波動(dòng)幅度0.1m。

3 模擬開(kāi)挖方案及方法

3.1 模擬工況及方法

試驗(yàn)現(xiàn)場(chǎng)為典型潮濕環(huán)境，池塘、河流、湖泊等地表水豐富，降水量大，導(dǎo)致地下水位高且補(bǔ)給源充足，為使模擬現(xiàn)場(chǎng)考古發(fā)掘作業(yè)正常進(jìn)行，在模擬探方開(kāi)挖過(guò)程中以及后期文化層信息提取收集過(guò)程中必須全程實(shí)施降排水措施，將地下水降至考古探方底部以下，減少地下水的不良影響。

由于實(shí)驗(yàn)現(xiàn)場(chǎng)土層多為滲透性極差的黏土組成，地下水滲透緩慢，滲透壓力對(duì)于模擬探方穩(wěn)定性影響相對(duì)較弱，同時(shí)為簡(jiǎn)化計(jì)算，本次模擬計(jì)算僅考慮地下水面以下的土體采用飽和狀態(tài)下的物理力學(xué)參數(shù)，地下水位面以上的土體采用天然狀態(tài)下的物理力學(xué)參數(shù)進(jìn)行模擬計(jì)算。

3.2 開(kāi)挖、邊界、網(wǎng)格

模擬開(kāi)挖可認(rèn)為是平面應(yīng)變問(wèn)題而且開(kāi)挖探方模型對(duì)稱，可以將模型簡(jiǎn)化為二維問(wèn)題并選取一半作為計(jì)算對(duì)象。模型的高度為12m，長(zhǎng)15m，在FLAC3D中對(duì)模型按現(xiàn)場(chǎng)土層情況進(jìn)行分組，模型由6360個(gè)單元（zones）和9882個(gè)節(jié)點(diǎn)（grid-

points）組成。

數(shù)值模擬計(jì)算的土體本構(gòu)模型采用Mohr-

Cloumb模型，除開(kāi)挖探方所形成的邊坡及坑底設(shè)為自由邊界外，模型底部為固定約束邊界，限制其在法向方向以及切向方向的位移，模型的四周方向?yàn)閱蜗蜻吔?，限制其在法向方向的位移?/p>

3.3 模擬開(kāi)挖方案

基于史前考古土遺址發(fā)掘?qū)嵺`及FLAC3D模擬軟件的可實(shí)施性，模擬開(kāi)挖方案為：開(kāi)挖分為四步進(jìn)行，每次開(kāi)挖0.6m，最終開(kāi)挖深度2.4m，開(kāi)挖情況如圖1所示。

4 模擬分析結(jié)果及穩(wěn)定性預(yù)判

4.1 位移云圖

位移監(jiān)測(cè)可以直接反應(yīng)出數(shù)值計(jì)算中模型監(jiān)測(cè)位置的變形情況，也可推演出現(xiàn)實(shí)情況下對(duì)應(yīng)位置的變形情況，結(jié)合現(xiàn)場(chǎng)地質(zhì)條件的認(rèn)識(shí)便可判別出對(duì)象的穩(wěn)定性以及潛在破壞的位置，而且在生產(chǎn)實(shí)踐中位移的監(jiān)測(cè)較應(yīng)力的監(jiān)測(cè)更加簡(jiǎn)便、精確。

由水平位移云圖2可知，水平橫向位移主要集中在坑壁附近，隨開(kāi)挖不斷進(jìn)行變形范圍逐漸增加，且位移量也逐漸增大；坡腳處也出現(xiàn)水平位移，且隨著開(kāi)挖進(jìn)行，出現(xiàn)水平位移的范圍逐漸增大。

圖3為每次開(kāi)挖結(jié)束后坑壁最終達(dá)到穩(wěn)定狀態(tài)時(shí)位移情況，Y軸為基坑深度，X軸為坑壁的水平位移值，負(fù)值表示位移方向向探方內(nèi)部位移，反之正方向表示向探方外側(cè)位移。由監(jiān)測(cè)結(jié)果可知，第一次開(kāi)挖，開(kāi)挖深度為0.6m相對(duì)較淺，水平位移量較小，在1mm左右；第二次開(kāi)挖0.6m，模擬探方深度為1.2m，由坡頂至坡底水平位移量逐漸增大，1.0m處水平位移量達(dá)到極值；第三次開(kāi)挖0.6m，模擬探方深度為1.8m，由坡頂至坡底水平位移量先增大后減小，在1.5m處時(shí)水平位移量達(dá)到極值1mm；第四次開(kāi)挖0.6m，模擬探方深度為2.4m，由坡頂至坡底水平位移量先增大后減小，在1.75m處時(shí)水平位移量達(dá)到極值3.9mm。

4.2 安全系數(shù)

FLAC3D中對(duì)于邊坡穩(wěn)定性的數(shù)值模擬常用方法為強(qiáng)度折減法，即通過(guò)降低巖土體的參數(shù)，使得邊坡剛好達(dá)到臨界破壞狀態(tài)，此時(shí)對(duì)于巖土體的抗剪強(qiáng)度進(jìn)行的折減程度即定義安全系數(shù)。安全系數(shù)為巖土體的實(shí)際抗剪強(qiáng)度與臨界破壞時(shí)的折減后剪切強(qiáng)度的比值。在數(shù)值計(jì)算中分別采用強(qiáng)度折減法計(jì)算每次開(kāi)挖后探方模型的安全系數(shù)，四次開(kāi)挖完成后的安全系數(shù)如圖4所示。

由圖可知隨開(kāi)挖深度不斷增加坑壁的安全系數(shù)逐漸降低，四次開(kāi)挖完成后的安全系數(shù)分別為14.71、4.86、2.93、2.37，第二次開(kāi)挖完成后，安全系數(shù)降低幅度最大，后兩次開(kāi)挖完成后安全系數(shù)變化值相對(duì)較小。

4.3 剪應(yīng)變?cè)隽?/p>

模擬開(kāi)挖完成后探方外土體的剪應(yīng)變?cè)隽吭茍D及位移矢量云圖如圖5所示。圖中箭頭方向代表對(duì)應(yīng)單元的位移方向，箭頭長(zhǎng)度代表位移值的相對(duì)大小，由圖可知在模型開(kāi)挖完成后，坑壁有向坑內(nèi)滑移的趨勢(shì)，坑壁處位移量最大，其后面土體的位移量逐漸降低。

由剪應(yīng)變?cè)隽吭茍D可知剪應(yīng)變?cè)隽糠植汲省皫А睜罘植肌Ｍ馏w的破壞主要有兩種：拉伸破壞與剪切破壞[13]。本次模擬中拉伸塑性區(qū)主要位于模擬探方外側(cè)后緣1.5m—2.0m處，沿垂向發(fā)育深度約為1.5m，此區(qū)域土體受水平向拉應(yīng)力作用，當(dāng)拉應(yīng)力超出土體的極限抗拉強(qiáng)度，此區(qū)域便會(huì)出現(xiàn)后緣拉裂縫。剪切塑性區(qū)主要位于坑壁2.0m—2.4m處，與坑壁呈45°夾角，在距頂部2.3m的坑壁處剪應(yīng)變?cè)隽孔畲蟆?/p>

5 分析與討論

5.1 位移云圖

隨著開(kāi)挖不斷進(jìn)行，變形量逐漸增大，坑壁的變形曲線逐漸呈“弓”型，水平位移由頂部至底部先增大到極值后減小，極值點(diǎn)附近即為此模擬探方的危險(xiǎn)帶。危險(xiǎn)帶位于距頂部1.75m左右，此處地層為淤泥質(zhì)粉質(zhì)黏土與黃色黏土，土體物理力學(xué)性質(zhì)極差，所以在開(kāi)挖完成后，由于開(kāi)挖卸荷出現(xiàn)臨空面。坑壁下部土體作為主動(dòng)區(qū)，失去了被開(kāi)挖掉的探方內(nèi)被動(dòng)區(qū)土體的支撐，所以此處將會(huì)產(chǎn)生較大變形。此處將是重點(diǎn)監(jiān)測(cè)、治理地帶。最終開(kāi)挖完成后的變形量最大值僅為3.9mm，約為直立探方開(kāi)挖深度的1.6‰，變形值相對(duì)較小，不會(huì)對(duì)考古探方內(nèi)的文化層信息造成破壞。

5.2 安全系數(shù)

第一次開(kāi)挖完成后，開(kāi)挖深度為0.6m模型安全系數(shù)14.71。第二次開(kāi)挖完成后，模型深度達(dá)1.2m，模型的安全系數(shù)降低至4.86，其穩(wěn)定性大幅降低。后面兩次開(kāi)挖完成后模型的安全系數(shù)變化量小，但安全系數(shù)相對(duì)較小。

最終模型的安全系數(shù)為2.37，安全儲(chǔ)備足。在降排水的條件下除發(fā)生極端事件，使初始條件發(fā)生明顯變化外，考古探方是可以達(dá)到穩(wěn)定狀態(tài)的，即在本次模擬的工況下進(jìn)行現(xiàn)場(chǎng)開(kāi)挖是不會(huì)出現(xiàn)失穩(wěn)破壞情況的。

5.3 剪應(yīng)變?cè)隽?/p>

巖土體的破壞是由于某一面上的剪應(yīng)力達(dá)到了巖土體的剪切強(qiáng)度，此時(shí)剪切面上必然發(fā)生了較大的剪切變形。通過(guò)剪應(yīng)變?cè)隽吭茍D的情況即可判斷出潛在滑移面的位置[14]，通過(guò)剪應(yīng)變?cè)隽康呢炌ㄒ部勺鳛檫吰伦冃纹茐牡臏?zhǔn)則[15]。此剪應(yīng)變?cè)隽吭茍D是在對(duì)土體參數(shù)進(jìn)行不斷折減的情況下獲取的，由此可以預(yù)知由于土體的力學(xué)性質(zhì)降低而引起的滑移破壞情況。

可看出潛在滑移帶所在位置，在開(kāi)挖初期由于開(kāi)挖形成臨空面，改變了初始穩(wěn)定狀態(tài)，臨空面出現(xiàn)向探方內(nèi)部位移趨勢(shì)，后部開(kāi)始出現(xiàn)拉張帶，約束臨空面的位移。隨開(kāi)挖深度增加，臨空面增大，拉張帶垂向延伸發(fā)育，坑壁底部出現(xiàn)剪切帶?？梢钥闯鰞烧呔诟髯韵鄳?yīng)方向發(fā)育，但兩者并未貫通，由此可知在浙江大學(xué)紫金港校區(qū)內(nèi)進(jìn)行的現(xiàn)場(chǎng)模擬史前考古發(fā)掘直立探方，在降排水措施下是不會(huì)出現(xiàn)失穩(wěn)破壞的；模擬探方潛在剪出口位于坑壁2.0m—2.4m處，此區(qū)域?yàn)橐灼茐膮^(qū)域，應(yīng)重點(diǎn)監(jiān)測(cè)其位移變化情況。

6 結(jié) 論

通過(guò)FLAC3D數(shù)值模擬軟件對(duì)浙江省杭州市浙江大學(xué)紫金港校區(qū)內(nèi)尺寸（長(zhǎng)×寬×深）為5m×5m×2.4m的直立式模擬考古探方的計(jì)算得出如下結(jié)論：

（1）隨著開(kāi)挖進(jìn)行，坑壁變形量逐漸增大，曲線形狀呈“弓”型，水平位移由頂部至底部先增大到極值后減小，危險(xiǎn)帶位于距頂部1.75m附近的坑壁處。

（2）在降排水的條件下，模擬探方開(kāi)挖完成后最終穩(wěn)定系數(shù)為2.37。

（3）得出了潛在滑移帶所在位置，坑壁剪切塑性區(qū)剪出口位于坑壁2.0m—2.4m處，此處為易破壞區(qū)域。

（4）據(jù)計(jì)算結(jié)果為現(xiàn)場(chǎng)考古模擬發(fā)掘試驗(yàn)提供參考，為直立式考古探方穩(wěn)定性預(yù)加固提供依據(jù)?？梢?jiàn)，通過(guò)采用數(shù)值模擬軟件FLAC3D，可以為我國(guó)東南方潮濕地區(qū)考古發(fā)掘探方穩(wěn)定性預(yù)判及預(yù)加固提供可靠依據(jù)。

參考文獻(xiàn)：

[1]劉杰，姚海林，任建喜.地鐵車站基坑圍護(hù)結(jié)構(gòu)變形監(jiān)測(cè)與數(shù)值模擬[J].巖土力學(xué)，2010，31（增2）：456-461.

[2]李磊，段寶福.地鐵車站深基坑工程的監(jiān)控量測(cè)與數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報(bào)，2013，32（增1）：2684-2691.

[3]周念清，唐益群，等.徐家匯地鐵站深基坑降水?dāng)?shù)值模擬與沉降控制[J].巖土工程學(xué)報(bào)，2011，33（12）：1950-1956.

[4]孔德森，門(mén)燕青，等.軟土深基坑坑底抗隆起穩(wěn)定性數(shù)值分析[J].巖土工程學(xué)報(bào)，2010，32（增1）：186-190.

[5]吳忠誠(chéng)，湯連生，廖志強(qiáng)，等.深基坑復(fù)合土釘墻支護(hù)FLAC-3D模擬及大型現(xiàn)場(chǎng)原位測(cè)試研究[J].巖土工程學(xué)報(bào)，2006，28（增）：1460-1465.

[6]蔡海波，吳順川，等.既有基坑延深開(kāi)挖穩(wěn)定性評(píng)價(jià)與支護(hù)方案確定[J].巖土力學(xué)，2011，32（11）：3306-3312.

[7]吳意謙，朱彥鵬.蘭州市濕陷性黃土地區(qū)地鐵車站深基坑變形規(guī)律監(jiān)測(cè)與數(shù)值模擬研究[J].巖土工程學(xué)報(bào)，2014，36（增2）：404-411.

[8]丁勇春，程澤坤，等.臨江基坑變形及受力性狀三維數(shù)值分析[J].巖土工程學(xué)報(bào)，2012，34（增）：243-247.

[9]胡安峰，張光建，等.地鐵換乘車站基坑圍護(hù)結(jié)構(gòu)變形監(jiān)測(cè)與數(shù)值模擬[J].巖土工程學(xué)報(bào)，2012，34（增）：77-81.

[10]寇曉東，周維垣，等.FLAC-3D進(jìn)行三峽船閘高邊坡穩(wěn)定分析[J].巖石力學(xué)與工程學(xué)報(bào)，2001，20（1）：610.

[11]于丹，郭舉興，等.某停車場(chǎng)深基坑支護(hù)方式的FLAC3D模擬分析[J].沈陽(yáng)建筑大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，31（1）：38-44.

[12]單仁亮，董洪國(guó)，等.中深部含軟弱夾層的深基坑土釘支護(hù)失穩(wěn)破壞數(shù)值模擬分析[J].巖土工程學(xué)報(bào)，2014，36（增2）：30-35.

[13]丁勇春，周順新，等.深基坑開(kāi)挖與土釘支護(hù)三維數(shù)值分析[J].上海交通大學(xué)學(xué)報(bào)，2011，45（4）：547-552.

[14]李劍，陳善雄，等.基于最大剪應(yīng)變?cè)隽康倪吰聺撛诨瑒?dòng)面搜索[J].巖土力學(xué)，2013，34（增1）：371-378.

[15]連鎮(zhèn)營(yíng)，韓國(guó)城，等.強(qiáng)度折減有限元法研究開(kāi)挖邊坡的穩(wěn)定性[J].巖土工程學(xué)報(bào)，2001，23（4）：407-411.