abaqus的Vgetvrm所提取的第一主应力为什么等于0? 20
具体程序段如下:
call vgetvrm( 'SP', rData, jData, cData, jStatus )
SIG1(K)=RDATA(K) 输出值为0
call vgetvrm( 'SP', rData, jData, cData, jStatus )
SIG1(K)=RDATA(2*NBLOCK+K) 有输出值
看完兵哥的视频后就在想是不是2022版本abaqus的主应力也是
S1 S1 S1......NBLOCK
S2 S2 S2......NBLOCK
S3 S3 S3......NBLOCK
或者是
S3 S3 S3......NBLOCK
S2 S2 S2......NBLOCK
S1 S1 S1......NBLOCK
试过了SIG1(K)=RDATA(K),但这样输出是0
如果通过提取主应力的那个子程序或者直接提取应力分量自己算,又会出现云图显示不全(大部分区域为黑色,少部分为彩色)的问题,abaqus2022我看帮助文件已经可以用vgetvrm提取了。。。。
最后附上能算的程序,但损伤值一直输出为0
NCL损伤模型
subroutine vusdfld(
c Read only variables -
1 nblock, nstatev, nfieldv, nprops, ndir, nshr,
2 jElem, kIntPt, kLayer, kSecPt,
3 stepTime, totalTime, dt, cmname,
4 coordMp, direct, T, charLength, props,
5 stateOld,
c Write only variables -
6 stateNew, field )
c
include 'vaba_param.inc'
c
dimension jElem(nblock), coordMp(nblock,*),
1 direct(nblock,3,3), T(nblock,3,3),
2 charLength(nblock), props(nprops),
3 stateOld(nblock,nstatev),
4 stateNew(nblock,nstatev),
5 field(nblock,nfieldv)
character*80 cmname
c
c Local arrays from vgetvrm are dimensioned to
c maximum block size (maxblk)
c
parameter( nrData=6 )
character*3 cData(maxblk*nrData)
dimension rData(maxblk*nrData), jData(maxblk*nrData)
real*8 s(NBLOCK,NDIR+NSHR),SIG1(NBLOCK)
real*8 eigVal(NBLOCK,3)
real*8 Xdamg(NBLOCK),XDD(NBLOCK),SEQ(NBLOCK),SEQE(NBLOCK)
REAL*8 DSEQE(NBLOCK)
c
do k = 1, nblock
C 提取米赛斯等效应力(亲测abaqus2022版本可以)
call vgetvrm( 'MISES', rData, jData, cData, jStatus )
SEQ(K)=rData(K)
C 提取主应力(最搞不清楚这个)
call vgetvrm( 'SP', rData, jData, cData, jStatus )
SIG1(K)=RDATA(K)
STATENEW(K,4)=SIG1(2*NBLOCK+K)
C 计算等效应变(vgetvrm无法提取总应变E,只能提取塑性应变或者对数应变(真实应变)帮助文档里没
C 写,但试了下有值,但为了避免出错还是采用LE)
call vgetvrm( 'LE', rData, jData, cData, jStatus )
EE1=rDATA(K)
EE2=rDATA(K+NBLOCK)
EE3=rDATA(K+2*NBLOCK)
EE4=rDATA(3*NBLOCK+K)
EE5=rDATA(4*NBLOCK+K)
EE6=rDATA(5*NBLOCK+K)
C 计算当前步等效应变
SEQE(K)=sqrt(((EE1-EE2)**2+(EE2-EE3)**2+ (EE1-EE3)**2+
1 6*(EE4**2+EE5**2+EE6**2))*(2.0d0/9.0D0))
C 计算等效应变增量,用分析步结束后的值-前的值
DSEQE(K)=SEQE(K)-STATEOLD(K,3)
C 存储当前步等效应变
STATENEW(K,3)=SEQE(K)
C 导入上一步损伤值
XDAMG(K)=STATEOLD(K,2)
C 若最大主应力小于等于0,则不增加损伤值。
if(SIG1(k).LE.0)then
xdamg(k)=STATEOLD(K,2)
else
C 当前步损伤值=上一步损伤值+最大主应力*等效应变增量/等效应力
XDAMG(K)=XDAMG(K)+SIG1(K)*DSEQE(K)/SEQ(K)
endif
C 更新损伤状态变量
STATENEW(K,2)=XDAMG(K)
C 若损伤状态变量大于等于0.58则删除单元
IF( STATENEW(K,2).GE.0.58)THEN
STATENEW(K,1)=0.0
ELSE
STATENEW(K,1)=1.0
ENDIF
ENDDO
c
return
end