一个子程序问题
大家好,想跟大家请教一下下面这个程序,我不论怎么在marc中运行,都不能够实现,正常的数据,好像计算不对.但是我又看不出原因,问题就在这个程序里面,期待有高手可以帮忙解决一下.谢谢大家了.
subroutine ueldam(m,n,nn,kc,inc,lovl,matus,timinc,cptim,toten,
2 deven,totend,totenv,surfc,surfd,dt,dtdl,
3 damd,ddamd)
c
include '../common/implicit'
dimensiondt(*),dtdl(*),matus(2)
c
c
c User subroutine to determine value of damage parameter
c
c Input:
c
c m = element number
c n = element/elsto number
c nn = integration point number
c kc = layer number
c inc = increment number
c lovl = 4 for assembly phase, 6 for stress recovery phase
c matus(1) = user material set number
c matus(2) = internal material set number
c timinc = time increment
c cptim = time at beginning of increment
c toten = total instantaneous strain energy at the end
c of the current step excluding damage
c deven = deviatoric part of the instantaneous strain energy at
c the end of the current step excluding damage
c totend = stored deviatoric energy at previous step (including damage)
c totenv = stored volumetric energy at previous step (including damage)
c surfc = current radius of continuous damage surface
c surfd = current radius of discontinuous damage surface
c dt = temperature
c dtdl = incremental temperature
c
c Required Output:
c
c damd = value of Kachanov deviatoric damage parameter
c ddamd = derivative of damage parameter with respect to maximum
c total strain energy
c
c
c include 'concom'
common/rubber/WW,vvv,fovinc(1600,9),fov(1600,9)
common/rubber/endom(1600,9)
vvv=toten
dinfalf=0.7719389
nualf=1.2
kc=1.D0
write(6,*) 'dans la routine'
call elmvar (1,m,nn,kc,strain)
fov(m,nn)=strain+1
if(abs(fovinc(m,nn)).gt.abs(fov(m,nn))) then
fov(m,nn)=fovinc(m,nn)
else
fov(m,nn)=fov(m,nn)
endif
damd=1-dinfalf*(1-exp((1-fov(m,nn))*nualf))
endom(m,nn)=damd
return
end 各位,有时间看看问题啊...很紧急,很重要....
页:
[1]