Abaqus 子程序开发日志03-umat之Orthogonally Anisotropic

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

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

由矩阵可以知道，正交各项异性材料仅需9个独立常数，只要知道这9个刚度系数，就可以构建出材料的本构方程。不懂去翻材料力学相关书籍。

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

材料属性就随便定了，详见下图umat说明。

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

分别用自带的正交各向异性与umat子程序进行计算。

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

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

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

D1111,D1122,D2222,D1133,D2233,D3333,D1212,D1313,D2323

(ABAQUS自带对照组的材料属性也对照此输入)

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

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

到此结束。cae文件如下：

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

正交各向异性umat：

共享库：