vumat子程序只用于显示分析步,用来定义用户自定义材料,这里首先对较为简单的各向同性材料本构模型进行构建。
vumat的工作原理与umat类似,只不过vumat不需要再提供雅可比矩阵,以矢量的形式进行储存运算。剪切应变不再是工程应变,而变成了张量的形式。
对于各向同性的刚度矩阵方程在umat部分已经提过,这里直接放方程式:
![Abaqus 子程序开发日志01-umat之isotropic的图1](static/picture/c51d726681224049973a801e0688ce46-formatwebpresizew_760.png)
找到变量之间的对应关系,再写入子程序中就行了。
本构方程是反应材料的固有的力学性能,不会受你用显示还是隐式求解发生改变,这里还是以E(杨氏模量)与NU(泊松比)对各向同性材料进行表征。
完成子程序的书写后先对其进行检验,abaqus中可以直接定义各向同性材料的材料属性,这里以其对照。
这里就随便建立一个长宽高为5;5;20的立方体对其进行施加一个单位为2的拉伸位移(较长方向拉伸)。
材料属性就随便设定为E=2000,v=0.2,密度(在显示分析中需要定义密度)为2e-09。
单元类型全部采用相同的完全积分单元。
分别用自带的各向同性材料属性与vumat子程序进行计算。
自带的材料属性的设定以及模型的建立、边界条件的设定不再说明,只对vumat的使用进行说明。
此次直接使用for文件,在材料属性的输入过程,在通用中找到用户自定义材料,输入设定的参数,第一个为杨氏模量,第二个为泊松比,密度另设。如图所示:
![Abaqus 子程序开发日志06-vumat之isotropic的图2](static/picture/0b7a773c2f7c4110a295605131104f8b-formatwebpresizew_760.png)
点击ok就好,其余的和平时操作一致,只需在最后提交job时提供用户子程序。
![Abaqus 子程序开发日志06-vumat之isotropic的图3](static/picture/afe469e0d2b84e748b32bc9506623844-formatwebp.png)
然后提交计算就好。结果如图所示:
![Abaqus 子程序开发日志06-vumat之isotropic的图4](static/picture/13d1c9c8551f44569b676e652236ae9a-formatwebpresizew_760.png)
到此结束。cae文件如下:
若选用缩减积分单元时要开启沙漏控制。
各向同性vumat:
这是两个对应的共享库
尾部加D的为双精度的,没的是单精度