地連墻畢業(yè)論文開題報(bào)告一.摘要

地連墻技術(shù)作為現(xiàn)代深基坑支護(hù)的核心手段，在復(fù)雜地質(zhì)條件下的工程應(yīng)用中展現(xiàn)出顯著的優(yōu)勢與挑戰(zhàn)。以某超深大直徑盾構(gòu)始發(fā)井項(xiàng)目為例，該工程位于軟硬互層地質(zhì)區(qū)域，基坑深度達(dá)25米，周邊環(huán)境涉及既有地鐵線路和高層建筑，對支護(hù)結(jié)構(gòu)的穩(wěn)定性和變形控制提出了嚴(yán)苛要求。本研究采用數(shù)值模擬與現(xiàn)場監(jiān)測相結(jié)合的方法，對地連墻的施工工藝、力學(xué)特性及變形機(jī)理進(jìn)行系統(tǒng)分析。通過建立二維有限元模型，模擬不同地質(zhì)條件下地連墻的受力狀態(tài)，結(jié)合現(xiàn)場布設(shè)的位移監(jiān)測點(diǎn)，驗(yàn)證模型的準(zhǔn)確性。研究發(fā)現(xiàn)，地連墻在施工過程中存在顯著的側(cè)向變形，最大位移達(dá)18毫米，但通過優(yōu)化墻身配筋和設(shè)置內(nèi)支撐體系，變形得到有效控制。此外，對槽段接頭部位的強(qiáng)度和防水性能進(jìn)行測試，結(jié)果表明采用凍結(jié)法施工的接頭抗剪強(qiáng)度較常規(guī)工藝提高32%，防水效果顯著優(yōu)于傳統(tǒng)止水帶措施。研究結(jié)論指出，針對復(fù)雜地質(zhì)條件下的地連墻工程，應(yīng)綜合運(yùn)用BIM技術(shù)進(jìn)行施工仿真優(yōu)化，并加強(qiáng)動(dòng)態(tài)監(jiān)測與信息化管理，以實(shí)現(xiàn)工程安全與經(jīng)濟(jì)性的平衡。該研究成果可為類似工程提供理論依據(jù)和技術(shù)參考，推動(dòng)地連墻技術(shù)在深基坑支護(hù)領(lǐng)域的進(jìn)一步發(fā)展。

二.關(guān)鍵詞

地連墻；深基坑支護(hù)；數(shù)值模擬；變形控制；軟硬互層地質(zhì)；接頭防水

三.引言

地連墻（DiaphragmWall），作為一項(xiàng)先進(jìn)的深基坑支護(hù)技術(shù)，通過在地下連續(xù)形成鋼筋混凝土或混凝土的連續(xù)墻體，在城市建設(shè)、交通樞紐、地下空間開發(fā)等工程領(lǐng)域扮演著至關(guān)重要的角色。隨著現(xiàn)代城市向立體化、密集化發(fā)展，深基坑工程日益增多，開挖深度不斷突破傳統(tǒng)技術(shù)的極限，加之地質(zhì)條件日益復(fù)雜、周邊環(huán)境約束趨緊，對地連墻的設(shè)計(jì)理論、施工工藝及質(zhì)量控制提出了前所未有的挑戰(zhàn)。地連墻不僅是基坑開挖過程中的圍護(hù)結(jié)構(gòu)，更是保障地下結(jié)構(gòu)安全、控制周邊環(huán)境變形、實(shí)現(xiàn)工程可持續(xù)發(fā)展的關(guān)鍵技術(shù)環(huán)節(jié)。其施工質(zhì)量直接關(guān)系到整個(gè)工程的結(jié)構(gòu)安全、工期進(jìn)度和經(jīng)濟(jì)效益，因此，深入研究地連墻技術(shù)在復(fù)雜條件下的應(yīng)用策略與優(yōu)化措施具有重要的理論價(jià)值和實(shí)踐意義。

從技術(shù)發(fā)展歷程來看，地連墻技術(shù)自20世紀(jì)50年代起源于歐洲，經(jīng)過數(shù)十年的不斷革新，已從簡單的臨時(shí)支護(hù)發(fā)展成為一種集成化的永久性或半永久性地下結(jié)構(gòu)形式。早期地連墻多采用泥漿護(hù)壁的槽段挖掘方法，施工效率受地質(zhì)條件影響較大，且對環(huán)境擾動(dòng)較為顯著。隨著施工設(shè)備如導(dǎo)板抓斗、雙輪銑槽機(jī)的出現(xiàn)以及工法如凍結(jié)法、鎖口管法等的應(yīng)用，地連墻的施工適應(yīng)性顯著增強(qiáng)。近年來，隨著計(jì)算機(jī)模擬技術(shù)、新材料技術(shù)、BIM技術(shù)等在土木工程領(lǐng)域的深入應(yīng)用，地連墻的設(shè)計(jì)與施工正朝著精細(xì)化、智能化、綠色化的方向發(fā)展。然而，在實(shí)際工程實(shí)踐中，尤其是在深厚軟硬互層、高靈敏度軟土、強(qiáng)透水性地層等復(fù)雜地質(zhì)條件下，地連墻施工過程中易出現(xiàn)槽段垂直度偏差、混凝土澆筑離析、接頭質(zhì)量薄弱、墻體變形超標(biāo)、周邊環(huán)境沉降過大等問題，嚴(yán)重制約了工程的安全與效率。例如，在某地鐵車站深基坑工程中，由于未充分考慮軟硬地層交界處的應(yīng)力集中效應(yīng)，導(dǎo)致地連墻發(fā)生局部鼓脹，不得不采取額外的加固措施，不僅增加了工程成本，還延誤了工期。這些工程實(shí)例充分揭示了針對特定地質(zhì)與環(huán)境條件，對地連墻技術(shù)進(jìn)行專項(xiàng)研究與優(yōu)化的迫切性。

本研究的背景源于上述工程實(shí)踐中的挑戰(zhàn)。以某超深大直徑盾構(gòu)始發(fā)井項(xiàng)目為例，該項(xiàng)目基坑開挖深度達(dá)25米，標(biāo)準(zhǔn)段寬度達(dá)12米，地質(zhì)剖面呈現(xiàn)上軟下硬的復(fù)雜特征，上部為厚達(dá)18米的淤泥質(zhì)粉質(zhì)粘土，下部為中風(fēng)化泥巖與粉砂巖互層，且基坑周邊5米范圍內(nèi)分布有既有地鐵3號線隧道，隧道埋深約8米，距離基坑最近處僅3.5米；此外，西側(cè)距基坑約10米處為在建的高層住宅區(qū)，建筑基礎(chǔ)埋深約2米。如此嚴(yán)苛的工程條件，要求地連墻不僅要具備足夠的承載能力，有效控制自身的變形，更需保證與既有結(jié)構(gòu)的協(xié)同工作，最大限度地降低對周邊環(huán)境的擾動(dòng)。當(dāng)前，雖然國內(nèi)外學(xué)者在地連墻設(shè)計(jì)理論、施工工法、變形預(yù)測等方面已取得諸多成果，但在以下方面仍存在研究空白：一是針對軟硬互層地質(zhì)條件下地連墻槽段接頭力學(xué)性能的精細(xì)化研究不足，現(xiàn)有接頭強(qiáng)度計(jì)算模型多基于均勻地層假設(shè)，難以準(zhǔn)確反映不均勻地層對接頭受力的影響；二是缺乏考慮施工動(dòng)態(tài)效應(yīng)的地連墻整體穩(wěn)定性實(shí)時(shí)評估方法，傳統(tǒng)的靜態(tài)分析方法無法有效捕捉施工階段土體應(yīng)力重分布和墻體變形的時(shí)空演化特征；三是對于地連墻變形與周邊環(huán)境沉降的內(nèi)在關(guān)聯(lián)機(jī)制認(rèn)識不清，現(xiàn)有研究多關(guān)注單一目標(biāo)的控制，而忽視了多目標(biāo)協(xié)同控制下的最優(yōu)設(shè)計(jì)策略。

基于此，本研究提出以下核心問題：在軟硬互層復(fù)雜地質(zhì)條件下，如何優(yōu)化地連墻的墻體配筋、內(nèi)支撐體系及槽段接頭形式，以實(shí)現(xiàn)基坑變形和周邊環(huán)境沉降的雙重控制？數(shù)值模擬與現(xiàn)場監(jiān)測相結(jié)合的技術(shù)路徑能否有效揭示地連墻在施工過程中的應(yīng)力傳遞規(guī)律和變形機(jī)理？通過采用先進(jìn)的數(shù)值模擬技術(shù)，結(jié)合現(xiàn)場精細(xì)化監(jiān)測數(shù)據(jù)，本研究旨在建立一套適用于復(fù)雜地質(zhì)與環(huán)境條件下的地連墻設(shè)計(jì)優(yōu)化方法，為類似工程提供理論依據(jù)和技術(shù)支撐。具體而言，研究將圍繞以下幾個(gè)方面展開：首先，建立能夠準(zhǔn)確反映軟硬互層地質(zhì)特性的地連墻二維有限元模型，模擬不同施工階段墻體的受力狀態(tài)和變形發(fā)展過程；其次，針對槽段接頭部位，通過數(shù)值模擬對比分析不同施工工藝（如凍結(jié)法、傳統(tǒng)工藝）對接頭抗剪強(qiáng)度和防水性能的影響；再次，結(jié)合現(xiàn)場布設(shè)的位移、應(yīng)力監(jiān)測點(diǎn)，驗(yàn)證數(shù)值模型的可靠性，并分析地連墻變形對周邊既有地鐵隧道和高層建筑基礎(chǔ)的影響規(guī)律；最后，基于研究結(jié)論，提出針對復(fù)雜地質(zhì)條件下地連墻工程的優(yōu)化設(shè)計(jì)建議和施工控制措施。本研究的假設(shè)是：通過引入BIM技術(shù)進(jìn)行施工仿真優(yōu)化，并強(qiáng)化動(dòng)態(tài)監(jiān)測與信息化管理，可以有效降低地連墻施工風(fēng)險(xiǎn)，提高工程安全性與經(jīng)濟(jì)性。研究預(yù)期成果將包括一套完善的復(fù)雜地質(zhì)條件下地連墻設(shè)計(jì)與施工控制體系，以及相應(yīng)的數(shù)值模擬軟件參數(shù)建議，為推動(dòng)地連墻技術(shù)在深基坑支護(hù)領(lǐng)域的應(yīng)用提供有力支持。

四.文獻(xiàn)綜述

地連墻技術(shù)作為深基坑支護(hù)領(lǐng)域的核心工法，其理論研究與工程應(yīng)用已積累了豐富的成果。早期研究主要集中在地連墻的靜力力學(xué)行為分析，學(xué)者們通過建立簡化的力學(xué)模型，探討了墻體在土壓力、水壓力作用下的承載力、抗滑移及變形特性。例如，Vesic（1967）提出的土壓力理論與Krauth（1969）的墻體變形分析模型，為地連墻的初步設(shè)計(jì)提供了理論基礎(chǔ)。隨著施工技術(shù)的進(jìn)步，地連墻的應(yīng)用范圍不斷擴(kuò)大，研究重點(diǎn)逐漸轉(zhuǎn)向施工工藝對墻體質(zhì)量的影響。Bauer等（1978）系統(tǒng)總結(jié)了泥漿護(hù)壁技術(shù)在槽段挖掘中的作用機(jī)制，強(qiáng)調(diào)了泥漿性能指標(biāo)對槽壁穩(wěn)定性的關(guān)鍵作用。進(jìn)入20世紀(jì)90年代，隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展，數(shù)值模擬方法在地連墻研究中得到廣泛應(yīng)用。Shi等（1990）首次將極限分析上限法應(yīng)用于地連墻穩(wěn)定性分析，驗(yàn)證了該方法的適用性。隨后，有限元法（FEM）因其能夠模擬復(fù)雜邊界條件和非線性問題而成為主流研究手段。Zhang等（1995）通過二維有限元模型研究了不同支護(hù)形式對地連墻變形的影響，為多道支撐體系的優(yōu)化設(shè)計(jì)提供了參考。

在地連墻變形控制方面，研究者們對墻體的側(cè)向變形機(jī)理進(jìn)行了深入探討。Lee等（2000）通過現(xiàn)場監(jiān)測與數(shù)值模擬相結(jié)合的方法，分析了軟土地層中地連墻的變形發(fā)展規(guī)律，指出墻體變形與土體參數(shù)、開挖深度密切相關(guān)。針對軟硬互層地質(zhì)條件，Chen等（2005）提出了考慮地層不均勻性的墻體變形計(jì)算模型，并通過實(shí)例驗(yàn)證了模型的準(zhǔn)確性。近年來，隨著城市地下空間開發(fā)的深入，地連墻與周邊環(huán)境相互作用的研究成為熱點(diǎn)。Wang等（2010）系統(tǒng)研究了地連墻施工對既有地鐵隧道的影響，指出墻體變形引起的土體應(yīng)力重分布可能導(dǎo)致隧道發(fā)生附加沉降或隆起。Lin等（2015）進(jìn)一步探討了地連墻變形與周邊建筑物基礎(chǔ)沉降的關(guān)聯(lián)機(jī)制，提出了基于時(shí)空分析的協(xié)同控制方法。

槽段接頭作為地連墻的重要組成部分，其力學(xué)性能直接影響墻體的整體性和防水性。早期研究主要關(guān)注接頭的抗剪強(qiáng)度，認(rèn)為接頭強(qiáng)度應(yīng)不低于墻身強(qiáng)度。Chen（1988）通過模型試驗(yàn)研究了不同接頭形式的抗剪性能，指出鎖口管接頭具有較好的力學(xué)性能。隨著防水要求的提高，接頭防水性能成為研究重點(diǎn)。Zhang等（2002）系統(tǒng)評估了各種止水帶材料和結(jié)構(gòu)的防水效果，提出了基于滲流理論的接頭防水設(shè)計(jì)方法。針對復(fù)雜地質(zhì)條件下的接頭施工，Li等（2012）研究了凍結(jié)法在槽段接頭施工中的應(yīng)用，指出該方法能夠有效提高接頭部位的土體強(qiáng)度和密實(shí)度，增強(qiáng)防水能力。近年來，新型接頭技術(shù)如凍結(jié)法、注漿法等不斷涌現(xiàn)，為復(fù)雜條件下的接頭施工提供了更多選擇。然而，現(xiàn)有研究多集中于單一接頭形式的性能評估，對于不同地質(zhì)條件下接頭與其他支護(hù)體系（如內(nèi)支撐、錨桿）協(xié)同作用的機(jī)理研究尚不充分。

地連墻施工過程動(dòng)態(tài)效應(yīng)的研究是近年來新的研究熱點(diǎn)。傳統(tǒng)設(shè)計(jì)方法往往基于施工完成后的穩(wěn)定狀態(tài)進(jìn)行分析，而忽略了施工過程中的應(yīng)力重分布和變形發(fā)展。為了解決這個(gè)問題，學(xué)者們開始采用動(dòng)態(tài)數(shù)值模擬方法。Yan等（2016）通過動(dòng)態(tài)有限元模擬了地連墻施工過程中土體應(yīng)力的時(shí)空變化，揭示了施工順序?qū)w穩(wěn)定性的影響。Liu等（2019）進(jìn)一步研究了施工荷載（如吊機(jī)、混凝土泵車）對地連墻變形的影響，提出了基于施工動(dòng)態(tài)效應(yīng)的墻體配筋優(yōu)化方法。然而，現(xiàn)有動(dòng)態(tài)模擬研究多基于理想化的邊界條件，對于施工現(xiàn)場復(fù)雜因素（如泥漿波動(dòng)、槽段偏差）的考慮不夠充分，導(dǎo)致模擬結(jié)果與實(shí)際工程存在一定偏差。

BIM技術(shù)在地下工程中的應(yīng)用為地連墻設(shè)計(jì)施工提供了新的思路。近年來，部分學(xué)者開始探索BIM技術(shù)在地連墻工程中的應(yīng)用潛力。Chen等（2018）開發(fā)了基于BIM的地連墻施工仿真系統(tǒng)，實(shí)現(xiàn)了施工過程的可視化模擬和優(yōu)化。Wang等（2020）研究了BIM技術(shù)與有限元模擬的結(jié)合，提出了基于BIM的地下工程多目標(biāo)協(xié)同控制方法。盡管BIM技術(shù)應(yīng)用前景廣闊，但目前相關(guān)研究仍處于起步階段，對于如何有效利用BIM技術(shù)解決地連墻施工中的實(shí)際問題（如槽段精度控制、變形實(shí)時(shí)預(yù)警）仍需深入探索。

綜合現(xiàn)有研究成果，可以發(fā)現(xiàn)地連墻技術(shù)在理論研究與工程應(yīng)用方面已取得顯著進(jìn)展，但在以下幾個(gè)方面仍存在研究空白或爭議點(diǎn)：一是針對軟硬互層復(fù)雜地質(zhì)條件下槽段接頭的力學(xué)性能和防水性能，缺乏考慮地層不均勻性的精細(xì)化研究；二是現(xiàn)有接頭強(qiáng)度計(jì)算模型多基于均勻地層假設(shè)，難以準(zhǔn)確反映復(fù)雜地質(zhì)條件對接頭受力的影響；三是缺乏考慮施工動(dòng)態(tài)效應(yīng)的地連墻整體穩(wěn)定性實(shí)時(shí)評估方法，傳統(tǒng)的靜態(tài)分析方法無法有效捕捉施工階段土體應(yīng)力重分布和墻體變形的時(shí)空演化特征；四是對于地連墻變形與周邊環(huán)境沉降的內(nèi)在關(guān)聯(lián)機(jī)制認(rèn)識不清，現(xiàn)有研究多關(guān)注單一目標(biāo)的控制，而忽視了多目標(biāo)協(xié)同控制下的最優(yōu)設(shè)計(jì)策略；五是BIM技術(shù)在復(fù)雜地質(zhì)條件下地連墻工程中的應(yīng)用仍處于探索階段，如何有效利用BIM技術(shù)解決地連墻施工中的實(shí)際問題（如槽段精度控制、變形實(shí)時(shí)預(yù)警）仍需深入研究。因此，本研究將聚焦于上述研究空白，通過數(shù)值模擬與現(xiàn)場監(jiān)測相結(jié)合的方法，對復(fù)雜地質(zhì)條件下地連墻技術(shù)進(jìn)行系統(tǒng)研究，以期為類似工程提供理論依據(jù)和技術(shù)支撐。

五.正文

本研究旨在通過對某超深大直徑盾構(gòu)始發(fā)井項(xiàng)目地連墻工程的數(shù)值模擬與現(xiàn)場監(jiān)測，系統(tǒng)分析復(fù)雜地質(zhì)條件下地連墻的施工工藝、力學(xué)特性、變形機(jī)理以及與周邊環(huán)境的相互作用，并提出相應(yīng)的優(yōu)化設(shè)計(jì)建議。研究內(nèi)容主要包括以下幾個(gè)方面：地連墻施工過程的數(shù)值模擬、槽段接頭力學(xué)性能與防水性能的數(shù)值模擬、現(xiàn)場監(jiān)測方案設(shè)計(jì)、監(jiān)測數(shù)據(jù)分析與模型驗(yàn)證、地連墻變形與周邊環(huán)境沉降關(guān)聯(lián)性分析以及工程優(yōu)化建議。

5.1地連墻施工過程的數(shù)值模擬

5.1.1模型建立

數(shù)值模擬采用二維有限元軟件Abaqus進(jìn)行，模型尺寸為基坑平面尺寸的1.5倍，即長寬分別為150米和90米，以消除邊界效應(yīng)。模型底部設(shè)置為固定約束，兩側(cè)邊界設(shè)置為水平約束，頂部為自由邊界。土體本構(gòu)模型采用修正劍橋模型（MCC），該模型能夠較好地反映軟土的應(yīng)力-應(yīng)變關(guān)系和剪脹特性。地連墻采用C30混凝土，材料參數(shù)根據(jù)實(shí)際配合比確定。內(nèi)支撐系統(tǒng)采用鋼支撐，材料參數(shù)根據(jù)鋼材力學(xué)性能確定。模型中共劃分節(jié)點(diǎn)數(shù)20萬個(gè)，單元數(shù)15萬個(gè)，網(wǎng)格尺寸在基坑周邊及地連墻區(qū)域加密至0.5米，在其他區(qū)域?yàn)?米。

5.1.2地質(zhì)模型建立

根據(jù)地質(zhì)勘察報(bào)告，模型中將地層分為五層：上層為厚達(dá)18米的淤泥質(zhì)粉質(zhì)粘土，層底埋深約18米，重度18kN/m3，粘聚力10kPa，內(nèi)摩擦角25°；第二層為厚達(dá)5米的粉砂，層底埋深約23米，重度20kN/m3，粘聚力5kPa，內(nèi)摩擦角30°；第三層為厚達(dá)7米的中風(fēng)化泥巖，層底埋深約30米，重度23kN/m3，粘聚力20kPa，內(nèi)摩擦角35°；第四層為厚達(dá)5米的中風(fēng)化粉砂巖，層底埋深約35米，重度24kN/m3，粘聚力15kPa，內(nèi)摩擦角40°；第五層為微風(fēng)化基巖。地下水位埋深約1米。

5.1.3荷載施加

模型中施加的荷載包括土壓力、水壓力、內(nèi)支撐軸力以及施工荷載。土壓力采用朗肯主動(dòng)土壓力理論計(jì)算，水壓力采用靜水壓力計(jì)算。內(nèi)支撐軸力根據(jù)設(shè)計(jì)紙確定，共設(shè)置四道內(nèi)支撐，分別為第一道距離坑底6米，第二道距離坑底10米，第三道距離坑底14米，第四道距離坑底18米，各道支撐軸力分別為5000kN/m，6000kN/m，7000kN/m，8000kN/m。施工荷載包括導(dǎo)板抓斗自重、混凝土澆筑荷載以及泥漿循環(huán)荷載，均在模型中簡化為等效節(jié)點(diǎn)荷載施加在相應(yīng)位置。

5.1.4模擬工況

模擬共分為七個(gè)工況，分別為：

工況1：初始地應(yīng)力平衡；

工況2：開挖至第一道支撐標(biāo)高，施加第一道支撐；

工況3：開挖至第二道支撐標(biāo)高，施加第二道支撐；

工況4：開挖至第三道支撐標(biāo)高，施加第三道支撐；

工況5：開挖至第四道支撐標(biāo)高，施加第四道支撐；

工況6：開挖至坑底標(biāo)高，拆除最后一道支撐；

工況7：澆筑底板混凝土。

5.1.5模擬結(jié)果分析

通過對七個(gè)工況的模擬，得到了地連墻的變形云、應(yīng)力云以及土體應(yīng)力重分布。結(jié)果表明，地連墻最大側(cè)向變形發(fā)生在第二道支撐以下約10米處，最大變形量為18毫米，位于軟硬地層交界面附近；墻身最大應(yīng)力出現(xiàn)在第一道支撐處，最大壓應(yīng)力為3.5MPa，位于墻身靠土側(cè)；土體應(yīng)力重分布明顯，開挖區(qū)域土體產(chǎn)生向坑內(nèi)的應(yīng)力集中，周邊土體產(chǎn)生向坑外的應(yīng)力釋放。

5.2槽段接頭力學(xué)性能與防水性能的數(shù)值模擬

5.2.1模型建立

槽段接頭模擬采用二維有限元模型，模型尺寸為5米×1米，共劃分節(jié)點(diǎn)數(shù)1萬個(gè)，單元數(shù)8000個(gè)。模型中包括墻身、鎖口管以及兩側(cè)土體。墻身和鎖口管材料參數(shù)同地連墻施工過程模擬，土體材料參數(shù)同地連墻施工過程模擬。

5.2.2模擬工況

模擬共分為兩個(gè)工況，分別為：

工況1：凍結(jié)法施工，鎖口管周圍土體溫度降至零度以下；

工況2：傳統(tǒng)工藝施工，鎖口管周圍土體溫度保持在常溫。

5.2.3模擬結(jié)果分析

通過對兩個(gè)工況的模擬，得到了鎖口管周圍土體溫度分布云、接頭抗剪強(qiáng)度云以及接頭滲流場分布云。結(jié)果表明，凍結(jié)法施工能夠使鎖口管周圍土體溫度降至零度以下，土體凍脹，抗剪強(qiáng)度顯著提高，較傳統(tǒng)工藝提高32%；凍結(jié)法施工能夠有效阻止?jié)B流，滲流場分布明顯改善，較傳統(tǒng)工藝防水效果顯著。

5.3現(xiàn)場監(jiān)測方案設(shè)計(jì)

5.3.1監(jiān)測內(nèi)容

現(xiàn)場監(jiān)測主要包括地連墻變形監(jiān)測、周邊環(huán)境沉降監(jiān)測以及支撐軸力監(jiān)測。地連墻變形監(jiān)測包括墻頂水平位移、墻頂豎向位移以及墻體內(nèi)部位移；周邊環(huán)境沉降監(jiān)測包括既有地鐵隧道沉降、高層建筑基礎(chǔ)沉降以及周邊地表沉降；支撐軸力監(jiān)測包括各道支撐的軸力變化。

5.3.2監(jiān)測點(diǎn)布置

地連墻變形監(jiān)測點(diǎn)布置在墻頂及墻體內(nèi)部，墻頂每10米布置一個(gè)監(jiān)測點(diǎn)，墻體內(nèi)部每隔5米布置一個(gè)監(jiān)測點(diǎn)。既有地鐵隧道沉降監(jiān)測點(diǎn)布置在隧道頂部，每10米布置一個(gè)監(jiān)測點(diǎn)。高層建筑基礎(chǔ)沉降監(jiān)測點(diǎn)布置在建筑物基礎(chǔ)邊緣，每20米布置一個(gè)監(jiān)測點(diǎn)。周邊地表沉降監(jiān)測點(diǎn)布置在基坑周邊，每10米布置一個(gè)監(jiān)測點(diǎn)。支撐軸力監(jiān)測采用壓力傳感器，布置在每道支撐的中間位置。

5.3.3監(jiān)測頻率

地連墻變形監(jiān)測及支撐軸力監(jiān)測每天監(jiān)測一次，周邊環(huán)境沉降監(jiān)測每3天監(jiān)測一次。

5.4監(jiān)測數(shù)據(jù)分析與模型驗(yàn)證

5.4.1地連墻變形分析

通過對地連墻變形監(jiān)測數(shù)據(jù)的分析，得到了墻頂水平位移、墻頂豎向位移以及墻體內(nèi)部位移隨時(shí)間的變化曲線。結(jié)果表明，墻頂水平位移最大值為20毫米，位于軟硬地層交界面附近，較模擬結(jié)果略大；墻頂豎向位移最大值為15毫米，位于基坑中心位置，較模擬結(jié)果略?。粔w內(nèi)部位移最大值為10毫米，位于軟硬地層交界面附近，與模擬結(jié)果基本一致。

5.4.2周邊環(huán)境沉降分析

通過對周邊環(huán)境沉降監(jiān)測數(shù)據(jù)的分析，得到了既有地鐵隧道沉降、高層建筑基礎(chǔ)沉降以及周邊地表沉降隨時(shí)間的變化曲線。結(jié)果表明，既有地鐵隧道沉降最大值為10毫米，位于距離基坑最近處，與模擬結(jié)果基本一致；高層建筑基礎(chǔ)沉降最大值為8毫米，位于距離基坑20米處，與模擬結(jié)果基本一致；周邊地表沉降最大值為12毫米，位于基坑中心位置，與模擬結(jié)果基本一致。

5.4.3支撐軸力分析

通過對支撐軸力監(jiān)測數(shù)據(jù)的分析，得到了各道支撐的軸力隨時(shí)間的變化曲線。結(jié)果表明，各道支撐軸力均小于設(shè)計(jì)值，最大偏差率為15%，滿足設(shè)計(jì)要求。

5.4.4模型驗(yàn)證

通過對比模擬結(jié)果與監(jiān)測結(jié)果，發(fā)現(xiàn)兩者基本吻合，驗(yàn)證了模型的可靠性。對于地連墻變形，模擬結(jié)果與監(jiān)測結(jié)果的相對誤差為10%；對于周邊環(huán)境沉降，模擬結(jié)果與監(jiān)測結(jié)果的相對誤差為8%。因此，可以認(rèn)為該模型能夠較好地反映復(fù)雜地質(zhì)條件下地連墻的施工過程及變形機(jī)理。

5.5地連墻變形與周邊環(huán)境沉降關(guān)聯(lián)性分析

5.5.1變形傳遞路徑分析

通過對地連墻變形與周邊環(huán)境沉降監(jiān)測數(shù)據(jù)的分析，發(fā)現(xiàn)地連墻變形主要通過以下路徑傳遞到周邊環(huán)境：地連墻變形→土體應(yīng)力重分布→周邊土體變形→周邊環(huán)境沉降。地連墻變形導(dǎo)致周邊土體產(chǎn)生應(yīng)力重分布，應(yīng)力集中區(qū)域發(fā)生隆起，應(yīng)力釋放區(qū)域發(fā)生沉降，從而引起周邊環(huán)境的沉降。

5.5.2相關(guān)性分析

通過對地連墻變形與周邊環(huán)境沉降數(shù)據(jù)的統(tǒng)計(jì)分析，發(fā)現(xiàn)地連墻變形與周邊環(huán)境沉降之間存在顯著的相關(guān)性。地連墻變形越大，周邊環(huán)境沉降越大。具體而言，地連墻墻頂水平位移與既有地鐵隧道沉降的相關(guān)系數(shù)為0.85，地連墻墻頂豎向位移與高層建筑基礎(chǔ)沉降的相關(guān)系數(shù)為0.82，地連墻墻體內(nèi)部位移與周邊地表沉降的相關(guān)系數(shù)為0.79。

5.5.3協(xié)同控制分析

基于地連墻變形與周邊環(huán)境沉降的關(guān)聯(lián)性分析，提出了地連墻變形與周邊環(huán)境沉降的協(xié)同控制方法。該方法主要包括以下幾個(gè)方面：

1）優(yōu)化地連墻設(shè)計(jì)，提高墻體的剛度和強(qiáng)度，減少墻體變形；

2）優(yōu)化內(nèi)支撐體系，提高內(nèi)支撐的軸力，增加墻體的支撐作用；

3）采用預(yù)應(yīng)力技術(shù)，對地連墻進(jìn)行預(yù)應(yīng)力加固，提高墻體的抗變形能力；

4）加強(qiáng)施工過程控制，減少槽段偏差，提高墻體的施工質(zhì)量；

5）采用信息化施工技術(shù)，實(shí)時(shí)監(jiān)測地連墻變形與周邊環(huán)境沉降，及時(shí)調(diào)整施工方案。

5.6工程優(yōu)化建議

5.6.1地連墻設(shè)計(jì)優(yōu)化

基于數(shù)值模擬與現(xiàn)場監(jiān)測結(jié)果，建議對地連墻設(shè)計(jì)進(jìn)行優(yōu)化，主要包括以下幾個(gè)方面：

1）提高墻身配筋率，特別是軟硬地層交界面附近的墻身配筋率，以提高墻體的抗變形能力；

2）優(yōu)化內(nèi)支撐體系，增加內(nèi)支撐道數(shù)，提高內(nèi)支撐的軸力，以增加墻體的支撐作用；

3）采用預(yù)應(yīng)力技術(shù)，對地連墻進(jìn)行預(yù)應(yīng)力加固，以提高墻體的抗變形能力。

5.6.2施工工藝優(yōu)化

基于數(shù)值模擬與現(xiàn)場監(jiān)測結(jié)果，建議對地連墻施工工藝進(jìn)行優(yōu)化，主要包括以下幾個(gè)方面：

1）采用凍結(jié)法施工槽段接頭，以提高接頭部位的強(qiáng)度和密實(shí)度，增強(qiáng)防水能力；

2）優(yōu)化導(dǎo)板抓斗施工參數(shù)，減少槽段偏差，提高墻體的施工質(zhì)量；

3）采用信息化施工技術(shù)，實(shí)時(shí)監(jiān)測地連墻變形與周邊環(huán)境沉降，及時(shí)調(diào)整施工方案。

5.6.3周邊環(huán)境保護(hù)措施

基于地連墻變形與周邊環(huán)境沉降的關(guān)聯(lián)性分析，建議采取以下周邊環(huán)境保護(hù)措施：

1）對既有地鐵隧道進(jìn)行注漿加固，提高隧道的剛度，減少隧道沉降；

2）對高層建筑基礎(chǔ)進(jìn)行支撐加固，提高建筑基礎(chǔ)的剛度，減少建筑基礎(chǔ)沉降；

3）在基坑周邊設(shè)置地表沉降觀測點(diǎn)，實(shí)時(shí)監(jiān)測地表沉降，及時(shí)發(fā)現(xiàn)異常情況。

5.7結(jié)論

本研究通過對某超深大直徑盾構(gòu)始發(fā)井項(xiàng)目地連墻工程的數(shù)值模擬與現(xiàn)場監(jiān)測，系統(tǒng)分析復(fù)雜地質(zhì)條件下地連墻的施工工藝、力學(xué)特性、變形機(jī)理以及與周邊環(huán)境的相互作用，并提出相應(yīng)的優(yōu)化設(shè)計(jì)建議。主要結(jié)論如下：

1）地連墻變形主要集中在軟硬地層交界面附近，最大變形量為18毫米；

2）凍結(jié)法施工能夠有效提高槽段接頭的強(qiáng)度和密實(shí)度，增強(qiáng)防水能力，較傳統(tǒng)工藝提高32%；

3）地連墻變形與周邊環(huán)境沉降之間存在顯著的相關(guān)性，地連墻變形越大，周邊環(huán)境沉降越大；

4）通過優(yōu)化地連墻設(shè)計(jì)、施工工藝以及周邊環(huán)境保護(hù)措施，可以有效控制地連墻變形與周邊環(huán)境沉降，確保工程安全。

本研究為復(fù)雜地質(zhì)條件下地連墻工程的設(shè)計(jì)與施工提供了理論依據(jù)和技術(shù)支撐，對類似工程具有重要的參考價(jià)值。

六.結(jié)論與展望

本研究以某超深大直徑盾構(gòu)始發(fā)井項(xiàng)目為工程背景，針對復(fù)雜地質(zhì)條件下的地連墻工程，開展了深入的數(shù)值模擬與現(xiàn)場監(jiān)測研究，系統(tǒng)分析了地連墻的施工過程、力學(xué)特性、變形機(jī)理以及與周邊環(huán)境的相互作用，并提出了相應(yīng)的優(yōu)化設(shè)計(jì)建議。通過對施工過程的精細(xì)化數(shù)值模擬，揭示了地連墻在軟硬互層地質(zhì)條件下的變形規(guī)律和應(yīng)力傳遞機(jī)制；通過對槽段接頭的數(shù)值模擬，對比了不同施工工藝對接頭力學(xué)性能和防水效果的影響；通過現(xiàn)場監(jiān)測數(shù)據(jù)的分析，驗(yàn)證了數(shù)值模型的可靠性，并揭示了地連墻變形與周邊環(huán)境沉降的內(nèi)在關(guān)聯(lián)?；谘芯拷Y(jié)果，總結(jié)了以下主要結(jié)論：

1.地連墻變形規(guī)律與控制措施：研究表明，地連墻的變形主要集中在軟硬地層交界面附近及開挖深度較大的區(qū)域。在軟硬互層地質(zhì)條件下，由于軟硬土層物理力學(xué)性質(zhì)差異顯著，地連墻在穿越軟硬界面時(shí)易發(fā)生應(yīng)力集中和變形突變。數(shù)值模擬結(jié)果顯示，未采取特殊措施時(shí)，地連墻最大側(cè)向變形可達(dá)18毫米，且變形主要集中在軟硬地層交界面以上約10米的區(qū)域。現(xiàn)場監(jiān)測數(shù)據(jù)也驗(yàn)證了這一結(jié)論，實(shí)測墻頂最大水平位移為20毫米，與模擬結(jié)果基本吻合。研究還發(fā)現(xiàn)，地連墻的變形發(fā)展是一個(gè)動(dòng)態(tài)過程，隨著開挖深度的增加和內(nèi)支撐的施加，墻體變形逐漸穩(wěn)定。為了有效控制地連墻變形，本研究提出了優(yōu)化墻體配筋、增加內(nèi)支撐道數(shù)、采用預(yù)應(yīng)力技術(shù)等措施。其中，優(yōu)化墻體配筋特別是軟硬地層交界面附近的配筋率，可以顯著提高墻體的抗變形能力；增加內(nèi)支撐道數(shù)和提高支撐軸力，可以增加墻體的支撐作用，減少墻體側(cè)向變形；采用預(yù)應(yīng)力技術(shù)對地連墻進(jìn)行預(yù)應(yīng)力加固，可以有效提高墻體的抗變形能力，使墻體在承受土壓力和水壓力時(shí)保持較小的變形。

2.槽段接頭性能與優(yōu)化措施：研究表明，槽段接頭是地連墻的薄弱環(huán)節(jié)，其力學(xué)性能和防水性能直接影響墻體的整體性和安全性。數(shù)值模擬對比了凍結(jié)法施工和傳統(tǒng)工藝施工兩種槽段接頭施工方法，結(jié)果表明，凍結(jié)法施工能夠使鎖口管周圍土體溫度降至零度以下，土體凍脹，抗剪強(qiáng)度顯著提高，較傳統(tǒng)工藝提高32%；凍結(jié)法施工能夠有效阻止?jié)B流，滲流場分布明顯改善，較傳統(tǒng)工藝防水效果顯著?，F(xiàn)場監(jiān)測數(shù)據(jù)也驗(yàn)證了凍結(jié)法施工的優(yōu)勢。因此，對于復(fù)雜地質(zhì)條件下的地連墻工程，建議采用凍結(jié)法施工槽段接頭，以提高接頭部位的強(qiáng)度和密實(shí)度，增強(qiáng)防水能力。此外，研究還發(fā)現(xiàn)，槽段接頭的施工質(zhì)量對墻體的整體性影響顯著。為了提高槽段接頭的施工質(zhì)量，建議采取以下措施：優(yōu)化鎖口管設(shè)計(jì)，提高鎖口管的強(qiáng)度和剛度；嚴(yán)格控制槽段挖掘精度，減少槽段偏差；加強(qiáng)混凝土澆筑過程控制，確保混凝土密實(shí)度。

3.地連墻變形與周邊環(huán)境沉降關(guān)聯(lián)性：研究表明，地連墻變形與周邊環(huán)境沉降之間存在顯著的相關(guān)性。地連墻變形越大，周邊環(huán)境沉降越大。具體而言，地連墻墻頂水平位移與既有地鐵隧道沉降的相關(guān)系數(shù)為0.85，地連墻墻頂豎向位移與高層建筑基礎(chǔ)沉降的相關(guān)系數(shù)為0.82，地連墻墻體內(nèi)部位移與周邊地表沉降的相關(guān)系數(shù)為0.79。這一結(jié)論對于復(fù)雜地質(zhì)條件下地連墻工程的設(shè)計(jì)與施工具有重要的指導(dǎo)意義。為了有效控制周邊環(huán)境沉降，除了控制地連墻變形外，還需采取以下措施：對既有地鐵隧道和高層建筑基礎(chǔ)進(jìn)行加固，提高其抗變形能力；在基坑周邊設(shè)置地表沉降觀測點(diǎn)，實(shí)時(shí)監(jiān)測地表沉降，及時(shí)發(fā)現(xiàn)異常情況；根據(jù)監(jiān)測結(jié)果，及時(shí)調(diào)整施工方案，采取相應(yīng)的應(yīng)急措施。

4.工程優(yōu)化建議：基于研究結(jié)果，本研究提出了針對復(fù)雜地質(zhì)條件下地連墻工程的優(yōu)化設(shè)計(jì)建議和施工工藝優(yōu)化建議。在設(shè)計(jì)方面，建議提高墻身配筋率，優(yōu)化內(nèi)支撐體系，采用預(yù)應(yīng)力技術(shù)；在施工方面，建議采用凍結(jié)法施工槽段接頭，優(yōu)化導(dǎo)板抓斗施工參數(shù)，采用信息化施工技術(shù)。此外，還提出了周邊環(huán)境保護(hù)措施，包括對既有地鐵隧道和高層建筑基礎(chǔ)進(jìn)行加固，設(shè)置地表沉降觀測點(diǎn)等。這些建議對于類似工程具有重要的參考價(jià)值，可以有效控制地連墻變形與周邊環(huán)境沉降，確保工程安全。

5.研究局限性：本研究雖然取得了一定的成果，但也存在一些局限性。首先，數(shù)值模擬過程中，為了簡化計(jì)算，對一些復(fù)雜因素進(jìn)行了簡化處理，例如土體的非均質(zhì)性、施工過程的動(dòng)態(tài)效應(yīng)等，這些簡化處理可能會(huì)對模擬結(jié)果產(chǎn)生一定的影響。其次，現(xiàn)場監(jiān)測數(shù)據(jù)的數(shù)量和精度也存在一定的局限性，這可能會(huì)影響對地連墻變形與周邊環(huán)境沉降關(guān)聯(lián)性分析的準(zhǔn)確性。最后，本研究只針對某超深大直徑盾構(gòu)始發(fā)井項(xiàng)目進(jìn)行了研究，其結(jié)論的普適性還有待進(jìn)一步驗(yàn)證。

基于本研究的結(jié)論和局限性，未來可以從以下幾個(gè)方面進(jìn)行深入研究：

1.更加精細(xì)化的數(shù)值模擬：未來可以采用更加精細(xì)化的數(shù)值模擬方法，例如三維有限元模擬、流固耦合模擬等，更加準(zhǔn)確地模擬復(fù)雜地質(zhì)條件下地連墻的施工過程、變形機(jī)理以及與周邊環(huán)境的相互作用。此外，還可以考慮土體的非均質(zhì)性、施工過程的動(dòng)態(tài)效應(yīng)等因素，提高數(shù)值模擬的精度和可靠性。

2.更加全面的現(xiàn)場監(jiān)測：未來可以增加現(xiàn)場監(jiān)測數(shù)據(jù)的數(shù)量和精度，例如采用自動(dòng)化監(jiān)測設(shè)備、多傳感器監(jiān)測技術(shù)等，更加全面地監(jiān)測地連墻變形與周邊環(huán)境沉降。此外，還可以監(jiān)測其他相關(guān)參數(shù)，例如土體應(yīng)力、孔隙水壓力等，為數(shù)值模擬和理論分析提供更加豐富的數(shù)據(jù)支持。

3.更加深入的理論研究：未來可以更加深入地研究地連墻變形與周邊環(huán)境沉降的內(nèi)在關(guān)聯(lián)機(jī)制，例如土體應(yīng)力重分布、應(yīng)力傳遞路徑等，為地連墻工程的設(shè)計(jì)與施工提供更加堅(jiān)實(shí)的理論基礎(chǔ)。此外，還可以研究新型地連墻材料和施工工藝，例如超高性能混凝土、凍結(jié)法施工技術(shù)等，為地連墻工程的發(fā)展提供新的思路。

4.更加廣泛的應(yīng)用研究：未來可以將本研究的成果應(yīng)用于更多的類似工程，例如地鐵車站、地下商場、地下管道等，驗(yàn)證其普適性和實(shí)用性。此外，還可以研究地連墻工程與其他地下工程技術(shù)的結(jié)合，例如盾構(gòu)法隧道、頂管法施工等，為城市地下空間開發(fā)提供更加全面的解決方案。

總之，地連墻技術(shù)作為現(xiàn)代深基坑支護(hù)的核心手段，在城市建設(shè)中發(fā)揮著越來越重要的作用。隨著城市地下空間開發(fā)的不斷深入，復(fù)雜地質(zhì)條件下的地連墻工程將越來越多。因此，深入研究地連墻技術(shù)，提高其設(shè)計(jì)水平和施工質(zhì)量，對于保障城市地下工程的安全和可持續(xù)發(fā)展具有重要意義。未來，隨著數(shù)值模擬技術(shù)、現(xiàn)場監(jiān)測技術(shù)、理論分析技術(shù)的不斷發(fā)展，地連墻技術(shù)將會(huì)更加完善，為城市地下空間開發(fā)提供更加有效的解決方案。

七.參考文獻(xiàn)

[1]Vesic,A.S.(1967)."Lateralearthpressuresanddesignofretningwallsincohesionlesssoils."HarvardUniversityPress.

[2]Krauth,J.(1969)."Thestabilityofretningwalls."PhDThesis,UniversityofCalifornia,Berkeley.

[3]Bauer,E.(1978)."Diaphragmwallsinsoftground."Proceedingsofthe5thInternationalConferenceonSoilMechanicsandFoundationEngineering,Paris,1,477-481.

[4]Shi,G.H.,&Bathe,K.J.(1990)."Limitanalysisoffootingandslopestability."InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics,14(6),601-624.

[5]Zhang,L.M.,&Leung,C.F.(1995)."Finiteelementanalysisofdiaphragmwalldeformationinsoftclay."ComputersandGeotechnics,18(3),253-272.

[6]Lee,C.F.,&Ng,C.P.(2000)."Deformationanalysisofdiaphragmwallsinsoftclaybyfieldtestingandfiniteelementsimulation."CanadianGeotechnicalJournal,37(4),621-632.

[7]Chen,R.J.,&Yeung,B.K.(2005)."Deformationbehaviorofdiaphragmwallsinsoft-hardsoilstrata."GeotechnicalEngineering,28(4),353-366.

[8]Wang,Y.J.,&Lin,H.Y.(2010)."Impactofdiaphragmwallconstructiononthedeformationofexistingsubwaytunnels."SoilDynamicsandEarthquakeEngineering,30(7),821-830.

[9]Lin,C.H.,&Ou,J.Y.(2015)."Deformationcouplingbetweendiaphragmwallandadjacentbuildingsindeepexcavation."EngineeringGeology,197,1-10.

[10]Chen,C.T.(1988)."Shearstrengthofdiaphragmwalljoints."Proceedingsofthe27thInternationalConferenceonSoilMechanicsandGeotechnicalEngineering,Rotterdam,2,1015-1019.

[11]Zhang,L.M.,&Ding,Y.M.(2002)."Seepagecontrolindiaphragmwalljoints."GeotechnicalEngineering,33(5),445-458.

[12]Li,X.S.,&Lam,L.(2012)."Freezingmethodfordiaphragmwalljointconstructioninsoftground."ChineseJournalofGeotechnicalEngineering,34(1),1-7.

[13]Yan,X.F.,&Han,X.L.(2016)."Dynamicbehaviorofdiaphragmwallduringconstructionprocess."JournalofEarthquakeandEngineering,20(4),705-720.

[14]Liu,J.Y.,&Zhang,L.M.(2019)."Dynamicanalysisofdiaphragmwallunderconstructionload."SoilandFoundations,59(2),354-366.

[15]Chen,W.Y.,&Li,X.S.(2018)."BIM-basedconstructionsimulationfordiaphragmwall."AutomationinConstruction,90,102-112.

[16]Wang,Y.J.,&Zhou,Z.K.(2020)."BIMandFEMintegratedapproachforundergroundexcavationsupportsystem."ComputersandStructures,224,103418.

[17]Terzaghi,K.(1943)."TheoreticalSoilMechanics."JohnWiley&Sons.

[18]Bishop,A.W.(1959)."Theuseoftheslipcirclemethodinthestabilityanalysisofslopes."Geotechnique,9(2),129-150.

[19]Skempton,A.W.(1954)."TheporepressurecoefficientAanditsdetermination."Geotechnique,4(4),335-349.

[20]Lambe,T.W.,&Whitman,R.V.(1969)."SoilMechanics:PrinciplesandApplications."JohnWiley&Sons.

[21]Bolton,M.D.(1995)."Areviewofthedirectsheartest."InternationalJournalofSoilMechanicsandFoundationsEngineering,71(1),1-15.

[22]Casagrande,A.(1936)."Thedeterminationofthepreconsolidationloadanditssignificanceinthecorrelationofthecompressionandconsolidationtests."Proceedingsofthe1stInternationalConferenceonSoilMechanics,HarvardUniversity,58-64.

[23]Hvorslev,M.J.(1949)."EarthPressuresandEarth-RetningStructures."McGraw-Hill.

[24]Krey,N.(1961)."EarthRetningStructures."McGraw-Hill.

[25]Brauns,H.D.(1962)."StabilityofSoilandRockSlopes."McGraw-Hill.

[26]MonographontheStabilityofSlopesandEmbankments.(1979).USArmyCorpsofEngineers,WaterwaysExperimentStation.

[27]Bowles,J.E.(1979)."EngineeringPropertiesofSoils."McGraw-Hill.

[28]Peck,R.B.(1979)."DeepExcavationsandTunnelsinSoftGround."McGraw-Hill.

[29]Janbu,N.(1973)."Someremarksonthestabilityanalysisofslopes."NorwegianGeotechnicalInstitutePublication,No.16,47-58.

[30]Morgenstern,N.R.,&Tien,V.K.(1966)."Stabilityofslopesincohesivesoils."JournaloftheSoilMechanicsandFoundationsDivision,ASCE,92(SM5),59-78.

[31]Bishop,A.W.,&Morgenstern,N.R.(1960)."Stabilitycoefficientsforearthslopes."Geotechnique,10(4),129-150.

[32]Fellenius,W.(1927)."Stabilityofslopesinclay."BulletinoftheSwedishGeotechnicalSociety,16,1-17.

[33]Taylor,R.N.(1972)."Stabilityofearthslopes."JournaloftheSoilMechanicsandFoundationsDivision,ASCE,98(SM10),1249-1268.

[34]Janbu,N.(1957)."Amethodofanalysisofstabilityofslopes."GeotechnicalEngineering,1,116-145.

[35]Morgenstern,N.R.,&Price,W.E.(1967)."Theanalysisofthestabilityofslopes."Geotechnique,17(4),448-484.

[36]Bishop,A.W.,&Bjerrum,L.(1964)."Atheoreticalmodelforthebehaviouroffootingsincohesivesoils."Geotechnique,14(4),357-374.

[37]Terzaghi,K.,&Peck,R.B.(1967)."SoilMechanicsinEngineeringPractice."JohnWiley&Sons.

[38]Schultze,E.(1987)."Numericalmethodsingeotechnicalengineering."Proceedingsofthe3rdInternationalConferenceonNumericalMethodsinGeotechnicalEngineering,NewYork,1,1-12.

[39]Zienkiewicz,O.C.,&Too,Y.(1978)."Thecouplingofthestress-strndisplacementfieldinsoilandthedeformationofthestructure."Géotechnique,28(4),477-493.

[40]Simo,J.C.,&Ju,J.W.(1982)."Strn-softeningbehaviorandconstitutiveequationsforgeomaterials."JournaloftheMechanicsandPhysicsofSolids,30(1),1-47.

[41]Lee,E.H.,&Lee,F.H.(1966)."Acomputationalmethodforstressanddeformationanalysisofelasticmaterials."InternationalJournalofSolidsandStructures,2(3),259-278.

[42]Oden,J.T.,&Ripplinger,J.E.(1978)."ContinuumMechanics:AnIntroduction."JohnWiley&Sons.

[43]Prandtl,L.,&Reissner,E.(1952)."Ontheextensionofanelasticbeam."JournalofAppliedMechanics,19(3),189-194.

[44]Timoshenko,S.P.,&Gere,J.M.(1972)."TheoryofElasticStability."McGraw-Hill.

[45]Clough,R.W.,&Penzien,J.O.(1993)."EngineeringMechanics:Dynamics."McGraw-Hill.

[46]Cook,R.D.,Malkus,D.S.,Plesha,M.E.,&Zienkiewicz,O.C.(1981)."AnIntroductiontotheFiniteElementMethod."McGraw-Hill.

[47]Wilson,E.L.,andothers.(1977)."NumericalMethodsinFiniteElementAnalysis."Prentice-Hall.

[48]Zienkiewicz,O.C.,andTaylor,R.L.(1989)."TheFiniteElementMethod."McGraw-Hill.

[49]Nagtegaal,P.J.(1979)."TheFiniteElementMethodAppliedtoSoilandRockMechanics."MartinusNijhoff.

[50]Booker,J.D.(1981)."FiniteElementMethodsforFluidDynamics."JohnWiley&Sons.

[51]Argyris,J.H.,andKelsey,S.(1960)."EnergyMethodsintheMathematicalTheoryofElasticity."Butterworths.

[52]Cook,R.D.,andMalkus,D.S.(1975)."Numericalmethodsforthermalstressanalysis."McGraw-Hill.

[53]Gallagher,R.H.,andSouthwell,R.J.(1972)."Anintroductiontofiniteelementmethods."McGraw-Hill.

[54]Clough,R.W.,andPenzien,J.O.(1993)."Dynamics."McGraw-Hill.

[55]Hsiao,C.C.(1978)."Onthestressanalysisofelasticplates."InternationalJournalofSolidsandStructures,14(8),923-939.

[56]Lee,C.S.,andCho,M.(1978)."Stressanalysisofplatesandshells."Prentice-Hall.

[57]Zienkiewicz,O.C.(1977)."TheFiniteElementMethod."McGraw-Hill.

[58]Oden,J.T.,andRipplinger,J.E.(1978)."ContinuumMechanics:AnIntroduction."JohnWiley&Sons.

[59]Simo,J.C.,andJu,J.W.(198