[2018年最新整理]5弹塑性本构关系简介.pptVIP

本构关系 有限元(二)的具体内容 材料(非线性)本构关系 固体力学大变形基础 非线性方程组的解法 材料非线性有限元分析 大变形有限元分析 边界元法－流体计算 弹塑性本构关系简介 1 弹塑性力学有关内容简介 2 几种常用弹塑性材料模型简介 3 弹塑性矩阵的建立步骤 固体力学大变形基本知识 1. 物体运动的物质描述 2. 欧拉、拉格朗日和克希荷夫应力 3. 大变形时平衡方程和虚位移原理 4. 大变形本构关系 非线性代数方程组的数值解法 1 直接迭代法 2 牛顿法和修正牛顿法 3 拟牛顿法 4 增量方法 5 增量弧长法 材料非线性有限元分析 弹塑性问题的有限单元法 大变形问题的有限单元法 1. 弹性大变形问题的有限元法 2. 物质描述大变形增量问题的T.L 、U.L法 弹塑性本构关系 应变空间 应力空间 例：等向强化-软化的米塞斯材料弹塑性矩阵的确定 在主应力状态下，第二不变量为 由薄壁圆筒的实验研究可得，这种材料的屈服面方程为 式中J2是应力偏张量的第二不变量 在单向拉伸状态下，J2=σ2/3。在纯剪状态下，J2=τ2。一般情况下，sij=σij -σkkδij /3,所以 . 屈服面式中χ(k)，是由单向应力状态的数据确定的屈服参数。在单向拉伸时为χ2=σB2/3。在纯剪状态下χ=τB。任何情况下χ都是硬化参数塑性功wp的函数。 根据屈服面表达式，可求得 因为 为了求A，需先由屈服面对wp的偏导数求M 式中Gp是 曲线的斜率。同理，对单向拉伸情况，-M=Ep/3。 纯剪 Gp是塑性剪切模量 单向拉伸 Ep是塑性拉伸模量 因为 由此可得A=G+Gp（或A=G+Ep/3），又因 由此可得弹塑性矩阵为 加、卸载准则为 统一的本构关系为 1) 根据实验研究建立材料合理的屈服条件f。 2) 由屈服条件求屈服面方程对应力或不变量等的导数。 等向强化-软化材料的屈服面方程为 3) 根据硬化参数的选取，计算M。 4) 由 计算A。 5) 由 计算塑性矩阵。 6) 计算 并由此判断H(l)。 7) 最后形成弹塑性矩阵 。 * 有限单元法II ——2005级硕士生课程 同济大学土木工程桥梁工程系 同济大学土木学院桥梁工程系 对线弹性介质在小变形情况下只有两个独立的弹性常数，但应力应变（本构）关系有多种表示形式： 用G和μ表示 用G和体积模量K表示 式中应力和应变偏张量分别为 如果用拉梅（Lame）常数表示，则有 弹性常数间有如下关系 利用上述关系，只要已知两个弹性常数就可写出有限元分析中的弹性矩阵(D)。 例如，当以G和μ表示时，以张量形式表示的本构关系为 由此可获得弹性张量Dijkl。其他可仿此写出。 弹性张量Dijkl 一般情况下本构关系可表为 在有限元分析中有两种应用形式：全量形和增量形本构关系。 一般是根据材料的力学试验通过拟合来得到的。 全量本构关系的表达形式和线性弹性情况相同，也即 全量形式本构关系 但其中的弹性系数Gs,μs不再是常数，它们是应变或应力的函数，分别称为割线弹性系数。可将它们看作与一定应力（或应变）水平对应的割线常数（割线剪切模量和割线泊松比）。 式中 为割线弹性张量，形式上它仍可表为 本构关系 例如对混凝土，Andenaes等依据实验给出，八面体正应力、切应力和八面体线应变、角应变间关系为 σoct εoct Ks Kt 并有 其中G、K分别为初始切线剪切和体积模量， 为混凝土单轴抗压强度，a、m、c和p为由试验确定的常数。 Gt Gs 增量形式本构关系 增量本构关系的表达形式为 但其中的弹性系数Gt,μt也不是常数，也是应变或应力的函数，分别称为切线弹性系数。可将它们看作与一定应力（或应变）水平对应的切线常数（切线剪切模量和切线泊松比）。 式中 为切线弹性张量，形式上仍可表为 塑性材料受外部作用的反应和变形的历史有关（可称为历史相关性或路径相关性），本构关系应写成增量关系。 韧性（塑性）金属材料单向拉伸试验曲线如下图示意 应力空间表述的弹塑性本构关系 弹性极限 屈服下限 屈服上限 强度极限 强化段 软化段 弹性变形 残余变形 卸载 包辛格效应 反向屈服点 卸载、反向加载 判断材料处于弹性还是塑性的准则，称为屈服条件或塑性条件。弹性和塑性区的分界面称为屈服面。 1) 屈服准