Télécharger equist.eso

Retour à la liste

Numérotation des lignes :

  1. C EQUIST SOURCE CHAT 05/01/12 23:44:40 5004
  2. c---------------------------------------------------------------------
  3. c
  4. SUBROUTINE EQUIST (STRES, TANGB, YOUNG, IFLAG, EQUIV)
  5. c
  6. c=====================================================================
  7. c =
  8. c This subroutine computes the equivalent strain in the eigen =
  9. c system. ----- =
  10. c ------ =
  11. c =
  12. c Variables: =
  13. c =
  14. c stres (3) desired part of stress tensor =
  15. c tangb tang(beta) (Drucker-Prager) =
  16. c young E =
  17. c equiv equivalent strain =
  18. c =
  19. c iflag: -1 - compression =
  20. c +1 - tension =
  21. c =
  22. c Author: Rui Faria Version: 93.04.20 =
  23. c -------- =
  24. c=====================================================================
  25. IMPLICIT INTEGER(I-N)
  26. integer*4 iflag
  27. real*8 equiv, tangb, young, stres (3)
  28. c
  29. real*8 oneth, root3, octno, octsh, argum
  30. c
  31. parameter (oneth = 0.33333333333333 d0,
  32. . root3 = 1.73205080756888 d0)
  33. c
  34. if (iflag .eq. -1) then
  35. octno = (stres (1) + stres (2) + stres (3)) * oneth
  36. c
  37. octsh = SQRT ((stres (1) - stres (2))**2 +
  38. & (stres (2) - stres (3))**2 +
  39. & (stres (3) - stres (1))**2) * oneth
  40. argum = root3 * (octno * tangb + octsh)
  41. if (argum .gt. 0.0 d0) then
  42. equiv = SQRT (argum)
  43. else
  44. equiv = 0.0 d0
  45. end if
  46. else
  47. c
  48. equiv = SQRT ((stres(1)**2+stres(2)**2+
  49. . stres(3)**2) / young)
  50. endif
  51. c
  52. return
  53. end
  54.  
  55.  
  56.  

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