ansys有限元程序在高路堤自身沉降变形分析中的应用.doc

PAGE PAGE 4 ANSYS有限元程序 在高路堤自身沉降变形分析中的应用 1 引言 随着西部地区高速公路建设的蓬勃发展，公路等级的不断提高，对公路建设的质量要求也越来越高。由于西部地区地形、地貌、地质条件的复杂，而公路又是线性结构物，路线长、占地广，高填深挖的路堤结构型式是不可避免的。在近几年来云、贵、川等西部地区的高等级公路建设中，大于20m填高的路段占有不少的比例，同时对建成通车的高填方路段的病害调查发现，多数路段在通车后不久均存在路面纵横向裂缝或不均匀沉降问题，因此对于在西部地区的高填方路基的修筑及沉降计算和工后沉降控制问题，已引起工程界的普遍重视。 路基沉降计算应包括地基沉降和路堤填料的自身沉降计算。在低填方路段，由于填方高度小，路堤填料在自重作用下的沉降相应较小，通常在沉降计算时不考虑填料的自身沉降问题，而只计算地基沉降。但对于高填方路基而言，不仅存在地基过量沉降问题，即由于填筑高度高而引起附加应力增大导致地基的沉降量过大，因此要求高填方地基应有具有足够的承载力；同时，对于填料在自重作用下自身压缩变形量也不容忽视。对于地基的沉降计算，目前已有大量的研究，对于高填路堤填料的自身沉降计算，公路、机场、水利部门的有关学者已进行了研究分析，主要采用的计算方法是对适用于地基沉降的理论分析方法进行修正，即改进的分层总和法、数值分析方法及根据现场实测数据进行拟合分析的方法。但总的来说，国内外对于高填路堤的沉降变形规律的研究还很少，还没有形成统一、完善的理论计算方法。 在应用数值分析法中的有限元法求解高路堤沉降变形时，需将研究区域离散成有限数据的区域单元，在每个单元内，假定位移场或应力场，应用变分原理建立代数方程组，然后结合初始条件和边界条件，以节点处的广义位移或广义应力作为未知量进行求解。有限元方法是一种近似的计算方法，应用有限元法计算路堤沉降时,需将地基和路堤作为一个整体来分析，对其进行网格划分，形成离散体结构，在荷载作用下算得任一时刻地基和路堤各点的位移和应力。其中基础底面的竖向位移就是所要求的沉降。有限元法可以将地基作为二维甚至三维问题来考虑，反映了侧向变形的影响。它可以考虑土体应力应变关系的非线性特性,采用非线性弹性的本构模型, 或者弹塑性本构 模型。目前使用最广泛的是邓肯—张双曲线模型。它不仅可以考虑应力历史对变形的影响（若应力低于或高于前期固结应力,则采用不同的弹性模量计算公式），还可以考虑土与结构的共同作用、复杂的边界条件、施工逐级加荷情况以及土层的各向异性等，同时对土体的固结计算,可以用Biot三维固结理论,避免了一维固结计算的许多弊端。从计算方法上来说，有限元法是一种较为完善的计算路堤沉降和变形的方法。它的缺点是计算工作量大，参数确定困难，要做三轴排水试验，目前主要用于重要工程、重点地段的计算【1】。因此，本文主要应用ANSYS大型有限元分析软件对高填方路堤的自身沉降变形规律进行分析探讨，寻找出高填方路堤沉降变形的基本趋势或规律，期望给修建高填方路基的工程技术人员提供一些有价值的参考意见。 2 ANSYS有限元分析软件概述 ANSYS是集结构、热、流体、电磁场、声场和耦合场分析于一体的大型通用有限元软件,广泛应 用于机械、土木、水利、机电、航天，冶金等领域。它具有强大的建模功能、非线性分析功能（可 进行几何非线性、材料非线性、及状态非线性分析）、智能网络划分功能及求解功能。ABSYS典型的分析过程由前处理（建立几何模型、定义单元类型、材料属性、对模型进行网络划分、施工荷载和约束）、求解计算（选择求解类型、进行求解选设定、求解）和后处理（从求解计算结果中读取 数据、对计算结果进行图形化显示、列表显示、进行各种后续分析）组成【2】。在用ANSYS进行路基沉降分析过程中，模型选取尤为重要，其中包括：有限单元的选取、有限元沉降计算方法和计算结果的收敛性判别。 3 高路堤的ANSYS有限元模型 在ANSYS有限元分析软件中，适用于土工材料如岩石、土壤等颗粒状材料分析计算的模型只有非金属材料模式中的Drucker-Prager屈服准则，它实际是对莫尔—库仑准则的近似。其流动准则可以使用相关流动准则，也可以使用不相关流动准则，其屈服面并不随材料的逐渐屈服而改变，因此没有强化准则。但是，它的屈服强度随侧限压力的增加而相应增加，其塑性行为被假定为理想弹塑性。该模型材料可以考虑由于屈服而引起的体积膨胀，但不考虑温度变化的影响。 3.1 Drucker-Prager屈服准则参数选取 Mohr-Coulomb准则是表征岩土剪切破坏的主要定律。事实上，如果将岩土体视为理想塑性材料，也可以把Mohr-Coulomb破坏准则视为屈服准则。Mohr-Coulomb准则以剪应力作为判断标准，但是这个剪应力并不是最大剪应力，而是相应于土