劉孝泰，顏慶智*，饒江

(1.中國(guó)石油大學(xué)(華東)儲(chǔ)運(yùn)與建筑工程學(xué)院，山東 青島 266580；2.勝利石油管理局供水公司，山東 東營(yíng) 257000)

黃島某地區(qū)綜合管廊應(yīng)力與變形的數(shù)值模擬

劉孝泰1，顏慶智1*，饒江2

(1.中國(guó)石油大學(xué)(華東)儲(chǔ)運(yùn)與建筑工程學(xué)院，山東 青島 266580；2.勝利石油管理局供水公司，山東 東營(yíng) 257000)

以黃島某段綜合管廊為分析對(duì)象，利用有限元軟件ANSYS建立其三維有限元模型，綜合考慮黃島該地區(qū)土體性質(zhì)、管廊施工方式、路面行人載荷和行車(chē)載荷，對(duì)管廊在正常使用階段的應(yīng)力分布和位移變化進(jìn)行數(shù)值模擬計(jì)算。分析結(jié)果表明，該段綜合管廊設(shè)計(jì)比較保守，材料沒(méi)有充分利用，還有進(jìn)一步優(yōu)化的空間；各艙室腋角處存在明顯的應(yīng)力集中現(xiàn)象，應(yīng)在結(jié)構(gòu)設(shè)計(jì)時(shí)予以重視，采取措施如增加腋筋配筋率，改善其受力狀態(tài)。所得結(jié)論可為今后類(lèi)似工程設(shè)計(jì)提供參考。

綜合管廊；應(yīng)力分布；位移；數(shù)值模擬

0 引言

綜合管廊是指城市地下用于集中敷設(shè)電力、通信、廣播電視、給水、熱力和燃?xì)獾仁姓芫€(xiàn)的公共隧道。推進(jìn)城市地下綜合管廊建設(shè)，統(tǒng)籌各類(lèi)市政管線(xiàn)規(guī)劃、建設(shè)和管理，解決反復(fù)開(kāi)挖路面、架空線(xiàn)網(wǎng)密集、管線(xiàn)事故頻發(fā)等問(wèn)題[1-8]。

黃島地區(qū)特殊性巖土主要是素填土、軟弱土和風(fēng)化基巖?？辈旖Y(jié)果表明，擬建場(chǎng)地及其影響范圍內(nèi)的周邊環(huán)境內(nèi)未發(fā)現(xiàn)影響場(chǎng)地穩(wěn)定性的巖溶、滑坡、危巖和崩塌、泥石流、采空區(qū)、地面沉降等不良地質(zhì)作用，不良地質(zhì)作用不發(fā)育。文中綜合管廊全長(zhǎng)12.8 km，管廊標(biāo)準(zhǔn)斷面采用4艙通行管溝結(jié)構(gòu)(分為污水艙、水力艙、電力艙和天然氣艙)，設(shè)計(jì)結(jié)構(gòu)標(biāo)準(zhǔn)斷面內(nèi)尺寸為13.85 m×4.6 m，如圖1所示。利用有限元軟件ANSYS建立綜合管廊的三維有限元模型，分析綜合管廊在正常使用階段的應(yīng)力與位移特點(diǎn)。

1 工況分析

1.1 材料參數(shù)

計(jì)算所涉及的材料主要是巖土、鋼筋、混凝土，鋼筋材料型號(hào)為HRB400、HPB300，混凝土強(qiáng)度等級(jí)為C45，混凝土材料參數(shù)取值見(jiàn)表2，鋼筋材料參數(shù)取值見(jiàn)表3，該段開(kāi)挖深度為8 m，土斷面示意圖如圖2所示。

土體參數(shù)根據(jù)地勘報(bào)告取值見(jiàn)表1，地下水位取地下0.5 m。

表1 土體參數(shù)Tab.1 Parameters of soil

圖1 綜合管廊橫斷面Fig.1 Cross section of utility tunnel

圖2 土層橫斷面Fig.2 Cross section of soil

C45混凝土容重/(kN·m-3)25彈性模量E/kPa3.35×107泊松比0.2抗拉強(qiáng)度設(shè)計(jì)值ft/kPa1.80×103抗壓強(qiáng)度設(shè)計(jì)值fc/kPa2.11×104

表3 鋼筋材料參數(shù)Tab.3 Parameters of reinforcement

1.2 荷載組合

荷載組合按照《建筑結(jié)構(gòu)荷載規(guī)范》GB50009-2012的相關(guān)規(guī)定計(jì)算。

γ0S=γ0(γGSGK+ψγQSQK)。

(1)

式中：γ0為結(jié)構(gòu)重要性系數(shù)，因?yàn)樵摻Y(jié)構(gòu)設(shè)計(jì)使用年限為100 a，故該系數(shù)取為1.1；S為承載能力極限狀態(tài)的荷載效應(yīng)(內(nèi)力)組合值；γG為永久荷載的分項(xiàng)系數(shù)，土壓力取1.27(有利時(shí)取1.0)，結(jié)構(gòu)自重取1.2(有利時(shí)取1.0)，土側(cè)壓力取1.4，水壓力取1.4；SGK為永久荷載標(biāo)準(zhǔn)值的效應(yīng)值，這里永久荷載包括土重力、土側(cè)壓力、結(jié)構(gòu)自重；ψ為可變荷載(人+車(chē)輛荷載)的組合系數(shù)，取為1.0；γQ為可變荷載(人+車(chē)輛荷載)作用效應(yīng)的分項(xiàng)系數(shù)，取為1.4；SQK為人+車(chē)輛荷載換算土重標(biāo)準(zhǔn)值的效應(yīng)[9-11]。

1.3 荷載的計(jì)算

(1)覆土壓力

p=h·γ。

(2)

式中：h為覆土厚度，m；γ為覆土重度，kN·m-3，地下水位以下的土層采用有效重度。

(2)結(jié)構(gòu)頂板自重

Q=25c。

(3)

式中：c為頂板厚度，混凝土重度25 kN/m3。

(3)土側(cè)壓力采用水土分算模型

結(jié)構(gòu)頂面土側(cè)壓：

q1=[γs+10(h-s)]K。

(4)

結(jié)構(gòu)底面土側(cè)壓力：

q2=[γs+10(H+h-s)]K。

(5)

結(jié)構(gòu)頂面水側(cè)壓力：

q3=10(h-s)]。

(6)

結(jié)構(gòu)頂面活載側(cè)壓力：

q4=P·K。

(7)

結(jié)構(gòu)底面處水側(cè)壓力：

q5=10(H+h-s)]。

(8)

式中：s為地下水位距離地面高度；h為頂板上覆土厚度為；H為結(jié)構(gòu)總高度；P為結(jié)構(gòu)頂面活載，該荷載中的行人荷載一般按4.0 kPa計(jì)算，行車(chē)荷載一般按20 kPa計(jì)算；K為靜止土壓力系數(shù)；K=1-sinφ(由表1取值)，采用水土分算模型，水的重度取10 kN/m3，在計(jì)算側(cè)壁活荷載時(shí)，汽車(chē)荷載不考慮動(dòng)力系數(shù)。

(4)水壓力

水壓力作用于結(jié)構(gòu)表面，計(jì)算公式如下：

pω=γω·hω。

(9)

式中：γω為水的重度，取10 kN/m3；hω為計(jì)算點(diǎn)距離地下水位面的深度，m[12-15]。

2 數(shù)值模擬

2.1 建立結(jié)構(gòu)計(jì)算模型

截取該工程標(biāo)準(zhǔn)段建立管廊模型，截面如圖1所示，管廊模型尺寸為13.85 m×4.6 m×20 m。管廊采用Solid單元建模，Solid65單元是專(zhuān)為混凝土、巖石等抗壓能力遠(yuǎn)大于抗拉能力的非均勻材料開(kāi)發(fā)的單元。可以較好的模擬混凝土中的配筋(或玻璃纖維、型鋼等)，以及材料的拉裂和壓潰現(xiàn)象[16-17]。管廊模型受力如圖3所示。

2.2 計(jì)算結(jié)果分析

經(jīng)數(shù)值模擬計(jì)算，管廊應(yīng)力分布如圖4所示，管廊最下層各艙整體受力均勻，各個(gè)艙室的下腋角處均存在應(yīng)力集中情況，其中污水管線(xiàn)艙和電力管線(xiàn)艙下腋角的應(yīng)力集中情況最明顯。

具體分析情況為：水力管線(xiàn)艙頂板應(yīng)力變化情況為中間大，沿左右兩邊方向逐漸減小，最大值為2.02 MPa，最小值為0.45 MPa。而管廊底板應(yīng)力變化呈相反情況，中間位置應(yīng)力最小，為0.2 KPa，最大應(yīng)力值出現(xiàn)在兩端腋角處，為1.35 MPa?；炷敛牧蠌?qiáng)度允許值為45 MPa，管廊的最大應(yīng)力值遠(yuǎn)小于規(guī)定的強(qiáng)度允許值。污水管線(xiàn)艙、天然氣管線(xiàn)艙和電力管線(xiàn)艙側(cè)壁自上往下應(yīng)力逐漸增大，最小應(yīng)力值為0.2 MPa，最大應(yīng)力值為1.12 MPa。

圖3 結(jié)構(gòu)荷載分布圖Fig.3 Distribution map of structure load

圖4 綜合管廊應(yīng)力分布云圖Fig.4 Stress distribution map of utility tunnel

圖5 綜合管廊位移云圖Fig.5 Deformation diagram of utility tunnel

由圖5分析可得：管廊整體變形趨勢(shì)呈凹字形，管廊底板不變形，頂板變形由中心向左右兩邊均勻減小，水力管線(xiàn)艙頂板中心位移最大，位移值為0.5 mm，最小位移值為零；電力管線(xiàn)艙和天然氣管線(xiàn)艙比較穩(wěn)定，整體位移接近于零；污水管線(xiàn)艙底板不變形，而側(cè)壁自下向上變形逐漸增大，變形范圍在0.2 mm左右。規(guī)范規(guī)定的撓度允許值是25 mm，管廊的最大位移值遠(yuǎn)小于規(guī)定允許值，結(jié)構(gòu)滿(mǎn)足要求。

3 結(jié)束語(yǔ)

經(jīng)過(guò)分析得出如下結(jié)論：

(1)管廊的最大位移和最大應(yīng)力均遠(yuǎn)小于設(shè)計(jì)允許值，說(shuō)明結(jié)構(gòu)設(shè)計(jì)比較保守，材料沒(méi)有充分利用，結(jié)構(gòu)還有進(jìn)一步優(yōu)化的空間。

(2)管廊各個(gè)艙室內(nèi)壁腋角處出現(xiàn)明顯的應(yīng)力集中現(xiàn)象，在結(jié)構(gòu)設(shè)計(jì)時(shí)應(yīng)適當(dāng)增加腋筋配筋率。

(3)管廊上部結(jié)構(gòu)比下部結(jié)構(gòu)變形大，設(shè)計(jì)時(shí)應(yīng)合理分配艙室大小，周?chē)馏w對(duì)管廊產(chǎn)生土壓力的同時(shí)也起到一定的約束作用。

[1] 孫厚超，王國(guó)良，楊平.考慮流固耦合的隧道開(kāi)挖引起地表沉降數(shù)值分析[J].森林工程，2015，31(3)：112-116.

[2] 朱云云，胡俊，衛(wèi)宏，等.地鐵聯(lián)絡(luò)通道凍結(jié)法加固溫度場(chǎng)數(shù)值優(yōu)化分析[J].森林工程，2017，33(1)：74-81.

[3] 薛學(xué)斌，殷吉彥，周洲，等.城市綜合管廊相關(guān)問(wèn)題探討[J].給水排水，2017，53(1)：137-142.

[4] 王恒棟.GB 50838—2015《城市綜合管廊工程技術(shù)規(guī)范》解讀[J].中國(guó)建筑防水，2016，34(14)：34-37.

[5] 劉文峰，尹力文.談南京河西新區(qū)江東南路綜合管廊設(shè)計(jì)[J].山西建筑，2014，40(36)：132-134.

[6] 曾慶紅.關(guān)于城市規(guī)劃建設(shè)指標(biāo)與地下綜合管廊、海綿城市建設(shè)銜接的探討[J].建設(shè)科技，2017，16(4)：73-74.

[7] 李軍.綜合管廊設(shè)計(jì)技術(shù)要點(diǎn)分析[J].建設(shè)科技，2017，16(5)：66-67.

[8] 周泓，姜青新.關(guān)于《關(guān)于推進(jìn)城市地下綜合管廊建設(shè)的指導(dǎo)意見(jiàn)》的初步解讀[J].環(huán)境與可持續(xù)發(fā)展，2015，40(5)：58-59.

[9] 中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.GB50009-2012，建筑結(jié)構(gòu)荷載規(guī)范[S].北京：中國(guó)建筑工業(yè)出版社，2012.

[10] 李杰，岳慶霞，陳雋.地下綜合管廊結(jié)構(gòu)振動(dòng)臺(tái)模型試驗(yàn)與有限元分析研究[J].地震工程與工程振動(dòng)，2009，29(4)：41-45.

[11] 郭健，夏鵬，殷俊.地下工程綜合管溝沉降的數(shù)值分析[J].武漢輕工大學(xué)學(xué)報(bào)，2015，34(4)：72-76.

[12] 滕振楠，田志宇，韋遠(yuǎn)飛，等.停機(jī)坪下穿隧道結(jié)構(gòu)荷載分析.西南公路[J]，2016，22(4)：73-76.

[13] 榮哲，孫玉品.城市綜合管廊設(shè)計(jì)與計(jì)算[J].工業(yè)建筑，2013，43(S1)：230-232.

[14] 喻卓爾，徐得力.城市綜合管廊工程重要節(jié)點(diǎn)設(shè)計(jì)研究[J].住宅與房產(chǎn)，2017，4(3)：255.

[15] 汪創(chuàng)，許夢(mèng)國(guó)，王平，等.車(chē)輛動(dòng)載作用下地下通道基坑穩(wěn)定性研究[J].科學(xué)技術(shù)與工程，2016，16(10)：216-220.

[16] 王魁，王浩，張達(dá).多層地下管廊使用階段受力與變形數(shù)值模擬[J].城市住宅，2017，24(1)：97-99.

[17] 姚守龍，王志鶴.關(guān)于城市地下綜合管廊的施工技術(shù)研究[J].綠色環(huán)保建材，2017，4(1)：111.

NumericalSimulationoftheStressandDeformationfortheUtilityTunnelinHuangdao

Liu Xiaoqin1，Yan Qingzhi1*，Rao Jiang2

(1.College of Pipeline and Civil Engineering，China University of Petroleum，Qingdao 266580； 2.Shengli Petroleum Management Bureau Water Supply Company，Dongying 257000)

In this paper，a period of utility tunnel in Huangdao is taken as the analysis object，the three-dimensional element model is established by using the finite element software ANSYS.Considering the soil nature in Huangdao，the construction method of utility tunnel，the pedestrian load and the traffic load，the stress distribution and the displacement change of the tunnel in the normal use stage are numerically simulated.The analysis results show that the design is conservative，the material is not fully utilized and there are room for further optimization.There are obvious stress concentration phenomena of axil angle in each cabin，which should be paid more attention to in structural design，and measures such as increasing the wail reinforcement ratio and improving the stress state should be taken.The conclusions can provide reference for similar engineering design in the future.

Utility tunnel；stress distribution；displacement；numerical simulation

TU 990.3

A

1001-005X(2017)06-0093-04

2017-05-23

國(guó)家自然科學(xué)基金(51308510)；山東省優(yōu)秀中青年科學(xué)家科研獎(jiǎng)勵(lì)基金(2014BSE28027)

劉孝泰，碩士研究生。研究方向：結(jié)構(gòu)可靠度。

*通信作者：顏慶智，博士，教授。研究方向：結(jié)構(gòu)可靠度。E-mail：yanqzhi@163.com

劉孝泰，顏慶智，饒江.黃島某地區(qū)綜合管廊應(yīng)力與變形的數(shù)值模擬[J].森林工程，2017，33(6)：93-96.