粘弹性状态下沥青混凝土路面层间剪力剖析.docVIP

粘弹性状态下沥青混凝土路面层间剪力剖析 　　【摘要】层间剪切滑移是山区及坡段公路半刚性基层沥青混凝土路面典型破坏形式之一。本文以数值模拟为主导，采用大型商用有限元软件ABAQUS建立三维路面模型对静载及半矩形波周期荷载作用下半刚性基层沥青路面的力学响应进行了三维粘弹性仿真分析，探讨了层间剪力随加载时间的变化规律。研究表明，沥青混合料的粘弹性能造成层间剪应力随时间而改变，柔性沥青面层间的剪应力随加载时间逐渐减小，沥青面层和半刚性基层间的剪应力随加载时间逐渐增大。对于半刚性基层沥青路面需要加强沥青面层与半刚性基层间的抗剪能力。 　　0 引言 　　我国普遍采用半刚性基层路面结构，半刚性基层是主要承载部分，沥青混凝土面层厚度相对较小，层间剪应力较大。对于山区公路及坡段公路半刚性基层沥青路面，层间剪切滑移是一种典型病害[1]。层间剪力是造成层间滑移的最直接原因。1989年，关昌余等人在测定层间抗剪模量时自行研制了直剪仪[3]，之后国内学者先后研制了各种室内试验设备对层间剪力进行了探讨。苏凯， 武建民等人在室内剪切试验的基础上应用计算软件分析了路面层间剪应力[1][2]。本文采用大型商用有限元软件ABAQUS分析了沥青面层处于粘弹性状态下半刚性基层沥青路面层间剪力响应，探讨了静载和半矩形波周期荷载下层间剪力随时间的变化规律。 　　1 ABAQUS粘弹性本构关系 　　ABAQUS的材料粘弹性本构关系采用了各向同性线性粘弹性应力与应变的遗传积分型式，该软件三维粘弹性关系表达式为[5][6]： 　　（1.1） 　　式中：、分别为偏应变和体应变张量；G和K分别表示剪切模量和体积模量；I为单位张量。 　　剪切松驰模量和体积松驰模量可表示为： 　　（1.2） 　　（1.3） 　　式中，和分别表示材料初始剪切模量和初始体积模量。 　　2 路面模型的建立与检验 　　选取路面结构如图1，由上至下各结构层依次为：上面层、下面层、基层、底基层、土基。上面层和下面层为沥青混凝土层，基层为半刚性基层。上面层、下面层、基层、底基层厚度分别取为：40mm、80mm、300mm、200mm。车辆正常行驶过程中，路面受到一向后的水平力，此水平力可表述为。路面平整度良好的情况下，车辆正常行驶时值为0.01～0.02[7]，取中值0.015。对路面结构施加标准轮压=0.7MPa，由式可得水平行车剪力0.0105MPa。双圆荷载半径δ=10.65cm，轮距31.95cm。周期性动载加载方式参照文献[8]，见图3。采用自动时间步长，荷载作用时间取30s。结构及荷载对称，取一半结构进行计算，见图2。忽略基层和土基土的粘弹性，假设基层以下材料均为线性材料，采用弹性模型。参照文献[9]，各层材料参数见表1。根据蠕变试验[10]得到试件各时刻的应变，如图4。在恒定加载应力下，根据材料瞬时本构关系 ，由应变曲线可得割线松弛模量曲线。假定泊松比为0.35，且松驰过程中不随时间改变，根据材料本构关系可以得到剪切松驰模量随时间变化的曲线。以初值为参照归一化处理之后可以得到μ的相对值随时间变化曲线，见图5。 　　图1 三维路面结构计算图示 图2 模型局部网格划分 　　图3 周期荷载 　　表1路面材料参数 　　上面层 下面层 基层 底基层 土基 　　模量（MPa） 180.55 120 1200 200 45 　　泊松比 0.35 0.35 0.2 0.3 0.4 　　图4应变-时间曲线 图5归一化剪切模量-时间曲线 　　3 计算结果分析 　　轮载作用位置附近上面层底面和下面层底面均受压，因此，纵向层间最大剪应力出现于水平力作用前方轮缘处，取前方轮缘处为纵向层间剪应力计算点；横向层间最大剪应力出现于车轮外轮缘处，见图1。 　　面层间最大剪应力随加载时间逐渐减小，基面层间最大剪应力随加载时间逐渐增大，见图6、7。加载初期，面层间剪应力要明显大于基面层间剪应力，随加载时间增长两条曲线差距在逐步缩小。纵向层间最大剪应力与横向层间最大剪应力相差不大。加载期间面层间纵向最大剪应力由0.154MPa减小到0.137MPa，减小了11.0%；基面层间纵向最大剪应力由0.11MPa增大为0.12MPa，增大了9%；面层间横向最大剪应力由0.149 MPa减小为0.132 MPa，减小了11.4%；基面层间横向最大剪应力由0.121MPa增大为0.131MPa，增大了8.3%。可见，沥青路面的粘弹性有利于面层间剪应力的松弛，但会增大基面层间剪应力。 　　上面层和下面层均为粘弹性材料，荷载作用下层间剪应力集中程度会随加载时间而减小，应力分布逐渐趋于均匀。结构体某一部分应力的减少必然伴随着其它部分应力的增长。剪应力在整个沥青混凝土面层内不断重分布，应力集中程度不断减小，当