Télécharger la plaquette Cast3M


fr en

Télécharger pstre3.eso

Retour à la liste

Numérotation des lignes : 

pstre3

  1. C PSTRE3    SOURCE    GOUNAND   26/08/19    21:15:02     12620          
  2. c---------------------------------------------------------------------
  3. c
  4.       SUBROUTINE PSTRE3 (SIGMA, SGPRI)
  5. c
  6. c=====================================================================
  7. c                                                                    =
  8. c   This routine calculates the eigen stress values of 3D            =
  9. c   'sigma (6)' and puts the results in 'sgpri (3)'.                 =
  10. c                                                                    =
  11. c    Note:   sigma = (Sxx, Syy, Szz, Sxy, Sxz, Syz)                  =
  12. c            sgpri = (S11, S22, S33)                                 =
  13. c                                                                    =
  14. C Gounand 2026/08 On renvoie sur jacob3.eso, plus robuste et precise =
  15. c=====================================================================
  16.       IMPLICIT REAL*8 (A-H,O-Z)
  17.       IMPLICIT INTEGER(I-N)
  18.       real*8    sigma (6), sgpri (3),a(3,3)
  19.  
  20.       a(1,1)=sigma(1)
  21.       a(2,2)=sigma(2)
  22.       a(3,3)=sigma(3)
  23.       a(1,2)=sigma(4)
  24.       a(2,1)=sigma(4)
  25.       a(1,3)=sigma(5)
  26.       a(3,1)=sigma(5)
  27.       a(2,3)=sigma(6)
  28.       a(3,2)=sigma(6)
  29.       call jacod3(a,3,sgpri)
  30.       return
  31.  
  32. * Old source
  33. c
  34.       c2=-SIGMA(1)-SIGMA(2)-SIGMA(3)
  35.       c1= (SIGMA(1)*SIGMA(2)+SIGMA(2)*SIGMA(3)+SIGMA(3)*SIGMA(1))
  36.      &   - SIGMA(5)**2 - SIGMA(4)**2 - SIGMA(6)**2
  37.       c0=-2.*SIGMA(4)*SIGMA(5)*SIGMA(6) + SIGMA(1)*SIGMA(6)**2
  38.      &   + SIGMA(2)*SIGMA(5)**2 + SIGMA(3)*SIGMA(4)**2
  39.      &   - SIGMA(1)*SIGMA(2)*SIGMA(3)
  40.       call degre3(c0,c1,c2,d1,XI1,d2,XI2,d3,XI3)
  41.  
  42.       sgpri (1) = D1
  43.       sgpri (2) = D2
  44.       sgpri (3) = D3
  45.  
  46.       return
  47.       end
  48.  
  49.  

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