摘 要：【目的】通過建立渠道混凝土襯砌結(jié)構(gòu)數(shù)值計(jì)算模型，深入了解渠道混凝土襯砌結(jié)構(gòu)和渠基土壤在凍脹條件下的溫度變化和變形特征。【方法】利用有限元軟件建立包含渠道襯砌結(jié)構(gòu)和渠基土壤的數(shù)值模型，并考慮渠基土壤溫度影響，模擬計(jì)算渠道底部溫度和襯砌結(jié)構(gòu)的凍脹位移。【結(jié)果】研究結(jié)果表明，無論是襯砌結(jié)構(gòu)和渠基土壤，考慮接觸還是未考慮接觸，渠基土壤溫度不同深度溫度變化均相同，在考慮接觸的襯砌結(jié)構(gòu)中，陰坡和陽坡的凍脹變位較渠底更大?！窘Y(jié)論】襯砌結(jié)構(gòu)接觸或未接觸條件下，渠基土壤溫度分布規(guī)律是一致的，襯砌結(jié)構(gòu)的陰坡和陽坡在頂部凍脹位移和變形存在明顯差異。研究成果對渠道襯砌結(jié)構(gòu)的抗凍脹能力方面的研究有一定的參考價(jià)值。

關(guān)鍵詞：渠道襯砌；數(shù)值模擬；凍脹特性

中圖分類號：TV698.26" "文獻(xiàn)標(biāo)志碼：A" " "文章編號：1003-5168（2024）11-0053-04

DOI：10.19968/j.cnki.hnkj.1003-5168.2024.11.011

Simulation Analysis of Frost Heare Characteristics of

Channel Lining Structures

ZHANG Yongling XIAO Rang FU Junke

（College of Civil Engineering，Hexi University，Research Institute of Water Resources Protection and

Utilization in Hexi Corridor，Zhangye 734000，China）

Abstract： [Purposes] This study aims to establish a numerical calculation model for the concrete lining structure of the channel， in order to gain a deep understanding of the temperature changes and deformation characteristics of the concrete lining structure and the soil of the channel foundation under frost heave conditions. [Methods] A numerical model was established using finite element software， which includes the channel lining structure and the soil of the channel foundation. The influence of soil temperature on the channel foundation was considered， and the temperature at the bottom of the channel and the frost heave displacement of the lining structure were simulated and calculated. [Findings] The" results show that regardless of whether the lining structure and the canal foundation soil consider contact or not， the temperature change of the canal foundation soil is the same at different depths. In the lining structure considering contact， the frost heave displacement of the shaded and sunny slopes is greater than that of the canal bottom. [Conclusions] The distribution pattern of soil temperature in the canal foundation is consistent under the contact or non-contact conditions of the lining structure， and there are significant differences in frost heave displacement and deformation at the top between the shaded and sunny slopes of the lining structure，which has certain reference value for improving the frost heave resistance of channel lining structures.

Keywords： channel lining; numerical modeing; frost heave characteristics

0 引言

渠道襯砌結(jié)構(gòu)是水利工程中常用的一種重要結(jié)構(gòu)形式，主要用于保證渠道的堅(jiān)固性和穩(wěn)定性。合理的襯砌結(jié)構(gòu)設(shè)計(jì)方案可以提高渠道的輸水效率，并減少渠道效率受到工程環(huán)境和輸水工況的影響。目前，越來越多的渠道工程在寒冷地區(qū)進(jìn)行建設(shè)，而在寒冷地區(qū)，凍脹是一個(gè)常見問題，會對渠道襯砌結(jié)構(gòu)的性能造成一定影響［1-3］。

針對凍土渠基環(huán)境下的渠道建設(shè)，特別需要關(guān)注襯砌結(jié)構(gòu)的防凍脹設(shè)計(jì)［4］。渠基底部凍土的存在會導(dǎo)致渠道附近土壤在凍融循環(huán)過程中土壓力產(chǎn)生變化，進(jìn)而引起渠道襯砌結(jié)構(gòu)產(chǎn)生凍脹應(yīng)力。劉超等［5］采用試驗(yàn)與模擬相結(jié)合的方法，在冬季停水渠道的研究中，利用有限差分軟件建立了渠基模型。研究結(jié)果顯示，渠基溫度與深度呈負(fù)相關(guān)，凍脹量在時(shí)間推移中經(jīng)歷了快速凍結(jié)階段、緩慢凍結(jié)階段和穩(wěn)定階段的變化。李彬斌［6］通過使用ANSYS建立了北方寒區(qū)的膜袋混凝土襯砌渠道有限元模型，并對具體工程進(jìn)行結(jié)構(gòu)優(yōu)化計(jì)算。郭富強(qiáng)等［7］建立了有限元渠道模型，并對鋪設(shè)不同厚度聚苯板的渠道進(jìn)行了溫度場和凍脹的模擬，并提出了預(yù)制混凝土渠道適宜的保溫模式。

因此，通過研究襯砌結(jié)構(gòu)的凍脹模型可以更好地了解襯砌厚度、土壤溫度等因素對渠道防凍脹的影響。研究成果可為輸水渠道襯砌結(jié)構(gòu)的設(shè)計(jì)和施工提供一定的參考。

1 研究方法

本研究采用有限元法，將待解區(qū)域逐個(gè)分割成有限數(shù)量的單元，并對其進(jìn)行離散化處理。這些單元的集合組成了一個(gè)離散的網(wǎng)格，在每個(gè)網(wǎng)格點(diǎn)上進(jìn)行線性插值，建立一個(gè)插值函數(shù)，用于描述待解問題的連續(xù)體［8］。本研究選擇了一個(gè)具體的渠道斷面作為研究對象，該斷面為U形50 mm厚預(yù)制混凝土板襯砌渠道，下部采用一定深度區(qū)域作為渠基墊層，幾何模型及網(wǎng)格如圖1所示，黑色區(qū)域?yàn)榛炷烈r砌結(jié)構(gòu)，淺黑色區(qū)域?yàn)榍?。?jīng)過有限元法的網(wǎng)格處理，該斷面被劃分為多個(gè)剖分的單元。在每個(gè)剖分單元中，利用離散網(wǎng)格點(diǎn)上的函數(shù)值進(jìn)行插值計(jì)算，從而建立了一個(gè)線性插值函數(shù)，為后續(xù)模擬計(jì)算做準(zhǔn)備。

在本次模擬中，使用特定的材料參數(shù)來描述凍土的性質(zhì)和行為。凍土的導(dǎo)熱性是決定溫度分布的關(guān)鍵參數(shù)之一。本次襯砌結(jié)構(gòu)和土壤等模擬采用的材料參數(shù)見表1。

在進(jìn)行溫度計(jì)算時(shí)，為了準(zhǔn)確模擬物體的溫度分布，應(yīng)在模型的邊界上施加溫度邊界條件。通過施加溫度邊界條件，可以確定模型中的未知溫度值，并計(jì)算整個(gè)溫度分布。在應(yīng)力應(yīng)變計(jì)算中，為了計(jì)算物體的應(yīng)力和應(yīng)變分布，可施加適當(dāng)?shù)奈灰七吔鐥l件。位移邊界條件可以包括固定的位移值或受到外部力作用產(chǎn)生的位移等。

2 渠基土壤溫度模擬分析

在進(jìn)行渠道襯砌結(jié)構(gòu)凍脹研究時(shí)，首先進(jìn)行渠基土壤溫度分析?？紤]接觸和未考慮接觸溫度分布如圖2所示。由圖2可知，無論考慮襯砌結(jié)構(gòu)接觸土壤與否，渠基土壤的溫度分布規(guī)律都是一致的。研究結(jié)果顯示，土壤溫度在不同深度范圍內(nèi)有所變化。渠底40 cm深度以下的土壤溫度變化快，而30～40 cm深度的變化相對緩慢。陰坡上部20～30 cm深度的變化速度略低于30～40 cm深度。陰坡下部的10～20 cm深度的變化梯度稍小于20～30 cm深度。陽坡下部30～40 cm深度的梯度稍小于40～50 cm深度。此外，在陽坡上部20 cm深度處，土壤溫度變化速度較快。

無論是否考慮襯砌結(jié)構(gòu)與渠基土壤接觸，渠基底部溫度分布具有高度相似性。在沿著渠基法線方向的不同深度范圍內(nèi)，土壤溫度變化呈現(xiàn)明顯的梯度差異。尤其在底部40 cm深度和陽坡上部20" cm深度的土壤溫度變化較快。相比之下，渠底30～40 cm深度的土壤溫度梯度較小。

渠基最大凍深溫度和最大融深溫度如圖3所示。在沿著襯砌結(jié)構(gòu)中間法線方向觀察時(shí)，渠基土壤的凍結(jié)深度每增加10 cm，對應(yīng)的溫度大約增加0.7～0.8℃。渠基土壤溫度計(jì)算對比結(jié)果如圖4所示。凍結(jié)深度和溫度兩個(gè)變量之間呈現(xiàn)近似線性遞增。

3 考慮接觸的襯砌結(jié)構(gòu)凍脹位移分析

襯砌結(jié)構(gòu)考慮接觸渠基土壤和未考慮接觸渠基土壤情況下的溫度分布基本相同。因此，在應(yīng)力應(yīng)變分析中，對考慮接觸的襯砌結(jié)構(gòu)凍脹位移進(jìn)行模擬計(jì)算后無須額外計(jì)算未考慮接觸情況下的位移變化。根據(jù)圖5所示的模擬結(jié)果顯示，凍脹變形后的砌板整體產(chǎn)生向上的位移，兩側(cè)受到了土壤擠壓，并出現(xiàn)微小偏轉(zhuǎn)變形。

考慮襯砌結(jié)構(gòu)接觸土壤情況下的計(jì)算結(jié)果顯示，襯砌結(jié)構(gòu)在陰坡和陽坡的凍脹位移分別為7.49 mm和7.20 mm，而渠底的凍脹位移為3.24 mm。這表明陰坡和陽坡的凍脹變位較大，而渠基的凍脹變位相對較小。

凍脹力分布和主應(yīng)力分布如圖6、圖7所示。由圖5和6可知，在分析接觸襯砌結(jié)構(gòu)時(shí)，發(fā)現(xiàn)陰坡和陽坡的凍脹力和位移分布存在顯著差異。陰坡的上部水平凍脹位移為3.50～5.53 mm，而中下部為1.47～3.50 mm。相比之下，陽坡的上部水平凍脹位移為2.60～4.64 mm，而中下部為0.57～2.60 mm。陰坡頂部的水平凍脹位移最大值為6.55 mm，而陽坡頂部為5.66 mm。

陰坡的中上部豎直凍脹位移為5.06～6.42 mm，而下部為4.56～4.99 mm。相比之下，陽坡的中上部豎直凍脹位移為4.38～5.74 mm，而下部為4.04～4.38 mm。此外，可以得到陰坡頂部的豎直凍脹位移最大值為7.11 mm，而陽坡頂部的最大值為6.77 mm。

綜上所述，在考慮接觸的襯砌結(jié)構(gòu)中，陰坡和陽坡的凍脹位移分布具有一定的差異。陰坡和陽坡上部的橫向變位較中下部更大，且陰坡頂部的豎向變位也較陽坡頂部大。

4 結(jié)論

本研究針對渠道襯砌結(jié)構(gòu)的凍脹特性進(jìn)行了詳細(xì)的模擬分析，建立了渠道混凝土襯砌結(jié)構(gòu)的數(shù)值計(jì)算模型，并分析了渠基土壤溫度的變化，主要得出以下結(jié)論。

①考慮襯砌結(jié)構(gòu)接觸土壤和未考慮襯砌結(jié)構(gòu)接觸土壤的渠基溫度分布規(guī)律一致。在不同深度范圍內(nèi)，渠基和陰坡、陽坡的土壤溫度變化有所不同。

②考慮襯砌結(jié)構(gòu)接觸土壤情況下，襯砌結(jié)構(gòu)的陰坡和陽坡的凍脹位移較大，渠基的凍脹位移較小，并且頂部位移存在明顯差異。其中考慮襯砌結(jié)構(gòu)接觸土壤后，襯砌結(jié)構(gòu)的陰坡和陽坡凍脹位移計(jì)算值為7～7.5 mm，渠基的凍脹位移計(jì)算值為3 mm左右。

參考文獻(xiàn)：

［1］張宇峰，王琦，蘇娟麗，等.幾種典型渠道斷面的實(shí)用經(jīng)濟(jì)性對比分析［J］.水利與建筑工程學(xué)報(bào)，2022，20（6）：46-49，142.

［2］張晨，王羿，韓孝峰，等.考慮接觸損傷效應(yīng)的襯砌渠道凍脹過程數(shù)值模擬方法［J］.巖土工程學(xué)報(bào)，2022，44（S2）：188-193.

［3］趙公亮.混凝土襯砌渠道的管理與維修養(yǎng)護(hù)淺析［J］.黑龍江水利科技，2022，50（6）：162-163，201.

［4］亞力坤.輸水渠道襯砌結(jié)構(gòu)凍脹特性影響變化仿真研究［J］.四川水利，2023，44（3）：17-21.

［5］劉超，郭妍秀.季凍區(qū)渠基凍脹模擬分析［J］.水利科學(xué)與寒區(qū)工程，2023，6（7）：1-5.

［6］李彬斌.寒區(qū)膜袋混凝土襯砌渠道結(jié)構(gòu)優(yōu)化［J］.水利科學(xué)與寒區(qū)工程，2023，6（9）：4-7.

［7］郭富強(qiáng)，李鋼鐵，霍軼珍，等.基于ADINA的河套灌區(qū)預(yù)制混凝土渠道保溫防凍脹效果數(shù)值模擬［J］.干旱區(qū)資源與環(huán)境，2023，37（10）：84-92.

［8］王海麗.混凝土襯砌渠道凍脹破壞有限元模擬分析［J］.水利科技與經(jīng)濟(jì)，2015，21（8）：27-29.

收稿日期：2023-11-15

基金項(xiàng)目：甘肅省高等學(xué)校創(chuàng)新基金項(xiàng)目“季節(jié)性凍土區(qū)農(nóng)田冬灌回歸水入滲對墊層混凝土渠道凍脹的誘發(fā)機(jī)理研究”（2021A-118）；甘肅省教育廳高等學(xué)校創(chuàng)新基金項(xiàng)目“不同斷面結(jié)構(gòu)明渠信息化精確量水關(guān)鍵技術(shù)研究”（2020B-212）；甘肅省級大學(xué)生創(chuàng)新訓(xùn)練計(jì)劃項(xiàng)目（S202110740046）；橫向課題“山丹縣馬營河大型灌區(qū)、寺溝河中型灌區(qū)農(nóng)田灌溉水有效利用系數(shù)測算分析”（H2023005）。

作者簡介：張永玲（1977—），女，碩士，副教授，研究方向：水利工程。

通信作者：肖讓（1979—），男，碩士，教授，研究方向：灌溉排水理論與技術(shù)。