沥青混合料蠕变试验数据处理与粘弹性计算
第37卷第6期
2007年11月
东南大学学报(自然科学版)
J0uRNAL
oF
V01.37No・6
soUT脏AsT
uNIvERs兀Y(NaⅢralsden∞Edi“on)
NoV.2007
沥青混合料蠕变试验数据处理与粘弹性计算
陈静云1
周长红1
王哲人1,2
(1大连理工大学土木水利学院,大连116023)(2哈尔滨工业大学交通科学与工程学院,哈尔滨150090)
摘要:为了获得沥青混合料粘弹性本构关系,并利用这种本构关系进行各种数值计算,结合贯入
试验中采集到的蠕变数据,采用Matlab软件对蠕变柔量进行拟台,得到了由广义Maxwell模型和Burge幅模型表示的粘弹性参数;针对ANsYS有限元软件的计算要求,推导了将其转化为
P唧ly级数形武的计算公式.通过实例计算表明:四参数Burge飓模型和六参数Ma】【well模型的
拟合相对误差分别小于0.7%和1.3%,利用Prony级数方法得到的计算结果与理论解误差不大干0.001%.Burgers模型比广义Ma】【well模型更能准确表达沥青混合料的本构关系;Prony级数的转化公式方法简单,计算精确.
关键词:沥青混合料;蠕变试验;粘弹性;ANsYs;Pmny级数中图分类号:u416.217
文献标识码:A
文章编号:lool一0505(2007)06一1091_05
DataprocessingandViscoelasticcomputationforcreeptest
of
ChenJingyuⅡ1
asphaltmixture
zhou
changhon91wallg
z11erenl・2
(‘scho“ofc~n棚嘶d娜“ch画懈d“g.D“锄u币V删ty缸儆h∞Iogy,D她116023.c虻柏)
(2sch∞lof&len∞柚dETl画ne州“g衄Co叫叫ni龃曲n,胁rbi玎h-sd加坼0fTkhnology.Harbjnl50090。0Ijna)
ob诅inthec0ⅡsⅡm虹veequationsof
is
Abstmct:Inorder
to
协num嘶al
SYS
asphalt“xtIlIe卸d
to
applymeme船ilyin.
byMatlabsoftware吐lroughanaIyzingcre印
of脚omod—dataacquiredillpenetrationtestofasphaltconcrete,ttlefea舳r,viscocla蛎c
els—BurgersaIldgeⅡemUzedMaxweUmodelareachjeved.A190,basedonmeDcquhmemofAN-
simul撕吼s,cfe印cornpliance
here
fi啦d
pa勰ters
software.卸alg嘶ⅡlIll
0f
isdeducedt0
tfansfo衄the
consdtIldveequacionsintoPronyseries.An
ex帅ple
Burgers
pene仃撕oⅡtes¨sofferedtoevaluatetlleerrorpercentage.Asa
model(fourparaI:[1ete巧)islesstIl柚0.7%,whneMaxwellmodel(sixpamme姆rs)lcss山en
no
result,relanVe唧巧。士
1.3%,扑dmesolunonerIⅢusingPronyseriesis
icalsolunon.So,it
moretllen0,00l%comparedwimthetheoret-
c舶k
proved血attheBurge飓modelismoreaccura【e
t11柚generaliz酣Ma】(well
fbnnulat0
modelinex口ressin叠Ⅱ1cviscoelasdcbchaviorofProny
asphalt
mixturc,aIld
tlle【ransf0加adon
s鲥esissimplea11defncient.
KeyworIls:asphaltconcrcte;cr∞ptest;viscoelastici可;ANSYS;Pronyse—es
目前,道路工程建设中,最主要的铺筑材料是沥青混凝土.它作为一种粘弹性材料,具有不同于其他弹性、弹塑性材料的本构关系,因而掌握其受力特点是较困难的.不过,随着粘弹性理论的不断
发展,以及ANsYS,ABAQus等有限元计算软件的不断完善,可以模拟和分析实际路面的受力状态‘1….
由于蠕变试验方法简单,控制方便,不少研究者做了大量的室内试验,以通过沥青混凝土材料的
收鞭日期:2007・04・16.
基金项目:国家自然科学基金资助项目(50578∞1).
作者简介:陈静云(1956一),女,教授.博士生导师,ch。哪y@dlul
甜u.印
应力应变关系,来揭示路面的真实受力状态.但是如何从这些应力应变数据中获得沥青混合料的粘弹性参数,进而利用计算机软件来模拟其受力特
1092
东南大学学报(自然科学版)
第37卷
点,仍较为困难.其原固有2条:一是数据拟合不够精确;二是拟合的数据公式与计算机模拟软件所需
要的数据公式不一致.这导致许多研究者有粘弹性
2
2.1
广义Maxwe¨模型参数
求解微分型本构方程系数
对式(1)一(3)施加零初始条件下的L印lace
数据反而采用弹性或弹塑性方法计算逼近路面的
受力状态.为能更好地利用这些数据,本文具体分
变换,经过整理可得如下递推公式H]:
析了适合于沥青混凝土的Bu蛳模型和广义
Maxwell模型的粘弹性参数拟合问题,并介绍了如何将这些参数转换为Prony级数,利用ANsYs对沥青混合料进行力学计算.
由于分析计算中计算量较大,本文采用计算效率高、操作简便的MaⅡab计算软件进行计算,得到了较好的效果.
当n=1时,“”=1,硝”=o,c5”=毛,cf”
=田I・
当n>1时,对于f=2,3,…,n有
Ajl。1’
Ap=
七=O
巫盟+∥)…‘
E.
t=1,2,…,f一2
刁,A::1’
E.
七=f—l
1一维粘弹性模型
1.1广义MaxwelJ模型
田。A5卜‘’+丑f卜”
七=1
广义Maxwell模型如图1(a)所示,其本构关口{f)=
系满足以下关系:
叩。A掣+争掣+丑…女=2,3,…,f一2
口=∑叽
矿l=Els+田l占
(1)(2)
可。AP+酗≯
d‘_1’
々;f.1
t=o
毋+》.=研g
f=2,3,…,n(3)
掣=
譬懈…t_1'z,…,H
延里
£
。:;
。
。
式中,n为Maxwell单元的个数.
于是可得微分方程
^一1
^
《
(^)广义地《w卸模型
岛
(b)踟。目∞模型
∑矶5‘孑=∑吼J‘F
式中,“=A:”(女=o,1,2,…,n一1);
,q”
女=0
(5)
—W卜刮.卜
图l一维粘弹性模型
1.2
qt={B:“’+c:“’
七=1,2。…,n一1(6)
t:。
【E一,
2.2
求解模型蠕变柔量由蠕变柔量的定义:
Burgers模型
此模型为四参数流体模型(见图1(b))其本
J(r)=£“l}l
Ls∑口。刘
Matlab软件求出.
2.3
『∑”‘]
(7)
构关系可表示为
盯+Pld+p2疗;鼋l占+日2言
(4)
…式中,£4[]算子表示Laplace逆变换,可以由
最小二乘法拟台粘弹性参数步骤
1)确定元件个数n,并选定弹性模量E及粘性系数饥的初始值.
2)代人式(6)和(7),利用Matlab中ilaplace函数,或采用LapJace逆变换的数值解法(Durbin
式中
”毒+警,舻毯i+T’p22珏
9-2
”扎”弩
第6期陈静云,等:沥青混合料蠕变试验教据处理与粘弹性计算
1093
算法‘“’)求得蠕变柔量,(f).
所以,由胡克定律和粘弹性问题对应原理,可3)调用MatIab优化工具箱中的lsqcurvefit函得单轴压缩时的Laplace解为
数,循环步骤1)、2).求得E。和吼的最优值.:一上。一皿丛!)!£∽3鱼:%2
F4n
3
Burgers模型参数
2———r一了5
由式(4).可以求出其蠕变柔量的表达式为Ⅲ
【手厂(s)+}歹(s)]吼
J(r)=击+去+乏(1一e也恤)
(8)
代入式(9)和(11)及x=耥G-得
最小二乘法拟合粘弹性参数步骤为:①选定
.,‘=[2,(r)+,(f)]/3=
弹性模量E:及粘性系数他的初始值;②调用可七而+砉+壶(卜e删’)
Matlab函数lsqcurvefit进行迭代求解,得到£。和
(12)
"。的最优值.为了区别一维粘性参数,用符号n代替町.对比4
ANSYS三维粘弹性计算方法
式(8)可见
巨
ql
由于ANsYs等工程软件在计算复杂的粘弹
“-2玎n而,”t
.
2了
性力学问题时本构模型采用Prony级数表示‘”,所
以本文只讨论此种形式.
G:=等,n:=等
假定体积应变为弹性,剪切应变呈流变性,则有
由式(10),Burgers模型的剪切模量,可表示口=P(f)_d£。=2缸口≥三;搿:2二乏”勺)协,
J
成Prony级数如下:
或
G(f)=G。+Go(口le叫71+02e嘶2)(13)
8#=,(f)}d5“
1
式中,G。=o(对沥青混合料),G0=G。;dl=
铲舶Ⅲ,=蠡j
p”
(G2/n2一a)/(口一a),rl=1/a;42=(G2/^2一声)/(a一芦),rt=1/p.
式中,%,%分别为应力应变的偏张量部分,口=
其中,口1,Ⅱ2,fl,r2分别对应于ANsYs中指定口√3,P。=s。/3分别为应力应变的球张量部分,
Prony本构模型时需输入的系数dI,d2,fI,屯.而¥号表示stiel日es卷积.
£。,p分别为各项同性材料初始弹性模量和泊松
对于Burgers模型有
比.
G∽=掣=南【(鲁一芦)‘
5算例
e虬(孚一。)e1】
沥青混凝土的单轴贯人试验被认为是评价沥\订,
u,
J
(-o)
青混合料抗剪切能力大小的有效方法”1.
“r)=击+去+去(卜一‰)(11)
本文在试验中采用的试件为直径150n1IIl,高
loo
mm,级配组成见表1”…,集料为石灰岩,沥青
式中
采用西太平洋产90号重交通道路石油沥青,沥青
“2——瓦~
p-+√p:一和2
质量分数为3,5%.利用MTs810试验机施加40卢:址≤业
kN恒载,压头直径50mm,有侧限,恒温45℃.取其中一组试验数据(见表2).
褒1试件混合料级配组成
孔径/mm19.016,0132
9.54.75236累积通过率/%
100
90.I
74
9
439
34.0
266
注I灏青质量分数为为3.5%;空隙率为3,7%
1094
东南大学学报(自然科学版)
表2单轴贯入蠕变试验数据拟合殛ANsYs计算结果(部分)
第37卷
5.1试验数据拟合
对蠕变试验数据分别进行广义Maxwell(元件
个数H=3)和Bu唱ers模型拟合,取泊松比为
0.25,拟合结果见表3.
表3试验数据拟台结果
I/B
图2蠕变柔量分别按2种模型拟合得到的曲线
用Burgers模型表示的Prony级数(取二阶)系数为
G0=665,580
6
图3
ANsYs计算模型(图中数宇为节点号)
25
m,nⅨY,1.0
1泊松比
口.=O.157288f1=7.42880142=0.842712
他,豫oNY,l,I,2,s肛AR!P∞“y级数TBTEm.45
TBDATA,,0.157288,7.42880l,O.842712.1065530942
f2=l065.530942
将计算所得的节点应变数据列在表2中,通过与Bu唱ers模型的理论解对比,可以看出:拟合数据与ANsYs计算结果完全一致,误差均在0.001%以下.
2种模型拟合曲线如图2所示,并分别计算相对误差(见表2).可以发现,这2种模型对于模拟沥青混合料的剪切变形都具有很高的精度,误差基本上在1%以内.通过比较还可以发现,四参数的
6结论
1)分析了沥青混合料蠕变试验数据的一般处
B“唱ers模型比六参数的Maxweu模型拟合精度更
高.
5.2
ANsYS计算验证
理过程,并对蠕变柔量采用广义M“well和B*
gers两种模型进行拟合,以获得对应的粘弹性参数.
选取可进行大变形计算的soLⅢ185单元,定
义如下模型(见图3),单元大小1m×lm×lm,共2个单元,12个节点,轴向均布荷载P=100N:
按照如下命令流定义材料性质:
2)结合ANsYs等工程软件计算粘弹性问题
时的本构关系,介绍了将粘弹性参数转化为Pmny级数的方法,并给出了ANsYs按此种方法计算的命令流.
Ⅷ,E)【,1,1663.9515
1初始模量
第6期陈静云,等:沥青混合料蠕变试验数据处理与粘弹性计算
3)通过沥青混合料的单轴贯入试验算例,验证了数据处理与模拟计算的正确性,并通过对比得出结论:Burgers模型对于模拟沥青混合料的蠕变柔量比用广义Maxwell使用参数少,而精度却更高.
参考文献(Referenc鹤)
[1]Lu
Y,w岫tPJ.Numcncal印PIo∞h0f“虻oel硒一
to巾l勰吐c帅alysis∞d姗m,”,1998,田(2):139一147.
f研∞phan
Jnix叫嚣[J],o唧“捃邶
[2]Erl(曲sSMJG,Ljux,scarp越A.3Dmitc
cl锄em
mDdelfbrQm曲,2002,2(3):305—330.
asphaltconc咒tc
re5P∞辨simul撕on[J].胁
J
[3]郭大智,任瑞渡.层状粘弹性体系力学[M].哈尔
滨:哈尔滨工业大学出版社,2001.[4]Durbjn
F.Nume^cal
il’vcrsi蚰0fLaPla∞们雎f0叽s:
如咖cient
iInprovement
t0
I沁bn玎and
Abatc’sⅡ圮删
[J].Tk07印Ⅳw如“"耐,1973,17(4):37l一376.
[5】张伟.六种数值L丑plas逆变换方法的评价[J].佳术
斯大学学报:自然科学版,2002,20(1):113一117.z11柚gWci.AIlevaluan加∞shmetIlodsbased叽
tIIenumericalinv盯sioⅡof
bplaceⅡ锄sfoⅡ∞n∞[J].
如HMlD,Jim吲踟fw耐J),:Ⅳ口n‘脚&im"副i“∞,
2002,∞(1):113一¨7.(inc撕肿辩)
[6]Pafk
s
w,schapery
R
A.M曲ods
0fin但ro优lvcrsion
betwe即1ine盯“sc∞㈦cmatedalfIlncdons.Panl—a
nuⅡ蛇一calme山0d咖MJ如“Ⅲ,0,勋f妇d硝&MⅡr目,1999,孙
based彻PI唧蛔jcs[J].血犯m4一(11):1653—1675.
[7]蒙上阳,唐国金,雷勇军,等.Bu增啪模型的参数获
取方法[J】固体火箭技术.2003,26(2):”一29.
Mengsh柚gy柚g,T8ngOuojin,LdYongjun,etalMe山0dofobcainingtl】ebu。ge阽modclpa咖nete倦[J].
m舢f可删耐胁胁跏胁fDgy,2003,26(2):27—
29.(in(=Iljnese)
[8]P0t、riⅡM曲Jo彘.comparis∞0fdnle-dom面n腼幛
el曲1皿tⅢodelli“g0fvi
sc_嘲枷c
snctIlresusing锄小
疗cient衄60nal
[D].Mon雠a1:蚴llVoigt—Kelv.mmodel
or
Pr衄y螂ics
uDiv啪时,2001.
[9]毕玉峰,孙立军.沥青混合料抗剪试验方法研究[J].
同济大学学报:自然科学版,2005,33(8):1036—
1040.
Bj
Y咖g,s曲圳岫.Rese黜h∞tcstm吐州of廿
pllalt
mixnl”’s
shea|ing
pr0Pcrdes[J].加"埘盯
r0喇f跏f"州母:胁M谢&把Me,2005,33(8):1036
—10柏.(inclline辩)
[10]曹建新.重载交通下级配碎石基层材料组成结构与
动力特性的研究[D].哈尔滨:哈尔滨工业大学交通科学与工程学院.2001.
沥青混合料蠕变试验数据处理与粘弹性计算
作者:作者单位:
陈静云, 周长红, 王哲人, Chen Jingyun, Zhou Changhong, Wang Zheren
陈静云,周长红,Chen Jingyun,Zhou Changhong(大连理工大学土木水利学院,大连,116023), 王哲人,Wang Zheren(大连理工大学土木水利学院,大连,116023;哈尔滨工业大学交通科学与工程学院,哈尔滨,150090)
东南大学学报(自然科学版)
JOURNAL OF SOUTHEAST UNIVERSITY(NATURAL SCIENCE EDITION)2007,37(6)9次
刊名:英文刊名:年,卷(期):被引用次数:
参考文献(10条)
1. Lu Y;Wright P J Numerical approach of viscoelastoplastic analysis for asphalt mixtures 1998(02)2. Erkens S M J G;Liu X;Scarpas A 3D finite element model for asphalt concrete response simulation2002(03)
3. 郭大智;任瑞波 层状粘弹性体系力学 2001
4. Durbin F Numerical inversion of Laplace transforms:an efficient improvement to Dubner and Abate'smethod 1973(04)
5. 张伟 六种数值Laplas逆变换方法的评价[期刊论文]-佳木斯大学学报(自然科学版) 2002(01)6. Park S W;Schapery R A Methods of interconversion between linear viscoelastic materialfunctions.Part Ⅰ-a numerical method based on Prony series 1999(11)
7. 蒙上阳;唐国金;雷勇军 Burgers模型的参数获取方法[期刊论文]-固体火箭技术 2003(02)
8. Potvin Marie-Josée Comparison of time-domain finite element modelling of viscoelastic structuresusing an efficient fractional Voigt-Kelvin model or Prony series 2001
9. 毕玉峰;孙立军 沥青混合料抗剪试验方法研究[期刊论文]-同济大学学报(自然科学版) 2005(08)10. 曹建新 重载交通下级配碎石基层材料组成结构与动力特性的研究[学位论文] 2001
本文读者也读过(9条)
1. 郑健龙. 吕松涛. 田小革 沥青混合料粘弹性参数及其应用[期刊论文]-郑州大学学报(工学版)2004,25(4)2. 张丽娟. 陈页开. Zhang Li-juan. Chen Ye-kai 重复荷载下沥青混合料变形的粘弹性有限元分析[期刊论文]-华南理工大学学报(自然科学版)2009,37(11)
3. 王锦河. WANG Jino-he 沥青混合料高温蠕变试验方法研究[期刊论文]-重庆交通大学学报(自然科学版)2009,28(1)
4. 王学军. 岳宏智. 周长红 沥青混凝土蠕变试验数据处理与粘弹性参数的获得[期刊论文]-山东交通科技2007(3)5. 吕松涛. 田小革. 郑健龙. LV Song-tao. TIAN Xiao-ge. ZHENG Jian-long 沥青混合料粘弹性参数的测定及其在本构模型中的应用[期刊论文]-长沙交通学院学报2005,21(1)
6. 周志刚. 钱国平. 郑健龙 沥青混合料粘弹性参数测定方法的研究[期刊论文]-长沙交通学院学报2001,17(4)7. 郑健龙. 吕松涛. 田小革. ZHENG Jian-long. LU Song-tao. TIAN Xiao-ge 基于蠕变试验的沥青粘弹性损伤特性[期刊论文]-工程力学2008,25(2)
8. 张熙颖 沥青混合料低温抗裂性能蠕变试验研究与分析[期刊论文]-中外公路2007,27(6)9. 魏建明 沥青混合料蠕变试验研究[期刊论文]-公路与汽运2007(4)
引证文献(9条)
2. 叶永. 杨新华. 陈传尧 沥青砂粘弹性模型参数的试验研究[期刊论文]-中外公路 2009(4)3. 刘广军 沥青路面面层剪应力的粘弹性有限元分析[期刊论文]-交通标准化 2009(6)
4. 张丽娟. 陈页开 重复荷载下沥青混合料变形的粘弹性有限元分析[期刊论文]-华南理工大学学报(自然科学版)2009(11)
5. 朱云升. 向会伦. 张谢东. 张彦 沥青路面结冰条件下抗滑性能[期刊论文]-武汉理工大学学报(交通科学与工程版) 2012(1)
6. 龚兰云. 高永毅. 张厚安 多孔材料的力学模型和粘弹性性能研究[期刊论文]-湖南科技大学学报(自然科学版)2010(1)
7. 应荣华. 侯昭光. 庞海峰 细集料粒径对沥青砂浆黏弹性参数影响程度的灰熵分析[期刊论文]-重庆交通大学学报(自然科学版) 2010(3)
8. 王东栋. 许建聪. 冯兆祥 泰州长江大桥群桩基础工后沉降数值分析[期刊论文]-解放军理工大学学报(自然科学版) 2009(6)
9. 高永毅. 龚兰云. 张厚安. 张翠娟 人体桡骨支架力学模型及承载性能[期刊论文]-力学与实践 2010(2)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_dndxxb200706030.aspx