Abaqus 子程序开发日志02-umat之Engineering Constants

umat子程序只用于通用静力分析步，使用用户自定义材料，这里对正交各向异性材料本构模型进行构建（使用工程常数）。

对于正交各向异性的刚度矩阵为：

由矩阵可以知道，正交各项异性材料仅需9个独立常数，用各项异性材料的9个工程常数进行表示，进而可以构建出材料的本构方程。不懂去翻材料力学相关书籍。

这里就随便建立一个长宽高为5；5；20的立方体对其进行施加一个单位为5的拉伸位移（较长方向拉伸）。

材料属性就随便找了一篇论文，按照上面的材料属性进行赋值（子程序已分享，可以自己尝试一些其他的工程常数）：

单元类型全部采用相同的完全积分单元。材料方向按全局坐标算（各向同性材料不需要指定方向，各向异性需要）

分别用自带的工程常数材料属性与umat子程序进行计算。

自带的材料属性的设定以及模型的建立、边界条件的设定等等不再说明，只对umat的使用进行说明。

此次直接使用for文件（二进制文件也一样），在材料属性的输入过程，在通用中找到用户自定义材料，输入设定的参数。如图所示：

这里要说明一下输入的材料参数，自上而下分别为：

E1,E2,E3,Nu23,Nu13,Nu12,G23,G13,G12

点击ok就好，其余的和平时操作一致，只需在最后提交job时提供用户子程序。

然后提交计算就好。结果如图所示：

到此结束。cae文件如下：

若选用缩减积分单元时要开启沙漏控制。

正交各向异性umat：

共享库：