Télécharger mat29.eso

Retour à la liste

Numérotation des lignes :

mat29
  1. C MAT29 SOURCE CB215821 20/07/29 21:15:44 10668
  2. SUBROUTINE MAT29(WR10,WRK9,INPLAS,IFOU,MFR)
  3. C
  4. C pour les materiaux endommageables de Lemaitre dependant ou non de la
  5. C temperature T :
  6. C YOG = tableau du module d'Young
  7. C YNU = tableau du coefficient de Poisson
  8. C YALFA = tableau du coefficient de dilatation thermique
  9. C YRHO = tableau de la masse volumique
  10. C YN, YM, YKK, YALFA1, YBETA1, YR, YA = tableaux des parametres non
  11. C lineaires du materiau 29
  12. C YN, YM, YKK = tableaux des parametres non lineaires du materiau 26
  13. C REMARQUE : si la dimension d'un tableau vaut 1, le parametre
  14. C correspondant est independant de T
  15. C sinon il contient T1 VAL1 T2 VAL2 ....
  16. C YKX = fonction k(X) tabulee en fonction de T pour le materiau 29
  17. C NKX = tableau de pointeurs sur les courbes k(X)
  18. C si la dimension de NKX est 1, k(X) est independant de T,NKX contient
  19. C le nbr. de pts. definissant k(X) et YKX contient X1 k1 X2 k2 ....
  20. C sinon NKX contient pour chaque courbe k(X) le nbr. de pts. qui la
  21. C definit et YKX contient T1 X11 k11 X12 k12 ...T2 X21 k21 X22 k22 ...
  22. C T3 ...
  23. C Pour le materiau 26 idem mais k(X) represente la courbe de traction
  24. C contraintes - deformations
  25. C
  26. C
  27. IMPLICIT INTEGER(I-N)
  28. IMPLICIT REAL*8(A-H,O-Z)
  29. SEGMENT WRK9
  30. REAL*8 YOG(NYOG),YNU(NYNU),YALFA(NYALFA),YSMAX(NYSMAX)
  31. REAL*8 YN(NYN),YM(NYM),YKK(NYKK),YALFA1(NYALF1)
  32. REAL*8 YBETA1(NYBET1),YR(NYR),YA(NYA),YKX(NYKX),YRHO(NYRHO)
  33. REAL*8 SIGY(NSIGY)
  34. INTEGER NKX(NNKX)
  35. ENDSEGMENT
  36. SEGMENT WR10
  37. INTEGER IABLO1(NTABO1)
  38. REAL*8 TABLO2(NTABO2)
  39. ENDSEGMENT
  40. C
  41. IF (INPLAS.EQ.26) GOTO 200
  42. C
  43. C cas du materiau viscoplastique endommageable de Lemaitre
  44. C
  45. DO 10 I=1,YOG(/1)
  46. 10 YOG(I)=TABLO2(I)
  47. IORIG=YOG(/1)
  48. DO 20 I=1,YNU(/1)
  49. 20 YNU(I)=TABLO2(IORIG+I)
  50. IORIG=IORIG+YNU(/1)
  51. DO 30 I=1,YALFA(/1)
  52. 30 YALFA(I)=TABLO2(IORIG+I)
  53. IORIG=IORIG+YALFA(/1)
  54. DO 40 I=1,YSMAX(/1)
  55. 40 YSMAX(I)=TABLO2(IORIG+I)
  56. IORIG=IORIG+YSMAX(/1)
  57. DO 50 I=1,YN(/1)
  58. 50 YN(I)=TABLO2(IORIG+I)
  59. IORIG=IORIG+YN(/1)
  60. DO 60 I=1,YM(/1)
  61. 60 YM(I)=TABLO2(IORIG+I)
  62. IORIG=IORIG+YM(/1)
  63. DO 70 I=1,YKK(/1)
  64. 70 YKK(I)=TABLO2(IORIG+I)
  65. IORIG=IORIG+YKK(/1)
  66. DO 80 I=1,YALFA1(/1)
  67. 80 YALFA1(I)=TABLO2(IORIG+I)
  68. IORIG=IORIG+YALFA1(/1)
  69. DO 90 I=1,YBETA1(/1)
  70. 90 YBETA1(I)=TABLO2(IORIG+I)
  71. IORIG=IORIG+YBETA1(/1)
  72. DO 100 I=1,YR(/1)
  73. 100 YR(I)=TABLO2(IORIG+I)
  74. IORIG=IORIG+YR(/1)
  75. DO 110 I=1,YA(/1)
  76. 110 YA(I)=TABLO2(IORIG+I)
  77. IORIG=IORIG+YA(/1)
  78. NNKX=NKX(/1)
  79. NYKX=YKX(/1)
  80. IF ((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.33.OR.MFR.EQ.5).
  81. + AND.IFOU.EQ.-2) THEN
  82. INTMAT=16
  83. *+TREF
  84. *+TALP
  85. ELSE
  86. INTMAT=15
  87. *+TREF
  88. *+TALP
  89. ENDIF
  90. IF (IABLO1(/1).EQ.INTMAT) THEN
  91. NKX(NNKX)=NYKX/2
  92. DO 120 I=1,NYKX
  93. 120 YKX(I)=TABLO2(IORIG+I)
  94. ELSE
  95. DO 130 I=1,NNKX
  96. 130 NKX(I)=IABLO1(12+I)
  97. DO 140 I=1,NYKX
  98. 140 YKX(I)=TABLO2(IORIG+I)
  99. ENDIF
  100. IORIG=IORIG+NYKX
  101. DO 150 I=1,YRHO(/1)
  102. 150 YRHO(I)=TABLO2(IORIG+I)
  103. SIGY(1)=0.D0
  104. GOTO 500
  105. C
  106. C cas du materiau plastique endommageable de Lemaitre
  107. C
  108. 200 CONTINUE
  109. DO 210 I=1,YOG(/1)
  110. 210 YOG(I)=TABLO2(I)
  111. IORIG=YOG(/1)
  112. DO 220 I=1,YNU(/1)
  113. 220 YNU(I)=TABLO2(IORIG+I)
  114. IORIG=IORIG+YNU(/1)
  115. DO 230 I=1,YRHO(/1)
  116. 230 YRHO(I)=TABLO2(IORIG+I)
  117. IORIG=IORIG+YRHO(/1)
  118. DO 240 I=1,YALFA(/1)
  119. 240 YALFA(I)=TABLO2(IORIG+I)
  120. IORIG=IORIG+YALFA(/1)
  121. NNKX=NKX(/1)
  122. NYKX=YKX(/1)
  123. IF ((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.33.OR.MFR.EQ.5).
  124. + AND.IFOU.EQ.-2) THEN
  125. *+DC INTMAT=9
  126. *+TREF
  127. *+TALP
  128. INTMAT=12
  129. ELSE
  130. *+DC INTMAT=8
  131. *+TREF
  132. *+TALP
  133. INTMAT=11
  134. ENDIF
  135.  
  136. IF (IABLO1(/1).EQ.INTMAT) THEN
  137. NKX(NNKX)=NYKX/2
  138. DO 250 I=1,NYKX
  139. 250 YKX(I)=TABLO2(IORIG+I)
  140. ELSE
  141. DO 260 I=1,NNKX
  142. 260 NKX(I)=IABLO1(5+I)
  143. DO 270 I=1,NYKX
  144. 270 YKX(I)=TABLO2(IORIG+I)
  145. ENDIF
  146. IORIG=IORIG+NYKX
  147. DO 280 I=1,YN(/1)
  148. 280 YN(I)=TABLO2(IORIG+I)
  149. IORIG=IORIG+YN(/1)
  150. DO 290 I=1,YM(/1)
  151. 290 YM(I)=TABLO2(IORIG+I)
  152. IORIG=IORIG+YM(/1)
  153. DO 300 I=1,YKK(/1)
  154. 300 YKK(I)=TABLO2(IORIG+I)
  155. C
  156. 500 CONTINUE
  157. RETURN
  158. END
  159.  
  160.  

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