空间8节点等参单元计算Drucker—Prager理想塑性材料.pdfVIP

维普资讯 天津城市建设学院学报 第 8卷 第 3期 2002年 9月 JournalofTianjinInstituteofUrbanConstruction Vo1．8 No．3 Sep．2002 空间8节点等参单元计算 Drucker—Prager理想塑性材料 唐柏鉴，侯朝胜 (天津大学 建筑工程学 院 ，天津 300072) 摘要 ：DruckerPrager模 型广泛应 用于混凝土 、岩石等材料 ，采用空 间 8节 点等参单元对其进 行增量迭代计算 ，每增量步 内NewtonRaphson法迭代 2次，4阶高斯积分 ，算例 结果与理论 解 吻合很 好 ． 关 键 词 ：有限元 ；非线性方程 ；Drucker—Prager模 型 中图分类号 ：TB301 文献标识码 ：A 文章编号 ：1006—6853(2002)03一O16904 1 Drucker—Prager理想塑性材料 2 有 限元法 Drucker—Prager模型 单 刚 ： 厂({ )一 aI + √J2一 K 一 0 (1) [七] JjJEB][]·EB]·IJI·ddd (6) 该模 型考虑体积力对 屈服 的影 响 ，因此广泛用于 体积力等效节点荷载 ： 混凝土 、岩石等材料 ；它改进 了Mohr—Coulomb破坏 曲 {Fo一…l-N] ·g·IJl·d钮 (7) 面 ，以六棱锥为 圆锥 ，因此便于数值分析 ． 内力列 阵 ： 由拉 、压单轴应力下 的破坏状态 ，求得 r 二 {o= …[曰]{ ·IJl·d钮 (8) a ： = — — ，，AK ： —===———一 L(2) v／3(厂 + f) v／3(厂 + f) 弹性 阶段式 (6)、(7)采用 1阶高斯积分 ；塑性 阶 理想塑性材料一致性条件 段式 (6)、(7)、(8)都采用 3阶高斯积分 ． {d ： 0 (3) 其 中．[ ]取决于单元应力状态 ，为 [。]=rED]+ (1一 r)[D ] 0≤ r≤ 1 由式 (1)、(3)可推 出Drucker—Prager材料的本构矩 阵 (9) [D ]= ED]一 [D] (4) 所谓 “单元应力 ”，笔者有两种处理手法 ：① 取 ( [D] l9K2a2{}{}+ = 0， ： 0， = 0)中心点代表单元应力