1 : C IDVAR4 SOURCE YURI 11/10/28 22:15:22 7181 2 : 3 : C======================================================================= 4 : * RECHERCHE DES NOMS DE COMPOSANTES DES VARIABLES INTERNES 5 : * -------------------------------------------------------- 6 : * 7 : * ENTREES : 8 : * --------- 9 : * MATEPL NUMERO DU MATERIAU PLASTIQUE 10 : * MFR NUMERO DE FORMULATION 11 : * IFOUR VALEUR DE IFOU DANS CCOPTIO 12 : * NPINT NOMBRE DE POINTS D'INTEGRATION DANS L'EPAISSEUR 13 : * POUR LES ELEMENTS DE COQUE INTEGRES 14 : * 15 : * SORTIES : 16 : * --------- 17 : * IPCOMP POINTEUR SUR LA LISTE DES NOMS DE COMPOSANTES 18 : * OBLIGATOIRES ET FACULTATIVES 19 : * NBROBL NOMBRE DE COMPOSANTES OBLIGATOIRES 20 : * NBRFAC NOMBRE DE COMPOSANTES FACTULTATIVES 21 : C======================================================================= 22 : 23 : SUBROUTINE IDVAR4 (MATEPL,MFR,IFOUR,NPINT,IPCOMP,NBROBL,NBRFAC) 24 : 25 : IMPLICIT INTEGER(I-N) 26 : SEGMENT NOMID 27 : CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC) 28 : ENDSEGMENT 29 : 30 : NBROBL=0 31 : NBRFAC=0 32 : NOMID=0 33 : 34 : C 1 - Modele de CHABOCHE et de OHNO 35 : C =================================== 36 : IF (MATEPL.EQ.25) THEN 37 : C ===== 38 : C 1.1 - Elements MASSIFs 39 : C ===== 40 : IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.33) THEN 41 : C= 1.1.1 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou 42 : C= des DEFORMATIONS PLANES GENERALISEES 43 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 44 : NBROBL=20 45 : SEGINI,NOMID 46 : LESOBL( 1)='X1XX' 47 : LESOBL( 2)='X1YY' 48 : LESOBL( 3)='X1ZZ' 49 : LESOBL( 4)='X1XY' 50 : LESOBL( 5)='X2XX' 51 : LESOBL( 6)='X2YY' 52 : LESOBL( 7)='X2ZZ' 53 : LESOBL( 8)='X2XY' 54 : LESOBL( 9)='GPXX' 55 : LESOBL(10)='GPYY' 56 : LESOBL(11)='GPZZ' 57 : LESOBL(12)='GPXY' 58 : LESOBL(13)='VIXX' 59 : LESOBL(14)='VIYY' 60 : LESOBL(15)='VIZZ' 61 : LESOBL(16)='VIXY' 62 : LESOBL(17)='EPSE' 63 : LESOBL(18)='RR ' 64 : LESOBL(19)='QQQ ' 65 : LESOBL(20)='QQ ' 66 : C= 1.1.2 - Cas AXISYMETRIQUE 67 : ELSE IF (IFOUR.EQ.0) THEN 68 : NBROBL=20 69 : SEGINI,NOMID 70 : LESOBL( 1)='X1RR' 71 : LESOBL( 2)='X1ZZ' 72 : LESOBL( 3)='X1TT' 73 : LESOBL( 4)='X1RZ' 74 : LESOBL( 5)='X2RR' 75 : LESOBL( 6)='X2ZZ' 76 : LESOBL( 7)='X2TT' 77 : LESOBL( 8)='X2RZ' 78 : LESOBL( 9)='GPRR' 79 : LESOBL(10)='GPZZ' 80 : LESOBL(11)='GPTT' 81 : LESOBL(12)='GPRZ' 82 : LESOBL(13)='VIRR' 83 : LESOBL(14)='VIZZ' 84 : LESOBL(15)='VITT' 85 : LESOBL(16)='VIRZ' 86 : LESOBL(17)='EPSE' 87 : LESOBL(18)='RR ' 88 : LESOBL(19)='QQQ ' 89 : LESOBL(20)='QQ ' 90 : C= 1.1.3 - Cas des series de FOURIER 91 : ELSE IF (IFOUR.EQ.1) THEN 92 : NBROBL=28 93 : SEGINI,NOMID 94 : LESOBL( 1)='X1RR' 95 : LESOBL( 2)='X1ZZ' 96 : LESOBL( 3)='X1TT' 97 : LESOBL( 4)='X1RZ' 98 : LESOBL( 5)='X1RT' 99 : LESOBL( 6)='X1ZT' 100 : LESOBL( 7)='X2RR' 101 : LESOBL( 8)='X2ZZ' 102 : LESOBL( 9)='X2TT' 103 : LESOBL(10)='X2RZ' 104 : LESOBL(11)='X2RT' 105 : LESOBL(12)='X2ZT' 106 : LESOBL(13)='GPRR' 107 : LESOBL(14)='GPZZ' 108 : LESOBL(15)='GPTT' 109 : LESOBL(16)='GPRZ' 110 : LESOBL(17)='GPRT' 111 : LESOBL(18)='GPZT' 112 : LESOBL(19)='VIRR' 113 : LESOBL(20)='VIZZ' 114 : LESOBL(21)='VITT' 115 : LESOBL(22)='VIRZ' 116 : LESOBL(23)='VIRT' 117 : LESOBL(24)='VIZT' 118 : LESOBL(25)='EPSE' 119 : LESOBL(26)='RR ' 120 : LESOBL(27)='QQQ ' 121 : LESOBL(28)='QQ ' 122 : C= 1.1.4 - Cas TRIDIMENSIONNEL 123 : ELSE IF (IFOUR.EQ.2) THEN 124 : NBROBL=28 125 : SEGINI,NOMID 126 : LESOBL( 1)='X1XX' 127 : LESOBL( 2)='X1YY' 128 : LESOBL( 3)='X1ZZ' 129 : LESOBL( 4)='X1XY' 130 : LESOBL( 5)='X1YZ' 131 : LESOBL( 6)='X1ZX' 132 : LESOBL( 7)='X2XX' 133 : LESOBL( 8)='X2YY' 134 : LESOBL( 9)='X2ZZ' 135 : LESOBL(10)='X2XY' 136 : LESOBL(11)='X2YZ' 137 : LESOBL(12)='X2ZX' 138 : LESOBL(13)='GPXX' 139 : LESOBL(14)='GPYY' 140 : LESOBL(15)='GPZZ' 141 : LESOBL(16)='GPXY' 142 : LESOBL(17)='GPYZ' 143 : LESOBL(18)='GPZX' 144 : LESOBL(19)='VIXX' 145 : LESOBL(20)='VIYY' 146 : LESOBL(21)='VIZZ' 147 : LESOBL(22)='VIXY' 148 : LESOBL(23)='VIYZ' 149 : LESOBL(24)='VIZX' 150 : LESOBL(25)='EPSE' 151 : LESOBL(26)='RR ' 152 : LESOBL(27)='QQQ ' 153 : LESOBL(28)='QQ ' 154 : C= 1.1.5 - Modes de calcul UNIDIMENSIONNELS (1D) PLAN 155 : ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN 156 : NBROBL=16 157 : SEGINI,NOMID 158 : LESOBL( 1)='X1XX' 159 : LESOBL( 2)='X1YY' 160 : LESOBL( 3)='X1ZZ' 161 : LESOBL( 4)='X2XX' 162 : LESOBL( 5)='X2YY' 163 : LESOBL( 6)='X2ZZ' 164 : LESOBL( 7)='GPXX' 165 : LESOBL( 8)='GPYY' 166 : LESOBL( 9)='GPZZ' 167 : LESOBL(10)='VIXX' 168 : LESOBL(11)='VIYY' 169 : LESOBL(12)='VIZZ' 170 : LESOBL(13)='EPSE' 171 : LESOBL(14)='RR ' 172 : LESOBL(15)='QQQ ' 173 : LESOBL(16)='QQ ' 174 : C= 1.1.6 - Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE 175 : ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN 176 : NBROBL=16 177 : SEGINI,NOMID 178 : LESOBL( 1)='X1RR' 179 : LESOBL( 2)='X1ZZ' 180 : LESOBL( 3)='X1TT' 181 : LESOBL( 4)='X2RR' 182 : LESOBL( 5)='X2ZZ' 183 : LESOBL( 6)='X2TT' 184 : LESOBL( 7)='GPRR' 185 : LESOBL( 8)='GPZZ' 186 : LESOBL( 9)='GPTT' 187 : LESOBL(10)='VIRR' 188 : LESOBL(11)='VIZZ' 189 : LESOBL(12)='VITT' 190 : LESOBL(13)='EPSE' 191 : LESOBL(14)='RR ' 192 : LESOBL(15)='QQQ ' 193 : LESOBL(16)='QQ ' 194 : ENDIF 195 : C ===== 196 : C 1.2 - Elements COQUEs 197 : C ===== 198 : ELSE IF (MFR.EQ.3) THEN 199 : IF (NPINT.EQ.0) THEN 200 : C= 1.2.1 - Modes de calcul TRIDIMENSIONNEL et FOURIER 201 : IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN 202 : NBROBL=36 203 : SEGINI,NOMID 204 : LESOBL( 1)='X1N1' 205 : LESOBL( 2)='X1N2' 206 : LESOBL( 3)='X1N3' 207 : LESOBL( 4)='X1N4' 208 : LESOBL( 5)='X1M1' 209 : LESOBL( 6)='X1M2' 210 : LESOBL( 7)='X1M3' 211 : LESOBL( 8)='X1M4' 212 : LESOBL( 9)='X2N1' 213 : LESOBL(10)='X2N2' 214 : LESOBL(11)='X2N3' 215 : LESOBL(12)='X2N4' 216 : LESOBL(13)='X2M1' 217 : LESOBL(14)='X2M2' 218 : LESOBL(15)='X2M3' 219 : LESOBL(16)='X2M4' 220 : LESOBL(17)='GPN1' 221 : LESOBL(18)='GPN2' 222 : LESOBL(19)='GPN3' 223 : LESOBL(20)='GPN4' 224 : LESOBL(21)='GPM1' 225 : LESOBL(22)='GPM2' 226 : LESOBL(23)='GPM3' 227 : LESOBL(24)='GPM4' 228 : LESOBL(25)='EPN1' 229 : LESOBL(26)='EPN2' 230 : LESOBL(27)='EPN3' 231 : LESOBL(28)='EPN4' 232 : LESOBL(29)='EPM1' 233 : LESOBL(30)='EPM2' 234 : LESOBL(31)='EPM3' 235 : LESOBL(32)='EPM4' 236 : LESOBL(33)='EPSE' 237 : LESOBL(34)='RR ' 238 : LESOBL(35)='QQQ ' 239 : LESOBL(36)='QQ ' 240 : C= 1.2.2 - Cas AXISYMETRIQUE 241 : ELSE IF (IFOUR.EQ.0) THEN 242 : NBROBL=28 243 : SEGINI,NOMID 244 : LESOBL( 1)='X1N1' 245 : LESOBL( 2)='X1N2' 246 : LESOBL( 3)='X1N3' 247 : LESOBL( 4)='X1M1' 248 : LESOBL( 5)='X1M2' 249 : LESOBL( 6)='X1M3' 250 : LESOBL( 7)='X2N1' 251 : LESOBL( 8)='X2N2' 252 : LESOBL( 9)='X2N3' 253 : LESOBL(10)='X2M1' 254 : LESOBL(11)='X2M2' 255 : LESOBL(12)='X2M3' 256 : LESOBL(13)='GPN1' 257 : LESOBL(14)='GPN2' 258 : LESOBL(15)='GPN3' 259 : LESOBL(16)='GPM1' 260 : LESOBL(17)='GPM2' 261 : LESOBL(18)='GPM3' 262 : LESOBL(19)='EPN1' 263 : LESOBL(20)='EPN2' 264 : LESOBL(21)='EPN3' 265 : LESOBL(22)='EPM1' 266 : LESOBL(23)='EPM2' 267 : LESOBL(24)='EPM3' 268 : LESOBL(25)='EPSE' 269 : LESOBL(26)='RR ' 270 : LESOBL(27)='QQQ ' 271 : LESOBL(28)='QQ ' 272 : C= 1.2.3 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou 273 : C= des DEFORMATIONS PLANES GENERALISEES 274 : ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-3) THEN 275 : NBROBL=28 276 : SEGINI,NOMID 277 : LESOBL( 1)='X1N1' 278 : LESOBL( 2)='X1N2' 279 : LESOBL( 3)='X1N3' 280 : LESOBL( 4)='X1M1' 281 : LESOBL( 5)='X1M2' 282 : LESOBL( 6)='X1M3' 283 : LESOBL( 7)='X2N1' 284 : LESOBL( 8)='X2N2' 285 : LESOBL( 9)='X2N3' 286 : LESOBL(10)='X2M1' 287 : LESOBL(11)='X2M2' 288 : LESOBL(12)='X2M3' 289 : LESOBL(13)='GPN1' 290 : LESOBL(14)='GPN2' 291 : LESOBL(15)='GPN3' 292 : LESOBL(16)='GPM1' 293 : LESOBL(17)='GPM2' 294 : LESOBL(18)='GPM3' 295 : LESOBL(19)='EPN1' 296 : LESOBL(20)='EPN2' 297 : LESOBL(21)='EPN3' 298 : LESOBL(22)='EPM1' 299 : LESOBL(23)='EPM2' 300 : LESOBL(24)='EPM3' 301 : LESOBL(25)='EPSE' 302 : LESOBL(26)='RR ' 303 : LESOBL(27)='QQQ ' 304 : LESOBL(28)='QQ ' 305 : ENDIF 306 : ELSE 307 : C= 1.2.4 - Cas des COQUES INTEGREES en TRIDIMENSIONNEL 308 : IF (IFOUR.EQ.2) THEN 309 : NBROBL=20 310 : SEGINI,NOMID 311 : LESOBL( 1)='X1SS' 312 : LESOBL( 2)='X1TT' 313 : LESOBL( 3)='X1NN' 314 : LESOBL( 4)='X1ST' 315 : LESOBL( 5)='X2SS' 316 : LESOBL( 6)='X2TT' 317 : LESOBL( 7)='X2NN' 318 : LESOBL( 8)='X2ST' 319 : LESOBL( 9)='GPSS' 320 : LESOBL(10)='GPTT' 321 : LESOBL(11)='GPNN' 322 : LESOBL(12)='GPST' 323 : LESOBL(13)='VISS' 324 : LESOBL(14)='VITT' 325 : LESOBL(15)='VINN' 326 : LESOBL(16)='VIST' 327 : LESOBL(17)='EPSE' 328 : LESOBL(18)='RR ' 329 : LESOBL(19)='QQQ ' 330 : LESOBL(20)='QQ ' 331 : ENDIF 332 : ENDIF 333 : C ===== 334 : C 1.3 - Element COQUE EPAISSE 335 : C ===== 336 : ELSE IF (MFR.EQ.5) THEN 337 : C= Mode de calcul TRIDIMENSIONNEL 338 : IF (IFOUR.EQ.2) THEN 339 : NBROBL=28 340 : SEGINI,NOMID 341 : LESOBL( 1)='X1SS' 342 : LESOBL( 2)='X1TT' 343 : LESOBL( 3)='X1NN' 344 : LESOBL( 4)='X1ST' 345 : LESOBL( 5)='X1SN' 346 : LESOBL( 6)='X1TN' 347 : LESOBL( 7)='X2SS' 348 : LESOBL( 8)='X2TT' 349 : LESOBL( 9)='X2NN' 350 : LESOBL(10)='X2ST' 351 : LESOBL(11)='X2SN' 352 : LESOBL(12)='X2TN' 353 : LESOBL(13)='GPSS' 354 : LESOBL(14)='GPTT' 355 : LESOBL(15)='GPNN' 356 : LESOBL(16)='GPST' 357 : LESOBL(17)='GPSN' 358 : LESOBL(18)='GPTN' 359 : LESOBL(19)='VISS' 360 : LESOBL(20)='VITT' 361 : LESOBL(21)='VINN' 362 : LESOBL(22)='VIST' 363 : LESOBL(23)='VISN' 364 : LESOBL(24)='VITN' 365 : LESOBL(25)='EPSE' 366 : LESOBL(26)='RR ' 367 : LESOBL(27)='QQQ ' 368 : LESOBL(28)='QQ ' 369 : ENDIF 370 : C ===== 371 : C 1.4 - Elements POUTRE et TUYAU 372 : C ===== 373 : ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN 374 : C= Mode de calcul TRIDIMENSIONNEL 375 : IF (IFOUR.EQ.2) THEN 376 : NBROBL=28 377 : SEGINI,NOMID 378 : LESOBL( 1)='X1FX' 379 : LESOBL( 2)='X1FY' 380 : LESOBL( 3)='X1FZ' 381 : LESOBL( 4)='X1MX' 382 : LESOBL( 5)='X1MY' 383 : LESOBL( 6)='X1MZ' 384 : LESOBL( 7)='X2FX' 385 : LESOBL( 8)='X2FY' 386 : LESOBL( 9)='X2FZ' 387 : LESOBL(10)='X2MX' 388 : LESOBL(11)='X2MY' 389 : LESOBL(12)='X2MZ' 390 : LESOBL(13)='GPFX' 391 : LESOBL(14)='GPFY' 392 : LESOBL(15)='GPFZ' 393 : LESOBL(16)='GPMX' 394 : LESOBL(17)='GPMY' 395 : LESOBL(18)='GPMZ' 396 : LESOBL(19)='EPFX' 397 : LESOBL(20)='EPFY' 398 : LESOBL(21)='EPFZ' 399 : LESOBL(22)='EPMX' 400 : LESOBL(23)='EPMY' 401 : LESOBL(24)='EPMZ' 402 : LESOBL(25)='EPSE' 403 : LESOBL(26)='RR ' 404 : LESOBL(27)='QQQ ' 405 : LESOBL(28)='QQ ' 406 : ENDIF 407 : C ===== 408 : C 1.5 - Element COQUE CISAILLEMENT TRANSVERSE 409 : C ===== 410 : ELSE IF (MFR.EQ.9) THEN 411 : C= Mode de calcul TRIDIMENSIONNEL 412 : IF (NPINT.EQ.0) THEN 413 : IF (IFOUR.EQ.2) THEN 414 : NBROBL=36 415 : SEGINI,NOMID 416 : LESOBL( 1)='X1N1' 417 : LESOBL( 2)='X1N2' 418 : LESOBL( 3)='X1N3' 419 : LESOBL( 4)='X1M1' 420 : LESOBL( 5)='X1M2' 421 : LESOBL( 6)='X1M3' 422 : LESOBL( 7)='X1V1' 423 : LESOBL( 8)='X1V2' 424 : LESOBL( 9)='X2N1' 425 : LESOBL(10)='X2N2' 426 : LESOBL(11)='X2N3' 427 : LESOBL(12)='X2M1' 428 : LESOBL(13)='X2M2' 429 : LESOBL(14)='X2M3' 430 : LESOBL(15)='X2V1' 431 : LESOBL(16)='X2V2' 432 : LESOBL(17)='GPN1' 433 : LESOBL(18)='GPN2' 434 : LESOBL(19)='GPN3' 435 : LESOBL(20)='GPM1' 436 : LESOBL(21)='GPM2' 437 : LESOBL(22)='GPM3' 438 : LESOBL(23)='GPV1' 439 : LESOBL(24)='GPV2' 440 : LESOBL(25)='EPN1' 441 : LESOBL(26)='EPN2' 442 : LESOBL(27)='EPN3' 443 : LESOBL(28)='EPM1' 444 : LESOBL(29)='EPM2' 445 : LESOBL(30)='EPM3' 446 : LESOBL(31)='EPV1' 447 : LESOBL(32)='EPV2' 448 : LESOBL(33)='EPSE' 449 : LESOBL(34)='RR ' 450 : LESOBL(35)='QQQ ' 451 : LESOBL(36)='QQ ' 452 : ENDIF 453 : ENDIF 454 : C ===== 455 : C 1.6 - Formulation MEMBRANE 456 : C ===== 457 : ELSE IF (MFR.EQ.25) THEN 458 : NBROBL=16 459 : SEGINI,NOMID 460 : LESOBL( 1)='X1SS' 461 : LESOBL( 2)='X1TT' 462 : LESOBL( 3)='X1ST' 463 : LESOBL( 4)='X2SS' 464 : LESOBL( 5)='X2TT' 465 : LESOBL( 6)='X2ST' 466 : LESOBL( 7)='GPSS' 467 : LESOBL( 8)='GPTT' 468 : LESOBL( 9)='GPST' 469 : LESOBL(10)='VISS' 470 : LESOBL(11)='VITT' 471 : LESOBL(12)='VIST' 472 : LESOBL(13)='EPSE' 473 : LESOBL(14)='RR ' 474 : LESOBL(15)='QQQ ' 475 : LESOBL(16)='QQ ' 476 : C ===== 477 : C 1.7 - Formulation UNIAXIALE 478 : C ===== 479 : ELSE IF (MFR.EQ.27) THEN 480 : NBROBL=8 481 : SEGINI,NOMID 482 : LESOBL(1)='X1FX' 483 : LESOBL(2)='X2FX' 484 : LESOBL(3)='GPFX' 485 : LESOBL(4)='EPFX' 486 : LESOBL(5)='EPSE' 487 : LESOBL(6)='RR ' 488 : LESOBL(7)='QQQ ' 489 : LESOBL(8)='QQ ' 490 : ENDIF 491 : 492 : C 2 - Modele VISCOPLASTIQUE PARFAIT 493 : C =================================== 494 : ELSE IF (MATEPL.EQ.43) THEN 495 : C= Formulation MASSIF : 496 : C= Modele non disponible en PLAN CONT et en FOURIER 497 : IF ( (MFR.EQ.1.OR.MFR.EQ.31) .AND. 498 : & (IFOUR.NE.-2.AND.IFOUR.NE.1) ) THEN 499 : NBROBL=1 500 : SEGINI,NOMID 501 : LESOBL(1)='EPSE' 502 : ENDIF 503 : 504 : C 3 - Modele a deux deformations inelastiques (DDI) 505 : C =================================================== 506 : ELSE IF (MATEPL.EQ.63) THEN 507 : C ===== 508 : C 3.1 - Elements MASSIFs 509 : C ===== 510 : IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.33) THEN 511 : C= 3.1.1 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou 512 : C= des DEFORMATIONS PLANES GENERALISEES 513 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 514 : NBROBL=45 515 : SEGINI,NOMID 516 : LESOBL( 1)='X1XX' 517 : LESOBL( 2)='X1YY' 518 : LESOBL( 3)='X1ZZ' 519 : LESOBL( 4)='X1XY' 520 : LESOBL( 5)='X2XX' 521 : LESOBL( 6)='X2YY' 522 : LESOBL( 7)='X2ZZ' 523 : LESOBL( 8)='X2XY' 524 : LESOBL( 9)='A1XX' 525 : LESOBL(10)='A1YY' 526 : LESOBL(11)='A1ZZ' 527 : LESOBL(12)='A1XY' 528 : LESOBL(13)='A2XX' 529 : LESOBL(14)='A2YY' 530 : LESOBL(15)='A2ZZ' 531 : LESOBL(16)='A2XY' 532 : LESOBL(17)='VIXX' 533 : LESOBL(18)='VIYY' 534 : LESOBL(19)='VIZZ' 535 : LESOBL(20)='VIXY' 536 : LESOBL(21)='Y1XX' 537 : LESOBL(22)='Y1YY' 538 : LESOBL(23)='Y1ZZ' 539 : LESOBL(24)='Y1XY' 540 : LESOBL(25)='Y2XX' 541 : LESOBL(26)='Y2YY' 542 : LESOBL(27)='Y2ZZ' 543 : LESOBL(28)='Y2XY' 544 : LESOBL(29)='B1XX' 545 : LESOBL(30)='B1YY' 546 : LESOBL(31)='B1ZZ' 547 : LESOBL(32)='B1XY' 548 : LESOBL(33)='B2XX' 549 : LESOBL(34)='B2YY' 550 : LESOBL(35)='B2ZZ' 551 : LESOBL(36)='B2XY' 552 : LESOBL(37)='EVXX' 553 : LESOBL(38)='EVYY' 554 : LESOBL(39)='EVZZ' 555 : LESOBL(40)='EVXY' 556 : LESOBL(41)='VIPW' 557 : LESOBL(42)='RP ' 558 : LESOBL(43)='V ' 559 : LESOBL(44)='RV ' 560 : LESOBL(45)='EPSE' 561 : C= 3.1.2 - Cas AXISYMETRIQUE 562 : ELSE IF (IFOUR.EQ.0) THEN 563 : NBROBL=45 564 : SEGINI,NOMID 565 : LESOBL( 1)='X1RR' 566 : LESOBL( 2)='X1ZZ' 567 : LESOBL( 3)='X1TT' 568 : LESOBL( 4)='X1RZ' 569 : LESOBL( 5)='X2RR' 570 : LESOBL( 6)='X2ZZ' 571 : LESOBL( 7)='X2TT' 572 : LESOBL( 8)='X2RZ' 573 : LESOBL( 9)='A1RR' 574 : LESOBL(10)='A1ZZ' 575 : LESOBL(11)='A1TT' 576 : LESOBL(12)='A1RZ' 577 : LESOBL(13)='A2RR' 578 : LESOBL(14)='A2ZZ' 579 : LESOBL(15)='A2TT' 580 : LESOBL(16)='A2RZ' 581 : LESOBL(17)='VIRR' 582 : LESOBL(18)='VIZZ' 583 : LESOBL(19)='VITT' 584 : LESOBL(20)='VIRZ' 585 : LESOBL(21)='Y1RR' 586 : LESOBL(22)='Y1ZZ' 587 : LESOBL(23)='Y1TT' 588 : LESOBL(24)='Y1RZ' 589 : LESOBL(25)='Y2RR' 590 : LESOBL(26)='Y2ZZ' 591 : LESOBL(27)='Y2TT' 592 : LESOBL(28)='Y2RZ' 593 : LESOBL(29)='B1RR' 594 : LESOBL(30)='B1ZZ' 595 : LESOBL(31)='B1TT' 596 : LESOBL(32)='B1RZ' 597 : LESOBL(33)='B2RR' 598 : LESOBL(34)='B2ZZ' 599 : LESOBL(35)='B2TT' 600 : LESOBL(36)='B2RZ' 601 : LESOBL(37)='EVRR' 602 : LESOBL(38)='EVZZ' 603 : LESOBL(39)='EVTT' 604 : LESOBL(40)='EVRZ' 605 : LESOBL(41)='VIPW' 606 : LESOBL(42)='RP ' 607 : LESOBL(43)='V ' 608 : LESOBL(44)='RV ' 609 : LESOBL(45)='EPSE' 610 : C= 3.1.3 - Cas des series de FOURIER 611 : ELSE IF (IFOUR.EQ.1) THEN 612 : NBROBL=65 613 : SEGINI,NOMID 614 : LESOBL( 1)='X1RR' 615 : LESOBL( 2)='X1ZZ' 616 : LESOBL( 3)='X1TT' 617 : LESOBL( 4)='X1RZ' 618 : LESOBL( 5)='X1RT' 619 : LESOBL( 6)='X1ZT' 620 : LESOBL( 7)='X2RR' 621 : LESOBL( 8)='X2ZZ' 622 : LESOBL( 9)='X2TT' 623 : LESOBL(10)='X2RZ' 624 : LESOBL(11)='X2RT' 625 : LESOBL(12)='X2ZT' 626 : LESOBL(13)='A1RR' 627 : LESOBL(14)='A1ZZ' 628 : LESOBL(15)='A1TT' 629 : LESOBL(16)='A1RZ' 630 : LESOBL(17)='A1RT' 631 : LESOBL(18)='A1ZT' 632 : LESOBL(19)='A2RR' 633 : LESOBL(20)='A2ZZ' 634 : LESOBL(21)='A2TT' 635 : LESOBL(22)='A2RZ' 636 : LESOBL(23)='A2RT' 637 : LESOBL(24)='A2ZT' 638 : LESOBL(25)='VIRR' 639 : LESOBL(26)='VIZZ' 640 : LESOBL(27)='VITT' 641 : LESOBL(28)='VIRZ' 642 : LESOBL(29)='VIRT' 643 : LESOBL(30)='VIZT' 644 : LESOBL(31)='Y1RR' 645 : LESOBL(32)='Y1ZZ' 646 : LESOBL(33)='Y1TT' 647 : LESOBL(34)='Y1RZ' 648 : LESOBL(35)='Y1RT' 649 : LESOBL(36)='Y1ZT' 650 : LESOBL(37)='Y2RR' 651 : LESOBL(38)='Y2ZZ' 652 : LESOBL(39)='Y2TT' 653 : LESOBL(40)='Y2RZ' 654 : LESOBL(41)='Y2RT' 655 : LESOBL(42)='Y2ZT' 656 : LESOBL(43)='B1RR' 657 : LESOBL(44)='B1ZZ' 658 : LESOBL(45)='B1TT' 659 : LESOBL(46)='B1RZ' 660 : LESOBL(47)='B1RT' 661 : LESOBL(48)='B1ZT' 662 : LESOBL(49)='B2RR' 663 : LESOBL(50)='B2ZZ' 664 : LESOBL(51)='B2TT' 665 : LESOBL(52)='B2RZ' 666 : LESOBL(53)='B2RT' 667 : LESOBL(54)='B2ZT' 668 : LESOBL(55)='EVRR' 669 : LESOBL(56)='EVZZ' 670 : LESOBL(57)='EVTT' 671 : LESOBL(58)='EVRZ' 672 : LESOBL(59)='EVRT' 673 : LESOBL(60)='EVZT' 674 : LESOBL(61)='VIPW' 675 : LESOBL(62)='RP ' 676 : LESOBL(63)='V ' 677 : LESOBL(64)='RV ' 678 : LESOBL(65)='EPSE' 679 : C= 3.1.4 - Cas TRIDIMENSIONNEL 680 : ELSE IF (IFOUR.EQ.2) THEN 681 : NBROBL=65 682 : SEGINI,NOMID 683 : LESOBL( 1)='X1XX' 684 : lESOBL( 2)='X1YY' 685 : LESOBL( 3)='X1ZZ' 686 : LESOBL( 4)='X1XY' 687 : LESOBL( 5)='X1YZ' 688 : LESOBL( 6)='X1ZX' 689 : LESOBL( 7)='X2XX' 690 : LESOBL( 8)='X2YY' 691 : LESOBL( 9)='X2ZZ' 692 : LESOBL(10)='X2XY' 693 : LESOBL(11)='X2YZ' 694 : LESOBL(12)='X2ZX' 695 : LESOBL(13)='A1XX' 696 : LESOBL(14)='A1YY' 697 : LESOBL(15)='A1ZZ' 698 : LESOBL(16)='A1XY' 699 : LESOBL(17)='A1YZ' 700 : LESOBL(18)='A1ZX' 701 : LESOBL(19)='A2XX' 702 : LESOBL(20)='A2YY' 703 : LESOBL(21)='A2ZZ' 704 : LESOBL(22)='A2XY' 705 : LESOBL(23)='A2YZ' 706 : LESOBL(24)='A2ZX' 707 : LESOBL(25)='VIXX' 708 : LESOBL(26)='VIYY' 709 : LESOBL(27)='VIZZ' 710 : LESOBL(28)='VIXY' 711 : LESOBL(29)='VIYZ' 712 : LESOBL(30)='VIZX' 713 : LESOBL(31)='Y1XX' 714 : LESOBL(32)='Y1YY' 715 : LESOBL(33)='Y1ZZ' 716 : LESOBL(34)='Y1XY' 717 : LESOBL(35)='Y1YZ' 718 : LESOBL(36)='Y1ZX' 719 : LESOBL(37)='Y2XX' 720 : LESOBL(38)='Y2YY' 721 : LESOBL(39)='Y2ZZ' 722 : LESOBL(40)='Y2XY' 723 : LESOBL(41)='Y2YZ' 724 : LESOBL(42)='Y2ZX' 725 : LESOBL(43)='B1XX' 726 : LESOBL(44)='B1YY' 727 : LESOBL(45)='B1ZZ' 728 : LESOBL(46)='B1XY' 729 : LESOBL(47)='B1YZ' 730 : LESOBL(48)='B1ZX' 731 : LESOBL(49)='B2XX' 732 : LESOBL(50)='B2YY' 733 : LESOBL(51)='B2ZZ' 734 : LESOBL(52)='B2XY' 735 : LESOBL(53)='B2YZ' 736 : LESOBL(54)='B2ZX' 737 : LESOBL(55)='EVXX' 738 : LESOBL(56)='EVYY' 739 : LESOBL(57)='EVZZ' 740 : LESOBL(58)='EVXY' 741 : LESOBL(59)='EVYZ' 742 : LESOBL(60)='EVZX' 743 : LESOBL(61)='VIPW' 744 : LESOBL(62)='RP ' 745 : LESOBL(63)='V ' 746 : LESOBL(64)='RV ' 747 : LESOBL(65)='EPSE' 748 : C= 3.1.5 - Modes de calcul UNIDIMENSIONNELS (1D) PLAN 749 : ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN 750 : NBROBL=35 751 : SEGINI,NOMID 752 : LESOBL( 1)='X1XX' 753 : LESOBL( 2)='X1YY' 754 : LESOBL( 3)='X1ZZ' 755 : LESOBL( 4)='X2XX' 756 : LESOBL( 5)='X2YY' 757 : LESOBL( 6)='X2ZZ' 758 : LESOBL( 7)='A1XX' 759 : LESOBL( 8)='A1YY' 760 : LESOBL( 9)='A1ZZ' 761 : LESOBL(10)='A2XX' 762 : LESOBL(11)='A2YY' 763 : LESOBL(12)='A2ZZ' 764 : LESOBL(13)='VIXX' 765 : LESOBL(14)='VIYY' 766 : LESOBL(15)='VIZZ' 767 : LESOBL(16)='Y1XX' 768 : LESOBL(17)='Y1YY' 769 : LESOBL(18)='Y1ZZ' 770 : LESOBL(19)='Y2XX' 771 : LESOBL(20)='Y2YY' 772 : LESOBL(21)='Y2ZZ' 773 : LESOBL(22)='B1XX' 774 : LESOBL(23)='B1YY' 775 : LESOBL(24)='B1ZZ' 776 : LESOBL(25)='B2XX' 777 : LESOBL(26)='B2YY' 778 : LESOBL(27)='B2ZZ' 779 : LESOBL(28)='EVXX' 780 : LESOBL(29)='EVYY' 781 : LESOBL(30)='EVZZ' 782 : LESOBL(31)='VIPW' 783 : LESOBL(32)='RP ' 784 : LESOBL(33)='V ' 785 : LESOBL(34)='RV ' 786 : LESOBL(35)='EPSE' 787 : C= 3.1.6 - Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE 788 : ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN 789 : NBROBL=35 790 : SEGINI,NOMID 791 : LESOBL( 1)='X1RR' 792 : LESOBL( 2)='X1ZZ' 793 : LESOBL( 3)='X1TT' 794 : LESOBL( 4)='X2RR' 795 : LESOBL( 5)='X2ZZ' 796 : LESOBL( 6)='X2TT' 797 : LESOBL( 7)='A1RR' 798 : LESOBL( 8)='A1ZZ' 799 : LESOBL( 9)='A1TT' 800 : LESOBL(10)='A2RR' 801 : LESOBL(11)='A2ZZ' 802 : LESOBL(12)='A2TT' 803 : LESOBL(13)='VIRR' 804 : LESOBL(14)='VIZZ' 805 : LESOBL(15)='VITT' 806 : LESOBL(16)='Y1RR' 807 : LESOBL(17)='Y1ZZ' 808 : LESOBL(18)='Y1TT' 809 : LESOBL(19)='Y2RR' 810 : LESOBL(20)='Y2ZZ' 811 : LESOBL(21)='Y2TT' 812 : LESOBL(22)='B1RR' 813 : LESOBL(23)='B1ZZ' 814 : LESOBL(24)='B1TT' 815 : LESOBL(25)='B2RR' 816 : LESOBL(26)='B2ZZ' 817 : LESOBL(27)='B2TT' 818 : LESOBL(28)='EVRR' 819 : LESOBL(29)='EVZZ' 820 : LESOBL(30)='EVTT' 821 : LESOBL(31)='VIPW' 822 : LESOBL(32)='RP ' 823 : LESOBL(33)='V ' 824 : LESOBL(34)='RV ' 825 : LESOBL(35)='EPSE' 826 : ENDIF 827 : C ===== 828 : C 3.2 - Elements COQUEs 829 : C ===== 830 : ELSE IF (MFR.EQ.3) THEN 831 : IF (NPINT.EQ.0) THEN 832 : C= 3.2.1 - Modes de calcul TRIDIMENSIONNEL et FOURIER 833 : IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN 834 : NBROBL=85 835 : SEGINI,NOMID 836 : LESOBL( 1)='X1N1' 837 : LESOBL( 2)='X1N2' 838 : LESOBL( 3)='X1N3' 839 : LESOBL( 4)='X1N4' 840 : LESOBL( 5)='X1M1' 841 : LESOBL( 6)='X1M2' 842 : LESOBL( 7)='X1M3' 843 : LESOBL( 8)='X1M4' 844 : LESOBL( 9)='X2N1' 845 : LESOBL(10)='X2N2' 846 : LESOBL(11)='X2N3' 847 : LESOBL(12)='X2N4' 848 : LESOBL(13)='X2M1' 849 : LESOBL(14)='X2M2' 850 : LESOBL(15)='X2M3' 851 : LESOBL(16)='X2M4' 852 : LESOBL(17)='A1N1' 853 : LESOBL(18)='A1N2' 854 : LESOBL(19)='A1N3' 855 : LESOBL(20)='A1N4' 856 : LESOBL(21)='A1M1' 857 : LESOBL(22)='A1M2' 858 : LESOBL(23)='A1M3' 859 : LESOBL(24)='A1M4' 860 : LESOBL(25)='A2N1' 861 : LESOBL(26)='A2N2' 862 : LESOBL(27)='A2N3' 863 : LESOBL(28)='A2N4' 864 : LESOBL(29)='A2M1' 865 : LESOBL(30)='A2M2' 866 : LESOBL(31)='A2M3' 867 : LESOBL(32)='A2M4' 868 : LESOBL(33)='EPN1' 869 : LESOBL(34)='EPN2' 870 : LESOBL(35)='EPN3' 871 : LESOBL(36)='EPN4' 872 : LESOBL(37)='EPM1' 873 : LESOBL(38)='EPM2' 874 : LESOBL(39)='EPM3' 875 : LESOBL(40)='EPM4' 876 : LESOBL(41)='Y1N1' 877 : LESOBL(42)='Y1N2' 878 : LESOBL(43)='Y1N3' 879 : LESOBL(44)='Y1N4' 880 : LESOBL(45)='Y1M1' 881 : LESOBL(46)='Y1M2' 882 : LESOBL(47)='Y1M3' 883 : LESOBL(48)='Y1M4' 884 : LESOBL(49)='Y2N1' 885 : LESOBL(50)='Y2N2' 886 : LESOBL(51)='Y2N3' 887 : LESOBL(52)='Y2N4' 888 : LESOBL(53)='Y2M1' 889 : LESOBL(54)='Y2M2' 890 : LESOBL(55)='Y2M3' 891 : LESOBL(56)='Y2M4' 892 : LESOBL(57)='B1N1' 893 : LESOBL(58)='B1N2' 894 : LESOBL(59)='B1N3' 895 : LESOBL(60)='B1N4' 896 : LESOBL(61)='B1M1' 897 : LESOBL(62)='B1M2' 898 : LESOBL(63)='B1M3' 899 : LESOBL(64)='B1M4' 900 : LESOBL(65)='B2N1' 901 : LESOBL(66)='B2N2' 902 : LESOBL(67)='B2N3' 903 : LESOBL(68)='B2N4' 904 : LESOBL(69)='B2M1' 905 : LESOBL(70)='B2M2' 906 : LESOBL(71)='B2M3' 907 : LESOBL(72)='A2M4' 908 : LESOBL(73)='EVN1' 909 : LESOBL(74)='EVN2' 910 : LESOBL(75)='EVN3' 911 : LESOBL(76)='EVN4' 912 : LESOBL(77)='EVM1' 913 : LESOBL(78)='EVM2' 914 : LESOBL(79)='EVM3' 915 : LESOBL(80)='EVM4' 916 : LESOBL(81)='VIPW' 917 : LESOBL(82)='RP ' 918 : LESOBL(83)='V ' 919 : LESOBL(84)='RV ' 920 : LESOBL(85)='EPSE' 921 : C= 3.2.2 - Cas AXISYMETRIQUE 922 : ELSE IF (IFOUR.EQ.0) THEN 923 : NBROBL=65 924 : SEGINI,NOMID 925 : LESOBL( 1)='X1N1' 926 : LESOBL( 2)='X1N2' 927 : LESOBL( 3)='X1N3' 928 : LESOBL( 4)='X1M1' 929 : LESOBL( 5)='X1M2' 930 : LESOBL( 6)='X1M3' 931 : LESOBL( 7)='X2N1' 932 : LESOBL( 8)='X2N2' 933 : LESOBL( 9)='X2N3' 934 : LESOBL(10)='X2M1' 935 : LESOBL(11)='X2M2' 936 : LESOBL(12)='X2M3' 937 : LESOBL(13)='A1N1' 938 : LESOBL(14)='A1N2' 939 : LESOBL(15)='A1N3' 940 : LESOBL(16)='A1M1' 941 : LESOBL(17)='A1M2' 942 : LESOBL(18)='A1M3' 943 : LESOBL(19)='A2N1' 944 : LESOBL(20)='A2N2' 945 : LESOBL(21)='A2N3' 946 : LESOBL(22)='A2M1' 947 : LESOBL(23)='A2M2' 948 : LESOBL(24)='A2M3' 949 : LESOBL(25)='EPN1' 950 : LESOBL(26)='EPN2' 951 : LESOBL(27)='EPN3' 952 : LESOBL(28)='EPM1' 953 : LESOBL(29)='EPM2' 954 : LESOBL(30)='EPM3' 955 : LESOBL(31)='Y1N1' 956 : LESOBL(32)='Y1N2' 957 : LESOBL(33)='Y1N3' 958 : LESOBL(34)='Y1M1' 959 : LESOBL(35)='Y1M2' 960 : LESOBL(36)='Y1M3' 961 : LESOBL(37)='Y2N1' 962 : LESOBL(38)='Y2N2' 963 : LESOBL(39)='Y2N3' 964 : LESOBL(40)='Y2M1' 965 : LESOBL(41)='Y2M2' 966 : LESOBL(42)='Y2M3' 967 : LESOBL(43)='B1N1' 968 : LESOBL(44)='B1N2' 969 : LESOBL(45)='B1N3' 970 : LESOBL(46)='B1M1' 971 : LESOBL(47)='B1M2' 972 : LESOBL(48)='B1M3' 973 : LESOBL(49)='B2N1' 974 : LESOBL(50)='B2N2' 975 : LESOBL(51)='B2N3' 976 : LESOBL(52)='B2M1' 977 : LESOBL(53)='B2M2' 978 : LESOBL(54)='B2M3' 979 : LESOBL(55)='EVN1' 980 : LESOBL(56)='EVN2' 981 : LESOBL(57)='EVN3' 982 : LESOBL(58)='EVM1' 983 : LESOBL(59)='EVM2' 984 : LESOBL(60)='EVM3' 985 : LESOBL(61)='VIPW' 986 : LESOBL(62)='RP ' 987 : LESOBL(63)='V ' 988 : LESOBL(64)='RV ' 989 : LESOBL(65)='EPSE' 990 : C= 3.2.3 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou 991 : C= des DEFORMATIONS PLANES GENERALISEES 992 : ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-3) THEN 993 : NBROBL=65 994 : SEGINI,NOMID 995 : LESOBL( 1)='X1N1' 996 : LESOBL( 2)='X1N2' 997 : LESOBL( 3)='X1N3' 998 : LESOBL( 4)='X1M1' 999 : LESOBL( 5)='X1M2' 1000 : LESOBL( 6)='X1M3' 1001 : LESOBL( 7)='X2N1' 1002 : LESOBL( 8)='X2N2' 1003 : LESOBL( 9)='X2N3' 1004 : LESOBL(10)='X2M1' 1005 : LESOBL(11)='X2M2' 1006 : LESOBL(12)='X2M3' 1007 : LESOBL(13)='A1N1' 1008 : LESOBL(14)='A1N2' 1009 : LESOBL(15)='A1N3' 1010 : LESOBL(16)='A1M1' 1011 : LESOBL(17)='A1M2' 1012 : LESOBL(18)='A1M3' 1013 : LESOBL(19)='A2N1' 1014 : LESOBL(20)='A2N2' 1015 : LESOBL(21)='A2N3' 1016 : LESOBL(22)='A2M1' 1017 : LESOBL(23)='A2M2' 1018 : LESOBL(24)='A2M3' 1019 : LESOBL(25)='EPN1' 1020 : LESOBL(26)='EPN2' 1021 : LESOBL(27)='EPN3' 1022 : LESOBL(28)='EPM1' 1023 : LESOBL(29)='EPM2' 1024 : LESOBL(30)='EPM3' 1025 : LESOBL(31)='Y1N1' 1026 : LESOBL(32)='Y1N2' 1027 : LESOBL(33)='Y1N3' 1028 : LESOBL(34)='Y1M1' 1029 : LESOBL(35)='Y1M2' 1030 : LESOBL(36)='Y1M3' 1031 : LESOBL(37)='Y2N1' 1032 : LESOBL(38)='Y2N2' 1033 : LESOBL(39)='Y2N3' 1034 : LESOBL(40)='Y2M1' 1035 : LESOBL(41)='Y2M2' 1036 : LESOBL(42)='Y2M3' 1037 : LESOBL(43)='B1N1' 1038 : LESOBL(44)='B1N2' 1039 : LESOBL(45)='B1N3' 1040 : LESOBL(46)='B1M1' 1041 : LESOBL(47)='B1M2' 1042 : LESOBL(48)='B1M3' 1043 : LESOBL(49)='B2N1' 1044 : LESOBL(50)='B2N2' 1045 : LESOBL(51)='B2N3' 1046 : LESOBL(52)='B2M1' 1047 : LESOBL(53)='B2M2' 1048 : LESOBL(54)='B2M3' 1049 : LESOBL(55)='EVN1' 1050 : LESOBL(56)='EVN2' 1051 : LESOBL(57)='EVN3' 1052 : LESOBL(58)='EVM1' 1053 : LESOBL(59)='EVM2' 1054 : LESOBL(60)='EVM3' 1055 : LESOBL(61)='VIPW' 1056 : LESOBL(62)='RP ' 1057 : LESOBL(63)='V ' 1058 : LESOBL(64)='RV ' 1059 : LESOBL(65)='EPSE' 1060 : ENDIF 1061 : ELSE 1062 : C= 3.2.4 - Cas des COQUES INTEGREES en TRIDIMENSIONNEL 1063 : IF (IFOUR.EQ.2) THEN 1064 : NBROBL=45 1065 : SEGINI,NOMID 1066 : LESOBL( 1)='X1SS' 1067 : LESOBL( 2)='X1TT' 1068 : LESOBL( 3)='X1NN' 1069 : LESOBL( 4)='X1ST' 1070 : LESOBL( 5)='X2SS' 1071 : LESOBL( 6)='X2TT' 1072 : LESOBL( 7)='X2NN' 1073 : LESOBL( 8)='X2ST' 1074 : LESOBL( 9)='A1SS' 1075 : LESOBL(10)='A1TT' 1076 : LESOBL(11)='A1NN' 1077 : LESOBL(12)='A1ST' 1078 : LESOBL(13)='A2SS' 1079 : LESOBL(14)='A2TT' 1080 : LESOBL(15)='A2NN' 1081 : LESOBL(16)='A2ST' 1082 : LESOBL(17)='VISS' 1083 : LESOBL(18)='VITT' 1084 : LESOBL(19)='VINN' 1085 : LESOBL(20)='VIST' 1086 : LESOBL(21)='Y1SS' 1087 : LESOBL(22)='Y1TT' 1088 : LESOBL(23)='Y1NN' 1089 : LESOBL(24)='Y1ST' 1090 : LESOBL(25)='Y2SS' 1091 : LESOBL(26)='Y2TT' 1092 : LESOBL(27)='Y2NN' 1093 : LESOBL(28)='Y2ST' 1094 : LESOBL(29)='B1SS' 1095 : LESOBL(30)='B1TT' 1096 : LESOBL(31)='B1NN' 1097 : LESOBL(32)='B1ST' 1098 : LESOBL(33)='B2SS' 1099 : LESOBL(34)='B2TT' 1100 : LESOBL(35)='B2NN' 1101 : LESOBL(36)='B2ST' 1102 : LESOBL(37)='EVSS' 1103 : LESOBL(38)='EVTT' 1104 : LESOBL(39)='EVNN' 1105 : LESOBL(40)='EVST' 1106 : LESOBL(41)='VIPW ' 1107 : LESOBL(42)='RP ' 1108 : LESOBL(43)='V ' 1109 : LESOBL(44)='RV ' 1110 : LESOBL(45)='EPSE' 1111 : ENDIF 1112 : ENDIF 1113 : C ===== 1114 : C 3.3 - Element COQUE EPAISSE 1115 : C ===== 1116 : ELSE IF (MFR.EQ.5) THEN 1117 : C= Mode de calcul TRIDIMENSIONNEL 1118 : NBROBL=65 1119 : SEGINI,NOMID 1120 : LESOBL( 1)='X1SS' 1121 : LESOBL( 2)='X1TT' 1122 : LESOBL( 3)='X1NN' 1123 : LESOBL( 4)='X1ST' 1124 : LESOBL( 5)='X1SN' 1125 : LESOBL( 6)='X1TN' 1126 : LESOBL( 7)='X2SS' 1127 : LESOBL( 8)='X2TT' 1128 : LESOBL( 9)='X2NN' 1129 : LESOBL(10)='X2ST' 1130 : LESOBL(11)='X2SN' 1131 : LESOBL(12)='X2TN' 1132 : LESOBL(13)='A1SS' 1133 : LESOBL(14)='A1TT' 1134 : LESOBL(15)='A1NN' 1135 : LESOBL(16)='A1ST' 1136 : LESOBL(17)='A1SN' 1137 : LESOBL(18)='A1TN' 1138 : LESOBL(19)='A2SS' 1139 : LESOBL(20)='A2TT' 1140 : LESOBL(21)='A2NN' 1141 : LESOBL(22)='A2ST' 1142 : LESOBL(23)='A2SN' 1143 : LESOBL(24)='A2TN' 1144 : LESOBL(25)='VISS' 1145 : LESOBL(26)='VITT' 1146 : LESOBL(27)='VINN' 1147 : LESOBL(28)='VIST' 1148 : LESOBL(29)='VISN' 1149 : LESOBL(30)='VITN' 1150 : LESOBL(31)='Y1SS' 1151 : LESOBL(32)='Y1TT' 1152 : LESOBL(33)='Y1NN' 1153 : LESOBL(34)='Y1ST' 1154 : LESOBL(35)='Y1SN' 1155 : LESOBL(36)='Y1TN' 1156 : LESOBL(37)='Y2SS' 1157 : LESOBL(38)='Y2TT' 1158 : LESOBL(39)='Y2NN' 1159 : LESOBL(40)='Y2ST' 1160 : LESOBL(41)='Y2SN' 1161 : LESOBL(42)='Y2TN' 1162 : LESOBL(43)='B1SS' 1163 : LESOBL(44)='B1TT' 1164 : LESOBL(45)='B1NN' 1165 : LESOBL(46)='B1ST' 1166 : LESOBL(47)='B1SN' 1167 : LESOBL(48)='B1TN' 1168 : LESOBL(49)='B2SS' 1169 : LESOBL(50)='B2TT' 1170 : LESOBL(51)='B2NN' 1171 : LESOBL(52)='B2ST' 1172 : LESOBL(53)='B2SN' 1173 : LESOBL(54)='B2TN' 1174 : LESOBL(55)='EVSS' 1175 : LESOBL(56)='EVTT' 1176 : LESOBL(57)='EVNN' 1177 : LESOBL(58)='EVST' 1178 : LESOBL(59)='EVSN' 1179 : LESOBL(60)='EVTN' 1180 : LESOBL(61)='VIPW' 1181 : LESOBL(62)='RP ' 1182 : LESOBL(63)='V ' 1183 : LESOBL(64)='RV ' 1184 : LESOBL(65)='EPSE' 1185 : ENDIF 1186 : 1187 : C 4 - Modele ELASTO-VISCOPLASTIQUE DE KOCKS 1188 : C =========================================== 1189 : ELSE IF (MATEPL.EQ.70) THEN 1190 : C= Formulation MASSIF : 1191 : C= Modele non disponible en PLAN CONT et en FOURIER 1192 : IF (MFR.EQ.1.AND.(IFOUR.NE.-2.AND.IFOUR.NE.1)) THEN 1193 : NBROBL=2 1194 : SEGINI,NOMID 1195 : LESOBL(1)='EPSE' 1196 : LESOBL(2)='S ' 1197 : ENDIF 1198 : 1199 : C 5 - Modele MISTRAL 1200 : C ==================== 1201 : ELSE IF (MATEPL.EQ.94) THEN 1202 : C ===== 1203 : C 5.1 - Elements MASSIFs (seule formulation disponible actuellement) 1204 : C ===== 1205 : C= Note : Les indices finaux 1, 2, 3, 4, 5, 6 1206 : C= correspondent respectivement aux composantes : 1207 : C= 11, 22, 33, 12, 13, 23 1208 : IF (MFR.EQ.1.OR.MFR.EQ.31) THEN 1209 : IF (IFOUR.NE.1) THEN 1210 : NBROBL=101 1211 : SEGINI,NOMID 1212 : LESOBL( 1)='ETH1' 1213 : LESOBL( 2)='ETH2' 1214 : LESOBL( 3)='ETH3' 1215 : LESOBL( 4)='ETH4' 1216 : LESOBL( 5)='ETH5' 1217 : LESOBL( 6)='ETH6' 1218 : LESOBL( 7)='EEL1' 1219 : LESOBL( 8)='EEL2' 1220 : LESOBL( 9)='EEL3' 1221 : LESOBL(10)='EEL4' 1222 : LESOBL(11)='EEL5' 1223 : LESOBL(12)='EEL6' 1224 : LESOBL(13)='EP01' 1225 : LESOBL(14)='EP02' 1226 : LESOBL(15)='EP03' 1227 : LESOBL(16)='EP04' 1228 : LESOBL(17)='EP05' 1229 : LESOBL(18)='EP06' 1230 : LESOBL(19)='EP11' 1231 : LESOBL(20)='EP12' 1232 : LESOBL(21)='EP13' 1233 : LESOBL(22)='EP14' 1234 : LESOBL(23)='EP15' 1235 : LESOBL(24)='EP16' 1236 : LESOBL(25)='EP21' 1237 : LESOBL(26)='EP22' 1238 : LESOBL(27)='EP23' 1239 : LESOBL(28)='EP24' 1240 : LESOBL(29)='EP25' 1241 : LESOBL(30)='EP26' 1242 : LESOBL(31)='EP31' 1243 : LESOBL(32)='EP32' 1244 : LESOBL(33)='EP33' 1245 : LESOBL(34)='EP34' 1246 : LESOBL(35)='EP35' 1247 : LESOBL(36)='EP36' 1248 : LESOBL(37)='EP41' 1249 : LESOBL(38)='EP42' 1250 : LESOBL(39)='EP43' 1251 : LESOBL(40)='EP44' 1252 : LESOBL(41)='EP45' 1253 : LESOBL(42)='EP46' 1254 : LESOBL(43)='ECR1' 1255 : LESOBL(44)='ECR2' 1256 : LESOBL(45)='ECR3' 1257 : LESOBL(46)='ECR4' 1258 : LESOBL(47)='ECR5' 1259 : LESOBL(48)='ECR6' 1260 : LESOBL(49)='EP0E' 1261 : LESOBL(50)='EP1E' 1262 : LESOBL(51)='EP2E' 1263 : LESOBL(52)='EP3E' 1264 : LESOBL(53)='EP4E' 1265 : LESOBL(54)='EP0Q' 1266 : LESOBL(55)='EP1Q' 1267 : LESOBL(56)='EP2Q' 1268 : LESOBL(57)='EP3Q' 1269 : LESOBL(58)='EP4Q' 1270 : LESOBL(59)='X11 ' 1271 : LESOBL(60)='X12 ' 1272 : LESOBL(61)='X13 ' 1273 : LESOBL(62)='X14 ' 1274 : LESOBL(63)='X15 ' 1275 : LESOBL(64)='X16 ' 1276 : LESOBL(65)='X21 ' 1277 : LESOBL(66)='X22 ' 1278 : LESOBL(67)='X23 ' 1279 : LESOBL(68)='X24 ' 1280 : LESOBL(69)='X25 ' 1281 : LESOBL(70)='X26 ' 1282 : LESOBL(71)='X31 ' 1283 : LESOBL(72)='X32 ' 1284 : LESOBL(73)='X33 ' 1285 : LESOBL(74)='X34 ' 1286 : LESOBL(75)='X35 ' 1287 : LESOBL(76)='X36 ' 1288 : LESOBL(77)='FIT ' 1289 : LESOBL(78)='PSI ' 1290 : LESOBL(79)='DTA ' 1291 : LESOBL(80)='SEA0' 1292 : LESOBL(81)='RRA0' 1293 : LESOBL(82)='VEA0' 1294 : LESOBL(83)='KPLA' 1295 : LESOBL(84)='EPSE' 1296 : LESOBL(85)='VET1' 1297 : LESOBL(86)='VET2' 1298 : LESOBL(87)='VET3' 1299 : LESOBL(88)='VET4' 1300 : LESOBL(89)='VET5' 1301 : LESOBL(90)='VET6' 1302 : LESOBL(91)='VEI1' 1303 : LESOBL(92)='VEI2' 1304 : LESOBL(93)='VEI3' 1305 : LESOBL(94)='VEI4' 1306 : LESOBL(95)='VEI5' 1307 : LESOBL(96)='VEI6' 1308 : LESOBL(97)='VP0E' 1309 : LESOBL(98)='VP1E' 1310 : LESOBL(99)='VP2E' 1311 : LESOBL(100)='VP3E' 1312 : LESOBL(101)='VP4E' 1313 : ENDIF 1314 : ENDIF 1315 : 1316 : C 6 - Modele de GATT_MONERIE 1317 : C ============================ 1318 : C= Modele disponible uniquement en formulation MASSIF 1319 : ELSE IF (MATEPL.EQ.107) THEN 1320 : IF (MFR.EQ.1.OR.MFR.EQ.31) THEN 1321 : IF (IFOUR.NE.1) THEN 1322 : NBROBL=6 1323 : SEGINI,NOMID 1324 : LESOBL(1)='EPSE' 1325 : LESOBL(2)='PORO' 1326 : LESOBL(3)='BU ' 1327 : LESOBL(4)='TETA' 1328 : LESOBL(5)='EPSD' 1329 : LESOBL(6)='EPSG' 1330 : ENDIF 1331 : ENDIF 1332 : 1333 : C 7 - Modele UO2 (= OTTOSEN + GATT_MONERIE) 1334 : C =========================================== 1335 : ELSE IF (MATEPL.EQ.108) THEN 1336 : C ===== 1337 : C 7.1 - Elements MASSIFs 1338 : C ===== 1339 : IF (MFR.EQ.1.OR.MFR.EQ.31) THEN 1340 : C= 7.1.1 - Cas TRIDIMENSIONNEL 1341 : IF (IFOUR.EQ.2) THEN 1342 : NBROBL=30 1343 : SEGINI,NOMID 1344 : LESOBL( 1)='EPSE' 1345 : LESOBL( 2)='W1MA' 1346 : LESOBL( 3)='W2MA' 1347 : LESOBL( 4)='W3MA' 1348 : LESOBL( 5)='W1 ' 1349 : LESOBL( 6)='W2 ' 1350 : LESOBL( 7)='W3 ' 1351 : LESOBL( 8)='VF1X' 1352 : LESOBL( 9)='VF1Y' 1353 : LESOBL(10)='VF1Z' 1354 : LESOBL(11)='VF2X' 1355 : LESOBL(12)='VF2Y' 1356 : LESOBL(13)='VF2Z' 1357 : LESOBL(14)='VF3X' 1358 : LESOBL(15)='VF3Y' 1359 : LESOBL(16)='VF3Z' 1360 : LESOBL(17)='IVL1' 1361 : LESOBL(18)='IVL2' 1362 : LESOBL(19)='IVL3' 1363 : LESOBL(20)='PREC' 1364 : LESOBL(21)='EVXX' 1365 : LESOBL(22)='EVYY' 1366 : LESOBL(23)='EVZZ' 1367 : LESOBL(24)='GVXY' 1368 : LESOBL(25)='GVXZ' 1369 : LESOBL(26)='GVYZ' 1370 : LESOBL(27)='PORO' 1371 : LESOBL(28)='BU ' 1372 : LESOBL(29)='EPSF' 1373 : LESOBL(30)='TETA' 1374 : C= 7.1.2 - Cas AXISYMETRIQUE 1375 : ELSE IF (IFOUR.EQ.0) THEN 1376 : NBROBL=23 1377 : SEGINI,NOMID 1378 : LESOBL( 1)='EPSE' 1379 : LESOBL( 2)='W1MA' 1380 : LESOBL( 3)='W2MA' 1381 : LESOBL( 4)='EPSM' 1382 : LESOBL( 5)='W1 ' 1383 : LESOBL( 6)='W2 ' 1384 : LESOBL( 7)='VEPS' 1385 : LESOBL( 8)='VF1X' 1386 : LESOBL( 9)='VF1Y' 1387 : LESOBL(10)='VF2X' 1388 : LESOBL(11)='VF2Y' 1389 : LESOBL(12)='IVL1' 1390 : LESOBL(13)='IVL2' 1391 : LESOBL(14)='IVL3' 1392 : LESOBL(15)='PREC' 1393 : LESOBL(16)='EVRR' 1394 : LESOBL(17)='EVZZ' 1395 : LESOBL(18)='EVTT' 1396 : LESOBL(19)='GVRZ' 1397 : LESOBL(20)='PORO' 1398 : LESOBL(21)='BU ' 1399 : LESOBL(22)='EPSF' 1400 : LESOBL(23)='TETA' 1401 : C= 7.1.3 - Cas des series de FOURIER 1402 : ELSE IF (IFOUR.EQ.1) THEN 1403 : NBROBL=25 1404 : SEGINI NOMID 1405 : LESOBL( 1)='EPSE' 1406 : LESOBL( 2)='W1MA' 1407 : LESOBL( 3)='W2MA' 1408 : LESOBL( 4)='EPSM' 1409 : LESOBL( 5)='W1 ' 1410 : LESOBL( 6)='W2 ' 1411 : LESOBL( 7)='VEPS' 1412 : LESOBL( 8)='VF1X' 1413 : LESOBL( 9)='VF1Y' 1414 : LESOBL(10)='VF2X' 1415 : LESOBL(11)='VF2Y' 1416 : LESOBL(12)='IVL1' 1417 : LESOBL(13)='IVL2' 1418 : LESOBL(14)='IVL3' 1419 : LESOBL(15)='PREC' 1420 : LESOBL(16)='EVRR' 1421 : LESOBL(17)='EVZZ' 1422 : LESOBL(18)='EVTT' 1423 : LESOBL(19)='GVRZ' 1424 : LESOBL(20)='GVRT' 1425 : LESOBL(21)='GVZT' 1426 : LESOBL(22)='PORO' 1427 : LESOBL(23)='BU ' 1428 : LESOBL(24)='EPSF' 1429 : LESOBL(25)='TETA' 1430 : C= 7.1.4 - Cas DEFOrmations PLANES (GENEralisees ou non) 1431 : ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 1432 : NBROBL=23 1433 : SEGINI NOMID 1434 : LESOBL( 1)='EPSE' 1435 : LESOBL( 2)='W1MA' 1436 : LESOBL( 3)='W2MA' 1437 : LESOBL( 4)='EPSM' 1438 : LESOBL( 5)='W1 ' 1439 : LESOBL( 6)='W2 ' 1440 : LESOBL( 7)='VEPS' 1441 : LESOBL( 8)='VF1X' 1442 : LESOBL( 9)='VF1Y' 1443 : LESOBL(10)='VF2X' 1444 : LESOBL(11)='VF2Y' 1445 : LESOBL(12)='IVL1' 1446 : LESOBL(13)='IVL2' 1447 : LESOBL(14)='IVL3' 1448 : LESOBL(15)='PREC' 1449 : LESOBL(16)='EVXX' 1450 : LESOBL(17)='EVYY' 1451 : LESOBL(18)='EVZZ' 1452 : LESOBL(19)='GVXY' 1453 : LESOBL(20)='PORO' 1454 : LESOBL(21)='BU ' 1455 : LESOBL(22)='EPSF' 1456 : LESOBL(23)='TETA' 1457 : C= 7.1.5 - Cas CONTraintes PLANes 1458 : ELSE IF (IFOUR.EQ.-2) THEN 1459 : NBROBL=20 1460 : SEGINI,NOMID 1461 : LESOBL( 1)='EPSE' 1462 : LESOBL( 2)='W1MA' 1463 : LESOBL( 3)='W2MA' 1464 : LESOBL( 4)='W1 ' 1465 : LESOBL( 5)='W2 ' 1466 : LESOBL( 6)='VF1X' 1467 : LESOBL( 7)='VF1Y' 1468 : LESOBL( 8)='VF2X' 1469 : LESOBL( 9)='VF2Y' 1470 : LESOBL(10)='IVL1' 1471 : LESOBL(11)='IVL2' 1472 : LESOBL(12)='PREC' 1473 : LESOBL(13)='EVXX' 1474 : LESOBL(14)='EVYY' 1475 : LESOBL(15)='EVZZ' 1476 : LESOBL(16)='GVXY' 1477 : LESOBL(17)='PORO' 1478 : LESOBL(18)='BU ' 1479 : LESOBL(19)='EPSF' 1480 : LESOBL(20)='TETA' 1481 : C= 7.1.6 - Modes de calcul UNIDIMENSIONNELs (1D) PLAN 1482 : ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN 1483 : C* NBROBL= 1484 : C* SEGINI,NOMID 1485 : C* LESOBL( 1)='EPSE' 1486 : C* LESOBL( 2)='W1MA' 1487 : C* LESOBL( 3)='W2MA' 1488 : C* LESOBL( 4)='W1 ' 1489 : C* LESOBL( 5)='W2 ' 1490 : C* LESOBL( 6)='VF1X' 1491 : C* LESOBL( 7)='VF1Y' 1492 : C* LESOBL( 8)='VF2X' 1493 : C* LESOBL( 9)='VF2Y' 1494 : C* LESOBL(10)='IVL1' 1495 : C* LESOBL(11)='IVL2' 1496 : C* LESOBL(12)='PREC' 1497 : C* LESOBL(13)='EVXX' 1498 : C* LESOBL(14)='EVYY' 1499 : C* LESOBL(15)='EVZZ' 1500 : C* LESOBL(16)='PORO' 1501 : C* LESOBL(17)='BU ' 1502 : C* LESOBL(18)='EPSF' 1503 : C* LESOBL(19)='TETA' 1504 : C= 7.1.7 - Modes de calcul UNIDIMENSIONNELs (1D) AXIS et SPHE 1505 : ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN 1506 : C* NBROBL= 1507 : C* SEGINI,NOMID 1508 : C* LESOBL( 1)='EPSE' 1509 : C* LESOBL( 2)='W1MA' 1510 : C* LESOBL( 3)='W2MA' 1511 : C* LESOBL( 4)='W1 ' 1512 : C* LESOBL( 5)='W2 ' 1513 : C* LESOBL( 6)='VF1X' 1514 : C* LESOBL( 7)='VF1Y' 1515 : C* LESOBL( 8)='VF2X' 1516 : C* LESOBL( 9)='VF2Y' 1517 : C* LESOBL(10)='IVL1' 1518 : C* LESOBL(11)='IVL2' 1519 : C* LESOBL(12)='PREC' 1520 : C* LESOBL(13)='EVRR' 1521 : C* LESOBL(14)='EVZZ' 1522 : C* LESOBL(15)='EVTT' 1523 : C* LESOBL(16)='PORO' 1524 : C* LESOBL(17)='BU ' 1525 : C* LESOBL(18)='EPSF' 1526 : C* LESOBL(19)='TETA' 1527 : ENDIF 1528 : C ===== 1529 : C 7.2 - Elements COQUEs MINCEs 1530 : C ===== 1531 : ELSE IF (MFR.EQ.3) THEN 1532 : C= 7.2.1 - Mode de calcul TRIDIMENSIONNEL 1533 : IF (IFOUR.EQ.2) THEN 1534 : NBROBL=22 1535 : SEGINI,NOMID 1536 : LESOBL( 1)='EPSE' 1537 : LESOBL( 2)='W1MA' 1538 : LESOBL( 3)='W2MA' 1539 : LESOBL( 4)='W1 ' 1540 : LESOBL( 5)='W2 ' 1541 : LESOBL( 6)='VF1X' 1542 : LESOBL( 7)='VF1Y' 1543 : LESOBL( 8)='VF2X' 1544 : LESOBL( 9)='VF2Y' 1545 : LESOBL(10)='IVL1' 1546 : LESOBL(11)='IVL2' 1547 : LESOBL(12)='PREC' 1548 : LESOBL(13)='EVSS' 1549 : LESOBL(14)='EVTT' 1550 : LESOBL(15)='GVST' 1551 : LESOBL(16)='RVSS' 1552 : LESOBL(17)='RVTT' 1553 : LESOBL(18)='RVST' 1554 : LESOBL(19)='PORO' 1555 : LESOBL(20)='BU ' 1556 : LESOBL(21)='EPSF' 1557 : LESOBL(22)='TETA' 1558 : ENDIF 1559 : C ===== 1560 : C 7.3 - Elements COQUEs en CISAILLEMENT TRANSVERSE$ 1561 : C ===== 1562 : ELSE IF (MFR.EQ.9) THEN 1563 : C= 7.3.1 - Mode de calcul TRIDIMENSIONNEL 1564 : IF (IFOUR.EQ.2) THEN 1565 : NBROBL=24 1566 : SEGINI,NOMID 1567 : LESOBL( 1)='EPSE' 1568 : LESOBL( 2)='W1MA' 1569 : LESOBL( 3)='W2MA' 1570 : LESOBL( 4)='W1 ' 1571 : LESOBL( 5)='W2 ' 1572 : LESOBL( 6)='VF1X' 1573 : LESOBL( 7)='VF1Y' 1574 : LESOBL( 8)='VF2X' 1575 : LESOBL( 9)='VF2Y' 1576 : LESOBL(10)='IVL1' 1577 : LESOBL(11)='IVL2' 1578 : LESOBL(12)='PREC' 1579 : LESOBL(13)='EVSS' 1580 : LESOBL(14)='EVTT' 1581 : LESOBL(15)='GVST' 1582 : LESOBL(16)='RVSS' 1583 : LESOBL(17)='RVTT' 1584 : LESOBL(18)='RVST' 1585 : LESOBL(19)='GVSN' 1586 : LESOBL(20)='GVTN' 1587 : LESOBL(21)='PORO' 1588 : LESOBL(22)='BU ' 1589 : LESOBL(23)='EPSF' 1590 : LESOBL(24)='TETA' 1591 : ENDIF 1592 : ENDIF 1593 : C 8 - Modele visqueux avec double endommagement 1594 : C ============================ 1595 : C= Modele disponible uniquement en formulation MASSIF 1596 : ELSE IF (MATEPL.EQ.130) THEN 1597 : IF (MFR.EQ.1.OR.MFR.EQ.31) THEN 1598 : IF (IFOUR.EQ.0) THEN 1599 : NBROBL=4 1600 : SEGINI NOMID 1601 : LESOBL(1)='R' 1602 : LESOBL(2)='EPSE' 1603 : LESOBL(3)='DD' 1604 : LESOBL(4)='DC' 1605 : IPCOMP=NOMID 1606 : SEGDES NOMID 1607 : ELSE IF (IFOUR.EQ.2) THEN 1608 : NBROBL=4 1609 : SEGINI NOMID 1610 : LESOBL(1)='R' 1611 : LESOBL(2)='EPSE' 1612 : LESOBL(3)='DD' 1613 : LESOBL(4)='DC' 1614 : IPCOMP=NOMID 1615 : SEGDES NOMID 1616 : ENDIF 1617 : ENDIF 1618 : C 9 - Modele visqueux CHAB_SINH_R et CHAB_NOR_R 1619 : C ============================================== 1620 : C= Modele disponible uniquement en formulation MASSIF 1621 : ELSE IF (MATEPL.EQ.136.OR.MATEPL.EQ.138) THEN 1622 : IF (MFR.EQ.1.OR.MFR.EQ.31) THEN 1623 : C--------AXIS 1624 : IF (IFOUR.EQ.0) THEN 1625 : NBROBL=2 1626 : SEGINI NOMID 1627 : LESOBL(1)='R' 1628 : LESOBL(2)='EPSE' 1629 : IPCOMP=NOMID 1630 : SEGDES NOMID 1631 : C--------3D 1632 : ELSE IF (IFOUR.EQ.2) THEN 1633 : NBROBL=2 1634 : SEGINI NOMID 1635 : LESOBL(1)='R' 1636 : LESOBL(2)='EPSE' 1637 : IPCOMP=NOMID 1638 : SEGDES NOMID 1639 : ENDIF 1640 : ENDIF 1641 : C 9 - Modele visqueux CHAB_SINH_X et CHAB_NOR_X 1642 : C ============================================== 1643 : C= Modele disponible uniquement en formulation MASSIF 1644 : ELSE IF (MATEPL.EQ.137.OR.MATEPL.EQ.139) THEN 1645 : IF (MFR.EQ.1.OR.MFR.EQ.31) THEN 1646 : C--------AXIS 1647 : IF (IFOUR.EQ.0) THEN 1648 : NBROBL=10 1649 : SEGINI NOMID 1650 : LESOBL( 1)='A1RR' 1651 : LESOBL( 2)='A1ZZ' 1652 : LESOBL( 3)='A1TT' 1653 : LESOBL( 4)='A1RZ' 1654 : LESOBL( 5)='A2RR' 1655 : LESOBL( 6)='A2ZZ' 1656 : LESOBL( 7)='A2TT' 1657 : LESOBL( 8)='A2RZ' 1658 : LESOBL(9)='r' 1659 : LESOBL(10)='EPSE' 1660 : IPCOMP=NOMID 1661 : SEGDES NOMID 1662 : C--------3D 1663 : ELSE IF (IFOUR.EQ.2) THEN 1664 : C write(6,*) 'on est en 3D',IFOUR 1665 : NBROBL=14 1666 : SEGINI NOMID 1667 : LESOBL( 1)='A1XX' 1668 : LESOBL( 2)='A1YY' 1669 : LESOBL( 3)='A1ZZ' 1670 : LESOBL( 4)='A1XY' 1671 : LESOBL( 5)='A1YZ' 1672 : LESOBL( 6)='A1ZX' 1673 : LESOBL( 7)='A2XX' 1674 : LESOBL( 8)='A2YY' 1675 : LESOBL( 9)='A2ZZ' 1676 : LESOBL(10)='A2XY' 1677 : LESOBL(11)='A2YZ' 1678 : LESOBL(12)='A2ZX' 1679 : LESOBL(13)='r' 1680 : LESOBL(14)='EPSE' 1681 : IPCOMP=NOMID 1682 : SEGDES NOMID 1683 : ENDIF 1684 : ENDIF 1685 : C 10 - Modeles de SYMONDS et COWPER (SYCO1 et SYCO2) 1686 : C =================================== 1687 : ELSE IF (MATEPL.EQ.153.OR.MATEPL.EQ.154) THEN 1688 : C element massif ou xfem 1689 : IF (MFR.EQ.1.OR.MFR.EQ.63) THEN 1690 : C Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou 1691 : C= tridimensionnel 1692 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.2) THEN 1693 : NBROBL=2 1694 : SEGINI,NOMID 1695 : LESOBL( 1)='EPSE' 1696 : LESOBL( 2)='VP' 1697 : ENDIF 1698 : ENDIF 1699 : ENDIF 1700 : 1701 : IPCOMP=NOMID 1702 : IF (IPCOMP.NE.0) THEN 1703 : SEGDES,NOMID 1704 : ELSE 1704 : CALL ERREUR(19) 1705 : ENDIF 1706 : 1707 : RETURN 1708 : END 1709 : 1710 : 1711 : 1712 : 1713 : 1714 : 1715 : 1716 : 1717 :
© Cast3M 2003 - All rights reserved.
Disclaimer