第23卷第3期 兰州工业学院学报 V01.23 No.3 Jun.2016 2016年6月 文章编号:1009—2269(2016)03—0036—04 沥青路面车辙数值模拟分析 田 静 (山西省交通科学研究院,山西太原030006) 摘要:为研究沥青路面车辙形变的规律,采用ANSYS有限元分析软件,结合Creep蠕变模型,对 沥青混合料中面层及全厚典型沥青路面的车辙变形进行数值计算,并与室内测试结果进行对比 分析.研究结果表明:单层板和全厚沥青路面的车辙变形数值计算结果与室内试验测试结果均随 着荷载作用的增加变形逐渐增大,但两者之间的初期变化差异性较大,差异性会随荷载作用的增 加而减小;ANSYS有限元中的蠕变模型能够有效地模拟沥青混合料路面车辙变形规律.研究成果 可为沥青路面结构设计及路面车辙病害防治提供技术支持. 关键词:沥青路面;车辙;蠕变模型;数值分析 中图分类号:U416.2 文献标志码:A 计算了沥青路面永久变形值.Erkens和Blab等 J 人则采用有限元三维方法分析路面力学性能.张久 鹏等 运用有限元建立了路面结构二维模型,分 0 引言 沥青路面的永久变形是高温和重载交通地区 所面临的主要问题之一,也是各国各地区所讨论的 析了沥青路面柔性基层的车辙变化.张丽娟 结合 D—P蠕变模型,利用粘弹性有限元对沥青混合料的 静载情况下的蠕变进行了分析.汪凡 采用John。 SOIl—Cook粘塑性模型分析了沥青混合料的温度效 研究热点之一L1。J.目前,中、美等国的沥青混合料设 计标准将路面结构的抗车辙能力作为衡量路面结构 性能的指标之一,由于影响沥青路面车辙的因素是 多样的,复杂的,并且沥青材料属于粘弹性材料,给 应、应变硬化等力学特征. 本次研究结合我国国内主要的沥青路面结构 形式,采用蠕变模型,运用数值计算方法对沥青路 面车辙进性模拟分析,并与试验和室外车辙结果进 沥青混合料车辙的相关研究带来了很多困难,研究 成果依然存在很多理想状态和不完善的情况. 沥青混合料性能与温度、时间、应力之间的关 系明显,属于非线性弹塑性材料.由于沥青混合料 的性能优劣差异,在高温及重、超载车辆的作用下 行对比分析,研究成果能够为沥青路面结构设计及 路面车辙病害的防治方法提供理论支持,也为相关 研究提供参考. 沥青路面产生不同程度的永久变形,从而出现车辙 痕迹,给路面行驶舒适性和安全性带来不利.因此, 沥青混合料永久变形的研究对路面结构设计及车 辙病害防治的意义非凡.目前在数值分析方面,路 1 模型建立 1.1 沥青结构层及模型参数选取 面车辙的相关研究成果较为丰富,例如wu_3 利用 ANSYS有限元软件,结合蠕变模型和(D—P)塑性模 型预测沥青混合料路面结构的永久变形.Hua 利 用ABAQUS有限元,将路面假设为平面应变问题, 收稿日期:2016—03-21 对于车辙的影响,不同的沥青结构面层所承担 的作用也不一样,在选取路面结构材料层的车辙情 况时,需要从各个层面的受力情况进行考虑,本文 所采用沥青路面结构层由上至下依次为AC一13, 作者简介:田静(1986一),女,山西长治人,工程师,硕士 ·38· 兰州 工业学 院学报 第23卷 研究结果与数值计算结果对比分析得到AC.20车 辙试验数据结果,见表2所示 表2室内试验与数值结果数据 l2 黑l o 8 。 数值分析结果的差异性较大,随着加载时间的增 加,室内试验结果与数值计算结果之间的差距逐渐 减小.初期结果差距较大的影响因素复杂,主要原 因在于模型、模型参数以及边界条件的选择都会对 试验结果产生较大的影响,尤其是时间参数的选 等效静载加栽时}f_】【/s 蒺4 誉2 择,时问越长,试验结果影响会越小,整体来看,室 图4单层板室内与数值车辙曲线 蔷 内试验与数值计算所得到的最终永久变形大小基 本一致. 。。未。20O 1 50 3全厚沥青路面对比分析 上述试验主要是针对中面层AC一20沥青混合 料,然而实际上单层板与全厚沥青路面两者的边界 条件不一样,尺寸大小和路面结构上也完全不同, 因此,为了更为真实的模拟路面车辙情况,对全厚 喜r<1。0。0 O 图5室内与数值车辙结果相对变异曲线 图4为室内车辙试验结果与有限元车辙数值 结果的对比曲线,从图4可以看出,运用ANSYS有 限元的Creep蠕变模型模拟得到的中面层AC一20 沥青路面进行车辙数值分析,并与室内相对车辙结 果进行对比. 本文选取ANSYS12.0中的时间硬化蠕变模型 的永久变形的计算结果与通过室内车辙板试验得 到的数值结果变化趋势具有一致性,随着等效静载 来分析全厚沥青路面车辙变形情况,由于建模只考 虑了路面面层结构,未考虑土基层的影响,将等效 静载作用时间折算成荷载作用次数.相对车辙的数 值计算结果与室内测试结构如图6所示. 加载时间的增加,变形逐步增加,曲线斜率逐渐减 小,最终有趋于稳定的趋势,曲线整体呈现出非线 性变化.图5为室内与数值车辙结果相对变异曲 线,相对变异率表示数值计算结果与室内试验车辙 试验结果的百分比值,相对变异率用来表示数据之 间的差异性.从图5可以看出,等效静载加载初期, 试验结果的相对变异率较大,高达300%左右,随 着等效静载加载时间超过100s后,相对变异率逐 渐降低,小于150%,并逐渐趋于100%.表明在等效 静载加载的初期阶段中,室内车辙试验变形结果与 筒载作翊次数/1 0。 图6全厚板室内与数值相对车辙曲线 从图6中可以看出,全厚沥青路面相对车辙的 第3期 田 静:沥青路面车辙数值模拟分析 ·39· 数值计算结果与实测结果曲线的变化趋势一致,相 参考文献: [1] 杨博.基于有限元方法的沥青路面车辙影响因素分析 及其应用研究[D].西安:长安大学,2010. [2] 王陆平,闫翠香.济青高速公路沥青路面维修后车辙 对车辙均随着作用次数的增加而增大,最终趋于稳 定,但有限元数值计算结果与室内实测结果之间存 在差距,数值计算结果普遍小于实测计算结果,但 随着作用次数超过2x 10 次之后,两者之间的差异 深度数值模拟[J].山西交通科技,2011(4):l0—12. [3] Wu Zhong.Finite element simulation of rutting on super— 逐渐减小.究其原因,在荷载作用的初始,由于沥青 混合料内部存在空隙率,在荷载的作用下,混合料进 pave pavements[D].Manhattan:Kansas State Universi— ty,2001. 行颗粒重排,内部空隙率减小,早期车辙变形增加, 但随着荷载作用次数的增加,这种变形速率逐渐减 小,并趋于稳定.此外,研究还表明ANSYS有限元软 [4]Hua Jian—feng,White T.A study of nonlinear tire contact pressure effects on HMA rutting[J].The International Journal of Geomechanics,2002,2(3):353—376. 件能够有效地模拟沥青混合料路面车辙变形规律. [5]Erkens SMJG,Liu X,Scavpas A.3D finite model for as— phah concrete respone simulation[J].The International Journal of Geomechanics,2002,2(3):305—330. 4 结语 本文利用ANSYS有限元分析软件,结合蠕变 模型,对沥青混合料中面层AC一20及全厚典型沥 [6] Blab R,Ha ̄ey J T.Modeling measured 3D tire contacts tress in a viscoelastic FE pavement model『J].The Inter— nationa1]ournal of Geomechanics,2002,2(3):271—290. 青路面的车辙变形进行数值计算,并与室内试验结 果进行对比分析.研究得到了以下主要结论: [7] 张久鹏,黄晓明,王晓磊.基于粘弹塑性理论的沥青路 面车辙分析[J].公路交通科技,2007(10):20—24. [8] 张丽娟,张肖宁,陈页开.沥青混合料变形的粘弹塑性 本构模型研究[J].武汉理工大学:科学与工程版, 2011(4):289—292. 1)单层板和全厚沥青路面的车辙变形数值计 算结果与室内试验测试结果的变化趋势一致,随着 荷载作用的增加变形逐渐增大,同时数值计算结果 与室内试验测试结果初期变化的差异性较大,差异 性随荷载作用次数增加逐渐减小. 2)ANSYS有限元中的蠕变模型能够有效地 模拟沥青混合料路面车辙变形规律,可为车辙变形 研究提供较好的技术手段. [9] 汪凡.基于流变学本构模型和动力有限元分析的沥青 路面车辙计算[D].重庆:重庆交通大学,2009. [10] 张久鹏,黄晓明,高英.沥青混合料非线性蠕变模型 及其参数确定[J].长安大学学报:自然科学版,2009 (2):24—27,55. 3)研究结果对于沥青路面结构设计及路面车 辙病害的防治方法提供理论支持,可进一步开展相 [11] 蔡红兵,曾志远,王福建.循环荷载作用下沥青路面 变形的粘弹塑性分析[J].公路交通科技:应用技术 版,2012(1O):50—52,55. 关研究,以细化数值计算的程度,扩大其应用价值. Numerical Analysis of Asphalt Pavement Rut TIAN Jing (Shanxi Transportation Science Research Institute,Taiyuan 030006,China) Abstract:In order to study the rutting deformation of asphalt pavement,the rutting deformation of asphalt mix— ture in the surface layer and the full thickness of the asphalt pavement was calculated by the ANSYS finite soft— ware and creep mode1.The results are compared with the indoor test results.Research results show that:rutting deformation numerical calculation results and laboratory test results both prove the deformation increases with more loads,and the difference of the change shows obvious at the beginning and becomes less gradually.The creep model can simulate the rutting deformation of asphalt mixture effectively.The research results can provide technical suppofl for asphalt pavement structure design and rut disease prevention. Key words:asphalt pavement;track;creep deformation model;numerical anylasis (责任编辑:曾贤灏)