Télécharger mat29.eso

Retour à la liste

Numérotation des lignes :

  1. C MAT29 SOURCE CHAT 05/01/13 01:33:35 5004
  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=14
  83. ELSE
  84. INTMAT=13
  85. ENDIF
  86. IF (IABLO1(/1).EQ.INTMAT) THEN
  87. NKX(NNKX)=NYKX/2
  88. DO 120 I=1,NYKX
  89. 120 YKX(I)=TABLO2(IORIG+I)
  90. ELSE
  91. DO 130 I=1,NNKX
  92. 130 NKX(I)=IABLO1(12+I)
  93. DO 140 I=1,NYKX
  94. 140 YKX(I)=TABLO2(IORIG+I)
  95. ENDIF
  96. IORIG=IORIG+NYKX
  97. DO 150 I=1,YRHO(/1)
  98. 150 YRHO(I)=TABLO2(IORIG+I)
  99. SIGY(1)=0.D0
  100. GOTO 500
  101. C
  102. C cas du materiau plastique endommageable de Lemaitre
  103. C
  104. 200 CONTINUE
  105. DO 210 I=1,YOG(/1)
  106. 210 YOG(I)=TABLO2(I)
  107. IORIG=YOG(/1)
  108. DO 220 I=1,YNU(/1)
  109. 220 YNU(I)=TABLO2(IORIG+I)
  110. IORIG=IORIG+YNU(/1)
  111. DO 230 I=1,YRHO(/1)
  112. 230 YRHO(I)=TABLO2(IORIG+I)
  113. IORIG=IORIG+YRHO(/1)
  114. DO 240 I=1,YALFA(/1)
  115. 240 YALFA(I)=TABLO2(IORIG+I)
  116. IORIG=IORIG+YALFA(/1)
  117. NNKX=NKX(/1)
  118. NYKX=YKX(/1)
  119. IF ((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.33.OR.MFR.EQ.5).
  120. + AND.IFOU.EQ.-2) THEN
  121. *+DC INTMAT=9
  122. INTMAT=10
  123. ELSE
  124. *+DC INTMAT=8
  125. INTMAT=9
  126. ENDIF
  127. IF (IABLO1(/1).EQ.INTMAT) THEN
  128. NKX(NNKX)=NYKX/2
  129. DO 250 I=1,NYKX
  130. 250 YKX(I)=TABLO2(IORIG+I)
  131. ELSE
  132. DO 260 I=1,NNKX
  133. 260 NKX(I)=IABLO1(5+I)
  134. DO 270 I=1,NYKX
  135. 270 YKX(I)=TABLO2(IORIG+I)
  136. ENDIF
  137. IORIG=IORIG+NYKX
  138. DO 280 I=1,YN(/1)
  139. 280 YN(I)=TABLO2(IORIG+I)
  140. IORIG=IORIG+YN(/1)
  141. DO 290 I=1,YM(/1)
  142. 290 YM(I)=TABLO2(IORIG+I)
  143. IORIG=IORIG+YM(/1)
  144. DO 300 I=1,YKK(/1)
  145. 300 YKK(I)=TABLO2(IORIG+I)
  146. C
  147. 500 CONTINUE
  148. RETURN
  149. END
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  

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