common

举个栗子

[Fortran] 纯文本查看 复制代码

REAL*8 AINTW(NNODE),XII(NNODE,2),XI(2),dNdxi(NNODE,2),
     1 VJACOB(2,2),dNdx(NNODE,2),VJABOBINV(2,2),AN(4),BP(2,NDOFEL),
     2 DP(2),SDV(NSTV),BB(3,NDOFEL),CMAT(3,3),EPS(3),STRESS(3),
     3 VNI(2,NDOFEL),ULOC(2),PHASENOD(NNODE)
     
       REAL*8 DTM,THCK,HIST,CLPAR,GCPAR,EMOD,ENU,PARK,ENG

       COMMON/KUSER/USRVAR(N_ELEM,NSTV,4)
C
C     ==================================================================
C     ******************************************************************
C     Constructing elemet TYPE 1 (phase)
C     ******************************************************************
C     ==================================================================
C      
       IF (JTYPE.EQ.ONE) THEN
C     ==================================================================
C     Time an iteration variables
C     ==================================================================
       TIMEZ=USRVAR(JELEM,17,1)
       IF (TIMEZ.LT.TIME(2)) THEN
        USRVAR(JELEM,17,1)=TIME(2)
        USRVAR(JELEM,18,1)=ZERO
       ELSE
        USRVAR(JELEM,18,1)=USRVAR(JELEM,18,1)+ONE
       ENDIF
       STEPITER=USRVAR(JELEM,18,1)

common是怎么使用？

vvt

看彭国伦的书第 8 章第 3 节 169页

举个栗子

vvt 发表于 2022-3-5 12:42
看彭国伦的书第 8 章第 3 节 169页

关键是后面的函数USRVSR(N_ELEM,NSTV,4)是怎么运用的；
而且这个函数全局中也没有声明；

vvt

USRVAR 是一个数组，USRVSR我没在你给出的代码里看到

举个栗子

vvt 发表于 2022-3-7 08:28
USRVAR 是一个数组，USRVSR我没在你给出的代码里看到

大佬，那这个数组是怎么使用的，位置就在common后面，这个数组里面有三个参数，到底是怎么使用的？

举个栗子

vvt 发表于 2022-3-7 08:28
USRVAR 是一个数组，USRVSR我没在你给出的代码里看到

REAL*8 DTM,THCK,HIST,CLPAR,GCPAR,EMOD,ENU,PARK,ENG

        COMMON / KUSER / USRVAR (N_ELEM,NSTV,4)


其中USRVAR (N_ELEM,NSTV,4)我感觉这是个函数，但是在代码中没有找到这个函数的用法，就只是在上面common这个位置说明一下，下文就直接调用了；这个for文件作为abaqus中的子程序来进行使用的。

vvt

是个数组，关于数组的用法，请参考彭国伦的书第 7 章

举个栗子

vvt 发表于 2022-3-9 08:14
是个数组，关于数组的用法，请参考彭国伦的书第 7 章

       COMMON/KUSER/USRVAR(N_ELEM,NSTV,4)
C
C     ==================================================================
C     ******************************************************************
C     Constructing elemet TYPE 1 (phase)
C     ******************************************************************
C     ==================================================================
C      
       IF (JTYPE.EQ.ONE) THEN
C     ==================================================================
C     Time an iteration variables
C     ==================================================================
      
       TIMEZ=USRVAR(JELEM,17,1)


       IF (TIMEZ.LT.TIME(2)) THEN
        USRVAR(JELEM,17,1)=TIME(2)
        USRVAR(JELEM,18,1)=ZERO
       ELSE
        USRVAR(JELEM,18,1)=USRVAR(JELEM,18,1)+ONE
       ENDIF
       STEPITER=USRVAR(JELEM,18,1)
若是数组的话，为啥上下还会变？

vvt

COMMON/KUSER/USRVAR(N_ELEM,NSTV,4)
这个意思是，USRVAR是个三维数组，每个维度的大小分别是 N_ELEM , NSTV , 4。即总大小 N_ELEM*NSTV*4。这里 N_ELEM  NSTV  4 是每个维度的大小。

TIMEZ=USRVAR(JELEM,17,1)
这个意思是，让 TIMEZ 的值，变成 USRVAR 这个三维数组的第 JELEM ，17，1 个元素的值。这里 JELEM,17,1 是索引下标。