Télécharger matsou.eso

Retour à la liste

Numérotation des lignes :

  1. C MATSOU SOURCE CHAT 05/01/13 01:35:34 5004
  2. C MATSOU SOURCE INSL 24/10/96
  3. SUBROUTINE MATSOU(D,EX,PX,NSTRS,IFOUR)
  4. C
  5. C FORMATION DE LA MATRICE DE SOUPLESSE
  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. D5(1,1)=1.D0/EX
  31. D5(1,2)=-PX/EX
  32. D5(1,3)=-PX/EX
  33. D5(2,1)=-PX/EX
  34. D5(2,2)=1.D0/EX
  35. D5(2,3)=-PX/EX
  36. D5(3,1)=-PX/EX
  37. D5(3,2)=-PX/EX
  38. D5(3,3)=1.D0/EX
  39. D5(4,4)=2.D0*(1.D0+PX)/EX
  40. D5(5,5)=1.2D0*D5(4,4)
  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. D5(1,1)=1.D0/EX
  50. D5(1,2)=-PX/EX
  51. D5(2,1)=-PX/EX
  52. D5(2,2)=1.D0/EX
  53. D5(3,3)=2.D0*(1.D0+PX)/EX
  54. D5(4,4)=D5(3,3)*1.2D0
  55. D5(5,5)=D5(4,4)
  56. C
  57. ELSE IF(IFOUR.EQ.-1) THEN
  58. C
  59. C DEFORMATIONS PLANES AVEC TOUS LES CISAILLEMENTS
  60. C
  61. E1=(1.D0+PX)/EX
  62. D5(1,1)=E1*(1.D0-PX)
  63. D5(1,2)=-PX*E1
  64. D5(2,1)=-PX*E1
  65. D5(2,2)=E1*(1.D0-PX)
  66. D5(3,3)=E1*2.D0
  67. D5(4,4)=D5(3,3)*1.2D0
  68. D5(5,5)=D5(4,4)
  69. ELSE
  70. GO TO 1000
  71. ENDIF
  72. GO TO 100
  73. C
  74. 40 CONTINUE
  75. IF(IFOUR.EQ.0.OR.IFOUR.EQ.-1) THEN
  76. C
  77. C CAS AXISYMETRIQUE OU DEFORMATION PLANES
  78. C
  79. D5(1,1)=1.D0/EX
  80. D5(1,2)=-PX/EX
  81. D5(1,3)=-PX/EX
  82. D5(2,1)=-PX/EX
  83. D5(2,2)=1.D0/EX
  84. D5(2,3)=-PX/EX
  85. D5(3,1)=-PX/EX
  86. D5(3,2)=-PX/EX
  87. D5(3,3)=1.D0/EX
  88. D5(4,4)=2.D0*(1.D0+PX)/EX
  89. C
  90. ELSE IF(IFOUR.EQ.-2) THEN
  91. C
  92. C CONTRAINTES PLANES
  93. C
  94. D5(1,1)=1.D0/EX
  95. D5(1,2)=-PX/EX
  96. D5(2,1)=-PX/EX
  97. D5(2,2)=1.D0/EX
  98. D5(3,3)=D5(1,1)*1.D+6
  99. D5(4,4)=2.D0*(1.D0+PX)/EX
  100. C
  101. ENDIF
  102. GO TO 100
  103. C
  104. 30 CONTINUE
  105. IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN
  106. C
  107. C CONTRAINTES PLANES SANS CISAILLEMENTS TRANSVERSAL
  108. C
  109. D5(1,1)=1.D0/EX
  110. D5(1,2)=-PX/EX
  111. D5(2,1)=-PX/EX
  112. D5(2,2)=1.D0/EX
  113. D5(3,3)=2.D0*(1.D0+PX)/EX
  114. C
  115. ELSE IF(IFOUR.EQ.-1) THEN
  116. C
  117. C DEFORMATIONS PLANES SANS CISAILLEMENTS TRANSVERSAL
  118. C
  119. E1=(1.D0+PX)/EX
  120. D5(1,1)=E1*(1.D0-PX)
  121. D5(1,2)=-PX*E1
  122. D5(2,1)=-PX*E1
  123. D5(2,2)=E1*(1.D0-PX)
  124. D5(3,3)=E1*2.D0
  125. C
  126. ELSE
  127. GO TO 1000
  128. ENDIF
  129. GO TO 100
  130. C
  131. 10 CONTINUE
  132. C
  133. C COMPORTEMENT UNIDIMENSIONNEL
  134. C
  135. D5(1,1)=1.D0/EX
  136. GO TO 100
  137. 1000 WRITE(*,200) IFOUR,NSTRS
  138. STOP
  139. 100 CONTINUE
  140. DO 300 I = 1,NSTRS
  141. DO 400 J = 1,NSTRS
  142. D(I,J) = D5(I,J)
  143. 400 CONTINUE
  144. 300 CONTINUE
  145. 200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS MATSOUP',/,
  146. * 10X,'IFOUR=',I2,' NSTRS=',I2,//)
  147. RETURN
  148. END
  149.  
  150.  
  151.  
  152.  

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