作者：2017級本科生 劉巍、戰(zhàn)朋禹、汪子涵 土木系
指導(dǎo)老師：林佳瑞 土木系
關(guān)鍵詞：有限元分析 粘彈性人工邊界 接觸面對法 地震響應(yīng)
摘要
地下管道網(wǎng)絡(luò)是城市基礎(chǔ)設(shè)施最重要的環(huán)節(jié)之一。在地震中，地下管道的直接破壞造成了巨大的公共財產(chǎn)損失，管道災(zāi)害等一系列次生災(zāi)害的伴生將損失進(jìn)一步擴(kuò)大。本作品旨在將地下管道的地震易損性因素引入管道設(shè)計軟件中，為地下管道設(shè)計提供可靠的分析工具。
設(shè)計插件207.png)
圖1 基于粘彈性人工邊界的地下管道地震動力學(xué)響應(yīng)項目流程簡圖
1 粘彈性人工邊界問題的參數(shù)選定
粘彈性人工邊界是針對無限區(qū)域或者半無限區(qū)域的有限元模擬中，如何從無限區(qū)域中取出一個有限區(qū)域并加之一個人工邊界條件，以盡可能真實準(zhǔn)確地反應(yīng)無限區(qū)域介質(zhì)的瞬態(tài)波動問題。我們采取的粘彈性人工邊界為彈簧-粘壺并聯(lián)約束邊界。
設(shè)計插件366.png)
圖2 粘彈性邊界條件
粘彈性人工邊界條件中的重要參數(shù)是彈簧的剛度K。根據(jù)推演出的公式，彈簧的剛度K與波源到這一點的距離R成反比。目前主流研究中采用的參數(shù)R為Rmin。但選取Rmin時，彈簧剛度K是這個平面上的最大值，這使得結(jié)構(gòu)的位移減小，使得分析結(jié)果低于實際值，存在低估管道地震下位移的風(fēng)險。
設(shè)計插件515.png)
圖3 不同R值下粘彈性邊界的位移時程
通過調(diào)整參數(shù)選擇方法，我們發(fā)現(xiàn)在設(shè)計插件576.png)
的物理參數(shù)下精度上升更加接近實際解，但是仍然有低估結(jié)構(gòu)位移量的問題；也發(fā)現(xiàn)當(dāng)選取設(shè)計插件642.png)
時模擬結(jié)果由低估轉(zhuǎn)化為高估結(jié)構(gòu)位移量，應(yīng)用此時模擬結(jié)果偏向于安全，但是存在精度較低的問題。因此在設(shè)計中采用，形成的區(qū)間作為參考值可以在保證設(shè)計安全性的同時提高精度。
2 動接觸問題
地下管道與土層的接觸面是一個動接觸問題。對于動接觸問題，我們采用接觸面對法，其相比于整體法和接觸單元法而言更加適用于強(qiáng)地震不連續(xù)約束現(xiàn)象。
接觸面對法有效地優(yōu)化了大變形接觸面的模擬，比如強(qiáng)地震作用時，地下結(jié)構(gòu)與周圍土體接觸面上的力學(xué)性態(tài)變動劇烈，常常存在開裂閉合等不連續(xù)約束現(xiàn)象，這將導(dǎo)致接觸在數(shù)值計算中往往不收斂。定義接觸面對可以較好地解決這種情形下的問題。
設(shè)計插件968.png)
圖4 接觸面對示意圖
3 管道地震響應(yīng)模擬
基于上述研究成果，在建模中采用接觸面對法作為動接觸面的邊界，通過采用粘彈性人工邊界的參數(shù)等，保證結(jié)果偏向于安全。
設(shè)計插件1050.png)
圖5 形變前Abaqus管道模型
對上述模型加上地震波進(jìn)行模擬，從變形前后對比的應(yīng)力紋云圖可以看出土層有明顯的下沉現(xiàn)象，并且在震源所在的豎直方向上下沉更加顯著，且在管道與土層表面之間的土體相較于其他土體而言擠壓更明顯。
設(shè)計插件1163.png)
圖6 變形前后對比Mises應(yīng)力云紋圖
從局部管道變形后應(yīng)力云紋圖中可以看出，節(jié)點相較于其他位置更容易受損，在管道節(jié)點處以及管道盡頭處的應(yīng)力更為集中，變形更加明顯，并且與土體脫離現(xiàn)象更加顯著；對于豎直方向放置的管道而言，其與土體的脫離現(xiàn)象相比于水平放置更加明顯；對于水平放置的管道其受擠壓而產(chǎn)生的變形更顯著。
設(shè)計插件1324.png)
圖7 局部管道變形后Mises應(yīng)力云紋圖
作者聯(lián)系方式：
清華大學(xué)土木工程系 劉巍 （北京 100084）
手機(jī)：18810730708
微信：lw729394910
歡迎進(jìn)一步交流指正！
