abaqus子程序用fortran语言编写，且子程序又通过abaqus调用，因而涉及到很多的力学及有限元理论知识，使得调试难度陡然增加。那么abaqus子程序常用的调试方法是什么呢？其实说出来也很简单，就是打印变量法，即把程序中的变量打印出来（一般打印结果在当前工作目录下的job-name.log的文件里），然后考察其合理性。
如在上述umat子程序中定义了一个save在内存的变量，用来记录进入umat的次数。
直接打印变量的方法很有用，但是当单元数增加后，众多的打印信息形成很多的干扰因素，这时候我们需要控制打印的频率，常用的方法为针对某个单元的某个积分点打印，如下面的一段程序所示，红色圆圈里内容的含义是当单元编号为1且积分点编号为1时才打印相关信息，这样调试更具针对性。
我们通过打印信息发现，一次迭代二次进入umat，第一次进入umat仅是给abaqus返回雅可比矩阵，第二次进入umat的目的是为了更新应力等信息。
视频教程有这个帖子的更详细解释，感兴趣的点击下面的链接观看
微信
QQ
微博