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    JK148537  26/08/27    21:15:08     12628          
  2.       SUBROUTINE MAT29(WR10,WRK9,INPLAS,IFOU,MFR,jmat)
  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. * jk148537 2026 a tester ?
  83.           INTMAT = jmat + 1
  84.       ELSE
  85.           INTMAT = jmat + 2
  86.       ENDIF
  87.       IF (IABLO1(/1).EQ.INTMAT) THEN
  88.                             NKX(NNKX)=NYKX/2
  89.                             DO 120 I=1,NYKX
  90.  120                        YKX(I)=TABLO2(IORIG+I)
  91.       ELSE
  92.                             DO 130 I=1,NNKX
  93.  130                        NKX(I)=IABLO1(12+I)
  94.                             DO 140 I=1,NYKX
  95.  140                        YKX(I)=TABLO2(IORIG+I)
  96.       ENDIF
  97.       IORIG=IORIG+NYKX
  98.       DO 150 I=1,YRHO(/1)
  99.  150   YRHO(I)=TABLO2(IORIG+I)
  100.       SIGY(1)=0.D0
  101.       GOTO 500
  102. C
  103. C  cas du materiau plastique endommageable de Lemaitre
  104. C
  105.  200  CONTINUE
  106.       DO 210 I=1,YOG(/1)
  107.  210  YOG(I)=TABLO2(I)
  108.       IORIG=YOG(/1)
  109.       DO 220 I=1,YNU(/1)
  110.  220  YNU(I)=TABLO2(IORIG+I)
  111.       IORIG=IORIG+YNU(/1)
  112.       DO 230 I=1,YRHO(/1)
  113.  230  YRHO(I)=TABLO2(IORIG+I)
  114.       IORIG=IORIG+YRHO(/1)
  115.       DO 240 I=1,YALFA(/1)
  116.  240  YALFA(I)=TABLO2(IORIG+I)
  117.       IORIG=IORIG+YALFA(/1)
  118.       NNKX=NKX(/1)
  119.       NYKX=YKX(/1)
  120.       IF ((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.33.OR.MFR.EQ.5).
  121.      +                        AND.IFOU.EQ.-2) THEN
  122. * jk148537 2026
  123.           INTMAT = jmat
  124.       ELSE
  125.           INTMAT = jmat
  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.  

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