易　富,朱鳳薇,楊宇婷

(遼寧工程技術(shù)大學土木與交通學院，阜新　123000)

?

溫度應(yīng)力耦合作用下瀝青路面應(yīng)力與變形分析

易富,朱鳳薇,楊宇婷

(遼寧工程技術(shù)大學土木與交通學院，阜新123000)

通過對瀝青路面內(nèi)部溫度場和應(yīng)力場耦合機理的分析，建立了熱應(yīng)力平衡方程，并且研究了溫度對材料力學性能的影響。應(yīng)用ADINA模擬軟件對溫度應(yīng)力耦合作用下的瀝青路面的應(yīng)力和變形變化規(guī)律進行有限元模擬。結(jié)果表明：在耦合作用下無論是位移還是應(yīng)力，在面層內(nèi)均不再是單調(diào)均勻變化，而是隨溫度場的變化出現(xiàn)了波動;并且把位移和應(yīng)力的變化規(guī)律進行分時段描述,得到了溫度與應(yīng)力耦合作用下瀝青路面面層應(yīng)力和變形波動變化的規(guī)律，同時位移和應(yīng)力的變化也是進行防車轍，防疲勞損傷考察的要素，對高等級瀝青路面防車轍防疲勞損傷具有一定的參考意義。

瀝青路面； 溫度應(yīng)力耦合； 應(yīng)力； 變形

1　引　言

完全處于自然環(huán)境中的瀝青路面常年受到車輛荷載、溫度應(yīng)力、濕度等因素的影響。炎熱的夏季，路面在車輛荷載的作用下會向四周推擠；而干冷的冬天，路面會發(fā)生凍脹和干裂[1]；北方地區(qū)的冷熱交替產(chǎn)生的溫度應(yīng)力周期性變化會使路面產(chǎn)生一定程度的損壞即疲勞破壞，會降低道路的使用壽命，對人民生命財產(chǎn)安全構(gòu)成重大威脅[2]。故為了減輕瀝青路面的疲勞損傷，延長道路的使用年限，深入的研究瀝青路面在溫度和應(yīng)力耦合作用[3]下的特殊規(guī)律是十分必要的。

國內(nèi)外對瀝青混合料進行了深入的研究并取得了大量的成果。馬正軍等[4]分析了熱力的參數(shù)對瀝青路面的影響，給出了在非標準熱力參數(shù)下瀝青路面溫度均值和標準差的修正公式；陳小燕等[5]用數(shù)值模擬的方法分析了夏日溫度全天變化下路面結(jié)構(gòu)層溫度場分布及不同降溫速率下瀝青路面溫度場分布；莊傳儀等[6]建立了基于晴天(或陰天)與雨天路表溫度的預估模型；Suh等[7]基于測試路面和加速路面車轍性能模型進行研究；Ozgan[8]確定了瀝青路面在不同溫度下的穩(wěn)定性進行了分析。以上研究成果多數(shù)依賴于數(shù)值模擬的方法，并且沒有考慮溫度場和應(yīng)力場耦合作用對瀝青路面的影響。實際上，溫度場的變化一方面以溫度應(yīng)力的形式與應(yīng)力場進行耦合[9]，另一方面影響了瀝青路面材料的力學性能。因此，有必要對溫度場應(yīng)力場耦合作用下瀝青路面應(yīng)力變形特性進行深入研究。

本文通過ADINA軟件建立有限元模型，分析了在溫度應(yīng)力耦合作用下瀝青路面的應(yīng)力和變形變化規(guī)律，可以為高等級瀝青路面進行車轍，疲勞損傷的防護提出一定的依據(jù)，延長了瀝青路面的使用年限。

2　溫度應(yīng)力耦合作用

由于瀝青路面常年暴露在復雜多變的環(huán)境中，周圍溫度的變化造成在路面內(nèi)部形成周期性變化的應(yīng)力。同時，具有粘彈性的瀝青材料的力學性質(zhì)也會對溫度的變化形成響應(yīng)[10]。

2.1溫度場

根據(jù)熱力學定律，熱分析的溫度場控制方程[11]為：

(1)

邊界條件為：

(2)

式中，qS=h(θe-θS)，θe為外界環(huán)境溫度，θS為土體表面溫度，h為熱對流系數(shù)；Kx、Ky、Kz為介質(zhì)三方向熱傳導率；qB為范圍內(nèi)熱源密度(單位體積熱生成率)；S1為上邊界溫度，S2為上邊界熱流量。

2.2溫度應(yīng)力耦合方程

(3)

對式(3)兩邊取跡有：

(4)

式中，εb=εxx+εyy+εzz為體積應(yīng)變；Pc=σkk/3=(σxx+σyy+σzz)/3為圍壓(平均正應(yīng)力)；K為宏觀體積模量；體積熱膨脹系數(shù)為3α。

由式(3)、(4)可以得到溫度應(yīng)力耦合作用下的應(yīng)力應(yīng)變方程：

σij=2Gεij+λ·tr(εij)δij+3αKδijθ

(5)

式中，λ為拉梅常數(shù)，λ=2Gν/(1-2ν)。

2.3變溫對力學參數(shù)的影響

瀝青路面在溫度場下流動變形是溫度和時間的函數(shù)。在試驗研究中發(fā)現(xiàn)，不同溫度和時間試驗測定得到的力學參數(shù)曲線具有大致相同的形狀。如果將溫度T1測得的曲線Er(T1,t)沿對數(shù)時間方向水平移動lgαT，則與溫度T0時測得的曲線Er(T0,t)大致相同，這種關(guān)系可以表示為：

Er(T1,t)=Er(T0,t/αT)

(6)

式中，αT為位移因子，僅是溫度的函數(shù)，αT=T1/T0。

由Williams、Landel和Ferry理論及試驗[13]分析，αT可以用式(7)求得：

(7)

式中，T0為參考溫度；取玻璃態(tài)脆化點溫度Tg

2.4外部荷載對瀝青路面的作用

荷載以一定速度在路面上行駛，忽略荷載受到路面凹凸的影響?？砂押奢d分為兩類：突加移動荷載和勻速移動荷載[15]。

突加移動荷載就是有一數(shù)值為的集中負荷突然作用在路面結(jié)構(gòu)上，然后又以恒定速度沿方向移動。其數(shù)學表達式[16]為：

F(t)=Pδ(x-νt)δ(y)H(t)

(8)

式中，δ(·)為狄拉克函數(shù)；H(t)為單位階躍函數(shù)。勻速移動荷載就是有一數(shù)值為的集中荷載在路面結(jié)構(gòu)上行駛，勻速移動速度為，其數(shù)學表達式為：

F(t)=P(t)δ(x-νt)δ(y)

(9)

2.5邊界條件

方程(6)所建立的溫度場與應(yīng)力場耦合模型中含有體積應(yīng)力和溫差兩個變量，若使該方程封閉，必須加上定解條件。上邊界條件為溫度按照正弦規(guī)律變化。外界氣溫無規(guī)律的波動對地溫影響的深度甚微，可不計。底基層底部溫度不受氣溫日變化影響。對于溫度場定解條件為：

(10)

對于應(yīng)力場，考慮如下邊界條件：

(11)

式中，▽為Laplus算子。

3　溫度應(yīng)力耦合有限元分析

3.1路面結(jié)構(gòu)參數(shù)

深度方向分別取18cm的瀝青混凝土面層，36cm的水泥穩(wěn)定碎石基層，20cm的低劑量水泥穩(wěn)定碎石底基層和5.26m的土基。對于各結(jié)構(gòu)層的材料參數(shù)，認為基層、地基層和土基為彈性體，瀝青混凝土面層為粘彈性體，即采用修正Burgers模型。取溫度為40 ℃，泊松比為0.4，瀝青路面材料性質(zhì)見表1。

對于瀝青混凝土，材料的泊松比隨溫度發(fā)生改變，其值介于0.15～0.45[17]，低溫時(小于0 ℃)用低值(0.15)，高溫時(50 ℃以上)取高值(0.45)。本文選擇瀝青混凝土40 ℃參數(shù)值，泊松比為0.4。

路面結(jié)構(gòu)體材料參數(shù)的如表2所示。

表1　面層瀝青混凝土的修正Burgers模型參數(shù)值(40 ℃)

表2　路面各結(jié)構(gòu)層材料參數(shù)

3.2有限元模型的建立

圖1　瀝青混凝土路面有限元模型Fig.1　Finite element model of Asphalt Concrete Pavement

行車荷載采用標準軸載0.7MPa，以5s為周期作用于路面上，作用時間為24h，分析瀝青混凝土路面的位移應(yīng)力變化規(guī)律。

對于邊界條件，結(jié)合實際情況，上邊界為自由邊界，兩側(cè)限制其側(cè)向位移而底邊固定。實際上，作用在行車道上的車輛荷載，路面兩側(cè)及底邊基本沒有響應(yīng)。瀝青路面的有限元模型如圖1。

3.3瀝青路面應(yīng)力與變形分析

將一天內(nèi)各個時刻對應(yīng)的豎向位移、水平應(yīng)力、剪應(yīng)力、豎向應(yīng)力分別分時段描述，其隨時間的變化關(guān)系如圖像2～5。

從圖2中可以看出，面層內(nèi)的豎向位移變化以非均勻的速率增加，而是以波動的形式增加，路表處尤其明顯。

圖2　路表豎向位移Fig.2　Vertical displacement of surface

圖3　底基層底面水平應(yīng)力Fig.3　Horizontal stress of sub-basing bottom

從圖3中可以看出，溫度場與應(yīng)力場耦合作用時，瀝青混凝土路面的最大拉應(yīng)力仍然出現(xiàn)在底基層底面。

圖4　面層地面剪應(yīng)力Fig.4　Shear stress of facing bottom

圖5　路表豎向應(yīng)力Fig.5　Vertical stress of surface

從圖4中可以看出，最大剪應(yīng)力出現(xiàn)在面層底面，最小剪應(yīng)力出現(xiàn)在底基層底面。面層內(nèi)剪應(yīng)力呈微小波動下降趨勢，路表的剪應(yīng)力波動最大但沒有出現(xiàn)下降現(xiàn)象，振幅小于0.02MPa。

從圖5中可以看出，最大的豎向應(yīng)力出現(xiàn)在路表，隨著深度的增加而逐漸減小，最小值位于底基層底面。面層內(nèi)豎向應(yīng)力呈波動遞減趨勢，路表處振幅最大，底基層底面的豎向應(yīng)力僅為路表的1/8。

為了更加深入的分析位移和應(yīng)力隨時間的變化規(guī)律，取波動比較明顯的路表處豎向位移和水平應(yīng)力分時段進行分析。路表處的位移和水平應(yīng)力變化規(guī)律可用表3～4來描述。

表3　路表豎向位移變化規(guī)律

表4　路表水平應(yīng)力變化規(guī)律

4　結(jié)　論

在溫度場和應(yīng)力場耦合作用下，無論是位移還是應(yīng)力均表現(xiàn)出了波動性，尤其是路表處表現(xiàn)的更為明顯。通過有限元分析了在溫度和應(yīng)力耦合作用下瀝青路面的響應(yīng)，分析結(jié)果如下：

(1)路表處的豎向位移增幅約為面層底面和底基層底面的2倍；底基層底面的剪應(yīng)力約為面層底面的1/8；底基層底面的豎向應(yīng)力約為路表的1/8；

(2)瀝青混凝土路面的面層對溫度敏感性較基層和底基層強烈。溫度波動的幅度自上而下呈衰減特性，因此溫度在道路的設(shè)計中是重點考慮的因素。溫度疲勞的最主要區(qū)域是在上面層，因其溫度梯度最大，所以在防治溫度疲勞裂縫和車轍等路面災(zāi)害的措施中要重點考察溫度疲勞；

(3)瀝青路面的主要承重部分為基 層及以下的結(jié)構(gòu)體，荷載大小對瀝青路面的變形與應(yīng)力的影響更為顯著，因此，在瀝青路面設(shè)計，使用時應(yīng)首先考慮荷載大小的因素。

[1] 楊建軍，鄭健龍，錢國平.路面凍脹開裂熱-應(yīng)力耦合分析[J].公路交通科技,2010,27(3):29-35.

[2] 不同類型瀝青混凝土凍融疲勞性能研究[J].硅酸鹽通報,2013,32(11):2385-2388.

[3] 易富，金艷，高健，等.耦合和非耦合效應(yīng)下瀝青路面受力對比分析研究[J].硅酸鹽通報2015,34(3):803-807.

[4] 馬正軍，談至明，錢晨，等.全國范圍瀝青混凝土路面路表溫度的分布特征研究[J].公路,2013,(3):1-5.

[5] 陳小燕，劉小文.瀝青混凝土路面溫度場數(shù)值分析[J].鐵道建筑,2011,(12):141-142.

[6] 莊傳儀，王林，申愛琴，等.瀝青路面路表溫度預估模型研究[J].公路交通科技,2010,27(3):39-44.

[7]SuhYC,ChoNH.Developmentofaruttingperformancemodelforasphaltconcretepavementbasedontestroadandacceleratedpavementtestdata[J].KSCE Journal of Civil Engineering, 2014, 18(1):165-171.

[8]ErcanO.Determiningthestabilityofasphaltconcreteatvaryingtemperaturesandexposuretimesusingdestructiveandnon-destructivemethods[J].Journal of Applied Sciences, 2007, 24(7):3870.

[9] 任德斌，蘇博.瀝青路面溫度應(yīng)力的有限元分析[J].沈陽建筑大學學報,2010,26(4):699-703.

[10] 張肖寧.瀝青與瀝青混合料的粘彈力學原理與應(yīng)用[M].人民交通出版社,2005.

[11] 易富，梁冰.瀝青混凝土路面溫度場和應(yīng)力場耦合模型[J].應(yīng)用基礎(chǔ)與工程科學學報,2009,4(8):597-604.

[12] 賀玉龍,楊立中,楊明.巖體溫度場與應(yīng)力場耦合作用的一種量化方法[J].西南交通大學學報，2002,37(1):10-13.

[13] 何曼君，張紅東，陳維孝，等.高分子物理[M].上海.復旦大學出版社，2007.

[14] 曹麗萍,譚憶秋,董澤蛟,等.應(yīng)用玻璃化轉(zhuǎn)變溫度評價SBS改性瀝青低溫性能[J].中國公路學報,2006,19(2):1-6.

[15] 陳俊,陳景雅,劉云,等.不同形式荷載下瀝青混凝土路面結(jié)構(gòu)力學響應(yīng)的對比分析[J].公路，2013,(8):138-134.

[16] 李澤光.狄拉克函數(shù)δ(f)的作用[J].大連大學學報,2012,33(3):14-18.

[17] 邵顯智，邵敏華，畢玉峰，等.瀝青混合料泊松比的測試方法[J].同濟大學學報(自然科學版),2006,34(11):1470-1474.

ActionofStressandDeformationAnalysisofAsphaltPavementunderTemperature-stressCoupling

YI Fu,ZHU Feng-wei,YANG Yu-ting

(InstituteofCivilEngineeringandTransportationLiaoning,TechnicalUniversity,Fuxin123000,China)

Throughtoanalysisoftheasphaltsurfaceinternalcouplingmechanismofthetemperaturefieldandstressfield,thethermalstressequilibriumequationisestablished,andstudiestheinfluenceofthetemperatureonmaterialmechanicsperformance.ApplingsoftwareADINAtoconductthefiniteelementsimulationofthetemperaturestressunderthecouplingeffectofthechangelawofstressanddeformationofasphaltsurfacepavement.Theresultsshowthatbothdisplacementandstressunderthecouplingeffect,inthesurfacelayerisnolongerdrabuniformvariation,butalongwiththechangeoftemperaturefieldofthevolatility.Andforthefirsttimeforperiodoftimetodescribethechangeruleofdisplacementandstressandtheruleofstressanddeformationofasphaltpavementsurfacefluctuationchangeunderthecoupledactionoftemperatureandstressisobtained,andatthesametime,thechangeofdisplacementandstressistopreventfurrow,antifatiguedamageoftheelements,andtohigh-gradeagainstruttingofasphaltpavementfatiguedamagehascertainreferencesignificance.

asphaltconcrete;temperature-stresscoupling;stress;deformation

遼寧省教育廳科研基金項目(L2013144);百千萬人才基金(2011921064)

易富(1978-)，男，教授，碩導.主要從事道路工程及環(huán)境巖土工程方面的研究.

U416

A

1001-1625(2016)01-0316-06