胡　俊,衛(wèi)　宏,劉　勇

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

?

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

胡俊1,2,衛(wèi)宏1,劉勇2

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

摘要：為解決現(xiàn)有人工凍結(jié)法施工后周圍地層產(chǎn)生凍脹融沉所引發(fā)的不良后果問題，設(shè)置加熱限位管對凍土帷幕的發(fā)展進行限制。運用有限元軟件分析在凍土帷幕主面上設(shè)置加熱限位管時對凍土帷幕溫度場發(fā)展的影響規(guī)律，主要得出：隨著加熱限位管鹽水溫度的升高，凍土帷幕厚度呈線性減小趨勢；限位管循環(huán)鹽水溫度越高，最終形成的凍土帷幕邊界平整性就越好，從而具有較好的限位效果；限位管應(yīng)與凍結(jié)管對齊設(shè)置在凍土帷幕主面上；限位管與凍結(jié)管距離由最終控制凍土帷幕的厚度決定；循環(huán)熱鹽水的溫度不宜過高，宜為5～10 ℃。

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

0引言

在城市地下工程中，采用人工凍結(jié)法施工后，會使周圍地層產(chǎn)生凍脹融沉現(xiàn)象，使得土的工程性質(zhì)和相鄰建筑物受到不良影響，造成地基失穩(wěn)，致使鄰近建筑物產(chǎn)生傾斜、裂縫，嚴(yán)重時會導(dǎo)致建筑物坍塌或使地下管線發(fā)生破壞等后果。為解決現(xiàn)有人工凍結(jié)法施工后周圍地層產(chǎn)生凍脹融沉所引發(fā)的不良后果問題，許多學(xué)者做了相關(guān)研究。商厚勝[1]以廣州地鐵6號線坦尾站—如意坊站區(qū)間對接段暗挖隧道工程為依托，揭示了淺覆土條件下人工凍土凍結(jié)溫度場、強制解凍溫度場的發(fā)展變化規(guī)律以及地層凍脹融沉變形規(guī)律，合理確定了跟蹤注漿控制地層融沉的技術(shù)參數(shù)；蔡海兵等[2-3]對地鐵隧道水平凍結(jié)工程地層凍脹融沉的預(yù)測方法進行了系統(tǒng)研究，建立了施工期地表凍脹融沉的歷時預(yù)測模型；胡俊等[4-5]公開了一種抑制杯型水平凍結(jié)凍脹融沉的施工方法，通過在杯型水平凍結(jié)壁上方設(shè)置2排卸壓孔兼注漿孔，可有效抑制杯型水平凍結(jié)時凍脹融沉對周圍地層及建筑物、交通和地下管線的影響，可將地層沉降量控制在8 mm，隆起量控制在5 mm；楊平等[6]公開了一種用水泥土加固法抑制水平凍結(jié)凍脹融沉的施工方法，使用水泥深層攪拌樁工藝在土體內(nèi)摻入一定量的水泥漿，從而改善土體的物理力學(xué)性質(zhì)，可抑制端頭水平凍結(jié)加固時凍脹融沉對周圍地層及建筑物的影響。

港珠澳大橋珠海連接線拱北隧道工程采用管幕凍結(jié)法施工，運用了“圓形凍結(jié)管+異形凍結(jié)管+加熱限位管”的凍結(jié)工藝，通過設(shè)置加熱限位管來達到抑制地層凍脹融沉的目的，這是國內(nèi)為數(shù)不多的實際工程中采用設(shè)置加熱限位管來抑制凍脹融沉的施工案例。文獻[7-11]為驗證拱北隧道工程限位管的效果和獲取合適的現(xiàn)場運行參數(shù)，通過采取數(shù)值模擬、室內(nèi)試驗、現(xiàn)場原位試驗等方法，對管幕凍結(jié)法中的限位管進行詳細(xì)的試驗研究，以此來證明采用限位管控制凍土帷幕的發(fā)展是可行的。為達到控制凍脹融沉的目的，可設(shè)置加熱限位管對凍土帷幕的發(fā)展進行限制，把凍土帷幕控制在一定的厚度范圍，這是一種新的抑制凍脹融沉的施工方法。本文運用有限元軟件分析在凍土帷幕主面上設(shè)置加熱限位管時對凍土帷幕溫度場發(fā)展的影響規(guī)律，以期為實際工程施工提供技術(shù)參考。

1溫度場計算相關(guān)理論

1.1熱分析的基本理論

熱分析遵循熱力學(xué)第一定律，即能量守恒定律。對于一個封閉系統(tǒng)(沒有能量的流入或流出)，有

Q-W=ΔU+ΔKE+ΔPE。

(1)

式中：Q為熱量；W為做功；ΔU為系統(tǒng)內(nèi)能；ΔKE為系統(tǒng)動能；ΔPE為系統(tǒng)勢能。

傳熱的基本形式有熱傳導(dǎo)、熱對流及熱輻射3種。對于沒有內(nèi)熱源的非穩(wěn)態(tài)三維傳熱過程，遵循以下能量控制方程[12]：

(2)

式中：ρ′為物體材料密度;c和k分別為材料的比熱容和導(dǎo)熱系數(shù)。

要求出具體的溫度場分布,需要給出具體的邊界條件，有以下3種形式[12]：

1)物體邊界上的溫度函數(shù)為

(3)

2)物體邊界上的熱流密度函數(shù)為

(4)

3)與物體接觸的流體介質(zhì)的溫度和對流換熱系數(shù)已知為

(5)

式(3)—(5)中：Γ為物體邊界;f (x，y，z，t)為已知溫度函數(shù);g(x，y，z，t)為熱流密度函數(shù);α為對流換熱系數(shù);Tf為流體介質(zhì)的溫度。

1.2溫度場控制微分方程

凍結(jié)溫度場是具有相變的傳熱問題，帶相變瞬態(tài)溫度場問題的熱量平衡控制微分方程為

(6)

在未凍區(qū)Ωu內(nèi)

(7)

式(6)—(7)中：f、u分別為凍、融狀態(tài)；Tf為正凍區(qū)Ωf內(nèi)巖土的溫度；Cf為正凍區(qū)Ωf內(nèi)巖土的體積比熱；λf為正凍區(qū)Ωf內(nèi)巖土的導(dǎo)熱系數(shù);帶有下標(biāo)u的參數(shù)為未凍區(qū)Ωf內(nèi)的相應(yīng)物理量。

2加熱限位管布設(shè)位置

加熱限位管一般可布設(shè)在凍土帷幕界面上(與凍結(jié)管錯位布置)或凍土帷幕主面上(與凍結(jié)管對齊布置)。本文將重點研究加熱限位管布設(shè)在凍土帷幕主面上的情況。

當(dāng)加熱限位管布設(shè)在凍土帷幕主面上時，影響凍土帷幕最終凍結(jié)效果的因素有很多，如限位管與凍結(jié)管的距離、循環(huán)鹽水溫度和凍結(jié)管間距等。本文考慮凍結(jié)管間距為800 mm的情況，在不同限位管與凍結(jié)管距離和不同鹽水加熱溫度條件下，研究最終凍土帷幕的限位效果。

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

圖1　模型幾何尺寸(單位：mm)

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

3.1計算基本假定

1)土層為1層，視為均質(zhì)、熱各向同性體，初始溫度取18 ℃(一般10 m以下地層恒溫帶溫度為15～20 ℃)；2)忽略水分遷移的影響。

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

建立二維溫度場數(shù)值模型，選取9節(jié)點網(wǎng)格劃分格式，網(wǎng)格劃分后的計算模型如圖2所示。依據(jù)相關(guān)報告及試驗[14-15]，模型的材料參數(shù)如表1所示。凍結(jié)前地層初始溫度取18 ℃，計算區(qū)域的外邊界看作是絕熱邊界。直接將溫度荷載施加到凍結(jié)管和加熱限位管管壁上，凍結(jié)管和限位管管壁為熱荷載邊界，以鹽水溫度作為邊界荷載，凍結(jié)期間凍結(jié)管鹽水降溫計劃見表2。根據(jù)降溫計劃，取凍結(jié)時間步為50 d，每步時間長為24 h。采用帶相變的瞬態(tài)導(dǎo)熱模型。

圖2　網(wǎng)格劃分后模型及研究路徑示意圖(單位：mm)

3.3研究路徑

為了更好地對比研究不同限位管與凍結(jié)管距離、不同限位管加熱溫度下凍土帷幕溫度場的差異，設(shè)置了1條路徑和在路徑上的11個分析點，如圖2所示。路徑1(1—11號分析點)設(shè)置在寬邊的中垂線上，兩凍結(jié)管中間為1號分析點，每隔100 mm設(shè)置1個分析點，11號分析點離凍結(jié)管1 000 mm。

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

4.1凍土帷幕邊界的平整性

加熱限位管布設(shè)在凍土帷幕主面上，平行設(shè)置在凍結(jié)管兩側(cè)，位置與凍結(jié)管對齊，與凍結(jié)管的距離分別取500、600、700、800 mm。數(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 ℃等值線。從圖3可以看出：雖然限位管與凍結(jié)管距離不同，但是隨著加熱限位管鹽水溫度的升高，凍土帷幕厚度呈線性減小趨勢；加熱限位管鹽水溫度每升高5 ℃，凍土帷幕厚度就減小約0.2 m；當(dāng)加熱限位管循環(huán)較低的鹽水溫度時(5、10、15 ℃)，凍土帷幕邊界的波浪形特征比較顯著，隨著鹽水溫度的升高(20、25 ℃)，凍土帷幕邊界變得平整，這一現(xiàn)象都存在于限位管與凍結(jié)管距離為500、600、700、800 mm的情況中，可見凍土帷幕邊界的平整性受限位管與凍結(jié)管距離的影響不大，主要是由限位管中循環(huán)鹽水的溫度而決定的，循環(huán)鹽水溫度越高，最終形成的凍土帷幕邊界平整性就越好；當(dāng)循環(huán)25 ℃的鹽水時，限位管與凍結(jié)管距離為500、600、700、800 mm的凍土帷幕最終厚度分別為0.6、0.7、0.8、1.0 m，限位管與凍結(jié)管距離每增加0.1 m，最終形成的凍土帷幕厚度就增加約0.1 m，可見限位管與凍結(jié)管距離是決定最終凍土帷幕厚度的主要因素。

表1　土體材料參數(shù)

表2　凍結(jié)管鹽水溫度降溫計劃

4.2凍土帷幕的平均溫度

單純解決限位問題而忽略凍土帷幕的力學(xué)性能問題可能會給工程帶來安全風(fēng)險。加熱限位管對凍土帷幕力學(xué)性能的削弱問題值得研究，選取2種工況進行對比：1)限位管離凍結(jié)管0.5 m，循環(huán)鹽水溫度5 ℃；2)限位管離凍結(jié)管0.8 m，循環(huán)鹽水溫度25 ℃。2種工況最終形成的凍土帷幕厚度都約為1 m，但是第1種工況凍土帷幕邊界呈波浪形，第2種工況凍土帷幕邊界比較平整(凍土帷幕邊界呈波浪形時，凍土帷幕厚度=(界面上厚度+主面上厚度)/2)。

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

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

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

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

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

圖4為2種工況下路徑1各點凍結(jié)50 d時的溫度空間分布圖。由圖4可以看出：路徑1上2種情況凍土帷幕最低溫度都接近-15 ℃，在1號分析點處；在距凍結(jié)管0.5 m處的6號分析點溫度都在0 ℃附近，故這2種情況凍土帷幕的平均溫度應(yīng)該相差不大(約-7.5 ℃)；在1號至4號分析點之間(Ⅰ區(qū))2種情況凍土帷幕的溫度幾乎一致(約-11 ℃)，4號至6號分析點(Ⅱ 區(qū))第1種工況凍土帷幕溫度略低(Ⅱ 區(qū)第1種工況平均溫度約為-5 ℃，第2種工況約為-3 ℃)。由于凍土帷幕的力學(xué)性能是與凍土帷幕的溫度密切相關(guān)的，凍土帷幕溫度越低，其力學(xué)性能也就越好，那么2種工況在Ⅰ區(qū)的凍土帷幕力學(xué)性能幾乎一致，在Ⅱ區(qū)，第1種工況凍土帷幕的力學(xué)性能要優(yōu)于第2種工況?？傮w上講，第1種工況凍土帷幕的力學(xué)性能要優(yōu)于第2種工況。

圖4　2種工況下各點不同時間溫度空間分布圖

第1種工況凍土帷幕邊界呈波浪形，第2種工況凍土帷幕邊界相對平整。波浪形邊界對整個凍土帷幕力學(xué)性能的影響還有待深入研究，但是就本文所模擬的情況來看，Ⅰ區(qū)平均溫度比Ⅱ區(qū)低得多，其強度也就比Ⅱ區(qū)高得多，在受力時真正起決定性作用的是在Ⅰ區(qū)。因此，單從凍土帷幕的平均溫度出發(fā)，第1種工況要比第2種工況限位后的凍結(jié)效果好。另外，第1種工況鹽水溫度為5 ℃，與第2種工況的25 ℃相比，節(jié)約了能源，是一種比較經(jīng)濟的方法。

4.3第1種工況詳細(xì)分析

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

(a)凍結(jié)29d (b)凍結(jié)30d (c)凍結(jié)31d (d)凍結(jié)32d

圖5開始循環(huán)5 ℃鹽水時凍土帷幕0 ℃等值線圖

Fig.50 ℃ contour maps of frozen soil curtain after using 5 ℃ saline water

圖6為第1種工況路徑1上各點溫度隨時間變化曲線。由圖6可以看出：路徑1上凍結(jié)管中間的1號分析點降溫最快，凍結(jié)10 d時溫度降到0 ℃，剩下各點離凍結(jié)管越遠(yuǎn)降溫越慢；各點在凍結(jié)30 d時均出現(xiàn)溫度明顯上升的現(xiàn)象，靠限位管越近溫度升高越多；隨著時間的推移，各點溫度趨于穩(wěn)定。

圖6　路徑1上各點溫度隨時間變化曲線

Fig.6Temperatures vs.time measured at different points in Path No.1

圖7為第1種工況各點不同時間的溫度空間分布曲線。由圖7可以看出：在加熱限位管未通熱鹽水之前，路徑1的降溫速度先快后慢，由鹽水降溫計劃所決定；凍結(jié)30 d后，由于加熱限位管的作用，路徑1溫度開始升高;凍結(jié)45 d和50 d時的曲線幾乎一致，說明凍結(jié)后期路徑1上各點溫度趨于穩(wěn)定。

圖7　路徑1上各點不同時間溫度空間分布曲線

Fig.7Temperature variations measured at different points in Path No.1

4.4限位管離凍結(jié)管0.5 m時不同熱鹽水溫度詳細(xì)分析

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

圖8　路徑1各點在凍結(jié)50 d時溫度空間分布曲線

Fig.8Temperature variations measured at different points in Path No.1 after freezing for 50 days

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

(a)1號分析點(b)5號分析點(c)8號分析點(d)11號分析點

圖9路徑1各點溫度隨時間變化曲線

Fig.9Temperatures vs.time measured at different points in Path No.1

5結(jié)論與建議

運用有限元軟件分析凍土帷幕主面上設(shè)置加熱限位管時對凍土帷幕溫度場發(fā)展的影響規(guī)律，考慮凍結(jié)管間距為800 mm的情況，在不同限位管與凍結(jié)管距離、不同鹽水加熱溫度條件下，研究最終凍土帷幕的限位效果。

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

2)隨著加熱限位管鹽水溫度的升高，凍土帷幕厚度呈線性減小趨勢；限位管循環(huán)鹽水溫度越高，最終形成的凍土帷幕邊界平整性就越好，說明具有較好的限位效果。

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

根據(jù)以上結(jié)論，為解決人工凍結(jié)法施工后地層凍脹融沉問題，可采用設(shè)置加熱限位管這一新方法；限位管設(shè)置在凍土帷幕主面上，與凍結(jié)管對齊布置；限位管與凍結(jié)管距離由最終控制凍土帷幕的厚度決定(若最終控制凍土帷幕厚度小于1 m，那么兩側(cè)限位管與凍結(jié)管的距離為500 mm)；循環(huán)熱鹽水的溫度不宜過高，一般為5～10 ℃。

參考文獻(References)：

[1]商厚勝.淺覆土人工凍結(jié)的凍脹融沉規(guī)律研究[D].北京:中國礦業(yè)大學(xué)，2015.(SHANG Housheng.Study of frost heave and thawing settlement of artificial freezing under shallow overburden[D].Beijing:China University of Mining &Technology,2015.(in Chinese))

[2]蔡海兵.地鐵隧道水平凍結(jié)工程地層凍脹融沉的預(yù)測方法及工程應(yīng)用[D].長沙：中南大學(xué)，2012.(CAI Haibing.Prediction method of thawing freezing frost formation engineering level and engineering application of subway tunnel [D].Changsha:Central South University,2012.(in Chinese))

[3]蔡海兵，程樺，姚直書,等.基于凍土正交各向異性凍脹變形的隧道凍結(jié)期地層位移數(shù)值分析[J].巖石力學(xué)與工程學(xué)報,2015,34 (8):1667-1676.(CAI Haibing,CHENG Hua,YAO Zhishu,et al.Numerical analysis of ground displacement due to orthotropic frost heave of frozen soil in freezing period of tunnel[J].Chinese Journal of Rock Mechanics and Engineering,2015,34 (8):1667-1676.(in Chinese))

[4]胡俊，李艷榮，肖天崟,等.一種抑制杯型水平凍結(jié)凍脹融沉的施工方法：中國，CN201410549588.6[P].2015-03-25.(HU Jun,LI Yanrong,XIAO Tianyin,et al.Horizontal ground freezing technique of cup-shape with controlling of frost heaving and thaw settlement：China，CN201410549588.6[P].2015-03-25.(in Chinese))

[5]胡俊,楊平.大直徑杯型凍土壁溫度場數(shù)值分析[J].巖土力學(xué),2015,36(2):523-531.(HU Jun,YANG Ping.Numerical analysis of temperature field within large-diameter cup-shaped frozen soil wall[J].Rock and Soil Mechanics，2015,36(2):523-531.(in Chinese))

[6]楊平，張翔宇，張婷，等.一種用水泥土加固法抑制水平凍結(jié)凍脹融沉的施工方法：中國，CN201210577408.6[P].2013-04-03.(YANG Ping,ZHANG Xiangyu,ZHANG Ting,et al.Horizontal ground freezing technique with cement-treated soils to control frost heaving and thaw settlement：China,CN201210577408.6[P].2013-04-03.(in Chinese))

[7]胡向東，任輝，陳錦，等.管幕凍結(jié)法積極凍結(jié)方案模型試驗研究[J].現(xiàn)代隧道技術(shù),2014,51 (5):92-98.(HU Xiangdong,REN Hui,CHEN Jin，et al.Model test study of the active freezing scheme for the combined pipe-roof and freezing method[J].Modern Tunnelling Technology，2014,51 (5):92-98.(in Chinese))

[8]李志宏，王文州，胡向東.施工熱擾動對管幕凍結(jié)止水帷幕影響研究[J].土木工程學(xué)報,2015,48 (增刊2):374-379.(LI Zhihong,WANG Wenzhou,HU Xiangdong.Research on influence of construction thermal disturbance on the freezing-sealing pipe roof[J].China Civil Engineering Journal，2015,48(S2):374-379.(in Chinese))

[9]李劍，李志宏，胡向東.管幕凍結(jié)暗挖工法凍結(jié)止水效果分析[J].地下空間與工程學(xué)報,2015,11(3):751-758.(LI Jian,LI Zhihong,HU Xiangdong.Analysis of water sealing effect of freezing-sealing pipe roof method[J].Chinese Journal of Underground Space and Engineering,2015,11(3):751-758.(in Chinese))

[10]HU Xiangdong,SHE Siyuan.Study of freezing scheme in freeze-sealing pipe roof method based on numerical simulation of temperature field[C]// Proceedings of the International Conference on Pipelines and Trenchless Technology.VA:ASCE,2012:1798-1805.

[11]張軍，胡向東，任輝.拱北隧道管幕凍結(jié)施工中限位管的凍結(jié)效果控制研究[J].隧道建設(shè)，2015,35(11):1157-1163.(ZHANG Jun,HU Xiangdong,REN Hui.Case study of control of freezing effect by installing limiting tubes in freezing-sealing pipe roof in Gongbei Tunnel[J].Tunnel Construction,2015,35(11):1157-1163.(in Chinese))

[12]Rohsenow W M，Hartnett J P，Gantic E N.傳熱學(xué)基礎(chǔ)手冊(上冊)[M].2版.齊欣，譯.北京：科學(xué)出版社，1992.(Rohsenow W M，Hartnett J P，Gantic E N.Handbook of heat transfer fundamentals[M].2nd.Translator:QI Xin.Beijing：Science Press，1992.(in Chinese))[13]張學(xué)富,蘇新民,賴遠(yuǎn)明,等.寒區(qū)隧道三維溫度場非線性分析[J].土木工程學(xué)報,2004,37(2):47 - 53.(ZHANG Xuefu,SU Xinmin,LAI Yuanming,et al.Non-linear analysis for three-dimensional temperature fields in cold-region tunnels[J].China Civil Engineering Journal,2004,37(2):47-53.(in Chinese))

[14]胡俊.水泥改良前后土體凍結(jié)溫度及力學(xué)特性試驗研究[J].鐵道建筑,2013,18(4):156-159.(HU Jun.Experimental study on freezing temperature and mechanical performance of soil before and after cement-improving[J].Railway Engineering,2013,18(4):156-159.(in Chinese))

[15]胡俊.高水壓砂性土層地鐵大直徑盾構(gòu)始發(fā)端頭加固方式研究[D].南京：南京林業(yè)大學(xué)，2012.(HU Jun.Study of end soil reinforcement methods of large-diameter shield launching in sandy clay with high water pressure[D].Nanjing:Nanjing Forestry University,2012.(in Chinese))

Numerical Analysis of Temperature Field of Frozen Soil Curtain with Thermal Limiting Pipes

HU Jun1,2,WEI Hong1,LIU Yong2

(1.College of Civil Engineering and Architecture,Hainan University,Haikou 570228,Hainan,China;2.Department of Civil &Environmental Engineering,National University of Singapore,Kent Ridge 117576,Singapore)

Abstract:The thermal limiting pipes are used to lay restrains on the development of frozen soil curtain.The influencing rules of thermal limiting pipe on development of frozen soil curtain are analyzed by means of finite-element method.The conclusions are drawn as follows:1) The frozen soil curtain thickness decreases linearly as the water temperature of thermal limiting pipe increases.2) A more even boundary of the frozen soil curtain can be obtained under a higher temperature of the water of thermal limiting pipes.3) The thermal limiting pipes should be aligned perpendicularly to the freezing pipes.4) The distance between thermal limiting pipe and freezing pipe is controlled by the thickness of frozen soil curtain.5) The temperature of thermal limiting pipes is recommended to be 5-10 ℃.

Keywords:thermal limiting pipe;frozen soil curtain;temperature field;ground freezing method;heaving thaw settlement;numerical simulation

收稿日期：2015-08-27;修回日期:2016-03-03

基金項目：中國博士后科學(xué)基金資助項目(2015M580559)；海南省科技項目(ZDXM2015117)；海南省重點研發(fā)計劃科技合作方向項目(ZDYF2016226)

第一作者簡介：胡俊(1983—)，男，四川樂山人，2012年畢業(yè)于南京林業(yè)大學(xué)，森林工程專業(yè)(環(huán)境巖土與地下工程方向)，博士后，副教授，主要從事隧道與地下工程方面的教學(xué)與研究工作。E-mail：hj7140477@hainu.edu.cn。

DOI:10.3973/j.issn.1672-741X.2016.06.006

中圖分類號：U 455

文獻標(biāo)志碼：A

文章編號：1672-741X(2016)06-0688-07