胡　俊，衛(wèi)　宏，劉　勇，趙聯(lián)楨

(1.海南大學(xué) 土木建筑工程學(xué)院，?？?570228；2.新加坡國立大學(xué) 土木與環(huán)境工程系，新加坡 肯特崗 117576)

?

凍土帷幕界面設(shè)置加熱限位管時溫度場數(shù)值分析

胡俊1，2，衛(wèi)宏1，劉勇2*，趙聯(lián)楨1

(1.海南大學(xué) 土木建筑工程學(xué)院，?？?570228；2.新加坡國立大學(xué) 土木與環(huán)境工程系，新加坡 肯特崗 117576)

摘要：為解決現(xiàn)有人工凍結(jié)法施工后周圍地層產(chǎn)生凍脹融沉所引發(fā)不良后果的問題，可設(shè)置加熱限位管對凍土帷幕的發(fā)展進行限制，從而達到控制凍脹融沉的目的。本文運用有限元軟件研究在凍土帷幕界面上設(shè)置加熱限位管時對凍土帷幕溫度場發(fā)展的影響規(guī)律，主要得出：加熱限位管鹽水溫度每升高5℃，凍土帷幕厚度就減小約0.2 m；在已凍結(jié)30 d，凍土厚度發(fā)展到1.6 m凍土上，加熱限位管循環(huán)5℃鹽水時，隨著時間的增加凍土厚度慢慢變小，加熱限位管作用效果明顯；在限位管開始循環(huán)熱水前期，各點溫度都有明顯上升，離限位管越近溫度所受影響越大；隨著時間的推移，各點溫度趨于穩(wěn)定，循環(huán)熱鹽水溫度越高，趨于穩(wěn)定的溫度值也越高。所得結(jié)果可為今后類似工程設(shè)計提供參考依據(jù)。

關(guān)鍵詞：加熱限位管；凍結(jié)法；凍脹融沉；數(shù)值模擬

0引言

人工凍結(jié)技術(shù)在城市地下工程建設(shè)過程中和施工后，會使周圍地層產(chǎn)生凍脹融沉現(xiàn)象：一是使鄰近建筑物產(chǎn)生傾斜、裂縫，嚴重時會導(dǎo)致建筑物坍塌等事故；二是使地下管線發(fā)生破壞等不良后果[1-6]。為了解決現(xiàn)有人工凍結(jié)法施工后周圍地層產(chǎn)生凍脹融沉所引發(fā)不良后果的問題，可設(shè)置加熱限位管對凍土帷幕的發(fā)展進行限制，把凍土帷幕控制在一定厚度范圍以內(nèi)，從而達到控制凍脹融沉的目的[7-12]。實際工程中，通過設(shè)置加熱限位管來主動地抑制凍脹融沉的方法在施工中還比較少見，關(guān)于加熱限位管對凍土帷幕溫度場發(fā)展的影響規(guī)律還缺乏深入的研究。本文運用有限元軟件分析在凍土帷幕界面上設(shè)置加熱限位管時對凍土帷幕溫度場發(fā)展的影響規(guī)律，以期對實際工程施工提供技術(shù)參考依據(jù)。

1溫度場數(shù)值模型的建立

1.1數(shù)值模擬情況簡介

本文采用二維數(shù)值模型進行分析計算，土層尺寸取長4000 mm、寬2400 mm的矩形平面，凍結(jié)管設(shè)置在長邊的中垂線上，間距800 mm；加熱限位管平行設(shè)置在凍結(jié)管兩側(cè)，位置設(shè)于兩個凍結(jié)管中間，距長邊中垂線500 mm(希望將凍土帷幕厚度控制在1 m以內(nèi))，間距也是800 mm；凍結(jié)管和加熱限位管直徑同為127 mm。在實施積極凍結(jié)29 d后，從凍結(jié)30 d開始在加熱限位管中循環(huán)5℃、10、15、20、25℃的熱鹽水，比較不同加熱溫度下凍土帷幕溫度場的差異，整個凍結(jié)時間取50 d。模型幾何尺寸及凍結(jié)管和限位管布置形式如圖1所示。

圖1　模型幾何尺寸(mm)Fig.1 Model size(mm)

1.2計算基本假定

假定土層具有均勻的初始溫度場，初始溫度取18℃(上海地區(qū)一般地層10 m以下恒溫帶溫度為15～20℃)；土層為一層，視為均質(zhì)、熱各向同性體；直接將溫度荷載施加到凍結(jié)管和加熱限位管管壁上；忽略水分遷移的影響。

1.3計算模型和參數(shù)選取

本文建立二維溫度場數(shù)值模型，選取了九節(jié)點網(wǎng)格劃分格式，網(wǎng)格劃分后的計算模型如圖2所示。

圖2　網(wǎng)格劃分后模型及研究路徑示意圖(mm)Fig.2 Model mesh and analysis paths(mm)表1　土體材料參數(shù) Tab.1 Parameters of soil materials

密 度/(kg·m-3)導(dǎo)熱系數(shù)/(kJ·m-1·d-1·℃-1)未凍土凍土比熱/(kJ·kg-1·℃-1)未凍土凍土相變潛熱/(×108J/m3)18801181791.531.611.20

表2　鹽水溫度降溫計劃Tab.2 Freezing time plan for brine temperature

依據(jù)相關(guān)報告及試驗[13-19]，模型的材料參數(shù)見表1。凍結(jié)前地層初始溫度取18℃，凍結(jié)管和限位管管壁為熱荷載邊界，以鹽水溫度作為邊界荷載，凍結(jié)期間凍結(jié)管鹽水降溫計劃見表2。根據(jù)降溫計劃，取凍結(jié)時間步為50 d，每步時間長為24 h。采用帶相變的瞬態(tài)導(dǎo)熱模型。

1.4研究路徑

為了更好地對比研究不同限位管加熱溫度下凍土帷幕溫度場的差異，分別設(shè)置了2條路徑和在路徑上的12個分析點，如圖2所示。路徑1(1～6號分析點)設(shè)置在寬邊的中垂線上，兩凍結(jié)管中間為1號點，6號點離上方限位管100 mm；路徑2(7～12號分析點)也垂直于寬邊設(shè)置，兩限位管中間為7號點，12號點離下方凍結(jié)管100 mm。

2溫度場計算結(jié)果與分析

2.1不同工況下的凍土帷幕情況

數(shù)值模擬了幾種不同的工況，一種是在整個凍結(jié)50 d內(nèi)，加熱限位管不起作用，無需循環(huán)熱鹽水；剩下幾種工況分別為：在實施積極凍結(jié)29 d后，從凍結(jié)30 d開始加熱限位管中循環(huán)5、10、15、20、25℃的熱鹽水。圖3為凍結(jié)50 d時不同工況下凍土帷幕0℃等值線，可以看出：當(dāng)加熱限位管不起作用時，間距800 mm的單排凍結(jié)管在凍結(jié)50 d時凍土帷幕厚度可以發(fā)展到將近2.4 m；當(dāng)從凍結(jié)30 d開始加熱限位管循環(huán)5、10、15、20、25℃的熱鹽水之后，凍結(jié)50 d時的凍土帷幕厚度約為1.4、1.2、1.0、0.8、0.6 m；加熱限位管鹽水溫度每升高5℃，凍土帷幕厚度就減小約0.2 m；隨著加熱限位管鹽水溫度的升高，凍土帷幕厚度呈線性減小。

(a)無需加熱

(b)循環(huán)5℃鹽水

(c)循環(huán)10℃鹽水

(d)循環(huán)15℃鹽水

(e)循環(huán)20℃鹽水

(f)循環(huán)25℃鹽水圖3　凍結(jié)50 d時不同工況下凍土帷幕0℃等值線圖Fig.3 Contour of 0℃ for 50 freezing days under different working conditions

2.2加熱限位管循環(huán)5℃鹽水時

圖4為加熱限位管開始循環(huán)5℃鹽水時凍土帷幕0℃等值線(凍結(jié)30～33 d)。

(a)凍結(jié)30 d

(b)凍結(jié)31 d

(c)凍結(jié)32 d

(d)凍結(jié)33 d圖4　開始循環(huán)5℃鹽水時凍土帷幕0℃等值線圖Fig.4 Contour of 0℃ with 5℃ brine at different freezing time

由圖4可以看出：在加熱限位管循環(huán)熱鹽水之前(凍結(jié)29 d)，凍土帷幕厚度發(fā)展到1.6 m，加熱限位管被凍土帷幕包圍，在凍土帷幕之中；凍結(jié)30 d開始，加熱限位管中循環(huán)5℃熱鹽水，限位管周圈的凍土帷幕開始解凍，解凍范圍以限位管中心為圓心呈同心圓向外發(fā)展，到凍結(jié)33 d時0℃等值線相交變成連續(xù)的封閉等值線，之后隨著凍結(jié)時間的增加，凍土帷幕厚度慢慢變小，加熱限位管起到了限制凍土帷幕發(fā)展的作用。

圖5為循環(huán)5℃鹽水時路徑上各點溫度隨時間變化曲線圖?？梢钥闯觯郝窂?上兩凍結(jié)管中間的1號分析點降溫最快，凍結(jié)10 d時溫度降到0℃，路徑2靠近凍結(jié)管的12號分析點降溫最快；剩下各點離凍結(jié)管越遠降溫越慢；各點在凍結(jié)30 d時均出現(xiàn)溫度明顯上升的現(xiàn)象，距離限位管最近的6號分析點溫度變化最大，溫度升高將近5℃，靠限位管越近溫度升高越多；隨著時間的推移，各點溫度趨于穩(wěn)定。

(a)路徑1

(b)路徑2圖5　循環(huán)5℃鹽水時路徑上各點溫度隨時間變化圖Fig.5 Temperature change at different paths with 5℃ brine

圖6為循環(huán)5℃鹽水時各點不同時間的溫度空間分布曲線?？梢钥闯觯涸诩訜嵯尬还芪赐猁}水之前，2條路徑的降溫速度先快后慢，由鹽水降溫計劃所決定；凍結(jié)30 d后，由于加熱限位管的作用，2條路徑溫度開始升高，凍結(jié)45 d和50 d時的曲線幾乎一致，說明凍結(jié)后期2條路徑上各點溫度趨于穩(wěn)定。

2.3分析點對比分析

圖7為凍結(jié)50 d時不同工況下各點溫度空間分布圖，可以看出：無論路徑1還是路徑2，限位管無需加熱時和循環(huán)5℃鹽水時的溫差很大，以1號分析點為例，限位管無需加熱時與循環(huán)5℃鹽水時的溫差約為6.9℃，而循環(huán)5℃鹽水時和循環(huán)25℃鹽水時的溫差約為5.7℃，說明加熱限位管對于凍土帷幕發(fā)展的限制作用效果十分顯著，離限位管越近溫度所受影響越大。

(a)路徑1

(b)路徑2圖6　循環(huán)5℃鹽水時各點不同時間溫度空間分布圖Fig.6 Temperature change at different points with 5℃ brine

(a)路徑1

(b)路徑2圖7　不同工況下各點在凍結(jié) 50天時溫度空間分布圖Fig.7 Spatial distribution of the temperature of each point at different working conditions with 50 freezing days

圖8為不同工況下各點溫度隨時間變化圖，可以看出：在限位管開始循環(huán)熱水前期，各點溫度都有明顯的上升，離限位管越近溫度所受影響越大，無需加熱和循環(huán)25℃鹽水時的溫差也越大；隨著時間的推移，各點溫度趨于穩(wěn)定，循環(huán)熱水溫度越高，趨于穩(wěn)定的溫度值也越高。

(a)1號分析點

(b)6號分析點

(c)7號分析點

(d)12號分析點圖8　不同工況下各點溫度隨時間變化圖Fig.8 Temperature change with freezing time at different points under different construction situation

3結(jié)論

本文運用有限元軟件，分析了在凍土帷幕界面上設(shè)置加熱限位管時對凍土帷幕溫度場發(fā)展的影響規(guī)律，通過在加熱限位管中循環(huán)不同的加熱鹽水溫度，數(shù)值對比分析了各加熱溫度下凍土帷幕溫度場的差異，主要得出：

(1)為解決現(xiàn)有人工凍結(jié)法施工后周圍地層產(chǎn)生凍脹融沉所引發(fā)不良后果的問題，可設(shè)置加熱限位管對凍土帷幕的發(fā)展進行限制，把凍土帷幕控制在一定厚度范圍內(nèi)，從而達到控制凍脹融沉的目的。

(2)加熱限位管鹽水溫度每升高5℃，凍土帷幕厚度就減小約0.2 m；隨著加熱限位管鹽水溫度的升高，凍土帷幕厚度呈線性減小。

(3)加熱限位管循環(huán)5℃鹽水時，在凍結(jié)30 d以前凍土帷幕厚度發(fā)展到1.6 m，加熱限位管被凍土帷幕包圍，之后隨著凍結(jié)時間的增加，凍土帷幕厚度慢慢變小，加熱限位管起到了限制凍土帷幕發(fā)展的作用。

(4)加熱限位管循環(huán)5℃鹽水時，各點在凍結(jié)30天均出現(xiàn)溫度明顯上升的現(xiàn)象，距離限位管最近的6號分析點溫度升高將近5℃，靠限位管越近溫度升高越多；隨著時間的推移，各點溫度趨于穩(wěn)定。

(5)在限位管開始循環(huán)熱水前期，各點溫度都有明顯的上升，離限位管越近溫度所受影響越大；隨著時間的推移，各點溫度趨于穩(wěn)定，循環(huán)熱鹽水溫度越高，趨于穩(wěn)定的溫度值也越高。

【參考文獻】

[1]胡俊，王效賓，袁云輝.盾構(gòu)隧道端頭杯型凍結(jié)壁溫度場發(fā)展與分布規(guī)律研究[M].北京：中國水利水電出版社，2015.

[2]胡俊，楊平.大直徑杯型凍土壁溫度場數(shù)值分析[J].巖土力學(xué)，2015，36(2)：523-531.

[3]Zhou J，Tang Y Q.Artificial ground freezing of fully saturated mucky clay：Thawing problem by centrifuge modeling[J].Cold Regions Science and Technology，2015，117：1-11.

[4]Zeng H，Hu J，Yang P.A numerical simulation study on the chemical reinforcement area at shield start shaft[A].2011 International Conference on Electric Technology and Civil Engineering(ICETCE)[C].IEEE，2011，4：29-34.

[5]Hu J，Yang P，Dong C，et al.Study on numerical simulation of cup-shaped horizontal freezing reinforcement project near shield launching[A].2011 International Conference on Electric Technology and Civil Engineering(ICETCE)[C].IEEE，2011，4：5522-5525.

[6]Tang Y Q，Zhou J，Hong J，et al.Quantitative analysis of the microstructure of Shanghai muddy clay before and after freezing[J].Bull.Eng.Geol.Environ.，2012，71：309-316.

[7]Hu X，She S.Study of freezing scheme in freeze-sealing pipe roof method based on numerical simulation of temperature field[A].Proceedings of the international conference on pipelines and trenchless technology[C].US：American Society of Civil Engineers，2012：1798-1805.

[8]李志宏，王文州，胡向東.施工熱擾動對管幕凍結(jié)止水帷幕影響研究[J].土木工程學(xué)報，2014，48(S2)：374-379.

[9]李劍，李志宏，胡向東.管幕凍結(jié)暗挖工法凍結(jié)止水效果分析[J].地下空間與工程學(xué)報，2015，11(3)：751-758.

[10]張軍，胡向東，任輝.拱北隧道管幕凍結(jié)施工中限位管的凍結(jié)效果控制研究[J].隧道建設(shè)，2015，35(11)：1157-1163.

[11]路耀邦，路耀平，劉洪震，等.盾構(gòu)法隧道球狀風(fēng)化體處理方法研究綜述[J].公路工程，2014，39(2)：139-142.

[12]胡向東，任輝，陳錦，等.管幕凍結(jié)法積極凍結(jié)方案模型試驗研究[J].現(xiàn)代隧道技術(shù)，2014，51(5)：92-98.

[13]胡俊，劉勇，李玉萍.凍結(jié)水泥土攪拌樁溫度場數(shù)值分析[J].森林工程，2015，31(5)：118-123.

[14]胡俊.水泥改良前后土體凍結(jié)溫度及力學(xué)特性試驗研究[J].鐵道建筑，2013，18(4)：156-159.

[15]Hu J，Zeng H，Wang X.Experimental research on the physi-mechanical performances of geosynthetics [J].Applied Mechanics and Materials，2013，341-342(5)：33-37.

[16]胡俊.高水壓砂性土層地鐵大直徑盾構(gòu)始發(fā)端頭加固方式研究[D].南京：南京林業(yè)大學(xué)，2012，6.

[17]胡俊，劉勇，曾暉.新型管幕凍結(jié)法不同管幕填充形式的溫度場數(shù)值對比分析[J].森林工程，2015，31(6)：135-141.

[18]林清輝，嚴佳佳.管幕超前支護淺埋暗挖隧道開挖方案對比研究[J].公路工程，2015，40(4)：106-110.

[19]胡俊，劉勇，梁乾乾.設(shè)置1-2根圓形凍結(jié)管時凍結(jié)水泥土攪拌樁溫度場數(shù)值對比分析[J].森林工程，2016，32(1)：77-82.

收稿日期：2015-10-28

基金項目：國家自然科學(xué)基金項目(51368017);海南省科技項目(ZDXM2015117);海南省重點研發(fā)計劃科技合作方向項目(ZDYF2016226);中國博士后科學(xué)基金資助項目(2015M580559)；海南省教育廳高等學(xué)校科研項目(Hnky2015-10)

第一作者簡介：胡俊，博士，副教授。研究方向：隧道及地下工程。 *通信作者：劉勇，博士，研究員。研究方向：隧道及地下工程方面的數(shù)值模擬。E-mail：ceeliuy@gmail.com

中圖分類號：TU 753

文獻標識碼：A

文章編號：1001-005X(2016)04-0070-05

Numerical Simulation of Temperature Field based on Frozen Soil Curtain Tube with Thermal Heat Limit

Hu Jun1，2，Wei Hong1，Liu Yong2*，Zhao Lianzhen1

(1.College of Civil Engineering and Architecture，Hainan University，Haikou 570228；2.Department of Civil & Environmental Engineering，National University of Singapore，Kentridge 117576，Singapore)

Abstract：In order to solve the problem of the existing artificial freezing method after construction that the surrounding formation may produce frost thaw settlement，which can cause adverse consequences，the heating limit can be used for the development of frozen soil wall pipe to limit thickness，so as to control frost thaw settlement.In this study，the finite element method was used to study the influence of temperature field of frozen soil curtain development when setting the heating limit pipe，the main conclusions were：the temperature increase of brine in the heating limit pipe by 5℃ will cause the reduction of frozen soil wall thickness by about 0.2 m；When the frost was lasting for 30 days and the permafrost thickness reached 1.6 m，the permafrost thickness gradually decreased as time increased with 5℃ brine，and the heating effect of limiting tube was obvious.Before circulating the brine，the temperature at each point increased significantly and the closer to the tube the greater the impact.As the time went on，the temperature at each point tended to stabilize，and the higher the circulating hot brine temperature the higher the stable temperature.The results can provide a theoretical reference for future similar projects.

Keywords：thermal limiting pipe；freezing method；heaving thaw settlement；numerical simulation

引文格式：胡俊，衛(wèi)宏，劉勇，等.凍土帷幕界面設(shè)置加熱限位管時溫度場數(shù)值分析[J].森林工程，2016，32(4)：70-74.