1 : C IDFORC SOURCE CHAT 12/06/07 21:15:21 7389 2 : 3 : C======================================================================= 4 : C= RECHERCHE DES NOMS DE COMPOSANTES DE FORCES = 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 IDFORC(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 : * formulation thermohydrique 33 : IF( MFR.EQ.65) THEN 34 : NBROBL=3 35 : SEGINI NOMID 36 : LESOBL(1)='QG ' 37 : LESOBL(2)='QC ' 38 : LESOBL(3)='Q ' 39 : * 40 : * 0/ MACRO ELEMENT 41 : * 42 : ELSEIF (MFR.EQ.61)THEN 43 : C 44 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN 45 : * 46 : * 0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES 47 : * 48 : NBROBL=5 49 : SEGINI NOMID 50 : LESOBL(1)='FX ' 51 : LESOBL(2)='FY ' 52 : LESOBL(3)='MZ ' 53 : LESOBL(4)='FM ' 54 : LESOBL(5)='MM ' 55 : ENDIF 56 : 57 : C 1 - Elements COQUE (3), COQUE EPAISSE (5), POUTRE (7), COQUE en 58 : C CISAILLEMENT TRANSVERSE (9), TUYAU (13), LINESPRING (15), TUYAU 59 : C FISSURE (17), Barre excentree BAEX (49), LIA2 (51), JOI1(75) 60 : C ==================================================================== 61 : ELSE IF (MFR.EQ. 3.OR.MFR.EQ. 5.OR.MFR.EQ. 7.OR.MFR.EQ. 9.OR. 62 : . MFR.EQ.13.OR.MFR.EQ.15.OR.MFR.EQ.17.OR.MFR.EQ.49.OR. 63 : . MFR.EQ.51.OR.MFR.EQ.75) THEN 64 : C ===== 65 : C 1.1 - Tridimensionnel 66 : C ===== 67 : IF (IFOU.EQ.2) THEN 68 : NBROBL=6 69 : SEGINI,NOMID 70 : LESOBL(1)='FX ' 71 : LESOBL(2)='FY ' 72 : LESOBL(3)='FZ ' 73 : LESOBL(4)='MX ' 74 : LESOBL(5)='MY ' 75 : LESOBL(6)='MZ ' 76 : C ===== 77 : C 1.2 - Fourier 78 : C ===== 79 : ELSE IF (IFOU.EQ.1) THEN 80 : NBROBL=4 81 : SEGINI,NOMID 82 : LESOBL(1)='FR ' 83 : LESOBL(2)='FZ ' 84 : LESOBL(3)='FT ' 85 : LESOBL(4)='MT ' 86 : C ===== 87 : C 1.3 - Axisymetrie 88 : C ===== 89 : ELSE IF (IFOU.EQ.0) THEN 90 : NBROBL=3 91 : SEGINI,NOMID 92 : LESOBL(1)='FR ' 93 : LESOBL(2)='FZ ' 94 : LESOBL(3)='MT ' 95 : C ===== 96 : C 1.4 - Bidimensionnel PLAN (CP/DP) 97 : C ===== 98 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN 99 : NBROBL=3 100 : SEGINI,NOMID 101 : LESOBL(1)='FX ' 102 : LESOBL(2)='FY ' 103 : LESOBL(3)='MZ ' 104 : C ===== 105 : C 1.5 - Bidimensionnel PLAN (CP/DP) 106 : C ===== 107 : ELSE IF (IFOU.EQ.-3) THEN 108 : NBROBL=6 109 : SEGINI,NOMID 110 : LESOBL(1)='FX ' 111 : LESOBL(2)='FY ' 112 : LESOBL(3)='MZ ' 113 : LESOBL(4)='FZ ' 114 : LESOBL(5)='MY ' 115 : LESOBL(6)='MX ' 116 : ENDIF 117 : 118 : C 2 - Elements LIQUIDE 119 : C ====================== 120 : ELSE IF (MFR.EQ.11) THEN 121 : NBROBL=2 122 : SEGINI,NOMID 123 : LESOBL(1)='FP ' 124 : LESOBL(2)='FPI ' 125 : 126 : C 3 - Elements TUYAU ACOUSTIQUE PUR 127 : C =================================== 128 : ELSE IF (MFR.EQ.41) THEN 129 : NBROBL=2 130 : SEGINI,NOMID 131 : LESOBL(1)='FPI ' 132 : LESOBL(2)='FP ' 133 : 134 : C 4 - Element de RACCORD LITU 135 : C ============================= 136 : ELSE IF (MFR.EQ.43) THEN 137 : NBROBL=5 138 : SEGINI,NOMID 139 : LESOBL(1)='FX ' 140 : LESOBL(2)='FY ' 141 : LESOBL(3)='FZ ' 142 : LESOBL(4)='FPI ' 143 : LESOBL(5)='FP ' 144 : 145 : C 5 - Elements HOMOGENEISE 146 : C ========================== 147 : ELSE IF (MFR.EQ.37) THEN 148 : C ===== 149 : C 5.1 - Fourier 150 : C ===== 151 : IF (IFOU.EQ.1) THEN 152 : NBROBL=6 153 : SEGINI,NOMID 154 : LESOBL(1)='FP ' 155 : LESOBL(2)='FPI ' 156 : LESOBL(3)='FR ' 157 : LESOBL(4)='MT ' 158 : LESOBL(5)='FT ' 159 : LESOBL(6)='MR ' 160 : C ===== 161 : C 5.2 - Axisymetrie 162 : C ===== 163 : ELSE IF (IFOU.EQ.0) THEN 164 : NBROBL=4 165 : SEGINI,NOMID 166 : LESOBL(1)='FP ' 167 : LESOBL(2)='FPI ' 168 : LESOBL(3)='FR ' 169 : LESOBL(4)='MT ' 170 : C ===== 171 : C 5.3 - Tridimensionnel 172 : C ===== 173 : ELSE IF (IFOU.EQ.2) THEN 174 : NBROBL=6 175 : SEGINI,NOMID 176 : LESOBL(1)='FP ' 177 : LESOBL(2)='FPI ' 178 : LESOBL(3)='FX ' 179 : LESOBL(4)='MY ' 180 : LESOBL(5)='FY ' 181 : LESOBL(6)='MX ' 182 : C ===== 183 : C 5.4 - Bidimensionnel PLAN (DP/CP/DPGE) 184 : C ===== 185 : ELSE 186 : NBROBL=4 187 : SEGINI,NOMID 188 : LESOBL(1)='FP ' 189 : LESOBL(2)='FPI ' 190 : LESOBL(3)='FX ' 191 : LESOBL(4)='FY ' 192 : ENDIF 193 : 194 : C 6 - Element de SURFACE LIBRE 195 : C ============================== 196 : ELSE IF (MFR.EQ.23) THEN 197 : NBROBL=3 198 : SEGINI,NOMID 199 : LESOBL(1)='FP ' 200 : LESOBL(2)='FPI ' 201 : LESOBL(3)='FZ ' 202 : 203 : C 7 - Element JOINT (35),JOINT CISAILLEMENT (53),JOINT GENERALISE (55) 204 : C ===================================================================== 205 : ELSE IF (MFR.EQ.35.OR.MFR.EQ.53.OR.MFR.EQ.55) THEN 206 : C ===== 207 : C 7.1 - Tridimensionnel 208 : C ===== 209 : IF (IFOU.EQ.2) THEN 210 : NBROBL=3 211 : SEGINI,NOMID 212 : LESOBL(1)='FX ' 213 : LESOBL(2)='FY ' 214 : LESOBL(3)='FZ ' 215 : C ===== 216 : C 7.2 - Bidimensionnel PLAN (CP/DP/DPGE) 217 : C ===== 218 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 219 : NBROBL=2 220 : SEGINI,NOMID 221 : LESOBL(1)='FX ' 222 : LESOBL(2)='FY ' 223 : C ===== 224 : C 7.3 - Axisymetrie 225 : C ===== 226 : ELSE IF (IFOU.EQ.0) THEN 227 : NBROBL=2 228 : SEGINI,NOMID 229 : LESOBL(1)='FR ' 230 : LESOBL(2)='FZ ' 231 : ENDIF 232 : 233 : C 8 - Elements MASSIFS, de MEMBRANE (25), UNIAXIALE (27) 234 : C ======================================================== 235 : ELSE IF (MFR.EQ.1.OR.MFR.EQ.25.OR.MFR.EQ.27.OR.MFR.EQ.31) THEN 236 : C ===== 237 : C 8.1 - Bidimensionnel PLAN (CP/DP) 238 : C ===== 239 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 240 : NBROBL=2 241 : SEGINI,NOMID 242 : LESOBL(1)='FX ' 243 : LESOBL(2)='FY ' 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)='FR ' 251 : LESOBL(2)='FZ ' 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)='FR ' 259 : LESOBL(2)='FZ ' 260 : LESOBL(3)='FT ' 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)='FX ' 268 : LESOBL(2)='FY ' 269 : LESOBL(3)='FZ ' 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)='FX ' 277 : LESOBL(2)='FY ' 278 : LESOBL(3)='FZ ' 279 : LESOBL(4)='MY ' 280 : LESOBL(5)='MX ' 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)='FX ' 289 : ELSE IF (IFOU.EQ.7.OR.IFOU.EQ.8) THEN 290 : NBROBL=2 291 : SEGINI,NOMID 292 : LESOBL(1)='FX ' 293 : LESOBL(2)='FY ' 294 : ELSE IF (IFOU.EQ.9.OR.IFOU.EQ.10) THEN 295 : NBROBL=2 296 : SEGINI,NOMID 297 : LESOBL(1)='FX ' 298 : LESOBL(2)='FZ ' 299 : ELSE IF (IFOU.EQ.11) THEN 300 : NBROBL=3 301 : SEGINI,NOMID 302 : LESOBL(1)='FX ' 303 : LESOBL(2)='FY ' 304 : LESOBL(3)='FZ ' 305 : ELSE IF (IFOU.EQ.12.OR.IFOU.EQ.13.OR.IFOU.EQ.15) THEN 306 : NBROBL=1 307 : SEGINI,NOMID 308 : LESOBL(1)='FR ' 309 : ELSE IF (IFOU.EQ.14) THEN 310 : NBROBL=2 311 : SEGINI,NOMID 312 : LESOBL(1)='FR ' 313 : LESOBL(2)='FZ ' 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)='Q ' 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)='FX ' 334 : LESOBL(2)='FY ' 335 : LESOBL(3)='FP ' 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)='FR ' 343 : LESOBL(2)='FZ ' 344 : LESOBL(3)='FP ' 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)='FR ' 352 : LESOBL(2)='FZ ' 353 : LESOBL(3)='FT ' 354 : LESOBL(4)='FP ' 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)='FX ' 362 : LESOBL(2)='FY ' 363 : LESOBL(3)='FZ ' 364 : LESOBL(4)='FP ' 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)='FX ' 377 : LESOBL(2)='FY ' 378 : LESOBL(3)='FP ' 379 : LESOBL(4)='FPQ ' 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)='FR ' 387 : LESOBL(2)='FZ ' 388 : LESOBL(3)='FP ' 389 : LESOBL(4)='FPQ ' 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)='FR ' 397 : LESOBL(2)='FZ ' 398 : LESOBL(3)='FT ' 399 : LESOBL(4)='FP ' 400 : LESOBL(5)='FPQ ' 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)='FX ' 408 : LESOBL(2)='FY ' 409 : LESOBL(3)='FZ ' 410 : LESOBL(4)='FP ' 411 : LESOBL(5)='FPQ ' 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)='FX ' 424 : LESOBL(2)='FY ' 425 : LESOBL(3)='FP ' 426 : LESOBL(4)='FPQ ' 427 : LESOBL(5)='FTP ' 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)='FR ' 435 : LESOBL(2)='FZ ' 436 : LESOBL(3)='FP ' 437 : LESOBL(4)='FPQ ' 438 : LESOBL(5)='FTP ' 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)='FR ' 446 : LESOBL(2)='FZ ' 447 : LESOBL(3)='FT ' 448 : LESOBL(4)='FP ' 449 : LESOBL(5)='FPQ ' 450 : LESOBL(6)='FTP ' 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)='FX ' 458 : LESOBL(2)='FY ' 459 : LESOBL(3)='FZ ' 460 : LESOBL(4)='FP ' 461 : LESOBL(5)='FPQ ' 462 : LESOBL(6)='FTP ' 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)='FP ' 471 : LESOBL(2)='FPI ' 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)='FX ' 483 : LESOBL(2)='FY ' 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)='FR ' 491 : LESOBL(2)='FZ ' 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)='FR ' 499 : LESOBL(2)='FZ ' 500 : LESOBL(3)='FT ' 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)='FX ' 508 : LESOBL(2)='FY ' 509 : LESOBL(3)='FZ ' 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)='FX ' 522 : LESOBL(2)='FY ' 523 : LESOBL(3)='MZ ' 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)='FR ' 531 : LESOBL(2)='FZ ' 532 : LESOBL(3)='MT ' 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)='FR ' 540 : LESOBL(2)='FZ ' 541 : LESOBL(3)='FT ' 542 : LESOBL(4)='MT ' 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)='FX ' 550 : LESOBL(2)='FY ' 551 : LESOBL(3)='FZ ' 552 : LESOBL(4)='MX ' 553 : LESOBL(5)='MY ' 554 : LESOBL(6)='MZ ' 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)='FX ' 563 : LESOBL(2)='FY ' 564 : LESOBL(3)='FZ ' 565 : LESOBL(4)='MX ' 566 : LESOBL(5)='MY ' 567 : LESOBL(6)='MZ ' 568 : LESOBL(7)='FW0 ' 569 : LESOBL(8)='FU1 ' 570 : LESOBL(9)='FU2 ' 571 : LESOBL(10)='FU3 ' 572 : LESOBL(11)='FU4 ' 573 : LESOBL(12)='FW1 ' 574 : LESOBL(13)='FW2 ' 575 : LESOBL(14)='FW3 ' 576 : LESOBL(15)='FW4 ' 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)='FALF' 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)='FBET' 589 : C 19 - Element XQ4R (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)='FX ' 598 : LESOBL(2)='FY ' 599 : LESFAC(1)='FAX ' 600 : LESFAC(2)='FAY ' 601 : LESFAC(3)='FB1X' 602 : LESFAC(4)='FB1Y' 603 : LESFAC(5)='FC1X' 604 : LESFAC(6)='FC1Y' 605 : LESFAC(7)='FD1X' 606 : LESFAC(8)='FD1Y' 607 : LESFAC(9)='FE1X' 608 : LESFAC(10)='FE1Y' 609 : LESFAC(11)='FB2X' 610 : LESFAC(12)='FB2Y' 611 : LESFAC(13)='FC2X' 612 : LESFAC(14)='FC2Y' 613 : LESFAC(15)='FD2X' 614 : LESFAC(16)='FD2Y' 615 : LESFAC(17)='FE2X' 616 : LESFAC(18)='FE2Y' 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)='FX ' 623 : LESOBL(2)='FY ' 624 : LESOBL(3)='FZ ' 625 : LESFAC(1)='FAX ' 626 : LESFAC(2)='FAY ' 627 : LESFAC(3)='FAZ ' 628 : LESFAC(4)='FB1X' 629 : LESFAC(5)='FB1Y' 630 : LESFAC(6)='FB1Z' 631 : LESFAC(7)='FC1X' 632 : LESFAC(8)='FC1Y' 633 : LESFAC(9)='FC1Z' 634 : LESFAC(10)='FD1X' 635 : LESFAC(11)='FD1Y' 636 : LESFAC(12)='FD1Z' 637 : LESFAC(13)='FE1X' 638 : LESFAC(14)='FE1Y' 639 : LESFAC(15)='FE1Z' 640 : LESFAC(16)='FB2X' 641 : LESFAC(17)='FB2Y' 642 : LESFAC(18)='FB2Z' 643 : LESFAC(19)='FC2X' 644 : LESFAC(20)='FC2Y' 645 : LESFAC(21)='FC2Z' 646 : LESFAC(22)='FD2X' 647 : LESFAC(23)='FD2Y' 648 : LESFAC(24)='FD2Z' 649 : LESFAC(25)='FE2X' 650 : LESFAC(26)='FE2Y' 651 : LESFAC(27)='FE2Z' 652 : ENDIF 653 : 654 : 655 : C 20 - Elements de zones cohesives 656 : C================================= 657 : ELSEIF(MFR.EQ.77) THEN 658 : C - Element ZCO2 (xfem meca rupture en 2D) 659 : IF(IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 660 : NBROBL=2 661 : NBRFAC=18 662 : SEGINI,NOMID 663 : LESOBL(1)='FAX ' 664 : LESOBL(2)='FAY ' 665 : LESFAC(1)='FX ' 666 : LESFAC(2)='FY ' 667 : LESFAC(3)='FB1X' 668 : LESFAC(4)='FB1Y' 669 : LESFAC(5)='FC1X' 670 : LESFAC(6)='FC1Y' 671 : LESFAC(7)='FD1X' 672 : LESFAC(8)='FD1Y' 673 : LESFAC(9)='FE1X' 674 : LESFAC(10)='FE1Y' 675 : LESFAC(11)='FB2X' 676 : LESFAC(12)='FB2Y' 677 : LESFAC(13)='FC2X' 678 : LESFAC(14)='FC2Y' 679 : LESFAC(15)='FD2X' 680 : LESFAC(16)='FD2Y' 681 : LESFAC(17)='FE2X' 682 : LESFAC(18)='FE2Y' 683 : C - Element ZCO3 et ZCO4 (xfem meca rupture en 3D) 684 : ELSE IF (IFOU.EQ.2) THEN 685 : NBROBL=3 686 : NBRFAC=27 687 : SEGINI,NOMID 688 : LESOBL(1)='FAX ' 689 : LESOBL(2)='FAY ' 690 : LESOBL(3)='FAZ ' 691 : LESFAC(1)='FX ' 692 : LESFAC(2)='FY ' 693 : LESFAC(3)='FZ ' 694 : LESFAC(4)='FB1X' 695 : LESFAC(5)='FB1Y' 696 : LESFAC(6)='FB1Z' 697 : LESFAC(7)='FC1X' 698 : LESFAC(8)='FC1Y' 699 : LESFAC(9)='FC1Z' 700 : LESFAC(10)='FD1X' 701 : LESFAC(11)='FD1Y' 702 : LESFAC(12)='FD1Z' 703 : LESFAC(13)='FE1X' 704 : LESFAC(14)='FE1Y' 705 : LESFAC(15)='FE1Z' 706 : LESFAC(16)='FB2X' 707 : LESFAC(17)='FB2Y' 708 : LESFAC(18)='FB2Z' 709 : LESFAC(19)='FC2X' 710 : LESFAC(20)='FC2Y' 711 : LESFAC(21)='FC2Z' 712 : LESFAC(22)='FD2X' 713 : LESFAC(23)='FD2Y' 714 : LESFAC(24)='FD2Z' 715 : LESFAC(25)='FE2X' 716 : LESFAC(26)='FE2Y' 717 : LESFAC(27)='FE2Z' 718 : ENDIF 719 : 720 : C 21 - Formulation ELECTROSTATIQUE (base MASSIF) 721 : C ================================== 722 : ELSE IF (MFR.EQ.71) THEN 723 : NBROBL=1 724 : SEGINI,NOMID 725 : LESOBL(1)='QEL ' 726 : 727 : C 22 - Formulation DIFFUSION (base MASSIF) 728 : C ============================ 729 : ELSE IF (MFR.EQ.73) THEN 730 : C* La grandeur DUALe contenue dans TYMODE(2) est ajoutee ailleurs 731 : NBROBL=1 732 : SEGINI,NOMID 733 : LESOBL(1)='Q.. ' 734 : ELSEIF(MFR.EQ.79) THEN 735 : NBROBL=1 736 : SEGINI,NOMID 737 : LESOBL(1)='Q ' 738 : 739 : ENDIF 740 : 741 : C Par DEFAUT : segment VIDE 742 : C =========================== 743 : IF (NOMID.EQ.0) THEN 744 : SEGINI,NOMID 745 : ELSE 746 : if (ifomod.eq.6) then 747 : nbrfa0 = nbrfac 748 : NBRFAC = NBROBL + (nbrfa0*2) 749 : segadj nomid 750 : do imo = 1,nbrobl 751 : lesfac(nbrfa0 + imo)(2:4) = lesobl(imo)(1:3) 752 : lesfac(nbrfa0 + imo)(1:1) = 'I' 753 : enddo 754 : do imo = 1,nbrfa0 755 : lesfac(nbrfa0+nbrobl+imo)(2:4) = lesfac(imo)(1:3) 756 : lesfac(nbrfa0+nbrobl+imo)(1:1) = 'I' 757 : enddo 758 : endif 759 : ENDIF 760 : 761 : SEGDES,NOMID 762 : IPNOMC=NOMID 763 : 764 : RETURN 765 : END 766 : 767 : 768 : 769 : 770 : 771 :
© Cast3M 2003 - All rights reserved.
Disclaimer