擠塑聚苯板在混凝土襯砌渠道凍脹防治中的應(yīng)用

鹿翔宇,劉仰玉,卞漢兵,3*,孫兆輝,王宸宇,邱秀梅*

擠塑聚苯板在混凝土襯砌渠道凍脹防治中的應(yīng)用

鹿翔宇1,劉仰玉2,卞漢兵1,3*,孫兆輝1,王宸宇1,邱秀梅1*

1. 山東農(nóng)業(yè)大學(xué) 水利土木工程學(xué)院, 山東 泰安 271000 2. 臨沂市恒泰安全科技有限公司, 山東 臨沂 276000 3. 法國里爾大學(xué) 綜合理工工程師學(xué)院, 法國

針對季節(jié)性凍土區(qū)山東省聊城市位山灌區(qū)混凝土襯砌渠道的凍融破壞問題，利用ANSYS有限元軟件建立鋪設(shè)不同厚度擠塑聚苯板的混凝土襯砌渠道模型，計算其瞬態(tài)溫度場分布和位移場分布，并結(jié)合相關(guān)實測數(shù)據(jù)進行分析。結(jié)果表明，經(jīng)一年持續(xù)觀測，在日最低氣溫隨時間變化的外界環(huán)境條件下，不論是否鋪設(shè)擠塑聚苯板，最大凍結(jié)深度出現(xiàn)的時間相對于年最低氣溫出現(xiàn)的時間來說均具有滯后性。通過對鋪設(shè)不同厚度擠塑聚苯板渠道產(chǎn)生的凍脹量和渠基土的凍結(jié)深度進行比對分析，得出擠塑聚苯板的保溫作用原理，推測選用5 cm擠塑聚苯板對該灌渠進行凍脹防治已滿足工程要求。本研究加入了瞬態(tài)溫度場計算，這一方案具有一定的先進性，可為今后混凝土襯砌渠道凍脹防治研究提供參考。

混凝土襯砌渠道; 凍脹防治; 數(shù)值模擬; 擠塑聚苯板（XPS）

我國是一個農(nóng)業(yè)大國，我國的灌區(qū)，特別是北方灌區(qū)大多處于季節(jié)性凍土區(qū)，灌溉渠道多發(fā)生凍融破壞[1-3]。對于混凝土渠道襯砌，渠基土的凍融作用致使渠道混凝土襯砌體發(fā)生破壞，輕者混凝土襯砌板表面出現(xiàn)裂縫或局部鼓起，重者襯砌板坍塌滑坡造成渠道不能正常輸水，每年都要花大量的人力、物力、財力進行維修維護。渠道混凝土襯砌結(jié)構(gòu)的凍脹過程是一個隨溫度變化，由表及里、由淺入深的逐步演變過程[4]。隨著外界溫度的逐漸降低，渠道下臥土層所產(chǎn)生的冰晶體越積越大[5-7]，土體凍脹，從而造成襯砌板鼓起、產(chǎn)生裂縫，最終導(dǎo)致襯砌板脫落發(fā)生凍融破壞[8,9]。針對我國灌渠凍脹破壞問題，國內(nèi)學(xué)者對西北、東北部灌渠進行了大量的分析研究[10-14]。

目前利用有限元軟件對渠道凍脹機理進行分析是一種常見方法，王正中等人[15]提出可以對渠道進行瞬態(tài)溫度場計算，并模擬了日內(nèi)氣溫的變化，對于抗凍脹設(shè)計來說，瞬態(tài)溫度場計算相對常用的穩(wěn)態(tài)溫度場計算更具有參考價值。

本文利用有限元軟件ANSYS，以山東省聊城市黃河流域位山灌區(qū)二干渠觀測段作為研究對象，參考相關(guān)資料[16-19]，建立仿真計算模型。在保溫性能方面，擠塑聚苯板（以下簡稱XPS）材料緊密的閉孔結(jié)構(gòu)決定了保溫的同時還可有效防滲，且具有導(dǎo)熱系數(shù)低、防腐蝕性好等優(yōu)勢，因此選用XPS作為該灌渠的保溫材料進行模擬研究。在有限元計算方面，由于瞬態(tài)溫度場計算相對于穩(wěn)態(tài)溫度場計算來說更具有先進性，且國內(nèi)針對完整凍結(jié)期內(nèi)灌渠凍深的模擬計算較少，因此對一年內(nèi)日最低氣溫條件下該灌渠的瞬態(tài)溫度場、位移場進行計算，分析和預(yù)測其凍深及凍脹量的變化規(guī)律。并將不同厚度XPS與混凝土襯砌板復(fù)合的渠道模型仿真計算結(jié)果與該混凝土襯砌渠道的仿真計算結(jié)果相比對，分析XPS對該灌渠的保溫作用，選取經(jīng)濟合理的防治措施。

1 工程背景

本文所做仿真計算原型為山東省聊城市位山灌區(qū)二干渠觀測段，該原型渠道如圖1，渠道邊坡比為1：2，坡長8.80 m，渠底寬1.80 m，左右水平邊界均為0.60 m。渠底不鋪設(shè)混凝土襯砌板，而在與渠坡轉(zhuǎn)折（坡腳）處設(shè)置0.50×0.50 mM7.5漿砌塊石，渠底回填粘性土壓實。其中直線1、2、3為試驗中選取的三條不同的觀測線。

圖 1 位山灌區(qū)二干渠渠道詳圖

2 試驗組與對照組選取

為將鋪設(shè)不同厚度XPS的渠道模型凍深及凍脹量與原渠道相比對，將不鋪設(shè)XPS的渠道模型作為對照組，鋪設(shè)0.01、0.03、0.05 mXPS的渠道模型分別作為試驗組1、試驗組2、試驗組3。

3 參數(shù)選取

由于瞬態(tài)溫度場計算時僅需要材料的密度、導(dǎo)熱系數(shù)和比熱容三項熱力學(xué)參數(shù)，因此參考該灌區(qū)相關(guān)試驗數(shù)據(jù)及相關(guān)資料[20,21]，整理渠道各個組成材料的熱力學(xué)相關(guān)參數(shù)見表1。

表 1 渠道結(jié)構(gòu)組成材料熱力學(xué)相關(guān)參數(shù)

由于凍土導(dǎo)熱系數(shù)在各個方向的差異不大，所以取坐標(biāo)軸方向上的導(dǎo)熱系數(shù)λ=λ=1.55。

位移場計算時需增加材料的彈性模量、泊松比和熱膨脹系數(shù)，因此根據(jù)瞬態(tài)溫度場計算結(jié)果，需要對渠基土的彈性模量和泊松比進行分層定義。

假定渠基土為各向同性的彈性材料，聊城市土質(zhì)為粉質(zhì)粘土，由于淤泥質(zhì)粉質(zhì)粘土的彈性模量與溫度變化是線性的，泊松比與溫度的關(guān)系呈近似線性關(guān)系，隨溫度的降低而減小，則其線性回歸關(guān)系分別為[22]：=19.82||-31.367 (1)

=-0.003||+0.2933 (2)

可推測聊城市位山灌區(qū)二干渠渠基土-3 ℃到-18 ℃的彈性模量和0 ℃到-18 ℃的泊松比。相關(guān)資料顯示粉質(zhì)粘土的彈性模量為7~20 MPa，0 ℃到-3 ℃的土層彈性模量可近似取20 MPa。則彈性模量及泊松比隨溫度的變化情況見圖2。

圖 2 渠基土彈性模量及泊松比隨溫度的變化曲線

取部分土樣進行室內(nèi)凍脹試驗，即取用位山灌區(qū)二干渠渠基土制備直徑100 mm，高65 mm的圓柱體試樣，先將試樣壓實后進行飽和[23]，放置于凍脹量監(jiān)測設(shè)備中?？刂圃嚇訙囟冉抵? ℃觀測其溫度和變形，隨后以-3 ℃/h的速度讓試樣降至-18 ℃，并保持-18 ℃持續(xù)觀測72 h后得到相關(guān)凍脹試驗數(shù)據(jù)見表2。

表 2 渠基土試塊凍脹試驗數(shù)據(jù)

將凍土的統(tǒng)一凍脹系數(shù)/（其中為凍脹率，為溫度）作為負(fù)膨脹系數(shù)進行熱力學(xué)計算。由于渠基土為粉質(zhì)黏土，凍結(jié)溫度接近于0 ℃，故以0 ℃為基準(zhǔn)。

混凝土襯砌板視為各向同性的彈塑性材料，其彈性模量為26000 MPa，泊松比為0.167，線膨脹系數(shù)為1.1×10-5。

擠塑聚苯板（XPS）的彈性模量可達6-8 Mpa[24-26]，泊松比近乎為0，線膨脹系數(shù)參考傳統(tǒng)聚苯板，取0.04。

4 瞬態(tài)溫度場計算

4.1 邊界條件

由相關(guān)觀測數(shù)據(jù)得到聊城市1951~2008年日最低氣溫-時間變化曲線見圖3，其中取部分?jǐn)?shù)據(jù)計算得年最低平均氣溫約為12攝氏度。為得到最大凍脹量，假設(shè)日最低氣溫-時間變化曲線為規(guī)律的正弦函數(shù)，并將其作為溫度荷載施加在渠道模型上邊界處。

圖 3 聊城市1951~2008年日最低氣溫變化曲線

上表面邊界溫度根據(jù)該環(huán)境溫度變化曲線可以將其簡化為正弦曲線形式，如下：

()=sin(+) (3)

其中：為溫度變化幅值，根據(jù)相關(guān)資料，可取18；為時間，單位s；

下邊界深度足夠達到穩(wěn)定溫度，取值為12 ℃，左右邊界設(shè)為絕熱條件。

4.2 計算模型

由于渠道是細(xì)長結(jié)構(gòu)，分析過程中可將其簡化為二維平面問題進行處理。幾何模型的建立需要考慮邊界條件等因素，當(dāng)土層達到足夠深度時，渠基土溫度可假設(shè)為恒溫，不產(chǎn)生位移變形及溫度變化。因此先對模型簡化進行一維計算。當(dāng)外界環(huán)境溫度為年最低氣溫時，實測聊城市凍土深度約為0.45 m左右。模型達足夠深時，下邊界可視為常年溫度不變層，因此假設(shè)該一維模型深度分別為3 m、10 m，其溫度隨深度的分布圖見圖4。

圖 4 環(huán)境溫度為年最低氣溫下的不同一維模型溫度隨深度的變化情況

由圖4可知，環(huán)境溫度為年最低氣溫時，渠基深度分別為3 m、10 m的有限元模型溫度隨深度的變化曲線近乎重合。根據(jù)計算數(shù)據(jù)得，當(dāng)深度為10 m，環(huán)境溫度為年最低氣溫時，其凍深（土溫為0 ℃時對應(yīng)的渠基土深度）為0.4451 m；當(dāng)深度為3 m，環(huán)境溫度為年最低氣溫時，其凍深為0.4447 m。因深度分別為10 m與3 m時的凍結(jié)深度相差僅0.0004 m，可知深度為3 m時，模型計算結(jié)果已穩(wěn)定。故渠基土下邊界向下取3.00 m，根據(jù)實際工況，左右邊界取0.60 m，混凝土襯砌板厚度為0.06 m。因具體計算研究的破壞面為圖1所示三條直線對應(yīng)的破壞面，故將斜坡面、坡頂?shù)幕炷烈r砌體分別視為不同的整體，得以更詳盡觀察其連接處的凍深變化。擠塑聚苯板鋪設(shè)在混凝土襯砌板下方，并假設(shè)渠基土與襯砌體之間無接觸面存在。各材料參數(shù)根據(jù)上文中數(shù)據(jù)導(dǎo)入。有限元模型見圖5、圖6。根據(jù)最大凍結(jié)深度的延遲作用和實地氣候原因，推測凍深應(yīng)大于0.45 m小于1 m，為使計算結(jié)果更為精確，網(wǎng)格采用上密下疏的劃分模式，即0~1 m深處的模型網(wǎng)格劃分應(yīng)更細(xì)密。單元選用PLANE55平面單元，共640個節(jié)點，1103個單元。

圖 5 進行瞬態(tài)溫度場計算時渠道（不鋪設(shè)XPS板）模型有限元網(wǎng)格劃分

圖6 進行瞬態(tài)溫度場計算時渠道（鋪設(shè)XPS板）模型有限元網(wǎng)格劃分

4.3 瞬態(tài)溫度場計算及分析

假設(shè)初始溫度場均勻，即先以年平均氣溫12攝氏度為溫度荷載進行穩(wěn)態(tài)熱分析計算，然后將計算結(jié)果作為瞬態(tài)熱分析的初始溫度場進行瞬態(tài)分析計算，得到其凍深隨時間的變化情況。為了更直觀的顯示仿真計算結(jié)果，對模型分別選取1.800、10.190、10.790即圖1所示的1、2、3三條直線上的點分別作為三組觀測數(shù)據(jù)，其布置位置分別為坡面底部與渠基土交界面處、坡面頂部、襯砌板右邊界處，根據(jù)觀測數(shù)據(jù)比較三組凍深隨時間的變化。

對模型按瞬態(tài)模式每24 h施加一次邊界溫度，模擬從5月25日~第二年5月25日的渠道溫度場分布，并分別繪制對照組和試驗組的渠道凍深隨時間的變化曲線以及渠基土上邊界溫度變化曲線。

4.3.1 對照組瞬態(tài)溫度場計算及分析通過對不施加XPS的渠道模型瞬態(tài)溫度場計算得到如下結(jié)論：

（1）三組觀測數(shù)據(jù)年最低溫度，即約為-18 ℃時對應(yīng)的凍深分別為0.3942 m、0.5254 m、0.4508 m。而因凍脹的延遲作用，隨時間變化所對應(yīng)的年最大凍深分別為0.5280 m、0.7552 m、0.6609 m，其平均最大凍深大約為0.6480 m。

（2）最低環(huán)境溫度大約出現(xiàn)在231 d，即1月11日左右，所對應(yīng)平均凍深約為0.4568 m，這與凍深的一維模型計算結(jié)果大致相符，但最低大氣溫度下所對應(yīng)的凍土深度并不是當(dāng)?shù)氐淖畲髢鐾辽疃龋鋬錾畹陌l(fā)展對于環(huán)境溫度來講有一定的滯后效應(yīng)，且根據(jù)相應(yīng)仿真計算可知，隨著氣溫的變化，渠基土?xí)霈F(xiàn)上部溫度高于0 ℃，中間低于0 ℃，下部高于0 ℃的現(xiàn)象。

（3）由圖7計算結(jié)果可知，直線2所對應(yīng)的凍深為三條直線中最大值，因此選取直線2與渠基土上邊界交界點，見圖8，作為渠基土上邊界溫度隨時間變化的觀測點。

圖 7 對照組凍深隨時間的變化曲線

圖 8 觀測基點

（4）不施加任何保溫措施的渠道模型，圖8所示觀測點溫度隨時間的變化曲線與環(huán)境溫度隨時間的變化曲線（圖9）是大致吻合的。故由于環(huán)境溫度的不斷變化，渠基土將隨時間發(fā)生往復(fù)的凍融循環(huán)，一旦因渠基土凍脹將襯砌板頂起造成襯砌板破壞，產(chǎn)生裂縫的襯砌板在滲漏和凍融的交替作用下極易加劇這種凍脹損害。

圖 9 環(huán)境溫度與對照組渠基土上邊界溫度隨時間的變化曲線

4.3.2 試驗組瞬態(tài)溫度場計算及分析通過對施加不同厚度XPS的渠道模型進行瞬態(tài)溫度場計算：

（1）鋪設(shè)0.01、0.03、0.05 mXPS的最大平均凍深分別為0.4843 m、0.2180 m、0.0333 m，其相對應(yīng)減少的凍深分別為0.1637 m、0.4300 m、0.6147 m，經(jīng)計算可知每0.01 mXPS平均可以減少0.1433 m的凍深。

（2）根據(jù)圖10的計算結(jié)果顯示，鋪設(shè)0.01、0.03、0.05 mXPS的渠道最大凍深均出現(xiàn)在直線2上，因此仍選用直線2與渠基土上邊界的交界點（與上文4.3.1（3）中提到的相同）作為渠基土上邊界溫度隨時間變化的觀測點。

圖 10 試驗組凍深隨時間的變化曲線

（3）由圖11可知，三組試驗組渠基土上邊界最低溫度分別出現(xiàn)在第246 d（1月26日），-12.482 ℃、第261 d（2月10日），-5.479 ℃和第269 d（2月18日），-1.589 ℃。

圖 11 環(huán)境溫度與試驗組渠基土上邊界溫度隨時間的變化曲線

（4）綜合（1）、（2）的結(jié)論分析，鋪設(shè)渠道XPS后0 ℃以下的天數(shù)逐漸減小，其中選用0.05 mXPS時的渠道抗凍效果最佳。與原模型仿真計算結(jié)果進行比對可知，鋪設(shè)XPS后，有效阻止了冷空氣進入土體，提高了渠基土的溫度，其內(nèi)部溫度受環(huán)境溫度的影響大大降低。

5 位移場計算

根據(jù)上文瞬態(tài)溫度場計算結(jié)果可知，位山灌區(qū)二干渠渠道年最大凍深約為0.6480 m，且由上文中一維計算（圖4）可知，當(dāng)環(huán)境溫度為年最低氣溫時，渠基土深度達1 m，土溫已穩(wěn)定。因此更改計算模型，進行位移場計算。

5.1 邊界條件

下邊界施加邊界溫度12 ℃，上邊界施加年最低氣溫-18 ℃。位移邊界條件為渠基土左右邊界無水平方向位移，下邊界無豎直方向位移，渠坡下邊界無水平和豎直位移。

5.2 計算模型

根據(jù)上文中提到的參數(shù)條件，需對土層進行分層建模，并分層定義材料的相關(guān)力學(xué)參數(shù)。更改有限元計算模型下邊界為1 m，混凝土襯砌板厚度仍為0.06 m，擠塑聚苯板鋪設(shè)在混凝土襯砌板下方，計算時將襯砌體與渠基土視為一個整體，僅考慮襯砌板自重。為更準(zhǔn)確反映凍脹量的分布，網(wǎng)格劃分仍采用上密下疏的模式，更改單元為PLANE183單元，共2689個節(jié)點，1198個單元。有限元網(wǎng)格劃分如圖。

圖 12 進行位移場計算時渠道（不鋪設(shè)XPS板）模型有限元網(wǎng)格劃分

5.3 位移場計算及分析

通過對對照組和試驗組分別進行位移場計算，得到其凍脹量變化情況如下：

5.3.1 對照組位移場計算及分析根據(jù)ANSYS計算結(jié)果顯示，此梯形模型的渠底處位移為最大值，近似為39.76 mm，平均凍脹量約為33.05 mm。將渠道視為軸對稱結(jié)構(gòu)，提取對稱軸右側(cè)（即渠底中心到渠底右邊界、右坡面和水平右邊界）凍脹量數(shù)據(jù)進行分析，計算結(jié)果如圖13所示。

圖 13 對照組凍脹位移數(shù)值模擬展開圖

由此計算結(jié)果分析可知，未加混凝土襯砌保護的渠底中部，凍脹量最大，其次為渠坡上邊界與水平邊界連接處，分析可知由于該點同時承受兩個方向的低溫荷載，故此點凍脹量相比其他部位較大。結(jié)合溫度場不同位置凍深計算結(jié)果可知渠底、渠坡、水平邊界的凍脹率分別為6.91%、4.22%、3.80%，即渠道的平均凍脹率約為4.98%。與室內(nèi)計算結(jié)果，見表2，絕對誤差0.14%，相對誤差2.9%，滿足工程精度要求。

5.3.2 試驗組位移場計算與分析如圖14所示，3組試驗組最大凍脹量分別為30.11 mm、15.73 mm、5.67 mm；平均凍脹量分別為25.54 mm、14.31 mm、4.34 mm。與對照組結(jié)果相比，三組試驗組分別減少了7.51 mm，18.74 mm、28.71 mm。故0.01 m擠塑聚苯板平均可減少6.50 mm的凍脹。

圖 14 試驗組凍脹位移數(shù)值模擬展開圖

加入擠塑聚苯板后，渠基土上邊界凍脹量明顯減小，同時隨著XPS厚度增加，凍脹量減小的更多，但減小的幅度有所減緩。

6 結(jié)論

（1）通過對渠道進行全年凍脹模擬可知，年最低氣溫出現(xiàn)在1月11日，而最大凍結(jié)深度出現(xiàn)在3月份左右，最大凍結(jié)深度滯后年最低氣溫約60 d，這與當(dāng)?shù)氐膶嶋H情況相符，說明瞬態(tài)溫度場計算結(jié)果是可靠且可行的。

（2）隨著XPS厚度的增加，渠基土上邊界最低溫度出現(xiàn)的時間與年最低氣溫的時間（1月11日）相比逐漸滯后。當(dāng)施加XPS時，保溫作用使得渠基土上邊界最低溫度出現(xiàn)的時間推遲，由此推測，假如XPS厚度足夠大，渠基土上邊界出現(xiàn)最低溫時，外界環(huán)境溫度可能在0 ℃以上。那么，在整個凍結(jié)期內(nèi)，渠基土上邊界有可能不出現(xiàn)負(fù)溫，即渠基土一直處于正溫環(huán)境中。

（3）通過對施加不同厚度的XPS對凍結(jié)深度以及凍脹量的減小比對可知，0.01 m的XPS平均可減少0.14 m的凍結(jié)深度和6.50 mm的凍脹量。實測數(shù)據(jù)顯示，聊城市環(huán)境溫度為年最低氣溫時，其最大凍結(jié)深度為0.45 m左右，小于上文中提到的直線2對應(yīng)的凍結(jié)深度，同時根據(jù)仿真計算結(jié)果產(chǎn)生的凍深延遲作用可推測，鋪設(shè)0.05 m的XPS已足夠滿足凍脹防治要求。

本文計算結(jié)果與實際工況相比有誤差的原因在于，模擬是在參考相關(guān)工程數(shù)據(jù)的基礎(chǔ)上進行的，如果測得本工程相關(guān)的力學(xué)參數(shù)，則模擬結(jié)果將更加精確。與先前所研究的不同在于，本試驗增加了瞬態(tài)溫度場的計算，即通過全年溫度的變化模擬出渠基土隨外界溫度的變化，總結(jié)其凍結(jié)深度的變化規(guī)律，可為后人的研究作參考。

[1] 宋立元,宗兆博,張欣.東港地區(qū)襯砌渠道凍脹有限元數(shù)值模擬研究[J].人民黃河,2012,34(11):99-103

[2] 汪魁峰,宋立元.興城海濱灌區(qū)襯砌渠道抗凍脹試驗與數(shù)值模擬[J].吉林水利,2016(11): 10-14

[3] 徐麗娟.南水北調(diào)中線工程總干渠渠道襯砌防凍脹試驗效果對比分析[J].水科學(xué)與工程技術(shù),2004(6):20-22

[4] Taber S. The mechanics of frost heaving[J]. The Journal of Geology, 1930,38(4):303-317

[5] Beskow G. Soil freezing and frost heaving with special application to roads and railroads[C]. The Swedish Geological Society Year book no.3, Northwestern University, 1935

[6] Everett DH. The thermodynamics of frost damage to porous solids[J]. Trans Faraday Soc, 1961,57:1541-1551

[7] Miller RD. Ice sandwich; Functional semipermeable membrane[J]. Science, 1970,169(3945):584-585

[8] 肖懷璽.渠道凍脹與防治探析[J].北京農(nóng)業(yè),2013(12):230-230

[9] 馬瑞忠.聚苯乙烯保溫板在渠道防凍中的應(yīng)用[J].水利規(guī)劃與設(shè)計,2016(1):96-97

[10] 何武全,鄭水蓉,沈長越.渠道防滲防凍脹復(fù)合型保溫塑料板的試驗研究[J].排灌機械工程學(xué)報,2012,30(5):553-557

[11] 劉瑞娟,劉振森.建筑新秀——復(fù)合式聚苯乙烯保溫板[J].房材與應(yīng)用,2000(2):95,97

[12] 郭利霞,王正中,李甲林.梯形與準(zhǔn)梯形渠道凍脹有限元分析[J].節(jié)水灌溉,2007(4):44-47

[13] 劉國華.保溫板在我市大型灌區(qū)渠道防凍脹中的應(yīng)用[J].內(nèi)蒙古水利,2018,2(186):74-75

[14] 劉志斌,劉建軍.苯板在渠道防凍脹中的應(yīng)用[J].低溫建筑技術(shù),2012,34(7):142-144

[15] 王正中,蘆琴,郭利霞,等.基于晝夜溫度變化的混凝土襯砌渠道凍脹有限元分析[J].農(nóng)業(yè)工程學(xué)報,2009,25(7):1-7

[16] 銀英姿,申向東,步豐湖,等.凍土區(qū)渠道保溫防凍措施對比研究[J].內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報,2010,31(4):210-213

[17] 劉志斌.寒區(qū)混凝土防滲渠道保溫防凍脹措施保溫層合理厚度研究[D].新疆:石河子大學(xué),2013

[18] 權(quán)曉龍.基于有限元的混凝土襯砌渠道凍脹性能研究[J].地下水,2014(1):153-155

[19] 葉爾肯·巴合提漢.基于ANSYS的混凝土襯砌渠道凍脹破壞有限元分析[J].水利科技與經(jīng)濟,2015,21(5):14-16

[20] 胡亮.擠塑聚苯乙烯保溫板在南水北調(diào)渠道工程中的應(yīng)用研究[J].南水北調(diào)與水利科技,2009,7(6):207-210

[21] 楊育紅,楊永來,孫恒.渠道防凍脹保溫板導(dǎo)熱系數(shù)及其測定因素探討[J].南水北調(diào)與水利科技,2014,12(3):191-194

[22] 賀俊,楊平,何文龍.蘇州地鐵典型土層凍土力學(xué)特性研究[J].水文地質(zhì)工程地質(zhì),2010,37(5):72-76

[23] Miller RD. Freezing and heaving of saturated and unsaturated soils[J]. Highway Research Board, 1972(393):1-11

[24] 沈嬌,李德英,介鵬飛.XPS板的性能分析及應(yīng)用研究[J].節(jié)能,2009,28(8):13-15

[25] 桑穎慧.改性聚苯板與傳統(tǒng)聚苯板的性能特點分析[J].山西建筑,2014(33):100-101

[26] 孫立國,杜成斌,王大鵬.EPS材料力學(xué)性能測試及裝飾構(gòu)件數(shù)值仿真[J].建筑材料學(xué)報,2013,16(4):637-641

The Application of Extruded Polystyrene Boards in Frost Heaving Prevention of Concrete Lining Channel

LU Xiang-yu1, LIU Yang-yu2, BIAN Han-bing1,3*, SUN Zhao-hui1, WANG Chen-yu1, QIU Xiu-mei1

1.2710002276000,3.

In view of freeze-thaw damage on concrete-lined channel in Weishan Irrigation District, Liaocheng, Shandong located in the seasonal permafrost region, ANSYS finite element software was utilized to establish a concrete-lined channel model laid with extruded polystyrene (XPS) board of varying thicknesses in order to figure out the transient temperature field distribution and displacement field distribution. Corresponding measured data were also combined to facilitate the analysis. The results indicate after one-year persistent observation and measurement, the appearance of maximum frozen depth lagged somewhat behind the annual lowest temperature arrival no matter whether XPS board was laid when the daily lowest temperature changed over the time. A comparative analysis of frozen volume from channel laid with varying-thickness XPS board as well as frozen depth of soil at the bottom of the channel revealed the insulating mechanism of XPS board. It was assumed that selection of 5cm-thick XPS board to deal with local frost heaving problem could meet the engineering requirements. In this study, transient temperature field calculation was added. This plan is of certain strengths and it can shed some light on future research about frost heaving prevention and control in concrete-lined channels.

Concrete lining channel; frost heaving prevention; numerical simulation; extruded polystyrene board (XPS)

TU757.4

A

1000-2324(2019)03-0460-08

10.3969/j.issn.1000-2324.2019.03.023

2018-05-20

2018-06-04

山東省重點研發(fā)計劃項目:渠道襯砌凍融破壞機理與工程防治措施研究(2017GSF16104)

鹿翔宇(1995-),女,碩士研究生,主要從事渠道襯砌凍融破壞研究. E-mail:lxy19951025@163.com

Author for correspondence. E-mail:hanbing·bian@univ-lorraine.fr; qxmxr@126.com