Télécharger la plaquette Cast3M


fr en

Télécharger equist.eso

Retour à la liste

Numérotation des lignes : 

equist

  1. C EQUIST    SOURCE    CB215821  17/07/21    21:15:08     9513           
  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   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