Télécharger pstre3.eso

Retour à la liste

Numérotation des lignes :

pstre3
  1. C PSTRE3 SOURCE PV 22/04/25 21:15:10 11344
  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=====================================================================
  15. IMPLICIT REAL*8 (A-H,O-Z)
  16. IMPLICIT INTEGER(I-N)
  17. real*8 sigma (6), sgpri (3)
  18. c
  19. c2=-SIGMA(1)-SIGMA(2)-SIGMA(3)
  20. c1= (SIGMA(1)*SIGMA(2)+SIGMA(2)*SIGMA(3)+SIGMA(3)*SIGMA(1))
  21. & - SIGMA(5)**2 - SIGMA(4)**2 - SIGMA(6)**2
  22. c0=-2.*SIGMA(4)*SIGMA(5)*SIGMA(6) + SIGMA(1)*SIGMA(6)**2
  23. & + SIGMA(2)*SIGMA(5)**2 + SIGMA(3)*SIGMA(4)**2
  24. & - SIGMA(1)*SIGMA(2)*SIGMA(3)
  25. call degre3(c0,c1,c2,d1,XI1,d2,XI2,d3,XI3)
  26.  
  27. sgpri (1) = D1
  28. sgpri (2) = D2
  29. sgpri (3) = D3
  30.  
  31. return
  32. end
  33.  
  34.  
  35.  
  36.  
  37.  

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