1 : C IDDEPL SOURCE AM 12/02/06 21:15:43 7281 2 : 3 : C======================================================================= 4 : C= RECHERCHE DES NOMS DE COMPOSANTES DE DEPLACEMENTS = 5 : C= ------------------------------------------------- = 6 : C= = 7 : C= Entrees : = 8 : C= MFR Numero de la FORMULATION element fini = 9 : C= IFOU valeur de IFOUR de CCOPTIO = 10 : C= = 11 : C= Sorties : = 12 : C= IPNOMC pointeur de type NOMID sur les listes de noms de = 13 : C= composantes OBLigatoires et FACultatives = 14 : C= NBROBL nombre de composantes OBLigatoires = 15 : C= NBRFAC nombre de composantes FACultatives = 16 : C======================================================================= 17 : 18 : SUBROUTINE IDDEPL (MFR,IFOU,IPNOMC,NBROBL,NBRFAC) 19 : 20 : IMPLICIT INTEGER(I-N) 21 : IMPLICIT REAL*8(A-H,O-Z) 22 : -INC CCOPTIO 23 : 24 : SEGMENT NOMID 25 : CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC) 26 : ENDSEGMENT 27 : 28 : NBROBL=0 29 : NBRFAC=0 30 : NOMID=0 31 : 32 : IF (MFR.EQ.65) THEN 33 : NBROBL=3 34 : SEGINI,NOMID 35 : LESOBL(1)='PG ' 36 : LESOBL(2)='PC ' 37 : LESOBL(3)='T ' 38 : * 39 : * 0/ MACRO ELEMENT 40 : * 41 : ELSEIF (MFR.EQ.61)THEN 42 : C 43 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN 44 : * 45 : * 0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES 46 : * 47 : NBROBL=5 48 : SEGINI NOMID 49 : LESOBL(1)='UX ' 50 : LESOBL(2)='UY ' 51 : LESOBL(3)='RZ ' 52 : LESOBL(4)='UM ' 53 : LESOBL(5)='RM ' 54 : ENDIF 55 : 56 : C 1 - Elements COQUE (3), COQUE EPAISSE (5), POUTRE (7), COQUE en 57 : C CISAILLEMENT TRANSVERSE (9), TUYAU (13), LINESPRING (15), TUYAU 58 : C FISSURE (17), Barre excentree BAEX (49), LIA2 (51), JOI1(75) 59 : C ==================================================================== 60 : ELSE IF (MFR.EQ. 3.OR.MFR.EQ. 5.OR.MFR.EQ. 7.OR.MFR.EQ. 9.OR. 61 : . MFR.EQ.13.OR.MFR.EQ.15.OR.MFR.EQ.17.OR.MFR.EQ.49.OR. 62 : . MFR.EQ.51.OR.MFR.EQ.75) THEN 63 : C ===== 64 : C 1.1 - Tridimensionnel 65 : C ===== 66 : IF (IFOU.EQ.2) THEN 67 : NBROBL=6 68 : SEGINI,NOMID 69 : LESOBL(1)='UX ' 70 : LESOBL(2)='UY ' 71 : LESOBL(3)='UZ ' 72 : LESOBL(4)='RX ' 73 : LESOBL(5)='RY ' 74 : LESOBL(6)='RZ ' 75 : C ===== 76 : C 1.2 - Fourier 77 : C ===== 78 : ELSE IF (IFOU.EQ.1) THEN 79 : NBROBL=4 80 : SEGINI,NOMID 81 : LESOBL(1)='UR ' 82 : LESOBL(2)='UZ ' 83 : LESOBL(3)='UT ' 84 : LESOBL(4)='RT ' 85 : C ===== 86 : C 1.3 - Axisymetrie 87 : C ===== 88 : ELSE IF (IFOU.EQ.0) THEN 89 : NBROBL=3 90 : SEGINI,NOMID 91 : LESOBL(1)='UR ' 92 : LESOBL(2)='UZ ' 93 : LESOBL(3)='RT ' 94 : C ===== 95 : C 1.4 - Bidimensionnel PLAN (CP/DP) 96 : C ===== 97 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN 98 : NBROBL=3 99 : SEGINI,NOMID 100 : LESOBL(1)='UX ' 101 : LESOBL(2)='UY ' 102 : LESOBL(3)='RZ ' 103 : C ===== 104 : C 1.5 - Bidimensionnel PLAN (CP/DP) 105 : C ===== 106 : ELSE IF (IFOU.EQ.-3) THEN 107 : NBROBL=6 108 : SEGINI,NOMID 109 : LESOBL(1)='UX ' 110 : LESOBL(2)='UY ' 111 : LESOBL(3)='RZ ' 112 : LESOBL(4)='UZ ' 113 : LESOBL(5)='RY ' 114 : LESOBL(6)='RX ' 115 : ENDIF 116 : 117 : C 2 - Elements LIQUIDE 118 : C ====================== 119 : ELSE IF (MFR.EQ.11) THEN 120 : NBROBL=2 121 : SEGINI,NOMID 122 : LESOBL(1)='P ' 123 : LESOBL(2)='PI ' 124 : 125 : C 3 - Elements TUYAU ACOUSTIQUE PUR 126 : C =================================== 127 : ELSE IF (MFR.EQ.41) THEN 128 : NBROBL=2 129 : SEGINI,NOMID 130 : LESOBL(1)='PI ' 131 : LESOBL(2)='P ' 132 : 133 : C 4 - Element de RACCORD LITU 134 : C ============================= 135 : ELSE IF (MFR.EQ.43) THEN 136 : NBROBL=5 137 : SEGINI,NOMID 138 : LESOBL(1)='UX ' 139 : LESOBL(2)='UY ' 140 : LESOBL(3)='UZ ' 141 : LESOBL(4)='PI ' 142 : LESOBL(5)='P ' 143 : 144 : C 5 - Elements HOMOGENEISE 145 : C ========================== 146 : ELSE IF (MFR.EQ.37) THEN 147 : C ===== 148 : C 5.1 - Fourier 149 : C ===== 150 : IF (IFOU.EQ.1) THEN 151 : NBROBL=6 152 : SEGINI,NOMID 153 : LESOBL(1)='P ' 154 : LESOBL(2)='PI ' 155 : LESOBL(3)='UR ' 156 : LESOBL(4)='RT ' 157 : LESOBL(5)='UT ' 158 : LESOBL(6)='RR ' 159 : C ===== 160 : C 5.2 - Axisymetrie 161 : C ===== 162 : ELSE IF (IFOU.EQ.0) THEN 163 : NBROBL=4 164 : SEGINI,NOMID 165 : LESOBL(1)='P ' 166 : LESOBL(2)='PI ' 167 : LESOBL(3)='UR ' 168 : LESOBL(4)='RT ' 169 : C ===== 170 : C 5.3 - Tridimensionnel 171 : C ===== 172 : ELSE IF (IFOU.EQ.2) THEN 173 : NBROBL=6 174 : SEGINI,NOMID 175 : LESOBL(1)='P ' 176 : LESOBL(2)='PI ' 177 : LESOBL(3)='UX ' 178 : LESOBL(4)='RY ' 179 : LESOBL(5)='UY ' 180 : LESOBL(6)='RX ' 181 : C ===== 182 : C 5.4 - Bidimensionnel PLAN (DP/CP/DPGE) 183 : C ===== 184 : ELSE 185 : NBROBL=4 186 : SEGINI,NOMID 187 : LESOBL(1)='P ' 188 : LESOBL(2)='PI ' 189 : LESOBL(3)='UX ' 190 : LESOBL(4)='UY ' 191 : ENDIF 192 : 193 : C 6 - Element de SURFACE LIBRE 194 : C ============================== 195 : ELSE IF (MFR.EQ.23) THEN 196 : NBROBL=3 197 : SEGINI,NOMID 198 : LESOBL(1)='P ' 199 : LESOBL(2)='PI ' 200 : LESOBL(3)='UZ ' 201 : 202 : C 7 - Element JOINT (35),JOINT CISAILLEMENT (53),JOINT GENERALISE (55) 203 : C ===================================================================== 204 : ELSE IF (MFR.EQ.35.OR.MFR.EQ.53.OR.MFR.EQ.55) THEN 205 : C ===== 206 : C 7.1 - Tridimensionnel 207 : C ===== 208 : IF (IFOU.EQ.2) THEN 209 : NBROBL=3 210 : SEGINI,NOMID 211 : LESOBL(1)='UX ' 212 : LESOBL(2)='UY ' 213 : LESOBL(3)='UZ ' 214 : C ===== 215 : C 7.2 - Bidimensionnel PLAN (CP/DP/DPGE) 216 : C ===== 217 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 218 : NBROBL=2 219 : SEGINI,NOMID 220 : LESOBL(1)='UX ' 221 : LESOBL(2)='UY ' 222 : C ===== 223 : C 7.3 - Axisymetrie 224 : C ===== 225 : ELSE IF (IFOU.EQ.0) THEN 226 : NBROBL=2 227 : SEGINI,NOMID 228 : LESOBL(1)='UR ' 229 : LESOBL(2)='UZ ' 230 : ENDIF 231 : 232 : C 8 - Elements MASSIFS, de MEMBRANE (25), UNIAXIALE (27) 233 : C ======================================================== 234 : ELSE IF (MFR.EQ.1.OR.MFR.EQ.25.OR.MFR.EQ.27.OR.MFR.EQ.31) THEN 235 : C ===== 236 : C 8.1 - Bidimensionnel PLAN (CP/DP) 237 : C ===== 238 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 239 : NBROBL=2 240 : SEGINI,NOMID 241 : LESOBL(1)='UX ' 242 : LESOBL(2)='UY ' 243 : 244 : C ===== 245 : C 8.2 - Axisymetrie 246 : C ===== 247 : ELSE IF (IFOU.EQ.0) THEN 248 : NBROBL=2 249 : SEGINI,NOMID 250 : LESOBL(1)='UR ' 251 : LESOBL(2)='UZ ' 252 : C ===== 253 : C 8.3 - Fourier 254 : C ===== 255 : ELSE IF (IFOU.EQ.1) THEN 256 : NBROBL=3 257 : SEGINI,NOMID 258 : LESOBL(1)='UR ' 259 : LESOBL(2)='UZ ' 260 : LESOBL(3)='UT ' 261 : C ===== 262 : C 8.4 - Tridimensionnel 263 : C ===== 264 : ELSE IF (IFOU.EQ.2) THEN 265 : NBROBL=3 266 : SEGINI,NOMID 267 : LESOBL(1)='UX ' 268 : LESOBL(2)='UY ' 269 : LESOBL(3)='UZ ' 270 : C ===== 271 : C 8.5 - Bidimensionnel PLAN DPGE 272 : C ===== 273 : ELSE IF (IFOU.EQ.-3) THEN 274 : NBROBL=5 275 : SEGINI,NOMID 276 : LESOBL(1)='UX ' 277 : LESOBL(2)='UY ' 278 : LESOBL(3)='UZ ' 279 : LESOBL(4)='RY ' 280 : LESOBL(5)='RX ' 281 : C ===== 282 : C 8.6 - Unidimensionnel (1D) 283 : C ===== 284 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN 285 : IF (IFOU.LE.6) THEN 286 : NBROBL=1 287 : SEGINI,NOMID 288 : LESOBL(1)='UX ' 289 : ELSE IF (IFOU.EQ.7.OR.IFOU.EQ.8) THEN 290 : NBROBL=2 291 : SEGINI,NOMID 292 : LESOBL(1)='UX ' 293 : LESOBL(2)='UY ' 294 : ELSE IF (IFOU.EQ.9.OR.IFOU.EQ.10) THEN 295 : NBROBL=2 296 : SEGINI,NOMID 297 : LESOBL(1)='UX ' 298 : LESOBL(2)='UZ ' 299 : ELSE IF (IFOU.EQ.11) THEN 300 : NBROBL=3 301 : SEGINI,NOMID 302 : LESOBL(1)='UX ' 303 : LESOBL(2)='UY ' 304 : LESOBL(3)='UZ ' 305 : ELSE IF (IFOU.EQ.12.OR.IFOU.EQ.13.OR.IFOU.EQ.15) THEN 306 : NBROBL=1 307 : SEGINI,NOMID 308 : LESOBL(1)='UR ' 309 : ELSE IF (IFOU.EQ.14) THEN 310 : NBROBL=2 311 : SEGINI,NOMID 312 : LESOBL(1)='UR ' 313 : LESOBL(2)='UZ ' 314 : ENDIF 315 : ENDIF 316 : 317 : C 9 - Formulation THERMIQUE 318 : C =========================== 319 : ELSE IF (MFR.EQ.29) THEN 320 : NBROBL=1 321 : SEGINI,NOMID 322 : LESOBL(1)='T ' 323 : 324 : C 10 - Formulation POREUX 325 : C ========================= 326 : ELSE IF (MFR.EQ.33) THEN 327 : C ====== 328 : C 10.1 - Bidimensionnel PLAN (CP/DP) 329 : C ====== 330 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 331 : NBROBL=3 332 : SEGINI,NOMID 333 : LESOBL(1)='UX ' 334 : LESOBL(2)='UY ' 335 : LESOBL(3)='P ' 336 : C ====== 337 : C 10.2 - Axisymetrie 338 : C ====== 339 : ELSE IF (IFOU.EQ.0) THEN 340 : NBROBL=3 341 : SEGINI,NOMID 342 : LESOBL(1)='UR ' 343 : LESOBL(2)='UZ ' 344 : LESOBL(3)='P ' 345 : C ====== 346 : C 10.3 - Fourier 347 : C ====== 348 : ELSE IF (IFOU.EQ.1) THEN 349 : NBROBL=4 350 : SEGINI,NOMID 351 : LESOBL(1)='UR ' 352 : LESOBL(2)='UZ ' 353 : LESOBL(3)='UT ' 354 : LESOBL(4)='P ' 355 : C ====== 356 : C 10.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 357 : C ====== 358 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 359 : NBROBL=4 360 : SEGINI,NOMID 361 : LESOBL(1)='UX ' 362 : LESOBL(2)='UY ' 363 : LESOBL(3)='UZ ' 364 : LESOBL(4)='P ' 365 : ENDIF 366 : 367 : C 11 - Formulation POREUX type Q 368 : C ================================ 369 : ELSE IF (MFR.EQ.57) THEN 370 : C ====== 371 : C 11.1 - Bidimensionnel PLAN (CP/DP) 372 : C ====== 373 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 374 : NBROBL=4 375 : SEGINI,NOMID 376 : LESOBL(1)='UX ' 377 : LESOBL(2)='UY ' 378 : LESOBL(3)='P ' 379 : LESOBL(4)='PQ ' 380 : C ====== 381 : C 11.2 - Axisymetrie 382 : C ====== 383 : ELSE IF (IFOU.EQ.0) THEN 384 : NBROBL=4 385 : SEGINI,NOMID 386 : LESOBL(1)='UR ' 387 : LESOBL(2)='UZ ' 388 : LESOBL(3)='P ' 389 : LESOBL(4)='PQ ' 390 : C ====== 391 : C 11.3 - Fourier 392 : C ====== 393 : ELSE IF (IFOU.EQ.1) THEN 394 : NBROBL=5 395 : SEGINI,NOMID 396 : LESOBL(1)='UR ' 397 : LESOBL(2)='UZ ' 398 : LESOBL(3)='UT ' 399 : LESOBL(4)='P ' 400 : LESOBL(5)='PQ ' 401 : C ====== 402 : C 11.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 403 : C ====== 404 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 405 : NBROBL=5 406 : SEGINI,NOMID 407 : LESOBL(1)='UX ' 408 : LESOBL(2)='UY ' 409 : LESOBL(3)='UZ ' 410 : LESOBL(4)='P ' 411 : LESOBL(5)='PQ ' 412 : ENDIF 413 : 414 : C 12 - Formulation POREUX type R 415 : C ================================ 416 : ELSE IF (MFR.EQ.59) THEN 417 : C ====== 418 : C 12.1 - Bidimensionnel PLAN (CP/DP) 419 : C ====== 420 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 421 : NBROBL=5 422 : SEGINI,NOMID 423 : LESOBL(1)='UX ' 424 : LESOBL(2)='UY ' 425 : LESOBL(3)='P ' 426 : LESOBL(4)='PQ ' 427 : LESOBL(5)='TP ' 428 : C ====== 429 : C 12.2 - Axisymetrie 430 : C ====== 431 : ELSE IF (IFOU.EQ.0) THEN 432 : NBROBL=5 433 : SEGINI,NOMID 434 : LESOBL(1)='UR ' 435 : LESOBL(2)='UZ ' 436 : LESOBL(3)='P ' 437 : LESOBL(4)='PQ ' 438 : LESOBL(5)='TP ' 439 : C ====== 440 : C 12.3 - Fourier 441 : C ====== 442 : ELSE IF (IFOU.EQ.1) THEN 443 : NBROBL=6 444 : SEGINI,NOMID 445 : LESOBL(1)='UR ' 446 : LESOBL(2)='UZ ' 447 : LESOBL(3)='UT ' 448 : LESOBL(4)='P ' 449 : LESOBL(5)='PQ ' 450 : LESOBL(6)='TP ' 451 : C ====== 452 : C 12.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 453 : C ====== 454 : ELSE IF(IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 455 : NBROBL=6 456 : SEGINI,NOMID 457 : LESOBL(1)='UX ' 458 : LESOBL(2)='UY ' 459 : LESOBL(3)='UZ ' 460 : LESOBL(4)='P ' 461 : LESOBL(5)='PQ ' 462 : LESOBL(6)='TP ' 463 : ENDIF 464 : 465 : C 13 - Elements de RACCORD 466 : C ========================== 467 : ELSE IF (MFR.EQ.19.OR.MFR.EQ.21) THEN 468 : NBROBL=2 469 : SEGINI,NOMID 470 : LESOBL(1)='P ' 471 : LESOBL(2)='PI ' 472 : 473 : C 14 - Element de RACCORD MASSIF (2e serie de composantes) 474 : C ========================================================== 475 : ELSE IF (MFR.EQ.1019) THEN 476 : C ====== 477 : C 14.1 - Bidimensionnel PLAN (CP/DP) 478 : C ====== 479 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 480 : NBROBL=2 481 : SEGINI,NOMID 482 : LESOBL(1)='UX ' 483 : LESOBL(2)='UY ' 484 : C ====== 485 : C 14.2 - Axisymetrie 486 : C ====== 487 : ELSE IF (IFOU.EQ.0) THEN 488 : NBROBL=2 489 : SEGINI,NOMID 490 : LESOBL(1)='UR ' 491 : LESOBL(2)='UZ ' 492 : C ====== 493 : C 14.3 - Fourier 494 : C ====== 495 : ELSE IF (IFOU.EQ.1) THEN 496 : NBROBL=3 497 : SEGINI,NOMID 498 : LESOBL(1)='UR ' 499 : LESOBL(2)='UZ ' 500 : LESOBL(3)='UT ' 501 : C ====== 502 : C 14.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 503 : C ====== 504 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 505 : NBROBL=3 506 : SEGINI,NOMID 507 : LESOBL(1)='UX ' 508 : LESOBL(2)='UY ' 509 : LESOBL(3)='UZ ' 510 : ENDIF 511 : 512 : C 15 - Element de RACCORD COQUE (2e serie de composantes) 513 : C ========================================================= 514 : ELSE IF (MFR.EQ.1021) THEN 515 : C ====== 516 : C 15.1 - Bidimensionnel PLAN (CP/DP) 517 : C ====== 518 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 519 : NBROBL=3 520 : SEGINI,NOMID 521 : LESOBL(1)='UX ' 522 : LESOBL(2)='UY ' 523 : LESOBL(3)='RZ ' 524 : C ====== 525 : C 15.2 - Axisymetrie 526 : C ====== 527 : ELSE IF (IFOU.EQ.0) THEN 528 : NBROBL=3 529 : SEGINI,NOMID 530 : LESOBL(1)='UR ' 531 : LESOBL(2)='UZ ' 532 : LESOBL(3)='RT ' 533 : C ====== 534 : C 15.3 - Fourier 535 : C ====== 536 : ELSE IF (IFOU.EQ.1) THEN 537 : NBROBL=4 538 : SEGINI,NOMID 539 : LESOBL(1)='UR ' 540 : LESOBL(2)='UZ ' 541 : LESOBL(3)='UT ' 542 : LESOBL(4)='RT ' 543 : C ====== 544 : C 15.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 545 : C ====== 546 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 547 : NBROBL=6 548 : SEGINI,NOMID 549 : LESOBL(1)='UX ' 550 : LESOBL(2)='UY ' 551 : LESOBL(3)='UZ ' 552 : LESOBL(4)='RX ' 553 : LESOBL(5)='RY ' 554 : LESOBL(6)='RZ ' 555 : ENDIF 556 : 557 : C 16 - Element TUYO 558 : C =================== 559 : ELSE IF (MFR.EQ.39) THEN 560 : NBROBL=6+9 561 : SEGINI,NOMID 562 : LESOBL(1)='UX ' 563 : LESOBL(2)='UY ' 564 : LESOBL(3)='UZ ' 565 : LESOBL(4)='RX ' 566 : LESOBL(5)='RY ' 567 : LESOBL(6)='RZ ' 568 : LESOBL(7)='W0 ' 569 : LESOBL(8)='U1 ' 570 : LESOBL(9)='U2 ' 571 : LESOBL(10)='U3 ' 572 : LESOBL(11)='U4 ' 573 : LESOBL(12)='W1 ' 574 : LESOBL(13)='W2 ' 575 : LESOBL(14)='W3 ' 576 : LESOBL(15)='W4 ' 577 : C 17 - Element POI1 materiau MODAL 578 : C =================== 579 : ELSE IF (MFR.EQ.26) THEN 580 : NBROBL=1 581 : SEGINI,NOMID 582 : LESOBL(1)='ALFA' 583 : C 18 - Element POI1 materiau STATIQUE 584 : C =================== 585 : ELSE IF (MFR.EQ.28) THEN 586 : NBROBL=1 587 : SEGINI,NOMID 588 : LESOBL(1)='BETA' 589 : C 19 - Element XFEM (xfem meca rupture) 590 : C =================== 591 : ELSEIF(MFR.EQ.63) THEN 592 : CTY - Element XQ4R (xfem meca rupture en 2D) 593 : IF(IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 594 : NBROBL=2 595 : NBRFAC=18 596 : SEGINI,NOMID 597 : LESOBL(1)='UX ' 598 : LESOBL(2)='UY ' 599 : LESFAC(1)='AX ' 600 : LESFAC(2)='AY ' 601 : LESFAC(3)='B1X ' 602 : LESFAC(4)='B1Y ' 603 : LESFAC(5)='C1X ' 604 : LESFAC(6)='C1Y ' 605 : LESFAC(7)='D1X ' 606 : LESFAC(8)='D1Y ' 607 : LESFAC(9)='E1X ' 608 : LESFAC(10)='E1Y ' 609 : LESFAC(11)='B2X ' 610 : LESFAC(12)='B2Y ' 611 : LESFAC(13)='C2X ' 612 : LESFAC(14)='C2Y ' 613 : LESFAC(15)='D2X ' 614 : LESFAC(16)='D2Y ' 615 : LESFAC(17)='E2X ' 616 : LESFAC(18)='E2Y ' 617 : CTY - Element XC8R (xfem meca rupture en 3D) 618 : ELSE IF (IFOU.EQ.2) THEN 619 : NBROBL=3 620 : NBRFAC=27 621 : SEGINI,NOMID 622 : LESOBL(1)='UX ' 623 : LESOBL(2)='UY ' 624 : LESOBL(3)='UZ ' 625 : LESFAC(1)='AX ' 626 : LESFAC(2)='AY ' 627 : LESFAC(3)='AZ ' 628 : LESFAC(4)='B1X ' 629 : LESFAC(5)='B1Y ' 630 : LESFAC(6)='B1Z ' 631 : LESFAC(7)='C1X ' 632 : LESFAC(8)='C1Y ' 633 : LESFAC(9)='C1Z ' 634 : LESFAC(10)='D1X ' 635 : LESFAC(11)='D1Y ' 636 : LESFAC(12)='D1Z ' 637 : LESFAC(13)='E1X ' 638 : LESFAC(14)='E1Y ' 639 : LESFAC(15)='E1Z ' 640 : LESFAC(16)='B2X ' 641 : LESFAC(17)='B2Y ' 642 : LESFAC(18)='B2Z ' 643 : LESFAC(19)='C2X ' 644 : LESFAC(20)='C2Y ' 645 : LESFAC(21)='C2Z ' 646 : LESFAC(22)='D2X ' 647 : LESFAC(23)='D2Y ' 648 : LESFAC(24)='D2Z ' 649 : LESFAC(25)='E2X ' 650 : LESFAC(26)='E2Y ' 651 : LESFAC(27)='E2Z ' 652 : ENDIF 653 : 654 : C 20 - Elements de zones cohesives 655 : C ================================== 656 : ELSEIF(MFR.EQ.77) THEN 657 : C - Element ZCO2 (xfem meca rupture en 2D) 658 : IF(IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 659 : NBROBL=2 660 : NBRFAC=18 661 : SEGINI,NOMID 662 : LESOBL(1)='AX ' 663 : LESOBL(2)='AY ' 664 : LESFAC(1)='UX ' 665 : LESFAC(2)='UY ' 666 : LESFAC(3)='B1X ' 667 : LESFAC(4)='B1Y ' 668 : LESFAC(5)='C1X ' 669 : LESFAC(6)='C1Y ' 670 : LESFAC(7)='D1X ' 671 : LESFAC(8)='D1Y ' 672 : LESFAC(9)='E1X ' 673 : LESFAC(10)='E1Y ' 674 : LESFAC(11)='B2X ' 675 : LESFAC(12)='B2Y ' 676 : LESFAC(13)='C2X ' 677 : LESFAC(14)='C2Y ' 678 : LESFAC(15)='D2X ' 679 : LESFAC(16)='D2Y ' 680 : LESFAC(17)='E2X ' 681 : LESFAC(18)='E2Y ' 682 : C - Elements ZCO3 et ZCO4(xfem meca rupture en 3D) 683 : ELSE IF (IFOU.EQ.2) THEN 684 : NBROBL=3 685 : NBRFAC=27 686 : SEGINI,NOMID 687 : LESOBL(1)='AX ' 688 : LESOBL(2)='AY ' 689 : LESOBL(3)='AZ ' 690 : LESFAC(1)='UX ' 691 : LESFAC(2)='UY ' 692 : LESFAC(3)='UZ ' 693 : LESFAC(4)='B1X ' 694 : LESFAC(5)='B1Y ' 695 : LESFAC(6)='B1Z ' 696 : LESFAC(7)='C1X ' 697 : LESFAC(8)='C1Y ' 698 : LESFAC(9)='C1Z ' 699 : LESFAC(10)='D1X ' 700 : LESFAC(11)='D1Y ' 701 : LESFAC(12)='D1Z ' 702 : LESFAC(13)='E1X ' 703 : LESFAC(14)='E1Y ' 704 : LESFAC(15)='E1Z ' 705 : LESFAC(16)='B2X ' 706 : LESFAC(17)='B2Y ' 707 : LESFAC(18)='B2Z ' 708 : LESFAC(19)='C2X ' 709 : LESFAC(20)='C2Y ' 710 : LESFAC(21)='C2Z ' 711 : LESFAC(22)='D2X ' 712 : LESFAC(23)='D2Y ' 713 : LESFAC(24)='D2Z ' 714 : LESFAC(25)='E2X ' 715 : LESFAC(26)='E2Y ' 716 : LESFAC(27)='E2Z ' 717 : ENDIF 718 : 719 : C 21 - Formulation ELECTROSTATIQUE (base MASSIF) 720 : C ================================== 721 : ELSE IF (MFR.EQ.71) THEN 722 : NBROBL=1 723 : SEGINI,NOMID 724 : LESOBL(1)='VEL ' 725 : 726 : C 22 - Formulation DIFFUSION (base MASSIF) 727 : C ============================ 728 : ELSE IF (MFR.EQ.73) THEN 729 : C* L'INCOnnue (DDL) contenue dans TYMODE(1) est ajoutee ailleurs 730 : NBROBL=1 731 : SEGINI,NOMID 732 : LESOBL(1)='.. ' 733 : ENDIF 734 : 735 : C Par DEFAUT : segment VIDE 736 : C =========================== 737 : IF (NOMID.EQ.0) THEN 738 : SEGINI,NOMID 739 : ELSE 740 : if (ifomod.eq.6) then 741 : nbrfa0 = nbrfac 742 : NBRFAC = NBROBL + (nbrfa0*2) 743 : segadj nomid 744 : do imo = 1,nbrobl 745 : lesfac(nbrfa0 + imo)(2:4) = lesobl(imo)(1:3) 746 : lesfac(nbrfa0 + imo)(1:1) = 'I' 747 : enddo 748 : do imo = 1,nbrfa0 749 : lesfac(nbrfa0+nbrobl+imo)(2:4) = lesfac(imo)(1:3) 750 : lesfac(nbrfa0+nbrobl+imo)(1:1) = 'I' 751 : enddo 752 : endif 753 : ENDIF 754 : 755 : SEGDES,NOMID 756 : IPNOMC=NOMID 757 : 758 : RETURN 759 : END 760 : 761 : 762 : 763 : 764 :
© Cast3M 2003 - All rights reserved.
Disclaimer