公路交通科技  2016, Vol. 31 Issue (4): 1-6

扩展功能

文章信息

颜可珍, 游凌云, 葛冬冬, 庾付磊
YAN Ke-zhen, YOU Ling-yun, GE Dong-dong, YU Fu-lei
横观各向同性沥青路面结构力学行为分析
Analysis of Structural Mechanical Behavior of Transverse Isotropic Asphalt Pavement
公路交通科技, 2016, Vol. 31 (4): 1-6
Journal of Highway and Transportation Research and Denelopment, 2016, Vol. 31 (4): 1-6
10.3969/j.issn.1002-0268.2016.04.001

文章历史

收稿日期: 2015-1-3
横观各向同性沥青路面结构力学行为分析
颜可珍1, 游凌云1, 葛冬冬1, 庾付磊2    
1. 湖南大学 土木工程学院, 湖南 长沙 410082;
2. 郑州市政工程勘测设计研究院, 河南 郑州 450052
摘要: 考虑到沥青混合料、碎石基层等材料横观各向同性特征对沥青路面结构力学行为的影响。基于现行大多数路面设计理论的综合分析,运用有限单元分析方法在弹性假设条件下引入材料水平模量与竖向模量的比值,建立三维有限元模型。提出更加符合现今公路路面受力情况的非均布荷载作用形式、材料参数、边界条件,对比分析非均布荷载作用下各向同性与不同横观各向同性条件下的各结构层应力、应变及路表弯沉的差异。结果表明,将路面结构材料考虑为各向同性,明显低估了路面结构永久变形和疲劳开裂的危害,在路面结构设计中应适当考虑材料的横观各向同性特征。
关键词: 道路工程     沥青路面     有限元     横观各向同性     非均布荷载    
Analysis of Structural Mechanical Behavior of Transverse Isotropic Asphalt Pavement
YAN Ke-zhen1, YOU Ling-yun1, GE Dong-dong1, YU Fu-lei2     
1. College of Civil Engineering, Hunan University, Changsha Hunan 410082, China;
2. Zhengzhou Municipal Engineering Survey and Design Institute, Zhengzhou Henan 450052, China
Abstract: Considering the effect of transverse isotropic characteristics of asphalt mixture and macadam on the structural mechanical behavior of asphalt pavement, based on the comprehensive analysis with most current pavement design theories, by introducing the ratio of material horizontal modulus to vertical modulus, a 3D finite element model is established using finite element numerical analysis under elastic assumption. The irregularly distributed load, material parameters, and boundary conditions which are more in line with the loads on current highway pavement are put forward. The differences of stress, strain in each structure layer and deflection of pavement surface under isotropy condition with non uniform load or under different transverse isotropic conditions are comparatively analysed. The result shows that pavement structural material based on isotropy obviously underestimates the permanent deformation and fatigue failure of pavement structure, the design of pavement structure should be appropriate to take the transversely isotropic characteristics of material into consideration.
Key words: road engineering     asphalt pavement     finite element     transverse isotropy     irregularly distributed load    
0 引言

现行大多数的路面设计理论是采用各向同性的线弹性体系,路面结构材料一般采用弹性或黏弹性模型,荷载一般采用均布静载模型。路面实际使用寿命大都明显低于设计使用寿命,应从路面结构设计和材料参数的角度重新考虑新的路面设计方法。实际路面结构材料如沥青混合料、碎石基层材料等由于其结构组成及压实工艺等不同而具有明显的各向异性特征[1, 2]。路面结构层材料的各向异性与各向同性具有明显不同的力学性质,基于各向同性的沥青路面设计将低估路面结构的剪应力及拉应力,导致路面出现早期的永久变形和疲劳开裂[1, 3, 4, 5, 6]。另外,实际的车辆轮胎接地压力呈现明显的非均布特征,而轮胎的接地压力大小与分布对沥青路面结构力学行为影响明显。现行均布静载模型无法准确反映出车辆荷载接地压力的非均布特性[7]。目前的沥青路面破坏机理分析中对这些影响因素的分析有所欠缺,因此很有必要针对沥青路面在非均布荷载作用下结构层材料横观各向同性特征深入研究其力学行为。

横观各向同性作为一种特殊的各向异性表征,能较好地表征路基路面结构材料的各向异性导致的路面结构行为差异[1, 8, 9, 10, 11, 12]。本文通过有限元软件建立非均布荷载作用下沥青路面结构模型,总结沥青路面结构响应规律;分析结构层材料横观各向同性对路面结构力学行为的影响,比较材料各向同性假设与横观各向同性假设对沥青路面结构力学行为的影响。

1 横观各向同性路面有限元模型建立 1.1 荷载模型

实际车轮荷载作用于路面不是圆形均布,而是非均匀分布。实际轮胎与路面的接触形状和轮胎的花纹会影响接触压力分布,一般情况下接触面上的压力分布亦是不均布。在一定荷载作用下,轮胎接地面积及每个轮胎的接触压力均可通过测试得到[13]。根据实际分析,本研究采用双组轮纵向花纹荷载模型(每个轮胎受力37.8 kN),荷载作用简化如图 1所示。

图 1 轮胎接地面积及压力分布简化图 Fig. 1 Simplified diagram of tire contact area and pressure distribution
1.2 有限元模型

本研究建立的有限元模型尺寸为5 m×5 m×10 m(xyz),x为道路横向,y为道路纵向,z为道路深度方向。非均布荷载作用于路面模型的中心区域,底面为固定面,无xyz向位移。网格划分在荷载受力区域并在面层及基层加密,其他区域作单精度处理。假定层间完全连续,单元类型为八节点线性六面体单元。模型边界条件和网格如图 2所示。

图 2 横观各向同性路面的三维有限元模型 Fig. 2 Three-dimensional finite element model of transverse isotropic pavement
1.3 路面结构及材料参数

有限元计算分析过程中采用典型的三层结构组合的路面结构形式,从上往下依次为沥青混合料面层、碎石材料基层、土路基,假定土路基为弹性半空间,路面各结构层间连续,分别考虑面层及基层材料的横观各向同性,并确定横观各向同性假设条件下的路面结构层材料参数[1, 6, 12, 14]

(1)假定沥青面层为横观各向同性,碎石基层和土路基层为各向同性。沥青面层竖向模量Ev=1 400 MPa,水平向与竖向泊松比相等(μv=μh=0.25),水平模量与竖向模量比nM分别取0.2,0.5,0.7,1.0。

(2)假定碎石基层为横观各向同性,沥青面层和土路基层为各向同性。碎石基层竖向模量Ev为500 MPa,水平向与竖向泊松比相等(μv=μh=0.25),水平模量与竖向模量比nJ分别取0.17,0.21,0.5和1.0,如表 1所示。

表 1 路面结构材料参数 Tab. 1 Material parameters of pavement structure
结构层 材料特性 弹性模量/MPa 泊松比 密度/(kg·m-3) 厚度/cm
沥青面层 各向同性 1 400 0.25 2 400 18
横观各向同性
碎石基层 各向同性 500 0.252 100 35
横观各向同性
土路基 各向同性 50 0.35 1 800
2 横观各向同性面层路面结构力学分析 2.1 横观各向同性面层路表弯沉分析

路表弯沉是从路面总体结构与宏观性能方面分析控制路面结构使用性能的关键指标。图 3是非均布荷载作用下考虑面层材料横观各向同性特征时路表弯沉的影响曲线。可以看出,荷载作用中心附近弯沉值最大,且距离中心位置越远弯沉值越小,随着面层材料水平模量与竖向模量比值nM的增大,路表弯沉逐渐减小。由图 3(a)可知,面层材料水平模量和竖向模量比值nM=0.2时路表弯沉值最大,峰值约790 μm;nM=1.0(面层材料各向同性)时路表弯沉值最小,峰值约为700 μm。由图 3(b)可知,轮隙中心弯沉与水平模量和竖向模量比值nM呈近似线性关系,且随nM增大弯沉逐渐减小。对比各测点弯沉,可以知道面层材料横观各向同性对路面变形影响相对较小。当采用弯沉作为控制指标进行路面结构设计及分析时,可以忽略面层材料横观各向同性的影响。

图 3 横观各向同性面层对路表弯沉的影响 Fig. 3 Influence of transverse isotropic surface layer on pavement surface deflection
2.2 横观各向同性面层路面结构层应力及应变分析

面层层底应力及应变是控制疲劳开裂和路面服务寿命的重要指标。图 4是非均布荷载作用下考虑面层材料横观各向同性特征时,面层层底应力及应变随X的变化曲线。由图 4(a)可知,随着距荷载中心位置距离的增大,面层层底拉应力先增大后减小;随着面层材料水平模量与竖向模量比nM的增大,面层层底拉应力增大。在距离荷载作用中心约x=0.17 m处,面层层底拉应力出现峰值(nM=1.0时峰值约为100 kPa,nM=0.2时峰值约为-80 kPa)。在荷载作用边缘x=0.5 m处,面层层底拉应力均开始趋近于零,且各曲线趋于重合。由图 4(b)可知,随着距荷载中心位置距离的增大面层层底拉应变先增大后减小,在荷载作用中心附近面层层底拉应变受面层材料水平模量与竖向模量比值nM的影响较小,而在荷载作用边缘约x=0.5 m处受面层材料水平模量与竖向模量比值nM的影响较大,随nM的增大而增大。因为面层拉应变是控制面层疲劳寿命的重要指标之一,从分析结果可以知道考虑面层材料横观各向同性时,在相同荷载作用位置路面服务寿命越短,面层越易出现疲劳开裂。因此,应在路面结构设计时考虑面层材料的横观各向同性特征。

图 4 考虑面层材料面层底部应力及应变 Fig. 4 Stresses and strains of surface bottom considering surface course material

基层层底应力及应变是沥青路面结构的重要设计指标,图 5是非均布荷载作用下考虑面层材料横观各向同性特征时基层层底应力及应变随X的变化曲线。可以看出,面层材料横观各向同性对基层层底应力及应变影响规律基本一致。由图 5可知,随着距荷载中心位置距离的增大,基层层底拉应力及拉应变均明显减小,而基层层底拉应力及拉应变受面层材料水平模量与竖向模量比值nM的影响相对较小,基本表现为nM越大,基层层底拉应力及拉应变越小。

图 5 考虑面层材料基层底部应力及应变 Fig. 5 Stresses and strains of base bottom considering surface course material

土路基顶压应变是路基永久变形和路面车辙的控制指标,图 6是非均布荷载作用下考虑面层材料横观各向同性特征时土路基顶压应力随X的变化曲线。可以看出,随着距荷载中心位置距离的增大,土路基顶压应变不断减小最终趋于零;随着面层材料水平模量与竖向模量之比值nM的增加,土路基顶压应变不断减小。在荷载作用中心位置,土路基顶压应变随水平模量与竖向模量nM的变化最为显著,面层材料水平模量与竖向模量之比nM=0.2时,土路基顶压应变峰值达到最大值641.4 μεnM=1.0(面层材料各向同性)时,土路基顶压应变峰值最小值508.1 με。考虑面层材料各向同性(nM=1.0)时,路面车辙将比实际偏小,设计偏不安全。当以车辙作为路面结构设计的控制指标时,应考虑面层材料的横观各向同性特征。

图 6 考虑面层材料土基顶部压应变 Fig. 6 Compressive strains of soil base top considering surface course material
3 横观各向同性基层路面结构力学分析 3.1 横观各向同性基层弯沉分析

图 7是非均布荷载作用下,碎石基层材料横观各向同性特征时路表弯沉影响曲线。由图 7(a)可知,随着距离荷载中心距离的增加,路表弯沉逐渐减小并最终趋于重合;随着基层材料水平模量与竖向模量比值nJ增大,路表弯沉减小。基层材料水平模量与竖向模量比值nJ=0.17时弯沉值最大,nJ=1.0(基层材料各向同性)时弯沉最小。在荷载作用中心附近,弯沉达到峰值,且越靠近作用中心,基层材料横观各向同性对路表弯沉影响越大。由图 7(b)可知,基层材料水平模量与竖向模量比值nJ越大,轮隙中心弯沉值越小,且随着nJ增大弯沉减小幅度降低。对比各测点弯沉可以知道,基层材料横观各向同性对路面变形影响较大,当采用弯沉作为路面结构设计控制指标时,应考虑基层材料的横观各向同性特征。

图 7 考虑基层材料横观各向同性基层对路表弯沉的影响 Fig. 7 Influence of transverse isotropic base course on pavement surface deflection Consideing base Course material
3.2 横观各向同性基层路面结构层应力及应变分析

图 8是非均布荷载作用下,考虑碎石基层材料横观各向同性特征时面层层底应力及应变随X的变化曲线。可以看出,基层材料横观各向同性对面层层底应力及应变影响规律基本一致。面层底部拉应力及拉应变随基层材料水平模量与竖向模量比值nJ的增加而减小,nJ=1.0时(基层材料各向同性)面层底部的拉应力及应变最小。面层底部应力及应变随X的增大,先增大后减小。在荷载作用中心附近约X=0.17 m处,面层底部拉应力及拉应变几乎同时出现峰值。当基层材料水平模量与竖向模量比值nJ=0.17时,面层底部拉应力及拉应变峰值最大(分别为240.5 kPa、194.3 με);nJ=1.0时,其峰值相对最小(分别为93.5 kPa、109.4 με)。面层拉应变是决定面层疲劳寿命的重要指标之一。从分析结果可以看出,考虑基层材料横观各向同性时,相同荷载作用位置路面结构疲劳寿命越短,面层越易出现疲劳开裂。因此,应在路面结构设计中考虑基层材料的横观各向同性特征。

图 8 考虑基层材料面层底部应力及应变 Fig. 8 Stresses and strains of surface bottom considering base course material
图 9是非均布荷载作用下,考虑碎石基层材料横观各向同性时,基层层底应力及应变随X的变化曲线。可以看出,基层底部拉应力及拉应变均随着距荷载中心距离增大相对减小;而基层层底拉应力随着基层材料水平模量与竖向模量比值nJ的增大不断增大;基层层底拉应变随着基层材料水平模量与竖向模量比值nJ的增大不断减小。基层底部拉应力及拉应变均在荷载作用中心达到峰值,且在nJ=1.0(基层材料各向同性)时基层层底拉应力最大,拉应变最小(分别为105.4 kPa、179.3 με);在nJ=0.17时基层层底拉应力最小,拉应变最大(分别为50.4 kPa、610.3 με)。基层拉应变是控制路面结构使用性能的关键指标之一。分析结果可以看出,考虑基层材料横观各向同性时,相同荷载作用位置路面结构使用性能相对降低(基层拉应变过大会导致面层较明显的反射裂缝)。因此,路面结构设计中应当考虑基层材料的横观各向同性特征。

图 9 考虑基层材料基层底部应力及应变 Fig. 9 Stresses and strains of base bottom considering base course material
图 10是非均布荷载作用下考虑碎石基层材料横观各向同性特征时土路基顶压应力随X的变化曲线。可以看出,从荷载作用中心到边缘, 土路基顶压应变逐渐减小并趋于平稳。在荷载作用中心附近,土路基顶压应变的变化最为显著,在X=0 m位置土路基顶压应变达到峰值。土路基顶压应变随着基层材料水平模量与竖向模量比值nJ的增大而减小。当nJ=0.17时,土路基顶压应变达最大值1 019.7 μεnJ=1.0(基层材料各向同性)时,土路基顶压应变峰值达到最小值508.1 με。考虑基层材料各向同性(nJ=1.0)时,路面车辙将会比实际偏小,设计偏危险。当以车辙为路面结构设计控制指标时,应当考虑基层材料的横观各向同性特征。

图 10 考虑基层材料土基顶部压应变 Fig. 10 Compressive strains of soil base top considering base course material
4 结论

(1)考虑面层材料横观各向同性时,面层材料的水平模量与竖向模量比值对路表弯沉、基层层底应力及应变影响相对较小;对面层层底应力及应变、土基顶压应变影响相对较大,且越是靠近荷载作用中心这种影响越明显。因此,在路面结构设计中应考虑面层材料的横观各向同性特征。

(2)考虑基层材料横观各向同性时,基层材料的水平模量与竖向模量比值对路表弯沉、面层和基层层底应力及应变、土基顶压应变影响均较明显,且越是靠近荷载作用中心这种影响越明显。因此,在路面结构设计中应考虑基层材料的横观各向同性特征。

(3)文中计算结果表明,考虑面层及基层材料横观各向同性特征时,路面结构使用性能基本在材料水平模量与竖向模量比值n=1.0(材料各向同性)时最佳,n<1.0(材料各向异性)时相对降低。因此,现行路面设计理论中,将结构层材料考虑为各向同性明显低估了路面结构永久变形及疲劳开裂等危害,在未来路面结构设计中应适当地将结构层材料横观各向同性特征考虑在内。

参考文献
[1] WANG L, HOYOS L R, WANG J, et al. Anisotropic Properties of Asphalt Concrete: Characterization and Implications for Pavement Design and Analysis[J].
[2] 彭余华,习红军,李芳芳,等.无核密度仪在沥青路面密度测定中的应用[J].筑路机械与施工机械化,2008,25(6):47-50.PENG Yu-hua, XI Hong-jun, LI Fang-fang, et al. Application of Pavement Quality Indicator in Asphalt Pavement Density Measuring [J].Road Machinery &Construction Mechanization, 2008,25(6):47-50.
[3] TUTUMLER E, LITTLE D, KIM S H. Validated Model for Predicting Field Performance of Aggregate Base Courses[J].
[4] 周晓青,李宇峙,应荣华,等. 沥青混合料拉伸疲劳试验下疲劳损伤特性研究[J]. 重庆建筑大学学报,2005,27(5):47-51.ZHOU Xiao-qing, LI Yu-zhi, YING Rong-hua, et al. Research on the Fatigue and Damage Characteristics of Asphalt Mixture Based on Uni-axial Fatigue Tension Testing at Low Temperature[J]. Journal of Chongqing Jianzhu University, 2005,27(5):47-51.
[5] 栗振锋,王秉刚,郭向云,等.基于横观各向同性倒装式沥青路面结构分析[J].重庆建筑大学学报,2007,29(1):65-69.LI Zhen-feng, WANG Bing-gang, GUO Xiang-yun, et al. Analysis of Inverted Asphalt Pavement Structure Based on Cross-anisotropy [J]. Journal of Chongqing Jianzhu University, 2007, 29(1):65-69.
[6] BOULBIBANE M, COLLINS I F, WEICHERT D, et al. Shakedown Analysis of Anisotropic Asphalt Concrete Pavements with Clay Subgrade[J].
[7] SEIBI A C. Development of Constitutive Relations for Asphalt Concrete under High Rates of Loading[D]. University Park: Pennsylvania State University,1993.
[8] OH J, FERNANDO E G, LYTTON R L. Evaluation of Damage Potential for Pavements due to Overweight Truck Traffic [J].
[9] OH J, LYTTON R L, FERNANDO E G. Modeling of Pavement Response Using Nonlinear Cross-Anisotropy Approach [J].
[10] CHEN J, PAN T, HUANG X. Numerical Investigation into the Stiffness Anisotropy of Asphalt Concrete from a Microstructural Perspective [J].
[11] WANG H, ALQADI I L. Importance of Nonlinear Anisotropic Modeling of Granular Base for Predicting Maximum Viscoelastic Pavement Responses under Moving Vehicular Loading [J].
[12] KIM S H, LITTLE D N, MASAD E, et al. Estimation of Level of Anisotropy in Unbound Granular Layers Considering Aggregate Physical Properties [J].
[13] ELSEIFI M A, AL-QADI I L, YOO P J, et al. Viscoelastic Modeling and Field Validation of Flexible Pavements[J].
[14] MASAD E, TASHMAN L, SOMEDAVAN N, et al. Micromechanics-based Analysis of Stiffness Anisotropy in Asphalt Mixtures[J].