1 : C MODEL2 SOURCE AM 12/02/06 21:16:00 7281 2 : 3 : C======================================================================= 4 : C= Objet MODELE - Cas de la formulation 'MECANIQUE' = 5 : C= = 6 : C= Definition des noms de modeles de MATERIAU = 7 : C= et des noms de types d'ELEMENTS FINIS = 8 : C======================================================================= 9 : 10 : SUBROUTINE MODEL2 (LESPRO,NBPRO,LESTEF,NBTEF,N1MAX,N2MAX) 11 : 12 : IMPLICIT INTEGER(I-N) 13 : IMPLICIT REAL*8 (A-H,O-Z) 14 : 15 : -INC CCOPTIO 16 : 17 : CHARACTER*16 MOMODL(200) 18 : CHARACTER*16 LESPRO(*) 19 : CHARACTER*4 LESTEF(*) 20 : 21 : C 1 - Noms des comportements LINEAIRES 22 : C -------------------------------------- 22 : CALL MODLIN(MOMODL,NMOD) 23 : DO i=1,NMOD 24 : LESPRO(i)=MOMODL(i) 25 : ENDDO 26 : NBPRO=NMOD 27 : 28 : C 2 - Noms des MATERIAUx pour le comportement LINEAIRE 29 : C ------------------------------------------------------ 29 : CALL MODELA(MOMODL,NMOD) 30 : DO i=NBPRO+1,NBPRO+NMOD 31 : LESPRO(i)=MOMODL(i-NBPRO) 32 : ENDDO 33 : NBPRO=NBPRO+NMOD 34 : 35 : C 3 - Noms des comportements NON LINEAIRES 36 : C ------------------------------------------ 36 : CALL MODNLI(MOMODL,NMOD) 37 : DO i=NBPRO+1,NBPRO+NMOD 38 : LESPRO(i)=MOMODL(i-NBPRO) 39 : ENDDO 40 : NBPRO=NBPRO+NMOD 41 : 42 : C 4 - Noms des modeles PLASTIQUES 43 : C --------------------------------- 43 : CALL MODPLA(MOMODL,NMOD) 44 : DO i=NBPRO+1,NBPRO+NMOD 45 : LESPRO(i)=MOMODL(i-NBPRO) 46 : ENDDO 47 : NBPRO=NBPRO+NMOD 48 : 49 : C 5 - Noms des modeles de FLUAGE 50 : C -------------------------------- 50 : CALL MODFLU(MOMODL,NMOD) 51 : DO i=NBPRO+1,NBPRO+NMOD 52 : LESPRO(i)=MOMODL(i-NBPRO) 53 : ENDDO 54 : NBPRO=NBPRO+NMOD 55 : 56 : C 6 - Noms des modeles VISCOPLASTIQUES 57 : C -------------------------------------- 57 : CALL MODVIS(MOMODL,NMOD) 58 : DO i=NBPRO+1,NBPRO+NMOD 59 : LESPRO(i)=MOMODL(i-NBPRO) 60 : ENDDO 61 : NBPRO=NBPRO+NMOD 62 : 63 : C 7 - Noms des modeles ENDOMMAGEABLES 64 : C ------------------------------------- 64 : CALL MODEND(MOMODL,NMOD) 65 : DO i=NBPRO+1,NBPRO+NMOD 66 : LESPRO(i)=MOMODL(i-NBPRO) 67 : ENDDO 68 : NBPRO=NBPRO+NMOD 69 : 70 : C 8 - Noms des modeles PLASTIQUES ENDOMMAGEABLES 71 : C ------------------------------------------------ 71 : CALL MODPLE(MOMODL,NMOD) 72 : DO i=NBPRO+1,NBPRO+NMOD 73 : LESPRO(i)=MOMODL(i-NBPRO) 74 : ENDDO 75 : NBPRO=NBPRO+NMOD 76 : 77 : C 9 - Noms des modeles ELASTIQUES NON LINEAIRES 78 : C ----------------------------------------------- 78 : CALL MODENL(MOMODL,NMOD) 79 : DO i=NBPRO+1,NBPRO+NMOD 80 : LESPRO(i)=MOMODL(i-NBPRO) 81 : ENDDO 82 : NBPRO=NBPRO+NMOD 83 : 84 : C 10 - Noms des modeles VISCOPLASTIQUES EXTERNES 85 : C ------------------------------------------------ 85 : CALL MODVIX(MOMODL,NMOD) 86 : DO i=NBPRO+1,NBPRO+NMOD 87 : LESPRO(i)=MOMODL(i-NBPRO) 88 : ENDDO 89 : NBPRO=NBPRO+NMOD 90 : 91 : C 12 - Noms des modeles LIAISONS 92 : C ------------------------------------------------ 92 : CALL MODLIA(MOMODL,NMOD) 93 : DO i=NBPRO+1,NBPRO+NMOD 94 : LESPRO(i)=MOMODL(i-NBPRO) 95 : ENDDO 96 : NBPRO=NBPRO+NMOD 97 : 98 : C 13 - Noms des modeles IMPEDANCE 99 : C ------------------------------------------------ 99 : CALL MODIMP(MOMODL,NMOD) 100 : DO i=NBPRO+1,NBPRO+NMOD 101 : LESPRO(i)=MOMODL(i-NBPRO) 102 : ENDDO 103 : NBPRO=NBPRO+NMOD 104 : 105 : IF (NBPRO.GT.N1MAX) GOTO 99 106 : 107 : C 11 - Noms de types d'ELEMENTS FINIS 108 : C ------------------------------------- 109 : C DIMENSION 1 : 110 : C --------------- 111 : IF (IDIM.EQ.1) THEN 112 : NBTEF=2 113 : IF (NBTEF.GT.N2MAX) GOTO 99 114 : LESTEF(1)='M1D2' 115 : LESTEF(2)='M1D3' 116 : C DIMENSION 2 : 117 : C --------------- 118 : C Les EF sont definis suivant la valeur de IFOUR 119 : ELSE IF (IDIM.EQ.2) THEN 120 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 121 : NBTEF=17 122 : IF (NBTEF.GT.N2MAX) GOTO 99 123 : LESTEF(1)='TRI3' 124 : LESTEF(2)='TRI6' 125 : LESTEF(3)='QUA4' 126 : LESTEF(4)='QUA8' 127 : LESTEF(5)='COQ2' 128 : LESTEF(6)='BARR' 129 : LESTEF(7)='TRIH' 130 : LESTEF(8)='CERC' 131 : LESTEF(9)='JOI2' 132 : LESTEF(10)='JOI3' 133 : LESTEF(11)='POLY' 134 : LESTEF(12)='BAR3' 135 : LESTEF(13)='QUAH' 136 : LESTEF(14)='Q4RI' 137 : LESTEF(15)='Q8RI' 138 : LESTEF(16)='XQ4R' 139 : LESTEF(17)='ZCO2' 140 : IF (IFOUR.EQ.0) THEN 141 : C Les elements incompressibles ne sont pas utilisables pour le calcul 142 : C en serie de FOURIER 143 : NBTEF=21 144 : IF (NBTEF.GT.N2MAX) GOTO 99 145 : LESTEF(18)='ICT3' 146 : LESTEF(19)='ICQ4' 147 : LESTEF(20)='ICT6' 148 : LESTEF(21)='ICQ8' 149 : ENDIF 150 : ELSE 151 : NBTEF=34 152 : IF (NBTEF.GT.N2MAX) GOTO 99 153 : LESTEF(1)='TRI3' 154 : LESTEF(2)='TRI6' 155 : LESTEF(3)='QUA4' 156 : LESTEF(4)='QUA8' 157 : LESTEF(5)='COQ2' 158 : LESTEF(6)='BARR' 159 : LESTEF(7)='TRIH' 160 : LESTEF(8)='JOI2' 161 : LESTEF(9)='POI1' 162 : LESTEF(10)='JOI3' 163 : LESTEF(11)='ICT3' 164 : LESTEF(12)='ICQ4' 165 : LESTEF(13)='ICT6' 166 : LESTEF(14)='ICQ8' 167 : LESTEF(15)='POLY' 168 : LESTEF(16)='BAR3' 169 : LESTEF(17)='QUAH' 170 : LESTEF(18)='TIMO' 171 : LESTEF(19)='POUT' 172 : LESTEF(20)='TRIS' 173 : LESTEF(21)='QUAS' 174 : LESTEF(22)='POIS' 175 : LESTEF(23)='TRH6' 176 : LESTEF(24)='SEGS' 177 : LESTEF(25)='POJS' 178 : LESTEF(26)='Q4RI' 179 : LESTEF(27)='Q8RI' 180 : LESTEF(28)='JCT3' 181 : LESTEF(29)='JCI4' 182 : LESTEF(30)='JGI2' 183 : LESTEF(31)='CIFL' 184 : LESTEF(32)='XQ4R' 185 : LESTEF(33)='JOI1' 186 : LESTEF(34)='ZCO2' 187 : ENDIF 188 : C DIMENSION 3 : 189 : C --------------- 190 : ELSE IF (IDIM.EQ.3) THEN 191 : NBTEF=43 192 : IF (NBTEF.GT.N2MAX) GOTO 99 193 : LESTEF( 1)='CUB8' 194 : LESTEF( 2)='CU20' 195 : LESTEF( 3)='PRI6' 196 : LESTEF( 4)='PR15' 197 : LESTEF( 5)='TET4' 198 : LESTEF( 6)='TE10' 199 : LESTEF( 7)='PYR5' 200 : LESTEF( 8)='PY13' 201 : LESTEF( 9)='COQ3' 202 : LESTEF(10)='COQ4' 203 : LESTEF(11)='COQ6' 204 : LESTEF(12)='COQ8' 205 : LESTEF(13)='DKT ' 206 : LESTEF(14)='BARR' 207 : LESTEF(15)='POUT' 208 : LESTEF(16)='TUYA' 209 : LESTEF(17)='LISP' 210 : LESTEF(18)='TUFI' 211 : LESTEF(19)='DST ' 212 : LESTEF(20)='JOI4' 213 : LESTEF(21)='TUYO' 214 : LESTEF(22)='TIMO' 215 : LESTEF(23)='TRIS' 216 : LESTEF(24)='QUAS' 217 : LESTEF(25)='POIS' 218 : LESTEF(26)='LISM' 219 : LESTEF(27)='JOT3' 220 : LESTEF(28)='BAR3' 221 : LESTEF(29)='BAEX' 222 : LESTEF(30)='LIA2' 223 : LESTEF(31)='CUBH' 224 : LESTEF(32)='POJS' 225 : LESTEF(33)='Q4RI' 226 : LESTEF(34)='Q8RI' 227 : LESTEF(35)='JGT3' 228 : LESTEF(36)='JGI4' 229 : LESTEF(37)='SEGS' 230 : LESTEF(38)='SHB8' 231 : LESTEF(39)='POI1' 232 : LESTEF(40)='XC8R' 233 : LESTEF(41)='JOI1' 234 : LESTEF(42)='ZCO3' 235 : LESTEF(43)='ZCO4' 236 : ENDIF 237 : RETURN 238 : 239 : C 12 - Traitement des ERREURS 240 : C ----------------------------- 240 : 99 CALL ERREUR(5) 241 : RETURN 242 : 243 : END 244 : 245 : 246 : 247 : 248 : 249 : 250 : 251 : 252 : 253 :
© Cast3M 2003 - All rights reserved.
Disclaimer