- 积分
- 0
- 注册时间
- 2015-3-25
- 仿真币
-
- 最后登录
- 1970-1-1
|
我自己做了一个简单的模型,变了一个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必须得在子程序中编辑材料参数吗?像我这样只用子程序来进行网格删除是不是不可行呢?
多谢各位!
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|