Télécharger la plaquette Cast3M


fr en

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