賈曉輝 劉愛文

?

跨斷層水泥管試驗(yàn)的有限元分析

賈曉輝 劉愛文

（中國(guó)地震局地球物理研究所，北京 100081）

在考慮管道的材料非線性和幾何非線性、管土相互作用的非線性和管道接口非線性的基礎(chǔ)上，建立了由管體梁?jiǎn)卧?、三向土彈簧單元和接口單元組成的埋地非連續(xù)管道在斷層位移作用下的有限元模型，并以美國(guó)密歇根大學(xué)Junhee等（2010）所做的跨斷層水泥管試驗(yàn)為原型進(jìn)行了模擬分析。有限元結(jié)果給出的水泥管最終變形、接口轉(zhuǎn)角、接口位移與實(shí)驗(yàn)結(jié)果基本一致，表明本文提出的跨斷層埋地非連續(xù)管道抗震計(jì)算的有限元分析方法具有一定的合理性。有限元結(jié)果和試驗(yàn)結(jié)果都表明，在逆沖斷層作用下，水泥管的破壞主要是因?yàn)樵诠艿澜涌谔幍妮S向壓力和彎矩的耦合作用，在斷層附近的管道接口承受了較大的轉(zhuǎn)動(dòng)和壓縮位移。本文所提出的分析方法可推廣到埋地非連續(xù)管道在其它永久地面變形作用下的有限元分析。

埋地管線 水泥管 斷層作用 有限元方法 數(shù)值分析

引言

埋地管道承擔(dān)著給水排水、輸油、供氣、供熱等服務(wù)功能，是現(xiàn)代城市生存和發(fā)展的命脈，被形象地稱為地下生命線工程。城鎮(zhèn)化程度越高，生產(chǎn)生活活動(dòng)對(duì)其依賴性越強(qiáng)。埋地管道一旦遭受破壞，不僅造成埋地管道本身結(jié)構(gòu)的毀壞和其承擔(dān)的服務(wù)功能的喪失，而且會(huì)更進(jìn)一步引發(fā)水、火等次生災(zāi)害。油氣輸送管道一般采用焊接和熱熔連接的連續(xù)鋼管和PE管，其整體性能較好；城市給排水管道一般使用鑄鐵管、混凝土管、鋼管、塑料管等，接口多采用承插連接或法蘭連接。

多次地震經(jīng)驗(yàn)表明，對(duì)埋地管道影響最大的是地震斷層，其次是地基液化，然后才是地震波。埋地管道通常要穿過廣大地區(qū)的線狀工程，考慮抗震設(shè)計(jì)時(shí)地面的永久變形往往比地震動(dòng)的峰值加速度更加重要。因此，埋地管道跨越錯(cuò)動(dòng)斷層的性能研究，一直是生命線工程研究領(lǐng)域的一項(xiàng)重要課題。目前，國(guó)內(nèi)外研究者在跨越斷層埋地連續(xù)管道的地震反應(yīng)方面做了大量的研究工作，取得了較多的研究成果。在理論解析方法方面，有經(jīng)典的解析方法（Newmark等，1975；Kennedy等，1977；Wang等，1985）以及其他學(xué)者提出的解析方法（張素靈，1999；劉愛文等，2002；王濱等，2011）。這些方法計(jì)算簡(jiǎn)便，易于工程應(yīng)用，但是無(wú)法準(zhǔn)確分析管道在斷層作用下，因?yàn)槭芗?、彎、拉、壓的綜合作用而發(fā)生的強(qiáng)度破壞、局部屈曲或整體失穩(wěn)。在數(shù)值模擬分析方面，國(guó)內(nèi)外學(xué)者以有限單元法為基礎(chǔ)，建立了跨斷層埋地連續(xù)管道的有限元計(jì)算模型，有梁?jiǎn)卧Ｐ停ü鳁澋龋?999）和殼單元模型（Takada，2001；劉愛文等，2005）等，有限元計(jì)算模型可以比較真實(shí)地模擬地震斷層作用下埋地連續(xù)管道的破壞過程。

以上研究主要針對(duì)埋地連續(xù)管道，而對(duì)于帶接口的埋地非連續(xù)管道的研究相對(duì)較少。原因之一是連續(xù)管道一般為長(zhǎng)距離輸油供氣的高壓能源管道，非連續(xù)管道多為用于給排水的市政管道，連續(xù)管道承擔(dān)的任務(wù)有很強(qiáng)的戰(zhàn)略意義；原因之二是連續(xù)管道的抗震研究不必考慮接口的破壞，而接口破壞機(jī)理的研究正是非連續(xù)管道抗震研究的關(guān)鍵問題。在管道的抗震設(shè)計(jì)中，對(duì)于埋地連續(xù)管道一般是通過計(jì)算管道的應(yīng)變或應(yīng)力來進(jìn)行抗震驗(yàn)算的；埋地非連續(xù)管道則與前者不同，一般是通過計(jì)算接口的位移變化和轉(zhuǎn)角變化來進(jìn)行抗震驗(yàn)算。近年來，國(guó)內(nèi)外學(xué)者對(duì)埋地非連續(xù)管道的研究方興未艾，主要是試驗(yàn)研究，如：對(duì)承插式埋地球墨鑄鐵供水管道進(jìn)行抗拉靜載試驗(yàn)（周靜海等，2009；陳春光等，2012），旨在分析管道接口的力學(xué)性能參數(shù)；美國(guó)密歇根大學(xué)進(jìn)行的跨斷層水泥管試驗(yàn)的研究（Junhee等，2010），旨在探討非連續(xù)管道在斷層作用下的破壞機(jī)制。

本文以Junhee等（2010）進(jìn)行的跨斷層水泥管原型試驗(yàn)研究為基礎(chǔ)，考慮了管土相互作用的非線性和管道接口的非線性，建立試驗(yàn)?zāi)Ｐ退鶎?duì)應(yīng)的有限元模型，模擬了水泥管在逆沖斷層作用下的破壞過程。

1 跨斷層水泥管試驗(yàn)

為研究分段水泥管在斷層作用下的變形行為，Junhee等（2010）在專業(yè)用于測(cè)試管道在永久地面變形作用的測(cè)試臺(tái)上，建立了由5根標(biāo)準(zhǔn)長(zhǎng)度和1截部分管段裝配而成的管道試驗(yàn)?zāi)Ｐ?，管道接口為承插口形式，接口處采用橡膠圈和水泥漿密封。位移傳感器密集排列在管道的相應(yīng)位置，實(shí)時(shí)提供管道在斷層作用下的位移響應(yīng)，并且分段水泥管接口部分的接口位移變化和接口轉(zhuǎn)角變化同時(shí)可以通過測(cè)量得到。測(cè)試臺(tái)中管道下伏土層的厚度為20cm，管道上覆土層的厚度為140cm。管段編號(hào)自左向右依次為Seg #1—Seg#5，接口編號(hào)自左向右依次為Joint #1—Joint #5。管道兩端采用鋼板鎖定的方式將管道的兩個(gè)端部邊界固定在測(cè)試臺(tái)上。斷層通過管道Seg #3的中心，并與管道呈50°交角，沿?cái)鄬泳€發(fā)生的斷層位移量最終為30.5cm，管道中心以南固定，管道中心以北的管段整體施加位移荷載，分12步加載完畢（圖1）。試驗(yàn)中所選用的水泥管試件的規(guī)格如表1所示。

表1 水泥管試件規(guī)格

2 有限元模型的建立

本文建立的跨斷層水泥管試驗(yàn)的有限元模型如圖2、圖3所示。其中，圖2為有限元模型示意圖；圖3為用有限元軟件建立的有限元模型。

埋地管道在斷層大位移作用下將進(jìn)入非線性反應(yīng)狀態(tài)，埋地管體本身的模擬涉及到材料非線性和幾何非線性的雙重非線性問題，這是模型建立的關(guān)鍵問題之一。管道周圍的土體在斷層作用下具有復(fù)雜的非彈性特征，土體在受斷層作用后是如何作用于管道并使管道發(fā)生變形的，以及如何實(shí)現(xiàn)合理的管土相互作用是模型建立的關(guān)鍵問題之二。震害調(diào)查表明，管道接口或其附近往往發(fā)生最頻繁和最嚴(yán)重的破壞，如何模擬斷層大位移作用下接口的非線性是模型建立的關(guān)鍵問題之三。

2.1 管體梁?jiǎn)卧?/p>

埋地非連續(xù)管道與連續(xù)管道的破壞不同，后者在斷層作用下常見屈曲破壞，所以多選用殼單元來實(shí)現(xiàn)殼體的屈曲模擬。前者在斷層作用下常見管道接口破壞，管體本身的破壞較少見，所以本文中管體選取能用于非線性分析的梁?jiǎn)卧M(jìn)行模擬。管體本身的材料非線性通過選用非線性本構(gòu)模型來實(shí)現(xiàn)，因?yàn)楣荏w的材料為混凝土，所以本文選取適用于混凝土的Drucker-Prager材料，材料常數(shù)見表2。幾何非線性通過選用有限元軟件中的非線性方程解法選項(xiàng)來實(shí)現(xiàn)。

表2 混凝土的參數(shù)

2.2 土彈簧單元

考慮土體的非線性，管土之間的相互作用通過三個(gè)定向非線性彈簧來模擬。管道每個(gè)梁?jiǎn)卧墓?jié)點(diǎn)上分別連接三個(gè)方向的定向土彈簧：管軸方向、水平方向和垂直方向。這三個(gè)彈簧分別用來考慮管軸方向的土摩擦力、水平方向及垂直方向的土壓力。在每個(gè)土彈簧的節(jié)點(diǎn)上施加斷層相應(yīng)方向的位移量，從而模擬管道所受的斷層作用。因此，在本文跨斷層水泥管試驗(yàn)的有限元模型中，斷層的位移是通過施加到土彈簧相應(yīng)的節(jié)點(diǎn)上來對(duì)水泥管進(jìn)行作用的。

試驗(yàn)所選用的土體參數(shù)見表3。土彈簧的力-位移關(guān)系為雙折線模型，如圖4所示。土彈簧的非線性模型中參數(shù)的計(jì)算參考《油氣輸送管道線路工程抗震技術(shù)規(guī)范》（中國(guó)石油天然氣集團(tuán)公司，2009）得到，所求參數(shù)結(jié)果見表4。

表3 土體介質(zhì)參數(shù)

表4 三個(gè)方向土彈簧參數(shù)

2.3 接口單元

埋地非連續(xù)管道的接口破壞是管道破壞的主要形式，在斷層作用下接口將受到拉伸、壓縮、彎曲、剪切等作用，相應(yīng)產(chǎn)生接口的拔脫、壓裂、轉(zhuǎn)動(dòng)、剪斷等破壞。依據(jù)接口的常見破壞現(xiàn)象和接口試驗(yàn)研究結(jié)果，接口采用軸向彈簧、彎曲彈簧和橫向彈簧來模擬，三個(gè)方向的彈簧分別用來模擬接口的伸縮、接口的彎曲和接口的剪切三個(gè)特性。

接口三個(gè)方向彈簧的廣義力-位移關(guān)系如圖5所示。其中，軸向彈簧的力-位移關(guān)系和彎曲彈簧的力矩-轉(zhuǎn)角關(guān)系均為雙折線模型；橫向彈簧的力-位移關(guān)系為線性模型，并且橫向采取大剛度進(jìn)行約束。圖5模型中參數(shù)的數(shù)值根據(jù)管道的接口試驗(yàn)所得，本文中所采取的參數(shù)的數(shù)值見表5。

這里還要特別提到，接口在壓縮力作用下，隨著接口壓縮位移量的增加，軸向壓力呈線性增長(zhǎng)變化，與拉伸情形下不同，最大作用力不能簡(jiǎn)化為水平線，如圖5中a所示，在壓縮區(qū)設(shè)定了彈簧剛度系數(shù)為的軸向壓縮彈簧。

表5 三個(gè)方向接口彈簧參數(shù)

2.4 邊界條件

在跨斷層水泥管試驗(yàn)中，水泥管的位移加載如圖1所示，以管道Seg #3的中心為界，中心以南固定，中心以北施加位移荷載。管道兩端在試驗(yàn)時(shí)均采取鋼板鎖定的方式，中心以北的管道（含管道北端）是隨著測(cè)試臺(tái)整體移動(dòng)的，因此本文中有限元模型的管道南端設(shè)置成完全固定邊界（即3個(gè)位移分量和3個(gè)轉(zhuǎn)動(dòng)分量均設(shè)定為零值），管道北端設(shè)置成只固定轉(zhuǎn)動(dòng)自由度的固定邊界（即3個(gè)轉(zhuǎn)動(dòng)分量設(shè)定為零值，3個(gè)位移分量與斷層的位錯(cuò)量一致）。

3 有限元結(jié)果分析

圖6（a）給出了跨斷層水泥管原型試驗(yàn)的最終管道變形，（b）和（c）給出了各接口最后的實(shí)際破壞狀態(tài)。

圖6 （a）測(cè)試后的管道變形；（b）斷層附近接口2、接口3在轉(zhuǎn)動(dòng)和壓裂作用下的綜合破壞；（c）在接口1、接口4、接口5處的壓裂破壞

Fig. 6 （a）The final shape of pipeline after the test;（b）Combined failure by rotation and contraction at two joints（#2 and #3）near the ground fault;（c）Compressive telescoping failures at the other joints（#1，#4，and #5）