均勻設(shè)計法在淤泥剛度參數(shù)確定中的應用

張智恒 龐小朝

鐵科院（深圳）研究設(shè)計院有限公司，廣東 深圳 518060

城市發(fā)展過程中，基坑工程不斷增多，對基坑安全性提出了更高要求?；影踩婕盎釉陂_挖過程中支護結(jié)構(gòu)橫向變形的發(fā)展，而橫向變形與土層信息及土體的主要剛度參數(shù)（參考割線模量Eref50、參考切線模量Erefoed、參考卸載再加載模量Erefur)有關(guān)。根據(jù)相關(guān)理論基礎(chǔ)和工程經(jīng)驗，硬化土本構(gòu)（Hardening Soil，HS）模型是反映基坑支護結(jié)構(gòu)橫向變形的常用土體本構(gòu)模型，把加載變形和卸載變形有所區(qū)分，體現(xiàn)了土體的應力-應變關(guān)系和土體的剪脹性。在HS模型基礎(chǔ)上通過有限元數(shù)值仿真、室內(nèi)土工試驗或合理參數(shù)反演手段得到土體的主要剛度參數(shù)。均勻設(shè)計法是一種參數(shù)反演分析方法，其本質(zhì)是通過數(shù)論中一致分布理論對區(qū)間范圍內(nèi)參數(shù)進行均勻篩選，取方差最小的參數(shù)或參數(shù)組合為待反演的參數(shù)或參數(shù)組合。楊蘭強等［1］通過室內(nèi)土工試驗得到土體應力-應變關(guān)系，根據(jù)模擬曲線確定了寧波地區(qū)典型土層的HS模型參數(shù)及其關(guān)系，驗證了HS模型及試驗獲取HS參數(shù)在基坑開挖數(shù)值計算中的適用性。劉書斌等［2］在對硬化土模型及模型參數(shù)進行理論分析的基礎(chǔ)上，基于室內(nèi)不同應力路徑試驗，確定了無錫地區(qū)典型土層HS參數(shù)，通過Z_Soil有限元數(shù)值模擬對比了基坑圍護結(jié)構(gòu)水平位移的計算值和實際測量值，表明土體的最大水平位移出現(xiàn)的區(qū)域相近，且計算值與實際測量值誤差較小，符合反演要求。王衛(wèi)東等［3］通過大量室內(nèi)三軸試驗確定了土體各土層的本構(gòu)模型參數(shù)，并確定了3個主要剛度參數(shù)與壓縮模量之間的比例關(guān)系，分析結(jié)果在上海軟土地區(qū)基坑工程應用中具有可靠性。劉暢［4］用HS模型進行有限元數(shù)值仿真，把土體的3個主要剛度參數(shù)與壓縮模量的關(guān)系進行了對比和擬合分析，反演得到的參數(shù)取值對于天津軟土地區(qū)HS模型的主要參數(shù)選取具有參考借鑒意義。劉關(guān)虎［5］結(jié)合珠海地區(qū)淤泥土層基坑工程，采用均勻設(shè)計法與Plaxis2D有限元軟件計算結(jié)合的方法，對淤泥土層的3個未知剛度參數(shù)進行反演分析，驗證了擬合參數(shù)的合理性。

深圳地區(qū)存在較多的淤泥土層，本文結(jié)合深圳地鐵12號線和平站基坑工程，土體本構(gòu)模型選用HS模型，采用Plaxis2D有限元軟件建立仿真模型并結(jié)合均勻設(shè)計法對淤泥的未知參數(shù)進行反演分析。

1 HS模型

HS模型（圖1）相比其他模型，在有限元數(shù)值模擬中得到的圍護結(jié)構(gòu)土體位移較為合理，其彈性部分合理地實現(xiàn)了加載和卸載模量的分別定義，塑性部分則采用了各向同性硬化準則和非相關(guān)聯(lián)流動法則，并在此基礎(chǔ)上滿足摩爾-庫倫破壞準則，在黏土和軟土地區(qū)有較大的發(fā)揮空間。

圖1 主應力空間中的HS模型屈服面

HS模型中剛度參數(shù)源于標準三軸排水試驗的雙曲線應力-應變關(guān)系（圖2）。HS模型共有11個參數(shù)，根據(jù)Plaxis2D可按以下標準取值：對于砂土，剪脹角ψ可取為φ-30°，對于黏土一般取0°，本次參數(shù)反演分析中淤泥土層也取0°；剛度應力水平相關(guān)冪指數(shù)m對于不同的土體稍有區(qū)分，本次分析m值取1；模量參考應力pref則按照理論定義取100 kPa；卸載再加載泊松比vur在Plaxis2D軟件中一般取0.2；靜止側(cè)壓力系數(shù)K0根據(jù)公式K0=1-sinφ'確定，φ'為土的參考內(nèi)摩擦角；在Plaxis2D軟件中Rf一般默認為0.9。模型3個剛度模量參數(shù)通過本次參數(shù)反演分析方法結(jié)合現(xiàn)場實測值反演確定。

圖2 三軸固結(jié)排水試驗中雙曲線應力-應變關(guān)系

HS模型中土的極限應力狀態(tài)是由土的黏聚力c、內(nèi)摩擦角φ及剪脹角ψ來簡化描述。不同的是，HS模型采用三個不同的輸入剛度可將土體剛度和土體變形描述得更為精確，從而使基坑開挖工程中的安全性更有保障。

2 均勻設(shè)計法

均勻設(shè)計法來源于數(shù)論中的一致分布理論，即所要反演的因素有若干個，有x個水平的試驗可供每個因素挑取。反演過程中，其選取的試驗點僅有x個，再通過x個試驗點和已知條件對比分析選取最符合反演結(jié)果的一組數(shù)據(jù)。均勻設(shè)計法在試驗次數(shù)相同或相近的情況下，將原來的試驗次數(shù)從x2次減少到x次，大幅優(yōu)化了參數(shù)篩選程序；由于所選的點在對應的參數(shù)區(qū)間里呈均勻分散分布，極大限度保證所選的點具有代表性和科學性。

均勻設(shè)計有特定的評價指標，一般通過查均勻設(shè)計表獲得均勻設(shè)計方法。每個均勻設(shè)計表是通過Un（nm）來體現(xiàn)，其中U表示均勻設(shè)計；n為試驗次數(shù)；m為最大因素數(shù)量，即要反演的因素數(shù)量。如U7（74）代表4因素7水平的均勻設(shè)計表，即有4個參數(shù)要反演，一共經(jīng)歷7次試驗。本次均勻設(shè)計法進行位移反演分析的本質(zhì)是尋找一組待反演的參數(shù)，通過將所得參數(shù)代入到有限元軟件進行數(shù)值仿真，使計算結(jié)果的位移值與實測位移逼近即可。目標函數(shù)F()

X為

式中：X為待反演的巖體參數(shù)；fi（X）為土體在第i個量測點對應的圍護結(jié)構(gòu)土體位移計算值；ui為相應的位移量的實測值；k為測點總數(shù)。

3 參數(shù)反演分析

3.1 工程概況

深圳地鐵12號線和平站基坑工程（圖3）位于松福大道與橋和路交叉路口，和平站沿橋和路下方呈東西向敷設(shè)，途經(jīng)建筑密集，與正在施工的穗莞深城際鐵路（高架）換乘，下穿穗莞深段凈高為14.1～14.6 m。車站為地下兩層12 m島式車站，全長230.5 m，標準段寬度為21.1 m，基坑深度16.8～18.7 m，下穿穗莞深城際鐵路段用1道混凝土支撐+4道伺服鋼支撐。

圖3 和平站平面示意

本工程中淤泥一般為灰黑、深灰等色，飽和，流塑。淤泥土層較厚，平均層厚9.1 m，而基坑深度在16.8～18.7 m，淤泥占比50%，遠超其他土層的層厚占比，因此通過本次有限元數(shù)值模擬與均勻設(shè)計法相結(jié)合反演的基坑工程土體剛度參數(shù)近似等于該地區(qū)淤泥土層的剛度參數(shù)。

3.2 Plaxis2D建模

模型（圖4）邊界長76 m、寬36 m，模型取車站區(qū)間含有伺服鋼支撐的基坑部分，基坑側(cè)壁有1道混凝土支撐和4道伺服鋼支撐。地連墻厚度1.2 m，抗拉剛度為3.6×107N，抗彎剛度為4.32×106Pa，抗扭截面系數(shù)為8 mm3，泊松比為0.15?；炷林紊疃葹?.5 m，第1道伺服鋼支撐在基坑中的深度為4.5 m，每道伺服鋼支撐間距為3.0 m?；炷林卧谀Ｐ椭杏缅^桿模擬，伺服鋼支撐的軸力在模型中用點荷載來模擬，軸力伺服系統(tǒng)的實現(xiàn)是依靠有限元模型中隨著土層的開挖地連墻上點荷載逐漸變化來實現(xiàn)的。案例分析中基坑土體為層狀分布，淤泥土層分布在上層，層厚9.1 m，深度為4.0～13.1 m，淤泥土體采用HS模型，其黏聚力為10 kPa，內(nèi)摩擦角為3°。其他土體除素填土采用摩爾-庫倫模型之外，均采用HS模型。其他土層參數(shù)的取值根據(jù)鐵科院（深圳）研究設(shè)計院多年對于深圳地區(qū)土層的參數(shù)反演經(jīng)驗取值得到，土層參考應力100 kPa，泊松比為0.2，見表1。

圖4 基坑建模Plaxis2D模型（單位：m）

表1 土體HS模型未參與反演的參數(shù)取值

3.3 剛度參數(shù)均勻設(shè)計及有限元計算

采用均勻設(shè)計法（UD）與Plaxis2D軟件計算相結(jié)合方法對基坑淤泥土層硬化土模型主要剛度參數(shù)進行參數(shù)反演分析。根據(jù)相關(guān)工程經(jīng)驗及地質(zhì)勘察資料可確定3個未知剛度參數(shù)的取值范圍為：。為使參數(shù)反演數(shù)據(jù)更加均勻合理，取均勻設(shè)計表U37（3712）來確定具體試驗設(shè)計方案，然后根據(jù)其對應的使用表來安排有限元計算，并輸入本構(gòu)參數(shù)進行計算。根據(jù)均勻設(shè)計表原理，選擇U37中列號為1、5、8的列為3個未知剛度參數(shù)的均勻設(shè)計方案，均勻度偏差最小，具有代表性和合理性，見表2和表3。

表2 3個未知剛度參數(shù)的均勻設(shè)計方案

表3 均勻設(shè)計表U37（3712）的均勻度偏差D

將表2中可能成立的組合數(shù)據(jù)代入Plaxis2D有限元軟件中進行計算，對比分析基坑地連墻水平位移的計算值與工程現(xiàn)場實測值，觀察模型結(jié)構(gòu)深度為4.0～13.1 m淤泥土層區(qū)間的地連墻水平位移，計算值與實測值最相近的一組參數(shù)組合視為反演出的主要剛度參數(shù)。經(jīng)計算，取=40.6 MPa時，位移最大值與現(xiàn)場實測位移最大值相差最小，見圖5及表4。

圖5 不同測點的位移

表4 測點位移監(jiān)測值與反演正算值對比

由表4可知，位移監(jiān)測值與反演正算值在測點ZQT15—ZQT17差距不大，均小于3.6 mm，且這三處測點的淤泥土層厚度均超過5 m，淤泥土層對地連墻水平位移的影響最大。測點ZQT13的位移監(jiān)測值與正算值差距達到6 mm，基于該測點淤泥土層厚度僅為3 m，因此實際反演影響不大，反演出的淤泥土層的3個未知剛度參數(shù)具有有效性。將計算值與實測值對比曲線中數(shù)據(jù)代入式（1）計算得到F（X）=72.66 mm2，整體差距較小，計算結(jié)果比較吻合。

Atkinson［6］通過測量剪切波速度得到土體在小應變下的剛度參數(shù)，并給出了土體剛度參數(shù)與土體強度指標的試驗規(guī)律。根據(jù)和平站淤泥土層十字板剪切試驗結(jié)果和反演得到的淤泥土層三個剛度參數(shù)（表5）對比可知，隨著深度的增加，原狀土和重塑土的抗剪強度逐漸增加；MLZ3-SHP-15試驗點埋深4.5～7.5 m為淤泥土層，且其頂層為土體抗剪強度的最小值，底層為土體抗剪強度的最大值。由此可推導與抗剪強度cu之間的大致關(guān)系為

表5 十字板剪切試驗結(jié)果

根據(jù)和平站標準貫入試驗成果表可以得到淤泥土層標準貫入擊數(shù)N最大值為2.0、最小值為1.0、平均值為1.1；不同測點對應的淤泥標貫擊數(shù)隨著淤泥層深度的變化基本保持穩(wěn)定。由此可以推導與標準貫入擊數(shù)N之間的一般關(guān)系公式為

4 結(jié)論

1）Plaxis2D有限元軟件計算值與實測值能夠很好地吻合，淤泥土層雖離散性較大，用HS模型也能得到預期的效果。HS模型除可以模擬硬土之外，對于模擬軟土也具有有效性。

2）采用均勻設(shè)計表法有效減少了參數(shù)反演分析中剛度參數(shù)取值次數(shù)和數(shù)值模擬分析次數(shù)，計算結(jié)果表明將均勻設(shè)計法和有限元數(shù)值模擬計算相結(jié)合的參數(shù)反演分析方法在巖土工程領(lǐng)域具有合理性。