应力更新算法概念.ppt

9 应力更新算法 本构方程率形式的积分算法称为应力更新算法（也称为本构更新算法），包括： 径向返回算法的一类图形返回算法， 算法模量与基本应力更新方案一致的概念， 大变形问题的增量客观应力更新方案， 基于弹性响应的应力更新方案，即自动满足客观性的超弹性势能。 给出描述本构模型的某些其它连续介质力学观点， 展示Eulerian，Lagrangian和两点拉伸的概念， 描述后拉、前推和Lie导数的运算， 材料框架客观性，材料的对称性，以本构行为的张量表示讨论了不变性的某些方面， 讨论由于热力学第二定律和某些附加的稳定性必要条件对材料行为的约束。 9 应力更新算法 对于积分率本构方程的数值算法称为本构积分算法或者应力更新算法。对于率无关和率相关材料提供了本构积分算法。 讨论简单的小应变塑性，将小应变算法扩展至大变形，将大变形分析的积分算法保持在基于本构方程客观性的基础上。 展示了关于大变形塑性的逐步客观积分算法。 讨论关于大变形超弹－塑性材料的应力更新算法，回避对应力率方程的积分。 描述了与本构积分算法相关的计算模量，采用隐式求解算法发展材料的切线刚度矩阵。 率无关塑性的图形返回算法 9 应力更新算法 小应变、率无关弹－塑性的本构方程 应力－应变反应与变形率无关的一种材料称为率无关；否则为率相关。 ， Kuhn-Tucker条件，上面第一个条件表明塑性率参数是非负的， 第二个条件表明当塑性加载时，应力状态必须位于或限制在塑性表面上， 最后条件也可以作为由已知一致性条件 的率形式。 塑性流动方向经常特指为 ，这里 称为塑性流动势 ）应力状态必须保持在屈服面 因此 。对于弹性加载或者卸载 ，没有塑性流动。 对于塑性加载（ 率无关塑性的图形返回算法 9 应力更新算法 上， 在时刻n 给出一组 和应变增量 本构积分算法的目的是计算 并满足加－卸载条件 在 时刻的应力给出为 求解的一致性条件给出 设想能够应用这个塑性参数值以提供更新的应力率、塑性应变率和内变量率，并且写出简单的向前Euler积分公式算法 率无关塑性的图形返回算法 9 应力更新算法 但在下一步，这些应力和内变量的更新值并不满足屈服条件，所以 由于解答从屈服表面漂移，常常导致不精确的结果，因此不受人青睐。公式也称为切线模量更新算法，形成了计算率无关塑性早期工作的基础。 率无关塑性的图形返回算法 9 应力更新算法 这导致考虑另外一些方法进行率本构方程的积分，目的之一是强化在时间步结束时的一致性，例如， 为避免离开屈服面的漂移。有许多不同的积分本构算法，这里主要关注一类方法－－返回图形算法，它是强健和精确的，被广泛应用。著名的von Mises塑性径向返回方法是返回图形算法的特例。 返回图形算法包括： 一个初始的弹性预测步，包含（在应力空间）对屈服表面的偏离， 以及塑性调整步使应力返回到更新后的屈服表面。 方法的两个组成部分是： 一个积分算法，它将一组本构方程转换为一组非线性代数方程， 一个对非线性代数方程的求解算法，该方法可基于不同的积分算法， 例如生成梯形法则，生成中点法则或者Runge-Kutta方法。 基于向后Euler算法，考虑一个完全隐式方法和一个半隐式方法。 完全隐式的图形返回算法 9 应力更新算法 在完全隐式的向后Euler方法中，在步骤结束时计算塑性应变和内变量的增量，同时强化屈服条件，这样，积分算法写成为 公式是一组关于求解 的非线性代数方程。注意到更新变量来自前一个时间步骤结束时的收敛值，这就避免了非物理意义的效果，例如当用不收敛的塑性应变和内变量值求解路径相关塑性方程时可能发生的伪卸载。 在时刻n 给出一组 和应变增量 通过方程系统的解答获得了应变 在时刻n ＋1， 完全隐式的图形返回算法 9 应力更新算法 如果解答过程是隐式的，可以理解应变 是在隐式解答算法的最后迭代后的总体应变。 塑性应变增量给出为 代入表达式 关联塑性的最近点投射方法 是弹性预测的试应力 是塑性修正量，它沿着一个方向，即规定为在结束点处塑性流动的方向， 返回或者投射试应力到适当更新的屈服表面（考虑硬化）。 而数值 完全隐式的图形返回算法 9 应力更新算法 由总体应变的增量驱动弹性预测状态，而由塑性参数的增量 驱动塑性修正状态。因此，在弹性预测阶段，塑性应变和内变量保持固定，而当塑性修正阶段，总体应变是不变的。在弹性预测阶段，由公式得到的结果为 关联塑性的最近点投射方法 其中 完全隐式的图形返回算法 9 应力更新算法 非线性代数方程组解答一般由Newton过程求解。基