粉煤灰路堤稳定性数值仿真分析.docVIP

PAGE PAGE 1 粉煤灰路堤稳定性数值仿真分析 　　摘要：粉煤灰填筑路堤是解决粉煤灰再利用的最有效途径,同时因为粉煤灰是一种轻质填料,对减轻土体结构自重,减少软土地基路堤沉降,提高土体抗剪强度,均有很好的效果。本文主要对粉煤灰路堤稳定性作数值仿真分析，希望对于粉煤灰路堤的设计以及施工提供一定的指导作用。 　　关键词：粉煤灰路堤；稳定性；数值分析 　　中图分类号：P618.11文献标识码：A文章编号： 　　0前言 　　边坡的稳定分析是岩土工程中的重要研究课题之一。边坡稳定分析的计算方法有很多，如条分法、数值分析方法、塑性极限方法和模糊数学方法等。目前，对于土质边坡来说，为人们所熟知和广泛讨论的方法有条分法和有限元方法。条分法因为其力学模型简单适用，可对边坡进行定量的评价，经过长期的工程实践和不断的完善和补充，己经成为边坡稳定分析的成熟方法。近期在将有限元计算成果与传统的安全系数结合起来，直接用于边坡设计的判断依据的研究获得了一定的成果，如建立在强度缩小有限元分析基础上的边坡稳定分析原理等，本文将利用这些理论成果，结合有限元分析软件ANSYS，进行高填方粉煤灰路堤的边坡稳定性分析。 　　1计算与分析 　　1.1计算模型的建立 　　土的力学性质非常复杂，影响其应力和变形的因素很多。例如粉煤灰的结构、孔隙、密度、应力历史、荷载特征、孔隙水及时间效应等，在如此之多的因素作用下要获得理论解几乎是不可能的。而用ANSYS可以很好地模拟粉煤灰的力学性能，本文计算采用的是Drucker-Prager弹性—理想塑性模型，在国外被称为德鲁克—普拉格准则（D-P准则）。 　　式（1） 　　式中，，分别为应力张量的第一不变量和应力偏张量的第二不变量。 　　有限元强度系数折减法的基本原理是将坡体强度参数：粘聚力C和内摩擦角Φ值同时除以一个折减系数F，得到一组新的C′、Φ′值，其中： 　　C′=C/F式（2） 　　Φ′=arctan(tanΦ/F)式（3） 　　然后作为新的资料参数输入，再进行试算，利用相应的稳定判断准则，程序可以自动根据弹塑性计算结果得到破坏滑动面，确定相应的F值为坡体的最小稳定安全系数，此时坡体达到极限状态，发生剪切破坏，同时又可得到坡体的破坏滑动面。上述屈服准则表示为 　　式（4） 　　式中为达到极限状态时的安全系数。 　　计算时，首先选取初始折减系数，折减土体强度参数，将折减后的参数作为输入参数，进行有限元计算，若程序收敛，则土体仍处于稳定状态，然后再增加折减系数，直到不收敛为止，此时的折减系数即位边坡的稳定安全系数，此时的滑移面即为实际滑移面，这种方法称为土体强度折减系数法。 　　1.2材料参数的选取 　　选取粉煤灰的材料参数为E=48.2MPa，泊松比u=0.3,容重。选取地基土材料参数c=18KPa,=28.85°，E=80MPa,泊松比u=0.25，容重。根据强度折减法的计算公式，列出下列计算参数、值，如表1所示。 　　表1强度折减法参数计算表 　　1.3ANSYS计算概述 　　①问题的描述 　　分析均值粉煤灰纯灰路堤边坡的稳定性，所取的材料参数以及模型的坡比都趋向于比较不利情况，能够模拟实际情况中粉煤灰边坡的稳定性。计算边坡结构层的最大水平应力、水平剪应力以及整个路堤的水平位移和竖向位移，得出等值线图，对比分析并总结出高填方粉煤灰路堤边坡的稳定性数值研究，并与离心试验和模型试验所得的数据互相印证，为施工规范的编制提供依据。 　　②建模及网格的划分、加载 　　建模采用从下向上的建模方式，即用低级别的图元如点、线、面来创建高级别的图元如线、面、体，这样建模可以创建出较为复杂的几何体，而且还可以较为方便的读取关键节点位置的应力、应变值。ANSYS中化分完后的横截面如图1所示。划分完网格以后，对模型施加约束和进行加载，在施加约束方面路堤两个侧面限制了x方向的位移，路堤底面限制x,y方向的位移，以模拟路堤实际的受力情况。本计算只考虑路堤自身的重力作用，在ANSYS中只需将整个模型置于重力场中即可。单元数总共600个。对模型施加约束和加载后的横截面如图2所示。 　　图1网格划分后的模型横截面图图2施加约束和荷载后的模型横截面 　　1.4计算结果分析 　　①路堤模型整个断面竖向位移和水平位移的计算与分析 　　针对不同的折减系数情况下，计算出各个对应抗滑安全系数下整个断面的竖向位移和水平位移，计算结果如表2所示。图3和图4给出了整个路堤水平位移的最小值和最大值以及整个路堤竖向位移的最大值与抗滑安全系数的关系。系数为1.4时，模型路堤的水平位移和竖向位移图。由此我们可以得到以下结论： 　　1）从路堤水平位移效果图可以看出，随着抗滑安全系数的增大，水平方向位移呈增大趋势，但总体来看，水平位移均不大，最大位移均出现在路堤坡角处，然