线性强化弹塑性umat子程序-umat在abaqus计算流程中的意义及调试方法

abaqus子程序用fortran语言编写，且子程序又通过abaqus调用，因而涉及到很多的力学及有限元理论知识，使得调试难度陡然增加。那么abaqus子程序常用的调试方法是什么呢？其实说出来也很简单，就是打印变量法，即把程序中的变量打印出来（一般打印结果在当前工作目录下的job-name.log的文件里），然后考察其合理性。

如在上述umat子程序中定义了一个save在内存的变量，用来记录进入umat的次数。

直接打印变量的方法很有用，但是当单元数增加后，众多的打印信息形成很多的干扰因素，这时候我们需要控制打印的频率，常用的方法为针对某个单元的某个积分点打印，如下面的一段程序所示，红色圆圈里内容的含义是当单元编号为1且积分点编号为1时才打印相关信息，这样调试更具针对性。

我们通过打印信息发现，一次迭代二次进入umat，第一次进入umat仅是给abaqus返回雅可比矩阵，第二次进入umat的目的是为了更新应力等信息。

视频教程有这个帖子的更详细解释，感兴趣的点击下面的链接观看

http://www.jishulink.com/college/video/c14014