肖 旻，王正中,2※，劉銓鴻，王 羿，葛建銳

（1. 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，楊凌 712100；2. 中國科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所凍土工程國家重點實驗室，蘭州730000）

·農(nóng)業(yè)水土工程·

考慮地下水位影響的現(xiàn)澆混凝土梯形渠道凍脹破壞力學(xué)分析

肖 旻1，王正中1,2※，劉銓鴻1，王 羿1，葛建銳1

（1. 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，楊凌 712100；2. 中國科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所凍土工程國家重點實驗室，蘭州730000）

地下水的補給與遷移是高地下水位渠道的凍脹破壞的主要影響因素。該文提出了一種考慮地下水位影響的梯形渠道襯砌凍脹力分布計算方法，推導(dǎo)出地下水位影響的渠道基土凍脹強度和凍結(jié)深度分布的計算公式，并得到現(xiàn)澆混凝土襯砌的截面最大彎矩和最易破壞截面位置的解析表達(dá)式。從整體與局部2個方面定量分析梯形渠道襯砌凍脹力分布的不均勻性，為渠道的抗凍性能評價和斷面優(yōu)化提供了定量指標(biāo)，結(jié)果表明：渠深越淺，坡板傾角越小，凍脹力分布越均勻，越不易發(fā)生破壞，揭示了寬淺式梯形渠道抗凍性能良好的原因。以塔里木灌區(qū)某梯形渠道為原型，對不同地下水埋深的渠道凍脹特征和受力進(jìn)行了分析，并與觀測資料進(jìn)行了對比，其中基土凍深的計算值與觀測值之間的最大相對誤差為3.5%，估算最大彎矩所在截面的位置為距離坡頂63.9%坡板長處，與灌區(qū)實地調(diào)查結(jié)果基本相符，表明了方法的實用性和合理性。最后，對高地下水位梯形渠道的凍害機理進(jìn)行了分析，該研究可為高地下水位現(xiàn)澆混凝土梯形渠道襯砌的抗凍設(shè)計和相關(guān)研究提供參考。

水分；渠道；凍土；凍脹；高地下水位；塔里木灌區(qū)；凍害機理

0 引 言

在中國廣大季節(jié)性凍土區(qū)，灌溉渠道經(jīng)常遭受嚴(yán)重凍害，破壞形式主要表現(xiàn)為襯砌的凍脹破壞[1]。新疆塔里木灌區(qū)以阿拉爾市為中心[2]，年最低氣溫?24～?29.3 ℃，已修渠道2 355 km，地表水豐沛，有阿克蘇新大河、和田河等五大河流貫穿，地下水為河流兩岸嵌入式水體，埋深僅2～3 m，渠道襯砌存在嚴(yán)重的凍脹破壞。由于雨量稀少且地下水埋深淺，引發(fā)凍脹的主要水分來源是地下水，這類地區(qū)渠道凍脹破壞規(guī)律的研究必須考慮地下水的影響。

王正中等[3-4]對現(xiàn)澆梯形、弧底梯形渠道，申向東等[5]對預(yù)制板梯形渠道，宋玲等[6]對冬季輸水條件下的梯形渠道進(jìn)行凍脹力學(xué)分析，建立了一系列凍脹破壞工程力學(xué)模型，為工程設(shè)計提供了簡明、實用的設(shè)計方法。但此類模型多基于工程經(jīng)驗預(yù)先假定凍脹力分布由靜力平衡條件求解，且未定量考慮地下水位的影響。渠道凍脹數(shù)值仿真可以考慮地下水位的影響[7-10]，但因其分析過程復(fù)雜繁瑣，不便于廣泛推廣應(yīng)用，工程實踐中迫切需要既能反映本質(zhì)，又簡捷、實用的渠道凍脹破壞設(shè)計計算方法。

眾多學(xué)者對考慮地下水影響的土體凍脹特性進(jìn)了研究[11-15]。中國北部如甘肅、新疆等大部分省區(qū)的水利、道路部門均設(shè)置大型凍脹試驗場觀測地下水對土體凍脹的影響[16-19]。Stephen[20]通過土體凍結(jié)試驗指出其他條件相同時，地下水埋深不同將引起土體不均勻凍脹；Michalowski等[21]基于孔隙率模型提出采用增長張量反映凍脹變形不均勻性的模型；王正中[3]認(rèn)為襯砌各點地下水位和凍結(jié)鋒面方向不同必然引起不均勻凍脹變形，導(dǎo)致結(jié)構(gòu)破壞；黃繼輝等[22]對寒區(qū)隧道的研究表明不均勻凍脹是產(chǎn)生凍脹應(yīng)力并導(dǎo)致結(jié)構(gòu)破壞的重要原因。

目前渠道凍脹的力學(xué)分析和數(shù)值模擬對表層環(huán)境因素（如太陽輻射、輸水條件等）的關(guān)注較多，考慮地下水影響的相對較少。工程中則僅按鄰近氣象站的地下水埋深進(jìn)行抗凍分析和設(shè)計，未考慮距地下水位較近的點由于地下水補給而導(dǎo)致凍脹的加劇，也未考慮各點至地下水位的距離不同導(dǎo)致的不均勻凍脹力和凍脹變形。這兩方面對高地下水位以及挖方渠道的影響尤為顯著，需更多地加以關(guān)注。

基于此，本文提出一種考慮地下水位影響的梯形渠道凍脹破壞簡明計算方法，并給出反映渠道基土凍脹強度和凍結(jié)深度不均勻分布的計算式，進(jìn)一步分析了襯砌的凍脹不均勻性。結(jié)合原型渠道對其在不同地下水位時的凍脹力學(xué)特性進(jìn)行了對比分析，最后闡明了高地下水位梯形渠道的凍害機理。

1 高地下水位渠道基土凍脹強度和凍深沿斷面的分布

Stephen等[23-24]對有無地下水補給的土體凍結(jié)進(jìn)行對比試驗，結(jié)果表明有補給時的土體凍脹強度比無補給時大得多。高地下水位渠道在特定地區(qū)的特定氣象、水分、土質(zhì)條件下，地下水補給強度是土體凍脹強度的主要影響因素；由于水分相變釋放潛熱，同時水分重分布和冰相的出現(xiàn)導(dǎo)致土體熱容量、導(dǎo)熱率等熱力學(xué)性質(zhì)發(fā)生改變，地下水補給強度對土體凍結(jié)深度也有顯著影響[25-26]。

結(jié)合已有研究成果和工程實踐[3-5]，作如下假設(shè)：1）襯砌渠道縱向尺寸遠(yuǎn)大于橫向尺寸，對渠道凍脹特征的分析簡化為二維平面應(yīng)變問題；2）本文重點針對高地下水位渠道研究考慮地下水位影響的襯砌凍脹破壞問題，地下水補給是主要影響因素，暫不考慮陰陽坡差異；3）渠道基土力學(xué)特性符合Winkler假設(shè)[27-28]，即特定地區(qū)特定氣象、水分、土質(zhì)條件下，襯砌各點對應(yīng)基土凍脹強度由其至地下水位的距離決定；4）計算均以坡板為例，這是因為渠道底板各點至地下水位的距離相同，凍脹力均勻分布，與已有文獻(xiàn)[3-4]一致，故不再贅述。

1.1 梯形渠道基土凍脹強度沿斷面分布規(guī)律

大量文獻(xiàn)和試驗研究表明[11-19]，土體凍脹強度（即凍脹率）與地下水位之間呈負(fù)指數(shù)關(guān)系。

式中η(z)為凍脹強度，%；z為計算點至地下水位的距離，cm；a、b為與特定地區(qū)特定氣象、水分、土質(zhì)條件有關(guān)的經(jīng)驗系數(shù)，常根據(jù)試驗數(shù)據(jù)由最小二乘法擬合。以甘肅省水利廳在張掖試驗場的試驗為例[18]，土質(zhì)為壤土?xí)r，a可取60.05，b可取0.0146。

把式（1）應(yīng)用于梯形襯砌渠道，在特定地區(qū)的特定氣象、水分、土質(zhì)條件下，系數(shù)a、b為定值，可得反映凍脹強度不均勻分布的計算公式如下

式中x為計算點沿坡面至坡頂?shù)木嚯x，cm；η(x)為各點基土凍脹強度，%；z(x)為各點至地下水位的距離，cm。

1.2 梯形渠道基土凍結(jié)深度沿斷面分布規(guī)律

根據(jù)相關(guān)專著和規(guī)范[16,29-30]，至地下水位的距離為z時的土體凍結(jié)深度H(z)可由式（3）計算。

式中H(z)為凍結(jié)深度，cm；H0為標(biāo)準(zhǔn)凍深，cm；φ1為土質(zhì)影響系數(shù)；φ2為考慮地域差異（如具體地點的經(jīng)度、緯度和海拔等）的影響系數(shù)；φ3(z)為與計算點的地下水位z有關(guān)的影響系數(shù)。φ3(z)由式（4）計算。

式中α為與土質(zhì)有關(guān)的經(jīng)驗系數(shù)，β為與地下水埋深（至渠頂）有關(guān)的經(jīng)驗系數(shù)。式（3）、式（4）中各影響系數(shù)可參照相關(guān)專著和規(guī)范取值[16,29-30]。

目前寒區(qū)工程抗凍設(shè)計在計算凍深時僅按鄰近氣象站地下水埋深考慮地下水的影響。在渠道工程中也沿用該方法，而在后文的分析中將指出，對地下水深埋的渠道該法的計算結(jié)果較準(zhǔn)確，而對高地下水位渠道將產(chǎn)生較大偏差。

在特定地區(qū)的特定氣象、水分、土質(zhì)條件下，襯砌各點基土凍結(jié)深度僅與各點的地下水位有關(guān)，即式中各影響系數(shù)沿斷面分布為定值。結(jié)合式（3）、式（4）可得基土凍結(jié)深度沿斷面的分布規(guī)律如下

式中H(x)為斷面各點基土凍結(jié)深度，cm。

2 高地下水位梯形渠道襯砌法向凍脹力分布和破壞類型

2.1 襯砌結(jié)構(gòu)法向凍脹力沿斷面分布規(guī)律

渠道襯砌發(fā)生凍脹位移的同時會一定程度上釋放作用在板上的凍脹力，其削減程度與凍脹位移的大小有關(guān)。為簡化計算，本文假定自由凍脹量Δh(x)完全約束，暫不考慮凍脹力削減，這是偏安全的。自由凍脹量Δh(x)（cm）由式（6）計算。

由木下誠一提出的凍脹力與凍脹率的線性函數(shù)關(guān)系[16-17,31]，結(jié)合式（2），可得襯砌所受法向凍脹力沿斷面的分布規(guī)律如下

式中q(x)為襯砌所受法向凍脹力分布，MPa；Ef為凍土彈性模量，MPa。式（7）具有通用性，對不同斷面型式的渠道，適當(dāng)建立坐標(biāo)系后把斷面上各點至地下水位的距離函數(shù)z(x)代入即可。僅考慮梯形渠道，視坡板為簡支梁[3-5]（圖1），坡頂受渠基土的法向凍結(jié)約束，坡腳則還與底板互相約束。由幾何關(guān)系有式（8）成立。

式中z0為渠道頂部至地下水位的距離，cm；θ為渠道坡板的傾角（圖3）。把式（8）代入式（7），得梯形渠道坡板法向凍脹力分布如下

圖1 渠道坡板受力的計算簡圖Fig.1 Calculating diagram of force situation of canal slope lining plate

2.2 襯砌結(jié)構(gòu)凍脹破壞特征和截面內(nèi)力

結(jié)合筆者對塔里木灌區(qū)渠道凍脹破壞狀況的調(diào)查發(fā)現(xiàn)[32]，局部彎矩過大導(dǎo)致結(jié)構(gòu)軸向拉裂、鼓脹或兩者共同作用導(dǎo)致襯砌板的折斷，是灌區(qū)渠道襯砌凍脹破壞的主要類型（圖2）。由此可知，對此類高地下水位地區(qū)，襯砌截面彎矩計算尤為重要。

圖2 局部彎矩過大導(dǎo)致的襯砌板凍脹破壞Fig.2 Frost heave damage of lining caused by large moment

通常襯砌某一截面發(fā)生破壞即意味著整塊板的凍脹破壞。因此，凍脹破壞驗算應(yīng)首先分析彎矩沿斷面的分布規(guī)律，進(jìn)而確定危險截面，并計算出最大彎矩。結(jié)合圖1，由材料力學(xué)方法計算渠道坡板各截面彎矩并化簡得其沿斷面的分布規(guī)律如式（10）所示。

式中f(bh)=(ebh?1)/(bh)、f(xbsinθ)=(exbsinθ?1)/(xbsinθ)，下同；k1=(l·Ef)/(b·h)；h為渠道深度，cm；M(x)為截面彎矩，kN·m?？紤]到h=l·sinθ，有M(0)=M(l)=0，這與假設(shè)坡板為簡支梁及已有文獻(xiàn)結(jié)果一致[3]。式中受特定地區(qū)特定土質(zhì)條件影響的經(jīng)驗系數(shù)a、b及凍土彈性模量Ef反映了基土力學(xué)特性對截面彎矩的影響。

由式（10）可知M(x)為連續(xù)函數(shù)，故其沿斷面分布必然存在最大值。由數(shù)學(xué)分析方法，可得最大彎矩作用截面（即最易破壞截面）xmax為

式中P(x)為截面剪力，MPa。

地下水位是本文考慮的主要因素，故暫不考慮陰陽坡差異。事實上，上述公式都針對單塊坡板進(jìn)行計算，只需由陰、陽坡太陽輻射差異導(dǎo)致的凍土層冬季最低溫度不同，適當(dāng)調(diào)整相關(guān)參數(shù)（如凍土彈性模量等）就可分別應(yīng)用于陰坡和陽坡的計算。

3 高地下水位梯形渠道凍脹不均勻性分析

渠道襯砌所受凍脹力的分布特征反映其凍脹適應(yīng)性的好壞，通常凍脹力分布越均勻，變化趨勢越平緩，表明凍脹適應(yīng)性越好。以下從整體、局部2個方面對高地下水位梯形渠道襯砌凍脹力分布的不均勻性即凍脹適應(yīng)性進(jìn)行分析。

3.1 整體凍脹受力不均勻性

均方差S(x)通常用來衡量渠道襯砌凍脹力分布的不均勻性。S(x)越小從而襯砌各點凍脹力相對平均值的偏離程度越小，表明凍脹力分布越均勻。由式（9）可知，渠道坡板所受凍脹力的平均值即凍脹力分布的數(shù)學(xué)期望E[q(x)]可由式（14）計算。

可見渠深h越大，地下水埋深z0越小，渠道坡板所受平均凍脹力越大，即整體凍脹強度越大。由上式，襯砌凍脹力分布的均方差可由式（15）表示。

為便于分析，把上式根據(jù)泰勒級數(shù)展開并略去高階項后簡化，得式（16）。

可見對地下水埋深z0固定的特定地區(qū)，渠深h越淺，均方差越小，渠道襯砌所受凍脹力整體分布越均勻。

3.2 局部凍脹受力不均勻性

渠道坡板上兩點間凍脹力（或凍脹量）的差值U可用于衡量凍脹力局部的不均勻性，由式（17）計算

式（17）中令Δx→0，可得渠道坡板某一特定點的凍脹不均勻性，可由式（18）計算。

對于z0固定的特定地區(qū)，坡板傾角θ越小即開口越寬，dq/dx越小，凍脹力隨x的變化趨勢越平緩。以上分析結(jié)果表明，渠深h越淺，渠道開口越寬（坡板傾角θ越小），襯砌所受凍脹力分布越均勻，變化趨勢越平緩，越不易受凍脹破壞，這正是工程實踐中寬淺式梯形渠道抗凍性能良好的原因。

如果考慮h→0 cm且θ→0的極端情況，渠道坡板退化為兩端無約束的有限長平直梁，法向凍脹力分布的均方差為0，分布均勻且無內(nèi)力產(chǎn)生，顯然不會發(fā)生凍脹破壞，與事實相符。此外，不僅新疆塔里木灌區(qū)，近年來對陜西石堡川、寶雞峽和涇恵渠等灌區(qū)渠道凍脹破壞狀況的調(diào)查結(jié)果也表明，寬淺式襯砌渠道的確具有良好的凍脹適應(yīng)性能。

需要指出的是，由于實際工程設(shè)計中襯砌渠道斷面型式需滿足一定的設(shè)計要求（如渠道邊坡穩(wěn)定、斷面設(shè)計流量和工程造價等），從而h和θ的取值是受約束的，其最優(yōu)取值是一個多目標(biāo)非線性的約束優(yōu)化問題，這為綜合考慮經(jīng)濟性和抗凍性的寒區(qū)襯砌渠道斷面優(yōu)化設(shè)計提供了參考。

4 原型渠道計算與結(jié)果分析

4.1 原型渠道概況

以新疆塔里木灌區(qū)某梯形渠道為原型（圖3，因?qū)ΨQ性僅繪制一側(cè)），C20混凝土襯砌，筆者于2010年－2011年越冬期進(jìn)行了原型觀測。凍深采用溫度剖面法即分層布設(shè)熱電偶測量溫度剖面，以0 ℃位置為凍深。凍土層冬季最低溫度約?12 ℃，板厚為8 cm，地下水埋深z0為3 m，土質(zhì)為壤土。

圖3 原型渠道尺寸與觀測點Fig.3 Cross-section dimension and observation points of trapezoidal canal

分別假定地下水埋深z0為2、2.5、3、3.5和4 m，對不同地下水埋深時梯形渠道的凍脹特征和受力進(jìn)行分析。

4.2 渠道基土凍深沿斷面的分布規(guī)律

依據(jù)文獻(xiàn)確定系數(shù)[16,29-30]：標(biāo)準(zhǔn)凍結(jié)深度取80 cm；渠基土質(zhì)為壤土，φ1可取1.2；考慮渠道走向和寬深比的修正后，φ2可取1.06；土質(zhì)為壤土?xí)rα取0.79，且當(dāng)z0≤3 m時取β=z0，當(dāng)z0>3 m時取β=3。綜上由式（5）可對不同地下水埋深z0時的渠道基土凍深沿渠坡的分布規(guī)律進(jìn)行計算（圖4）。

圖4 不同地下水埋深時渠坡凍結(jié)深度分布圖Fig.4 Distribution of frozen depth of canal slope with different groundwater tables

如圖4，渠頂受地下水影響較小，不同地下水埋深z0時該處凍深相差不大，而在渠坡中下部，z0對各點凍深有顯著影響。對z0固定的特定渠道，渠坡各點凍深存在明顯橫向差異，且隨z0增大，差異減小且逐漸趨于一致，即地下水埋深z0越大，其對渠道基土凍結(jié)過程的影響越不顯著，這與事實相符。由此可見，對高地下水位渠道如僅按鄰近氣象臺地下水位（相當(dāng)于渠頂?shù)叵滤唬┯嬎銉錾?，將與實際產(chǎn)生較大的偏差，而當(dāng)?shù)叵滤盥駮r偏差不顯著。

4.3 法向凍脹力沿斷面的分布規(guī)律

為確定式（9）中各參數(shù)，條件具備時應(yīng)通過現(xiàn)場試驗獲取數(shù)據(jù)并由最小二乘法擬合。本文由文獻(xiàn)確定參數(shù)：Ef按凍土層在冬季最低溫度時的凍土彈性模量取為3.82 MPa[16]，這是偏安全的；土質(zhì)為壤土?xí)r[2,13,16]，a可取44.33，b可取0.011。綜上由式（9）可對不同地下水埋深z0時的渠道基土凍結(jié)深度沿渠坡的分布規(guī)律進(jìn)行計算（圖5）。

圖5 不同地下水埋深時渠坡法向凍脹力分布圖Fig.5 Distribution of normal frost heave force of canal slope with different groundwater table

如圖5，地下水埋深z0越小，襯砌整體受力越大，且分布的橫向差異越顯著；當(dāng)z0越大，各點凍脹力變化趨于平緩。對于z0固定的特定渠道，法向凍脹力沿坡板呈指數(shù)規(guī)律增長，且距地下水位越近時趨于線性規(guī)律，這種凍脹力趨于線性分布的現(xiàn)象，在土體分散性越強時表現(xiàn)得越明顯[16-18]。實際上，把式（9）按泰勒級數(shù)展開取一級近似，即得凍脹力分布的線性分布規(guī)律，這也與文獻(xiàn)[3-4]中由工程實踐經(jīng)驗預(yù)先假定的凍脹力分布規(guī)律相一致。

4.4 截面彎矩沿斷面的分布規(guī)律與最大彎矩

由式（10）可對襯砌截面彎矩沿渠坡的分布規(guī)律進(jìn)行計算（圖6）。由圖6可知，不同地下水埋深z0對截面彎矩變化的總體趨勢影響較小，但對截面彎矩大小尤其是最大彎矩的量值影響顯著。隨著地下水埋深越淺，截面最大彎矩迅速增大。為由式（12）可知，對于特定斷面的梯形渠道，最大彎矩隨地下水埋深z0越小，呈指數(shù)規(guī)律增大，即寒區(qū)高地下水位渠道極易遭受凍脹破壞，與事實相符。

圖6 不同地下水埋深時渠坡截面彎矩分布圖Fig.6 Distribution of section bending moment of canal slope with different groundwater table

4.5 原型渠道驗證與誤差分析

如圖3所示，在實際地下水埋深時（即z0=3 m），各觀測點基土凍深計算值與觀測值基本相符，表明方法的合理性。在對基土做分層級配分析時，發(fā)現(xiàn)土層中夾雜部分砂土、砂壤土層，而本文均按壤土層進(jìn)行參數(shù)選取和計算，導(dǎo)致了一定誤差（即計算值相對觀測值普遍偏?。?，但誤差不顯著，本文方法的計算值與觀測值之間的最大相對誤差為3.5%。

結(jié)合式（11），可對渠道襯砌最易拉裂、鼓脹和折斷的部位即最大彎矩截面進(jìn)行估算。

將相關(guān)參數(shù)取值和原型斷面尺寸代入式（19），得xmax/l約為63.9%。調(diào)查表明灌區(qū)渠道坡板的凍脹破壞多發(fā)生在中下部距離坡頂約60%～80%坡板長處，少量發(fā)生上半部，計算結(jié)果與工程實際相符。

由于凍土物理力學(xué)性質(zhì)和水分遷移、相變的復(fù)雜性，目前尚缺乏完備的室內(nèi)模型試驗數(shù)據(jù)可供驗證，筆者結(jié)合塔里木灌區(qū)渠道凍脹破壞情況的實地調(diào)查及灌區(qū)某梯形渠道的原型觀測，對部分計算結(jié)果進(jìn)行對比驗證，結(jié)果與原型觀測基本相符，表明本文方法具有一定的實用性和合理性。更深入系統(tǒng)地分析和驗證還有待進(jìn)一步研究。

5 高地下水位梯形渠道襯砌的凍害機理

有研究表明[33-34]，土體凍結(jié)較慢且水分來源充足時，將引發(fā)劇烈凍脹。在中國廣大季節(jié)性凍土區(qū)，冬季漫長而氣溫下降緩慢，渠道基土凍結(jié)速率也較慢，水分有充足的時間遷移和補給；同時高地下水位渠道地下水埋深淺，又正好為凍結(jié)過程提供了充足的水分來源。充足的水分遷移時間和補給來源是高地下水位渠道襯砌凍脹破壞的主要原因之一。

由于渠道特殊的斷面特性，襯砌各點至地下水位的距離不同導(dǎo)致地下水補給強度沿斷面不均勻分布，引起基土凍脹強度和凍結(jié)深度沿斷面的橫向差異，進(jìn)而導(dǎo)致襯砌所受凍脹力和凍脹變形沿斷面的不均勻分布。渠道襯砌凍脹力分布不均勻是高地下水位渠道襯砌凍脹破壞的另一個主要原因。圖7簡要概括了高地下水位梯形渠道的凍害機理。

圖7 高地下水位梯形渠道襯砌凍害機理Fig.7 Mechanism of frost heave damage of trapezoidal canal with high groundwater table

6 結(jié)論和討論

1）該文提出一種考慮地下水位影響的梯形渠道襯砌凍脹力分布計算方法，推導(dǎo)出地下水位影響的渠道基土凍脹強度和凍結(jié)深度不均勻分布的計算公式，并得到截面最大彎矩和最易破壞截面位置的解析表達(dá)式。

2）從整體與局部2個方面定量分析梯形渠道襯砌凍脹力分布的不均勻性，為渠道的抗凍性能評價和斷面優(yōu)化提供了定量指標(biāo)。分析指出，渠深越淺，坡板傾角越小，凍脹力分布越均勻，變化趨勢越平緩，即越不易發(fā)生凍脹破壞，這正是寬淺式梯形渠道抗凍性能良好的主要原因，與實際相符。

3）以塔里木灌區(qū)某梯形渠道為原型，對不同地下水埋深的渠道凍脹特征和受力進(jìn)行了分析，并與觀測資料進(jìn)行對比，其中基土凍深的計算值與觀測值之間的最大相對誤差為3.5%，估算最大彎矩截面的位置為距離坡頂63.9%坡板長處，與灌區(qū)實地調(diào)查結(jié)果基本相符，表明了方法的實用性和合理性。此外，分析結(jié)果還表明，地下水埋深對高地下水位襯砌渠道的凍脹特征和受力有顯著影響，地下水埋深越小渠坡襯砌凍脹力越趨于線性分布，截面最大彎矩則隨地下水埋深越小呈指數(shù)規(guī)律增大。由此可見，高地下水位渠道極易產(chǎn)生嚴(yán)重凍脹破壞。

本文假設(shè)自由凍脹量被完全約束，暫不考慮襯砌凍脹位移引起的凍脹力削減效應(yīng)（這是偏安全的），如需考慮襯砌與凍土間的動態(tài)相互作用及由此引起的凍脹力削減效應(yīng)，可采用彈性地基梁理論進(jìn)行渠道凍脹破壞分析，仍有待進(jìn)一步的研究。

[1] 周幼吾，郭東信，邱國慶，等. 中國凍土[M]. 北京：科學(xué)出版社，2000.

[2] 塔里木灌區(qū)水利管理處史志編撰委員會.塔里木灌區(qū)水利管理處志[M]. 向詩漁，主編. 烏魯木齊：新疆人民出版社，2001.

[3] 王正中. 梯形渠道砼襯砌凍脹破壞的力學(xué)模型研究[J]. 農(nóng)業(yè)工程學(xué)報，2004，20(3)：24－29. Wang Zhengzhong. Establishment and application of mechanics models of frost heaving damage of concrete lining trapezoidal open canal[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2004, 20(3): 24－29. (in Chinese with English abstract)

[4] 王正中，李甲林，陳濤，等. 弧底梯形渠道砼襯砌凍脹破壞的力學(xué)模型研究[J]. 農(nóng)業(yè)工程學(xué)報，2008，24(1)：18－23. Wang Zhengzhong, Li Jialin, Chen Tao, et al. Mechanics models of frost heaving damage of concrete lining trapezoidal canal with arc-bottom[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2008, 24(1): 18－23. (in Chinese with English abstract)

[5] 申向東，鄭玉佩，王麗萍. 混凝土預(yù)制板襯砌梯形斷面渠道的凍脹破壞受力分析[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(16)：80－85. Shen Xiangdong, Zheng Yupei, Wang Liping.Stress analysis of frost heave for precast concrete panel lining trapezoidalcross-section channel[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2012, 28(16): 80－85. (in Chinese with English abstract)

[6] 宋玲，歐陽輝，余書超. 混凝土防滲渠道冬季輸水運行中凍脹與抗凍脹力驗算[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(18)：114－120. Song Ling, Ouyang Hui, Yu Shuchao. Frozen heaving and capacity of frozen heaving resistance of trapezoidal concrete lining canal with water in winter[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2015, 31(18): 114－120. (in Chinese with English abstract)

[7] 王正中，沙際德，蔣允靜，等. 正交各向異性凍土與建筑物相互作用的非線性有限元分析[J]. 土木工程學(xué)報，1999，32(3)：55－60. Wang Zhengzhong, Sha Jide, Jiang Yunjing, et al. Nonlinear finite element analysis of interaction of orthotropic frozen ground and construction[J]. China Civil Engineering Journal, 1999, 32(3): 55－60. (in Chinese with English abstract)

[8] 陳立杰，王正中，劉旭東，等. 高地下水位灌排渠道襯砌結(jié)構(gòu)抗凍脹數(shù)值模擬[J]. 長江科學(xué)院院報，2009，26(9)：66－70. Chen Lijie, Wang Zhengzhong, Liu Xudong, et al. Simulation analysis of channel frost heave under high groundwater lever[J]. Jouranl of Yangtze River Scientific Research Institute, 2009, 26(9): 66－70. (in Chinese with English abstract)

[9] 安鵬，邢義川，張愛軍. 基于部分保溫法的渠道保溫板厚度計算與數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報，2013，29(17)：54－62. An Peng, Xing Yichuan, Zhang Aijun. Thickness calculation and numerical simulation of insulation board for canal using partial insulation method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2013, 29(17): 54－62. (in Chinese with English abstract)

[10] Li Zhuo, Liu Sihong, Feng Youting, et al. Numerical and study on the effect of frost heave prevention with different canal lining structure in seasonally frozen ground regions[J]. Cold Regions Science and Technology, 2013, 85(1): 242－249.

[11] 王希堯. 不同地下水埋深和不同土壤條件下的凍結(jié)和凍脹試驗研究[J]. 冰川凍土，1980，3(2)：40－45. Wang Xiyao. Experimental research of frost heave under different level of groundwater in various soils[J]. Journal of Glaciology and Geocryology, 1980, 3(2): 40－45. (in Chinese with English abstract)

[12] 陳肖柏. 土凍結(jié)作用研究近況[J]. 力學(xué)進(jìn)展，1991，21(2)：226－235.

[13] 李安國，李浩，陳清華. 渠道基土凍脹預(yù)報的研究[J]. 西北水資源與水工程學(xué)報，1993，4(3)：17－23. Li Anguo, Li Hao, Chen Qinghua. Study on the prediction of frost heave in the bedsoil of canals[J]. Water Resources & Water Engineering, 1993, 4(3): 17－23. (in Chinese with English abstract)

[14] 田亞護(hù). 動、靜荷載作用下細(xì)粒土凍結(jié)時水分遷移與凍脹特性實驗研究[D]. 北京：北京交通大學(xué)，2008. Tian Yahu. The Experimental Study on Moisture Migration and Frost action of Fine-grained Soils during Freezing under Dynamic and Static Loading[D]. Beijing: Beijing Jiaotong University, 2008. (in Chinese with English abstract)

[15] 盛岱超，張升，賀佐躍. 土體凍脹敏感性評價[J]. 巖石力學(xué)與工程學(xué)報，2014，33(3)：594－605. Sheng Daichao, Zhang Sheng, He Zuoyue. Assessing frost susceptibility of soils[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(3): 594－605. (in Chinese with English abstract)

[16] 陳肖柏，劉建坤，劉鴻緒，等. 土的凍結(jié)作用與地基[M].北京：科學(xué)出版社，2006.

[17] 李甲林，王正中.渠道襯砌凍脹破壞力學(xué)模型及防凍脹結(jié)構(gòu)[M]. 北京：中國水利水電出版社，2013.

[18] 甘肅省渠道防滲抗凍試驗研究小組. 甘肅省渠道防滲抗凍試驗研究報告[R]. 甘肅省水利廳，1985.

[19] 山西省渠道防滲工程技術(shù)手冊編委會. 山西省渠道防滲工程技術(shù)手冊[M]. 太原：山西科技出版社，2003.

[20] Stephen T. Frost heaving[J].Journal of Geology, 1929, 37(5): 428－461.

[21] Michalowski R L, Zhu Min. Frost heave modelling using porosity rate function[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2006, 30(8): 703－722.

[22] 黃繼輝，夏才初，韓常領(lǐng)，等. 考慮圍巖不均勻凍脹的寒區(qū)隧道凍脹力解析解[J]. 巖石力學(xué)與工程學(xué)報，2015，34(增2)：3766－3774. Huang Jihui, Xia Caichu, Han Changling, et al. Analytical solution of frost heave force action on cold-region tunnel liner considering anisotropy frost heave of surrounding rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(Supp.2): 3766－3774. (in Chinese with English abstract)

[23] Stephen T. The mechanics of frost heaving[J]. Journal of Geology, 1930, 38(4): 303－317.

[24] Leonid B. The modelling of the freezing process in fine-grained porous media: Application to the frost heave estimation[J]. Cold Regions Science and Technology, 2009, 56: 120－134.

[25] 馬巍，王大雁. 凍土力學(xué)[M]. 北京：科學(xué)出版社，2014.

[26] 鄭秀清，樊貴盛，邢述彥. 水分在季節(jié)性非飽和凍融土壤中的運動[M]. 北京：地質(zhì)出版社，2002.

[27] Selvadural A P S. 土與基礎(chǔ)相互作用的彈性分析[M]. 范文田等，譯. 北京：中國鐵道出版社，1984.

[28] 肖旻，王正中，劉銓鴻，等. 開放系統(tǒng)預(yù)制混凝土梯形渠道凍脹破壞力學(xué)模型及驗證[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(19)：100－105. Xiao Min, Wang Zhengzhong, Liu Quanhong, et al. Mechanical model and validation of frost heave damage of precast concrete slab lining trapezoidal canal in open system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CASE), 2016, 32(19): 100－105. (in Chinese with English abstract)

[29] 中華人民共和國水利部. 渠系工程抗凍脹設(shè)計規(guī)范：SL23-2006[S]. 北京：中國水利水電出版社，2006.

[30] 國家能源局. 水工建筑物抗冰凍設(shè)計規(guī)范：NB/T 35024-2014[S]. 北京：中國電力出版社，2015.

[31] 木下誠一著. 凍土物理學(xué)[M]. 王異，張志權(quán)，譯. 長春：吉林科學(xué)技術(shù)出版社，1985.

[32] 肖旻. 塔里木灌區(qū)防滲渠道破壞機理及綜合防治措施研究[D]. 阿拉爾：塔里木大學(xué)，2011. Xiao Min. Studies on Channel Damage Mechanism and its Prevention Measures in Tarim Irrigation[D]. Alar: Tarim University, 2011. (in Chinese with English abstract).

[33] Peppin S S L, Style T W. The physics of frost heave and ice-lens growth[R]. Oxford Centre for Collaborative Applied Mathematics, 2012.

[34] Leonid B, Regina B. Modeling frost heave in freezing soils[J]. Cold Regions Science and Technology, 2010, 61(12): 43－46.

Mechanical analysis on frost heave damage of cast-in-place concrete trapezoidal canal considering influence of groundwater table

Xiao Min1, Wang Zhengzhong1,2※, Liu Quanhong1, Wang Yi1, Ge Jianrui1
(1. College of Water Resources and Architectural Engineering, Northwest A&F University, Yangling 712100, China; 2. State Key Laboratory of Frozen Soil Engineering, Cold and Arid Region Environmental and Engineering Research Institue, CAS, Lanzhou 730000, China)

Groundwater migration and supply is the main influence factor causing frost heave damage of concrete lining trapezoidal canals with high groundwater table. Considering the influence of groundwater table on frost heave of canal foundation soil, the article puts forward a method to calculate the distribution of frost heave force acting on concrete trapezoidal canal lining, presents the design formulas to calculate the distribution of frost heave intensity and frozen depth of canal foundation soil, and then proposes the analytical expression to determine the maximum bending moment of cast-in-place concrete trapezoidal canal lining and the position on canal lining plate which is the easiest to crack. The inhomogeneity of distribution of frost heave force is quantitatively analyzed both locally and integrally, and a new indicator to judge frost heave resistance of canal lining and a method for cross-section pattern optimization are provided. For specific region with fixed groundwater table, the smaller the canal cross-section depth, the more homogeneous the distribution of frost heave force integrally. For specific region with fixed groundwater table from canal top, the smaller the inclination of canal slope plate, the smoother the variation tendency of frost heave force locally. In a word, when the canal cross-section depth is shallower and the inclination of canal slope lining plate is smaller, the distribution of frost heave force is more homogeneous, and the canal lining is harder to be damaged by frost heave force, therefore revealing the reason that the frost heave resistance of trapezoidal canal with wide and shallow cross-section is favorable. By taking a trapezoidal concrete lining canal in Tarim irrigation area as prototype, the frost heave characteristics of canal lining are analyzed, as well as the distribution of frost heave force acting on trapezoidal canal lining with different groundwater table, and the rationality of the method is confirmed by comparing the calculated value with the measured value. The results of comparison show that the maximum relative error of frozen depth of canal foundation soil is 3.5%, and the position of the section undergoing the maximum bending moment estimated is 63.9% of the length of canal slope lining plate from canal top. In addition, the analysis results also show that different groundwater table remarkably influences frost heave characteristics of canal lining and distribution of frost heave force acting on canal lining with high groundwater table. The tendency of variation of frost heave force becomes more linear when groundwater table is higher and the maximum bending moment on canal slope lining plate increases following exponential law. Therefore, the canal with high groundwater table is extraordinarily easy to suffer serious frost heave damage. Eventually, the mechanism of frost heave damage of trapezoidal canal with high groundwater table is analyzed. Because of long and severe winter in North China, in addition to slowly decreasing of air temperature, frozen rate of canal foundation soil is small, and then soil water has sufficient time to migrate and supply. While groundwater table is shallow, groundwater migration and supply is enough to cause serious frost heave damage to canal lining. On account of specific characteristics of canal cross-section, different distances from each position of canal slope plate to groundwater table result in inhomogeneous distribution of the intensity of groundwater supply, and then cause transverse difference of intensity of frost heave and frozen depth of canal foundation soil, and therefore the distribution of frost heave force and the deformation of canal lining can be inhomogeneous. In fact, inhomogeneous distribution of frost heave force acting on canal lining is also a main reason causing serious frost heave damage to canal lining. The research can provide theoretical basis for relevant study and frost heave resistance design of cast-in-place concrete trapezoidal canal with high groundwater table.

moisture; canals; frozen soils; frost heave; high groundwater level; Tarim irrigation area; frost heave damage mechanism

10.11975/j.issn.1002-6819.2017.11.012

S277; TV67

A

1002-6819(2017)-11-0091-07

肖 旻，王正中，劉銓鴻，王 羿，葛建銳. 考慮地下水位影響的現(xiàn)澆混凝土梯形渠道凍脹破壞力學(xué)分析[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(11)：91－97.

10.11975/j.issn.1002-6819.2017.11.012 http://www.tcsae.org

Xiao Min, Wang Zhengzhong, Liu Quanhong, Wang Yi, Ge Jianrui. Mechanical analysis on frost heave damage of cast-in-place concrete trapezoidal canal considering influence of groundwater table[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(11): 91－97. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.11.012 http://www.tcsae.org

2016-11-03

2017-05-02

國家自然科學(xué)基金（51279168）；國家科技支撐計劃（2012BAD10B02）；博士點基金項目（2012020411100254）

肖 旻，男，博士生，主要從事凍土工程及渠道防滲抗凍脹研究。楊凌 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，712100。

Email：xmhdts@163.com

※通信作者：王正中，男，教授，博士生導(dǎo)師，主要從事水工結(jié)構(gòu)工程及凍土工程學(xué)科研究。楊凌 西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，712100。Email：wangzz0910@163.com