非线性有限元及弹塑性力学讲解课件_1.ppt

2000.3 哈尔滨建筑大学 王焕定教授制作 第二章 弹塑性本构关系简介 1. 弹性介质本构关系 2. 弹塑性力学有关内容简介 3. 几种常用弹塑性材料模型简介 4. 弹塑性矩阵的建立步骤 1. 弹性介质本构关系 2. 弹塑性力学有关内容简介 3. 几种常用弹塑性材料模型简介 4. 弹塑性矩阵的建立步骤 3.1 等向强化-软化的米塞斯（Mises）材料 由薄壁圆筒的实验研究可得，这种材料的屈服面方程为 在主应力状态下，第二不变量为 式中J2是应力偏张量的第二不变量 ，由于偏张量第一不变量等于零 ，因此 在单向拉伸状态下，J2=σ2/3。在纯剪状态下，J2=τ2。一般情况下，sij=σij -σkkδij /3,所以 . 屈服面式中χ(k)，是由单向应力状态的数据确定的屈服参数。在单向拉伸时为χ2=σB2/3。在纯剪状态下χ=τB。任何情况下χ都是硬化参数塑性功wp的函数。 根据屈服面表达式，可求得 因为 因此 为了求A，需先由屈服面对wp的偏导数求M 式中Gp是 曲线的斜率。同理，对单向拉伸情况，-M=Ep/3。 由此可得A=G+Gp（或A=G+Ep/3），又因 因为 纯剪 Gp是塑性剪切模量 单向拉伸 Ep是塑性拉伸模量 由此可得 弹塑性矩阵为 加、卸载准则为 统一的本构关系为 3.2 随动强化的米塞斯材料 这种材料的屈服面方程为 由此屈服面方程出发，求导可得 式中α是与内变量有关的量， 称为应力迁移张量，由它可以确定屈服面在应力空间的位置。 是一个常量，表示屈服面形状不变。 如上推导即可得到应变空间的本构关系和流动准则。米塞斯屈服准则主要适用于金属材料。 3.3 岩土工程广泛使用的Mohr-Coulomb材料 这种材料的屈服准则为 式中τ是破坏面上切应力，σn是破坏面的正压力。c是材料粘性系数，φ是内摩擦角，c和φ是两个材料常数。在主应力空间，此屈服准则可用下图示意 * 对线弹性介质只有两个独立的弹性常数，但应力应变（本构）关系有多种表示形式： 用G和μ表示 用G和体积模量K表示 1.1 线性弹性小变形 式中应力和应变偏张量分别为 如果用拉梅（Lame）常数表示，则有 弹性常数间有如下关系 利用上述关系，只要已知两个弹性常数就可写出有限元分析中的弹性矩阵(D)。 例如，当以G和μ表示时，以张量形式表示的本构关系为 由此可获得弹性张量Dijkl。其他可仿此写出。 非线性弹性介质的本构关系，一般是根据材料的力学试验通过拟合来得到的。例如金属材料单向拉伸Romberg-Osgood模型的关系为 式中k和n为拟合的实验参数，E为初始弹性模量。一般情况下本构关系可表为 1.2 非线性弹性小变形 在有限元分析中有两种应用形式：全量形和增量形本构关系。 全量本构关系的表达形式和线性弹性情况相同，也即 1.2.1 全量形式本构关系 但其中的弹性系数Gs,μs不再是常数，它们是应变或应力的函数，分别称为割线弹性系数。可将它们看作与一定应力（或应变）水平对应的割线常数（割线剪切模量和割线泊松比）。 式中 为割线弹性张量，形式上它仍可表为 例如对混凝土，Andenaes等依据实验给出，八面体正应力、切应力和八面体线应变、角应变间关系为 σoct εoct Ks Kt 并有 其中G、K分别为初始切线剪切和体积模量， 为混凝土单轴抗压强度，a、m、c和p为由试验确定的常数。 Gt Gs 1.2.2 增量形式本构关系 增量本构关系的表达形式为 但其中的弹性系数Gt,μt也不是常数，也是应变或应力的函数，分别称为切线弹性系数。可将它们看作与一定应力（或应变）水平对应的切线常数（切线剪切模量和切线泊松比）。 式中 为切线弹性张量，形式上仍可表为 上面介绍的是哥西方法，讲义上还简述了格林方法，大家可自行阅读。 韧性（塑性）金属材料单向拉伸试验曲线如下图示意 2.1 应力空间表述的弹塑性本构关系 弹性极限 屈服下限 屈服上限 强度极限 强化段 软化段 弹性变形 残余变形 卸载 包辛格效应 反向屈服点 卸载、反向加载 由单向拉伸曲线可见，弹塑性材料受外部作用的反应和变形的历史有关（可称为历史相关性或路径相关性），因此本构关系应写成增量关系。又因弹塑性状态下加载和卸载有不同的规律，所以其本构关系的表述要比非线性弹性情况复杂。 以应力为坐标，其每点代表一个应力状态，如此的空间称为应力空间。 判断材料处于弹性还是塑性的