非线性有限元下册.pdfVIP

非 线性 有 限单 元 法 下册 （程序部分） 任青文 编著 河海大学工程力学系 第五章 弹塑性有限元程序设计 本章将给出一个二维的弹塑性有限元程序 NFA2D ，它可以解决平面应力、平面应变和轴 对称问题。鉴于不同类型的材料显示出不同的弹塑性性质，程序中有四种屈服准则可供选择， 其中 Tresca 和 Von Mises 准则适用于金属材料，Mohr-Coulomb 和 Drucker-Prager 准则适用于 混凝土、岩石和土。 5.1 程序中的弹塑性模型 5.1.1 屈服准则 屈服准则用来确定塑性变形发生时的应力水平，第三章已经给出其一般形式及 Mises 模型和 Drucker—Prager 模型屈服准则的具体表达式。本程序共有四类屈服准则可供使用。 Tresca 准则 该准则认为最大剪应力达到一定数值时材料产生塑性屈服。若一点的主应力为 σ1、σ 和 2 σ ，且 σ≥σ≥σ ，则其表达式为 3 1 2 3 f =( - )- 0 (5.1) 1 3 s 如果不知道 σ 、σ 和 σ 的大小，则上式可改写为 1 2 3  2 2 2 2 2 2  f          0 (5.2)  1 2 s 2 3 s 3 1 s  式中，σ 为单向应力状态下，材料的屈服应力，它是硬化参数 的函数。式（5.2 ）表示屈服面 s 为应力空间内的正六棱柱（图 5.1 ），它的六个面分别平行于三个主应力。它在平面上的投 影是一个正六边形（图 5.2a ）。由于这种模型与静水压力无关，屈服函数只与应力偏量的第 二不变量 J 和第三不变量 J 有关，因而可表示为f (σ -σ ，σ -σ)，其二维的图像如图 5.2 （b ） 2 3 1 3 2 3 所示。 σ1 = σ2 = σ3 Tresca Mises σ3 σ2 σ1 π平面 σ1 + σ2 +σ3 0