Esope source files

   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 :