ADINA、ANSYS、FLAC的强度折减法.docxVIP

PAGE 1 PAGE 1 ADINA、ANSYS、FLAC的强度折减法  与传统极限平衡法相比，强度折减法有众多优势。计算程序和屈服准则是影响强度折减求安全系数的重要因素。以某复杂边坡为例，通过ANSYS—ADINA和ANSYS—FLAC接口程序，建立ADINA、ANSYS、FLAC同套网格模型，克服网格对安全系数的影响。计算中分别选用M-C和D-P准则，分析得出：相对于D-P准则，M-C准则求得的安全系数偏于保守；对于一系列的D-P准则，外角点外接圆D-P和内切圆D-P求得的安全系数分别对应上限和下限；对于同种屈服准则，有限差分法比有限元法求得的安全系数略大。  在边坡稳定性分析中常用的方法是传统的刚体极限平衡法，但极限平衡法具有忽视边坡岩土体本身是变形体、边坡的破坏过程是其内部应力分布和变形不断调整的过程这一缺点。此外，滑动面必需事先假定的这个前提也严重限制了上述方法的应用范围。目前强度折减法得到普遍重视，但绝大部分研究集中于较为简洁的边坡模型，而且选用的计算程序繁多，计算结果的牢靠性尚未得到证明和比较。另外，屈服准则是影响强度折减求安全系数的一个非常重要的因素，对于岩土工程，Mohr-Coulumb(M-C)准则和Drucker-Prager(D-P)准则最为常见，D-P准则又有多种拟合方式，是一系列D-P准则的合称。对于屈服准则对强度折减求安全系数的影响，较为全面的分析和研究尚不多见。  鉴于上述几点，本文基于ADINA、ANSYS、FLAC，选取某工程边坡，通过子程序接口，建立ADINA、ANSYS、FLAC网格模型，在同等条件下利用强度折减法分别计算安全系数。选用ADINA、ANSYS、FLAC，基于以下考虑：  (1)ADINA、ANSYS为有限元软件，FLAC为有限差分软件，可以进行有限元强度折减和有限差分强度折减的比较；  (2)ADINA、FLAC内置M-C模型，ANSYS内置D-P模型，可以进行M-C准则强度折减和D-P准则强度折减的比较。  1强度折减法基本原理  对于边坡潜在滑动面为曲面时，边坡的稳定安全系数满意下列关系式：  强度折减法分析的实质是不断地按下式调整岩土材料的c、ψ值(为保证折减效率和安全系数的精度，折减间隔为o．02)，直至边坡失稳：  式中：c、ψ——粘聚力(MPa)和内摩擦角(°)；  c＇、ψ＇——经过折减后的粘聚力(MPa)和内摩擦角(°)；  F——折减系数。  强度折减法的失稳判据主要有3种：①以收敛性判定失稳状态；②通过计算域内塑性区是否贯穿判定失稳状态；③依据计算域内关键点位移与折减系数之间关系曲线变化特征判定失稳状态。说明如下：对于判据①，收敛性与收敛精度亲密相关，文中的收敛精度严格遵循各计算软件的默认收敛精度；对于判据②，FLAC采用塑性剪应变增量的方式；对于判据③，关键点通常选在简单发生塑性流淌的区域，对于不同计算软件，塑性区推断方式不同，关键点难于统。  2屈服准则  2．1M．C准则  M-C准则可表示为：  式中δ1——岩体破坏时的最大主应力(MPa)；  δ3——岩体破坏时的最小主应力(MPa)。  2．2D-P准则  D-P屈服准则是假定材料的破坏取决于第一应力不变量和第二应力不变量，其破坏准则可表示为：  式中：a、k——与岩体抗剪强度参数c、ψ有关的系数。  D-P准则与M-C准则有多种不同的拟合方法。表1给出了部分常见的D-P准则的系数，图1为各屈服准则在π平面上的曲线。 表1D-P准则的系数图1各屈服准则在π平面上的曲线 为了说明各种D-P准则对安全系数的影响，选取具有代表性的外角点外接圆D-P和内切圆D-P，分别代表半径最大的圆和半径最小的圆。由于ANSYS未内置内切圆D-P，可以通过以下参数转换方法间接地运用内切圆D-P准则。  设外角点外接圆D-P准则中粘聚力、内摩擦角分别为c、ψ，内切圆D—P准则中的转化粘聚力、内摩擦角(ANSYS输入参数)为c1、ψ1，二者的转换关系如下:  从而可求得内切圆D—P准则的ANSYS输入参数c1、ψ1的值。  3有限元模型及安全系数  计算剖面取自向家坝水电站某边坡工程，首先利用ANSYS建立有限元模型(图2)；将ANSYS节点和单元文件通过ANSY-ADINA接口转换，建立相同网格的ADINA有限元模型(图3)；类似的，通过ANSYS-FLAC接口建立FLAC有限差分模型(图4)。ADINA、ANSYS、FLAC底部均采用竖向约束，两侧采用水平向约束。 图2ANSYS网格