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