找回密码
 注册
Simdroid-非首页
查看: 308|回复: 7

[子程序] ABAQUS帮助文档例题中UMAT的应力更新

[复制链接]
发表于 2018-5-18 09:26:56 | 显示全部楼层 |阅读模式 来自 湖北襄阳
在学习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程序解读中的注释,做了下图,请大家帮忙看看,是否存在误解

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
 楼主| 发表于 2018-5-18 10:00:26 | 显示全部楼层 来自 湖北襄阳
Simdroid开发平台
看到程序里面在应力更新之前,有这段应力清零的指令,说明这种应变全量的应力更新方法是正确的。      
IF (NDI .EQ. 3) THEN
         DO I = 1, NTENS
            STRESS(I)=ZERO
现在疑问是用应变全量更新与应变增量更新结果不同,还是以单向加载为例画个图
应变增量来更新应力的结果((1)式)与采用应变全量更新的结果((2)式)有不同,这里哪种是对的呢
应变增量更新应力的程序我也看到过,两种都有使用。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 3 不支持 0

使用道具 举报

发表于 2018-5-18 10:09:32 | 显示全部楼层 来自 辽宁盘锦
支持,支持楼主!
回复 不支持

使用道具 举报

发表于 2018-11-13 17:09:43 | 显示全部楼层 来自 湖北武汉
增量的迭代过程是要按照物理事实拉回的,拉回到屈服面,一维下面对应的是红虚线那个终点吧
回复 不支持

使用道具 举报

发表于 2018-12-12 09:44:48 | 显示全部楼层 来自 陕西汉中
正好进来学习一下
回复 不支持

使用道具 举报

发表于 2019-4-13 16:05:52 | 显示全部楼层 来自 大连理工大学西山生活区
楼主您好,请问您搞明白了吗?哪种应力更新的方式是对的?
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-25 19:34 , Processed in 0.032635 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表