基于路面材料模量应力和应变依赖模型的沥青路面结构分析当量方法与流程

本发明主要涉及道路工程领域，特别涉及由于路面材料非线性特性所引起的沥青路面结构受力分析问题。
背景技术：
：荷载作用下路面结构的真实力学响应，一直是路面结构分析理论与设计方法研究中的基础理论问题。由于组成材料的多样性和复杂性，路面结构在车辆荷载和自然环境作用下的力学响应，多表现出较为明显的非线性特性。多年来，为了揭示路面结构服役性能的演变机理、建立更加符合实际路面结构和材料力学响应行为的设计理论，国内外研究人员一直在努力探索荷载作用下路面结构的真实力学响应，但由于路面材料的非线性特性，常不能得到满意的结果。现阶段，基于线弹性假设的路面结构分析方法，由于物理方程不能正确描述沥青路面材料的非线性本构关系，加上结构层模量试验方法的多样性和取值的不确定性，导致难以准确、客观地评价路面结构的真实受力状态，不利于指导新建路面设计和在役路面寿命评估。按照经典的线弹性力学理论，材料的模量是其固有属性，与环境、外力无关，一般为定值。然而，大量的室内外试验表明，路面材料的模量与应力、应变、强度等指标密切相关，会随着荷载水平的变化而变化，而且不同应力路径、不同试验条件下的模量数值差异较大，表现出十分明显的非线性特性，路面材料的模量并不是定值，而是一个与应力、应变或强度有关的函数表达式，即应力(应变)依赖模型，为了获得路面结构的真实力学响应，必须考虑材料非线性问题对路面力学计算所带来的影响。技术实现要素：本发明为基于路面材料模量应力(应变)依赖模型的沥青路面结构分析当量方法，根据路面材料的非线性特性，采用路面材料模量应力(应变)依赖模型表征沥青路面结构各层材料模量，以沥青路面结构各层受力最不利点位作为当量计算点，提出一种全新沥青路面结构力学分析方法，该方法考虑了材料非线性问题对路面力学计算所带来的影响，使沥青路面结构受力分析更趋合理和完善。基于路面材料模量应力和应变依赖模型的沥青路面结构分析当量方法，包括如下步骤：1)选择拟分析的沥青路面结构形式和材料类型，确定计算荷载类型、层间结合条件、各层厚度hi和泊松比μi；其中：参数i代表沥青路面结构层数量，为正整数；2)根据拟分析的沥青路面结构形式和材料类型，通过室内试验确定沥青混合料复模量应变依赖模型、半刚性基层和路基土回弹模量应力依赖模型，以该组模型作为沥青路面结构分析时各层的模量取值依据，该组模型采用如下步骤获得：(1)沥青混合料复模量应变依赖模型确定方法①.在不同应变水平、不同加载频率、不同试验温度下进行沥青混合料复模量试验；②.按照下述方法获得以某一频率为基准频率，不同应变水平下、基于温度参数的沥青混合料复模量主曲线；(i)在某应变水平下，采用boltzmann函数表征某一频率时温度与沥青混合料复模量对数的相关关系，(ii)以(i)中的boltzmann函数关系为基准，根据时温等效原理引入温度移位因子α，将该应变水平下、其他频率时温度与复模量的试验数据进行平移处理，得到一组新的平移后温度-复模量试验，(iii)采用式(1)对(ii)中的数据进行回归分析，可得到该应变水平下，以某一频率作为基准频率基于温度参数的沥青混合料复模量主曲线，函数表达式为式(1)，式中：t——温度，单位为℃；e——以10为底的沥青混合料复模量对数；a1、a2、x0、dx——回归参数，(iv)采用(i)、(ii)、(iii)中的方法，可得到其他应变水平下，以某一频率为基准频率、基于温度参数的沥青混合料复模量主曲线，③.将步骤②中不同应变水平下的平移后温度-复模量试验数据，整理成自变量为温度和应变水平、因变量为复模量对数的数据格式，采用式(2)对其进行回归分析，可得到基于温度参数的沥青混合料复模量应变依赖模型；式中：t——温度，单位为℃；e——以10为底的沥青混合料复模量对数；ε——应变，单位为1×102με；a、b、c、x0、dx——回归参数；(2)半刚性基层回弹模量应力依赖模型确定方法[1].根据拟分析的沥青路面结构形式和材料类型，确定半刚性基层的无机结合料类型、被稳定材料类型和试件养生龄期，通过调整无机结合料剂量和被稳定材料级配组成，得到不同强度水平的半刚性基层材料；[2].采用步骤[1]中不同强度水平的半刚性基层材料，分别进行压缩受力模式的强度及回弹模量试验、弯拉受力模式的强度及回弹模量试验；[3].根据受力模式，将步骤[2]中的试验数据整理成自变量为强度和应力水平，因变量为回弹模量的数据格式；[4].采用式(3)的二元二次模型对步骤[3]中压缩受力模式下的数据进行回归分析，采用式(4)的幂函数模型对步骤[3]中的弯拉受力模式下的数据进行回归分析，可得到以式(3)和式(4)表达的基于强度水平的半刚性基层回弹模量应力依赖模型；e＝f(σ,r)＝a·σ2+b·σ+c·r+d式(3)e＝f(σ,r)＝a·rc·(σ+1)b式(4)式中：e——压缩或弯拉回弹模量，单位为mpa；σ——应力水平，单位为mpa；r——压缩或弯拉强度，单位为mpa；a、b、c、d——回归参数；(3)路基土回弹模量应力依赖模型确定方法<1>.在不同围压(σ3)和偏应力(σd)下进行路基土动态三轴回弹模量试验；<2>.在步骤<1>的不同围压(σ3)下进行路基土三轴破坏强度试验；<3>.将步骤<1>和步骤<2>中的试验数据整理成自变量为偏应力和破坏强度，因变量为动态回弹模量的数据格式；<4>.采用式(5)对步骤<3>中的数据进行回归分析，可得到以式(5)表达的基于强度指标的路基土回弹模量应力依赖模型。式中：e——动态三轴回弹模量，单位为mpa；σd——三轴试验偏应力，单位为mpa；pa——大气压力，取0.10138mpa；r——三轴破坏强度，单位为mpa；a、b、k1、k2、k3——回归参数；3)根据拟分析的沥青路面结构形式和材料类型，给定一组各层材料模量初始值ei0，结合步骤1)中的厚度hi和泊松比μi，以沥青路面结构各层受力最不利点位作为当量计算点，采用弹性层状体系理论计算各当量计算点的mises等效应力和等效应变；4)将步骤3)中各沥青层当量计算点的mises等效应变代入沥青混合料复模量应变依赖模型式(2)、各半刚性基层当量计算点的mises等效应力代入半刚性基层回弹模量应力依赖模型式(3)或式(4)、路基层当量计算点的mises等效应力代入路基土回弹模量应力依赖模型式(5)，可得到一组新的各层材料模量值5)比较和当二者相等或相对误差达到最小值时，计算结束；否则，重复步骤3)和步骤4)进行迭代计算，直至满足收敛精度要求为止；6)采用步骤5)中满足收敛精度要求时的各层模量值作为各结构层的最终模量ei；7)采用步骤6)中各结构层的最终模量，根据弹性层状体系理论计算沥青路面结构内部任意一点的应力、应变和位移，从而得到基于路面材料模量应力(应变)依赖模型和当量方法的沥青路面结构分析结果。该方法从路面材料的非线性特性出发，采用路面材料模量应力(应变)依赖模型表征沥青路面结构各层模量，以沥青路面结构各层受力最不利点位作为当量计算点，通过迭代计算获得各结构层的最终模量，在此基础上根据弹性层状体系理论计算沥青路面结构内部任意一点的应力、应变和位移。由于该方法考虑了材料非线性问题对路面力学计算所带来的影响，使沥青路面结构受力分析更趋合理和完善。所述步骤2)中沥青混合料复模量应变依赖模型确定方法的应变水平通常为30、60、90、120、150με，加载频率通常为10、15、20、25、30、35、40hz，试验温度通常为0、5、10、15、20、25、30、35、40、45℃，所述基准频率为10hz。所述步骤2)中半刚性基层回弹模量应力依赖模型确定方法的无机结合料类型为水泥、石灰或粉煤灰，被稳定材料类型为碎石、砾石、砂砾、砂或土，试件养生龄期通常为7d、28d、90d或180d，强度水平通常为5种以上，压缩受力模式指轴向无侧限压缩试验，弯拉受力模式指三点或四点梁式弯拉试验，强度及回弹模量数值均为代表值。所述步骤2)中路基土回弹模量应力依赖模型确定方法的不同围压(σ3)通常为3-5种水平，压力范围0-0.5mpa；不同偏应力(σd)通常为5种水平，压力范围0-1.0mpa。所述步骤3)和步骤4)中的沥青路面结构各层受力最不利点位(当量计算点)是指：沥青层和半刚性基层为荷载中心线与各层层底的交点、路基层为荷载中心线与路基顶面以下40cm位置的交点。所述步骤5)和步骤6)中的收敛精度要求为各层第n次与第n-1次模量相等或相对误差最小。本发明从路面材料的非线性特性出发，以沥青路面结构各层受力最不利点位作为当量计算点，提出了一种基于路面材料模量应力(应变)依赖模型的全新沥青路面结构力学分析方法，使沥青路面结构受力分析更趋合理和完善。附图说明图1计算分析体系，图2频率10hz时温度与复模量的boltzmann函数关系(应变水平30με)，图3不同频率时温度与复模量试验数据平移结果(应变水平30με、10hz为基准曲线)，图4基于温度参数的沥青混合料复模量主曲线(应变水平30με、10hz为基准曲线)，图5全部应变水平下基于温度参数的沥青混合料复模量主曲线(10hz为基准频率)，图6z轴上各点最大主应力计算结果图，图7z轴上各点最大主应变计算结果图，图8z轴上各点竖向位移计算结果图。具体实施方式以某沥青路面结构为例，说明基于路面材料模量应力(应变)依赖模型的沥青路面结构分析当量方法，具体步骤如下：步骤1.选定拟分析的沥青路面结构形式、材料类型和计算分析体系。(1)选择如图1所示的半刚性基层长寿命沥青路面结构作为拟分析的路面结构。该结构共7层，路面厚度和材料自上而下分别为：4cmsbs-ac13沥青层、8cmhmac20沥青层、20cm水泥稳定碎石cbg25半刚性基层、20cm水泥稳定碎石cbg25半刚性基层、20cm水泥稳定土cs半刚性基层、20cm水泥稳定土cs半刚性基层、路基层。(2)sbs-ac13和水泥稳定土cs泊松比取0.25，hmac20和水泥稳定碎石cbg25泊松比取0.20，路基层泊松比取0.35。(3)计算荷载类型为单圆均布荷载，荷载大小p＝0.7mpa，荷载圆半径＝15cm。以单圆荷载中心点为坐标系原点o，水平向右为x轴正向，垂直纸面向外为y轴正向，竖直向下为z轴正向，建立三维坐标系，如图1所示。(4)路面各结构层的层间结合条件为完全连续。步骤2.通过试验确定沥青路面各层模量的应力(应变)依赖模型。(1)确定sbs-ac13沥青层复模量应变依赖模型①.开展沥青混合料sbs-ac13的复模量试验，应变水平为30、60、90、120、150με，加载频率为10、15、20、25、30、35、40hz，试验温度为0、5、10、15、20、25、30、35、40、45℃，所得到的沥青混合料复模量试验结果见表1。对表1中的复模量数据取对数，结果见表2。表1沥青混合料sbs-ac13复模量试验结果表2沥青混合料sbs-ac13复模量试验结果(复模量取对数)②.根据下述步骤获得应变水平为30με时，以10hz为基准频率、基于温度参数的沥青混合料复模量主曲线。(i)以表2中数据为基础，采用boltzmann函数，建立应变水平为30με、频率10hz时温度与沥青混合料复模量对数的关系，如图2所示。(ii)以频率10hz时温度与复模量的boltzmann函数关系为基准曲线，根据时温等效原理计算每个加载频率下的温度移位因子αi,j(结果见表3)，将不同频率时温度与复模量的试验数据进行平移，平移后的结果见图3和表4。表3不同频率下的温度移位因子(应变水平30με)表4不同频率下平移后温度与复模量试验数据(应变水平30με)(iii)采用式(1)对表3数据进行回归分析，可得到应变水平为30με时、基于温度参数的沥青混合料复模量主曲线(如图4所示)，主曲线表达式见式(6)。③.根据步骤②的方法可以获得应变水平为60、90、120、150με时，以10hz为基准频率、基于温度参数的沥青混合料复模量主曲线，如图5所示。主曲线中的各个回归参数数值见表5。表5不同应变水平下温度与复模量主曲线关系中的回归参数回归参数30με60με90με120με150μεa14.48544.48294.50134.49594.5656a22.37332.16771.98681.98651.7653x029.507130.969132.058031.256232.2689dx14.172314.007314.584914.400216.0970④.将步骤③中不同应变水平下温度与复模量试验数据的平移结果，整理成自变量为温度和应变水平、因变量为复模量对数的数据格式，见表6。表6不同应变水平下温度与复模量试验数据平移结果整理表⑤.采用式(2)对表5中数据进行回归分析，可得到沥青混合料sbs-ac13基于温度参数的沥青混合料复模量应变依赖模型，见式(7)。(2)确定hmac20沥青层复模量应变依赖模型采用步骤(1)中的方法，同理可确定hmac20沥青层复模量应变依赖模型，见式(8)。(3)确定水泥稳定碎石cbg25半刚性基层回弹模量应力依赖模型由图1所示的沥青路面结构形式和材料类型可知，该结构中半刚性基层的弯拉受力模式比压缩受力模式更为不利，因此选择弯拉受力模式确定水泥稳定碎石回弹模量的应力依赖模型。①.水泥稳定碎石cbg-25结合料类型为水泥、被稳定材料为碎石、试件养生龄期为90d。通过调整无机结合料剂量和被稳定材料级配组成，得到了弯拉受力模式下5种不同强度水平的材料。②.采用步骤①中的水泥稳定碎石cbg-25，进行弯拉受力模式的强度及回弹模量试验，结果见表6。③.采用式(4)对表7中弯拉受力模式下的数据进行回归分析，可得到式(9)。e＝11272.5×r0.84×(σ+1)-0.29式(9)式(9)即为水泥稳定碎石cbg-25不同受力模式下基于强度水平的无机结合料稳定材料回弹模量应力依赖模型。表7水泥稳定碎石cbg-25弯拉受力模式下强度及回弹模量试验结果(4)确定水泥稳定土cs半刚性基层回弹模量应力依赖模型采用步骤(3)中的方法，同理可确定水泥稳定土cs回弹模量应力依赖模型，见式(10)。e＝7951.2×r0.53×(σ+1)-0.36式(10)(5)确定路基土回弹模量应力依赖模型①.进行不同围压(σ3)和偏应力(σd)下路基土动态三轴回弹模量试验。围压水平(σ3)为3种：0.14mpa、028mpa、042mpa，偏应力(σd)水平为5种：0.2mpa、0.3mpa、0.4mpa、0.5mpa、0.6mpa，试验结果见表8。表8路基土动态三轴回弹模量试验结果②.采用步骤①中的围压水平，进行不同围压(σ3)下路基土三轴破坏强度试验，结果见表9。表9路基土三轴破坏强度试验结果σ3/mpa0.140.280.42r/mpa0.0760.1030.130③.将步骤①和步骤②中的试验数据整理成自变量为偏应力和破坏强度，因变量为动态回弹模量的数据格式，见表10。表10路基土动态三轴回弹模量及破坏强度试验结果汇总表④.采用式(5)对表10中的数据进行回归分析，各回归参数数值见表11，代入式(5)可得到以式(11)表达的某路基土基于强度指标的路基土回弹模量应力依赖模型。e＝f(σ,r)＝157.2274·(9.8639·σd+1)-1.3721(r)0.4189(11)表11路基土三轴试验模式下应力依赖模型的回归参数回归参数abk1k2k3相关系数r2数值3.363310.245984.04214-1.372111.702980.90(6)路面结构各层材料模量应力(应变)依赖模型汇总至此，可以得到拟分析路面结构各层材料模量的应力(应变)依赖模型，见表12。表12路面材料模量应力(应变)依赖模型步骤3.计算路面材料模量为初始值时，各当量计算点的mises等效应力和等效应变。(1)给定一组各层材料模量初始值结合步骤1可确定用于沥青路面结构分析的材料参数，见表13。表13沥青路面各层材料模量初始值(2)确定当量计算点位置以沥青路面结构各层受力最不利点位作为当量计算点，其中：沥青层和半刚性基层为荷载中心线z轴与各层层底的交点、路基层为荷载中心线z轴与路基顶面以下40cm位置的交点，即图1中的点a1至a7。在三维坐标系中，各当量计算点的坐标分别为(0,0,4)、(0,0,12)、(0,0,32)、(0,0,52)、(0,0,72)、(0,0,92)、(0,0,132)，单位为cm。(3)采用弹性层状体系理论计算各当量计算点的mises等效应力和等效应变，结果见表14。表14模量初始值时各当量计算点的mises等效应力和等效应变步骤4.将表14中的mises等效应力或等效应变代入表12中的路面材料模量应力(应变)依赖模型，计算各层材料新的模量值见表14。表15路面各层材料新的模量值步骤5.确定各结构层的最终模量ei。(1)比较和计算二者相对误差，见表16。由表16可知，二者相对误差较大，不满足收敛精度要求，故需要重复步骤3和步骤4进行迭代计算，直至满足收敛精度要求为止。表16和相对误差(2)迭代计算过程见表17。由表17可知，当迭代计算至第7次时，已满足收敛精度要求时，此时的各层模量值即为最终模量ei。表17迭代计算结果步骤6.采用步骤5中的各结构层最终模量ei，根据弹性层状体系理论计算沥青路面结构内部任意一点的应力、应变和位移，至此便得到基于路面材料模量应力(应变)依赖模型和当量方法的沥青路面结构分析结果。此处仅列出z轴上各点的最大主应力、最大主应变和竖向位移计算结果，见图6至图8。当前第1页12