Numérotation des lignes :

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