Esope source files

   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 :