5.本构模型-应力更新专题-UMAT和VUMAT分解.pptx

计算固体力学应力更新专题 柳占立 庄茁 liuzhanli@tsinghua.edu.cn 2017年4月3日 导出有限元方程（完全，更新、共旋） 不同用户要求： 知道有限元基本流程； 编写本构关系、节点内力 开发有限元程序 第4-6章 求解有限元方程（隐式和显式） 有限元基本知识 有限元高级知识 单元技术、结构单元、接触等 第7-10章 应力更新 什么是应力更新 几种客观率 前推后拉及Lie导数 常用客观率应力更新流程 ABAQUS用户自定义材料 理论 应用 有限元程序 赋予力学以活力 1.应力更新 ？ 应力增量 2.几种客观率 Truesdell率： 前推 Lie导数 t= Js 2.几种客观率 Green-Naghdi率： 前推 Lie导数 共轴旋转应力率 2.几种客观率 ？ Jaumann率： W = W Rigid body rotation or The principal directions of the strain remaining constant 如何得到 2.几种客观率 Truesdell率： Green-Naghdi率： Jaumann率： 应力更新: 本构关系 例4.1 考虑处于剪切状态的一个单元，如图所示。对于次弹性各向同性材料，应用Jaumann，Truesdell和Green－Naghdi率求出剪切应力 x(t) = X + k t Y , y(t) = Y , z(t) = Z 2.几种客观率 对于不同客观率采用了相同的材料常数，其差别是非常大的。事实上，这是误用了材料模型。材料模型必须根据不同的率转换。这是变形体，若是刚体转动，Jaumann率与Truesdell率是一致的。 例4.1 2.几种客观率 Green-Naghdi率 Jaumann率 ？ 发生有限剪切时 慎用Jaumann率 切线模量之间的关系 对于各向同性材料Jaumann率的切线模量为 次弹性本构关系共同应用的形式为 对于同一种材料，切线模量不同，材料反应的率形式不同，如 如果 是常数， 不是常数。 切线模量 证明见第5.4.5节，推导复杂 2.几种客观率的关系 如何得到正确的结果？ 2.几种客观率的关系 如何得到正确的结果？ 假设CSE已知： 如何处理各向异性材料？ 2.几种客观率的关系 如何得到正确的结果？ 如何处理各向异性材料？ PK2和共轴旋转应力 横观各向同性 2.几种客观率的关系 什么时候可以不区分几种客观率？ Assume F ≈ R, ignoring the stretch component of F Assume W = W ？ 退化 简化 Comparison of different objective stress rate Truesdell Difficult to implement Not used in commercial software Green-Naghdi Kinematically consistent with the rate of Cauchy stress Must accurately determine the rotation tensor R Jaumann Relatively easy to implement Produces symmetric tangent moduli 2.几种客观率的关系 总结一： 2.几种客观率的关系 总结二： 当参考构型和当前构型重合时： F=I 当旋转构型和当前构型重合时： R=I 使TL和UL的互相转换变得非常方便！ Solver Element Type Objective Rate Abaqus/Standard Solid (Continuum) Jaumann Structural (Shells, Membranes, Beams, Trusses) Green-Naghdi Abaqus/Explicit Solid (Continuum) Jaumann Green-Naghdi Structural (Shells, Membranes, Beams, Trusses) Green-Naghdi 2.几种客观率的关系 总结三： 后拉、前推 Euler张量 Lagrangian张量 两点张量 Green应变张量E PK2应力 线单元dX 线单元dx Cauchy应力 速度梯度 L＝D＋W Lagrangian矢量dX和Eulerian矢量dx定义的二阶张量 可以由后拉和前推运算给出E-L张量之间映射的统一描述。 例如，L矢量dX由F前推到当前构形给出E矢量dx E矢