请问怎么实现双曲线

yejet

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收集的资料