《路面设计原理》讲稿--沥青路面应力分析

第五章沥青路面应力分析一古典设计方法1麻省公式图51古典公式示意图年，美国麻省道路委员会第八次年会上发表了世界上第一个路面设计的公式。它假定汽车是一个集中荷载P，荷载以角通过碎石基层分布于边长为碎石层厚2倍的正方形面积的土基上，所以（51）QH212）（载荷中集度强载承基土中式PQ2公式年，对麻省公式进行修正，认为荷载在路面层内的传布与垂直方向成某一分布角的圆锥上，所以传到路面的顶面时，压力分布于一个圆形的面积上而不是正方形，但他仍假定汽车荷载为集中荷载。据此图52古典公式改进（52）PHTGQ20564载荷中集度强载承基土中式PQ3公式1934年、认为由于汽车荷载轮胎接触路面由一个面积，所以不应当假定汽车荷载为集中荷载，而应当假定汽车荷载为圆形均布荷载，并设轮载接地圆形面积的半径为A，即（53）PHTGAQ（）（）210564载荷中集度强载承基土中式PQ4评述古典理论公式是假定路面只要起分布荷载的作用，采用简单的分布角的概念，这个朴素思想的路面力学理论应予解决的问题；从各公式得知，路面厚度主要取决于土基承载力得大小，这就是土基强度得问题。但初期没有提出土基参数的测定问题；古典公式以轮载作为交通荷载，它不能反映交通量的因素，这在当时轻交通时代可能矛盾不突出，但随着交通得发展，不考虑交通量是无法使用的解决的办法就是在土基承载力取值上应根据交通量的大小采取不同的安全系数。二弹性半空间体解答过程布辛尼斯克得到完整的解答，方法是采用半逆解法。年势能法得到了解答。采用路面力学中的方法，同样可以得到解答。解轮隙弯沉的计算及应用采用以上公式2101212/1EPAWZARZ）（时，当）（时，当）（）（64202715RARRRFF时，当三多层体系解答过程1945年，得到理论解19451955研究层状体系的工程应用1955，希夫曼得到非轴对称的解计算方法采用查诺模图法采用程序计算法四计算程序沥青路面通常是多层体系。自从本世纪四十年代以来无论在理论分析，还是在数值计算方面，都取得很大进展，特别是计算机科学的发展及其在工程技术中广泛应用，使层状体系理论的研究的日趋完善，其中有波米斯特DMBURMISTER1945年及英因福克斯LFOX、阿堪姆WEACUM、苏联科岗KORAH及英国琼斯AJONES等所作的贡献。在荷载形式方面，包括轴对称均布荷载与非轴对称单向水平荷载，都可直接进行数值计算，在层次结构方面，由双层体系、三层体系发展到多层体系。在计算机程序方面，有壳牌公司编制的BISAR程序，雪弗隆公司编制的CHEVRON程序，美国地沥青学会所采用的DAMA程序。基本图式与基本假定多层体系在圆形均布垂直荷载作用下的计算图式如图1所示。层状体系基本假定1各层都是由均质、各向同性的弹性材料组成，这种材料的力学性能服从虎克定律；2假定土基在水平方向和向下的深度方向均为无限，其上的各层厚度均为有限，但水平方向仍为无限；3上层表示作用着轴对称圆形均布垂直荷载，同时在上层无限深度处及水平无限远处应力和应变都是零；4层间接触面假定完全连续。Z1,E1,1,H1Z2,E2,2,H2Z3,E3,3,H3ZN1,EN1,N1,HN1EN,N1RPR图1计算图式基本原理根据弹性理论，对于轴对称空间体，其几何方程为RZZRURWU；1其物理方程为RRZE1ZZRZRZRE212式中G为剪数模量，GE1为弹性体的泊桑比轴对称空间课题微分单元的平衡微分方程为RZRRZZ03从式1式3看出，三式中共有十个变量，并且已有十个方程式，结合边界条件即可解出未知量值。但这种解法相当困难，甚至不可能得到应力分量。因此一般采用应力函数求解。研究物体的变形一般是针对物体内部割出的一块微分单元体，显然各相邻单元体的变形应是谐调的。所以物体在变形前是一个连续体，在变形后也应是一个连续体。消去位移分量，可得变形连续方程为22104RRR第一应力不变量，Z变形连续方程又称相容条件，是由圣维南BDESAINTVENANT于1864年提出的。实际上式4应有四个相容条件，但确是等效的。采用应力函数法求解轴对称课题主要有LOVE函数法及SOUTHWELL函数法，这里介绍LOVE函数法。设应力函数R,Z并给定RZR222ZZ（）2ZR（）152将式5代入平衡微分方程式3和变形连续方程式4，除平衡微分方程中第一个恒等于零外，其余全部转化为重调和方程，即206这就是说，如果应力函数是重调和方程的解，则能满足平衡微分方程的变形连续方程。并可由式5求得应力分量，再由物理方程求得应变分量。位移分量可由下式求得。UERZWZ（）127222重调和方程的求解可采用分离变量法。对于多层体系中某一层J，可以给定应力函数为JJGJR08代入重调和方程可以得出DZZ229令，解以上方程式，可得应力函数为HZR；101203JJJJJJJJJDCEBEAJ式中参数；J第一类零阶贝塞尔函数；无量纲系数；JJZA,B,C,D为积分常数，可由每一层的边界条件和层间结合条件等确定。下标J从1到N，表示同该层次相应的计算参数。将应力函数式（10）代入洛夫应力函数与应力关系式。ZJJJJJJACEBDEJJ0122以上表达的各项分量并非因荷载所引起，而是由JR所引起。通过HANKEL变换等推演过程，可以求得（12）PJD0假设PJ0（13）D0假设式（11）的各项分量为Y，而由实际荷载P产生的各项应力、位移分量为Y，这两者有以下关系（14）YP0将式（11）代入15212101DEDBECAJJJJJJJJJJZ由Q这个实际荷载引起的各项分量。积分常数计算对N层体系具有4N个积分常数。由以上算式可以看出，多层体系的应力应变计算的关键是要确定对应于各个层次的积分常数，然后通过贝塞尔函数的无穷积分计算，便可完成全部计算分析工作。确定积分常数，可以根据相应的边界条件与层间结合条件来进行。在多层体系顶面J1,0具有以下边界条件）J0（6）在第J层与第J1层之间的结合面上J，若这两层是完全连续的，则具有以下连续条件（）（）U）UWW（7）此外，在地基的无限深处，应力与位移皆满足，,U,WR0（8）则要求应力函数0，因E，即A0（9）则对于N层体系，还有4N2个待定积分系数，而根据边界条件可以建立4N2个方程式，因此全部积分系数均可以求解。确定待定积分系数，用矩阵法非常简单，便于使用计算机分析计算。为此可将应力和位移中包含有A,B,C,D的系数写成矩阵形式ZJRJJJJJJJJUWMEHCD，20式中HM44的矩阵根据连续条件，可以写成EHABCDEHABCDJJJJJJJJJJ，11112由式21可以看出，第J层积分常数可由第J1层的积分常数求得。通过逐层计算，可以将第一层的积分常数与第N层的积分常数联系起来。并利用下式可得ABCDNBDCJNJNNNN11102由多层体系顶面的边界条件代入（11）得AEEHH11112B120则11111220111DCBAFEEHH003FCBDN因此，在计算积分常数时，可按以下步骤进行计算1形成矩阵2形成矩阵3计算B，D4由下而上逐层计算各层的积分常数。在积分常数确定之后，通过贝塞尔函数及无穷积分数值可计算应力分量及位移分量。数值积分在进行应力分量及位移分量计算时，可以归纳为以下的形式EFD0245贝塞尔函数的计算贝塞尔方程（25）XDYXNY220贝塞尔函数的解JKNNK1120N0JXXXK024211N1JXXK13221函数特性当时，（），其它各阶均为零是衰减函数，趣近于正弦函数贝塞尔函数的计算当小于时JXABXIII007214（）IAIBIIAIBI0124044425840177758313999999450070925300236617239999973266666616000767720002206931777756088888470000501400001290当大于时JXPTXTQXIII（）44025025COSSNTTTXII1105105334（）式中03989423001246690398942300374008000175310000456400029218000063900001734000008700002232000010650000048800000342000005810000039900000174000001420000020100000162000000370000003200000042000000376贝塞尔函数的无穷积分高斯积分根据高斯积分公式，在，上取个插值点，计算积分时有（）（）其中插值点为勒朗德多项式的个根，称为高斯型点，为高斯系数。111PFTDFTDAFTXTPAKKPKK如果积分区间不是，而是，则可进行变换，高斯系数也应作相应的变换，12XTBAKKK高斯系数高斯型点的选择，可以根据计算精度确定。T057735026921077459666920555555555560088888888889086113631160347854845103399810436065211451549090617984590236926885105384693101047862867050058888888893贝塞尔函数的积分为避免因贝塞尔函数的波动性而引起积分误差，在积分时采用贝塞尔函数零点积分的办法，即首先找到贝塞尔函数的零点，然后在该范围内进行高斯积分，则JXDFEKKNNK01控制精度对于两个贝塞尔函数的积分，方法同前，只是零点应包括两个贝塞尔函数的零点。计算程序五层状体系计算程序层状体系应力，应变及位移的计算时主要涉及的是贝塞尔函数的无穷积分，因此无穷积分的精度将直接影响计算精度。因为贝塞尔函数是波动衰减函数，如果采用第二章所述的高斯积分，必须合理选取高斯积分段。在一般的数值积分法中主要采用等值增长的办法，使计算相对误差达到规定的精度。但由于贝塞尔函数是波动衰减函数，如果高斯积分区段一端函数值为正，另一端函数值为负，高斯积分点或为正，或为负，那么计算结果误差则比较大，为了防止以上这种情况的发生，在程序中采用零点分段的办法。由于贝塞尔函数的零点为已知，那么零点与其它任何数相乘均为零，那么积分时选用的分段区间为相邻两零点之间，则积分结果精度较高。6求解N层体系积分常数AI，BI，CI，DI子程序SOL积分常数计算顺序由下而上进行，即由第N层的积分常数AN，BN，CN，DN计算AN1,BN1,CN1,DN1，然后逐层向上，直到希望计算的某一层。程序中，积分常数计算的主要任务是确定系数阵F及C。其执行程序为SOL。程序中符号说明TM相当于；HNH每层的结构厚度取总厚度的相对值；ZNH各层界面的竖向坐标取总厚度的相对值；NH结构层数不包括第N层土基，NHN1；PRN各层次泊桑比；EENH相当于REJJJ1VV4存放每个层次的积分常数AJ，BJ，CJ，DJ；LPTNH层次结构的顺序；A4，4系数矩阵J；C4，4系数矩阵C；F2，4系数矩阵F。7积分计算子程序同前8余项计算在某一区段内无穷积分积的执行程序可见多层体系计算程序。在程序中余项值的计算由控制变量INTT控制执行。当INTT值为零时，表明计算点在层状体系顶面Z0，则要求计算余顶值，当INTT为1时，表明计算点在层状体系顶面以下点Z0，则不要计算余项。余项计算值代表符VS11，RS11，TS11及W11分别表示Z，R，ZR及在有限积分段的余项值。9多层体系应力、应变及位移计算程序说明本程序适用于多层体系结构，对双圆或多圆均布荷载可采用应力迭加原理得到。执行程序如下程序中主要变量说明INTT计算点位置符INTT0，表示要求余项；IC积分次数；VS，RS，TS，分别存放每次积分后Z，,W的对应值DEL积分的相对精度WT，D分别为高斯积分宽度和高斯积分点。其它变量名同子程序SOL六、双圆或多圆荷载应力的计算现行的路面设计规范多用双圆均匀荷载，利用本程序进行应力计算时，必须对源程序进行修正。修正的基本方法是根据应力迭加原理。要求荷轮隙、轮印中心及轮印中心两侧的应力或位移，可利用计算两点位移或应力迭加的方法。七弹性多层体系应力、位移分析程序AP01程序功能弹性多层体系应力，位移分析程序适用于N层组成的多层结构体系，具有如下功能1适用于多层弹性体系，层数不限，在此最大值定为L6；2每个层次的弹性模量和泊桑比不受限制；3适用于计算各层体系任意一点的应力，位移计算，可同时算出多个点的应力及位移，计算点最大值为24点；4荷载为单园垂直均布荷载，作用于上层顶面；5对双圆荷载，则利用单圆荷载进行应力迭加。程序的输入输出说明输入变量NL层状体系的层数NL6；NS应力、应变及位移计算点数NS24；NIC积分最多次数NIC40；INTT计算点状态参数INTT0，表面点，INTT0，内部点；DEL近似积分的精度，常用00001；CR荷载圆的半径CM；CP荷载的单位接触压力KG/CM；RNS每个计算点离荷载中心的径向坐标值CM；ZNS每个计算点离表面的垂直坐标值CM；ENL每一层的弹性模量，KG/CM；PRNL每一层的泊桑比；HANH每一结构层的厚度NHNL1CM；输出变量STRESSZZ方向正应力KG/CM；STRESSRR方向正应力KG/CM；STRESSTT方向正应力KG/CM；STRESSZRZR方向剪应力KG/CM；STRAINZZ方向应变；STRAINRR方向应变；STRIINTT方向应变；DISPZZ方向位移，CM。3计算实例例1、某一弹性三层体系，有关参数如下3,10,1880,1,100001,707,10650,10,20,30,50,80,150,250,0,080,10,50025,025,03520000,10000,50010,404计算结果STRESSES,STRAINSVERTICALDISPLACEMENTSINANELASTICULTILAYERSYSTEMUNDERASURFACECIRCLEUNIFORMLOADNUMBEROFLAYERS3NUMBEROFPOINTS10CONTACTRADIUSOFLOAD707000CMCONTACTPRESSUREOFLOAD1065000KSCTOLERANCEFORINTEGRATION00010RADIUSCOORDINATESANDVERTICALCOORDINATESFOREACHCONPUTINGPOINTNORCMZCM1000000210000000320000000430000000550000000680000000715000000082500000009000100001000050000PARAMETERSFOREACHLAYERNOPREKSCHCM12502000010000225010000400003350500NUMBEROFITERATIONSIC5RETURN1STRESSE,STRAINVERTICALDISPLACEMENTFOREACHCOMPUTINGPOINTNOSTRESSZSTRESSRSTRESSTSTRESSZR1106492297697997697900000200000401682615010000030000001439108695000004000001100270087000005000001205843806000006000000076125757000007000000688908904000008000000559002340000009402456153290153290000001006498618276182700000NOSTRAINZSTRAINRSTRAINTDISPZ128822E0323325E0323325E0321370E01227667E0452772E0413577E0316424E01313407E0414306E0454527E0414085E01410136E0432601E0533668E0412942E01569829E0555306E0620396E0411447E01633147E0528392E0512783E0496647E02725193E0645575E0553133E0565413E02840623E0630877E0518690E0539669E02923955E0310779E0310779E0318071E011037411E0447994E0447994E0413806E01八计算结果分析（1）路基应力铺设路面结构层的主要目的是扩散车轮荷载，以减少传给路基的应力值，因为过大的应力值使路基出现剪切破坏或出现塑性变形，从而使路面结构破坏。图56是相对刚度不同的双层体系，沿荷载截面中轴上路基竖向应力系数随深度而变化的情况。图中可以明显Z地看出，在路面厚度不变的情况下，随路面材料刚度的增长E1/E0，路基的应力急剧减少，特别是路基顶面处的应力值降得更快。例如，在两层分界面处，按均质半无限体E1/E01计算所得的约为竖向应力的68，而设置模量增大9倍的面层后，约为竖向应力的ZZ30。利用三层体系的数值解，可以分析基层或面层的厚度和刚度对路基顶面竖向应力的影响。面层和路基的刚度不变时，竖向应力系数随基层刚度和Z厚度而变化的情况下（如图57A）所示。可以看出，随层Z刚度和厚度的增加而减少。面层刚度的影响如57B，路基应力随面层刚度的增加而减少。面层刚度很大时，基层厚度对路基应力的影响很小。由此可见，为把路基应力降到某一容许值，可以采用增加面层或基层厚度或刚度办法，其中增加刚度比增加厚度效果大。这个规律对于设计沥青路面的基层有重要的意义。采用粒料基层时，由于本身的模量值很低，只能通过增加厚度来减少路基应力；而采用刚度较大的稳定类基层，则可明显减少路基应力，并且在相同图56路基竖向应力系数随深度、E1/E0变化情况基层的厚度和刚度面层的厚度和刚度图57基层或面层的厚度和刚度对路基顶面竖向应力的影响的路基类型和容许应力（弯沉）条件下，其厚度可比粒料基层减少很多。（2）路面弯沉路面弯沉是路基和路面结构不同深度处竖向应变的总和。对于等级不太高的路面来说，其中7095由路基提供，各点的应变是三向应力状态的函数，因此，影响路基应力的诸因数也会影响路面弯沉。图58给出了三层体系荷载面中轴处的表面弯沉系数随层厚和模量而变化的情况。0W增加面层或基层的厚度都可减少路面弯沉；但在面层或基层厚度较薄时，增加厚度对降低弯沉量的影响比层后大时显著得多。也可通过增加路基、基层或面层的刚度使路面弯沉量降低。对比图中曲线变化可以看出，在路基刚度低时，路基刚度对弯沉量的影响要比基层和面层的影响明显得多。（3）基层底面的拉应力采用刚度较大基层将提高荷载扩散能力，使路基的应力和弯沉量减少。但是随着基层相对刚度的增大，基层底面的拉应力增大。此拉应力如果超过材料的抗拉强度，基层将会断裂，并导致路面破坏。图59给出了在面层相对刚度和厚度不变时，基层底面拉应力系数随基层相对刚度和厚度而变化的情况。可以2R看出，增加基层的相对刚度，将导致增大，而在2R基层较薄时，刚度对的影响要比厚基层严重得多。2R因此，为降低路基的应力或路面弯沉值而选用相对刚度较大的基层时，应验算基层底面的拉应力，使材料的抗拉强度与之相适应。基层底面最大拉应力位置一般在荷载作用面中轴处；在双圆荷载作用下，则出现在其中一个荷载作用面的中轴处。（4）面层的底面的拉应力垂直荷载作用下，面层底面的径向应力并非都是拉应力（如图510）。面层较薄而相对刚度较小时（E2/E1035，H1/01，H1/025），可能出现拉应力。面层较厚和刚度较大时，面层底面便出现拉应力。它随面层相对刚度的增大而增大，特别是面层的相对刚度很大时（E2/E103），拉应力随刚度的增大而急剧增大。底面最大拉应力的位置，一般在荷载面的中轴处；双圆荷载时，最大拉应力一般出现在某一荷载面的中轴处，但在面层很厚时，随层厚增增大而移向双圆图58A）厚度影响（E2/E1E0/E202）B路基和面层刚度的影响（H1/05H2/10）图59基层底面拉应力系数（E2/E102，H1/05）图510面层底面拉应力系数（E2/E010，H2/2）荷载面的对称轴处。在圆形均布的单向水平荷载作用下，面层内会出现较大的径向拉应力。特别在路面荷载作用面的边缘处，其数值很大（图511）。面层较薄时，其底面也会出现较大的径向拉应力（512）。（5）剪应力增加上层的刚度，还将导致层内剪应力的增加。垂直荷载作用下，面层内任一水平面上的最大剪应力一般出现在通过荷载作用面边缘的垂ZR直线上（图513）。在面层相对刚度增大时，最大剪应力随深度而变化的情况最大应力出现在ZR面层中部，并随面层刚度的增大而增大。但在面基的分界面上，最大剪应力最大剪应力随面层刚度的增加而减少。面层的厚度对最大剪ZR应力也由很大的影响。在面层和基层的相对ZR刚度不变的情况下，随面层厚度的减少，最大剪应力增加。其最大值出现的位置逐渐上移。ZR由二分点上升到三分点附近。因此，面层相对刚度很大而厚度较薄时，垂直荷载将产生较大的剪应力。应采取措施以防止面层出现出现较大的应力。面层受到圆形均布的单向水平荷载作用时，面层内各水平面上受到的最大剪应力随深度的增ZR加而衰减的很快（514）。面基曾分界处，最大值已下降大水平力的不到10。而在基层底面，最大剪应力已经可以忽略不计。ZR垂直荷载和水平荷载共同作用下，面层内最大剪应力也随深度的增加而减少，并随面ZR层相对刚度和厚度增加而增加。路面的最大剪应力出现在荷载面边缘处。其值主要受水平大小的影响，同时也受面层厚度和刚度的影