地下管廊與土體接觸的模擬與試驗研究

丁天 裴曉峰 孫虎 項長生

摘 要：文章通過Midas GTS利用有限元模擬對地下管廊和土體間的接觸進行分析，通過彈塑性支撐來搭建接觸界面，利用反作用系數來對地基支撐剛度進行改變，以此來分析地下復合管廊和土體接觸面間的關聯性，同時通過對土體的內摩擦角作出改變以此來對地下管廓受力情況進行研究。得出土體內摩擦角的改變對地下綜合管廊—土體的界面剛度影響很小，界面摩擦角和土體內摩擦角呈正向反饋，擴大土體內的摩擦角將對地下復合廊的強度和變形產生積極影響。根據上述研究表明，可以針對不同質地的土壤通過相對應的措施，對地下復合管廊的沉降和變形達到控制的作用。

關鍵詞：地下綜合管廊;土體;沉降變形;接觸;Midas GTS

中圖分類號：U354 文獻標識碼：A 文章編號：1001-5922（2021）07-0153-06

Simulation and Experimental Study on the Contact between Underground Pipe Gallery and Soil

Ding Tian1， Pei Xiaofeng1， Sun Hu2， Xiang Changsheng2

（1.China Gansu International Corporation for Economic and Technical Cooperation， Lanzhou 730000，China;

2. School of Civil Engineering， Lanzhou University of Technology，Lanzhou 730000，China）

Abstract：In this paper， Midas GTS uses finite element simulation to analyze the contact between the underground pipe gallery and the soil， and builds the contact interface through elastoplastic support， use the reaction coefficient to change the stiffness of the foundation support to analyze the correlation between the underground composite pipe gallery and the contact surface of the soil， at the same time， the internal friction angle of the soil is changed to study the force of the underground pipe profile. It is concluded that the change of friction angle in soil has little influence on the interfacial rigidity of underground composite pipe gallery-soil， and the magnitude of interfacial friction angle is positively related to the magnitude of friction angle in soil. Expanding the friction angle in soil will positively affect the strength and deformation of underground composite pipe gallery.According to the above research， the settlement and deformation of underground composite pipe gallery can be controlled by corresponding measures for different soil textures.

Key words：underground integrated pipe gallery; soil mass; settlement deformation; contact; Midas GTS

0 引言

地下管道走廊始于1833年，建于巴黎的下水道中。那時，構建良好的衛生系統，對防止城市疾病擴散起著重要作用，并且地下管道的修建和開發也應運而生。西方國家經過一個世紀的發展壯大，所構建的地下復合走廊已經十分成熟，其中法國走在隊伍前列，建造的地下管廊公里數位于首位[1]，其由于特殊的管道結構被鋪設在地下，而又因地下管廊發生受力和變形時不易被人察覺，因此，社會時刻高度重視地下走廊的安全穩定性。對土體地下管廊力學特性和沉降變形的研究，使得其未來運行和維護能夠得到更好地技術支撐，以更有效地促使發展綜合地下管廊。Desai等在文獻[2-4]上討論了具有接觸損傷的Desai薄膜元件的厚度，而盧廷浩通過對法向以及切向耦合進行測定，對此進行了深入研究，并對薄層單元模型進行了推廣[5];唐世斌基于損傷力學構建了接觸面模型，考慮了其不均勻性和易損性，得出了剪切損傷是由于高剪應力造成的結論，不均勻性也會在非靜態調整過程中產生破壞性應力[6]。鄧健采取的方式是通過模擬接觸面在演化時彈塑性損傷，結論認為通過模擬接觸面的非線性，其精度能夠達到工業需求[7-8]。國內目前針對土體等接觸做了很多直接剪切實驗[9-10]，接觸性破壞也在土體中有所體現，而直接剪切實驗結果對混凝土或鋼板具有較大剛性的損失準確率不高，后續對此還需做進一步研究。

文章通過對地下管廊與土體接觸的數值進行模擬可得，土體管廊遭受沉降是受到其內摩擦角和基床系數的影響，針對土體的差異采取不能的解決方式，從而對地下綜合管廊的沉降進行控制。

1 理論分析

1.1 接觸界面單元

土體和結構之間發生接觸是經常存在的，土體的剛度和混凝土的有所區別，二者有著較大差異，即使材料不一致，其接觸在有限元軟件通過滑動接觸元件來完成，Goodman單元被廣泛采用[11]。四節點接觸如圖1所示。

σ為接觸單元的正應力，τ為剪應力，二者有著如下關系：

本構關系矩陣D為：

其中kn，ks、分別為軸向剛度和切向剛度，形函數（Ni）如下：

利用參形函數表示單元內任意點位移的公式如下：

接觸單元的形函數為：

單元坐標系的位移公式為：

三角函數值可由式（6）微分得到：

定義三角函數為：

應變增量公式為：

其中ue表示所以單元節點位移向量的集合。

B是整體坐標系形函數的微分，對單元坐標系形函數進行微分得到：

則整體坐標系下：

將雅克比矩陣值用于積分轉換通過得到：

二維單元的剛度矩陣單元坐標系定義為：

由位移反推應變應力，即可推導出基于摩爾-庫倫屈服準則的剪應力為：

1.2 基床系數

基床系數是指假設中基底任隨意某點的壓強P與該點沉降量S的比值，通常也表示土體的剛度，基床系數Kv通過載荷試驗P～S曲線獲得，標準承壓板的變成一般設置為30cm，而利用非標準承壓板所做的載荷試驗獲取的基床系數需進行數值修正。

黏土：

砂土：

式中，b表示為非標準承壓板的邊長，Kv1為在非標準承壓板試驗下進行載荷試驗獲取的基床系數。

但是現場測試的弊端日益突出，并且測試誤差率一直較高，無法滿足現場測試需求，室內測試包括3軸和集成測試長時間內也沒有參考量，因此基床系數可靠性較低。地基的剛度大多數情況下也可用土體的壓縮模量進行表示，壓縮模量和基床系數呈正比關系，由于土體基床系數獲取的測量值與實際值之間存在很大差異，因此通常使用壓縮模量來表示土體剛度。在項目實施進程中，壓縮模量根體土體的壓縮模量的熟知而設定，從而對土體結構進行分析研究。

2 數值模擬

2.1 有限元模型建立

本文基于無為武威市地下復合管廊，對地下結構和周圍土體進行力學分析，將地下管廊埋在3m深處，復合管廊的大小和深度如圖2所示。

混凝土強度設置為C30，彈性模量設置為3×104MPa，重度設為25kN/m3，標準段綜合管廊寬為8.1m，兩外墻厚度取值為0.35m，兩內墻墻厚度取值為0.25m，3個艙室寬度從左到右分別為2.9m、2.3m、1.7m;標準段綜合管廊高設為4.15m，頂板厚設置為0.35m，底板厚度取值為0.4m，三艙室凈高設置為3.4m。綜合管廊的底部包括一個10cm的素混凝土防水墊，并從綜合管廊的兩個外壁延伸10cm。為了對建模過程進行簡化，省略了10cm素防水混凝土墊層。

將地下管道放置在半空間主體中，主體大小為20m×15m，將管道土體采用單一土體來代替，通過測量出的內摩擦角和層系數，將其進行對比，同時將接觸單元放置在混凝土以及土壤間，具體如圖3示意圖所示，表1為單一土體參數。

通常，簡單的直接剪切測試無法模擬結構與土壤之間的接觸行為，由于在現場操作中無法完成簡單加載或者循環加載，但是，結構與土體的接觸行為對現場直接剪切試驗、室內三軸試驗和固結試驗都適用。接觸效應可以通過在混凝土結構和土體間建立接觸單元來模擬。

利用Midas GTS NX軟件能夠完成對地下復合管廊中應力點實施不間斷監控，將7個監視點從左到右在其底部依次放置，監視點位于垂直墻和房間中央，圖4表示特定監視點所在坐標，圖5表示Midas GTS NX中的接觸式接口設備的顯示，圖6表示土體中半無限體模型。

2.2 壓縮模量數值分析

表1反應了Midas GTS NX中單個土體的壓縮模量變化，調整基床系數，控制土體參數不變，其內摩擦角φ為20°;當前，在地下廊道底面的接觸部位上，可以監視7個監察點的垂直力、切向力、垂直位移和切向位移，圖7表示垂直力的監控數據;圖8表示切向力的監控數據;圖9表示法向位移監控數據。

從圖7我們可以看出，在7個監視點中，編號1、7垂直力檢測數據大致相同，其余5個編號與之相比偏差較大，當土體壓縮量從7.45 MPa變為19.45 MPa時，上述監察點的垂直力值幾乎保持不變。然而編號1、2、6、7的垂直力呈下降狀態，編號3、4、5與之相反，呈上升狀態，隨著土體壓縮模量的上升，這加強了結構中間部分的剛度，并且得出了土體反應模量的上升對結構幾乎沒有影響。

如圖8所示。切向力可為正方向或者是負方向，在7個監視點中，編號2、3方向相同，其余5個編號相同?？梢钥闯龇ㄏ蛄Φ谋O測值要遠遠大于同一位置的切向力，并且當土體的壓縮模量從7.45 MPa上升到19.45 MPa時，由于每個位置處的切向力都在不斷減小，并且降低的值都不是很大，因此土體壓縮系數的上升表明對結構幾乎沒有影響。

在圖9中，7個監測點的位移是協同的，幾乎不變，并且當土體的壓縮模量從7.45MPa上升到19.45MPa時，法向位移標量發生了明顯變化，當其從7.45MPa上升到19.45MPa時，垂直位移總體減小了約29.5%;當其從10.45MPa上升到13.45MPa時，垂直位移總體減小了約22.0%，與壓縮量是7.45MPa時相比，減小了約44.6%;當其從13.45MPa上升到16.45MPa時，垂直位移總體減小了約18.2%，與壓縮量是7.45MPa時相比，減小了約54.6%;當其從16.45MPa上升到19.45MPa時，垂直位移總體減小了約15.4%，與壓縮量是7.45MPa時相比，減小了約61.8%;當其值由小變大時，垂直位移總體減少，而減小量則從大到小。地下管廊結構的剛度沒有發生改變，接觸界面的剛度與土體呈正比關系，土體抵抗結構沉降隨著土體的壓縮模量增大而變強。該分析表明，地下結構的沉降和變形受堅實的土體影響。

圖10表明，切向力可為正方向或者是負方向，每個監測點切向位移數值絕對值均比較小，當壓縮模量為7.45 MPa時，監測點2的切向位移相對于其它6個點數值偏大，這是因為此時土體壓縮模量過小，界面接觸剛度過低，2號監測點所在艙室空間大，整體剛度也較低，故而數值較大;當壓縮模量提高至10.45 MPa以后，2號監測點的數值穩定保持較低水平;總體而言，地下管道通道的觀察點的切向位移在土體的每個壓縮量都是一致的，切向位移的絕對值和地基土體的壓縮模量呈反比，土體的壓縮模量越大，切向位移越小，這在地下綜合管廊的總應力中發揮著積極作用。

2.3 內摩擦角數值分析

控制土體其他參數不變，通過改變表1中的內摩擦角來改變結構—土體界面摩擦角，土體的壓縮模量為7.45MPa;在地下廊道底面的接觸部位上，可以監視7個監察點的垂直力、切向力、垂直位移和切向位移，圖11表示法向力的監測數據;圖12表示切向力的監測數據;圖13表示垂直位移監測數據;圖14表示切向位移監測數據。

圖11表明，順著土體管道走廊的順序排列的7個觀測點的值顯示出中間比兩邊要大的態勢;當其內摩擦角從20°變成40°時，每個觀察點的垂直位移值不會發生顯著變化，具體而言，特定觀察點的值變化很小，改變土體摩擦角大小不會顯著影響下管道結構的垂直受力形態。

如圖12所示。切向力可為正方向或者是負方向，在7個監視點中，編號2、3方向相同，其余5個編號相同，可以看出法向力與切向力都很小，法向力值要遠遠大于同一位置監測的切向力，當其內摩擦角從20°變成40°時，每個觀察點的垂直位移值不會發生顯著變化，具體而言，特定觀察點的值變化很小，改變土體摩擦角大小不會顯著影響下管道結構的水平向受力形態。

圖13表明，順著土體管道走廊的順序排列的7個觀測點的值顯示出中間比兩邊要大的態勢，當土體的內摩擦角從20°增加到40°時，每個觀察點的法向位移值不會發生顯著變化，然而，總體態勢是，土體管道底部的應力條件隨著土壤的摩擦角越大變得越強，垂直位移值會略有降低，改變土體摩擦角大小不會顯著影響下管道結構的受力形態，但是增加內部摩擦角對地下結構的強度具有正面影響。

圖14表明順著土體管道走廊的順序排列的7個觀測點的垂直位移大體一致，可為正方向或者是負方向，切向位移的絕對值很小，當內摩擦角從20°增加到40°時，每個監視點的切向位移沒有明顯改變，因此土體內摩擦角對地下綜合管廊的切向位移影響可忽略不計。

3 結論

文章通過Midas GTS利用有限元模擬對地下管廊和土體間的接觸進行分析，將7個監視點管廊底部底部依次放置，監視點位于垂直墻和房間中央，分別監視垂直力、切向力、垂直位移和切向位移。利用有限元模型，改變土體內摩擦角及壓縮模量，對監測結果進行分析，從而得出其對接觸面受力的影響，結論如下：

（1）對4節點接觸單元進行描述，并給出了界面單元詳細的力學計算公式;通過對基床系數和土體壓縮模量的簡單描述，同時解釋了可用壓縮模量替代基床系數，并能保證其可靠性;采用地下復合管廊參數在Midas GTS NX中搭建有限元模型，并把7個監控點放置在地下綜合管廊與土體間的界面接觸單元中。

（2）土壤壓縮系數的變化對地下管道的機械狀態機會無影響，對切向位移的也影響甚微，法向位移則會受到明顯影響，法向位移和壓縮模量呈反比關系，增加土體壓縮模量對綜合管廊的強度和變形具有正面影響;改變土體摩擦角大小不會顯著影響綜合管廊的受力形態，對法向位移、切向位移數值的也影響甚微，但是增加內摩擦角對綜合管廊的強度和變形具有正面影響。

參考文獻

[1]于晨龍，張作慧.國內外城市地下綜合管廊的發展歷程及現狀[J].建設科技，2015，179（12）：49-51.

[2]C.S.Desai，B.K.Nagaraj. Modeling for cyclic normal and shear behavior of interfaces[J]. Journal of Engineering Mechanics，1988，114（7）：1198-1217.

[3]Brandt J R T.Behavior of soil-concrete interfaces[M].Canda： The University of Alberta， 1985.

[4]K.G.Sharma，C.S.Desai.Analysis and implementation of thin-layer element for interfaces and joints[J].Journal of Engineering Mechanics，1992，118（12）：2442-2462.

[5]盧廷浩，鮑伏波.接觸面單元薄層耦合本構模型[J].水利學報，2000，02：71-75.

[6]唐世斌，唐春安，梁正召.接觸面損傷演化過程的數值模擬[J].計算力學學報，2011，28（1）： 146-151.

[7]鄧建，肖明，謝冰冰.三維接觸面單元彈塑性損傷數值模擬分析方法[J].中南大學學報（自然科學版）， 2016，47（7）：2383-2389.

[8]Fredlund D G，Rahadjo H.Soil mechanics for unsaturated soil mechanics[M].New York：John Wiley and Sons，1993.

[9]Evgin E，Fakharian K.Effect of stress paths on the behavior of sand-steel interfaces[J]. Canadian Geotechnical Journal，1996，33（33）：853-865.

[10]Uesugi M，Kishida H.Influential factors of between steel and dry sands[J].Soils and Foundations，1986，26（2）：33-46.

[11]胡黎明，濮家騮.土與結構物接觸面損傷本構模型[J].巖土力學，2002，23（1）：6-11.