Télécharger eigglo.eso

Retour à la liste

Numérotation des lignes :

eigglo
  1. C EIGGLO SOURCE CB215821 17/07/21 21:15:08 9513
  2. c---------------------------------------------------------------------
  3. c
  4. SUBROUTINE EIGGLO (NDIME, NSTR1, NTENS,
  5. . betam, prosi, sgten, soxyz)
  6. c
  7. c=====================================================================
  8. c =
  9. c This subroutine performs rotation of tension stress matrix, =
  10. c from eigen to global system (optimized). =
  11. c =
  12. c Input variable: =
  13. c =
  14. c ntens no. of tension componentes =
  15. c betam rotation matrix =
  16. c sgten (i) (i=1,3) principal tension stress tensor =
  17. c =
  18. c Output variable: =
  19. c =
  20. c sgten tension stress tensor contribution =
  21. c =
  22. c Author: Rui Faria Version: 96.09.18 =
  23. c -------- =
  24. c=====================================================================
  25. IMPLICIT INTEGER(I-N)
  26. integer ndime, nstr1, ntens
  27. real*8 sgten (nstr1),
  28. . betam (ndime,ndime), prosi (ndime,ndime),
  29. . soxyz (ndime,ndime)
  30. c
  31. integer istrs, jstrs
  32. real*8 provi
  33. c
  34. if (ntens .gt. 0) then
  35. do 10 istrs = 1, ndime
  36. do 20 jstrs = 1, ndime
  37. prosi (istrs,jstrs) = betam (jstrs,istrs) *
  38. . sgten (jstrs)
  39. 20 continue
  40. 10 continue
  41. c
  42. call proma1 (soxyz, prosi, betam, ndime, ndime, ndime)
  43. c
  44. if (nstr1 .eq. 6) then
  45. sgten (1) = soxyz (1,1)
  46. sgten (2) = soxyz (2,2)
  47. sgten (3) = soxyz (3,3)
  48. sgten (4) = soxyz (1,2)
  49. sgten (5) = soxyz (1,3)
  50. sgten (6) = soxyz (2,3)
  51. else
  52. sgten (1) = soxyz (1,1)
  53. sgten (2) = soxyz (2,2)
  54. provi = sgten (3)
  55. sgten (3) = provi
  56. sgten (4) = soxyz (1,2)
  57. end if
  58. else
  59. call zeror (sgten, nstr1)
  60. end if
  61. c
  62. return
  63. end
  64.  
  65.  
  66.  
  67.  

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