考慮地下水位影響的現澆混凝土梯形渠道凍脹破壞力學分析

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

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

·農業水土工程·

考慮地下水位影響的現澆混凝土梯形渠道凍脹破壞力學分析

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

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

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

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

0 引 言

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

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

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

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

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

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

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

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

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

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

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

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

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

1.2 梯形渠道基土凍結深度沿斷面分布規律

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

式中H(z)為凍結深度，cm；H0為標準凍深，cm；φ1為土質影響系數；φ2為考慮地域差異（如具體地點的經度、緯度和海拔等）的影響系數；φ3(z)為與計算點的地下水位z有關的影響系數。φ3(z)由式（4）計算。

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

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

在特定地區的特定氣象、水分、土質條件下，襯砌各點基土凍結深度僅與各點的地下水位有關，即式中各影響系數沿斷面分布為定值。結合式（3）、式（4）可得基土凍結深度沿斷面的分布規律如下

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

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

2.1 襯砌結構法向凍脹力沿斷面分布規律

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

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

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

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

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

2.2 襯砌結構凍脹破壞特征和截面內力

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

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

通常襯砌某一截面發生破壞即意味著整塊板的凍脹破壞。因此，凍脹破壞驗算應首先分析彎矩沿斷面的分布規律，進而確定危險截面，并計算出最大彎矩。結合圖1，由材料力學方法計算渠道坡板各截面彎矩并化簡得其沿斷面的分布規律如式（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，這與假設坡板為簡支梁及已有文獻結果一致[3]。式中受特定地區特定土質條件影響的經驗系數a、b及凍土彈性模量Ef反映了基土力學特性對截面彎矩的影響。

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

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

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

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

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

3.1 整體凍脹受力不均勻性

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

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

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

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

3.2 局部凍脹受力不均勻性

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

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

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

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

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

4 原型渠道計算與結果分析

4.1 原型渠道概況

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

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

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

4.2 渠道基土凍深沿斷面的分布規律

依據文獻確定系數[16,29-30]：標準凍結深度取80 cm；渠基土質為壤土，φ1可取1.2；考慮渠道走向和寬深比的修正后，φ2可取1.06；土質為壤土時α取0.79，且當z0≤3 m時取β=z0，當z0>3 m時取β=3。綜上由式（5）可對不同地下水埋深z0時的渠道基土凍深沿渠坡的分布規律進行計算（圖4）。

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

如圖4，渠頂受地下水影響較小，不同地下水埋深z0時該處凍深相差不大，而在渠坡中下部，z0對各點凍深有顯著影響。對z0固定的特定渠道，渠坡各點凍深存在明顯橫向差異，且隨z0增大，差異減小且逐漸趨于一致，即地下水埋深z0越大，其對渠道基土凍結過程的影響越不顯著，這與事實相符。由此可見，對高地下水位渠道如僅按鄰近氣象臺地下水位（相當于渠頂地下水位）計算凍深，將與實際產生較大的偏差，而當地下水深埋時偏差不顯著。

4.3 法向凍脹力沿斷面的分布規律

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

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

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

4.4 截面彎矩沿斷面的分布規律與最大彎矩

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

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

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

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

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

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

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

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

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

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

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

6 結論和討論

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

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

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

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

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

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

[3] 王正中. 梯形渠道砼襯砌凍脹破壞的力學模型研究[J]. 農業工程學報，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] 王正中，李甲林，陳濤，等. 弧底梯形渠道砼襯砌凍脹破壞的力學模型研究[J]. 農業工程學報，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] 申向東，鄭玉佩，王麗萍. 混凝土預制板襯砌梯形斷面渠道的凍脹破壞受力分析[J]. 農業工程學報，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]. 農業工程學報，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]. 土木工程學報，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] 陳立杰，王正中，劉旭東，等. 高地下水位灌排渠道襯砌結構抗凍脹數值模擬[J]. 長江科學院院報，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] 安鵬，邢義川，張愛軍. 基于部分保溫法的渠道保溫板厚度計算與數值模擬[J]. 農業工程學報，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] 王希堯. 不同地下水埋深和不同土壤條件下的凍結和凍脹試驗研究[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] 陳肖柏. 土凍結作用研究近況[J]. 力學進展，1991，21(2)：226－235.

[13] 李安國，李浩，陳清華. 渠道基土凍脹預報的研究[J]. 西北水資源與水工程學報，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] 田亞護. 動、靜荷載作用下細粒土凍結時水分遷移與凍脹特性實驗研究[D]. 北京：北京交通大學，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]. 巖石力學與工程學報，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] 陳肖柏，劉建坤，劉鴻緒，等. 土的凍結作用與地基[M].北京：科學出版社，2006.

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

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

[19] 山西省渠道防滲工程技術手冊編委會. 山西省渠道防滲工程技術手冊[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] 黃繼輝，夏才初，韓常領，等. 考慮圍巖不均勻凍脹的寒區隧道凍脹力解析解[J]. 巖石力學與工程學報，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] 馬巍，王大雁. 凍土力學[M]. 北京：科學出版社，2014.

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

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

[28] 肖旻，王正中，劉銓鴻，等. 開放系統預制混凝土梯形渠道凍脹破壞力學模型及驗證[J]. 農業工程學報，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] 中華人民共和國水利部. 渠系工程抗凍脹設計規范：SL23-2006[S]. 北京：中國水利水電出版社，2006.

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

[31] 木下誠一著. 凍土物理學[M]. 王異，張志權，譯. 長春：吉林科學技術出版社，1985.

[32] 肖旻. 塔里木灌區防滲渠道破壞機理及綜合防治措施研究[D]. 阿拉爾：塔里木大學，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

肖 旻，王正中，劉銓鴻，王 羿，葛建銳. 考慮地下水位影響的現澆混凝土梯形渠道凍脹破壞力學分析[J]. 農業工程學報，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

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

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

Email：xmhdts@163.com

※通信作者：王正中，男，教授，博士生導師，主要從事水工結構工程及凍土工程學科研究。楊凌 西北農林科技大學水利與建筑工程學院，712100。Email：wangzz0910@163.com