Télécharger mathoo.eso

Retour à la liste

Numérotation des lignes :

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

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