Télécharger dfhardcc.eso

Retour à la liste

Numérotation des lignes :

  1. C DFHARDCC SOURCE CB215821 16/04/21 21:16:27 8920
  2. FUNCTION dfhardCC(rkappac,fc,fc0,rkappa1,ac,bc)
  3.  
  4. c This subroutine calculates the derivative of the compressive
  5. c hardening function with respect to rkappac
  6.  
  7. IMPLICIT REAL*8 (A-B,D-H,O-Z)
  8. implicit integer (I-K,M,N)
  9. implicit logical (L)
  10. implicit character*10 (C)
  11.  
  12. r0 = 0.
  13. hc = bc
  14.  
  15. if (rkappac.le.rkappa1) then
  16. rn = (rkappac/rkappa1)
  17. rloc = 2.*(fc-fc0)*rn/(1.0d0+(rn*rn))
  18. ry = (EXP(ac*rkappac))*(fc0+rloc)
  19. rloc2 = 2.*(fc-fc0)*(1.0d0-(rn*rn))/
  20. . (rkappa1*(1.0d0+(rn*rn))*(1.0d0+(rn*rn)))
  21. dry = ac*ry + (EXP(ac*rkappac))*rloc2
  22. else
  23. rloc = (1.0d0 + hc*(rkappac-rkappa1)) *
  24. . (EXP(-bc*(rkappac-rkappa1)))
  25. ry = (EXP(ac*rkappac))*fc*rloc
  26. rloc2 = fc * (hc-bc-hc*bc*(rkappac-rkappa1))
  27. dry = ac*ry + (EXP((ac-bc)*rkappac + bc*rkappa1)) * rloc2
  28. endif
  29.  
  30. dfhardcc = dry
  31.  
  32. RETURN
  33. END
  34.  
  35.  
  36.  
  37.  
  38.  
  39.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales