ANSYS在溷凝土基层沥青路面结构计算中的应用

会计学1ANSYS在溷凝土基层沥青路面结构计算中的应用1有限元法和ANSYS简介

ANSYS软件是融结构、流体、电场、磁场、声场和耦合场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS公司开发，是美国机械工程师协会、美国核安全局及近20种专业技术协会认证的标准分析软件。ANSYS软件主要包括三个部分:前处理模块,分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具可以方便地构造有限元模型。在加载求解过程中，包括位移、力、温度在内的任何荷载均可以直接施加在任意几何实体或者有限元实体上，载荷可以是具体数值，也可以是与时间或者坐标有关的任意函数。后处理模块可将计算结果以彩色等值线显示、梯度显示、立体切片显示等图形方式显示出来,也可将计算结果以图表、曲线形式显示或输出。软件提供了100种以上的单元类型,用来模拟工程中的各种结构和材料。ANSYS不仅具有结构静力分析、结构动力学分析、结构非线性分析、动力学分析、热分析等，而且还具有优化设计、建立子结构模型等高级应用功能。第1页/共35页2混凝土基层沥青路面的应力分析混凝土基层沥青路面以普通水泥混凝土基层和贫混凝土基层作为路面基层，上面加铺沥青混凝土面层，这种路面结构可弥补柔性路面和水泥混凝土路面的不足，大大改善路面的使用性能，延长路面的使用寿命。沥青面层为表面功能层，主要是提高路面表面的平整度、耐磨和抗滑性能，并扩散车轮荷载，减少冲击振动，防止反射裂缝的产生和发展。基层是车辆荷载和温度应力的主要承载层。在行车荷载和温度变化的反复作用下，混凝土基层内会产生荷载疲劳应力和温度疲劳应力。当行车荷载和温度变化产生的荷载疲劳应力和温度疲劳应力之和超过混凝土基层的弯拉强度时，基层板将发生断裂，继而整个路面结构破坏。因此，应计算混凝土基层的荷载疲劳应力和温度疲劳应力，且两者之和不应超过混凝土的设计弯拉强度，以控制因基层板断裂而造成的结构破坏。路面结构不但受到行车荷载的反复作用，而且在自然环境下受到太阳辐射、气温等因素不断变化的影响。因此路面结构内的温度是随着气温日变化、年变化而变化。当温度变化时，混凝土基层在温度作用下产生的应力可分为：由于基层顶与基层底温度不一致而产生的翘曲应力；由于温度均匀上升或下降引起基层的热胀冷缩而在基层内产生的胀缩应力。胀缩应力对温度应力的影响不大，且可在接缝设计时考虑，本文主要考虑温度翘曲应力。第2页/共35页3路面结构的建模和计算3.1几何模型和计算参数混凝土基层沥青路面结构可视为由沥青面层、混凝土基层及地基组成的弹性层状体系。沥青加铺层、混凝土基层和地基础结构层采用8节点实体单元。上下层间的接触状况为，混凝土基层与地基之间假设为绝对光滑，与沥青加铺层之间假设为完全连续。模型的边界条件为：（1）地基为弹性半空间体，根据收敛性分析，取为有限尺寸。基层为有限尺寸；（2）地基的底面为自由度完全约束的固定面，基层和地基的其它面均为自由悬空面。收敛性分析时采用标准轴载BZZ-100，轮胎内压0.7MPa，单个轮压作用范围18.9cm×18.9cm，接触面积357.21cm2，双轮间距32cm，两侧轮隙间距182cm（图1）。荷载作用位置为临界荷位，即纵向边缘中部，轮载作用如图2所示。第3页/共35页

图1车轮布载平面图（单位：cm）图2临界荷位处轮载作用示意图第4页/共35页

计算参数范围的确定：对混凝土基层沥青路面进行有限元分析，先要确定计算参数的大小和范围。主要计算参数如表1所示。表1主要计算参数结构层厚度（cm）弹性模量（MPa）泊松比导热系数（w/m℃）线膨胀系数（1/℃）沥青混凝土面层5～23600~22000.251.22.1×10-5混凝土基层12～308500~355000.151.51.0×10-5地基－50~9500.301.00.5×10-5第5页/共35页通过收敛性分析，结合实际应用情况，确定混凝土板块尺寸为10.0m×4.5m。此外，已有研究表明，当地基尺寸比混凝土板长宽方向均扩大3m，深度取为9m时，混凝土板底的荷载应力早已收敛，因此地基的尺寸取为16×10.5×9m。通过上述分析，最终确定混凝土基层沥青路面进行有限元分析计算的模型如图3所示。与之对应的路面整体有限元网格划分如图4所示。第6页/共35页图3结构计算模型（单位：cm）第7页/共35页图4路面结构有限元网格划分图第8页/共35页表2列举了沥青层厚度ha为0～20cm,基层厚度hc为22cm时，各自然区划内最大温度梯度变化范围。表2不同自然区划最大温度梯度推荐值(℃/cm)ha(cm)自然区划ⅡⅢⅣⅤⅥⅦ00.88～0.940.95～1.000.93～0.980.90～0.960.92～0.970.98～1.0440.62～0.660.65～0.710.64～0.700.63～0.680.64～0.690.69～0.7480.43～0.460.48～0.510.46～0.500.44～0.470.45～0.490.50～0.53120.30～0.320.32～0.340.31～0.330.30～0.320.30～0.330.33～0.35160.30～0.310.32～0.330.31～0.320.30～0.310.30～0.310.33～0.34200.30.320.310.30.30.33第9页/共35页当基层厚度变化时，可按表3进行修正。表3不同基层厚度最大温度梯度修正系数在进行热与结构的耦合分析中,主要考虑温度相对稳定时的翘曲应力,将各层接触面温度作为约束条件,不考虑空气的对流和路面的辐射,通过计算可求得各节点的温度分布,另外也可以输入辐射量、对流系数等参数直接求出温度分布,但需考虑一天中温度近似于正弦分布,较复杂,可操作性不强。ANSYS软件中对于热应力分析问题推荐采用顺序耦合分析,先进行热分析计算,将热分析所得的节点温度施加到结构单元的节点上,再求解温度应力,施加汽车荷载后可单独求汽车荷载应力或汽车和温度共同作用下产生的应力。为便于与规范计算值对照,先用ANSYS分别求出温度应力、荷载应力,再探讨温度场和行车荷载对路面所形成的复合应力。基层厚度（cm）1618202224262830修正系数1.241.151.0810.930.870.810.76第10页/共35页材料的几种强度理论（1）最大拉应力理论这一理论认为最大拉应力是引起材料脆性断裂破坏的主要因素。即不论材料处于什么应力状态，只要最大拉应力达到了材料在简单拉伸破坏时的极限应力，就会发生脆性断裂破坏。（2）最大剪应力理论这一理论认为最大剪应力是引起材料塑性屈服破坏的主要因素。即不论材料处于什么应力状态，只要危险点最大剪应力达到了材料在简单拉伸屈服时的极限剪应力，就会发生屈服破坏。（3）形状改变比能理论这一理论认为形状改变比能是引起材料塑性屈服破坏的主要因素。即不论材料处于什么应力状态，只要危险点处形状改变比能达到简单拉伸屈服时的形状改变比能，材料就会发生屈服破坏。第11页/共35页3.2具体计算步骤3.2.1热分析1)首先选定单元类型,本分析建立的是三维模型,热单元选用Thermalsolid70。2)在材料参数设置时,将材料的回弹模量、泊松比及导热系数输入。3)建立三维模型。对于层间光滑接触状态，引入正交各向异性单元模拟。选取正交各向异性接触单元的参数如表4所示。对于层间完全连续,用体与体的粘结(glue)。表4正交各向异性单元的参数取值说明：正交各向异性单元的厚度取为0.1cm参数取值弹性模量（MPa）泊松比剪切模量（MPa）ExEyEzμxyμyzμxzGxyGyzGxz0.00010.00011000000000.00010.00010.0001第12页/共35页对于层间光滑的情况，主要应考虑E33取值尽可能大，G12尽可能小。但为保证接触单元的本构关系矩阵为正定矩阵，上述参数的取值需满足根据余应变能理论推导出的模型约束公式：上述条件是理论分析后得出的应该满足的x方向各参数间的相互关系，y方向上参数间的关系可类似地导出。第13页/共35页对接触的更好的模拟方法：采用面面接触分析，用目标单元Targe170来模拟混凝土基层的上表面，用Conta173来模拟沥青面层下表面。在摩擦类型中选用基本的库仑摩擦模型，在这个模型中，两个接触面在开始相互滑动之前，在它们的界面上会有达到某一大小的剪应力产生，这种状态称为粘合状态。库仑摩擦模型定义了一个等效剪应力τ，在某一法向压应力p作用下剪应力达到此值时表面开始滑动（τ=μp+COHE，其中μ是摩擦系数，而COHE是粘聚力）。根据已有的文献摩擦系数μ采用1.8，粘聚力COHE为0.07Mpa。一旦剪应力超过此值后，两个表面之间将开始相互滑动。这种状态，叫做滑动状态。粘合、滑动计算决定什么时候一个点从粘合状态到滑动状态，或从滑动状态变到粘合状态。第14页/共35页4)把材料属性施加到体上,然后划分网格。5)按温度梯度施加温度荷载。按正温度梯度设置,即板顶温度比板底温度高,板边向下翘曲,板底面受拉应力。6)进行稳态热分析求解，观察结果。3.2.2

温度应力和荷载应力耦合分析1)将热单元转换为相应的结构单元。2)在材料特性中将热膨胀系数加入。3)读入热分析结果。5)施加约束,定义边界条件，模型中设置地基底面自由度完全约束。6)求解。如果单独求汽车荷载产生的应力,可直接建立结构体单元施荷加约束求得。如果要求温度和汽车的复合应力,可在3.2.2计算中加入施加汽车荷载的步骤。本模型中都采用线弹性材料,对于非线性材料,ANSYS也有强大的求解功能。第15页/共35页3.2.3后处理ANSYS有两个后处理器,通用后处理器(POST1)和时间历程后处理器(POST26)。对于静力学分析,可在POST1中图象显示结果或列表显示结果,对自己感兴趣的位置进行应力、位移分析。也可以指定特定的路径,对路径上各点的应力、位移等指标的变换规律进行研究。时间历程后处理器可用于研究指定点的分析结果与时间的函数关系。它有强大的分析能力,可进行图形显示、列表显示、微分、响应谱生成等操作。在动力学分析中,可通过定义变量,提取各结点、单元的位移、应力随时间的变化关系。为了加快后处理的速度,可以用ANSYS的命令编制命令流文件,通过选择、组合命令提取特定位置的结点或单元,将其随时间变化的数值映射到数组中,另存为自己定制的文件中。第16页/共35页选取以下的参数对混凝土基层进行有限元计算：面层厚度为14cm，模量为1800MPa，基层厚度为21cm，模量为25000MPa，地基模量为600MPa。温度梯度为0.3℃/cm。3.3荷载应力分析第17页/共35页基层横向受力图第18页/共35页第19页/共35页基层底部沿纵向边缘受力图第20页/共35页第21页/共35页3.4温度应力分析第22页/共35页基层横向受力图第23页/共35页第24页/共35页基层底部沿纵向边缘受力图第25页/共35页第26页/共35页3.5耦合应力分析第27页/共35页基层横向受力图第28页/共35页第29页/共35页基层底部沿纵向边缘受力图第30页/共35页第31页/共35页4结论1)利用大型有限元软件ANSYS,通过层面温度约束可方便快捷地求出结构层各节点的温度分布,通过其耦合功能,可求出最大温度应力及其位置。在施加汽车荷载作用下,可单独求出最大汽车荷载应力或汽车和温度场的复合应力。利用ANSYS的强大后处理功能,可观察一定横截面上的应力分布图。进一步可分析面层厚度、基层厚度、基层模量、地基模量和温度梯度等参数对应力的影响，并回归出荷载和温度应力的实用计算公式，用于路面结构设计。2)利用本文所建模型求出的温度应力和荷载应力与规范计算值有一定的出入,这主要是由网格划分、约束条件、参数选取、模型尺寸等因素造成的,可以通过施加约束方程、调整边界条件、建立接触单元等办法使解更接近真实值