摘 要：文章結(jié)合哈爾濱到大連客運(yùn)專(zhuān)線路基防凍脹施工技術(shù)進(jìn)行分析通用有限元軟件ANSYS模擬路基溫度場(chǎng)的變化規(guī)律，將溫度場(chǎng)與變形場(chǎng)進(jìn)行耦合，分析熱作用對(duì)路基變形的影響，預(yù)測(cè)多年后路基中溫度場(chǎng)的變化情況及路基的變形趨勢(shì)，了解凍脹可能引起的危害，及時(shí)提出預(yù)防措施。

關(guān)鍵詞：路基凍脹；變形；預(yù)防措施

1 引言

在季凍區(qū)，由于溫度、水分的存在，使得溫度場(chǎng)發(fā)生改變，導(dǎo)致土體中發(fā)生冰-水相變，這對(duì)土體的力學(xué)性質(zhì)影響較大。通用有限元軟件ANSYS模擬路基溫度場(chǎng)的變化規(guī)律，將溫度場(chǎng)與變形場(chǎng)進(jìn)行耦合，分析熱作用對(duì)路基變形的影響，預(yù)測(cè)多年后路基中溫度場(chǎng)的變化情況及路基的變形趨勢(shì)，了解凍脹可能引起的危害，及時(shí)提出預(yù)防措施。

2 數(shù)學(xué)模型的建立

一般來(lái)說(shuō)凍脹發(fā)生會(huì)伴隨著路基內(nèi)熱（溫度），流（水）的遷移與重分布，因此路基內(nèi)水-熱作用的機(jī)理問(wèn)題可歸結(jié)為伴有相變的、考慮水分遷移的傳熱傳質(zhì)問(wèn)題。

2.1 傳熱傳質(zhì)數(shù)學(xué)模型

取連續(xù)介質(zhì)微元體，表面積 ，體積V，建立能量守恒方程：

（1）

其中，qi是熱流矢量在i方向的分量，T為溫度，t為時(shí)間，？籽為密度（即單位體積的質(zhì)量），Q為單位體積的熱生成率，c是比熱。由傅立葉定律，熱流矢量和溫度梯度存在如下關(guān)系：

（2）

其中，？姿ij是材料的導(dǎo)熱系數(shù)，對(duì)于各向同性材料，導(dǎo)熱系數(shù)為常數(shù)。

由式1和式2整理得微元體積V內(nèi)的導(dǎo)熱微分方程：

（3）

式中Q為溫度變化過(guò)程中的熱生成（即熱源），是溫度的函數(shù)，其它熱物理參數(shù)也是溫度的函數(shù)，如對(duì)于凍土而言，其導(dǎo)熱系數(shù)在正負(fù)溫度條件下差別極大，可用下式表示：

？姿=f（T） （4）

2.2 邊界條件

2.2.1 三類(lèi)邊界條件

在物體外表面？祝上通常有三類(lèi)不同的邊界條件：（1）第一類(lèi)邊界條件，表面？祝T上的溫度T是給定的，溫度隨空間位置和時(shí)間的不同而變化，滿(mǎn)足函數(shù)T=T（xi，t）。（2）第二類(lèi)邊界條件，在？祝q邊界上的熱流強(qiáng)度q是給定，熱流隨空間位置和時(shí)間變化滿(mǎn)足函數(shù)■（xi，t）=？姿n■，其中■表示？祝q的外法線方向。

（3）第三類(lèi)邊界條件，在？祝邊界上給定邊界面周?chē)黧w溫度和邊界面與流體之間的對(duì)流換熱系數(shù)，表達(dá)式為：q=h（T？祝-T∞）。

其中，T？祝是表面溫度，h是表面對(duì)流換熱系數(shù)，T∞是外界環(huán)境介質(zhì)溫度。對(duì)流換熱系數(shù)h通常是流動(dòng)的速度、流體換熱表面的結(jié)構(gòu)尺寸、環(huán)境介質(zhì)的溫度、材料熱性能系數(shù)、表面特征等的復(fù)雜函數(shù)。

2.2.2 附面層原理

路基表面由于收到太陽(yáng)輻射、表面水的蒸發(fā)和凝結(jié)、地面的紊流作用等，路基的換熱過(guò)程極其復(fù)雜，因此路基的表面溫度不能作為路基熱傳導(dǎo)計(jì)算的上邊界溫度，所以需要引入“附面層原理”。所謂“附面層”是指附著于地基上的具有一定厚度的物理層，這一物理層的溫度、濕度與上附空氣層和下附土層的溫度、濕度在大氣及太陽(yáng)輻射作用下具有不連續(xù)的突變現(xiàn)象，其中與空氣接觸的稱(chēng)為上附面層，與土壤相接的稱(chēng)為下附面層。“附面層原理”以變化穩(wěn)定的附面層底作為上邊界條件，摒棄了干擾因素眾多的土表面層，提高了路基傳熱傳質(zhì)計(jì)算的可靠性與精度。

2.2.3 上邊界條件

根據(jù)試驗(yàn)段實(shí)測(cè)的溫度，取路基的上邊界條件為表面下0.4m深度處的地溫函數(shù)。結(jié)合青藏高原多年凍土的大量研究，一般將地溫模擬為正玄函數(shù)。同時(shí)基于全球氣候變暖的趨勢(shì)，在未來(lái)20年，氣溫在逐漸升高，這一現(xiàn)象也在青藏高原凍土地區(qū)得到驗(yàn)證，隨著氣候變暖，凍土區(qū)氣溫升高，出現(xiàn)了凍土融化現(xiàn)象，因此，在東北季凍區(qū)也考慮了氣候變暖的影響，假設(shè)試驗(yàn)段地溫滿(mǎn)足標(biāo)準(zhǔn)正玄函數(shù)，形式為：

（5）

式中，T為溫度；t為時(shí)間，天；a、b、c、d和e為待定參數(shù)。

根據(jù)現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)，采用最小二乘法對(duì)實(shí)測(cè)表面下0.4m處的地溫曲線進(jìn)行擬合分析，求出待定參數(shù)。通過(guò)擬合，上邊界條件為以下形式：

（6）

3 求解方法

利用通用有限元軟件ANSYS進(jìn)行求解。由于路基縱向很長(zhǎng)，屬于平面應(yīng)變問(wèn)題，因此采用平面四邊形8節(jié)點(diǎn)單元進(jìn)行網(wǎng)格劃分。定義不同土體單元的材料特性，計(jì)算區(qū)域的邊界上施加邊界條件，邊界條件只考慮溫度條件和熱流密度條件，熱流率條件、對(duì)流條件和熱生成率條件不予考慮。計(jì)算時(shí)，首先讓原地表邊界的溫度條件在逐年不升溫的情況下計(jì)算2年，得到地基內(nèi)土體各點(diǎn)的初始溫度，然后刪除地表邊界的溫度約束，在填筑路基后的新邊界上施加逐年升高的溫度邊界，然后再計(jì)算20年。

4 模擬結(jié)果及分析

4.1 地溫線分析

不同年份相同日期地溫曲線見(jiàn)圖1。

（a）2013年1月26日 （b）2030年1月26日

（c）2013年4月16日 （d）2030年4月16日

圖1 各階段溫度場(chǎng)分布（單位：℃）

圖1（a）是2013年1月底的溫度分布圖，從現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)可知1月底是地溫負(fù)溫最大的階段，可以看出路基表面溫度低，越往下越高，熱量在向上傳遞，負(fù)溫在向下傳遞，說(shuō)明此時(shí)盡管表面溫度最低，但并不是凍深最大的時(shí)刻。由1（a）和（b）的對(duì)比可以看出，隨著氣候變暖和溫度周期循環(huán)作用，從2011年到2030年地溫各個(gè)深度的地溫都有所增加，溫度場(chǎng)穩(wěn)定，各個(gè)深度的溫度場(chǎng)形式基本保持不變。

圖1（c）是2013年4月中旬的溫度分布圖，可以看出盡管表面溫度已經(jīng)大于0℃，但內(nèi)部的地溫仍然有負(fù)溫存在，此階段負(fù)溫值很小，正在向正溫轉(zhuǎn)變，此時(shí)熱量向下傳遞，說(shuō)明這一階段凍深基本達(dá)到最大。此時(shí)的凍深為1.8m，與實(shí)測(cè)的1.66m基本接近。圖1（d）為2030年4月中旬的溫度場(chǎng)，可以看出，隨著溫度循環(huán)交替以及氣候變暖的作用，該階段的溫度場(chǎng)已經(jīng)完全處于正溫狀態(tài)，說(shuō)明路基凍結(jié)的時(shí)間縮短，最大凍深減小。

圖1（a）、（b）、（c）都有路基填料部分凍結(jié)深度小于原地基處凍結(jié)深度的規(guī)律，說(shuō)明填料有良好的抗凍性，這也與實(shí)測(cè)的結(jié)果一致。

4.2 路基變形數(shù)值模擬及分析

哈大客專(zhuān)對(duì)沉降要求嚴(yán)格，具體要求見(jiàn)表1，而季凍區(qū)土的凍脹融沉變化必然會(huì)引起路基的沉降變形，因此，控制路基沉降是哈大客專(zhuān)建設(shè)的關(guān)鍵問(wèn)題之一。因此，采用數(shù)值模擬對(duì)路基的沉降進(jìn)行預(yù)測(cè)是必要的，可以通過(guò)分析定性地了解路基沉降的發(fā)展趨勢(shì)，便于提前采取有效措施，防患于未然。

表1無(wú)碴軌道路基工后沉降控制值

參考文獻(xiàn)

[1]孟凡松，劉建平，劉永智.黑北公路凍土路基設(shè)計(jì)原則及病害特征[J].冰川凍土，2001，23（3）：307-311.

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

[3]TZ212-2005客運(yùn)專(zhuān)線鐵路路基工程施工技術(shù)指南[Z].

[4]TB10035-2002，J158-2002.鐵路特殊路基設(shè)計(jì)規(guī)范[S].

作者簡(jiǎn)介：李寶昌（1981-），男，黑龍江哈爾濱人，碩士，從事土木工程。