DF秽土转生 发表于 2018-5-18 09:26:56

ABAQUS帮助文档例题中UMAT的应力更新

在学习Abaqus Example Problems 1.4.6 Failure of blunt notched fiber metal laminates
论坛里也看到有几个UMAT程序详解的帖子在讲这段程序。
该程序中最后部分求解应力，写着
C   UPDATE THE STRESS STATE IF 3D CASE
C
      IF (NDI .EQ. 3) THEN
         DO I = 1, NTENS
            STRESS(I)=ZERO
            DO J = 1, NTENS
               STRESS(I)=STRESS(I)+CDFULL(I,J) * STRANT(J)
            END DO
         END DO
请问大家，STRESS(I)=STRESS(I)+CDFULL(I,J) * STRANT(J)这一句，为什么不是STRESS(I)=STRESS(I)+CDFULL(I,J) * DSTRAN(J)
程序中STRANT为当前迭代步初始应变量，DSTRAN为上次迭代步结束时的应变增量。
我根据UMAT程序解读中的注释，做了下图，请大家帮忙看看，是否存在误解

DF秽土转生 发表于 2018-5-18 10:00:26

看到程序里面在应力更新之前，有这段应力清零的指令，说明这种应变全量的应力更新方法是正确的。      
IF (NDI .EQ. 3) THEN
         DO I = 1, NTENS
            STRESS(I)=ZERO
现在疑问是用应变全量更新与应变增量更新结果不同，还是以单向加载为例画个图
应变增量来更新应力的结果（（1）式）与采用应变全量更新的结果（（2）式）有不同，这里哪种是对的呢
应变增量更新应力的程序我也看到过，两种都有使用。

墨梦呀 发表于 2018-5-18 10:09:32

支持，支持楼主！

blswjs 发表于 2018-11-13 17:09:43

增量的迭代过程是要按照物理事实拉回的，拉回到屈服面，一维下面对应的是红虚线那个终点吧

张大毛666 发表于 2018-12-12 09:44:48

正好进来学习一下

tyut1001 发表于 2019-1-20 17:10:34

支持楼主

qqcyhappy 发表于 2019-4-13 16:05:52

楼主您好，请问您搞明白了吗？哪种应力更新的方式是对的？

莫与陌与沫 发表于 2019-4-18 10:13:16

学习学习
页: [1]
查看完整版本: ABAQUS帮助文档例题中UMAT的应力更新