不好意思各位大佬,上次写的UMAT案例代码似乎不大好,发现和abaqus的结果对不上,所以这次写了个相对标准的elastic*的code。
先上代码
模型尺寸是一个1:1:1的正方体,材料性质如下:
边界条件和加载条件是一端FIX,另外一端通过couple一个reference point施加0.1单位的位移,所以算是单向受压吧。
mesh 的尺寸是0.5,记得在element type里面取消掉reduced integration, 不然会有hourglass 的问题
跑完看结果,和abaqus内置的elastic模型做对比
搞定!