基于ABAQUS对DruckerPrager系列准则的讨论.pdfVIP

第28卷第11期 水电能源科学 V01．28No．11 Water 201 RCSOtIrCCSandPower Nov．20l0 0年11月 文章编号：1000—7709(2010)I卜0070—03 张艳山 潘玉珍 (长江三峡勘测研究院有限公司(武汉)，湖北武汉430074) 摘要：分析了Drucker—Prager(D—P)系列准则的应用条件及相互问关系，导出了各准则间的转换公式，结合算 例基于强度折减法，采用位移突变判据和塑性区贯通判据，分别计算分析了M—C和D—P准则下的安全系数， 并探讨了各准则的选取对计算结果的影响。 关键词：Drucker-Prager系列准则；Mohr-Coulomb；强度折减法；摩擦角I黏聚力 中图分类号：TV642．4；TU457 文献标志码：A 为更好地描述岩土材料的强度特性，在Mohr-忌及各准则的物理意义，见表1。 1．2 Coulomb(M—C)准则的基础上，出现了Drucker—D—P系列准则间的等价变换 Prager(D—P)系列屈服准则，且在岩土工程领域有限元软件ABAQus仅提供了M—C准则 有限元软件(ABAQUS、ANSYS、FLAC等)中得外角点外接圆的D—P。屈服准则。由于D—P系歹lj 到广泛应用。但D—P准则的参数选取不当，有可 准则在7(平面上为一系列的同心圆，其半径为： 能导致计算结果误差较大，甚至错误。文献[1]研 r=瓜一拉(五一aI，) (1) 究了边坡稳定有限元分析的影响，表明D—P的选 式中，J。为应力偏张量第二不变量；It为应力张 取直接影响边坡稳定有限元分析结果，不同D-P 量第一不变量；a、忌为与岩土材料内摩擦角钆粘 计算的边坡位移场和应力场差异显著。鉴此，本 聚力C有关的参数。 文分析了D-P系列准则的基本特性，导出了各准 r与足、a、I。有关，对特定的应力状态下，f， 则间的转换公式，基于强度折减法，采用位移突变 不变，因此可通过调整愚、口实现D-P准则间的等 判据n]和塑性区贯通判据[3]，结合典型的边坡算价变换[4|。以D—P。、D—P。准则为例：①对一定的 例分别计算了M—C和D—P系列准则下的安全系岩土材料参数C、9，若利用D—Pz准则计算，由表 数，揭示了应用各准则计算误差原因，为数值计算 1可知a。、忌。表达式；②若利用D—P．准则计算， 参数正确选取提供理论依据。 由表1可知口。、矗，表达式。③令口。=cr2、忌t=如， 即可求出fz、伫，并代入D—P，准则计算，与用f、妒 1 D—P系列准则 采用D—P。准则计算等价。因此，在仅提供D—P。 屈服准则的软件ABAQUS中用Cz、仇便实现了 1．1 M—C与D-P准则间参数换算D—P2准则的有限元计算。④同理，由表1计算可 求得弘、ct(i。1，2，…，5)；⑤将仍、ci代入表1对 文献[4]对比分析了M—C和D—P准则，得出 应力罗德角为两准则相匹配的对应参数，由此可 应的a．k，由此可实现D—P系列准则的转换。D—P 导出与M—C匹配的D—P系列准则的对应参数a、系列准则间的转化详见文献[4]。 表1 各准则参数换算及使用条件 anduseconditions Tab．1 ofdifferentcriterions