Télécharger la plaquette Cast3M


fr en

Télécharger pstres.eso

Retour à la liste

Numérotation des lignes : 

pstres

  1. C PSTRES    SOURCE    GOUNAND   26/08/19    21:15:03     12620          
  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 Gounand 2026/08 On renvoie sur jacob3.eso, plus robuste et precise =
  17. c=====================================================================
  18.       IMPLICIT INTEGER(I-N)
  19.       IMPLICIT REAL*8 (A-H,O-Z)
  20.       integer   nstr1
  21.       real*8    sigma (nstr1), sgpri (3),A(3,3)
  22.       real*8    avalu, bvalu
  23. c
  24.       if      (nstr1 .eq. 4)  then
  25.          a(1,1)=sigma(1)
  26.          a(2,2)=sigma(2)
  27.          a(1,2)=sigma(4)
  28.          a(2,1)=sigma(4)
  29.          a(3,3)=sigma (3)
  30.          call jacod2(a,sgpri)
  31.          return
  32. * Old source
  33.         avalu = (sigma (1) + sigma (2)) * 0.5 d0
  34.         bvalu = SQRT (0.25d0*(sigma(1)-sigma(2))**2+sigma(4)**2)
  35. c
  36.         sgpri (1) = avalu + bvalu
  37.         sgpri (2) = avalu - bvalu
  38.         sgpri (3) = sigma (3)
  39. c        write(6,*)'sgpri',(sgpri(ie1),ie1=1,3)
  40.       else if (nstr1 .eq. 6)  then
  41.         call pstre3 (sigma, sgpri)
  42.       end if
  43. c
  44.       return
  45.       end
  46.  
  47.  

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