[Fortran] 纯文本查看 复制代码
C USER INPUT FOR ADAPTIVE MESH CONSTRAINT
SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK,LSMOOTH)
INCLUDE 'ABA_PARAM.INC'
C USER DEFINED DIMENSION STATEMENTS
DIMENSION ULOCAL(*)
DIMENSION ALOCAL(NDIM,*)
DIMENSION JGVBLOCK(*),JMATYP(*)
DIMENSION TIME(2)
DIMENSION ARRAY(15)
DIMENSION JELEMLIST(10),JELEMTYPE(10)
double precision k,PEEQ
C
NELEMS=10
k=1.0
C
CALL GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,JRCD,JGVBLOCK)
C
CALL GETVRMAVGATNODE(NODE,JTYP,'PE',ARRAY,JRCD,JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
C
PEEQ=ARRAY(7)
C
if(PEEQ.GT.0.0) then
ULOCAL(NDIM)=-k*PEEQ
end if
RETURN
END