摘要：不規(guī)則超大深基坑的數(shù)值仿真時一方面建模難度較高，另一方面涉及存在大量的結(jié)構(gòu)物與巖、土相互作用的結(jié)構(gòu)面，且存在大量接觸面隨開挖“消失”或“出現(xiàn)”的過程，由此導(dǎo)致單一手段無法準(zhǔn)確模擬不規(guī)則深基坑的位移場和應(yīng)力場。針對這一難點，采用聯(lián)合ANSYS及ABAQUS對不規(guī)則地鐵車站基坑進(jìn)行數(shù)值仿真。利用ANSYS前期建模準(zhǔn)確及ABAQUS非線性處理能力強(qiáng)的優(yōu)勢，獲得不規(guī)則超大深基坑的位移場和應(yīng)力場。數(shù)值分析與實測對比的結(jié)果說明，這一仿真思路值得在大型復(fù)雜工程的計算分析中推廣與應(yīng)用。

關(guān)鍵詞：建模；仿真計算；深基坑；ANSYS；ABAQUS

0" "引言

三維深基坑的仿真涉及眾多墻－土單元的非線性接觸，基坑開挖模擬過程比較復(fù)雜[1-3]，隨著基坑的開挖，會有大量土體單元及接觸單元的“殺死”以及支撐桿單元的生成。針對不規(guī)則超大深基坑，可利用數(shù)值分析手段，對基坑的位移場和應(yīng)力場進(jìn)行分析計算，找到基坑潛在的風(fēng)險區(qū)域，并采取針對性的處理措施，優(yōu)化原有基坑支護(hù)方案，以保證基坑施工的安全性和降低成本。

如何獲得更接近實際的數(shù)值計算結(jié)果，也成超大深基坑數(shù)值仿真計算的難點之一。一方面是因為不規(guī)則超大深基坑建模難度較高，有限元網(wǎng)格剖分較為困難。一般的有限元前處理軟件采用先繪制模型再剖分單元格的方式，很難精準(zhǔn)剖分出形狀適宜的單元格，易出現(xiàn)畸形網(wǎng)格，導(dǎo)致數(shù)值計算難以收斂。且往往需要較高的計算成本，使得數(shù)值計算無法完成迭代[4-5]。另一方面是因為不規(guī)則超大深基坑存在大量的結(jié)構(gòu)物與巖、土相互作用的結(jié)構(gòu)面，且存在大量接觸面隨開挖過程的“消失”或“出現(xiàn)”。這就需要數(shù)值分析軟件能提供合適的指令，模擬這一動態(tài)施工過程，且能提供豐富的接觸面本構(gòu)關(guān)系、巖土體本構(gòu)關(guān)系，來準(zhǔn)確模擬基坑工作性狀[6-7]。

鑒于此，本文綜合ANSYS、ABAQUS兩款通用有限元軟件的優(yōu)點。利用ANSYS靈活的建模能力，應(yīng)用APDL語言編制相應(yīng)的程序，可方便快捷地構(gòu)建不規(guī)則深基坑的模型。ANSYS完成建模后，再用APDL的功能函數(shù)*GET，將不規(guī)則深基坑的模型信息改寫成ABAQUS的命令流格式。ABAQUS比較ANSYS在處理非線性問題方面的優(yōu)點在于計算收斂快，特別是對于基坑開挖求解問題，其對初始地應(yīng)力場的平衡效果是ANSYS軟件所無法比擬的[8-9]。

ABAQUS自帶許多適宜模擬土體變形的本構(gòu)模型，可為模擬巖土變形問題及巖土-結(jié)構(gòu)物接觸問題提供多方位選擇，這也是ANSYS或者其他非線性有限元軟件沒有的功能。此外ABAQUS還提供用戶自定義材料模型的子程序UMAT，方便用戶添加自己的本構(gòu)模型。ABAQUS非線性求解能力，能準(zhǔn)確模擬不規(guī)則深大基坑的位移場和應(yīng)力場，找到基坑三維深基坑潛在的風(fēng)險區(qū)域，以供基坑支護(hù)方案的優(yōu)化與調(diào)整。此方案結(jié)合了ANSYS、ABAQUS兩款通用有限元軟件的優(yōu)點，可在巖土工程設(shè)計與分析計算中推廣與應(yīng)用。

1" "ANSYS及ABAQUS用于基坑有限元前處理

1.1" "應(yīng)用ANSYS方法對基坑建模

對于基坑問題，在定義墻與土體的接觸分析時，往往采用NODE TO SURFACE類型。在CAE中要一個個選擇節(jié)點，當(dāng)節(jié)點很多的時候不是很方便且費時，而ANSYS中，則可十分方便地通過NSEL命令選擇從面上的節(jié)點，還可以通過CM命令將節(jié)點集合、組合等一系列操作，從而避免ABAQUS自帶CAE操作中選擇節(jié)點費時的缺點。另外，也可將后處理中需要用到的節(jié)點單元集導(dǎo)入到ABAQUS中，以方便后處理。例如可以定義斜測管處的節(jié)點為一個節(jié)點集，在INP文件中添加后處理命令輸出該節(jié)點集的斜測值，從而提高后處理效率。

ANSYS主要包括掃瓊劃分、映射劃分、自由劃分和自適應(yīng)劃分等建模方式，也可直接生成有限元模型[10]。復(fù)雜的三維基坑模型需綜合上述各種建模方法。通常模擬土體六面體單元采用掃瓊劃分，模擬支撐桿的三維殼單元采用映射劃分。

1.2" "ANSYS格式輸出命令

利用ANSYS建好的基坑模型并不能直接為ABAQUS軟件所用。第一步需用APDL中的功能函數(shù)*GET獲取完整的三維基坑有限元模型信息，包括了節(jié)點信息、單元信息、邊界條件信息、荷載信息、材料信息、截面特性等；第二步在相應(yīng)位置再用APDL命令寫入ABAQUS關(guān)鍵詞；第三步按照ABAQUS的INP文件所需格式類型。

復(fù)雜深基坑的建模網(wǎng)格密度難以一次剖分到位，利用APDL語言可以輕松修改模型網(wǎng)格密度，并實現(xiàn)模型信息的快速更改和提取。由APDL命令建立的有限元分析模型數(shù)據(jù)文件，往往僅為CAE數(shù)據(jù)文件的1%，不僅便于保存和攜帶，在網(wǎng)上或平常的交流也很方便。

1.3" "ABAQUS方法用于基坑工程分析

聯(lián)合ANSYS、ABAQUS兩款通用有限元軟件，進(jìn)行不規(guī)則深大基坑建模的流程圖如圖1所示。用APDL提出ABAQUS能識別的模型數(shù)據(jù)后，將其存儲為.inp文件，在三維基坑模型數(shù)據(jù)信息后，寫入ABAQUS數(shù)值軟件能識別的分析語言。在.inp文件中需加入以下分析步驟。

第一類分析步驟是在整體基坑模型建成、裝配、并剖分好單元格后，在基坑單元格信息后面，進(jìn)行基坑初始應(yīng)力場的平衡。此時基坑初始應(yīng)力場的平衡采用*Geostatic（地應(yīng)力平衡）分析步實現(xiàn)。

第二類分析步驟是隨著基坑的土體的開挖，分層架設(shè)支撐，對基坑土體開挖和鋼支撐添加分別采用*Model Change、Remove、*Model Change、Add來模擬。

第三類分析步驟是不規(guī)則深基坑采用橫向支撐為預(yù)應(yīng)力支撐，要給支撐施加預(yù)應(yīng)力時，可采用初始條件語句*Initial conditions Type=Stress來設(shè)置。

2" "工程概況

某工程為地下二層島式車站基坑，基坑?xùn)|西向布置。車站總長163.3m，標(biāo)準(zhǔn)段寬度18.9m。地鐵車站基坑為不規(guī)則長條形，如圖2所示?；娱_挖時分東、西兩區(qū)，并進(jìn)行分區(qū)開挖。西部擴(kuò)大段最大寬度37.4m。車站底板埋深約16.5m，基坑兩端深約16.8m，車站凈高13.55m。車站基坑整體呈兩端寬、中間窄的結(jié)構(gòu)形式。

圖3為基坑橫剖面示意圖。由圖3可知，車站基坑共四道橫向支撐，基坑開挖深度范圍內(nèi)以粉質(zhì)黏土和淤泥質(zhì)粉質(zhì)黏土為主。針對這種開挖深度大、土層條件不優(yōu)越的不規(guī)則大型基坑，開挖前進(jìn)行基坑數(shù)值仿真分析，獲得基坑整體的位移場和應(yīng)力場，以優(yōu)化和調(diào)整基坑支護(hù)方案十分有必要。

3" "基坑有限元仿真

按照圖1所示的基坑模擬流程圖，聯(lián)合ANSYS及ABAQUS兩種大型通用有限元計算軟件，對車站基坑支護(hù)和開挖進(jìn)行數(shù)值模擬。采用*Geostatic（地應(yīng)力平衡）分析步驟模擬基坑開挖前的初始地應(yīng)力場。采用*Model Change，Remove模擬土體分層開挖，采用*Model Change、Add來模擬支撐體系架設(shè)。采用分區(qū)域“殺死”、“激活”的方式分步驟，模擬基坑開挖過程及支撐架設(shè)過程。計算分析步驟結(jié)束后，選取圖2所示的CX1、CX2、CX3、CX4等位置基坑位移值及支護(hù)結(jié)構(gòu)應(yīng)力計算值，同實測數(shù)據(jù)對比，以反映計算結(jié)果的準(zhǔn)確性。

3.1" "基坑模型簡化

假設(shè)土體初始位移為零，各層土體采用摩爾庫倫本構(gòu)模型?；咏７治鲞^程中，根據(jù)實際情況假定，模型兩側(cè)邊界為水平位移為零，豎直位移無約束，模型底邊位無位移，模型頂部為自由端、無約束?；拥慕３叽鐬榭傞L163.3m，深度為16.5m。寬度依據(jù)基坑形狀略有改變，分別選擇37.4m（Z1區(qū)寬）、33.5m（Z1、Z2區(qū)寬）及18.9m（標(biāo)準(zhǔn)段寬），未考慮東區(qū)擴(kuò)大寬度。

3.2" "模型參數(shù)選取

三維基坑模型仿真參數(shù)如表1和表2所示。

3.3" "仿真初始模型

有限元模型共有單元46700個單元，其中土體實體單元個數(shù)為36922個，接觸面單元數(shù)為6252個，連續(xù)墻殼單元個數(shù)為3126個，鋼支撐梁單元數(shù)為400個。有限元計算范圍長×寬×高為666m×200m×64m。撐、地下連續(xù)墻、基坑被動區(qū)的墻高為32m，開挖深度為16.5m。三維基坑仿真模擬模型如圖4所示。

4" "仿真結(jié)果分析

4.1" "仿真基坑變形云圖

因模型比較大，選擇地下連續(xù)墻殼單元的位移情況來反映基坑整體變形性狀。圖5a為地應(yīng)力平衡效果圖，由圖5a可以看出，場地豎向位移達(dá)到10-14m級，表明地應(yīng)力平衡效果十分理想。圖5b、圖5c分別為地下連續(xù)墻U1、U2方向的位移云圖。由圖5b、圖5c可知，不管是U1還是U2方向的側(cè)移，在最終開挖步連續(xù)墻的最大水平位移均出現(xiàn)在基坑開挖面附近。

圖5d為最終開挖步下基坑底部位移變化云圖（U3方向位移）。由圖5d可知，坑底最大隆起位移首先出現(xiàn)在基坑擴(kuò)大段，隨后整個基坑底部都慢慢出現(xiàn)較大隆起量?；娱_挖完畢后，坑底最大位移隆起量為100.0mm，與基坑立柱的實測隆起量平均值95.6mm接近，說明了數(shù)值仿真模擬技術(shù)的正確性。

4.2" "計算結(jié)果分析

4.2.1" "墻身位移

地下連續(xù)墻墻身位移曲線及與實測值對比如圖6所示。由圖6可知，最終開挖步下，墻身的最大側(cè)向位移值出現(xiàn)在開挖面附近出現(xiàn)，基坑實際監(jiān)測也得到類似的規(guī)律。基坑連續(xù)墻上的4#測斜管，對應(yīng)的位置為圖6中的CX1，12#測斜管對應(yīng)的位置為圖6中的CX3。將有限元計算得到的地連墻水平位移和測斜管實測值進(jìn)行對比，發(fā)現(xiàn)地下連續(xù)墻測斜位移計算值與實測值擬合效果十分好，進(jìn)一步說明采用該方法對車站基坑數(shù)值模擬結(jié)果的準(zhǔn)確性。

4.2.2" "坑底隆起計算結(jié)果分析

5個開挖步完畢時，西區(qū)基坑擴(kuò)大段底部土體的隆起位移變化曲線如圖7所示。因為兩端有墻體的側(cè)向約束作用，故坑底隆起位移值為中間大兩端小，這一曲線形態(tài)與實測結(jié)果一致。位移曲線顯示出隨著土體開挖深度變深，卸荷量變大，土體回彈量也變大，土體的隆起位移增大。5個開挖步結(jié)束時，坑底最終的最大隆起位移量為100.0mm，與立柱實測的最大隆起位移值105.3mm比較接近?；觾蓚?cè)土體由于墻－土接觸摩擦作用，土體隆起位移增量比較均勻，平均值為10mm左右，最大隆起位移值為43.5mm，可見由于墻－土間的摩擦作用對基坑土體的隆起位移起到了一定，約束作用。

由數(shù)值分析的坑底隆起量較實際監(jiān)測值要偏大，是因為實際開挖過程中，基坑分區(qū)開挖至坑底時，會及時進(jìn)行底板澆筑及施作地下結(jié)構(gòu)。數(shù)值分析時，未考慮坑底土體隆起的加固措施，也未考慮實際的分區(qū)開挖步驟，故數(shù)值計算結(jié)果會比實際工況的隆起量要大。計算結(jié)果的差異性，說明模擬基坑開挖時，進(jìn)一步考量基坑分區(qū)開挖的時空效應(yīng)，才會得到更接近真實情況的坑底土體隆起位移量。

5" "結(jié)束語

本文針對不規(guī)則地鐵深基坑建模難度高及計算分析過程復(fù)雜，分別取ANSYS前期建模準(zhǔn)確及ABAQUS非線性處理能力強(qiáng)的優(yōu)勢，提出聯(lián)合ANSYS及ABAQUS對不規(guī)則地鐵車站基坑進(jìn)行數(shù)值仿真。數(shù)值模擬結(jié)果與實測值的對比，說明了這一思路的可行，值得在大型復(fù)雜工程的計算分析中推廣與應(yīng)用。

參考文獻(xiàn)

[1] 李興磊，周堅，宗鑫.基于ANSYS軟件高層建筑深基坑施

工風(fēng)險分析[J].建筑技術(shù)開發(fā)，2023，50（9）：59-61.

[2] 匡文壯，邵金，張玉峰，等.軟土基坑中長短樁支護(hù)的數(shù)值

模擬研究[J].土工基礎(chǔ)，2023（12）：856-864.

[3] 余庭毒，符必昌.基于FLAC軟件對某基坑穩(wěn)定性的數(shù)值

模擬研究[J].地質(zhì)災(zāi)害與環(huán)境保護(hù)，2020：68-72.

[4] 陳江，孫志浩，徐長節(jié).軟土地區(qū)深大基坑支護(hù)方案優(yōu)化三

維數(shù)值模擬分析[J].科學(xué)技術(shù)與工程，2019，19 （31） ：319-326.

[5] 徐冉.地下連續(xù)墻支護(hù)深基坑開挖過程模擬分析[J].水利

科技與經(jīng)濟(jì)， 2022，28（5）：55- 58.

[6] 王志杰，周飛聰，周平，等.基于強(qiáng)近接大型基坑單側(cè)開挖

卸載既有車站變形理論研究[J].巖石力學(xué)與工程學(xué)報，2020

（10）：2131-2147.

[7] CHEN A X，WANG Q，CHEN Z D， et al.Investigating pile

anchorsupport system for deep foundation pitina

congested area of Changchun[J].Bulletin of Engineering

Geology and the Environment.2021.80 （2）： 1125-1136.

[8] 楊波，孫大齊，周記名.地鐵車站深基坑動態(tài)施工過程變形

規(guī)律與數(shù)值仿真分析[J].河北工程大學(xué)學(xué)報（自然科學(xué)版），

2020，37（2）：38-44.

[9] 胡軍學(xué)，黃子豐，胡裕琛，等. 深大基坑異步分區(qū)開挖的三

維數(shù)值模擬研究[C].2021年全國工程建設(shè)行業(yè)施工技術(shù)交

流會論文集：194-196.

[10] 吳清平，時偉，戚鍥鐘，等.超大深基坑BIM施工全過程

模擬與分析研究[J.工程建設(shè)，2013，45（5）： 20-23.