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

[子程序] VUMAT 网格删除问题

[复制链接]
发表于 2016-7-27 22:11:25 | 显示全部楼层 |阅读模式 来自 陕西西安
我自己做了一个简单的模型,变了一个VUMAT子程序,主要是用来测试网格是否能顺利删除。
CAE模型为悬臂梁一段施加力,如下图所示:
材料参数我是直接在CAE中写入的,如下图所示:
分析步中勾选了SDV和STATUS,如下图所示:
提交作业之后运算,没有删除网格,如此下图所示:
附上我的简单的子程序代码。
C**********call user subroutine VUMAT*******************
C**********write fixed head file here*******************
      subroutine vumat(
C Read only -
     1  nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,
     2  stepTime, totalTime, dt, cmname, coordMp, charLength,
     3  props, density, strainInc, relSpinInc,
     4  tempOld, stretchOld, defgradOld, fieldOld,
     5  stressOld, stateOld, enerInternOld, enerInelasOld,
     6  tempNew, stretchNew, defgradNew, fieldNew,
C Write only -
     7  stressNew, stateNew, enerInternNew, enerInelasNew)
C
      include 'vaba_param.inc'

C All arrays dimensioned by (*) are not used in this algorithm
      dimension props(nprops), density(nblock),
     1  coordMp(nblock,*),
     2  charLength(*), strainInc(nblock,ndir+nshr),
     3  relSpinInc(*), tempOld(*),
     4  stretchOld(*), defgradOld(*),
     5  fieldOld(*), stressOld(nblock,ndir+nshr),
     6  stateOld(nblock,nstatev), enerInternOld(nblock),
     7  enerInelasOld(nblock), tempNew(*),
     8  stretchNew(*), defgradNew(*), fieldNew(*),
     9  stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev),
     1  enerInternNew(nblock), enerInelasNew(*)

C*******************************************************

      character*80 cmname
      PARAMETER (ZERO=0.D0,ONE=1.D0,TWO=2.D0)

C****************head file ends here********************

C*************user code begins here*********************

C*************define state******************************

C      DO i = 1,nblock

C*************CHECK TOTALTIME***************************     

      IF(totalTime.GT.0.0005)THEN
         DO j=100,150,1
            stateNew(j,7)=ZERO
         END DO
      ELSE
            stateNew(i,7)=ONE
      END IF

C*******************************************************

C User code ends here

C*******************************************************  

C      end DO

      return
      end
C*******************************************************

分析步有0.003秒,想设置在totaTtime超过0.001秒之后删除一些网格,结果总是不行。

想咨询各位:
1)我的CAE中是不是有一些设置不对或者漏掉?
2)我的子程序没有材料定义,只有网格状态设置,不知道这样是不是abaqus就不会调用子程序,所以才不删除网格?
3)调用VUMAT必须得在子程序中编辑材料参数吗?像我这样只用子程序来进行网格删除是不是不可行呢?

多谢各位!

本帖子中包含更多资源

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

×
发表于 2016-9-5 15:06:41 | 显示全部楼层 来自 陕西西安
Simdroid开发平台
非常感谢,正在学习
回复 不支持

使用道具 举报

发表于 2018-10-28 23:17:05 来自手机 | 显示全部楼层 来自 北京
问:想设置在totaTtime超过0.001秒之后删除一些网格
回复 不支持

使用道具 举报

发表于 2018-10-28 23:18:25 来自手机 | 显示全部楼层 来自 北京
问:想设置在totaTtime超过0.001秒之后删除一些网格。答:你会把所有网格都删除掉的。
回复 不支持

使用道具 举报

发表于 2019-2-28 21:22:19 | 显示全部楼层 来自 陕西西安
同问,我也出现了网格无法删除的现象,强行删除也不行
回复 不支持

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-3-28 23:28 , Processed in 0.048324 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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