Télécharger parinv.eso

Retour à la liste

Numérotation des lignes :

parinv
  1. C PARINV SOURCE CHAT 05/01/13 02:09:11 5004
  2. C PARINV SOURCE INSL 24/10/96
  3. SUBROUTINE PARINV(S,SX,SY,SZ,SXY,SXZ,SYZ,NSTRS)
  4. C ****************************************************
  5. C * CALCUL DES IVARIANT : I1, J2, J3 et COS(3TETA) *
  6. C * ( MODEL D'OTTOSEN ) *
  7. C ****************************************************
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. DIMENSION S1(6),S(NSTRS)
  11. C
  12. C INITIALISATION
  13. C
  14. CALL ZERO(S1,6,1)
  15. DO 11 I=1,NSTRS
  16. S1(I)=S(I)
  17. 11 CONTINUE
  18. C
  19. SX=0.D0
  20. SY=0.D0
  21. SZ=0.D0
  22. SXY=0.D0
  23. SXZ=0.D0
  24. SYZ=0.D0
  25. C
  26. GO TO (10,1000,30,40,50,60),NSTRS
  27. GOTO 1000
  28. C
  29. 60 CONTINUE
  30. C
  31. C COMPORTEMENT TRIDIMENSIONNEL
  32. C
  33. SX =S1(1)
  34. SY =S1(2)
  35. SZ =S1(3)
  36. SXY=S1(4)
  37. SXZ=S1(5)
  38. SYZ=S1(6)
  39. GO TO 100
  40. C
  41. 50 CONTINUE
  42. C
  43. C CONTRAINTES OU DEFORMATIONS PLANES
  44. C AVEC TOUS LES CISAILLEMENTS
  45. C
  46. SX =S1(1)
  47. SY =S1(2)
  48. SXY=S1(3)
  49. SXZ=S1(4)
  50. SYZ=S1(5)
  51. GO TO 100
  52. C
  53. 40 CONTINUE
  54. C
  55. C CAS AXISYMETRIQUE
  56. C
  57. SX =S1(1)
  58. SY =S1(2)
  59. SZ =S1(3)
  60. SXY=S1(4)
  61. GO TO 100
  62. C
  63. 30 CONTINUE
  64. C
  65. C CONTRAINTES OU DEFORMATIONS PLANES
  66. C SANS CISAILLEMENT TRANSVERSAL
  67. C
  68. SX =S1(1)
  69. SY =S1(2)
  70. SXY=S1(3)
  71. GO TO 100
  72. C
  73. 10 CONTINUE
  74. C
  75. C COMPORTEMENT UNIDIMENSIONNEL
  76. C
  77. SX =S1(1)
  78. GO TO 100
  79. C
  80. 1000 WRITE(*,200) NSTRS
  81. STOP
  82. 200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS PARINV',/,
  83. * 20X,' NSTRS=',I2,//)
  84. 100 CONTINUE
  85. C
  86. RETURN
  87. END
  88.  
  89.  
  90.  
  91.  

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