Abaqus 子程序开发日志04-umat之Transversely Isotropic

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

对于横观各向同性的刚度矩阵为：

由矩阵可以知道，横观各向同性材料仅需5个独立常数，用横观各向同性的5个参数进行表示，进而可以构建出材料的本构方程。不懂去翻材料力学相关书籍。

这里就随便建立一个高20，半径为1的圆柱体对其进行施加一个单位为1的拉伸位移（较长方向拉伸）。

材料属性就随便找了一篇论文，按照上面的材料属性进行赋值：

为了对比，G23改取3183.1MPa（还是会有一些小误差）。

单元类型全部采用相同的完全积分单元。材料方向指定圆柱体的轴向。

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

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

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

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

Et(轴/横向弹性模量）、Ep（面内弹性模量）、Nutp（面外泊松比，也就是说假设1为轴向，那么Nu12/13的值为Nutp）、Nup（面内泊松比）、Gt（轴/横向剪切模量）

（这里的轴/横向对abaqus来说是局部方向的1方向）

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

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

到此结束。cae文件如下：

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

横观各向同性umat：

共享库：