鄭 源，湯 驊，姜海波

(石河子大學水利建筑工程學院，新疆石河子 832000)

我國對水資源的需求非常大，而水資源緊缺、浪費現(xiàn)象又十分嚴重，導致我國面臨越來越嚴峻的農(nóng)業(yè)用水危機，節(jié)水措施已成為保持農(nóng)業(yè)持續(xù)發(fā)展所需解決的重要問題［1］。渠道襯砌防滲技術(shù)是我國應(yīng)用最普遍的節(jié)水措施之一，其減少的滲流量高達90%［2］。但是在我國季節(jié)凍土區(qū)存在著土體凍融現(xiàn)象，導致渠道襯砌體因土體凍脹而出現(xiàn)破壞現(xiàn)象，從而降低甚至達不到防滲節(jié)水效果。因此，要重點進行渠道防凍脹研究，關(guān)鍵是要研究渠道襯砌體在凍脹時的受力情況，從力學角度對渠道的凍脹破壞進行分析，揭示渠道凍脹破壞的機理，為渠道防滲、防凍脹技術(shù)提供科學依據(jù)。

研究表明，混凝土襯砌結(jié)構(gòu)產(chǎn)生的凍脹破壞是由襯砌材料抗拉強度低的內(nèi)部原因和凍脹變形被約束的外部原因共同造成的［3-5］。有研究人員分析了混凝土襯砌渠道的力學特性并進行了內(nèi)力計算，但并沒有考慮復合襯砌渠道中土工膜與渠基土壤的相互作用［6-10］。在襯砌結(jié)構(gòu)"可以有一定凍脹位移量"的設(shè)計標準和"削減或消除凍脹"的防凍脹設(shè)計原則上加入復合土工膜［11-12］，在凍脹作用下襯砌體將會在復合土工膜和渠基土壤接觸面產(chǎn)生一定的凍脹位移量，使得襯砌體的變形約束減小，從而達到渠道襯砌抗凍脹破壞的目的。筆者曾分析過U形復合襯砌渠道凍脹破壞的力學特征，但將襯砌坡板看作為簡支梁，沒有考慮到坡板與底板整體性強的特點［13］。在剛性混凝土襯砌體凍脹力學模型基礎(chǔ)之上，建立了復合襯砌渠道弧底梯形斷面的凍脹破壞力學模型，進行內(nèi)力計算，對復合土工膜在襯砌體結(jié)構(gòu)破壞中的影響進行了探討，分析了影響渠道凍脹破壞的因素，揭示了季節(jié)性寒冷凍土區(qū)渠道的凍脹機理。

1 力學模型的建立

在有水分補給的條件下，季節(jié)性寒冷凍土區(qū)的土體將在溫度梯度和水分遷移的影響下發(fā)生凍結(jié)和凍脹，體積隨之增大并產(chǎn)生凍脹量。若無外界約束，凍脹量將能夠自由釋放。但土體的自由凍脹被襯砌體約束，反之土體將會對襯砌體產(chǎn)生約束反力。垂直于襯砌體的凍脹力稱作法向凍脹力，平行于襯砌體的力稱作切向凍結(jié)力。復合襯砌結(jié)構(gòu)在產(chǎn)生一定程度的凍脹位移量后，作用在襯砌體上的有效凍脹力與有效凍結(jié)力將在土工膜與土壤摩擦力的作用下減小，再與重力、底板與坡板相互約束力共同作用達到平衡。根據(jù)以上凍脹破壞的特征及原因，在建立力學模型時作出相應(yīng)的假設(shè)和簡化。

1．1 模型基本假設(shè)及簡化

a．將渠基凍結(jié)土壤和混凝土襯砌板均視為彈性材料，且混凝土板的彈性模量遠大于凍土的彈性模量。

b．渠基土壤的固結(jié)已經(jīng)全部完成，不考慮未凍結(jié)凍土壤的壓縮效應(yīng)。

c．渠道襯砌體的切向凍結(jié)力在弧形底板上線性分布，在弧形底板底部的值為零;渠道襯砌體的法向凍脹力在弧形底板由陰坡向陽坡減小。切向凍結(jié)力和法向凍脹力在坡板上沿著坡板軸線線性分布，在坡板與弧形底板連接處達到最大值，在頂部為零。

d．渠道襯砌體在坡底處坡板與弧形底板相互約束，在坡頂處呈自由變形，將襯砌坡板看作懸臂梁。

e．將混凝土板和復合土工膜看成黏結(jié)的整體，在底板頂推力作用下，其與凍結(jié)土壤之間產(chǎn)生凍脹位移量時存在摩擦力。摩擦力在弧形底板底部的值為零，沿坡板線性分布。

1．2 渠道斷面圖、受力圖、坡板及底板計算簡圖

1．2．1 斷面圖

弧底梯形復合襯砌渠道斷面如圖1所示(設(shè)渠道坡板長為L，弧底中心角為2α，弧底半徑為R，襯砌板厚為b)。

圖1 弧底梯形復合襯砌渠道斷面

1．2．2 受力圖

渠道凍結(jié)后，襯砌體在法向凍脹力、凍結(jié)力、摩擦力、襯砌板的相互約束力及襯砌體的重力作用下保持平衡。渠道坡板陰坡與陽坡的計算簡圖相同，只是其上作用荷載數(shù)值不同(對陽坡的有關(guān)外力及反力在相關(guān)變量符號右上角加撇以示區(qū)別)。渠道受力情況如圖2所示(設(shè)最大法向凍脹力為q0，最大切向凍結(jié)力為τ0，摩擦力為τ1)。根據(jù)分析和假設(shè)得

圖2 弧底梯形渠道受力情況示意圖

渠道在達到凍脹極限平衡狀態(tài)時，渠道陰坡襯砌板與渠基凍結(jié)土壤之間的凍結(jié)力達到最大值τ0，其大小由溫度、土壤含水量、土質(zhì)等因素決定，屬于已知反力。最大摩擦力與最大法向凍脹力線性相關(guān)。根據(jù)假設(shè)并由受力圖建立靜力平衡方程。

渠道受力圖上的未知力只有q0，列出豎向受力平衡方程如下:

其中 τ1= μq0

式中μ為土工膜與土壤摩擦因數(shù)。由式(1)～(3)聯(lián)立求出渠道在凍脹破壞作用下達到極限平衡狀態(tài)時的 q0、τ1，進一步求解各內(nèi)力。

1．2．3 坡板計算簡圖

根據(jù)以上分析，將坡板視為懸臂梁(如圖3所示)，其受力作用包括沿切向線性分布的凍結(jié)力、沿坡板線性分布的法向凍脹力、重力、沿切向分布的摩擦力，以及弧形底板與坡板之間的約束力。RBy、RBx、MB表示弧形底板對坡板的約束力。

圖3 渠道坡板計算簡圖

1．2．4 弧形底板計算簡圖

渠道弧形底板兩端的約束反力因陰坡陽坡往往不相等，這里選取陰坡弧形底板即中心線以左的部分(圖4)進行計算。圖4中N、M、Q表示坡板對底板的約束力，N0、M0、Q0表示求解截面內(nèi)力。

圖4 弧形底板計算簡圖

2 力學模型求解

2．1 渠道坡板內(nèi)力計算

坐標系如圖3所示。

a．支座處反力計算式如下:

b．任意截面的軸壓力計算式如下:

c．任意截面的彎矩計算式如下:

最大彎矩截面

d．任意截面剪力計算式如下:

各內(nèi)力分布如圖5所示。圖5結(jié)果表明:坡板軸力單調(diào)函數(shù)在坡板底部產(chǎn)生最大值;坡板彎矩和剪力在坡板內(nèi)取得極值，極值位置與參數(shù)的取值有關(guān)。圖中規(guī)定正負號按軸力拉正壓負、彎矩下側(cè)受拉為正、剪力以隔離體順時針轉(zhuǎn)動為正。

2．2 渠道底板內(nèi)力計算

以弧底中心處(見圖4)為坐標原點，根據(jù)靜力平衡條件，由∑X=0可得與中心線成θ角處的軸力為

圖5 渠道坡板內(nèi)力示意圖

由∑M=0，可得與中心線成θ角處的彎矩為

由∑Y=0，可得與中心線成θ角處的剪力為

弧形底板內(nèi)力如圖6所示。由圖6可以看出:

圖6 弧形底板內(nèi)力示意圖

弧形底板軸力和剪力為單調(diào)函數(shù);彎矩受拉側(cè)在弧底板內(nèi)表面，體現(xiàn)了弧形底板的反拱作用。圖中規(guī)定正負號按軸力拉正壓負、彎矩下側(cè)受拉為正、剪力以隔離體順時針轉(zhuǎn)動為正。

2．3 混凝土襯砌板厚度驗算

混凝土襯砌板最大彎矩處的最大拉應(yīng)變是否超過允許值決定了襯砌板是否脹裂;在此過程中，將坡板和弧形底板視為壓彎構(gòu)件，一般不考慮剪力的影響。

a．渠道襯砌坡板:坡板最大彎矩處拉應(yīng)力最大，其最大拉應(yīng)力計算式如下:

其中x0按式(9)計算?？沽褩l件驗算:

式中Ec、εt分別為C20混凝土的極限拉應(yīng)變及彈性模量，可查相關(guān)手冊得到。

b．渠道襯砌弧形底板:弧形底板最大彎矩處拉應(yīng)力值最大，其最大拉應(yīng)力計算式如下:

抗裂條件同式(16)。

3 應(yīng)用實例

某弧底梯形渠道，由C20混凝土構(gòu)成，材料密度ρ=2400 kg/m3，邊坡板長L=3.5 m，弧底半徑R=3.24m，襯砌板厚b=0.15m，坡角α=33.69°。渠床土壤為粉質(zhì)壤土，土工膜與土壤摩擦角為5.33°，陰坡、陽坡凍土層最低溫度分別為－15℃、－12℃，需判斷該襯砌結(jié)構(gòu)是否會發(fā)生凍脹破壞。C20混凝土的極限拉應(yīng)變εt=1.0×10－4，彈性模量Ec=260GPa。

3．1 復合襯砌形式結(jié)構(gòu)計算

步驟1 最大切向凍結(jié)力計算。由經(jīng)驗公式得τ0=c+mt= －9.4 kPa(c、m 與土質(zhì)、水分等有關(guān)，c取0.4，m 取0.6)。

步驟2 復合襯砌渠道最大法向凍脹力、摩擦力計算。由式(1)～(3)可得 q0=7.08 kPa，τ1=0.66 kPa。

步驟3 復合襯砌渠道坡板內(nèi)力計算及凍脹破壞判斷。按式(9)計算得x0=3.15 m，按式(7)計算得N(x0)=18.63 kN，按式(8)計算得 M(x0)=5.05 kN·m。按式(15)計算得 σ0=1.22 MPa，而混凝土板的極限拉應(yīng)力 σt=εtEc=2.6 MPa，σ0＜ σt，所以坡板不會發(fā)生凍脹破壞。

步驟4 復合襯砌渠道弧底板內(nèi)力計算及凍脹破壞判斷。經(jīng)計算，最大彎矩發(fā)生在弧形底部，θ0=0°，N(θ0)=30.95 kN，M(θ0)=10.29 kN·m;同理，最大拉應(yīng)力 σ'0=2.43 MPa，σ0＜ σt，所以弧形底板不會發(fā)生凍脹破壞。

3．2 素混凝土板襯砌形式結(jié)構(gòu)計算

步驟1 最大法向凍脹力計算。素混凝土襯砌體在不考慮摩擦力情況下對應(yīng)的凍脹力q01=7.27kN。

步驟2 坡板內(nèi)力計算及凍脹破壞判斷。計算得:最大彎矩截面x01=3.08 m，N(x01)=18.87 kN，M(x01)=5.25 kN·m;最大拉應(yīng)力σ01=1.27 MPa，σ01＜σt，襯砌坡板不會發(fā)生凍脹破壞，與實際情況相符。

步驟3 弧形底板內(nèi)力計算及凍脹破壞判斷。最大彎矩發(fā)生在弧形底部，θ01=0°，N(θ01)=32.76 kN，M(θ01)=9.89 kN·m;同理，最大拉應(yīng)力σ'01=2.53 MPa，σ'01＜σt，弧形底板不會發(fā)生凍脹破壞，與實際情況相符。

3．3 計算結(jié)果分析

計算結(jié)果表明，坡板彎矩最大值在坡板頂部往下2/3處偏下，弧形底板彎矩最大值在中心線附近，則襯砌板容易在此處產(chǎn)生凍脹破壞現(xiàn)象，這與實際工程情況吻合。

復合襯砌渠道與素混凝土板渠道相比，作用于襯砌體的最大有效凍脹力和有效切向約束力分別減小了2.6%和7%，坡板與弧形底板抗拉強度分別提高了3.9%和4%，該結(jié)果表明復合襯砌渠道通過允許一定的凍脹位移量使作用于襯砌體的有效凍脹力和有效切向約束力減小，從而提高抗拉強度，達到防凍脹的目的。

4 結(jié)論

a．渠道的凍脹破壞是在土質(zhì)、水分、溫度等因素共同作用下的復雜過程。通過適當假設(shè)，根據(jù)渠道凍脹破壞的特征和規(guī)律，建立弧底梯形復合襯砌渠道凍脹力學模型，將渠道凍脹破壞表達為最大切向凍結(jié)力的函數(shù)，切向凍結(jié)力可根據(jù)經(jīng)驗或試驗確定，將渠道凍脹破壞過程簡單定量化。

b．加入復合土工膜后，襯砌體在凍脹力作用下將向上滑動，從而產(chǎn)生一定的凍脹位移量，使渠基土壤對襯砌體的有效凍脹力和有效凍結(jié)約束減小，襯砌體內(nèi)力值減小，抗裂強度提高，達到復合襯砌渠道防凍脹破壞的目的。當渠坡板與弧底板發(fā)生較大凍脹錯位時，渠坡板可能發(fā)生滑動失穩(wěn)。在實際工程中應(yīng)加強復合防滲襯砌渠道坡板和弧形底板之間的伸縮縫管理，或使用整體現(xiàn)澆結(jié)構(gòu)，以防止滑動失穩(wěn)。

c．弧底梯形復合襯砌渠道的弧形底板反拱作用明顯，與梯形斷面渠道相比，其變形恢復能力與抗凍脹破壞能力明顯增強。

［1］肖旻，賀興宏，李壽寧．梯形渠道混凝土襯砌凍脹破壞的力學模型探討［J］．東北水利水電，2010(10):58-61．(XIAO Min，HE Xinghong，LI Shouning．Discussion on mechanical model of frost swelling destroyfor concrete lining of trapezoid channel［J］．Water Resources ＆Hydropower of Northeast， 2010(10):58-61．(in Chinese))

［2］李安國．我國渠道防滲工程技術(shù)綜述［J］．防滲技術(shù)，2000，6(1):1-4，16．(LI Anguo．The reviewed of canal seepage control engineering technology in China［J］．Technique of Seepage Control，2000，6(1):1-4，16．(in Chinese))

［3］王正中．梯形渠道混凝土襯砌凍脹破壞的力學模型研究［J］．農(nóng)業(yè)工程學報，2004，20(3):24-29．(WANG Zhengzhong．Establishment and application of mechanics models of frost heavingdamage of concrete lining trapezoidal open canal［J］．Transactions of the CSAE，2004，20(3):24-29．(in Chinese))

［4］劉旭東．混凝土襯砌渠道抗凍脹技術(shù)措施及其機理研究［D］．楊凌:西北農(nóng)林科技大學，2011．

［5］艾樹衡，陳清芬，余祥全，等．混凝土及膜料復合渠道防滲技術(shù)的研究［J］．防滲技術(shù)，1996，2(3):9-12．(AI Shuheng，CHEN Qingfen，YU Xiangquan，et al．Study of technique for canal seepage prevention by composite of concrete and membrane［J］．Technique of Seepage Prevention，1996，2(3):9-12．(in Chinese))

［6］李甲林．渠道襯砌凍脹破壞力學模型及防凍脹結(jié)構(gòu)研究［D］．楊凌:西北農(nóng)林科技大學，2008．

［7］郭利霞．渠道凍脹力學模型及有限元分析［D］．楊凌:西北農(nóng)林科技大學，2007．

［8］李超．弧底梯形混凝土襯砌渠道凍脹力學模型及數(shù)值模擬［D］．石河子:石河子大學，2010．

［9］王正中，李甲林，陳濤，等．弧底梯形渠道混凝土襯砌凍脹破壞的力學模型研究［J］．農(nóng)業(yè)工程學報，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 CSAE，2008，24(1):18-23．(in Chinese))

［10］劉旭東，王正中，閆長城，等．基于數(shù)值模擬的雙層薄膜防滲襯砌渠道抗凍脹機理探討［J］．農(nóng)業(yè)工程學報，2011，27(1):29-35．(LIU Xudong，WANG Zhengzhong，YAN Changcheng，et al．Exploration on anti-frost heave mechanism of lining canal with double films based on computer simulation［J］．Transactions of the CSAE，2011，27(1):29-35．(in Chinese))

［11］姜海波，侍克斌．壩坡復合土工膜防滲體的抗滑穩(wěn)定分析［J］．水資源與水工程學報，2010，21(6):15-18．(JIANG Haibo，SHI Kebin．Stability analysis against slide of anti-seepage structure indam slope with composite geomembrane［J］．Journal of Water Resources＆ Water Engineering，2010，21(6):15-18．(in Chinese))

［12］姜海波，侍克斌，劉亮． 復合土工膜與粗粒料的摩擦特性試驗研究［J］． 中國農(nóng)村水利水電，2011( 3) : 86-89．( JIANG Haibo，SHI Kebin，LIU Liang． Experimental research on friction characteristics between compsite geomembrane and coarse-grained materials ［J］． China Rural Water and Hydropower，2011 ( 3 ) : 86-89． ( in Chinese) )

［13］鄭源，湯驊，姜海波．U形復合襯砌渠道凍脹破壞力學模型研究［J］．中國農(nóng)村水利水電，2013(11):80-84．(ZHENG Yuan，TANG Hua，JIANG Hai-bo．Research on mechanics models of forst heaving damage of composite linling U-shape canals［J］．China Rural Water and Hydropower，2013(11):80-84．(in Chinese))