1、摘要：對埋地管道結構分析，考慮管道與土的相互作用問題是非常必要的。其相互作用問題包括管土界面滑移、分離現象及界面處接觸應力的確定。采用非線性接觸的處理方法建立了管土相互作用分析模型，管道及其周圍土體分別以空間薄殼單元和實體單元進行離散，采用有限元手段分析了斷層錯動下管土界面的非線性反應情況。通過算例，分析了位錯量、埋深、管道徑厚比及場地條件對管土界面非線性反應的影響；此外對管土界面接觸應力進行了分析，發現近斷層處管土發生強烈的非線性反應，在距離斷層較遠處，管土切向摩擦應力的是以一種“拋物線”式的方式遞減，而非以線性遞減的方式減少。關鍵詞：埋地管道，非線性接觸，斷層運動，接觸應力 1.引言 埋地

2、管線是重要的生命線工程，影響人們日常的生產生活。西氣東輸和川氣東送同三峽工程、青藏鐵路、南水北調并稱為中國五大工程，建成后將形成覆蓋全國的天然氣管道網絡。這些重大的天然氣管道工程遍及全國各地的自來水、供燃氣、排水等管道工程，這些管道工程構成了維系現代社會正常運轉的基石。管線工程多埋于地下，形成的網絡系統龐大復雜，需要覆蓋廣闊的地域，對于一些不利甚至危險的地域也必須經過。例如，西氣東輸工程跨越了多個地震活動地區，需要穿越100多個斷層。埋地管線擔負著輸送生產和生活物質的功能，它需要覆蓋廣闊的區域，遇到不利甚至危險的地質條件(斷層、液化、震陷、滑坡等)也無法避讓而只能穿過。這些地質災害將導致地面產

3、生永久變形，嚴重威脅著埋地管線工程的安全性，有可能造成管線系統服務功能的下降甚至喪失，對社會生產、生活產生影響。更為嚴重的是，油、氣、水等泄露不僅經濟損失巨大，而且可能導致環境污染，甚至引發火災等次生災害。因此，研究在地面永久變形下埋地管道的工作性狀，是不得不面臨和需要解決的重要問題之一。埋地管道在上覆土層和外荷載等作用下，管道將因受力而變形，管周土體既是作用在管道上的荷載，同時又是增強管道強度和剛度的一種介質。因此在研究地埋管道的工作機理時，必須把管道周圍一定范圍內的土體作為結構的一部分加以考慮，即考慮管土相互作用問題。地面永久變形下管土相互作用問題主要包括管土接觸界面的相對滑移、分離及閉合

4、等非線性反應現象以及界面處分布接觸應力的確定1-3。地面大變形下埋地管道的反應分析研究從最初的索及梁模型發展到殼模型，是一個巨大的改進；但是目前處理管土相互作用時，多將土體對管道的作用簡化為一系列彈簧參數4-6來替代，而土體的強非線性也導致了彈簧參數的難確定。本文采用非線性接觸7-12的處理方法建立了管土相互作用模型，通過有限元計算分析了位錯量、埋深、管道徑厚比及場地條件對管土界面非線性反應的影響，并對管土界面上接觸應力及反應變化規律進行了簡單的分析研究。2.有限元模型及參數選擇 大量的動力和靜力實驗及研究分析表明管線一般情況下受震和受到斷層作用運動形式和原因是不同的，跨越斷層埋地管線運動主要

5、是由于地表的永久性位移產生，其動力放大作用微小，可以忽略管道慣性力的影響，因此對管線的地震反應可以采用擬靜力方法進行研究。建立了地面永久變形下管土相互作用的三維有限元分析模型，如圖1所示。管道及周圍土體從半無限地球介質中取出，采用非線性接觸單元來模擬管土之間的相對滑移、分離及閉合等非線性現象。由模型可知，計算區域包括管道結構及周圍土體介質，由6面邊界條件組成，除了地表為自由面外，其他5個面為人工邊界條件。邊界條件為：地表為自由面，左右端面、模型前后兩面為水平約束，左側底面采用豎向約束，右側底面為位移輸入邊界；管道兩側為自由邊界。圓柱形中空薄壁管道采用四節點薄殼單元，土體采用八節點三維實體單元，

6、考慮初始應力場的影響，采用線性位移加載方式來表征斷層位錯作用，對管土相互作用體系進行非線性分析，研究管土交界面上非線性接觸的反應過程。取管徑為0.4m的鋼質管道作為案例分析對象，上、下盤土體縱向長度均取為160m，即為管徑的400倍長度；土體寬度取為8m（管徑的20倍長度）；管道計算長度取為324m。在斷層兩側各20m區域內網格進行細分。模型中鋼質管道采用三折線本構模型，其彈性階段彈性模量MPa，當應力達到418MPa管道開始進入彈塑性階段，當應力與塑性應變分別達到516MPa、0.016時完全進入塑性狀態，泊松比52.110E=0.3=。土體采用Mohr-Coulomb模型。 3.管土界面反

7、應現象分析 為了了解認識地震導致的地面永久變形下管土交界面上非線性接觸反應現象，主要包括管土界面的滑移及管土相對分離現象，對正斷層作用下管土相互作用反應進行了參數影響分析。以接觸面下側的豎向位移為指標來衡量管土間的相對分離反應，以接觸面上側的軸向位移為指標來表征管土間的相對滑移行為。依據模型計算結果，從土體位錯量、埋設深度、管子徑厚比及場地條件等方面來分析其對管土相互作用反應的影響。不同位錯情況下埋地管道與周圍土體接觸界面上產生的相對分離及相對滑移的分布中可以看出隨著位錯量的增加，斷層處管道下側豎向位移隨之增大，距斷層越近分離越明顯，且管土分離區域也慢慢隨位錯增大而增大；當位錯達到0.4m以后

8、，管土分離區域基本不變，分離區域約6m(12D)，這里管道發生明顯的橢圓化變形甚至“鉸式”屈曲。管土界面上側的相對滑移的反應情況隨著位錯量增加管子靠近斷層處的軸向位移隨之增大，且管子與土體產生相對滑動的區域不斷增加，在距離斷層一定遠處，管土之間近乎無相對滑移，可以將其看成是“錨固點”，給有效計算長度的選取提供了一個理論依據。同一加載水平下不同埋深對管土相對分離及相對滑移行為的影響中不難看出，隨著埋設深度的增大，管道豎向位移峰值減小，且管土分離的區域也隨之減小，但是由于相同管徑及相同位錯的緣故，管土分離的區域近似限制在距離斷層約6m處。可以知道管道埋設深度越淺，管土間相對存在滑移的區域越大，從而

9、在地震斷層作用下能夠更多的吸收地震能量，即管道抵抗斷層的能力與埋深成比，淺埋時管道易于變形，益于抗震。同一加載水平下管道徑厚比不同時導致的接觸界面上管土之間相對分離及滑移的變化中可以看出管土分離區域隨著徑厚比的增大而呈減小趨勢，但是在斷層處管道的豎向位移峰值相近；圖4(b)是界面上側管土相對滑移情況圖，明顯的是隨著徑厚比的增大，存在相對滑移的區域在減小，這可以從側面說明當管徑大小不變時，隨壁厚增加，滑移區域越大，管道在土中越易于活動變形，從而利于抗震。土質材料不同情況時對管土相對分離及相對滑移的影響看出，隨著土體硬度增加，管土相對分離區域以及斷層處豎向位移峰值隨之減小；從中發現管土界面存在相對

10、滑移的區域隨著土體硬度的增加而減小，同時管道上側軸向位移峰值也隨土體變硬而減小。因此，當管道埋設在松散的或中等密實的土體中時，管道在斷層作用下容易變形，能耗散較大的地震能力。4.接觸界面接觸應力分析 斷層左側固定部分管土接觸面斷層處接觸應力隨斷層位錯變化而變化的反應過程中可以看出管土接觸界面上側節點249及周圍臨近節點248、250先隨著位錯增大而接觸壓力和摩擦應力隨之劇增，然后慢慢減小，當位錯量達到0.4m時節點248及250與管道發生分離，管土之間轉化為自由邊界，沒有力的傳遞。對于其它的下側的節點，由于斷層錯動的產生，很早就產生管土的分離現象，不進行力學傳遞，成為自由邊界。此外，不難發現，

11、斷層附近處管土非線性接觸現象非常強烈，管道與土體之間隨著斷層錯動量的增加不斷的改變其接觸狀態。距離斷層界面80m處(節點611-620)管土接觸面上接觸應力隨著斷層位錯變化而變化的過程中可以看出，接觸界面上各節點的接觸壓力在斷層位錯量增大時近乎保持恒值，這是因為距離斷層較遠處管道與周圍土體之間沒有發生分離現象，其接觸壓力值為初始應力狀態下管土接觸界面上的相互作用值。而從中可以看出隨著斷層錯動量的增加，管土界面上軸向摩擦應力值隨之呈現“拋物線”式的增大趨勢，所以認為遠離斷層截面處的接觸界面上摩擦應力隨位錯量增大呈線性增大的說法是不正確的；也可以看出管土接觸面上兩側的軸向摩擦力呈對稱分布。當然，從

12、圖中的曲線變化來看，距離斷層較遠時管土相互作用相對斷層附近處相互作用要弱很多。斷層左側固定部分距離斷層30m-160m段內管土界面上側節點摩擦應力隨位錯量的變化關系中橫向坐標“0”代表的是左側土體的端部界面上節點，坐標“130”代表的是距離斷層30m處管土接觸面上側土體節點位置。之所以取距離斷層處30m之外的界面上的摩擦應力進行分析，是因為30m外已經是大變形之外了，這段范圍內管道及周圍土體應力和變形基本處于彈性范圍之內。從中我們發現，隨著位錯量的增加管土界面切向摩擦應力隨之增大。此外，從中還可以看出切向摩擦應力的是以一種“拋物線”式的遞減方式，逐漸減少，到距離端部40m處變為零。所以假設管土

13、系統反應大變形之外的管土界面摩擦應力以一種線性方式而減少的論斷是不合理的。5.結語 采用非線性接觸問題的處理方法建立了管土相互作用分析模型，研究了交界面上管土相對分離、滑移的非線性接觸行為，對影響管土接觸界面反應的幾種因素作了分析。通過參數分析，得到如下結論：(1)位錯量增加時，管土相對分離區域及管道下側豎向位移隨之增大，但由于管徑因素，其相對分離的區域受到限制；同時其相對滑移區域及管道的軸向滑移量也隨位錯增大而增大。(2)埋深越大，其相對分離區域及管子豎向位移越小，且存在相對滑移的區域也隨埋深增大而減小；埋深越大，管道不易在土體中活動，對抗震不利。(3)管道徑厚比越大管土相對分離區域越小，存在相對滑移的區域越小；管徑大小不變時，隨壁厚增加，滑移區域越大，管道在土中越易于活動變形，耗散更多地震能量，利于抗震。(4)從軟土到硬土，管土相對分離區域