程 鵬,張淑菡,陳艷金,趙潔冰
(張家口翰得交通公路勘察設計有限責任公司,河北 張家口 075000)
目前,國內橋面鋪裝設計壽命一般為15~20年,但在服務期限內都會或多或少地出現裂縫、波浪推移、局部擁包(瀝青包、高溫氣包)及車轍等一種或多種破壞現象。其中車轍破壞最為常見,同時是使道路工作者擔心的路面病害之一。這是由于車轍嚴重降低了路面的行車安全性,從而影響了高等級道路的投資效益。對車轍的控制和防治是目前迫切需要解決的問題,其中容許車轍深度是最為關鍵的指標之一。
本文采用理論分析與試驗研究相結合的方法,利用ANSYS有限元軟件對通泰大橋車轍進行分析和計算。
1.1.1 靜態蠕變試驗
為了預估瀝青混凝土路面的車轍,需采用恰當的實驗方法來確定本構方程中的材料參數,才能較為準確地模擬實際瀝青混凝土路面的應力和變形狀態。因此需要通過有效的實驗手段來對建立的本構關系相對應的材料參數進行確定。通過長期分析發現,蠕變試驗能夠較好的模擬瀝青混合料的粘彈性應變—時間關系,也能很好的對Burgers模型的材料參數進行求解。
采用澳大利亞UTM-100動態伺服液壓瀝青混合料試驗系統進行試驗,試驗中電腦自動采集數據,將蠕變試驗結果繪于以應變為縱軸,時間為橫軸的坐標系內,通過蠕變曲線能很直觀地了解瀝青混合料的粘彈性變形特性,分析瀝青混合料的變形隨時間的變化關系。蠕變試驗曲線如圖1所示。

圖1 SMA13的蠕變曲線
1.1.2 基于粒子群算法的Burgers模型參數擬合

圖2 粒子群算法計算流程
粒子群算法的基本思想是用于優化問題的每一個解稱為一個粒子,定義一個符合度函數來衡量每個粒子解的優越程度。每個粒子根據自己和其他粒子的“飛行經驗”游動,從而達到從全空間搜索最優解的目的。參數擬合流程如圖2所示。由靜態蠕變試驗所得到蠕變曲線關系圖以及大量數據,根據Burgers模型加載、卸載過程中的蠕變關系方程,利用粒子群算法原理進行編程分析計算,輸入蠕變實驗數據(應變ε及時間t)并根據已有文獻參考初值即可得到材料的Burgers的模型參數(E1、η1、E2、η2)。 參數擬合結果如表1所示。


表1 材料的Burgers模型參數
通泰大橋設計行車速度為60km/h,根據徐世法等人對輪載作用瀝青路面上的“代表時間”的研究[1],對該行駛速度范圍內的車輛,統一將車輛作用時間簡化定為0.015s以方便計算。確定后,每一點的累計加載時間就決定于車道輪載的橫向分布系數。車道內的橫向分布系數參考文獻[2],對于雙向六車道道路的主車道,車道內輪跡的橫向分布系數為0.33~0.36,本文按最嚴重的情況取推薦范圍上限0.36。由此便可以得到輪載對路面的累加作用時間為:

式中:N——累計當量軸次。
因此可計算出不同累計當量軸次下輪載對路面的累加作用時間(如表2所示)。
將由路面不平整度產生的附加動荷載進行簡化,使其形式更加簡單以便于進行瀝青路面車轍的計算。
現行路面設計方法采用豎向靜載作用下彈性多層體系的理論模型。對于實際道路上行駛的車輛,由于路面不平整而產生振動,路面受到的車輛荷載隨時間和路表特性改變而變化,是一個典型的動力荷載。根據汽車動力學的研究成果,汽車振動荷載的幅值一般為靜載的0.05~0.6倍,這主要與車輛懸架特性、車速和路面不平整度等有關。為了使計算簡便并使其結果具有代表性,也通過上述計算分析驗證,取汽車振動荷載幅值為靜載的0.3倍,并根據上述附加動荷載的建立得出了如下汽車動載模型方程:

式中:F(t)——汽車荷載;
P0——車輪靜載;
ω——振動頻率;
Pf′(t)——汽車振動荷載幅值,且

具體取值為:根據BZZ-100的標準,取單邊單軸靜載輪重,當車速為80km/h時,輪載作用時間取0.015s,由此可算ω=,T=0.015s。
以上取值經計算可得:

由于當量圓半徑δ=0.1065m,式(5)所得汽車荷載換算為接地壓強可得:

盡管該汽車動載模型形式簡單,但對汽車特性、幾何線形、路況及車速等多方面因素進行了綜合考慮,計算簡便,可以用來模擬汽車荷載。汽車荷載與時間的關系如圖3所示。

圖3 汽車荷載與時間的關系
通過ANSYS有限元分析軟件中APDL(ANSYS Parametric Design Language參數化設計語言[3])對汽車動載模型函數進行輸入即可得到與圖3汽車荷載與時間關系圖一致的荷載典線。該車輛荷載模型與已有的模型相比,具有形式簡單、便于計算的優點,且與幾何曲線、路況及車速關系密切,可以用來模擬汽車荷載。對瀝青路面進行加載,以便于進行車轍的分析。
車轍是瀝青路面高溫季節的產物,它對溫度有極高的敏感性,經研究表明,在溫度低于20℃時,瀝青混合料產生的永久變形可以忽略,所以不會產生車轍,因此,在車轍計算中可以只考慮路面溫度高于20℃的月份。
雖然路面內溫度場的一維熱傳導方程可以應用在不同邊界條件和方法來求解該偏微分方程,得到溫度場的解析式或直接算得不同時刻在不同深度處的溫度值。但在實際的應用過程中,并不需要確切知道路面內任意處的溫度。因此AI提出了確定路面代表溫度的經驗式,根據得出的溫度值,來選擇瀝青以及其他材料。
AI認為瀝青層的典型溫度由下式決定:

式中:T——路面月平均代表溫度,℉;
T0——大氣月平均溫度,℉;
Z——瀝青面層厚度的1/3厚度,英寸。
為了計算通泰大橋瀝青路面的代表溫度,統計了張家口地區1971~2000年的平均氣溫情況(如表3所示),并用AI方法計算了張家口地區不同月份瀝青路面的代表溫度。

表3 張家口平均氣溫統計(1971~2000年)
從表3可知,張家口地區4~9月路面代表溫度高于20℃,代表溫度區間23.3~31.9℃,取其平均溫度作本文計算時的代表溫度為27.9℃。
有了代表溫度便可以通過插值得到在該代表溫度下的前面通過蠕變試驗得到的瀝青混合料Burgers模型的參數,由此將溫度問題轉換為模型參數中在代表溫度下的模量及Burgers模型參數來進行計算。
利用ANSYS有限元計算程序,建立橋面二維計算模型。考慮到結構及荷載作用的整體性,可建立瀝青路面如結構B的實體模型進行分析。
模型建立后,計算單元瀝青層采用PLANE182,聚丙烯纖維混凝土及正交異性鋼橋面板則采用PLANE42。瀝青面層為粘彈性材料,材料的本構模型采用Burgers模型在ANSYS內使用其的Prony級數形式,因此要進行輸入必須對其進行相應的轉換。聚丙烯纖維混凝土及正交異性鋼橋面板則認為是完全彈性體系,材料本構模型采用廣義虎克定律。選用路面結構的各層材料參數如表4所示,路面有限元模型如圖4所示。

表4 通泰大橋計算模型材料參數

圖4 ANSYS中橋面模型
采用加載次數10萬次~2000萬次對通泰大橋車轍深度即永久變形進行了計算,計算結果如圖5所示。

圖5 不同加載次數時的車轍
圖5所示的是軸載作用次數與路面車轍之間的關系。隨著軸載作用次數的增加瀝青路面的車轍深度逐漸增大。在軸載作用次數的前期(<100萬次),路面車轍深度增長相對緩慢,當軸載作用次數達到100萬次以后車轍增加速度有所增加。總體而言,此瀝青橋面鋪裝抗車轍性能較好,當軸載次數達到2000萬時,其車轍深度僅為3.978mm。
為比較不同鋼橋面板厚度對瀝青鋪裝車轍的影響,計算了12mm、16mm、20mm三種鋼板厚度情況下的瀝青鋪裝車轍。計算結果如圖6所示。

圖6 不同鋼橋面板厚度時的車轍
從以上圖表看出,橋面瀝青鋪裝車轍深度隨著鋼橋面板厚度增加表現出先增大后略有減小的趨勢。但是,從增長的幅度來看,當鋼橋面板厚度達到16mm,車轍深度最大為3.978mm,當鋼橋面厚度為12mm時,車轍最大深度為3.727mm,前者相對于后者僅增加了6.7%。同時,3.978mm的車轍深度仍然處于很低的水平,遠小于文獻[4]中車轍的養護標準規定的15mm車轍深度。這表明,總體而言通泰大橋瀝青鋪裝抗車轍能力強,而鋼橋面板厚度對車轍深度的影響較小。
為比較不同瀝青鋪裝厚度對瀝青鋪裝車轍的影響,計算了50mm、70mm、90mm三種瀝青鋪裝厚度情況下的車轍。計算結果如圖7所示。

圖7 不同瀝青鋪裝厚度時的車轍
從以上圖表看出,當瀝青鋪裝厚度為70mm時,車轍深度達到5.580mm,相對于50mm瀝青鋪裝結構增加了1.602mm,增加幅度達到40.3%。當瀝青鋪裝厚度為90mm時,車轍深度達到6.992mm,相對于70mm瀝青鋪裝結構增加了1.412mm,增加幅度為25.3%。這表明通泰大橋瀝青鋪裝車轍深度隨著瀝青鋪裝厚度增加而快速增加,但增加的趨勢逐步放緩。從減小車轍的方面考慮,通泰大橋采用50mm單層鋪裝,其抗車轍效果好于采用雙層瀝青鋪裝的更厚的瀝青鋪裝結構。
對影響瀝青路面車轍的主要因素進行了分析,討論了計算車轍時各種因素和車轍深度之間的關系,并對通泰大橋的車轍進行了計算,得到了以下幾個方面的結論:
a)針對通泰大橋進行的單軸靜載蠕變試驗,通過粒子群算法擬合得到了不同溫度下材料的Burgers模型參數(E1、η1、E2、η2);
b)將影響車轍變形深度的道路交通條件用交通量大小以及荷載重復作用次數來控制,通過汽車作用時間的研究得到了車轍計算時的交通條件參數——累加作用時間;
c)根據當地實地溫度資料統計出各月的平均氣溫,并以此為依據計算了當地的代表氣溫和路面代表溫度;
d)通過對通泰大橋的車轍進行計算,結果表明當軸載作用次數達到100萬次以后車轍增加速度有所增加;鋼橋面板厚度對瀝青鋪裝車轍深度影響較小,瀝青鋪裝厚度對車轍深度影響大。
[1]徐世法,朱照宏.按粘彈性理論預估瀝青路面車轍[J].同濟大學學報,1990,(9):299-305.
[2]湯文,孫立軍.車道荷載橫向分布影響下的路面車轍變形分析[J].上海公路,2011,(4):14-17.
[3]陳志華,劉紅波,周婷,等.空間鋼結構APDL參數化計算與分析[M].北京:中國水利水電出版社,2009.
[4]JTJ 073.2—2001,公路瀝青路面養護技術規范[S].