Télécharger pstres.eso

Retour à la liste

Numérotation des lignes :

pstres
  1. C PSTRES SOURCE CB215821 17/07/21 21:15:29 9513
  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 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
  32. call pstre3 (sigma, sgpri)
  33. end if
  34. c
  35. return
  36. end
  37.  
  38.  
  39.  

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