动载作用下典型沥青路面剪应力特性分析Word下载.docx

动载作用下典型沥青路面剪应力特性分析Word下载.docx

《动载作用下典型沥青路面剪应力特性分析Word下载.docx》由会员分享，可在线阅读，更多相关《动载作用下典型沥青路面剪应力特性分析Word下载.docx（7页珍藏版）》请在冰豆网上搜索。

道路工程；

沥青路面；

三维有限元分析；

剪应力；

动荷载0　引言随着国民经济日益发展，公路运输对于高速、重载的要求与日俱增，高等级沥青路面在使用过程中暴露出寿命不足、早期病害严重等缺陷。

其中，由于路面抗剪切能力不足所产生的路面推移、拥包以及流动性或压密性车辙问题[1]，是沥青路面使用过程中的显著病害。

因此，研究沥青路面的剪应力特性，对合理解释路面损坏原因，保障路面使用寿命具有重要意义。

我国现行沥青路面设计规范[2]将车辆荷载视为均布竖向静态荷载，与实际行车荷载间存在较大差异，此外在控制和验算指标选取时，未对剪应力指标作出明确规定；

文献[3-4]采用路面结构有限元法模拟了剪应力在沥青路面结构中的分布情况;

文献[5]考虑荷载影响作用，研究了超载状态下路面剪应力的变化规律;

文献[6]针对半刚性基层沥青路面研究了最大剪切应力的大小和位置。

但以上研究都是基于静力学理论的“理想化”分析结果，与动载模式下的实际力学响应存在较大差异。

文献[7]分析了半正弦荷载作用下，温度及深度因素对沥青路面剪应力的影响规律，但其研究过程未考虑水平荷载的影响作用。

文献[8]对沥青路面的动态剪切响应进行了研究，得出了刹车反复作用下沥青路面的力学行为特征，但在模型构建时将路面各结构层视为完全连续，忽略了层间接触状态对路面抗剪性能的影响。

鉴于此，本研究以商业通用有限元软件ABAQUS为平台，研究了半刚性（S1）、倒装式（S2）和组合式（S3）3种典型沥青路面内部动态剪应力的空间分布特性，并通过动、静剪应力对比，明确了路面动力行为分析的重要性。

在此基础上，进一步分析了路面结构参数及车辆荷载参数对剪应力的影响规律，以期为探究路面损坏机理提供参考。

1　沥青路面结构动力有限元模型1.1　路面结构基本模型参考文献调研结论并结合试算验证，以计算结果收敛稳定为依据，拟定路面结构的平面尺寸为6m×

6m，路基计算深度取6m。

路面结构和其荷载形式具有对称性，为提高运算效率，取1/4模型进行分析。

边界条件为：

路基底面（z方向）固定、各结构层对称面施加对称约束、侧面（x和y方向）分别施加水平约束[9-10]，其中x和y方向分别表示道路横向和纵向，有限元模型中各结构层采用八节点六面体完全积分单元C3D8模拟。

路面结构有限元基本模型如图1所示。

图1　路面结构有限元模型（1/4模型）

Fig.1　Finiteelementmodelofpavementstructure（1/4model）1.2　荷载参数考虑到行车荷载模式的复杂性，以固定位置（车胎与路表接触面）的冲击荷载作为路面简化荷载场，将轮印区域视为当量圆，并假设接触压力均匀分布。

车辆静止时，轴载采用恒定静载模式；

对于动态荷载，其荷载集度变化参考YangH.Huang在KENLAYER程序中的简易方法[11]，荷载强度随时间呈半正弦函数变化（图2），具体荷载模型见式

（1）。

（1）式中，q为荷载强度峰值，标准轴载取0.7MPa；

T为动荷载作用周期；

R为轮胎与地面的接地当量圆半径，标准轴载取0.1065m。

图2　动荷载强度变化

Fig.2　Variationofdynamicloadstrength对于动态荷载模式，其荷载作用时间取决于行车速度v和轮胎接触面半径R，根据KENLAYER计算程序假设[11]，荷载作用周期T的计算见式

（2）：

（2）式中符号同上。

在具体计算中，为降低网格划分的难度，并提高网格的质量，将轮印面积简化为与当量圆面积相等的矩形面积[11]（图3）。

以单轴双轮组100kN荷载作为标准轴载，根据图3所示矩形的长宽比，每个矩形轮印尺寸为：

长a=22.78cm（行车方向），宽b=15.68cm，如图4所示。

图3　轮印面积简化图形

Fig.3　Simplifiedgraphicsoftireimpressionarea图4　标准轴载荷载分布图示（单位：

cm）

Fig.4　Diagramofstandardaxleloaddistribution（unit:

cm）1.3　路面结构组合材料参数参考国内外沥青路面结构的主导形式[12]，以半刚性路面（S1）、倒装式路面（S2）和组合式路面（S3）3类典型沥青路面结构作为研究对象开展动力分析，如表1所示。

表1　路面结构组合Tab.1　Pavementstructurecombinations结构层位路面结构S1（半刚性路面）路面结构S2（倒装式路面）路面结构S3（组合式路面）上面层SMA-13SMA-13SMA-13中面层AC-20CAC-20CAC-20C下面层AC-25CAC-25CAC-25C基层水泥稳定碎石级配碎石沥青稳定碎石底基层水泥稳定碎石水泥稳定碎石水泥稳定碎石路基土基对于路面材料的动态模量而言，相同材料的模量值在不同荷载、温度、湿度等因素影响下的变化范围较大，本次计算参考沥青路面材料动态力学参数的取值范围[12]，以车速80km/h（作用时间0.058s）为标准条件，从中选取代表值作为模型计算参数。

此外，在实际路面结构中，由于工序间隔、施工污染及各结构层材料属性差异等因素，沥青路面各结构层的层间界面之间很难做到完全黏结，即通常处于部分连续状态，因此采用以层间摩阻系数μ表征的接触模型更有助于获取路面结构内部的真实应力传递规律[13]。

具体计算参数设置见表2。

表2　路面结构计算参数Tab.2　Calculationparametersofpavementstructure结构层位材料类型模量/MPa泊松比密度/（kg·

m-3）层间状态[14]上面层中面层下面层基层底基层路基SMA-13110000.252400AC-20140000.252400AC-25130000.252400水泥稳定碎石150000.202200级配碎石3500.352000沥青稳定碎石100000.252400水泥稳定碎石130000.202100土基1000.351800μ=0.7μ=0.7μ=0.5μ=0.5Tie绑定注：

表中μ值分别对应不同结构层之间的摩阻系数，其中上面层与中面层间、中面层与下面层间，μ取0.7；

下面层与基层间、基层与底基层间，μ取0.5；

底基层与路基之间为绑定连接（Tie）。

分析过程考虑阻尼影响作用，首先通过模态分析，获取结构模型在主振方向（z方向）的振型及自振频率，然后将路面结构阻尼比ζ设定为0.05[15]，通过第一、二阶自振频率ω1，ω2与结构阻尼ζ确定相应的Rayleigh阻尼分项系数α=5.159和β=3.622×

10-4，如式（3）所示：

（3）2　剪应力特性分析2.1　剪应力空间分布特性计算中车辆轴重取标准轴载BZZ-100，行车速度80km/h（对应的荷载作用周期为0.058s），空间坐标（图5）以轮轴中心为原点O。

道路横向为x轴，行车方向（纵向）为y轴，厚度方向为z轴。

考虑到面层内部最大剪应力的方向变化复杂，不宜用直观方式表示，因此难以准确评价剪应力所对应的路面破坏模式，本文取水平纵向（y方向）剪应力分析计算路面剪切状态，对应评价沥青路面由于纵向流动剪切引起的推移、拥包及车辙病害。

图5　空间坐标系

Fig.5　Spacecoordinatesystem图7　剪应力三维分布图

Fig.7　Three-dimensionaldistributionofshearstresses图6为路面内部剪应力在内侧轮印中心位置处随道路深度的变化情况。

3种路面结构的剪应力最大值均出现在上面层与中面层之间（40mm深度），中面层顶面位置。

其中，S2结构的剪应力最大值是S1和S3结构的2倍以上，因此更容易出现推移、拥包等破坏现象。

此外，由于S2结构采用级配碎石作为基层材料，与面层沥青混合料材料性能差异较大，剪应力在面层与基层之间（下面层底面与基层顶面间）发生急剧突变。

图6　剪切应力沿路面厚度的分布

Fig.6　Distributionofshearstressesalongpavementthickness图7为剪应力在沿路面厚度方向的最不利位置——中面层顶面的分布情况，剪应力主要集中在车轮荷载作用处，轮印区域的剪应力较轮隙处存在突变，剪应力大小随着与荷载作用点的纵、横向距离增加而急剧减小。

图8为剪应力沿道路纵向及横向的分布情况，通过对比发现，剪应力的最大值并非出现在轮印中心位置，往往出现在轮印边缘位置处。

3种路面结构的剪应力大小关系为：

S2结构的剪应力最大，而S1与S3结构的剪应力较为接近，S3略大于S1结构。

图8　中面层顶部的剪应力分布情况

Fig.8　Distributionofshearstressesattopofmid-surface2.2　动、静剪应力对比基于前文建立的沥青路面结构3D有限元模型，针对3种典型沥青路面开展动、静剪应力对比分析，如图9所示。

其中，柱状图分别对应不同荷载下的静力与动力响应量，图中带数据标记的折线表示路面结构动、静剪应力比值随轴重的变化。

图9　路面动、静剪应力响应对比

Fig.9　Comparisonbetweendynamicandstaticshearstressresponsesinpavement图10　模量对剪应力的影响

Fig.10　Influenceofmodulusonshearstress由图9可以看出，路面结构的动、静剪应力均随轴重的增加逐渐增大，但其大小差异明显：

相比于静态荷载，动载作用下的路面结构在力学性质上表现出模量值增大，其剪应力的动力响应明显大于静力响应。

若以动力修正系数μd（μd=Rd/Rs）表示各力学指标动力响应Rd与静力响应Rs的定量比值关系，剪应力的动、静力响应比值范围为1.2～2.3，其中S1与S3接近，比值范围为1.2～1.7；

S2为1.6～2.3。

在相同路面结构中，动、静剪应力之比随轴重的增加呈逐渐增大的趋势，表明随着轴重水平的提高，动态与静态荷载模式所引起的剪应力差异有所增大。

沥青路面动、静剪应力的差异性表明，为了更加有效地抑制路面早期剪切病害，并对其使用性能的及时修复做出合理决策，应采用反映实际交通荷载下路面结构力学行为的动力分析方法，开展路面分析与设计工作。

2.3　路面结构参数对剪应力的影响2.3.1　材料模量对动力响应指标的影响参考典型路面结构层的经验厚度（面层18cm+基层30cm+底基层20cm），结合各结构层材料动态模量值的变化范围，以面层模量（4水平）与基层模量（5水平）为变化参数（如表3所示），沥青层剪应力随路面材料模量的变化规律如图10所示。

表3　材料模量组合Tab.3　Combinationofmaterialmoduli组合因素路面结构类型半刚性基层路面倒装式路面组合式路面面层模量/MPa2000～200004水平（2000,8000,14000,20000）基层模量/MPa3000～360005水平（3000,6000,12000,24000,36000）200～6005水平（200,300,400,500,600）5水平（与面层模量的比值:

0.5,0.75,1.0,1.25,1.5）注：

组合式路面的面层与基层均为沥青类结构层，其基层模量水平的选取通过基层与面层的模量比值进行控制。

由图10可知3种路面具有相似的变化规律：

剪应力随基层模量增加逐渐减小，随面层模量的增大逐渐增大。

比较不同面层或基层模量水平下剪应力值的变化幅度发现：

随面层模量水平提高，S1与S3中剪应力随基层模量增大而减小的幅度增大，而S2中的变化幅度有所减小；

随基层模量水平提高，S1与S3中剪应力随面层模量增加而增大的幅度有所减小，而S2中的变化幅度有所增大。

这表明，在S1和S3中，面层模量水平越高，通过提高基层模量减小路面剪应力的效果越明显，而基层模量水平越低，通过降低面层模量减小路面剪应力的效果越明显；

与之相反，在S2中，面层模量水平越低，提高基层模量对减小剪应力的效果越明显，基层模量水平越高，降低面层模量对减小剪应力的效果越明显。

2.3.2　基层厚度对动力响应指标的影响在沥青路面结构设计过程中，结构厚度是一项重要的设计参数。

考虑到沥青面层材料要求严、造价高的特点，分析中根据经验采用典型面层厚度，而将基层厚度作为路面结构厚度变量。

以半刚性路面为代表性结构，基层厚度分别取10，20，30，40cm及50cm，剪应力随基层厚度的变化情况如图11所示。

由图11可知，剪应力随基层厚度增加逐渐减小，而随着基层厚度的增大，其递减速率呈逐渐减小的趋势。

以基层模量取8000MPa为例，当基层厚度由10cm增加至20cm时，沥青层剪应力减小了15.5%，随着基层厚度水平逐渐升高，当其厚度由40cm增至50cm时，剪应力仅减小了5.0%。

这表明，当基层厚度较大（>

40cm）之后，通过提高基层厚度来改善路面抗剪性能的效果明显减弱。

图11　基层厚度对剪应力的影响

Fig.11　Influenceofbasethicknessonshearstress2.4　车辆荷载参数对剪应力的影响2.4.1　车速、轴重及轴距对剪应力的影响以半刚性路面（S1）为标准模型，不同行车速度、轴重及轴距条件下，剪应力的时程曲线分别如图12（a）～（c）所示。

图12　车速、轴重与轴距对剪应力的影响

Fig.12　Influenceofvehiclespeed,axleloadandwheelbaseonshearstress由图12（a）可知，车速对剪切应力有较大影响，随车速提高，剪应力逐渐增大，车速120km/h时的剪应力较20km/h增加了20.46%。

由此可推断，在车辆运行过程中对超速现象进行控制，是避免沥青路面在车辆荷载作用下产生剪切类病害的有效途径。

由图12（b）可知，随着荷载轴重的增加，车辆作用与路面上的动荷载效应逐渐增大，在整个加载时程内，路面剪应力有显著提高，当轴重超载160%（荷载260kN）时，剪应力峰值增加90.2%。

由图12（c）可知，以1.1m轴距为分界点，剪应力随轴距增加呈先增大后减小的趋势，轴距取1.1m时，剪应力出现最大值，对沥青层抗剪最为不利。

当轴距≥1.4m之后，剪应力响应与单轴标准轴载基本相等，且轴距变化对剪应力的影响作用较小。

2.4.2　考虑水平荷载效应的剪应力分析在车辆与路面相互作用过程中，不仅产生竖向荷载，还由于轮胎与路面间的滚动摩阻而产生水平切向作用力[16]。

车胎与路面间的摩阻作用通常较小，水平荷载可忽略不计。

但在车辆频繁启动、制动以及陡坡急弯路段，水平荷载效应趋于明显。

本节重点针对车辆制动过程中水平荷载所引起的路面动力响应进行相关数值计算，车辆启动时的情况可通过相同方法进行等效分析。

为简化计算过程并体现代表性，计算中对车辆制动过程及其荷载作如下假设[8]：

（1）车辆在制动过程作匀减速直线运动，车速在制动距离内迅速衰减为0；

（2）车辆对路面施加竖向与水平综合荷载作用，水平及竖向荷载在荷载接触面上均呈均匀分布，且作用周期相同；

（3）水平荷载仅考虑车胎与路面间的滑动摩擦作用，其荷载峰值可通过牛顿第二定律确定。

为更加全面地评价水平荷载对于路面结构剪切特性的影响，以半刚性基层路面结构为代表，采用单轴双轮组荷载，分别考虑轴重取100kN（标准轴载）、180kN（超载），车速取20km/h（低速）、80km/h（高速），对应的制动加速度取10m/s2（紧急制动）、5m/s2（普通制动）、2m/s2（轻微制动）以及0m/s2（正常行驶），以对比不同工况下的动态剪切应力。

其中水平荷载峰值结合牛顿第二定律以及匀减速直线运动方程，由轴重、行车速度、加速度以及荷载作用面积确定，计算结果见表4。

表4　不同工况计算方案及结果Tab.4　Calculationschemesandresultsunderdifferentconditions轴载/kN车速/（km·

h-1）制动加速度/（m·

s-2）水平荷载峰值/kPa剪应力/kPa1002000144.22140.0166.85350.0200.710699.9257.11008000169.12140.0192.55350.0227.510699.9285.91802000198.62170.1228.25425.3272.510850.5346.61808000246.82170.1272.25425.3318.110850.5394.6图13直观地反映出路面结构在车辆荷载不同制动情况下的剪切动响应状态。

由表4与图13可知，在不同轴重及车速条件下，制动加速度对面层剪应力的影响均表现为：

随加速度a（m/s2）的增加0（正常行驶）→2（轻微制动）→5（普通制动）→10（紧急制动），峰值剪应力基本呈线性关系增大。

车速越高，相同制动加速度所引起的剪应力越大，具体而言：

对于紧急制动（a=10m/s2），车速80km/h时产生的剪应力较20km/h提高了11.2%（100kN轴载）和13.8%（180kN轴载）。

同样，轴载水平越高，相同制动加速度引起的剪应力越大，紧急制动（a=10m/s2）时，轴载180kN产生的剪应力较100kN增加了34.8%（20km/h）和38.0%（80km/h）。

车速的提高与轴重的增加，都会使路面在车辆制动过程中的剪应力增大。

轴载100kN、以20km/h车速匀速行驶的车辆荷载产生的剪应力为144.2kPa；

而轴载180kN、车速80km/h、制动加速度10m/s时产生的路面剪应力为394.6kPa。

后者为前者的2.7倍，即当重载车辆在高速运行过程中面临紧急制动时，对于路面剪切受力极为不利。

因此，为避免路面出现剪切类早期病害，应在控制车辆超速、超载的基础上，提高道路设计的合理性，保障车辆的制动距离，并加强频繁启动、制动路段的材料抗剪性能。

图13　不同工况的剪应力对比

Fig.13　Comparisonofshearstressesunderdifferentconditions3　结论

（1）沿路面深度方向，3种路面结构的剪应力最大值均出现在上面层与中面层之间（40mm深度），即中面层顶面位置；

剪应力在水平面的分布主要集中在车轮荷载作用处，其响应量随着与荷载作用点的纵、横向距离增加而急剧减小，最大值出现在轮印边缘位置处。

3种路面的剪应力大小为：

倒装式路面（S2）的剪应力最大，而半刚性路面（S1）与组合式路面（S3）的剪应力较为接近，S3略大于S1。

（2）相比于静态荷载，动荷载作用下的剪应力值更大，并且轴重水平越高，动、静剪应力的差异越显著。

剪应力随基层模量增加而减小，随面层模量的增大而增大。

剪应力随基层厚度增加逐渐减小，而随着基层厚度的增大，其递减速率呈逐渐减小的趋势，当基层厚度较大（≥40cm）之后，通过提高基层厚度来改善路面抗剪性能的效果明显减弱。

（3）随着车速和轴重的提高，路面剪应力均逐渐增大；

剪应力随轴距增加呈先增大后减小的趋势，在轴距取1.1m时，剪应力出现最大值，而当轴距≥1.4m之后，剪应力响应与单轴标准轴载基本相等，且轴距变化对剪应力的影响作用较小。

（4）随着行车荷载加速度的增加，路面剪应力基本呈线性关系增大。

车速越高，相同制动加速度所引起的剪应力越大；

轴载水平越高，相同制动加速度引起的剪应力越大。

重载车辆在高速运行过程中的紧急、频繁制动，是造成路面剪切损坏的重要原因。

References:

[1]　沈金安，李福普，陈景.高速公路沥青路面早期损坏分析与防治对策[M].北京：

人民交通出版社，2004.SHENJin-an,LIFu-pu,CHENJing.AnalysisandPreventiveTechniquesofPrematureDamageofAsphaltPavementinExpressway[M].Beijing:

ChinaCommunicationsPress,2004.[2]　JTGD50—2006，公路沥青路面设计规范[S].JTGD40—2006，SpecificationsforDesignofHighwayAsphaltPavement[S].[3]　聂忆华，张起森.高等级公路沥青路面剪应力分析与应用[J].中南大学学报：

自然科学版，2007，38（6）：

1232-1238.NIEYi-hua,ZHANGQi-sen.AnalysisandApplicationofShearStressinHigh-typeAsphaltPavement[J].JournalofCentralSouthUniversity:

ScienceandTechnologyEditon,2007,38（6）:

1232-1238.[4]　聂忆华，张起森，徐阳，等.高等级公路沥青路面剪应力分布研究[J].重庆建筑大学学报，2007，29（5）：

85-90.NIEYi-hua,ZHANGQi-sen,XUYang,etal.DistributionRegularitiesandInfluencingFactorsofShearStressinHighGradeAsphaltPavementLayers[J].JournalofChongqingJianzhuUniversity,2007,29（5）:

85-90.[5]　张静，宁晓霞，王耀军.超重车辆荷载作用下沥青混凝土路面剪应力分析[J].公路，2013，1

（1）：

165-170.ZHANGJing,NINGXiao-xia,WANGYao-jun.SShearStressAnalysisofAsphaltConcretePavementUnderSuper-heavyTraffic[J].Highway，2013，1

（1）:

165-170.[6]　陈宝，吴德军，刘崭.沥青路面结构的最大剪应力[J].长安大学学报：

自然科学版，2010，30（6）：

20-23.CHENBao,WUDe-jun,LIUZhan.