abaqus在岩土工程中的应用-边坡稳定分析

高等土力学边坡稳定分析专业：岩土工程姓名：xxx指导老师：xxx学号：xxx1 .前言边坡稳定分析是边坡设计的前提，它决定着边坡是否失稳以及边坡失稳时存在多大推力，以便为支护结构设计提供科学依据。然而这个问题至今仍未得到妥善解决，因为解决这一问题必须先要查清坡体的地质状况及其强度参数，同时又要有科学合理的分析方法1。对于均质土坡，传统方法主要有：极限平衡法，极限分析法，滑移线场法等，就目前工程应用而言，主要还是极限平衡法，但需要事先知道滑动面位置和形状。对于均质土坡，可以通过各种优化方法来搜索危险滑动面，但是对于岩质边坡，由于实际岩体中含有大量不同构造、产状和特性的不连续结构面，传统极限平衡方法尚不能搜索出危险滑动面以及相应的稳定安全系数。边坡稳定分析涉及复杂的地质地形边界条件、材料的应力-应变的非线性行为、初始地应力、水压力、地震荷载的耦合分析等等，多数情况下不能获得解析解。在计算机和计算方法不断发展的背景下，以有限元为代表的数值分析方法在20世纪70年代已逐步在岩土工程中推广应用，并发展成为一种强有力的计算分析工具。然而传统的数值分析方法一般只是得出边坡应力、位移、塑性区等，不能直接与边坡稳定建立定量关系。随着计算机技术的发展，有限元强度折减法近来在国内外受到关注212,对于均质土坡已经得到了较好的结论，但尚未在工程中实用，本文采用有限元强度折减法，对均质土坡进行了系统分析，证实了其实用于工程的可行性，得到了节理岩质边坡坡体的危险滑动面和相应的稳定安全系数。该方法可以对贯通和非贯通的节理岩质边坡进行稳定分析，同时可以考虑地下水、施工过程对边坡稳定性的影响，可以考虑各种支挡结构与岩土材料的共同作用，为边坡稳定分析开辟了新的途径。2 .有限元强度折减法原理c=c/,tan甲=tan邛/与这种方法早在70年代就提出来了。1975年Zienkiewize就利用有限元进行边坡稳定分析，但是由于受计算条件的限制，此法一直没有流行起来。近年来，这种方法随着计算机软件和硬件的发展又有了新的发展。3 .有限元强度折减法判断依据有限元强度折减法分析边坡稳定性的一个关键问题是如何根据有限元计算结果来判别边坡是否处于破坏状态。目前的失稳判据主要有两类:(1)在有限元计算过程中采用力和位移的不收敛作为边坡失稳的标志。(2)以广义塑性应变或者等效塑性应变从坡脚到坡顶贯通作为边坡破坏的标志。以上两种判据得到的安全系数相差不大4 .案例分析4.1 .工程概况3有一局H二13m坡脚为45度的均质边坡，土体容重20KN/m,粘聚力c=12.38Kpa,摩擦角20度。按平衡极限分析法，求土坡的稳定安全系数。4.2 模型的建立与求解4.2.1 建立部件在Part模块中建立，如下所示的部件下所示的4.2.2 设置材料及截面特性。4.2.3 装配部件4.2.4 定义分析步在step模块中建立load、reduce分析步4.2.5 定义载荷、边界条件。在Load模块中，定义载荷和边界条件，如下所示4.2.6 划分网格，在mesh模块中划分网格，如下所示图1网格图4.2.7 修改模型文件4.2.8 提交作业4.3 结果分析本算例在第二个分析步的0.23526时无法收敛，计算终止图2FV1随U1的变化关系强度折减到某一程度后，土体就已经失稳，既就是在t=0.23526时土坡失稳破坏。我们由图可以看出，以数值计算不收敛作为土坡稳定的评价标准对应的FV1=0.86,既就是安全系数为0.86。由于该图没有明显的该店所以不能以位移拐点作为评价标准。观察塑性区的发展，下图绘制出了时刻塑性区，我们可以看出在一开始土坡坡脚先出现屈服，然后向上延伸，直到二级坡的塑性区贯图3t=0.125的塑性区图4t=0.2328的塑性区从图可以发现滑动面的大致位置，与平衡极限法分析得到的一样，大致呈圆弧状，通过坡脚。用增量位移判断滑动面635555555a10:灯勺。灯灯1。5./-7.7有-1.-1,-1.-i三三4.4 理正岩土复核|1K8.LiLiLi10.LILILI12.80651.340022.0000.0002.0000.000033.0003.0004.24345,000043.LiLiLi0.LILILI3.0000.0000图6参数界面辘助功能-说明|倾斜土层土坡稳定计算2基本坡面土层水面加筋坡面线段数晋投感竖回投翱朋线城籍加角举曩官月长加)和沁隹加)(JW)不薮图7参数界面地功能说明倾斜土层土坡稳定计算21原点以上土国A1燃点即土国ft1重度粘结强度（h,）抗期附113.QW20,00012DODOJ市值12.360基本坛面土层水囱加防髓爵重度dtH/ma)W粘结强偿OtFfc）抗期雨粘果,婚L3.00D20,0001200D0C、中值12.3S0图8参数界面C-3.200,19.300)-195G3K=0840(-3.200/19.340)-19.603K=0840C-3.200,19.300)-19.64ZK=08如C-3.200,19一420)一19.662K=0.340C-3.200,19460)19.T21K=0.340(-3,200,19.500)-19.761K=0840图9运行计算界面与结果输出界面由图可以看出理正岩土计算的的安全系数为0.86,与有限元强度折减法求得的安全系数0.84基本接近，这就说明有限元计算结果的有效性。5 .结论有限元强度折减法不需要对滑动面的形状和位置进行假定，通过强度折减达到边坡的不稳定状态时，非线性有限元计算不收敛，此时折减系数就是安全系数，同时可得到边坡破坏时的滑动面。有限元强度折减法相对传统计算方法的优势；能够对具有复杂地貌，地质的边坡计算；考虑了土体非线性弹塑性本构关系，以及变形对应力的影响。