Abaqus 子程序开发日志06-vumat之isotropic

vumat子程序只用于显示分析步，用来定义用户自定义材料，这里首先对较为简单的各向同性材料本构模型进行构建。

vumat的工作原理与umat类似，只不过vumat不需要再提供雅可比矩阵，以矢量的形式进行储存运算。剪切应变不再是工程应变，而变成了张量的形式。

对于各向同性的刚度矩阵方程在umat部分已经提过，这里直接放方程式：

找到变量之间的对应关系，再写入子程序中就行了。

本构方程是反应材料的固有的力学性能，不会受你用显示还是隐式求解发生改变（好像有点错了），这里还是以E（杨氏模量）与NU（泊松比）对各向同性材料进行表征。

完成子程序的书写后先对其进行检验，abaqus中可以直接定义各向同性材料的材料属性，这里以其对照。

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

材料属性就随便设定为E=2000，v=0.2，密度（在显示分析中需要定义密度）为2e-09。

单元类型全部采用相同的完全积分单元。

分别用自带的各向同性材料属性与vumat子程序进行计算。

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

此次直接使用for文件，在材料属性的输入过程，在通用中找到用户自定义材料，输入设定的参数，第一个为杨氏模量，第二个为泊松比，密度另设。如图所示：

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

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

到此结束。cae文件如下：

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

各向同性vumat：

这是两个对应的共享库

尾部加D的为双精度的，没的是单精度