abaqus的Vgetvrm所提取的第一主应力为什么等于0? 20

浏览:1867 回答:1

具体程序段如下:

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


邀请回答 我来回答

全部回答

(1)
默认 最新
民_2382
大佬我不会您的问题,但是我想问一下您如何确认输出的为0,确认输出的值是什么。难道都要等待漫长的计算之后在从后处理云图里看吗
2023年12月20日
评论 点赞

没解决?试试专家一对一服务

换一批