马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
subroutine fric(lm,tau,ddtddg,ddtddp,dslip,sed,spd,
1 ddtddt,pnewdt,statev,dgam,taulm,press,dpress,ddpddh,
2 slip,kstep,kinc,time,dtime,noel,ciname,slname,
3 msname,npt,node,npatch,coords,rcoord,drot,temp,
4 predef,nfdir,mcrd,npred,nstatv,chrlngth,props,nprops)
c
include 'aba_param.inc'
c
character*80 ciname,slname,msname
dimension tau(nfdir),ddtddg(nfdir,nfdir),ddtddp(nfdir),
1 dslip(nfdir),ddtddt(nfdir,2),statev(*),
2 dgam(nfdir),taulm(nfdir),slip(nfdir),time(2),
3 coords(mcrd),rcoord(mcrd),drot(2,2),temp(2),
4 predef(2,*),props(nprops)
c
parameter(zero=0.0D0,two=2.0D0)
c
if (lm .eq. 2) return
c
a=props(1)
b=props(2)
c
lm=0
tau(1)=dslip(1)/(a+b*dslip(1))
tau(2)=dslip(2)/(a+b*dslip(2))
ddtddg(1,1)=tau(1)/slip(1)
ddtddg(2,2)=tau(2)/slip(2)
ddtddg(1,2)=zero
ddtddg(2,1)=zero
dslip(1)=dgam(1)
dslip(2)=dgam(2)
return
end
来自:likuiming收集的资料 |