一起学习UMAT课件.pdfVIP

一起学习UMAT ZHANG chunyu 1、什么时候用用户定义材料（User-defined material, UMAT）？ 很简单，当ABAQUS 没有提供我们需要的材料模型时。所以，在决定自己定义一 种新的材料模型之前，最好对 ABAQUS 已经提供的模型心中有数，并且尽量使用现 有的模型，因为这些模型已经经过详细的验证，并被广泛接受。 2、好学吗？需要哪些基础知识？ 先看一下ABAQUS 手册（ABAQUS Analysis Users Manual）里的一段话： Warning: The use of this option generally requires considerable expertise. The user is cautioned that the implementation of any realistic constitutive model requires extensive development and testing. Initial testing on a single element model with prescribed traction loading is strongly recommended. 但这并不意味着非力学专业，或 力学基础知识不很丰富 （就如我本人J） 就只能望洋兴叹，因为我们的任务不是开发一套完整的有限元软件，而只是提供一 个描述材料力学性能的本构方程 （Constitutive equation）而已。当然，最基本 的一些概念和知识还是要具备的，比如 应力(stress),应变 （strain）及其分量； volumetric part 和 deviatoric part ；模量（modulus ） 泊松比 (Poisson’s ratio) 拉美常数 (Lame constant)；矩阵的加减乘除甚至求逆；还有一些高等数学知识如积分 微分等。 3 UMAT 的基本任务？ 我们知道，有限元计算 （增量方法）的基本问题是： 已知第 n 步的结果（应力，应变等） ，ε ； 然后给出一个应变增量dε , n n n 1 计算新的应力 n1 。UMAT 要完成这一计算，并要计算Jacobian 矩阵DDSDDE(I,J) = 。 是应力增量矩阵 （张量或许更合适）， 是应变增量矩阵。 DDSDDE(I,J) 定义了第 J 个应变分量的微小变化对第 I 个应力分量带来的变化。 该矩阵只影响收敛速度，不影响计算结果的准确性 （当然，不收敛自然得不到结 果）。 4、UMAT 到底起什么作用？为什么Jacobian 矩阵不影响计算结果的准确性却影响收 敛速度？ 有限元计算的中心问题就是求得节点的位移 （进而应变、应力），以使内力和 外力达到平衡： Fint d Fext (1) d 是节点位移矩阵，黑体字表示矩阵或矢量。除了小变形、线弹性问题，方程 （1）是非 性的，要用迭代的方法解出： K d Fext Fint di (2) T n 1 n 1 di 1 di d (3) n1 n1 i 表示一个增量步内的第 i 次迭代，n 表示第 n 个增量步。KT 是切线刚度，由材料 的 Jacobian 矩阵结合单元计算组装而得。刚度矩阵其实就是力对位移的梯度。要 想快速收敛，位移增量应沿该梯度方向变化，也就是说，如果 Jacobian 矩阵不是 那么准确，自然KT 也不怎么准确，那么满足 （1）式的位移被找到的速度也就变 慢，甚至发散，根本找不到。但收敛速度无论