P02_Chap_01_材料本构模型.doc

第二部分 材料库midas FEA 第二部分 材料库 midas FEA 第一章 | 材料本构模型midas 第一章 | 材料本构模型 midas FEA We Analyze and Design the Future PAGE We Analyze and Design the Future PAGE 208 PAGE 207 PAGE 207 We Analyze and Design the Future 第一章 材料本构模型1-1 概要 第一章 材料本构模型 1-1- 材料的塑性效应与弹性效应不同，即使卸去荷载仍然会存在永久变形。为了描述这种特性，将应变区分为弹性应变和塑性应变两个部分来考虑。 (1.1.1 其中： ：总应变 ：弹性应变 (elastic strain) ：塑性应变 (plastic strain) 虎克定律(Hook’s law)可定义弹性范围内的应变与应力的关系，将其适用于式(1.1.1) (1.1.2 其中： : 应力向量 : 材料刚度矩阵(material stiffness matrix) 荷载的结构内任意一点的应力大小可作为该点处于弹性状态还是塑性状态的判断标准。这种判断标准对于钢材或者混凝土，根据材料的特性会有不同的定义，我们把它叫做屈服标准。材料的屈服标准是通过对不同应力状态的大量的构件的实验来决定的，导致发生塑性流的应力值可按应力空间形式使用下列函数来表示。 (1.1.3 其中： : 屈服函数 (yield function) : 硬化参数 (hardening parameter) 屈服函数比0小时不发生塑性流动，比0大时发生塑性流动。 1-1-2 塑性流动法则( 材料的屈服会导致塑性流动发生，为了保持材料的平衡状态，这种塑性流动会使应力重新得到分配。塑性流动的计算表现为非线性，因此需要使用增量的方式表述。在计算塑性流动的各项参数中主要参数有应力增量方向参数和塑性应变增量的方向参数，其中应力增量方向如下： (1.1.4) 其中：是沿应力空间屈服面的法线方向，是指屈服函数的个数。 塑性应变增量根据Koiter’s law将大小和方向分开表示如下： (1.1.5 其中：为塑性势函数(plastic potential function)，一般是应力和硬化参数的函数，通过材料实验能够得到。为塑性乘数(plastic multiplier)，它的具体数值和材料硬化法则有关，需满足如下的Kuhn-Tucker条件。 ，， (1.1.6) 根据上述条件，屈服函数比0小时，总是为0，不会发生塑性流动。 式(1.1.5)中的是描述塑性增量方向的向量。此时可以对塑性势函数取和屈服函数相同的形式，即用来定义塑性增量的方向，称之为关联塑性流动法则 (associated flow rule)；使用塑性势函数的来定义塑性增量方向的方法称为非关联塑性流动法则(non-associated flow rule)。对于von Mises或者Tresca等本构，其应力空间上沿着等压轴(hydrostatic axis)有恒定的偏差应力分布，对于这种材料本构一般使用关联塑性流动法则。 对于Mohr-Coulomb和Drucker-Prager等本构其应力空间上沿着等压轴有变化的偏差应力分布，对于这种材料本构则需要使用非关联塑性流动法则。此时由于其应力增量方向和应变增量方向不一致，可以起到抑制过分的体积膨胀的效果。不过因为刚度矩阵为非对称，所以需要使用非对称求解器求解，导致其收敛速度会非常慢。因此一般还是使用关联塑性流动法则。 但是当混凝土受箍筋、钢管或者钢箱的约束效果很大时，建议使用非关联塑性流动法则。 1-1-3 在弹塑性模型中定义硬化特性的方法有应变硬化假设和塑性功硬化假设两种。应变硬化假设是假设硬化过程随等效塑性应变的增加而进展，塑性功硬化假设是假设硬化过程随塑性功的增加而进展。 应变硬化中的硬化参数使用的是无量刚的等效塑性应变，可按下式定义。 (1.1.7 其中： , 此时使用主塑性应变来计算硬化参数很便利，因此较为普遍地在使用，其计算式如下： (1.1.8 其中： 塑性功硬化是将应力单位的塑性功(plastic work)作为硬化参数来使用，其硬化行为表述如下： (1.1.9 上述两种形式的硬化都是对时间进行积分，计算到当前荷载步骤的硬化参数。 (1.1.10 利用硬化参数