main.f90

Go to the documentation of this file.

PROGRAM test_general_umat
use,intrinsic :: iso_fortran_env
include 'ABA_PARAM.INC'
include 'PARAM_UMAT.INC'

!C     ADD COMMON BLOCKS HERE IF NEEDED ()
!C      COMMON /KBLOCK/KBLOCK

parameter(ntens = 6, nstatev = nsdv, nprops = 6, ndi=3, nshr=3)
parameter(noel = 1, npt = 8)
!
CHARACTER*8 CMNAME
dimension stress(ntens),statev(nstatev),ddsdde(ntens,ntens),ddsddt(ntens),      &
drplde(ntens),stran(ntens),dstran(ntens),time(2),predef(1),dpred(1),            &
props(nprops),coords(3),drot(3,3),dfgrd0(3,3),dfgrd1(3,3)
!
i=1.0d0
j=1.0d0
DO i=1,ntens
    DO j=1,ntens
        ddsdde(i,j)=0.0d0
    ENDDO
    stress(i)=0.0d0
ENDDO
!
! DEFORMATION GRADIENT
dfgrd1(1,1)= 1.1d0
dfgrd1(1,2)= 0.0d0
dfgrd1(1,3)= 0.0d0
dfgrd1(2,1)= 0.0d0
dfgrd1(2,2)= 1.0d0/dfgrd1(1,1)
dfgrd1(2,3)= 0.0d0
dfgrd1(3,1)= 0.0d0
dfgrd1(3,2)= 0.0d0
dfgrd1(3,3)= 1.0d0/dfgrd1(1,1)
!
time(1)=0.d0
time(2)=0.d0
call uexternaldb(0,0,time,0.d0,0,0)
!
! MATERIAL PROPERTIES
!
! k PENALTY PARAMETER
props(1)=2.d0/100000.d0
! ISOTROPIC MATRIX
! C10=
props(2)=1.00d0
! C01
props(3)=1.00d0
!k1
props(4)=1.00d0
!k2
props(5)=0.1d0
!kappa
props(6)=0.1d0
!
! !viscous parameters - maxwell
! ! v - number of dashpots
! PROPS(17)=0
! !tau1 %
! PROPS(18)=2.0d0
! !teta1
! PROPS(19)=0.835d0
! !tau2 %
! PROPS(20)=1.2d0
! !teta2
! PROPS(21)=7.0d0
! !tau3 %
! PROPS(22)=12.d0
! !teta3
! PROPS(23)=2.0d0
! !
statev=0.d0
!
erf=0.d0
rho=0.d0
!
!
 dfgrd1(1,1)=  1.0d0
 dfgrd1(1,2)=  0.0d0
 dfgrd1(1,3)=  0.0d0
 dfgrd1(2,1)=  0.0d0
 dfgrd1(2,2)=  1.0d0
 dfgrd1(2,3)=  0.0d0
 dfgrd1(3,1)=  0.0d0
 dfgrd1(3,2)=  0.0d0
 dfgrd1(3,3)=  1.0d0
!
!################################################################################################!
!!     TENSILE MONOTONIC LOAD TEST
 dfgrd1(1,1)=  1.2d0
 dfgrd1(1,2)=  0.0d0
 dfgrd1(1,3)=  0.0d0
 dfgrd1(2,1)=  0.0d0
 dfgrd1(2,2)=  1/sqrt(dfgrd1(1,1))
 dfgrd1(2,3)=  0.0d0
 dfgrd1(3,1)=  0.0d0
 dfgrd1(3,2)=  0.0d0
 dfgrd1(3,3)=  1/sqrt(dfgrd1(1,1))
!
 CALL umat(stress,statev,ddsdde,sse,spd,scd,rpl,ddsddt, drplde,drpldt,stran,     &
dstran,time,dtime,temp,dtemp,predef,dpred,cmname,ndi,nshr,ntens,nstatev,props,  &
nprops,coords,drot,pnewdt,celent,dfgrd0,dfgrd1,noel,npt,layer,kspt,kstep,kinc)

 write(*,*) stress
 write(*,*)
 write(*,*) ddsdde
close(150)
!################################################################################################!
!
END PROGRAM