土体三维弹塑性分析中的体积闭锁分析及单元选择-中国科技论文在线.PDFVIP

中国科技论文在线 第17卷第 6期 工 程 力 学 Vol.17 No.6 2000年 12 月 ENGINEERING MECHANICS Dec. 2000 文章编号 1000-4750(2000)06-006-08 土体三维弹塑性分析中 的体积闭锁分析及单元选择 杨 军 宋二祥 陈肇元 清华大学土木工程系 北京 100084 摘 要 在材料不可压缩或胀/缩塑性流动情况下 传统低阶单元有体积闭锁问题 以摩尔 库仑模型为例 推导了塑性剪切应变和塑性体积应变的关系 揭示闭锁产生的原因 分析8节 点等参元 Wilson非协调元 EAS单元和14节点单元的闭锁性态 表明8节点单元有严重闭锁 性 Wilson非协调元也有闭锁性 EAS单元和采用降阶积分的14节点单元能克服闭锁 单元测 试和方形基础的承载力计算两个数值算例证实了分析的结果 为土体三维分析中选择有效可 靠的单元提供依据 关键词 体积闭锁 有限元 三维分析 弹塑性分析 中图分类号 TU435 文献标识码 A 1 引言 随着工程的复杂化和计算手段的进步 人们已经不满足于对土体的二维分析 而在土 体三维有限元分析中 单元形式选择的重要性和复杂性较之二维分析更为突出 低阶单元 (如标准8节点等参元)便于网格生成 易于与其它单元类型(如板单元 壳单元 梁单元 界面单元等)联合使用 且生成的总体刚度矩阵带宽小 高阶单元(如20节点等参元) 由于自 由度太多 且使总体刚度矩阵的带宽过大 在实际问题中难以应用 但是 低阶单元往往有体积闭锁的缺陷 这可以被看做是离散体系不能正确反映由荷 载 边界条件和本构关系引起的连续性变形 研究表明 对于不可压缩 塑性体积应变为 零或胀/缩塑性流动的材料都有体积闭锁问题 为了解决上述低阶单元的闭锁问题 已提出不少方法 如用降阶积分和采用较高阶的 单元 但降阶积分可能会使单元刚度矩阵的秩降低 出现零能变形模式 用高阶单元将增 加计算量 有学者提出假定应变单元 其单元应变插值与位移是独立的 如B-bar方法[1] 最近 Simo提出一种叫做改进假定应变(EAS)的单元 单元应变假定为满足协调条件的应 变和附加应变两部分之和[2] [3] 此外 Wilson为提高单元精度 提出一种非协调元 Smith提出一类14节点单元并用 [4] 于土体的三维分析 本文将首先讨论体积闭锁及其原因 接着对8节点等参元 Wilson非 协调元 EAS单元和14节点单元的闭锁性态进行分析 最后用方形基础的承载力计算来考 收稿日期 1999-04-29 修订日期 1999-05-30 基金项目 国家自然科学基金资助项目 作者简介 杨军(1974) 男 四川邻水人 博士 从事饱和土动力数值分析 转载 中国科技论文在线 土体三维弹塑性分析中的体积闭锁分析及单元选择 7 察各种单元的性态 从而对土体三维分析中选择有效可靠的单元提供依据 2 体积闭锁及其原因概述 为说明闭锁产生的原因 我们来看用两个三节点三角形单元分析土体变形的一个简单 例子(图1) 假设为平面应变问题 很明显 任何使材料不可压缩的本构关系 如不可压缩 弹性或Von Mises完全塑性 均要求这两个三角形单元的面积保持不变 当把这种运动约束 加在右下单元时 自由节点A 的位移被约束在水平线上 同理 左上单元约束