李銳鐸， 樂金朝， 李瑞霞， 姚立陽

(1.鄭州大學(xué) 水利與環(huán)境學(xué)院，鄭州　450002； 2.河南城建學(xué)院 交通運輸工程學(xué)院，河南 平頂山　467036)

?

基于分?jǐn)?shù)階導(dǎo)數(shù)的瀝青膠漿動態(tài)力學(xué)特性

李銳鐸1,2， 樂金朝1， 李瑞霞1， 姚立陽2

(1.鄭州大學(xué) 水利與環(huán)境學(xué)院，鄭州450002； 2.河南城建學(xué)院 交通運輸工程學(xué)院，河南 平頂山467036)

通過對粉膠比為0.62、0.82、1.02、1.22和1.42的瀝青膠漿在20℃、30℃、40℃和50℃條件下進(jìn)行動態(tài)頻率掃描試驗，研究了不同粉膠比及試驗溫度條件下瀝青膠漿復(fù)模量、抗車轍因子和相位角的變化規(guī)律。基于分?jǐn)?shù)階導(dǎo)數(shù)理論，建立了Nutting蠕變方程與經(jīng)典分?jǐn)?shù)階導(dǎo)數(shù)Abel黏壺蠕變模型之間的關(guān)系，從而明確了Nutting蠕變方程各參數(shù)的物理意義。對分?jǐn)?shù)階Riemann-Liouville算子黏彈性蠕變本構(gòu)模型的動態(tài)力學(xué)響應(yīng)進(jìn)行分析，提出了利用瀝青膠漿動態(tài)頻率掃描試驗結(jié)果確定蠕變本構(gòu)模型中參數(shù)A值和γ值的新方法。

瀝青膠漿;粉膠比;動態(tài)剪切流變試驗;分?jǐn)?shù)階導(dǎo)數(shù);Nutting蠕變方程;變形因子;黏彈因子

按照膠漿理論，瀝青膠漿是構(gòu)成瀝青混合料的微分散系，主要起填充細(xì)集料空隙和黏結(jié)作用[1]，其組成結(jié)構(gòu)決定了瀝青混合料的高低溫性能。瀝青膠漿是典型的黏彈性材料，常用經(jīng)典的元件模型描述其力學(xué)性能[2]。元件模型主要由彈簧元件和黏壺元件構(gòu)成，不同的元件數(shù)量及組合形式構(gòu)成了 Maxwell、Kelvin、Kelvin-Voigt、廣義Maxwell、廣義Kelvin、廣義Kelvin-Voigt、Burgers等模型。近幾十年來，國外許多學(xué)者[3-6]以及國內(nèi)的劉麗等[7-10]先后采用傳統(tǒng)的力學(xué)理論及評價指標(biāo)對瀝青膠漿的性能進(jìn)行了深入研究。但傳統(tǒng)的力學(xué)模型存在參數(shù)多、參數(shù)確定不便和不能精確描述材料的力學(xué)特性等缺點。

近年來，分?jǐn)?shù)階導(dǎo)數(shù)模型因為能夠克服經(jīng)典整數(shù)階微分模型與試驗結(jié)果吻合不好的缺點而被廣泛應(yīng)用于各個工程領(lǐng)域[11-13]。Celauro等[14-15]依據(jù)冪函數(shù)經(jīng)驗方程構(gòu)建了以分?jǐn)?shù)階Riemann-Liouville算子描述黏彈性材料黏彈性應(yīng)變的蠕變本構(gòu)模型，并對純?yōu)r青的力學(xué)特性進(jìn)行了研究。李銳鐸等[16]提出了改進(jìn)的分?jǐn)?shù)階導(dǎo)數(shù)經(jīng)驗蠕變模型并利用該模型來描述瀝青膠砂的蠕變特性。

鑒于瀝青膠漿組成的微分散系對瀝青混合料的性質(zhì)起著重要影響，本文結(jié)合分?jǐn)?shù)階導(dǎo)數(shù)模型優(yōu)勢建立了Nutting蠕變方程[17]與經(jīng)典分?jǐn)?shù)階導(dǎo)數(shù)Abel黏壺蠕變模型之間的關(guān)系。以不同粉膠比的瀝青膠漿為研究對象，研究不同溫度條件下的瀝青膠漿動態(tài)力學(xué)特性，探討Nutting蠕變方程各參數(shù)表征的物理意義，并提出利用動態(tài)剪切流變試驗結(jié)果進(jìn)行一元線性回歸分析來確定蠕變方程的參數(shù)，為進(jìn)一步深入研究瀝青膠漿在不同溫度條件下的動態(tài)流變性能提供借鑒。

1　分?jǐn)?shù)階導(dǎo)數(shù)黏彈性本構(gòu)模型

黏彈性高聚物材料的蠕變服從Nutting公式[17]：

lgε(t)=lgA+γlgt+lgσ,0<γ<1

(1)

其蠕變?nèi)崃?/p>

J(t)=Atγ,0<γ<1

(2)

式中：A、γ為反映材料性能的材料參數(shù)，可通過試驗確定。

Celauro等[18]依據(jù)冪函數(shù)經(jīng)驗方程構(gòu)建了以分?jǐn)?shù)階Riemann-Liouville算子描述黏彈性材料黏彈性應(yīng)變的蠕變本構(gòu)模型如式(3)所示

(3)

Koeller等[19]以Abel黏壺元件模型構(gòu)建的黏彈性材料的蠕變?nèi)崃繛?/p>

0≤γ≤1

(4)

式中：Γ(1+γ)為1+γ的gamma函數(shù)；φ(t)為蠕變?nèi)崃繂挝浑A躍函數(shù)，即

當(dāng)γ=1時，該元件為Newton黏壺，代表理想流體，γ=0時，該元件變?yōu)閺椈稍?，代表理想固體。

2　分?jǐn)?shù)階導(dǎo)數(shù)黏彈性本構(gòu)模型動態(tài)力學(xué)特性

在交變應(yīng)力荷載作用下，基于分?jǐn)?shù)階導(dǎo)數(shù)的黏彈性材料的本構(gòu)方程[15]為

(5)

式中：G*(ω)=G1(ω)+iG2(ω)；G*(iω)為動態(tài)剪切模量，G1(ω)為存儲模量，G2(ω)為損耗模量。

將cγ=AγΓ(γ)代入式(5)并取絕對值可得

(6)

3　試驗材料及試驗條件確定

3.1試驗材料

瀝青為重交70號石油瀝青，基本技術(shù)性質(zhì)測試結(jié)果如表1所示。填料為石灰石礦粉，對礦粉試樣進(jìn)行X衍射分析(如圖1所示)，主要包含了CaC03、CaMg(C03)2和SiO2衍射峰。其密度為2.814 g/cm3，親水系數(shù)為0.86。參照實際工程瀝青混合料配合比中粉膠比，使用高速剪切乳化機制備粉膠比分別為0.62、0.82、1.02、1.22和1.42的瀝青膠漿。

表1　基質(zhì)瀝青基本技術(shù)性質(zhì)

圖1　礦粉XRD衍射圖譜Fig.1 XRD diffraction pattern of mineral powder

3.2線黏彈性范圍確定

剪切流變儀測試原理是基于流變學(xué)線黏彈性假設(shè)的基礎(chǔ)上，對基質(zhì)瀝青及粉膠比為1.42的瀝青膠漿進(jìn)行應(yīng)變掃描試驗，結(jié)果如圖2所示。可以看出：隨著粉膠比的增大，線黏彈性應(yīng)變范圍逐漸降低，同一粉膠比條件下，溫度越低，線黏彈性范圍也越小。因此，為保證所有試樣均能處于線黏彈性范圍內(nèi)，并考慮試驗數(shù)據(jù)分析方便，最終確定在進(jìn)行頻率掃描時的應(yīng)變幅值為0.1%。

4　試驗結(jié)果及分析

4.1瀝青膠漿動態(tài)頻率掃描試驗結(jié)果及分析

瀝青膠漿頻率掃描試驗結(jié)果如圖3所示?？梢钥闯觯孩?在雙對數(shù)坐標(biāo)中，同一溫度條件下，任一粉膠比瀝青膠漿的復(fù)數(shù)剪切模量隨著角頻率的增加而線性增加；② 在不同溫度條件下，同一粉膠比瀝青膠漿的復(fù)數(shù)剪切模量均隨著溫度的增加而減小，說明其抗變形能力隨著溫度的升高而降低。而對于不同粉膠比，復(fù)數(shù)剪切模量隨著粉膠比的增大而增大。說明礦粉含量的增大提高了瀝青膠漿的抗變形能力。

圖2　瀝青膠漿應(yīng)變掃描曲線圖Fig.2 Asphalt mortar dynamic shear strain scan

瀝青膠漿的抗車轍因子(G*(ω)/sinδ)[20]與粉膠比和溫度關(guān)系如圖4所示?？梢钥闯觯孩?除40℃條件下粉膠比為1.42的瀝青膠漿的抗車轍因子比粉膠比為1.22的瀝青膠漿略低外，其它不同溫度條件下，隨著粉膠比的增大，瀝青膠漿的抗車轍因子均有所增大。說明從瀝青混合料的高溫穩(wěn)定性考慮，宜選用較大的粉膠比，但并非粉膠比越大越好，過大的礦粉含量不但會增大瀝青用量，增加工程造價，同時礦粉也不宜均勻分散在瀝青混合料中，影響瀝青混合料整體結(jié)構(gòu)的高溫穩(wěn)定性；② 當(dāng)粉膠比增大，隨著溫度的降低，抗車轍因子急劇增大。說明溫度對瀝青膠漿的抗車轍因子有著重要影響。

圖3　瀝青膠漿動態(tài)頻率掃描曲線Fig.3 Dynamic frequency scanning curve of asphalt mortar

圖4　不同粉膠比不同溫度瀝青膠漿抗車轍因子Fig.4 The rutting factor of asphalt mortar at at different temperature and ratio

動態(tài)頻率掃描試驗得到的瀝青膠漿的相位角(振蕩角頻率為10 rad·s-1)如圖5所示?？梢钥闯觯孩?對同一粉膠比，相位角隨著溫度的升高逐漸增大，說明瀝青膠漿的黏性特性越來越強；② 隨著粉膠比的增大，各溫度條件下的相位角均逐漸降低，說明礦粉含量的增加增強了瀝青膠漿的彈性固體特性。

圖5　不同粉膠比不同溫度瀝青膠漿相位角Fig.5 The phase angle of asphalt mortar at different temperature and ratio

4.2基于分?jǐn)?shù)階導(dǎo)數(shù)瀝青膠漿動態(tài)力學(xué)特性分析

基于分?jǐn)?shù)階導(dǎo)數(shù)的黏彈性材料的本構(gòu)方程可以表示為式(6)，其中A和γ均為材料參數(shù)，對特定材料其為常數(shù)。為了確定材料參數(shù)A和γ，可將式(6)改寫為，

y=a+bx

(7)

圖6　不同粉膠比和不同溫度條件下擬合得到的A值Fig.6 The value of parameter A at different temperature and ratio

圖7　不同粉膠比和不同溫度條件下擬合得到的γ值Fig.7 The value of parameterγ at different temperature and ratio

由圖6可知，A值隨粉膠比的增大而減小；對任一粉膠比，隨著溫度的增高，A值逐漸增大，且礦粉含量越低，A值增加的幅度越大，說明A值對粉膠比和溫度的變化均比較敏感。因此，同溫度或同粉膠比條件下，A可以作為瀝青膠漿的高溫性能評價指標(biāo)，稱A為變形因子，A值越大，材料的高溫穩(wěn)定性越差。

由圖7可知，同一溫度條件下，粉膠比越大，說明瀝青膠漿的彈性成分越大，而γ值越小；同一粉膠比條件下，隨著溫度的升高，瀝青膠漿表現(xiàn)出來的黏性性能逐漸增加，越來越接近牛頓流體。因此，同溫度或同粉膠比條件下，可以稱γ為黏彈因子，γ值越大，瀝青膠漿的黏性性能越強，而彈性越弱。

5　結(jié)　論

(1) 在雙對數(shù)坐標(biāo)中，對任一粉膠比的瀝青膠漿同一溫度條件下的復(fù)模量隨著角頻率的增加而線性增加；在不同溫度條件下，瀝青膠漿的復(fù)模量均隨著溫度的增加而減??；對不同粉膠比的瀝青膠漿的復(fù)模量隨著粉膠比的增大逐漸增大。

(2) 瀝青膠漿的抗車轍因子隨著粉膠比的增大逐漸增大，且隨著溫度的降低急劇增大，而相位角則逐漸降低，說明溫度對瀝青膠漿的抗變形能力有著重要的影響；同一粉膠比條件下，隨著溫度的升高，抗車轍因子逐漸降低，而相位角則逐漸增大。

(3) 基于分?jǐn)?shù)階導(dǎo)數(shù)理論，對Nutting蠕變方程和Abel黏壺蠕變方程進(jìn)一步推導(dǎo)，建立了Nutting蠕變方程與經(jīng)典分?jǐn)?shù)階導(dǎo)數(shù)Abel黏壺蠕變模型之間的關(guān)系，并進(jìn)一步明確了Nutting蠕變方程的各參數(shù)的物理意義。提出A為變形因子，A值越大，材料的抗變形能力越差，而γ為黏彈因子，γ值越大，材料表現(xiàn)的黏性越強。對不同溫度下的瀝青膠漿的動態(tài)剪切流變試驗結(jié)果進(jìn)行了擬合分析，確定了Nutting蠕變方程中的變形因子和黏彈因子，提出一種確定Nutting蠕變方程參數(shù)的新方法。

[1] Puzinauskas V P. A study of mineral fillers in asphalt paving mixture[J]. Proceeding of the Association of Asphalt Paving Technologists,1961,30:439-440.

[2] 張肖寧. 瀝青與瀝青混合料的黏彈力學(xué)原理及應(yīng)用[M]. 北京: 人民交通出版社,2006:31-42.

[3] Van der Poel C. A general system describing the visco elastic properties of bitumens and its relation to routine test data[J]. Journal Appl Chem., 1954, 4(5): 221-236.

[4] Anderson D A, Goetz W H. Mechanical behavior and reinforcement of mineral filler-asphalt mixtures[J]. Proc., Association of Asphalt Paving Technologists, 1973, 42: 37-66.

[5] Palade L I, Attané P, Camaro S. Linear viscoelastic behavior of asphalt and asphalt based mastic[J]. Rheol Acta, 2000, 39: 180-190.

[6] Kim Y R, Little D N. Linear viscoelastic analysis of asphalt mastic[J]. Journal of Materials in civil Engineering?ASCE/MARCH/APRIL, 2004, 16: 122-132.

[7] 劉麗. 瀝青膠漿技術(shù)性能及評價方法研究[D]. 西安: 長安大學(xué), 2004.

[8] 張爭奇, 王永財. 瀝青膠漿對瀝青混合料高低溫性能的影響[J]. 長安大學(xué)學(xué)報:自然科學(xué)版, 2006, 26(2): 1-5.

ZHANG Zheng-qi, WANG Yong-cai. Influence of asphalt mortar on hot mix asphalt performance at high and low temperature[J]. Journal of Chang’an University:Natural of Science, 2006, 26(2):1-5.

[9] 耿韓, 趙永利. 瀝青膠漿路用性能的流變試驗研究[J]. 公路工程, 2009, 34(1): 84-89.

GENG Han, ZHAO Yong-li. Pavement performance study of asphalt-mineral filler mastic by rheology test[J]. Highway Engineering,2009,34(1): 84-89.

[10] 樊亮, 魏建明, 張玉貞, 等. 礦粉對瀝青膠漿的性質(zhì)影響及作用機理[J]. 建筑材料學(xué)報, 2014, 17(6): 1096-1101.

FAN Liang, WEI Jian-ming, ZHANG Yu-zhen,et al. Acting mechanism and performance of asphalt mortars by mineral filler[J]. Journal of Building Materials,2014,17(6):1096-1101.

[11] Bagley R L. Power law and fractional calculus model of viscoelasticity[J]. AIAA Journal,1989,27(10):1412-1417.

[12] 陳文, 孫紅廣, 李西成. 力學(xué)與工程問題的分?jǐn)?shù)階導(dǎo)數(shù)建模[M]. 北京: 科學(xué)出版社, 2010: 7-10.

[13] 李軍強, 劉宏昭,王忠民. 線性黏彈性本構(gòu)方程及其動力學(xué)應(yīng)用研究綜述[J]. 振動與沖擊,2005,24(2):116-122.

LI Jun-qiang, LIU Hong-zhao, WANG Zhong-min. Review on the linear the constitutive equation and its dynamics applicaions to visoelastic materials[J]. Journal of Vibration and Shock,2005,24(2):116-122.

[14] Celauro C, Fecarotti C, Pirrotta A, et al. Experimental validation of a fractional for creep/recovery testing of asphalt mixtures[J]. Construction and Building Materials, 2012, 36: 458-466.

[15] Fecarotti C, Celauro C, Pirrotta A. Linear viscoelastic behaviour of bitumen via fractional models[J]. Social and Behavioural Sciences, 2012(53): 450-461.

[16] 李銳鐸, 樂金朝, 馮師蓉, 等. 瀝青膠砂改進(jìn)分?jǐn)?shù)階導(dǎo)數(shù)冪函數(shù)經(jīng)驗蠕變本構(gòu)模型[J]. 建筑材料學(xué)報,2015,18(2):237-242.

LI Rui-duo, YUE Jin-chao, FENG Shi-rong, et al. Improved fractional order derivative empirical power function model for describing the creep of asphalt mortar[J]. Journal of Building Materials, 2015, 18(2): 237-242.

[17] 何平笙. 高聚物的力學(xué)性能[M]. 第2版. 合肥: 中國科學(xué)技術(shù)大學(xué)出版社, 2008: 104-105.

[18] Celauro C, Di Paolam M, Lo Presti D, et al. Modeling of the viscoelastic behavior of paving bitumen fractional derivatives[J]. Meccanica dei Materiali e delle Strutture, 2009, 1(2): 38-51．

[19] Koeller R C. Applications of fractional calculus to the theory of viscoelasticity[J]. Journal of Applied Mechanics, 1984, 51: 229-307.

[20] Anderson D A, Bahia H U, Christensen D W. Binder characterization and evaluation-Volume 3: Physical properties.Strategic Highways Research Program, Rep.No.SHRP-A-369[R], National Research Council, Washington, D.C. 1994: 43-79.

Dynamic performances of asphalt mortar based on fractional derivatives

LI Ruiduo1,2, YUE Jinchao1, LI Ruixia1, YAO Liyang2

(1. School of Water Conservancy and Environment, Zhengzhou Univ., Zhengzhou 450002, China;2. School of Transportation Engineering, Henan University of Urban Construction, Pingdingshan 467036, China)

Asphalt mortars were prepared with different ratios of filler bitumen of 0.62, 0.82, 1.02, 1.22 and 1.44. Then the dynamic behaviors of asphalt mortars at 20℃, 30℃, 40℃ and 50℃ were studied with a dynamic shear rheometer. The relationships among complex modulus, rutting resistance factor and phase angle udner different ratios of filler bitumen and test temperatures were studied. The relationships between Nutting creep equation and Abel sticky pot creep model of classic fractional order derivatives were established based on the theory of fractional derivatives. The dynamic responses of fractional order Riemann-Liouville operator viscoelastic creep constitutive model were analyzed. The model parametersAandγwere determined using the results of dynamic shear rheological tests.

asphalt mortar; ratio of filler bitumen; dynamic shear rheological tests; fractional derivative;Nutting creep equation; deformation factor; viscoelastic factor

河南省科技攻關(guān)計劃(112102310602)；教育部“長江學(xué)者和創(chuàng)新團(tuán)隊發(fā)展計劃”創(chuàng)新團(tuán)隊(IRT1186)；水利與交通基礎(chǔ)設(shè)施安全防護(hù)河南省協(xié)同創(chuàng)新中心(2012-6)

2015-09-21修改稿收到日期：2016-01-09

李銳鐸 男，博士生，講師，1980年生

樂金朝 男，博士，教授，博士生導(dǎo)師，1965年生

E-mail: yuejc@zzu.edu.cn

U414

A

10.13465/j.cnki.jvs.2016.15.008