基于黏彈模型的交叉口路面力學(xué)響應(yīng)分析

孫思林 叢 林

(同濟(jì)大學(xué)道路與交通工程教育部重點(diǎn)實(shí)驗(yàn)室 上海 201804)

在城市道路交叉口信號(hào)控制區(qū)域，小轎車會(huì)頻繁的剎車和啟動(dòng)，產(chǎn)生較大的水平荷載，大大增加路面的拉應(yīng)力和剪應(yīng)力水平，使路表產(chǎn)生剪切推移變形，導(dǎo)致路面容易產(chǎn)生推移和車轍變形[1]。如何從車轍行成機(jī)理的角度，提出車轍的防治措施成為學(xué)者關(guān)注的問題。黃佳等[2]從荷載工況、溫度、車速、水平力作用等方面對(duì)路面受力響應(yīng)進(jìn)行分析。結(jié)果表明瀝青面層內(nèi)部剪應(yīng)力和豎向壓應(yīng)力是導(dǎo)致車轍病害的主要原因。

交叉口方面的已有研究多集中在交叉口材料的改性上[3]，朱志強(qiáng)等[4]開發(fā)了交叉口用高性能抗車轍混合料，并對(duì)混合料的性能進(jìn)行了驗(yàn)證。Liao[5]研發(fā)了一種適合交叉口路段使用的nail-reinforce瀝青混合料，用來抵抗剪切流動(dòng)變形、增強(qiáng)抗車轍效果；另一些學(xué)者通過有限元的角度，對(duì)路面的荷載特征進(jìn)行了分析，張洪偉等[6]建立了線彈性有限元模型，探究了不同行駛狀態(tài)時(shí)，路面最大剪應(yīng)力的分布，梁悅等[7]利用ANSYS有限元軟件建立了線彈性模型，探究了重載對(duì)交叉口路面剪應(yīng)力、剪應(yīng)變的影響。綜上，已有研究，未能充分考慮交叉口處水平制動(dòng)力較大的荷載特征，且采用的力學(xué)分析模型為線彈性模型，而瀝青混合料在溫度較高時(shí)主要表現(xiàn)為黏彈性。故本研究擬利用室內(nèi)動(dòng)態(tài)模量試驗(yàn)對(duì)黏彈有限元模型的參數(shù)進(jìn)行確定，探究水平力、重載、溫度對(duì)路面剪應(yīng)力的影響，從車轍產(chǎn)生機(jī)理的角度，為交叉口路面車轍的防治提供參考。

1 黏彈性參數(shù)的確定

1.1 動(dòng)態(tài)模量試驗(yàn)

依據(jù)ASHTTO TP 62-03規(guī)范的要求，在試驗(yàn)溫度為-10，4.4，21.1，37.8，54.4 ℃，正弦荷載頻率為25，10，5，1，0.5，0.1 Hz條件下，測(cè)試了SMA-13、AC-16、AC-20的動(dòng)態(tài)模量和相位角數(shù)據(jù)。對(duì)頻率采用對(duì)數(shù)坐標(biāo)，可得動(dòng)態(tài)模量和頻率的半對(duì)數(shù)坐標(biāo)圖，以SMA-13瀝青混合料為例，對(duì)瀝青混合料的黏彈性參數(shù)進(jìn)行介紹，圖1和圖2分別為SMA-13瀝青混合料的動(dòng)態(tài)模量和相位角的半對(duì)數(shù)坐標(biāo)圖。

圖1 SMA-13 動(dòng)態(tài)模量半對(duì)數(shù)坐標(biāo)圖

圖2 SMA-13 相位角半對(duì)數(shù)坐標(biāo)

由圖1可見，SMA-13動(dòng)態(tài)模量隨著頻率的增加而增大，不同溫度下瀝青混合料動(dòng)態(tài)模量隨頻率的變化趨勢(shì)相同，低溫時(shí)瀝青混合料的動(dòng)態(tài)模量大于高溫時(shí)瀝青混合料的動(dòng)態(tài)模量。由圖2可見，SMA-13瀝青混合料的相位角隨著頻率的增加而降低，隨著溫度的升高而升高。瀝青混合料動(dòng)態(tài)模量和相位角隨溫度頻率變化的特點(diǎn)充分反映了其黏彈性。

1.2 松弛模量主曲線

利用動(dòng)態(tài)模量的試驗(yàn)數(shù)據(jù)，結(jié)合相關(guān)公式對(duì)松弛模量進(jìn)行計(jì)算。

通過動(dòng)態(tài)模量和相位角計(jì)算E′(f)。

E′(f)=|E*|cosφ

(1)

式中：E′(f)為儲(chǔ)能模量，MPa；E*為動(dòng)態(tài)模量，MPa；φ為相位角，(°)。

對(duì)E′(f)進(jìn)行正規(guī)化，計(jì)算f和E′(f)的對(duì)數(shù)，得到lgf和lgE′(f)關(guān)系，在每個(gè)頻率點(diǎn)處，lgf對(duì) lgE′(f)進(jìn)行微分求解，計(jì)算出n。

n=[d lgE′(f)]/[d lg(f)]

(2)

按式(3)計(jì)算調(diào)整函數(shù)，先計(jì)算伽馬函數(shù)，再帶入公式求解調(diào)整函數(shù)λ′。

λ′=Γ(1-n)·cos (nπ/2)

(3)

進(jìn)而計(jì)算松弛模量E(t)

E(t)=E′(f)/λ′

(4)

對(duì)于黏彈性材料，作用時(shí)間t(或頻率f)和溫度θ可以相互轉(zhuǎn)化，低溫可以轉(zhuǎn)換為高頻。高溫可以轉(zhuǎn)換為低頻。將不同溫度下曲線移動(dòng)到參考溫度下即形成主曲線。表征各個(gè)溫度相對(duì)參考溫度的平移距離的參數(shù)稱為移位因子aθ。

aθ=t/ξ

(5)

式中：ξ為減縮時(shí)間。

利用西歌摩德(Sigmoidal)s型函數(shù)模型可對(duì)松弛模量的主曲線進(jìn)行擬合。通過回歸確定出不同溫度θ下的移位因子aθ。以21.1 ℃為參考溫度，將其他溫度下松弛模量經(jīng)平移，可得到21.1 ℃時(shí)，SMA-13松弛模量的主曲線，結(jié)果如圖3。

圖3 SMA-13松弛模量主曲線

(6)

式中：δ,α,β,γ為回歸參數(shù)。

相對(duì)于參考溫度的其他溫度的移位因子的計(jì)算，可以利用半理論半經(jīng)驗(yàn)公式Wilams-Landel-Ferry(WLF)實(shí)現(xiàn)，公式如式(7)。

(7)

式中：C1、C2為回歸系數(shù)，可由溫度θ下的移位因子aθ回歸確定；θref為參考溫度。

1.3 Prony級(jí)數(shù)

在ABAQUS有限元軟件中， Prony級(jí)數(shù)表示黏彈性材料的時(shí)間依賴性，輸入的參數(shù)包括一組Prony級(jí)數(shù)參數(shù)及與之對(duì)應(yīng)的延遲時(shí)間。剪切模量G(t)和剪切模量比g(t)的計(jì)算方法如式(8)、(9)所示。

(8)

式中：G(t)為剪切模量，MPa；ν為泊松比。

(9)

式中：N為Prony級(jí)數(shù)項(xiàng)數(shù)；gi為Prony級(jí)數(shù)參數(shù)；τi為延遲時(shí)間。

結(jié)合計(jì)算的剪切模量數(shù)據(jù)，采用5項(xiàng)Prony級(jí)數(shù)進(jìn)行參數(shù)擬合，擬合21.1 ℃的SMA-13、SMA-16、SMA-20的Prony級(jí)數(shù)數(shù)據(jù)見表1～表3。

表1 SMA-13級(jí)數(shù)

表2 AC-16級(jí)數(shù)

表3 AC-20級(jí)數(shù)

利用表1擬合的SMA-13混合料的5項(xiàng)Prony級(jí)數(shù)，對(duì)剪切模量比進(jìn)行預(yù)測(cè)，擬合及實(shí)測(cè)的剪切模量比見圖4。由圖4可見，預(yù)測(cè)數(shù)據(jù)和實(shí)際剪切模量比數(shù)據(jù)較為吻合，獲得的Prony級(jí)數(shù)較為準(zhǔn)確。

圖4 實(shí)測(cè)和Prony預(yù)測(cè)SMA-13瀝青混合料剪切模量比

2 有限元模型

參照J(rèn)TG D50-2017《公路瀝青路面設(shè)計(jì)規(guī)范》[8]規(guī)定的軸載設(shè)計(jì)參數(shù)，采用標(biāo)準(zhǔn)軸載為100 kN的單軸雙輪軸載，單輪接地當(dāng)量圓直徑為0.213 m,兩輪中心距為0.319 5 m, 輪胎接地壓強(qiáng)為0.7 MPa。

考慮路面荷載的作用形式為矩形加上2個(gè)半圓，進(jìn)一步將標(biāo)準(zhǔn)軸載簡化為2個(gè)0.192 m×0.184 m的矩形荷載，荷載間距為0.135 m,荷載大小為0.7 MPa。對(duì)于交叉口路面，考慮豎直荷載的同時(shí)，也需要考慮車輛制動(dòng)產(chǎn)生的水平荷載，下文中水平荷載用H表示，豎向荷載用V表示。

本文采用三維模型進(jìn)行分析，綜合考慮模型的計(jì)算精度和效率，模型尺寸選擇為6 m×3 m×3 m。其有限元模型見圖5。

圖5 路面有限元模型圖

路面結(jié)構(gòu)層為多層，層與層之間的連接是連續(xù)的，面層采用3層結(jié)構(gòu)。路面的結(jié)構(gòu)圖見圖6。

圖6 路面結(jié)構(gòu)圖(單位:cm)

結(jié)構(gòu)層材料參數(shù)見表4。

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

基于該模型對(duì)不同水平荷載、豎直荷載、溫度作用下，路面沿深度方向的最大剪應(yīng)力進(jìn)行探究，其中深度方向指圖5中通過K點(diǎn)的深度方向。

3 最大剪應(yīng)力影響因素分析

3.1 水平荷載

圖7為不同水平荷載作用下路面深度方向最大剪應(yīng)力曲線。

圖7 不同水平荷載作用下最大剪應(yīng)力深度方向分布

由圖7可見，在路面深度方向，當(dāng)施加0.7 MPa的豎向荷載時(shí)，最大剪應(yīng)力沿深度方向呈先增后減的趨勢(shì)，路面的最大剪應(yīng)力在路表以下8 cm附近達(dá)到最大；當(dāng)施加豎向荷載和水平力作用時(shí)，路面的最大剪應(yīng)力沿深度方向呈遞減趨勢(shì)，最大剪應(yīng)力在路表取得最大值,隨著水平荷載的增加，路面的最大剪應(yīng)力也在不斷增大。

有、無水平力作用下最大剪應(yīng)力在路表的分布情況見圖8、圖9。

圖8 有水平力時(shí)路表最大剪應(yīng)力

圖9 無水平力時(shí)路表最大剪應(yīng)力

由圖8、圖9可見，在無水平力作用時(shí)，最大剪應(yīng)力在路表呈現(xiàn)對(duì)稱分布，而在有水平荷載作用時(shí)，最大剪應(yīng)力在路表縱向呈非對(duì)稱分布，最大剪應(yīng)力的位置偏離荷載中心，在有水平荷載作用的交叉口路面，路面最大剪應(yīng)力的最不利位置位于圖5的G點(diǎn)，在瀝青混合料設(shè)計(jì)的過程中，可增加該點(diǎn)為剪應(yīng)力最不利位置驗(yàn)算點(diǎn)。

3.2 豎直荷載

豎向荷載對(duì)路面最大剪應(yīng)力的影響見圖10。

圖10 不同豎向荷載作用下最大剪應(yīng)力深度方向分布曲線

由圖10可見，隨著豎向荷載的增加，路面的最大剪應(yīng)力隨之增大，路面的最大剪應(yīng)力沿深度方向的變化趨勢(shì)相同，當(dāng)豎向荷載為0.84 MPa時(shí)，即超載為20%時(shí)，路面的最大剪應(yīng)力增加約15%，當(dāng)超載為50%時(shí)，最大剪應(yīng)力增加約43%。故對(duì)交叉口路段，應(yīng)對(duì)車輛的超載進(jìn)行嚴(yán)格限制以減小路面的最大剪應(yīng)力，達(dá)到防治車轍的效果。

3.3 溫度

溫度對(duì)最大剪應(yīng)力的影響見圖11。

圖11 不同溫度最大剪應(yīng)力深度方向分布曲線

由圖11可見，在21.1和54.4 ℃下，路面最大剪應(yīng)力沿深度方向的分布趨勢(shì)基本相同，最大剪應(yīng)力集中分布在路面的上面層和中面層；路面的上中面層對(duì)溫度更為敏感，隨著溫度的升高，上面層和中面層的最大剪應(yīng)力逐漸增大，54.4 ℃時(shí)的最大剪應(yīng)力為21.1 ℃時(shí)的1.3倍，而中面層以下部分最大剪應(yīng)力基本保持不變，其主要因?yàn)樯稀⒅忻鎸訛轲椥圆牧希瑢?duì)溫度變化較為敏感。故夏季高溫時(shí)段可采取灑水降溫等方式，降低路面溫度，減小最大剪應(yīng)力，從而防治車轍。

4 結(jié)語

在無水平荷載作用時(shí)，路面最大剪應(yīng)力的最大值位于路面的中面層，當(dāng)有水平荷載作用時(shí)，路面最大剪應(yīng)力的最大值位于荷載圓沿道路縱向的最前方G點(diǎn)，在城市路面設(shè)計(jì)階段，對(duì)瀝青面層進(jìn)行剪應(yīng)力最大值驗(yàn)算時(shí)，可考慮增加該點(diǎn)為剪應(yīng)力控制點(diǎn)。

升高溫度會(huì)使路面的最大剪應(yīng)力增大，且上、中面層的最大剪應(yīng)力對(duì)溫度更加敏感。對(duì)于高溫重載的交叉口路段，在進(jìn)行上、中面層材料設(shè)計(jì)時(shí)，可采用添加抗車轍劑和使用改性材料等措施提高其抗車轍性能，同時(shí)，夏季高溫時(shí)段可采取灑水降溫等方式，降低路面溫度，防治車轍。

豎向荷載大小對(duì)路面的最大剪應(yīng)力有著顯著的影響，可通過加大對(duì)超載車輛的查處力度,嚴(yán)格限制超載等措施,實(shí)現(xiàn)交叉口車轍的防治。