umat.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	umat (stress, statev, ddsdde, sse, spd, scd, rpl, ddsddt, drplde, drpldt, stran, dstran, time, dtime, temp, dtemp, predef, dpred, cmname, ndi, nshr, ntens, nstatv, props, nprops, coords, drot, pnewdt, celent, dfgrd0, dfgrd1, noel, npt, layer, kspt, kstep, kinc)

Function/Subroutine Documentation

umat()

subroutine umat	(	real*8, dimension(ntens)	stress,
		real*8, dimension(nstatv)	statev,
		real*8, dimension(ntens,ntens)	ddsdde,
		real*8	sse,
		real*8	spd,
		real*8	scd,
		real*8	rpl,
		real*8, dimension(ntens)	ddsddt,
		real*8, dimension(ntens)	drplde,
		real*8	drpldt,
		real*8, dimension(ntens)	stran,
		real*8, dimension(ntens)	dstran,
		real*8, dimension(2)	time,
		real*8	dtime,
		real*8	temp,
		real*8	dtemp,
		real*8	predef,
		real*8	dpred,
		character*80	cmname,
		integer	ndi,
		integer	nshr,
		integer	ntens,
		integer	nstatv,
		real*8, dimension(nprops)	props,
		integer	nprops,
		real*8, dimension(3)	coords,
		real*8, dimension(3,3)	drot,
		real*8	pnewdt,
		real*8	celent,
		real*8, dimension(3,3)	dfgrd0,
		real*8, dimension(3,3)	dfgrd1,
		integer	noel,
		integer	npt,
		integer	layer,
		integer	kspt,
		integer	kstep,
		integer	kinc
	)

!
!     here, an ABAQUS umat routine can be inserted
!
!     note that reals should be double precision (REAL*8)
!
      implicit none
!
      character*80 cmname
!
      integer ndi,nshr,ntens,nstatv,nprops,noel,npt,layer,kspt,
     &  kstep,kinc
!
      real*8 stress(ntens),statev(nstatv),
     &  ddsdde(ntens,ntens),ddsddt(ntens),drplde(ntens),
     &  stran(ntens),dstran(ntens),time(2),celent,
     &  props(nprops),coords(3),drot(3,3),dfgrd0(3,3),dfgrd1(3,3),
     &  sse,spd,scd,rpl,drpldt,dtime,temp,dtemp,predef,dpred,
     &  pnewdt
!
!     START EXAMPLE LINEAR ELASTIC MATERIAL
!
      integer i,j
      real*8 e,un,al,um,am1,am2
!
c      write(*,*) 'noel,npt ',noel,npt
c      write(*,*) 'stress ',(stress(i),i=1,6)
c      write(*,*) 'stran ',(stran(i),i=1,6)
c      write(*,*) 'dstran ',(dstran(i),i=1,6)
c      write(*,*) 'drot ',((drot(i,j),i=1,3),j=1,3)
      e=props(1)
      un=props(2)
      al=un*e/(1.d0+un)/(1.d0-2.d0*un)
      um=e/2.d0/(1.d0+un)
      am1=al+2.d0*um
      am2=um
!
!     stress
!      
      stress(1)=stress(1)+am1*dstran(1)+al*(dstran(2)+dstran(3))
      stress(2)=stress(2)+am1*dstran(2)+al*(dstran(1)+dstran(3))
      stress(3)=stress(3)+am1*dstran(3)+al*(dstran(1)+dstran(2))
      stress(4)=stress(4)+am2*dstran(4)
      stress(5)=stress(5)+am2*dstran(5)
      stress(6)=stress(6)+am2*dstran(6)
!
!     stiffness
!
      do i=1,6
         do j=1,6
            ddsdde(i,j)=0.d0
         enddo
      enddo
      ddsdde(1,1)=al+2.d0*um
      ddsdde(1,2)=al
      ddsdde(2,1)=al
      ddsdde(2,2)=al+2.d0*um
      ddsdde(1,3)=al
      ddsdde(3,1)=al
      ddsdde(2,3)=al
      ddsdde(3,2)=al
      ddsdde(3,3)=al+2.d0*um
      ddsdde(4,4)=um
      ddsdde(5,5)=um
      ddsdde(6,6)=um
!
!     END EXAMPLE LINEAR ELASTIC MATERIAL
!
      return