李南生,蘆 赫,王雨舟
(同濟(jì)大學(xué) 土木工程學(xué)院,上海 200092)
我國超過一半的面積處于凍土區(qū)域,其中多年凍土和季節(jié)性凍土占國土面積70%以上,是世界上第三大凍土國家,尤其是我國北部和青藏高原地區(qū)分布大量季節(jié)性凍土[1].隨著“一帶一路”戰(zhàn)略的進(jìn)一步實施,基礎(chǔ)建設(shè)向中亞地區(qū)發(fā)展,沿線區(qū)域工程建設(shè)將遇到前所未有的發(fā)展良機(jī),工程中遇到的凍融災(zāi)害問題應(yīng)得到更多的關(guān)注.
當(dāng)溫度低于0 ℃時,土體中的液相水轉(zhuǎn)化為固相冰,氣溫回暖后融化,多年循環(huán).在這一過程中由于水分凍融循環(huán)過程會引起土體的體積變化,宏觀表現(xiàn)為凍脹和融沉,對寒區(qū)工程土體有嚴(yán)重危害[2].從內(nèi)在機(jī)理上來看,土體的凍融循環(huán)過程是溫度場、水分場及力學(xué)場相互作用的復(fù)雜滲流及熱力學(xué)綜合作用過程[3].
針對這一問題,國內(nèi)外經(jīng)過大量理論及試驗研究提出了凍脹及融沉理論,并建立了基于溫度場、水分場和應(yīng)力場的耦合方程[4].Harlan[5]認(rèn)為在凍結(jié)過程中,由于熱水耦合作用引起水分的遷移,因此建立了水動力模型理論;Konard和Morgenstern[6]引入分凝勢概念,提出了水分遷移量和溫度梯度的關(guān)系.本文基于耦合理論,通過時域差分法建立碾壓土本構(gòu)方程,將冰水相變方程引入耦合方程,分析預(yù)測多年凍融循環(huán)作用下耦合場變化規(guī)律,以此得到水分遷移、應(yīng)力變化以及位移變化情況,并能預(yù)測多年凍融循環(huán)下路基沉降以及凍結(jié)深度.
寒區(qū)道路工程常采用碾壓土體作為基層材料,土體經(jīng)碾壓后可以有效降低土體孔隙率,提高土體密實度,進(jìn)而提高土體強(qiáng)度.因此作為基層材料,為了提高其初始強(qiáng)度,一般都要經(jīng)過碾壓處理.
在寒區(qū)的特殊環(huán)境下,一方面土體中冰水相變伴隨著熱量的轉(zhuǎn)移,另一方面含水率的變化會使得土水勢變化造成水分遷移,這些溫度場和水分場的變化引起重力重分布,往往是導(dǎo)致地質(zhì)災(zāi)害的主要原因.
針對工程土體的環(huán)境特點,在凍融循環(huán)過程中會引起熱、水、力三場耦合作用,因此,本文以連續(xù)介質(zhì)力學(xué)和熱力學(xué)為理論基礎(chǔ),分別建立了三場控制方程.
凍土地區(qū)不同于一般地區(qū),環(huán)境溫度會影響土中冰水相變,因此在建立數(shù)學(xué)模型時,基于土體非穩(wěn)態(tài)溫度場引入了相變方程[7].對于空間中每一點對應(yīng)一個溫度標(biāo)量T(x,y,z),凍土相變溫度場的溫度梯度為
(1)
假定土質(zhì)分層均勻,質(zhì)量遷移僅通過液相遷移和冰水相變,采用顯熱容法對相變區(qū)處理建立相變溫度場數(shù)學(xué)模型,即
(2)
C=Csθs+Cwθw+Ciθi+Caθa
(3)
λ=λsλwλiλa
(4)
式中:T為瞬態(tài)溫度;t為過程進(jìn)行的時間;λ為導(dǎo)熱系數(shù);ρw為水的容重;V為流體的相對速度矢量;C為熱容量;Cs、Cw、Ci、Ca分別為土顆粒熱容量、水熱容量、冰熱容量及空氣熱容量;θs、θw、θi、θa分別為體積土顆粒含量、體積含水量、體積含冰量及體積空氣含量;λs、λw、λi、λa分別為土的導(dǎo)熱系數(shù)、水的導(dǎo)熱系數(shù)、冰的導(dǎo)熱系數(shù)及空氣的導(dǎo)熱系數(shù).
由表6可以看出,赴菲游客在住宿方式上更偏好入住星級酒店,對經(jīng)濟(jì)型酒店和度假村風(fēng)險感知無較大差異。在餐飲類型上,赴菲游客更偏好選擇高檔酒店或者主題餐廳,游客更看中餐廳就餐環(huán)境、服務(wù)以及食品安全。在交通工具的選擇上,赴菲游客認(rèn)為乘坐旅行社統(tǒng)一安排的旅行中巴或者商務(wù)小巴更加安全可靠,而乘坐Tri-bike、螃蟹船和三輪摩托具有一定風(fēng)險。跳島游時需要乘坐螃蟹船,涉及水上活動感知風(fēng)險增加。島上自由活動時可能會選用Tri-bike或三輪摩托,而摩托的危險性更大,因而感知風(fēng)險越大。
(5)
式中:φ為總水頭;ψ為基質(zhì)熱或毛細(xì)壓力;K為有效導(dǎo)水系數(shù);ΔS為單位時間單位體積冰的變化.
一維穩(wěn)定或非穩(wěn)定對流傳熱的熱傳導(dǎo)方程為
(6)
(7)
結(jié)合總水頭計算公式和導(dǎo)熱系數(shù)計算公式得到Harlan水熱耦合模型[9],即
(8)
σ=Dε
(9)
式中:σ為應(yīng)力矩陣;D為剛度矩陣;ε為應(yīng)變矩陣.
考慮相變的應(yīng)力場控制方程為
σ=σg+DεV
(10)
式中:σg為由碾壓土體自身重力造成的引力;εV為土體凍結(jié)膨脹變形.
(11)
選取青藏高原某段道路路基作為研究對象,取路基中間豎直方向土體為研究對象,通過耦合模型構(gòu)建數(shù)值計算模型.
研究路段路基填土為粉質(zhì)黏土及弱風(fēng)化巖,其參數(shù)如表1、2所示.
表1 冰、水熱力學(xué)參數(shù)Tab.1 Thermomechanics parameters for ice and water
天然地表的溫度計算表達(dá)式為
(12)
圖1為第1,2,5,10年冬夏溫度場分布.根據(jù)土層溫度分布可知,地表溫度隨氣溫迅速降低,由于熱傳導(dǎo)效應(yīng)使得淺層土溫度劇烈變化,達(dá)到一定深度后變化幅度明顯減小,原因在于深層土熱量補(bǔ)給充分.可以發(fā)現(xiàn)在約0 ℃區(qū)域相變劇烈,由于耦合作用使得溫度傳導(dǎo)過程明顯減慢,體現(xiàn)于溫度場斜率達(dá)到0 ℃后明顯降低并趨于0.對比第2年、第5年、第10年結(jié)果可知,1.0 m深度土層處于相變劇烈區(qū),經(jīng)過多年凍融循環(huán)后趨于一個穩(wěn)定狀態(tài),尤其在第10年出現(xiàn)近1.0 m恒溫區(qū),阻礙溫度向更深處土層傳導(dǎo).
表2 土熱力學(xué)參數(shù)Tab.2 Thermomechanics parameters for soil
注:ρd為土體干密度;+、-分別為對應(yīng)土層的上下表層。
圖1 第1,2,5,10年冬夏溫度場分布Fig.1 Distribution of temperature field in summer and winter in year of 1,2,5 and 10
圖2為1 m深度處溫度變化.取深度1 m土溫度變化為研究對象,可知在冬季(凍融循環(huán)過程)出現(xiàn)溫度變化放緩區(qū),隨著凍融循環(huán)過程不斷重復(fù),這一區(qū)域的溫度場趨于一個穩(wěn)定值,溫度變化不再明顯.
圖2 1 m深度處溫度變化Fig.2 Temperature variation at depth of 1 m
圖3為第2,5,10年冬夏含水量分布.根據(jù)水分場分布可知,土層中水分在一定深度處急劇變化,這一深度與溫度場所反應(yīng)出的相變劇烈區(qū)相對應(yīng),隨著凍融循環(huán)過程的往復(fù),這一聚變區(qū)不斷加深并擴(kuò)張[10].表層土得到補(bǔ)給,含水量變化明顯,深層土與溫度場相似,達(dá)到一定深度區(qū)域穩(wěn)定值后,并不會因表層或相變引起水分波動而變化.
圖3 第2,5,10年冬夏含水量分布Fig.3 Distribution of water content in summer and winter in year of 2,5 and 10
圖4、5為1、2 m深度處含水量變化.對比1、2 m深度土層水分場可知,1 m深度水分場呈現(xiàn)明顯的周期性特點,當(dāng)超過2 m后周期性特點明顯降低,耦合作用效果降低,水分場波動多來源于水分補(bǔ)給或土層因素.
圖4 1 m深度處含水量變化Fig.4 Variation of water content at depth of 1 m
圖5 2 m深度處含水量變化Fig.5 Variation of water content at depth of 2 m
圖6為第2,5,10年冬夏應(yīng)力增量場.對比應(yīng)力增量可知,土層凍結(jié)區(qū)域出現(xiàn)一個峰值,由于土層中多場耦合及溫度變化引起相變和水分場變化,通過凍結(jié)力與應(yīng)力場產(chǎn)生耦合作用,造成了凍結(jié)區(qū)的應(yīng)力峰值.隨著多年凍融循環(huán),應(yīng)力增量僅表現(xiàn)在凍結(jié)區(qū),對深處土層影響減弱,與溫度場和水分場變化協(xié)調(diào),符合工程實際.
通過上述計算可以得到土層的位移和凍結(jié)深度變化,如圖7、8所示.隨著凍結(jié)次數(shù)增加,凍結(jié)深度有所增加,在達(dá)到第六次凍結(jié)后凍結(jié)深度不再發(fā)生變化,即凍結(jié)循環(huán)過程趨于穩(wěn)定.在本次計算中,最大凍結(jié)深度可達(dá)到1.8 m.
由圖7可知,正值表示地基表面升高,負(fù)值表示低于初始值,初始狀態(tài)為第一年冬天.隨著溫度降低,土體凍結(jié)導(dǎo)致體積膨脹,地基表層出現(xiàn)隆起.隨著固相冰向液相轉(zhuǎn)化,土體表面下沉,即融沉效應(yīng),同時淺層土體融化的液相水會導(dǎo)致局部水分增加,減弱土體有效應(yīng)力,加劇融沉,因此,土體下沉常表現(xiàn)為短時間急劇下降.隨著時間增加,凍融循環(huán)作用往復(fù),形成較為穩(wěn)定的凍結(jié)帶,融沉深度趨于穩(wěn)定.由于凍結(jié)深度緩慢加深,凍脹作用隨之不斷增加,長期來看對于道路工程是一種危害.
圖6 第2,5,10年冬夏應(yīng)力增量場Fig.6 Stress increment field in summer and winter in year of 2,5 and 10
圖7 路基表面位移變化Fig.7 Variation of surface displacement of highway foundation
圖8 凍結(jié)深度變化Fig.8 Variation of frozen depth
本文基于Duncan-Chang模型,考慮溫度場、水分場以及應(yīng)力場耦合條件建立了寒區(qū)凍土的控制方程.通過時域差分法,選取研究實例,對多年往復(fù)凍融循環(huán)進(jìn)行計算模擬,得到以下結(jié)論:
2) 本文在一維條件下將碾壓土體分層,通過時域差分法建立耦合計算模型,得到三場分布及變化規(guī)律,模型條件和結(jié)果符合實際情況.該計算模型能夠在一定精度下簡化計算過程,具有較高的計算效率.
3) 計算模型精確預(yù)測了土層劇烈應(yīng)變區(qū)域以及凍結(jié)深度,與實際變化規(guī)律一致,可用于對土層分析預(yù)測.本文結(jié)果反映了溫度場、水分場和應(yīng)力場的緊密耦合關(guān)系,尤其是在明顯的相變劇烈區(qū)域.