沥青混凝土弯曲疲劳试验疲劳损伤分析.docx
- 文档编号:26291237
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:30
- 大小:44.32KB
沥青混凝土弯曲疲劳试验疲劳损伤分析.docx
《沥青混凝土弯曲疲劳试验疲劳损伤分析.docx》由会员分享,可在线阅读,更多相关《沥青混凝土弯曲疲劳试验疲劳损伤分析.docx(30页珍藏版)》请在冰豆网上搜索。
沥青混凝土弯曲疲劳试验疲劳损伤分析
第42卷第6期
2011年6月中南大学学报(自然科学版)JournalofCentralSouthV01.42No.6University(ScienceandTechnology)June201l
沥青混凝土弯曲疲劳试验疲劳损伤分析
周志刚,张清平,袁秀湘
(长沙理工大学道路结构与材料交通行业重点实验室,湖南长沙,410004)
摘要:
针对现有沥青混合料疲劳损伤试验时直接依据刚度或模量计算损伤值而无法真实反映材料微观损伤特性的不足,提出应用反分析方法以得到沥青混凝土的疲劳损伤演化规律。
其步骤是:
根据沥青混凝土疲劳损伤模型,得到三点弯曲梁试件疲劳损伤力学理论经典解:
根据疲劳试验数据确定沥青混凝土材料疲劳损伤模型中的特征参数。
应用非线性有限元方法,模拟计算并分析试件疲劳损伤过程中的弯拉应力、挠度、损伤变量、裂纹扩展速率等的变化规律,并预测沥青混凝土试件疲劳寿命及失稳断裂时的裂纹长度。
研究结果表明:
疲劳寿命的数值模拟结果与试验结果较吻合,证明了所提出的沥青混凝土疲劳损伤模型的合理性和有效性。
关键词:
沥青混凝土;弯曲疲劳试验:
非线性疲劳损伤;仿真
中图分类号:
U416.217
-・’文献标志码:
A’・文章编号:
1672—7207(2011)06-1743-09^’_・・●1一■。
attguedamageanalysisolaspahltmixturespecimen
duringcyclicbendingtest
ZHOUZhi—gang,ZHANGQing-ping,YUANXiu-xiang
(KeyLaboratoryofRoadStructureandMaterial,MinistryofTransport,
ChangshaUniversityofScienceandTechnology,Changsha410004,China)
Abstract:
Basedonthefactthatdamageindexiscalculatedusingstiffness
toormodulusdirectlyinfatiguedamageexperimentforasphaltmixtureanditisunable
back—analysis
classicrevealtherealityofasphaltmixture’Smicro—damagecharacteristic,themethodwassuggestedtodeterminethefatiguedamagecharacteristicofasphaltmixture.Atfirst,thebeamspecimenduringthreepointsbendingtestfatiguedamagetheoricalsolutionofwasdeducedusingthefatiguedamage
basedonmodelsuggestedforasphaltmixture,andthecharacteristicparametresinthemodelweredeterminedvariationsofbendingtheresultsoffatiguetest.Thestress,deflection,damageindex,crackexpandingvelocity,etc.inthespecimenduringfatiguedamageweresimulatedbyusingnonlinearfilliteelementmethod,andthefatiguelifeandthecracklengthwhilebreakingwerepredi【cted.Theresultsshowthatthenumercalfatiguelifeiswellidenticalwithexperimentalresults,whichprovesthereasonablenessandeffectivenessofthefatiguedamagemodelsuggestedforasphaltmixture.
Keywords:
asphaltmixtnre;cyclicbendingtest;nonlinearfatiguedamage;simulation
沥青路面的疲劳性能是表征沥青路面耐久性的重
要性质之一。
鉴于沥青混合料以及路面结构的疲劳破
坏实质上是材料局部出现劣化损伤并逐渐扩展、发展
成为裂缝直至断裂破坏的过程,为此,需要在试验基
收稿日期:
2010--07-22:
修回日期:
2010--10-09础上进一步应用疲劳损伤力学理论与方法揭示材料及结构疲劳损伤破坏的演化过程,认识疲劳破坏的机制。
关于沥青混合料疲劳损伤,一般基于2类损伤理论,即弹性损伤理论【卜51和黏弹性损伤理论[6--11l。
这些研究
基金项目:
国家自然科学基金资助项111(50778027):
湖南省自然科学基金资助项I|(04JJ30331
通信作者:
周志刚(1966-),男,湖南长沙人,博士,教授,博士生导师。
从事道路工程结构与材料研究;电话;0731.85259256:
E-mail:
zzgcs@yahoo.c,01n011
万方数据
1744中南大学学报(自然科学版)
第42卷
中尽管应用损伤力学理论与方法,但损伤值一般用试件的刚度或模量的变化来定义,该参量具有平均意义,并不能真实地反映试件局部损伤及其发展过程。
例如弯曲疲劳试验研究中采用根据梁式试件挠度与弯矩计算得到的截面等效弯曲刚度来定义计算梁的损伤值,显然,这表示梁在整个截面平均的损伤值,忽视了同一断面上不同受力状态区域损伤的差异性,由此得到的疲劳损伤演化关系及其有关参数不能真实地反映材料微观损伤特性及其损伤演化过程。
其原因主要是:
在这些试验中,未能获取试件局部应变、微观结构变化等与经典损伤定义直接相关的信息。
为此,本文作者结合沥青混凝土梁式试件弯曲疲劳试验,采用反分析的方式,以提出真实体现沥青混合料损伤演化的模型。
其步骤是:
首先提出沥青混合料的非线性疲劳损伤演化模型;然后,应用解析与数值仿真分析的方法演绎疲劳试验过程,将所得结果与疲劳试验结果进行对比,以论证模型的合理性和有效性。
1
梁式试件疲劳损伤力学理论的经
典解
1.1疲劳损伤方程
由于一般疲劳试验过程只采集循环荷载及其次数、梁的挠度等,而无法直接获取损伤信息,故采用以下非线性疲劳损伤演化方程【12】来进行反分析:
竺=adN‘(忐1——2
l—ID卜D广
Il一∥}1
IlI
I—J、
7
(1)一
式中:
a’,P和g为材料的损伤特性参数,与载荷循环特征R及温度丁有关;损伤变量D采用Kachanov的经典定义,将经典的弹性损伤本构方程o=E(1-D)£代入式(1),并令a=a’F,有:
—dD—:
a占p(1一D)一9——=占‘ll—U)1
(2)tZ,
dⅣ
1.2梁式试件的拉压分区及有关力学量计算
假定沥青混合料梁式试件损伤前拉压弹性模量没有区别,即具有各向同性性质,但由于其损伤破坏体现为拉伸形式,故出现在梁式试件的受拉一侧。
假定受压一侧不会出现因损伤引起的力学性质的变化,故沥青混合料梁式试件x断面的拉压分界点zo(x)由断面轴向受力平衡确定:
其中:
乒-z/h:
毒o=Zo(x)/h:
啪;d黜。
J:
:
(1-t)0一D)dt=吉(古一1)2
(3)
工断面内的拉压应力可按下式计算:
万方数据
嘶)-[1_晰)】【zo(炉z)】等
(4)
在受压区,
嘶)_【zo(圹z】等
(5)
其中:
E7(x)为材料损伤后x断面的等效弯曲刚度;
7(x)为x断面的等效惯性矩。
m)铜(州扣f)2【1_∞力仲+詈(赤嘲
(6)
x断面内受拉区任一点的损伤值可根据最危险点
M的损伤值进行计算:
百dD(x,z)《型MM・警Zo・斟,[寄卜峨【
●,-
J(圳【l一%j~
1.3疲劳裂纹形成与扩展的疲劳寿命一=‘一・:
=———一,l—I
M
1.3.1疲劳裂纹形成阶段的疲劳寿命
试件未断裂前以M点表示最危险断面(梁跨中断面)内受拉区距中性轴最远的点(试件表面),锄和Du分别表示M点的正应变幅值与损伤值。
由式(2),M
点达到损伤值巩时的疲劳寿命(疲劳裂纹形成寿
命地为:
%=J:
:
1i1。
面Mzo)叩(1一%)9峨
(8)
记名义应埘:
竿:
器,如:
击砘记名义应力嘞’=—}2老》,如。
击醐3,
则
Ⅳ盯=C・(Oo’)一p
(9)
其中:
~世赤枷训9蛾](10,
材料特性损伤参数a与荷载、损伤无关,从而,
a=C’/C
(11)
其中:
c,一Eplc古c扣・训9蛾](12,
计算疲劳裂纹形成阶段的疲劳寿命的基本步骤见文献【12】。
1.3.2疲劳裂纹扩展阶段的疲劳寿命
当循环次数达到M,时,疲劳裂纹形成,试件进
第6期
周志刚,等:
沥青混凝土弯曲疲劳试验疲劳损伤分析1745
入疲劳裂纹扩展阶段。
采用与疲劳裂纹形成阶段一样的分段方式,即在计算疲劳裂纹扩展阶段的疲劳寿命时,将计算断面等分为m2段(含m2+l点),每段长度为△^勘砌2。
当循环加载次数Ⅳ达到^0时,试件表面点1的损伤值达到1,此时,相邻点2的损伤值达到D5n。
随着Ⅳ进一步增加,点2的损伤值达到D52)=1,产生损伤失效。
此时,所对应的循环加载次数增加量为:
AN2
2珐俨-l(兹九1_D{D)9dD2
(13)
其中:
z:
为产生裂纹长度al砀1l后试件有效断面内
在以点2为坐标原点的新坐标系下的拉压分界点坐标;J’为此时的等效惯性矩。
依此类推,第f-1点失效时,对于第i点,其损伤值D}1)为前f-l段开裂过程的累积损伤结果。
当其损伤值达到1时,所对应的循环加载增加量为:
ANi
2岛一铥九1一蚋9dDi(14)
疲劳裂纹扩展寿命AN与疲劳全寿命M分别为:
△Ⅳ=AN2+AN3+…+AⅣm+I
(15)Nf=Nc【+厶N
(16)
2
沥青混合料弯曲疲劳试验及疲劳损伤特性参数
采用三点弯曲疲劳试验测试方法。
分别选用3种
70号沥青制作中粒式密级配矿料沥青混凝土,试件采用轮碾成型后切割成规定尺寸的制作方法。
进行疲劳试验时,采用10Hz的正弦波作为加载波形,加载时
间为16nl¥,相邻波形间没有间歇时间。
试验所用的3种沥青混凝土的抗弯拉强度和弯拉回弹模量见表l。
根据弯曲疲劳试验结果,疲劳方程^r-国矽中的疲劳参数C和P回归值见表2。
表1
沥青混凝土的抗弯拉强度和弯拉回弹模量
Table1
Bendingstrenghthandmodulusofasphaltmixtures
注:
在测试抗弯强度实验中,A,B和C的品种数量分别为6,6和5个;在测试抗拉回弹模量实验中,A,B和C的品种数量分别为4,4和3个。
万方数据
表2疲劳试验结果
Table2
Resultsoffatiguetests
取表2中3种沥青混凝土疲劳试验结果的平均值,
即C=40467,p=3.45。
根据疲劳损伤裂纹形成阶段的
有关计算及式(11),得到沥青混凝土疲劳损伤特征参
数a=1.525x10-5。
3
三点弯曲试件疲劳损伤的有限元
分析
3.1疲劳损伤力学一有限元方法
采用平面应变条件下的八节点等参单元有限元方法,考虑损伤时弹性矩阵S为:
1-/z∥0s:
墨!
!
二望
(1+∥)(1—2/z)
∥
1-g
0
(17)
0
O—1-—2//
2
式(1)中应力改为等效应力变程crc=1层%2。
其
’二
中:
勋为应力偏量分量的变程,J旷=毋一F6。
盯√3。
在进行疲劳损伤分析与裂纹形成及扩展寿命预估时,应用文献【13】中所述步骤进行计算。
在模拟计算过程中,当某点累积损伤值达到1.0时即认为该点开裂,当连续多点损伤值均达到1.0时即认为它们形成了1条裂纹。
利用自编的有限元程序计算损伤不同增长步长AD为0.01,0.02,0.03,0.04和O.05时的疲劳损伤,结果表明均能取得较高的精度,故本文的有限元计算中,均取AD=0.02的增长量步长。
分别取5种集中荷载,使各个荷载初次作用于沥青混合料梁式试件时在其跨中断面底部产生的最大弯拉应力ao分别为0.4,0.6,0.8,1.0和1.2MPa。
然后,计算分析沥青混合料梁式试件在这5种集中荷载循环作用下的疲劳损伤开裂及其断裂过程,并与经典损伤理论分析和试验结果进行对比。
三点弯曲试件的长X宽X高为250
mmX30mfn
×35
mm,跨中长度为200mrn。
有限元网格划分为
756个八节点等参单元,共2405个节点。
1746中南大学学报(自然科学版)第42卷
3.2弯拉应力
3.2.1
梁底面弯拉应力
试件开裂前其底面弯拉应力吼底分布随疲劳荷载
作用的变化曲线见图1。
从图1可见:
随着疲劳损伤变量JD材的增大,梁式试件底面产生最大弯拉应力的位置逐渐地从跨中向两侧转移;当Djl4=1.0时,最大弯拉应力点的位置距跨中M点约为1.0cm。
并且梁式试件底面最大弯拉应力也由0.40MPa降低为O.25MPa;同时,M点的弯拉应力由最大值逐渐衰减,当DM=I.0时变为0MPa,见图2。
从图2可见:
其衰减的速度随损伤变量DM的增大而加快;此外,对于跨
0.4
O.3山
西
善0.2
蜓b
O.1
O125
75
25
25
75
125
x/mm
0.4
O.3
_(∽
山
蒌0.2蜓
b
0.1
O
125
:
/75
25
\
25
75
125
x/mm
O.4
O.3
.(c)
罡
善0.2
蜓b
0.1
O
125
:
/75
25
\
25
75
125
x/mm
上k:
(a)O;(b)0.5;(c)1.0
图1
a0=0.4
MPa时开裂前试件底面弯拉应力分布
Fig.1
Distributionofbendingstress
at
specimen
bouom
beforecrackingwhen00=0.4
MPa
万方数据
l一有限元;2一经典
图2
a0=0.4
MPa时试件跨中M点弯拉应力随损伤的变化
Fig.2
VariationofMpoint’Sbending
sUess
withdamagein
specimenmiddlewhen
00=0.4MPa
中M点的弯拉应力,其有限元计算结果(实线)与经典当开裂以后,随着裂纹长度的增加,梁底面最大降低,见图3。
并且裂缝附近弯拉应力零值区稍扩大。
O4
O3
.(b)
O2
-
叠善世0
OJ
0
125
。
门75
25
八
25
75
125
x/mm
裂纹长度/ram:
(a)5;(b)10
图3开裂后试件底面弯拉应力分布(cr0=o.4MPa)
Fig.3
Distributionofbending
stressat
specimen
bottomafter
cracking(oo=0.4MPa)
理论结果(虚线)很接近。
弯拉应力点位置继续向两侧转移,最大弯拉应力有所
第6期周志刚,等:
沥青混凝土弯曲疲劳试验疲劳损伤分析
1747
3.2.2跨中断面弯拉应力
在出现损伤前,跨中断面大部分截面弯拉应力几乎呈线性,只是在梁式试件项面集中荷载作用区附近出现异常现象,即压应力急剧增大(见图4)。
这与集中荷载区的应力集中现象有关。
该影响一直存在于损伤、开裂全过程。
随着梁底面产生损伤,M点的弯拉应力逐渐减小,断面上产生最大弯拉应力的位置逐渐向上转移;开裂以后,由于有效断面减小,弯拉刚度降低,在跨中断面上产生的最大弯拉应力和压应力均增大。
在裂纹出现、扩展全过程中,拉压分区点位置逐渐向上发展。
g
墨
(a)起裂前;(b)开裂后
图4跨中断面弯拉应力分布(ao=O.4MPa,有限元方法)
Fig.4
Distributionofbendingstressalong
cross
sectionof
specimenmiddle(cro=0.4MPa,FEM)
跨中断面弯拉应力经典理论值见图5。
对比图4和图5可见:
经典损伤理论解与有限元方法计算结果较接近,主要差别在于集中荷载附近受压区的应力分布。
由于经典损伤理论解未考虑集中荷载附近应力集中,故该受压区应力与有限元结算结果相差较大。
3.3梁的挠度
在梁底面开裂前,梁的跨中挠度矿随着疲劳损伤的发展而增大,但增大的速度逐渐趋缓(见图6(a))。
这
万方数据
g
戛
Z)M:
l—O;2—旬.5;3一1.0
图5跨中断面弯拉应力分布(aoffi0.4MPa,经典损伤理论解)
Fig.5
Distributionofbendingstressalongcrosssectionof
specimenmiddle(oo=0.4MPa,Classictheoricalsolution)
g
毒
瑙端
(a)开裂前;(b)开裂后
图6啊o---o.4MPa时跨中挠度随损伤.开裂的变化
Fig.6
Variationofdeflectionin
specimen
middlewith
damagearidcrackingwhen啊o--o.4MPa
是由于梁开始出现损伤时,损伤区主要在跨中,对该处弯拉刚度影响较大;随着损伤的进一步加深,损伤区范围扩大,受损区弯拉刚度降低对梁跨中挠度的影响趋于平缓。
当梁底面开裂后,跨中挠度随着裂纹的扩展增大,且增大的趋势逐渐加剧,当裂纹长度达到25—30mm
十自★{{报(自#H}*l
第42%
时,出现失稳断裂(见图6(b))。
3A损伤变量
开裂前损伤区主要位于粱底部和跨中两侧,损伤最大的位于跨中底,将从此处产生疲劳裂纹(见图7)。
起裂后,疲劳裂纹持续沿着跨中断面向上发展(见
图81。
_
酽
J・—■—●■●_一
-
●
+
,,=n7
图7
cro=04
MPa时开裂前拯伤场
Fig.7
Distribution
ofdamagebeforecracking
whenao=04MPa
:
-i
图8
d神4MPa和庐20nun时开裂后损伤场
Fig.8
Dis研bulionofdamageaftercracking
when㈣4MPaandF20ram
在粱上部及支点之间存在着大片无损和轻微损伤区域.这是由于在疲劳损伤中只考虑了受拉和剪切损伤,认为受压不会引起损伤。
3.5疲劳损伤寿命
在粱式试件底面开裂前,根据有限元方法计算结果绘制的疲劳曲线(即疲劳寿命与作用荷载应力关系曲线)与试验曲线非常接近(见图9)。
有限元方法与试验回归所得疲劳寿命见表3。
由表3可见:
二者最大
误差在177%以内;而疲劳寿命预估误差在30%以内,
均是可以接受的。
表3有限元方法与试验回归的疲劳寿命比较Comparmgtheth∞ritlcalfatigueliveswith
experimentalresults
方法
!
竺里塑竺竺1
04
06
08
l0
12
当粱式试件底面开裂后,在同一裂纹长度时,荷载应力大,则疲劳寿命短;当裂纹长度达到约20mm后.疲劳寿命趋于稳定.即将进入失稳扩展阶段(见图
万方数据
l
o)。
开裂前后疲劳寿命比较结果见表4。
由表4可见:
结合前面粱的挠度分析,若以裂纹长度为25mm作为失稳扩展分界,则失稳时的疲劳寿命比产生裂纹时的
疲劳寿命延长约O74倍。
lgⅣ
—试验结果;卜有限元计算结果
固9疲劳曲线
Fig.9
Fatiguecurve
1窟Ⅳ
doPMPa:
1--04:
2--08:
3一l2
图10疲劳寿命随裂纹长度的变化]ngao
Variationoffariguelifewithcrackl即gm
表4开裂前后疲劳寿命比较
Table4
Comparisonofmefatiguelivesbefomcracking
and
af时cmcking次
荷载应力n皿a
阶段
04
06
08
10
12
3.6裂纹扩展速率
根据不同荷载应力下三点弯曲粱式试件疲劳损伤
第6期周志刚,等:
沥青混凝土弯曲疲劳试验疲劳损伤分析1749
开裂、断裂过程的有限元计算,可以求出不同荷载应力、不同裂纹长度时的裂纹扩展速率见图ll。
从图ll可见:
随着裂纹长度的增大,裂纹扩展速率逐渐加大,
当裂纹扩展至20mm以后呈现出急剧变化趋势。
这结
其扩展寿命。
针对本文所分析的沥青混合料三点弯曲梁式试件,根据裂纹扩展速率的有限元计算结果,下面计算回归Paris公式中的参数彳和m。
da/dN=AK“(18、
其中:
da/dN为裂纹扩展速率,ram/次;K为应力强度
因子,
MPa・mmu2。
果与前面根据梁的挠度、疲劳寿命与裂纹长度关系推测结果相一致。
并且对于同一裂纹长度,在荷载应力较大时,其裂纹扩展速率相应增大。
3.7
为此,首先根据文献[14]提供的公式,计算各个荷载应力下不同裂纹长度时的I型应力强度因子K:
KI=0.720-owl/2[(1一a/w)弓一(1一a/w)3】172
(19)
Paris公式参数
在疲劳裂纹扩展分析中经常采用Paris公式预测
gg
专
已
董
已
g
董
已
§
喜
3
cro/MPa:
(a)O.4;(b)O.6;(c)0.8;(d)1.O:
(e)1.2
图ll裂纹扩展速率与裂纹长度关系
Fig.11
Relationship
be“㈣crack
expandingvelocityandcracklength
万方数据
1750
中南大学学报(自然科学版)第42卷
力:
j1・15—60(alw)2,o≤∥瞩o・05^2‘
(20)
I二V,
【1.0,a/w>l
o.05
一
其中:
口为裂纹长度;w为梁式试件高度;(70为荷载应力。
裂纹扩展速率和应力强度因子关系曲线见图12。
从图12可见:
当产生裂纹后,初始阶段裂纹扩展速度较小;当裂纹长度发展到6.25ranl后,进入稳定扩展阶段;当裂纹长度接近20111In时,裂纹扩展速度加快,将进入失稳扩展阶段。
由于Paris公式是反映稳定扩展阶段的疲劳裂纹扩展规律,因此,在采用双对数坐标下的最小二乘法回归(Paris公式,式(18))中疲劳断裂参数时,剔除裂纹长度大于20toni的点。
回归结果为:
A=I.832×10-5,m=1.938:
相关系数为0.870。
若剔除裂纹扩展初始阶段和20mill附近少量点,则可以获得相关性更强的回归结果。
疲劳裂纹扩散速率与应力强度因子的关系见图12。
从图12可见:
在稳定扩展阶段,裂纹扩展速率与应力强度因子具有良好的双对数线性关系,可以采用Paris公式预测裂纹扩展寿命。
a0/MPa:
●--0.4;◆—旬.6;;●--0.8:
▲一1.O;V—1.2
图12疲劳裂纹扩展速率与应力强度因子的关系
Fill2
Relationship
bctw∞nfatiguecrackexpanding
velocityandstresssHenghthindex
(1)根据疲劳损伤力学理论,推导出疲劳过程中的演化方程,以及疲劳裂纹形成阶段和扩展阶段疲劳万方数据
(2)应用疲劳损伤力学有限元方法,模拟计算了三
点弯曲梁疲劳损伤断裂过程,分析了试件疲劳损伤过程中的弯拉应力、跨中位移、损伤变量、裂纹扩展速率等的变化规律,并预测了沥青混合料试件失稳断裂时的裂纹长度。
疲劳寿命的数值模拟结果与试验结果较吻合,证明了在梁式试件三点弯曲疲劳试验情况下所提出的沥青混合料疲劳损伤模型的合理性和有效
性
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 沥青 混凝土 弯曲 疲劳 试验 损伤 分析