- 积分
- 0
- 注册时间
- 2011-1-27
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-2-19 10:00:08
|
显示全部楼层
来自 安徽合肥
我把整个subroutine都按照f90改了(自由格式)
subroutine vumat(&
nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, stepTime, totalTime, dt, cmname, &
coordMp, charLength, props, density, strainInc, relSpinInc, &
tempOld, stretchOld, defgradOld, fieldOld, stressOld, stateOld, enerInternOld, enerInelasOld, &
tempNew, stretchNew, defgradNew, fieldNew, stressNew, stateNew, enerInternNew, enerInelasNew )
implicit none
integer,parameter::precision = 4 !精度
integer i
!只读
character(len=80)::cmname
integer::nblock,ndir,nshr,nstatev,nfieldv, nprops,lanneal
real(kind=precision)::stepTime ,totalTime ,dt
real(kind=precision)::props(nprops) ,density(nblock) ,strainInc(nblock,ndir+nshr) ,relSpinInc(nblock,nshr) ,coordMp(nblock,2) ,charLength(nblock)
real(kind=precision)::tempOld(nblock) ,stretchOld(nblock,ndir+nshr) ,defgradOld(nblock,ndir+nshr+nshr) ,fieldOld(nblock,nfieldv)
real(kind=precision)::tempNew(nblock) ,stretchNew(nblock,ndir+nshr) ,defgradNew(nblock,ndir+nshr+nshr) ,fieldNew(nblock,nfieldv)
real(kind=precision)::stressOld(nblock,ndir+nshr) ,stateOld(nblock,nstatev), enerInternOld(nblock) ,enerInelasOld(nblock)
!只写
real(kind=precision)::stressNew(nblock,ndir+nshr) ,stateNew(nblock,nstatev) ,enerInternNew(nblock) ,enerInelasNew(nblock)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do i = 1,nblock
end do
return
end |
|