車(chē)忠新

(撫順市新賓滿(mǎn)族自治縣水務(wù)局，遼寧 新賓 113200)

科學(xué)研究及工程設(shè)計(jì)

小型U形玻璃鋼渠道抗凍脹數(shù)值模擬

車(chē)忠新

(撫順市新賓滿(mǎn)族自治縣水務(wù)局，遼寧 新賓 113200)

玻璃鋼渠道在我國(guó)北方部分灌溉渠道建設(shè)中得到應(yīng)用，其抗凍脹效果得到普遍認(rèn)可。由于地區(qū)差異大，其性能表現(xiàn)也可能存在差異。本文以遼寧新賓縣渠道工程為例，通過(guò)ANSYS有限元分析，對(duì)小型U形玻璃鋼渠道模型、混凝土渠道模型的抗凍脹性能進(jìn)行對(duì)比，結(jié)果表明：兩類(lèi)渠道模型保溫性能差異很??；與混凝土渠道相比，小型U形玻璃鋼渠道韌性更強(qiáng)、應(yīng)力分散作用更顯著，小型U形玻璃鋼渠道可以更好地維持渠道襯砌板穩(wěn)定性。

U形玻璃鋼渠道；抗凍脹性能；數(shù)值模擬

我國(guó)東北地區(qū)耕地多為粉質(zhì)黏土，當(dāng)氣溫降低至冰點(diǎn)后，若土壤含水量達(dá)到一定程度，則極易出現(xiàn)凍脹現(xiàn)象。目前這些地區(qū)的灌溉渠道多采用混凝土結(jié)構(gòu)，一旦出現(xiàn)基土凍脹，可能導(dǎo)致渠道變形、塌陷，從而影響渠道使用穩(wěn)定性，還可能引發(fā)嚴(yán)重事故，危害農(nóng)田及生產(chǎn)者安全。為有效應(yīng)對(duì)渠道基土凍脹危害，黑龍江、寧夏、內(nèi)蒙古和青海等地嘗試推廣玻璃鋼渠道[1-2]。事實(shí)上在水文和地質(zhì)條件一定時(shí)，溫度是影響基土凍脹的關(guān)鍵因素[3-4]。在我國(guó)各地的溫度差異十分顯著，在引入玻璃鋼渠道時(shí)，需要對(duì)溫度變化與抗凍脹性能的關(guān)系進(jìn)行研究。本文以新賓縣農(nóng)業(yè)灌溉渠道為例，通過(guò)建立小型U形玻璃鋼渠道模型，采用ANSYS軟件對(duì)其抗凍脹性能進(jìn)行模擬。通過(guò)模擬分析，希望對(duì)推動(dòng)小型U形玻璃鋼渠道的應(yīng)用有積極作用。

1 小型U形玻璃鋼渠道原型概況

小型U形玻璃鋼渠道原型，為遼寧省新賓縣農(nóng)業(yè)灌溉工程示范區(qū)。通過(guò)對(duì)該區(qū)域進(jìn)行地質(zhì)測(cè)量，渠道基土呈顯著粉質(zhì)特征，屬于凍脹性土質(zhì)。該渠道斷面為U形。根據(jù)工程設(shè)計(jì)，整體呈南北走向，總長(zhǎng)度為1.50km，縱比降為1/8000；渠道寬度為0.98m，深度為0.83m；渠底圓心角為2α=155°，對(duì)應(yīng)半徑R為0.45m；斜坡直線段長(zhǎng)度L為0.50m，其垂直傾斜度β=13°。在完成渠道開(kāi)挖后，有兩種材質(zhì)備選方案：?鋪設(shè)混凝土，厚度為0.06m；?鋪設(shè)玻璃鋼，厚度為0.005m(玻璃鋼試塊見(jiàn)圖1)。

2 小型U形玻璃鋼渠道有限元模型

2.1 基本假定

在對(duì)玻璃鋼渠道進(jìn)行建模分析時(shí)，不可能對(duì)濕度、溫度、含水量、基土土質(zhì)及襯砌體剛度等進(jìn)行全面、精確模擬，因此，需對(duì)關(guān)鍵要素作如下假定：?基土各個(gè)方向的性質(zhì)均勻分布；?在含水量、土質(zhì)等條件既定時(shí)，僅將溫度作為抗凍脹性能的影響因素；?凍結(jié)過(guò)程中，基土封閉，不再有水體滲入或者滲出；?襯砌體對(duì)渠道抗凍脹不產(chǎn)生影響，渠道各個(gè)方向的抗凍脹性能不存在差異。

2.2 有限元基本原理

2.2.1 熱傳導(dǎo)方程

遼寧地區(qū)的渠道凍結(jié)過(guò)程較長(zhǎng)，與熱傳導(dǎo)過(guò)程相同。因此，可根據(jù)傅里葉導(dǎo)熱定律建立如下穩(wěn)態(tài)二維熱傳導(dǎo)微分方程[5]：

(1)

式中T——溫度，℃；

λ——坐標(biāo)軸不同方向的導(dǎo)熱系數(shù)，W/(m·℃)；

ρ——渠道容重，kg/m3；

Q——熱源強(qiáng)度(內(nèi)部)，W/kg；

t——持續(xù)時(shí)間，h；

c——比熱，J/(kg·℃)。

通過(guò)式(1)，可以建立有限元分析所需的溫度場(chǎng)。然而，在研究渠道導(dǎo)熱問(wèn)題時(shí)，還需確定初始條件、邊界條件。其中，初始條件表示渠道開(kāi)始導(dǎo)熱時(shí)(t=0)的溫度分布特征，其表達(dá)式為

(2)

邊界條件是指渠道基土邊緣的換熱特征，其表達(dá)式為

(3)

在建立起穩(wěn)態(tài)溫度場(chǎng)后，即默認(rèn)渠道基土內(nèi)無(wú)熱源。此時(shí)，Q=0，且各測(cè)試點(diǎn)溫度恒定，即

(4)

于是，可將式(1)作如下簡(jiǎn)化：

(5)

因?yàn)橐呀?jīng)假設(shè)溫度在各個(gè)方向無(wú)差異，即

(6)

于是，式(5)可進(jìn)一步簡(jiǎn)化為

(7)

2.2.2 應(yīng)力應(yīng)變方程

當(dāng)環(huán)境溫度變化時(shí)，渠道基土、襯砌體等會(huì)發(fā)生凍脹變形，從而出現(xiàn)結(jié)構(gòu)應(yīng)變。依據(jù)熱應(yīng)力理論，此時(shí)的渠道結(jié)構(gòu)應(yīng)力應(yīng)變方程為

(8)

式中γxy——剪應(yīng)變；

ε——正應(yīng)變；

τxy——剪應(yīng)力；

E——彈性模量；

α——襯砌體自由膨脹系數(shù)；

Δt——溫度變化值；

μ——泊松比。

2.3 參數(shù)選取

在采集渠道工程數(shù)據(jù)后，根據(jù)式(1)～式(8)，可得渠道模型主要參數(shù)，結(jié)果見(jiàn)表1。同時(shí)，凍土材料彈性模量隨著溫度變化有所不同：當(dāng)溫度為-1℃時(shí)，彈性模量為10MPa；當(dāng)溫度為-2℃時(shí)，彈性模量為26MPa；當(dāng)溫度為-3℃時(shí)，彈性模量為33MPa；當(dāng)溫度為-5℃時(shí)，彈性模量為46MPa。

表1 玻璃鋼與混凝土材料參數(shù)選取

2.4 有限元模型建立

根據(jù)“中國(guó)季節(jié)性?xún)鐾翗?biāo)準(zhǔn)凍深線圖”，結(jié)合所選渠道工程實(shí)際情況，將凍深最大值確定為1.20m。由圖2可見(jiàn)，小型U形玻璃鋼渠道基土，最大凍深線沿著渠頂板、渠底板下方1.20m分布。以最大凍深線為界線，其上為凍土區(qū)(標(biāo)記為Ⅰ)，其下為非凍土區(qū)(標(biāo)記為Ⅱ)。

圖2 渠道基土結(jié)構(gòu)(單位：m)

2.5 有限元計(jì)算結(jié)果分析

在計(jì)算有限元模型熱邊界條件時(shí)，其上邊界條件取凍融期最低氣溫，其下邊界條件取凍融期最高氣溫。由表2可知，本研究中的有限元模型上邊界、下邊界溫度分別為-15.05℃、10.40℃。由于渠道模型位移邊界不受約束，因而屬于自由凍脹。

表2 2012—2016年凍融周期月平均氣溫 單位：℃

2.5.1 溫度場(chǎng)分析

溫度是影響渠道基土凍脹性能的最重要參數(shù)，在很大程度上決定著渠道的穩(wěn)定性。通過(guò)將邊界溫度施加于渠道模型中，可得其溫度場(chǎng)分布云圖(見(jiàn)圖3)。在凍土區(qū)，溫度變化十分明顯，溫差梯度也十分明顯；在非凍土區(qū)，溫度線則十分平整。可見(jiàn)，渠道基土深層受到邊界溫度的影響十分微弱。同時(shí)，通過(guò)比較渠頂部、斜坡以及底板處溫度變化情況可以發(fā)現(xiàn)，渠底板下凍土深度明顯更小。這表明底板可以起到隔絕冷空氣的作用，有利于保護(hù)渠底穩(wěn)定性。

圖3 渠道溫度場(chǎng)云圖(單位：℃)

2.5.2 位移場(chǎng)分析

渠道基土發(fā)生凍脹后，會(huì)產(chǎn)生垂直方向的凍脹力，同時(shí)，渠道基土與襯砌板之間會(huì)形成凍結(jié)力。在以上兩種外力作用下，渠道及其襯砌體會(huì)向上位移。為了更加直觀地觀察渠道位移結(jié)果，將其變位數(shù)值展開(kāi)(見(jiàn)圖4)。渠道整體存在上抬特征。

圖4 渠道襯砌板變位數(shù)值模擬展開(kāi)圖

由圖4可知，兩類(lèi)渠道凍脹位移呈現(xiàn)如下幾個(gè)特征：與混凝土渠道位移值相比，陽(yáng)坡板、陰坡板以及渠底板處，小型U形玻璃鋼位移值更小。不過(guò)，二者的位移值變化特征十分類(lèi)似，最大值均出現(xiàn)在兩側(cè)(陽(yáng)坡板略大于陰坡板)；最小值均出現(xiàn)在渠底板。最大位移值之所以會(huì)出現(xiàn)在頂板處，是因?yàn)轫敯寤敛皇芗s束，在自由凍脹狀態(tài)下的凍脹量也就較大。

在渠底，混凝土渠道位移值為4.14cm，小型U形玻璃鋼位移值為4.54cm，前者要小于后者；在兩側(cè)，混凝土渠道位移值為2.32cm，小型U形玻璃鋼位移值為3.26cm，前者要小于后者。

可見(jiàn)，與混凝土渠道相比，小型U形玻璃鋼位移值更大。然而，這并不意味著后者安全性和穩(wěn)定性更低。這是因?yàn)?，小型U形玻璃鋼的彈性模量、厚度均更小，基土凍脹分布更加均勻，對(duì)渠道的凍脹力作用也更均勻；同時(shí)，小型U形玻璃鋼的韌性更佳，當(dāng)凍土融化后，渠道能夠快速恢復(fù)至原狀。

2.5.3 應(yīng)力場(chǎng)分析

在進(jìn)行應(yīng)力場(chǎng)模擬時(shí)，須從法向凍脹力、切向凍脹力這兩方面進(jìn)行。

當(dāng)基土凍脹后，會(huì)對(duì)襯砌板形成法向作用力，即為法向凍脹力。當(dāng)溫度設(shè)定在-1℃后，基土開(kāi)始出現(xiàn)凍脹現(xiàn)象?；炷燎滥Ｐ团c小型U形玻璃鋼渠道模型的法向凍脹力最大值均出現(xiàn)在頂板處。這是因?yàn)椋敯逄幨艿剿胶涂v向雙重凍結(jié)，加之其約束呈剛性，因而應(yīng)力最為集中?；炷燎滥Ｐ偷姆ㄏ騼雒浟ψ畲笾禐?.32MPa，小型U形玻璃鋼渠道模型的法向凍脹力最大值為6.05MPa?？梢?jiàn)，采用玻璃鋼渠道可以將法向凍脹力最大值降低4.46%，這可以更有效地防止渠道局部應(yīng)力過(guò)于集中，避免影響其使用安全性。

當(dāng)基土出現(xiàn)凍脹后，為了阻止渠道襯砌板上抬，凍土?xí)纬上孪蜃饔昧?，即為切向凍脹力。在凍土區(qū)，兩類(lèi)模型的切向凍脹力最大值均出現(xiàn)在渠坡與底板相切處左側(cè)。這是因?yàn)?，渠坡板受到的約束多，且剛性特征顯著，故應(yīng)力在此最為集中。與最大值分布相對(duì)應(yīng)的是，切向凍脹力最小值出現(xiàn)在陽(yáng)坡一側(cè)。這是因?yàn)椋?yáng)坡溫度相對(duì)更高，基土凍脹量更小，加之底板約束小，因而形成較小的切向凍脹力。混凝土渠道模型所對(duì)應(yīng)的切向凍脹力分布在區(qū)間[4.69MPa，8.44MPa]，小型U形玻璃鋼渠道模型所對(duì)應(yīng)的切向凍脹力分布在區(qū)間[4.43MPa，8.02MPa]。

可見(jiàn)，小型U形玻璃鋼渠道可以有效地降低切向凍脹力，使基土凍脹所引起的渠道變形可以得到更好的釋放，從而增加渠道穩(wěn)定性。

3 結(jié) 語(yǔ)

本文以新賓縣灌溉工程示范區(qū)的渠道為原型，分別構(gòu)建混凝土渠道模型、小型U形玻璃鋼渠道模型；在完成有限元原理討論后，選定主要參數(shù)，并運(yùn)用ANYSYS軟件構(gòu)建有限元模型。經(jīng)過(guò)有限元計(jì)算，繪制出兩類(lèi)渠道模型的溫度場(chǎng)、位移場(chǎng)等云圖。

通過(guò)分析和討論，得出如下幾點(diǎn)結(jié)論：與混凝土渠道模型相比，小型U形玻璃鋼渠道模型在韌性、應(yīng)力分散以及位移控制等方面有更好表現(xiàn)，同時(shí)，小型U形玻璃鋼渠道具有更好的防滲效果，可以防止凍脹過(guò)程中補(bǔ)入含水量而加劇襯砌板變形；兩類(lèi)模型的保溫效果不存在顯著差異；小型U形玻璃鋼渠道可以有效地降低法向凍脹力、切向凍脹力，提升渠道抗壓/抗拉性能，從而有效釋放渠道變形，提升渠道運(yùn)行安全性和穩(wěn)定性。

[1] 鄭輝玲. 甘州區(qū)盈科干渠工程凍脹破壞的原因及防治措施[J].水利建設(shè)與管理, 2011(6):48-49.

[2] 柴振軍. 弧底梯形渠道無(wú)喉道量水槽水動(dòng)力學(xué)性能分析[J].中國(guó)水能及電氣化, 2016(10):37-40.

[3] 李根,何武全,宋清林,等. 混凝土襯砌渠道置換防凍脹標(biāo)準(zhǔn)化技術(shù)模式[J].節(jié)水灌溉,2015(4):79-83.

[4] 唐少容,王紅雨,李存云.考慮襯砌板跨度影響的U形渠道熱——力耦合研究[J].灌溉排水學(xué)報(bào),2015,(11):65-69.

[5] 程傳勝,田軍倉(cāng),王斌,等. 旱寒地區(qū)U形玻璃鋼和混凝土渠道抗凍脹性能的研究[J]. 水資源與水工程學(xué)報(bào),2015(6):173-177.

Numerical Simulation on Freezing Risistance in Small-sized U-shaped Glass Steel Channels

CHE Zhongxin

(FushunXinbinManchuAutonomousWaterSuppliesBureau,Xinbin113200,China)

Glass fiber reinforced plastics have been applied in building of irrigation channels in North China and its freezing risistance effects have obtained high approval. However, due to big regional variance, their performance might vary. This article puts channel project in Xinbin County, Liaoning as an example. Through ANSYS finite element analysis, it carries out comparison on small-scaled U-shaped glass steel channel model and concrete channel model. Results show that the two channel models show minor difference in property of thermal conservation. Compared with concrete channel, small-scaled U-shape glass steel channels are better in tenacity and more significant in stress dispersion effect, so stability of canal lining can be better maintained.

U-shape glass steel channel; freezing risistance performance; numerical simulation

10.16617/j.cnki.11-5543/TK.2017.04.013

TV672+.2

A

1673-8241(2017)04- 0049- 04