Esope source files

   1 : C PSTRES    SOURCE    CHAT      05/01/13    02:38:22     5004           
   2 : c---------------------------------------------------------------------
   3 : c
   4 :       SUBROUTINE PSTRES (NSTR1, SIGMA, SGPRI)
   5 : c
   6 : c=====================================================================
   7 : c                                                                    =
   8 : c   This routine calculates the principal stress values of           =
   9 : c   sigma (nstr1) and puts the results in sgpri (3).                 =
  10 : c                                                                    =
  11 : c   Note:    sigma = (Sxx, Syy, Szz, Sxy, Sxz, Syz)   for NSTR1=6    =
  12 : c            sgpri = (S11, S22, S33)                                 =
  13 : c            sigma = (Sxx, Syy, Szz, Sxy)             for NSTR1=4    =
  14 : c            sgpri = (S11, S22, Szz)                                 =
  15 : c                                                                    =
  16 : c=====================================================================
  17 :       IMPLICIT INTEGER(I-N)
  18 :       integer*4 nstr1
  19 :       real*8    sigma (nstr1), sgpri (3)
  20 : c
  21 :       real*8    avalu, bvalu
  22 : c
  23 :       if      (nstr1 .eq. 4)  then
  24 :         avalu = (sigma (1) + sigma (2)) * 0.5 d0
  25 :         bvalu = SQRT (0.25d0*(sigma(1)-sigma(2))**2+sigma(4)**2)
  26 : c
  27 :         sgpri (1) = avalu + bvalu
  28 :         sgpri (2) = avalu - bvalu
  29 :         sgpri (3) = sigma (3)
  30 : c        write(6,*)'sgpri',(sgpri(ie1),ie1=1,3)
  31 :       else if (nstr1 .eq. 6)  then
  31 :         call pstre3 (sigma, sgpri)
  32 :       end if
  33 : c
  34 :       return
  35 :       end   
  36 :  
  37 :