ABAQUS子程序UMAT里弹塑本构的实现

409.89949494 411.06797181 412.12435565 413.09580085 414

414.8492424 415.65247584 416.41645516 417.1464282 417.8465683 418.52025918

0.0002 0.00025 0.0003 0.00035 0.0004 0.00045 0.0005 0.00055 0.0006 0.00065 0.0007

421.57544901 422.13594362 422.68259244 423.21637353 423.73815494 424.24871131 424.74873734 425.23885893 425.7196423 426.19160171 426.65520587

0.00095 0.001 0.00105 0.0011 0.00115 0.0012 0.00125 0.0013 0.00135 0.0014 0.00145

采用ABAQUS自带的弹塑性材料应力计算结果如图6-2所示

图6-2 自带弹塑性材料Mises等效应力

图6-2表明Mises等效应力值介于308.3MPa—421.3MPa之间。

应力应变关系如图6-3所示:

43

图6-3 自带弹塑性材料应力应变关系

图6-3表明，材料的屈服点为400Mpa，塑性阶段按幂函数增长，这个结果与参考材料的应力应变关系完全一致。

6.4. 常刚度法的UMAT验证

模型与前文的6.2一致，材料定义如下： 状态变量数：1

用户自定义材料见表6-2：

表6-2 用户自定义材料表

弹性模量E 泊松比Mu

200000 0.3

屈服应力Yield0 400 A B C

采用常刚度法应力计算结果如图6-4所示

700 0.5 400

44