清华大学计算固体力学第五次_本构模型详解.ppt

6 多轴塑性 Tresca屈服准则 Mises屈服准则 在有限元程序中一般应用哪种屈服准则？为什么？ 摩擦滑移屈服表面 6 多轴塑性 Mohr-Coulomb本构模型 滑移方向（塑性流动）是水平的（沿Q 的方向）而不是垂直屈服面。这是非关联塑性流动的例子。对于连续体和多轴应力－应变状态的行为，M-C准则具有普适性。它应用于模拟土壤和岩石。 M-C准则是基于这样的概念，即当任意面上的切应力和平均法向应力达到临界组合时在材料中发生屈服 c是内聚力，通过 定义内摩擦角 6 多轴塑性 Mohr-Coulomb屈服行为 Mohr-Coulomb屈服表面 Drucker-Prager屈服表面 在Mohr平面上的两条直线代表了方程式，它们是Mohr圆的包络 并称为Mohr破坏或者失效包络。假设主应力 应力状态 屈服准则 6 多轴塑性 考虑 的特殊情况并让 ，代表剪切屈服强度， 上式成为 即为Tresca准则。 在Tresca和M-C屈服表面上的直线线段便于塑性问题的解析处理。 然而，从计算的观点看，夹角使得本构方程难以建立(例如，计算屈服面的法线)。通过改进von Mises屈服准则结合压力的影响，Drucker-Prager屈服准则避免了与夹角有关的问题： 这是一个光滑圆锥的方程， 为等效Cauchy应力，选择常数有 D-P屈服表面通过了M-C屈服表面上的内部或者外部顶点（取加号对应于内部顶点，而取减号对应于外部顶点）。 9 应力更新算法 本构方程率形式的积分算法称为应力更新算法（也称为本构更新算法），包括： 径向返回算法的一类图形返回算法， 算法模量与基本应力更新方案一致的概念， 大变形问题的增量客观应力更新方案， 基于弹性响应的应力更新方案，即自动满足客观性的超弹性势能。 给出描述本构模型的某些其它连续介质力学观点， 展示Eulerian，Lagrangian和两点拉伸的概念， 描述后拉、前推和Lie导数的运算， 材料框架客观性，材料的对称性，以本构行为的张量表示讨论了不变性的某些方面， 讨论由于热力学第二定律和某些附加的稳定性必要条件对材料行为的约束。 9 应力更新算法 对于积分率本构方程的数值算法称为本构积分算法或者应力更新算法。对于率无关和率相关材料提供了本构积分算法。 讨论简单的小应变塑性，将小应变算法扩展至大变形，将大变形分析的积分算法保持在基于本构方程客观性的基础上。 展示了关于大变形塑性的逐步客观积分算法。 讨论关于大变形超弹－塑性材料的应力更新算法，回避对应力率方程的积分。 描述了与本构积分算法相关的计算模量，采用隐式求解算法发展材料的切线刚度矩阵。 率无关塑性的图形返回算法 9 应力更新算法 小应变、率无关弹－塑性的本构方程 应力－应变反应与变形率无关的一种材料称为率无关；否则为率相关。 ， Kuhn-Tucker条件，上面第一个条件表明塑性率参数是非负的， 第二个条件表明当塑性加载时，应力状态必须位于或限制在塑性表面上， 最后条件也可以作为由已知一致性条件 的率形式。 塑性流动方向经常特指为 ，这里 称为塑性流动势 屈服条件 是标量塑性流动率， 是塑性流动方向 h 塑性模量 q 内变量 ）应力状态必须保持在屈服面 因此 。对于弹性加载或者卸载 ，没有塑性流动。 对于塑性加载（ 率无关塑性的图形返回算法 9 应力更新算法 上， 在时刻n 给出一组 和应变增量 本构积分算法的目的是计算 并满足加－卸载条件 在 时刻的应力给出为 求解的一致性条件给出 设想能够应用这个塑性参数值以提供更新的应力率、塑性应变率和内变量率，并且写出简单的向前Euler积分公式算法 率无关塑性的图形返回算法 9 应力更新算法 但在下一步，这些应力和内变量的更新值并不满足屈服条件，所以 由于解答从屈服表面漂移，常常导致不精确的结果，因此不受人青睐。公式也称为切线模量更新算法，形成了计算率无关塑性早期工作的基础。 率无关塑性的图形返回算法 9 应力更新算法 这导致考虑另外一些方法进行率本构方程的积分，目的之一是强化在时间步结束时的一致性，例如， 为避免离开屈服面的漂移。有许多不同的积分本构算法，这里主要关注一类方法－－返回图形算法，它是强健和精确的，被广泛应用。著名的von Mises塑性径向返回方法是返回图形算法的特例。 返回图形算法包括： 一个初始的弹性预测步，包含（在应力空间）对屈服表面的偏离， 以及塑性调整步使应力返回到更新后的屈服表面。