Télécharger cremat.eso

Retour à la liste

Numérotation des lignes :

cremat
  1. C CREMAT SOURCE CHAT 05/01/12 22:29:25 5004
  2. SUBROUTINE CREMAT(D,EX,PX,NSTRS,IFOUR)
  3. C
  4. C FORMATION DE LA MATRICE D ELASTICITE LINEAIRE
  5. C-----------------------------------------------------------------------
  6. C VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU ET NECOU
  7. C
  8. C IFOUR INDICE DU TYPE DE PROBLEME
  9. C -2 CONTRAINTES PLANES
  10. C -1 DEFORMATIONS PLANES
  11. C 0 AXISYMETRIQUE
  12. C 1 SERIE DE FOURIER
  13. C 2 TRIDIMENSIONNEL
  14. C-----------------------------------------------------------------------
  15. IMPLICIT INTEGER(I-N)
  16. IMPLICIT REAL*8(A-H,O-Z)
  17. DIMENSION D5(6,6),D(4,4)
  18. C
  19. CALL ZERO(D5,6,6)
  20. CALL ZERO(D,4,4)
  21. C
  22. GO TO (10,1000,30,40,50,60),NSTRS
  23. GOTO 1000
  24. C
  25. 60 CONTINUE
  26. C
  27. C COMPORTEMENT TRIDIMENSIONNEL
  28. C
  29. E1=EX/(1.D0+PX)/(1.D0-2.D0*PX)
  30. D5(1,1)=E1*(1.D0-PX)
  31. D5(1,2)=E1*PX
  32. D5(1,3)=D5(1,2)
  33. D5(2,1)=D5(1,2)
  34. D5(2,2)=D5(1,1)
  35. D5(2,3)=D5(1,3)
  36. D5(3,1)=D5(1,2)
  37. D5(3,2)=D5(1,2)
  38. D5(3,3)=D5(1,1)
  39. D5(4,4)=E1*(1.D0-2.D0*PX)/2.D0
  40. D5(5,5)=D5(4,4)/1.2D0
  41. D5(6,6)=D5(5,5)
  42. GO TO 100
  43. C
  44. 50 CONTINUE
  45. IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN
  46. C
  47. C CONTRAINTES PLANES AVEC TOUS LES CISAILLEMENTS
  48. C
  49. E1=EX/(1.D0-PX*PX)
  50. D5(1,1)=E1
  51. D5(1,2)=PX*E1
  52. D5(2,1)=PX*E1
  53. D5(2,2)=E1
  54. D5(3,3)=E1*(1.D0-PX)/2.D0
  55. D5(4,4)=D5(3,3)/1.2D0
  56. D5(5,5)=D5(4,4)
  57. C
  58. ELSE IF(IFOUR.EQ.-1) THEN
  59. C
  60. C DEFORMATIONS PLANES AVEC TOUS LES CISAILLEMENTS
  61. C
  62. E1=EX/(1.D0+PX)/(1.D0-2.D0*PX)
  63. D5(1,1)=E1*(1.D0-PX)
  64. D5(1,2)=PX*E1
  65. D5(2,1)=PX*E1
  66. D5(2,2)=E1*(1.D0-PX)
  67. D5(3,3)=E1*(1.D0-2.D0*PX)/2.D0
  68. D5(4,4)=D5(3,3)/1.2D0
  69. D5(5,5)=D5(4,4)
  70. ELSE
  71. GO TO 1000
  72. ENDIF
  73. GO TO 100
  74. C
  75. 40 CONTINUE
  76. IF(IFOUR.EQ.0.OR.IFOUR.EQ.-1) THEN
  77. C
  78. C CAS AXISYMETRIQUE OU DEFORMATION PLANES
  79. C
  80. E1=EX/(1.D0+PX)/(1.D0-2.D0*PX)
  81. D5(1,1)=E1*(1.D0-PX)
  82. D5(1,2)=E1*PX
  83. D5(1,3)=E1*PX
  84. D5(2,1)=E1*PX
  85. D5(2,2)=E1*(1.D0-PX)
  86. D5(2,3)=E1*PX
  87. D5(3,1)=E1*PX
  88. D5(3,2)=E1*PX
  89. D5(3,3)=E1*(1.D0-PX)
  90. D5(4,4)=E1*(1.D0-2.D0*PX)/2.D0
  91. C
  92. ELSE IF (IFOUR.EQ.-2) THEN
  93. C
  94. C CONTRAINTES PLANES
  95. C
  96. E1=EX/(1.D0-PX*PX)
  97. D5(1,1)=E1
  98. D5(1,2)=PX*E1
  99. D5(2,1)=PX*E1
  100. D5(2,2)=E1
  101. D5(3,3)=D5(1,1)*1.D-6
  102. D5(4,4)=E1*(1.D0-PX)/2.D0
  103. C
  104. ENDIF
  105. GO TO 100
  106. C
  107. 30 CONTINUE
  108. IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN
  109. C
  110. C CONTRAINTES PLANES SANS CISAILLEMENTS TRANSVERSAL
  111. C
  112. E1=EX/(1.D0-PX*PX)
  113. D5(1,1)=E1
  114. D5(1,2)=PX*E1
  115. D5(2,1)=PX*E1
  116. D5(2,2)=E1
  117. D5(3,3)=E1*(1.D0-PX)/2.D0
  118. C
  119. ELSE IF(IFOUR.EQ.-1) THEN
  120. C
  121. C DEFORMATIONS PLANES SANS CISAILLEMENTS TRANSVERSAL
  122. C
  123. E1=EX/(1.D0+PX)/(1.D0-2.D0*PX)
  124. D5(1,1)=E1*(1.D0-PX)
  125. D5(1,2)=PX*E1
  126. D5(2,1)=PX*E1
  127. D5(2,2)=E1*(1.D0-PX)
  128. D5(3,3)=E1*(1.D0-2.D0*PX)/2.D0
  129. ELSE
  130. GO TO 1000
  131. ENDIF
  132. GO TO 100
  133. C
  134. 10 CONTINUE
  135. C
  136. C COMPORTEMENT UNIDIMENSIONNEL
  137. C
  138. D5(1,1)=EX
  139. GO TO 100
  140. 1000 WRITE(*,200) IFOUR,NSTRS
  141. STOP
  142. 100 CONTINUE
  143. C
  144. DO 300 I = 1,NSTRS
  145. DO 400 J = 1,NSTRS
  146. D(I,J) = D5(I,J)
  147. 400 CONTINUE
  148. 300 CONTINUE
  149. 200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS CREMAT',/,
  150. * 10X,'IFOUR=',I2,' NSTRS=',I2,//)
  151. RETURN
  152. END
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  

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