-
Notifications
You must be signed in to change notification settings - Fork 5
/
nhinterface.f90
26 lines (26 loc) · 1.21 KB
/
nhinterface.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
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)
use numerichyper, only: update_umat
use umatutils, only: dp
implicit none
! This is a hack. The content of the 'aba_param.inc' is simply the
! Following two lines. I commented the first one, and modified the
! second one.
! include 'aba_param.inc'
! implicit real*8(a-h,o-z)
! parameter (nprecd=2)
integer, parameter :: nprecd=2
character*80, intent(in) :: cmname
real(dp), intent(in) :: stran(ntens),dstran(ntens),time(2),predef(1),dpred(1),&
props(nprops),coords(3),drot(3,3),dfgrd0(3,3),dfgrd1(3,3), dtime, temp, &
dtemp, celent
integer, intent(in) :: ndi, nshr, ntens, nstatv, nprops, noel, npt, layer,&
kspt, kstep, kinc
real(dp), intent(inout) :: stress(ntens), statev(nstatv), sse, spd, scd,&
rpl, ddsdde(ntens, ntens), ddsddt(ntens), drplde(ntens), drpldt, pnewdt
! Update the stress and ddsdde
call update_umat(props, dfgrd1, ntens, stress, ddsdde, statev)
end subroutine umat