边坡稳定的剪切带计算朱以文.pdfVIP

第24卷第4期 　 计 算 力 学 学 报 　 Vol . 2 4 ,No . 4 2007 年 8 月 Chinese Journal of Computational Mechanics Augu st 2007 文章编号 (2007) 04044 106 边坡稳定的剪切带计算 1 1 ,2 1 1 ,3 朱以文 , 　徐 　晗 , 　蔡元奇 , 　朱方敏 ( 1. 武汉大学 土木建筑工程学院 ,武汉 430072 ;2 . 长江科学院 土工研究所 ,武汉 4300 10 ; 3 . 广东省水利电力勘测设计研究院 ,广州 5107 10) 摘 　要 :为了解决边坡稳定分析中剪切带有限元网格的依赖性问题 ,采用梯度塑性理论 ,从本构关系中引入特征 长度入手 ,建立计算模型 。提出了一种 8 节点缩减积分的梯度塑性单元 ,并采用梯度塑性理论推导了Drucker Prager 屈服准则的软化模型的有限元格式 ,在 ABA QU S 中进行了二次开发 ,嵌入了本文提出的 8 节点单元和本 构模型 ,并用 ABA QU S 软件进行了边坡剪切带的计算 。计算结果表明 ,本文提出的方法消除了经典有限元计算 的网格依赖性问题 ,可以得到与单元剖分无关的稳定的剪切带宽度 。本文所提出的方法可适用于其他场合的剪 切带计算 。 关键词 :梯度塑性 ;应变局部化 ;剪切带 ;边坡稳定 中图分类号 : TU43 1 　　　文献标识码 :A ( ) 1 　引　言 或软化 参数不在局部意义上确定 ,而是由积分点 及其邻域在非局部意义上确定 ,并在本构模型中引 应变局部化现象在岩土工程中大量存在 ,最为 入了材料内在特征长度参数 l , 即梯度塑性理论是 典型的就是土体失稳时形成的剪切带 。剪切带形 从本构关系入手来解决网格依赖性问题 ,它在应力 成的研究对于评价土工结构物的安全性和稳定性 应变关系中加入了塑性应变的梯度项 ,其屈服函数 等问题具有重要意义 。剪切带现象的本质是材料 的一般形式为 的不稳定性 ,材料不稳定研究中的一个十分重要的 2 2 ( σ ) ( ) f = f { } , k , l k = 0 1 问题是多尺度和标度律的问题[ 1 ] 。在材料发生不 式中 k 为硬化参数, 在局部点上定义, 可以用塑性 稳定性时 ,不同的物理过程 ,不同的微细观结构 , 2 功或者是等效塑性应变表示; k 为塑性应变的 在某个时间尺度内,将以其特有的动力学行为 ,在 梯度项, 梯度塑性模型中的梯度项补充了局部点邻 某个特征尺度上表现自