Télécharger fhardcc.eso

Retour à la liste

Numérotation des lignes :

fhardcc
  1. C FHARDCC SOURCE PV 21/12/18 07:15:03 11240
  2. FUNCTION fhardCC(rkappac,fc,fc0,rkappa1,ac,bc,lerror)
  3.  
  4. c hardening in effective stress!
  5. c strain hardening inspired from concrete ETC curve in uniaxial compression
  6. c for the ascending branch, i.e. until peak stress and peak stress strain
  7. c then the descending branch is an exponential
  8.  
  9. c Elastic from (0;0) to (0;fc0)
  10. c The curve begins on (0;fc0),
  11. c then ascending branch based on ETC until peak
  12. c then descending branch exponential
  13. c and at the junction, the tangent is horizontal
  14.  
  15. IMPLICIT REAL*8 (A-B,D-H,O-Z)
  16. implicit integer (I-K,M,N)
  17. implicit logical (L)
  18. implicit character*10 (C)
  19.  
  20. hc = bc
  21.  
  22. if (rkappac.le.rkappa1) then
  23. rn = (rkappac/rkappa1)
  24. rloc = 2.*(fc-fc0)*rn/(1.+(rn*rn))
  25. ry = (EXP(ac*rkappac))*(fc0+rloc)
  26. else
  27. rloc = (1. + hc*(rkappac-rkappa1)) *
  28. . (EXP(-bc*(rkappac-rkappa1)))
  29. ry = (EXP(ac*rkappac))*fc*rloc
  30. endif
  31.  
  32. if (ABS(ry).gt.10.E+25) then
  33. lerror = .true.
  34. return
  35. endif
  36.  
  37. fhardCC= ry
  38.  
  39. RETURN
  40. END
  41.  
  42.  
  43.  
  44.  
  45.  
  46.  

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