Télécharger partit.eso

Retour à la liste

Numérotation des lignes :

  1. C PARTIT SOURCE CB215821 17/07/21 21:15:22 9513
  2. c---------------------------------------------------------------------
  3. c
  4. SUBROUTINE PARTIT (NDIME, NSTR1, NTENS, NCOMP,
  5. . betam, sgeff, sgten, sgcom)
  6. c
  7. c=====================================================================
  8. c =
  9. c This subroutine determines principal stresses and directions and =
  10. c performs split of stress tensor into tensile and compressive =
  11. c contributions (eigen system). =
  12. c ------------ =
  13. c Variables: =
  14. c =
  15. c ntens no. of tension componentes =
  16. c ncomp " " compre. " =
  17. c =
  18. c betam rotation matrix =
  19. c sgeff effect. stress tensor to be splitted =
  20. c sgten tension stress tensor contribution =
  21. c sgcom compre. " " " =
  22. c =
  23. c stpri eigen stresses =
  24. c =
  25. c Author: Rui Faria Version: 93.08.25 =
  26. c -------- =
  27. c=====================================================================
  28. IMPLICIT INTEGER(I-N)
  29. integer ndime, nstr1, ntens, ncomp
  30. real*8 sgeff (nstr1), sgten (nstr1), sgcom (nstr1),
  31. . betam (ndime,ndime)
  32. c
  33. integer istrs
  34. real*8 eigen, stpri (3)
  35. c=====----------------------------------------------------------------
  36. c determination of principal stresses and directions |
  37. c=====----------------------------------------------------------------
  38. C+CCd
  39. if (nstr1 .eq. 6) then
  40. call pstres (6, sgeff, stpri)
  41. call pdirec (3, nstr1, sgeff, stpri, betam)
  42. else
  43. call pdirec (2, nstr1, sgeff, stpri, betam)
  44. end if
  45. C CCd
  46. c=====----------------------------------------------------------------
  47. c partition of stress tensor |
  48. c=====----------------------------------------------------------------
  49. ntens = 0
  50. ncomp = 0
  51. c
  52. do 10 istrs = 1, 3
  53. eigen = stpri (istrs)
  54. if (eigen .ge. 0.0 d0) then
  55. ntens = ntens + 1
  56. sgten (istrs) = eigen
  57. sgcom (istrs) = 0.0 d0
  58. else
  59. ncomp = ncomp + 1
  60. sgten (istrs) = 0.0 d0
  61. sgcom (istrs) = eigen
  62. end if
  63. 10 continue
  64. c
  65. return
  66. end
  67.  
  68.  
  69.  
  70.  

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