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

劉孝泰，顏慶智*，饒江

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

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

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

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

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

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

0 引言

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

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

1 工況分析

1.1 材料參數(shù)

計算所涉及的材料主要是巖土、鋼筋、混凝土，鋼筋材料型號為HRB400、HPB300，混凝土強度等級為C45，混凝土材料參數(shù)取值見表2，鋼筋材料參數(shù)取值見表3，該段開挖深度為8 m，土斷面示意圖如圖2所示。

土體參數(shù)根據(jù)地勘報告取值見表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抗拉強度設(shè)計值ft/kPa1.80×103抗壓強度設(shè)計值fc/kPa2.11×104

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

1.2 荷載組合

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

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

(1)

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

1.3 荷載的計算

(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計算，行車荷載一般按20 kPa計算；K為靜止土壓力系數(shù)；K=1-sinφ(由表1取值)，采用水土分算模型，水的重度取10 kN/m3，在計算側(cè)壁活荷載時，汽車荷載不考慮動力系數(shù)。

(4)水壓力

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

pω=γω·hω。

(9)

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

2 數(shù)值模擬

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

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

2.2 計算結(jié)果分析

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

具體分析情況為：水力管線艙頂板應(yīng)力變化情況為中間大，沿左右兩邊方向逐漸減小，最大值為2.02 MPa，最小值為0.45 MPa。而管廊底板應(yīng)力變化呈相反情況，中間位置應(yīng)力最小，為0.2 KPa，最大應(yīng)力值出現(xiàn)在兩端腋角處，為1.35 MPa?；炷敛牧蠌姸仍试S值為45 MPa，管廊的最大應(yīng)力值遠小于規(guī)定的強度允許值。污水管線艙、天然氣管線艙和電力管線艙側(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分析可得：管廊整體變形趨勢呈凹字形，管廊底板不變形，頂板變形由中心向左右兩邊均勻減小，水力管線艙頂板中心位移最大，位移值為0.5 mm，最小位移值為零；電力管線艙和天然氣管線艙比較穩(wěn)定，整體位移接近于零；污水管線艙底板不變形，而側(cè)壁自下向上變形逐漸增大，變形范圍在0.2 mm左右。規(guī)范規(guī)定的撓度允許值是25 mm，管廊的最大位移值遠小于規(guī)定允許值，結(jié)構(gòu)滿足要求。

3 結(jié)束語

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[16] 王魁，王浩，張達.多層地下管廊使用階段受力與變形數(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

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

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

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

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