1 : C IDDEIN SOURCE AM 12/02/28 00:28:26 7288 2 : 3 : C======================================================================= 4 : C= RECHERCHE DES NOMS DE COMPOSANTES DE DEFORMATIONS INELASTIQUES = 5 : C= -------------------------------------------------------------- = 6 : C= = 7 : C= Entrees : = 8 : C= IPMODL pointeur ACTIF sur objet IMODEL = 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 IDDEIN (IPMODL,IFOU,IPNOMC,NBROBL,NBRFAC) 19 : 20 : IMPLICIT INTEGER(I-N) 21 : IMPLICIT REAL*8(A-H,O-Z) 22 : 23 : -INC CCOPTIO 24 : -INC SMMODEL 25 : 26 : NBROBL=0 27 : NBRFAC=0 28 : NOMID=0 29 : 30 : IMODEL=IPMODL 31 : MELE=NEFMOD 32 : MFR=NUMMFR(MELE) 33 : NFOR=FORMOD(/2) 34 : 35 : mfr2=0 35 : CALL PLACE(FORMOD,NFOR,ITHEHY,'THERMOHYDRIQUE') 36 : IF(ITHEHY.EQ.1) mfr2=65 36 : CALL PLACE(FORMOD,NFOR,ITTHE,'CONVECTION') 36 : CALL PLACE(FORMOD,NFOR,ITHER,'THERMIQUE') 37 : ITHER=ITHER+ITTHE 38 : if (ITHER.EQ.1) mfr2=29 38 : CALL PLACE(FORMOD,NFOR,IMAGN,'MAGNETODYNAMIQUE') 39 : if (IMAGN.EQ.1) mfr2=69 39 : CALL PLACE(FORMOD,NFOR,IELEC,'ELECTROSTATIQUE') 40 : IF (IELEC.NE.0) mfr2=71 40 : CALL PLACE(FORMOD,NFOR,IDIFF,'DIFFUSION') 41 : IF (IDIFF.NE.0) mfr2=73 42 : IF (mfr2.NE.0) GOTO 100 43 : 44 : C Coque integree ou pas ? 45 : NPINT=0 46 : IF (INFMOD(/1).NE.0) NPINT=INFMOD(1) 47 : * 48 : * -1/ MACRO ELEMENT 49 : * 50 : IF (MFR.EQ.61)THEN 51 : C 52 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN 53 : * 54 : * 0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES 55 : * 56 : NBROBL=5 57 : SEGINI NOMID 58 : LESOBL(1)='EIS1' 59 : LESOBL(2)='CIZ1' 60 : LESOBL(3)='GIXY' 61 : LESOBL(4)='EIS2' 62 : LESOBL(5)='CIZ2' 63 : ENDIF 64 : 65 : C 1 - Element JOINT 66 : C =================== 67 : ELSE IF (MFR.EQ.35) THEN 68 : C ===== 69 : C 1.1 - Tridimensionnel 70 : C ===== 71 : IF (IFOU.EQ.2) THEN 72 : NBROBL=3 73 : SEGINI,NOMID 74 : LESOBL(1)='FRS1' 75 : LESOBL(2)='FRS2' 76 : LESOBL(3)='FRN ' 77 : C ===== 78 : C 1.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie 79 : C ===== 80 : ELSE IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR. 81 : . (IFOU.EQ.0)) THEN 82 : NBROBL=2 83 : SEGINI,NOMID 84 : LESOBL(1)='FRSN' 85 : LESOBL(2)='FRN ' 86 : ENDIF 87 : 88 : C 2 - Element JOINT CISAILLEMENT (2D) 89 : C Bidimensionnel PLAN (CP/DP/DPGE) 90 : C ====================================== 91 : ELSE IF (MFR.EQ.53) THEN 92 : IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1)) THEN 93 : NBROBL=2 94 : SEGINI,NOMID 95 : LESOBL(1)='FRS1' 96 : LESOBL(2)='FRS2' 97 : ENDIF 98 : 99 : C 3 - Element JOINT GENERALISE 100 : C ============================== 101 : ELSE IF (MFR.EQ.55) THEN 102 : C ===== 103 : C 3.1 - Bidimensionnel PLAN (CP/DP/DPGE) 104 : C ===== 105 : IF ((IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.(IFOU.EQ.-3)) THEN 106 : NBROBL=4 107 : SEGINI,NOMID 108 : LESOBL(1)='FRSN' 109 : LESOBL(2)='FRN ' 110 : LESOBL(3)='EITT' 111 : LESOBL(4)='EIZZ' 112 : C ===== 113 : C 3.2 - Tridimensionnel 114 : C ===== 115 : ELSE IF (IFOU.EQ.2) THEN 116 : NBROBL=6 117 : SEGINI,NOMID 118 : LESOBL(1)='FRS1' 119 : LESOBL(2)='FRS2' 120 : LESOBL(3)='FRN ' 121 : LESOBL(4)='EITT' 122 : LESOBL(5)='EIUU' 123 : LESOBL(6)='GITU' 124 : ENDIF 125 : 126 : C 4 - Elements MASSIFS 127 : C ====================== 128 : ELSE IF (MFR.EQ.1.OR.MFR.EQ.31.or.MFR.EQ.63) THEN 129 : C ===== 130 : C 4.1 - Bidimensionnel PLAN (CP/DP/DPGE) 131 : C ===== 132 : IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 133 : NBROBL=4 134 : SEGINI,NOMID 135 : LESOBL(1)='EIXX' 136 : LESOBL(2)='EIYY' 137 : LESOBL(3)='EIZZ' 138 : LESOBL(4)='GIXY' 139 : C ===== 140 : C 4.2 - Axisymetrie 141 : C ===== 142 : ELSE IF (IFOU.EQ.0) THEN 143 : NBROBL=4 144 : SEGINI,NOMID 145 : LESOBL(1)='EIRR' 146 : LESOBL(2)='EIZZ' 147 : LESOBL(3)='EITT' 148 : LESOBL(4)='GIRZ' 149 : C ===== 150 : C 4.3 - Fourier 151 : C ===== 152 : ELSE IF (IFOU.EQ.1) THEN 153 : NBROBL=6 154 : SEGINI,NOMID 155 : LESOBL(1)='EIRR' 156 : LESOBL(2)='EIZZ' 157 : LESOBL(3)='EITT' 158 : LESOBL(4)='GIRZ' 159 : LESOBL(5)='GIRT' 160 : LESOBL(6)='GIZT' 161 : C ===== 162 : C 4.4 - Tridimensionnel 163 : C ===== 164 : ELSE IF (IFOU.EQ.2) THEN 165 : NBROBL=6 166 : SEGINI,NOMID 167 : LESOBL(1)='EIXX' 168 : LESOBL(2)='EIYY' 169 : LESOBL(3)='EIZZ' 170 : LESOBL(4)='GIXY' 171 : LESOBL(5)='GIXZ' 172 : LESOBL(6)='GIYZ' 173 : C ===== 174 : C 4.5 - Unidimensionnel (1D) 175 : C ===== 176 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN 177 : NBROBL=3 178 : SEGINI,NOMID 179 : IF (IFOU.LE.11) THEN 180 : LESOBL(1)='EIXX' 181 : LESOBL(2)='EIYY' 182 : LESOBL(3)='EIZZ' 183 : ELSE IF (IFOU.GE.12.AND.IFOU.LE.14) THEN 184 : LESOBL(1)='EIRR' 185 : LESOBL(2)='EIZZ' 186 : LESOBL(3)='EITT' 187 : ELSE IF (IFOU.EQ.15) THEN 188 : LESOBL(1)='EIRR' 189 : C Voir TC pour LESOBL(2)='EIFF' 190 : LESOBL(2)='EIZZ' 191 : LESOBL(3)='EITT' 192 : ENDIF 193 : ENDIF 194 : 195 : C 5 - Elements COQUES 196 : C ===================== 197 : ELSE IF (MFR.EQ.3.OR.MFR.EQ.39) THEN 198 : IF (NPINT.EQ.0) THEN 199 : C ===== 200 : C 5.1 - Tridimensionnel et Fourier 201 : C ===== 202 : IF (IFOU.EQ.2.OR.IFOU.EQ.1) THEN 203 : NBROBL=6 204 : SEGINI,NOMID 205 : LESOBL(1)='EISS' 206 : LESOBL(2)='EITT' 207 : LESOBL(3)='GIST' 208 : LESOBL(4)='RISS' 209 : LESOBL(5)='RITT' 210 : LESOBL(6)='RIST' 211 : C ===== 212 : C 5.2 - Axisymetrie 213 : C ===== 214 : ELSE IF (IFOU.EQ.0) THEN 215 : NBROBL=4 216 : SEGINI,NOMID 217 : LESOBL(1)='EISS' 218 : LESOBL(2)='EITT' 219 : LESOBL(3)='RISS' 220 : LESOBL(4)='RITT' 221 : C ===== 222 : C 5.3 - Bidimensionnel PLAN (CP/DP/DPGE) 223 : C ===== 224 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 225 : NBROBL=4 226 : SEGINI,NOMID 227 : LESOBL(1)='EISS' 228 : LESOBL(2)='EIZZ' 229 : LESOBL(3)='RISS' 230 : LESOBL(4)='RIZZ' 231 : ENDIF 232 : ELSE 233 : C ===== 234 : C 5.4 - Tridimensionnel - Cas particulier 235 : C ===== 236 : IF (MFR.EQ.3.AND.IFOU.EQ.2) THEN 237 : NBROBL=4 238 : SEGINI,NOMID 239 : LESOBL(1)='EISS' 240 : LESOBL(2)='EITT' 241 : LESOBL(3)='EINN' 242 : LESOBL(4)='GIST' 243 : ENDIF 244 : ENDIF 245 : 246 : C 6 - Elements de COQUE EPAISSE 247 : C Tridimensionnel, axisymetrie et Fourier 248 : C ============================================= 249 : ELSE IF (MFR.EQ.5) THEN 250 : NBROBL=5 251 : SEGINI,NOMID 252 : LESOBL(1)='EISS' 253 : LESOBL(2)='EITT' 254 : LESOBL(3)='GIST' 255 : LESOBL(4)='GISN' 256 : LESOBL(5)='GITN' 257 : 258 : C 7 - LIA2 : Element 3D de liaison a 2 noeuds 259 : C ============================================= 260 : ELSE IF (MFR.EQ.51) THEN 261 : NBROBL=6 262 : SEGINI,NOMID 263 : LESOBL(1)='FEPS' 264 : LESOBL(2)='FGXY' 265 : LESOBL(3)='FGXZ' 266 : LESOBL(4)='FFCX' 267 : LESOBL(5)='FFCY' 268 : LESOBL(6)='FFCZ' 269 : C 7bis - JOI1: Element 3D de liaison a 2 noeuds 270 : C ============================================= 271 : ELSE IF (MFR.EQ.75) THEN 272 : IF (IFOU.EQ.2) THEN 273 : NBROBL=6 274 : SEGINI,NOMID 275 : LESOBL(1)='FEXX ' 276 : LESOBL(2)='FGXY ' 277 : LESOBL(3)='FGXZ ' 278 : LESOBL(4)='FCXX ' 279 : LESOBL(5)='FCXY ' 280 : LESOBL(6)='FCXZ ' 281 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 282 : NBROBL=3 283 : SEGINI,NOMID 284 : LESOBL(1)='FEXX ' 285 : LESOBL(2)='FGXY ' 286 : LESOBL(3)='FCXZ ' 287 : ENDIF 288 : 289 : C 8 - BAEX : BARRE EXCENTREE 290 : C ============================ 291 : ELSE IF (MFR.EQ.49) THEN 292 : NBROBL=1 293 : SEGINI,NOMID 294 : LESOBL(1)='EISS' 295 : 296 : C 9 - Elements POUTRE et TUYAU 297 : C ============================== 298 : ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN 299 : C ===== 300 : C 9.1 - Tridimensionnel 301 : C ===== 302 : IF (IFOU.EQ.2) THEN 303 : NBROBL=6 304 : SEGINI,NOMID 305 : LESOBL(1)='FEPS' 306 : LESOBL(2)='FGXY' 307 : LESOBL(3)='FGXZ' 308 : LESOBL(4)='FFCX' 309 : LESOBL(5)='FFCY' 310 : LESOBL(6)='FFCZ' 311 : C ===== 312 : C 9.2 - Bidimensionnel PLAN (CP/DP/DPGE) 313 : C ===== 314 : ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 315 : NBROBL=3 316 : SEGINI,NOMID 317 : LESOBL(1)='FEPS' 318 : LESOBL(2)='FGXY' 319 : LESOBL(3)='FFCZ' 320 : ENDIF 321 : 322 : C 10 - Element de COQUE en CISAILLEMENT TRANSVERSE 323 : C ================================================== 324 : ELSE IF (MFR.EQ.9) THEN 325 : IF (NPINT.EQ.0) THEN 326 : C ====== 327 : C 10.1 - Tridimensionnel 328 : C ====== 329 : IF (IFOU.EQ.2) THEN 330 : NBROBL=8 331 : SEGINI,NOMID 332 : LESOBL(1)='EISS' 333 : LESOBL(2)='EITT' 334 : LESOBL(3)='GIST' 335 : LESOBL(4)='RISS' 336 : LESOBL(5)='RITT' 337 : LESOBL(6)='RIST' 338 : LESOBL(7)='GISN' 339 : LESOBL(8)='GITN' 340 : ENDIF 341 : ENDIF 342 : 343 : C 11 - Elements LINESPRING (tridimensionnel) 344 : C ============================================ 345 : ELSE IF (MFR.EQ.15) THEN 346 : IF (IFOU.EQ.2) THEN 347 : NBROBL=6 348 : SEGINI,NOMID 349 : LESOBL(1)='EIZZ' 350 : LESOBL(2)='GIXZ' 351 : LESOBL(3)='GIYZ' 352 : LESOBL(4)='RIXX' 353 : LESOBL(5)='RIZZ' 354 : LESOBL(6)='FDJP' 355 : ENDIF 356 : 357 : C 12 - Elements de MEMBRANE 358 : C =========================== 359 : ELSE IF (MFR.EQ.25) THEN 360 : NBROBL=3 361 : SEGINI,NOMID 362 : LESOBL(1)='EISS' 363 : LESOBL(2)='EITT' 364 : LESOBL(3)='GIST' 365 : 366 : C 13 - Formulation UNIAXIALE 367 : C ============================= 368 : ELSE IF (MFR.EQ.27) THEN 369 : NBROBL=1 370 : SEGINI,NOMID 371 : LESOBL(1)='EISS' 372 : 373 : C 14 - Formulation THERMIQUE 374 : C ============================ 375 : ELSE IF (MFR.EQ.29) THEN 376 : C ====== 377 : C 14.1 - Tridimensionnel 378 : C ====== 379 : IF (IFOU.EQ.2) THEN 380 : NBROBL=3 381 : SEGINI,NOMID 382 : LESOBL(1)='FRDX' 383 : LESOBL(2)='FRDY' 384 : LESOBL(3)='FRDZ' 385 : C ====== 386 : C 14.2 - Bidimensionnel PLAN (CP/DP/DPGE) 387 : C ====== 388 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 389 : NBROBL=2 390 : SEGINI,NOMID 391 : LESOBL(1)='FRDX' 392 : LESOBL(2)='FRDY' 393 : C ====== 394 : C 14.3 - Axisymetrie et Fourier 395 : C ====== 396 : ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN 397 : NBROBL=2 398 : SEGINI,NOMID 399 : LESOBL(1)='FRDX' 400 : LESOBL(2)='FRDY' 401 : C*OF LESOBL(1)='FRDR' 402 : C*OF LESOBL(2)='FRDZ' 403 : C ====== 404 : C 14.4 - Unidimensionnel (1D) 405 : C ====== 406 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN 407 : NBROBL=1 408 : SEGINI,NOMID 409 : C*OF IF (IFOU.LE.11) THEN 410 : LESOBL(1)='FRDX' 411 : C*OF ELSE 412 : C*OF LESOBL(1)='FRDR' 413 : C*OF ENDIF 414 : ENDIF 415 : 416 : C 15 - Formulation POREUX 417 : C ========================= 418 : ELSE IF (MFR.EQ.33) THEN 419 : C ===== 420 : C 15.A - Elements POREUX 421 : C ====== 422 : IF (MELE.GE.79.AND.MELE.LE.83) THEN 423 : C ======== 424 : C 15.A.1 - Bidimensionnel PLAN (CP/DP) 425 : C ======== 426 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 427 : NBROBL=5 428 : SEGINI,NOMID 429 : LESOBL(1)='EIXX' 430 : LESOBL(2)='EIYY' 431 : LESOBL(3)='EIZZ' 432 : LESOBL(4)='GIXY' 433 : LESOBL(5)='PIDE' 434 : C ======== 435 : C 15.A.2 - Axisymetrie 436 : C ======== 437 : ELSE IF (IFOU.EQ.0) THEN 438 : NBROBL=5 439 : SEGINI,NOMID 440 : LESOBL(1)='EIRR' 441 : LESOBL(2)='EIZZ' 442 : LESOBL(3)='EITT' 443 : LESOBL(4)='GIRZ' 444 : LESOBL(5)='PIDE' 445 : C ======== 446 : C 15.A.3 - Fourier 447 : C ======== 448 : ELSE IF (IFOU.EQ.1) THEN 449 : NBROBL=7 450 : SEGINI,NOMID 451 : LESOBL(1)='EIRR' 452 : LESOBL(2)='EIZZ' 453 : LESOBL(3)='EITT' 454 : LESOBL(4)='GIRZ' 455 : LESOBL(5)='GIRT' 456 : LESOBL(6)='GIZT' 457 : LESOBL(7)='PIDE' 458 : C ======== 459 : C 15.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 460 : C ======== 461 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 462 : NBROBL=7 463 : SEGINI,NOMID 464 : LESOBL(1)='EIXX' 465 : LESOBL(2)='EIYY' 466 : LESOBL(3)='EIZZ' 467 : LESOBL(4)='GIXY' 468 : LESOBL(5)='GIXZ' 469 : LESOBL(6)='GIYZ' 470 : LESOBL(7)='PIDE' 471 : ENDIF 472 : C ====== 473 : C 15.B - Elements JOINT POREUX 474 : C ====== 475 : ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN 476 : C ======== 477 : C 15.B.1 - Tridimensionnel 478 : C ======== 479 : IF (IFOU.EQ.2) THEN 480 : NBROBL=4 481 : SEGINI,NOMID 482 : LESOBL(1)='FRS1' 483 : LESOBL(2)='FRS2' 484 : LESOBL(3)='FRN ' 485 : LESOBL(4)='PIDE' 486 : C ======== 487 : C 15.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie 488 : C ======== 489 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN 490 : NBROBL=3 491 : SEGINI,NOMID 492 : LESOBL(1)='FRSN' 493 : LESOBL(2)='FRN ' 494 : LESOBL(3)='PIDE' 495 : ENDIF 496 : ENDIF 497 : 498 : C 16 - Formulation POREUX type Q 499 : C ================================ 500 : ELSE IF (MFR.EQ.57) THEN 501 : C ===== 502 : C 16.A - Elements POREUX Q 503 : C ====== 504 : IF (MELE.GE.173.AND.MELE.LE.177) THEN 505 : C ======== 506 : C 16.A.1 - Bidimensionnel PLAN (CP/DP) 507 : C ======== 508 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 509 : NBROBL=6 510 : SEGINI,NOMID 511 : LESOBL(1)='EIXX' 512 : LESOBL(2)='EIYY' 513 : LESOBL(3)='EIZZ' 514 : LESOBL(4)='GIXY' 515 : LESOBL(5)='PIDE' 516 : LESOBL(6)='QIDE' 517 : C ======== 518 : C 16.A.2 - Axisymetrie 519 : C ======== 520 : ELSE IF (IFOU.EQ.0) THEN 521 : NBROBL=6 522 : SEGINI,NOMID 523 : LESOBL(1)='EIRR' 524 : LESOBL(2)='EIZZ' 525 : LESOBL(3)='EITT' 526 : LESOBL(4)='GIRZ' 527 : LESOBL(5)='PIDE' 528 : LESOBL(6)='QIDE' 529 : C ======== 530 : C 16.A.3 - Fourier 531 : C ======== 532 : ELSE IF (IFOU.EQ.1) THEN 533 : NBROBL=8 534 : SEGINI,NOMID 535 : LESOBL(1)='EIRR' 536 : LESOBL(2)='EIZZ' 537 : LESOBL(3)='EITT' 538 : LESOBL(4)='GIRZ' 539 : LESOBL(5)='GIRT' 540 : LESOBL(6)='GIZT' 541 : LESOBL(7)='PIDE' 542 : LESOBL(8)='QIDE' 543 : C ======== 544 : C 16.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 545 : C ======== 546 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 547 : NBROBL=8 548 : SEGINI,NOMID 549 : LESOBL(1)='EIXX' 550 : LESOBL(2)='EIYY' 551 : LESOBL(3)='EIZZ' 552 : LESOBL(4)='GIXY' 553 : LESOBL(5)='GIXZ' 554 : LESOBL(6)='GIYZ' 555 : LESOBL(7)='PIDE' 556 : LESOBL(8)='QIDE' 557 : ENDIF 558 : C ====== 559 : C 16.B - Elements JOINT POREUX Q 560 : C ====== 561 : ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN 562 : C ======== 563 : C 16.B.1 - Tridimensionnel 564 : C ======== 565 : IF (IFOU.EQ.2) THEN 566 : NBROBL=5 567 : SEGINI,NOMID 568 : LESOBL(1)='FRS1' 569 : LESOBL(2)='FRS2' 570 : LESOBL(3)='FRN ' 571 : LESOBL(4)='PIDE' 572 : LESOBL(5)='QIDE' 573 : C ======== 574 : C 16.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie 575 : C ======== 576 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN 577 : NBROBL=4 578 : SEGINI,NOMID 579 : LESOBL(1)='FRSN' 580 : LESOBL(2)='FRN ' 581 : LESOBL(3)='PIDE' 582 : LESOBL(4)='QIDE' 583 : ENDIF 584 : ENDIF 585 : 586 : C 17 - Formulation POREUX type R 587 : C ================================ 588 : ELSE IF (MFR.EQ.59) THEN 589 : C ===== 590 : C 17.A - Elements POREUX R 591 : C ====== 592 : IF (MELE.GE.178.AND.MELE.LE.182) THEN 593 : C ======== 594 : C 17.A.1 - Bidimensionnel PLAN (CP/DP) 595 : C ======== 596 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 597 : NBROBL=7 598 : SEGINI,NOMID 599 : LESOBL(1)='EIXX' 600 : LESOBL(2)='EIYY' 601 : LESOBL(3)='EIZZ' 602 : LESOBL(4)='GIXY' 603 : LESOBL(5)='PIDE' 604 : LESOBL(6)='QIDE' 605 : LESOBL(7)='TIDE' 606 : C ======== 607 : C 17.A.2 - Axisymetrie 608 : C ======== 609 : ELSE IF (IFOU.EQ.0) THEN 610 : NBROBL=7 611 : SEGINI,NOMID 612 : LESOBL(1)='EIRR' 613 : LESOBL(2)='EIZZ' 614 : LESOBL(3)='EITT' 615 : LESOBL(4)='GIRZ' 616 : LESOBL(5)='PIDE' 617 : LESOBL(6)='QIDE' 618 : LESOBL(7)='TIDE' 619 : C ======== 620 : C 17.A.3 - Fourier 621 : C ======== 622 : ELSE IF (IFOU.EQ.1) THEN 623 : NBROBL=9 624 : SEGINI,NOMID 625 : LESOBL(1)='EIRR' 626 : LESOBL(2)='EIZZ' 627 : LESOBL(3)='EITT' 628 : LESOBL(4)='GIRZ' 629 : LESOBL(5)='GIRT' 630 : LESOBL(6)='GIZT' 631 : LESOBL(7)='PIDE' 632 : LESOBL(8)='QIDE' 633 : LESOBL(9)='TIDE' 634 : C ======== 635 : C 17.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 636 : C ======== 637 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 638 : NBROBL=9 639 : SEGINI,NOMID 640 : LESOBL(1)='EIXX' 641 : LESOBL(2)='EIYY' 642 : LESOBL(3)='EIZZ' 643 : LESOBL(4)='GIXY' 644 : LESOBL(5)='GIXZ' 645 : LESOBL(6)='GIYZ' 646 : LESOBL(7)='PIDE' 647 : LESOBL(8)='QIDE' 648 : LESOBL(9)='TIDE' 649 : ENDIF 650 : C ====== 651 : C 17.B - Elements JOINT POREUX R 652 : C ====== 653 : ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN 654 : C ======== 655 : C 17.B.1 - Tridimensionnel 656 : C ======== 657 : IF (IFOU.EQ.2) THEN 658 : NBROBL=6 659 : SEGINI,NOMID 660 : LESOBL(1)='FRS1' 661 : LESOBL(2)='FRS2' 662 : LESOBL(3)='FRN ' 663 : LESOBL(4)='PIDE' 664 : LESOBL(5)='QIDE' 665 : LESOBL(6)='TIDE' 666 : C ======== 667 : C 17.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie 668 : C ======== 669 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN 670 : NBROBL=5 671 : SEGINI,NOMID 672 : LESOBL(1)='FRSN' 673 : LESOBL(2)='FRN ' 674 : LESOBL(3)='PIDE' 675 : LESOBL(4)='QIDE' 676 : LESOBL(5)='TIDE' 677 : ENDIF 678 : ENDIF 679 : 680 : C 18 - Elements TUYAU FISSURE (tridimensionnel) 681 : C =============================================== 682 : ELSE IF (MFR.EQ.17) THEN 683 : IF (IFOU.EQ.2) THEN 684 : NBROBL=8 685 : SEGINI,NOMID 686 : LESOBL(1)='FEPS' 687 : LESOBL(2)='FGXY' 688 : LESOBL(3)='FGXZ' 689 : LESOBL(4)='FFCX' 690 : LESOBL(5)='FFCY' 691 : LESOBL(6)='FFCZ' 692 : LESOBL(7)='FPS7' 693 : LESOBL(8)='FPS8' 694 : ENDIF 695 : 696 : C 19 - Elements de SECTION/FIBRE 697 : C ================================ 698 : ELSE IF (MFR.EQ.47) THEN 699 : C ====== 700 : C 19.1 - Bidimensionnel PLAN (CP/DP/DPGE) 701 : C ====== 702 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 703 : NBROBL=2 704 : SEGINI,NOMID 705 : LESOBL(1)='EIXX' 706 : LESOBL(2)='GIXY' 707 : C*OF LESOBL(3)='GIXZ' 708 : C ====== 709 : C 19.2 - Tridimensionnel 710 : C ====== 711 : ELSE IF (IFOU.EQ.2) THEN 712 : NBROBL=3 713 : SEGINI,NOMID 714 : LESOBL(1)='EIXX' 715 : LESOBL(2)='GIXY' 716 : LESOBL(3)='GIXZ' 717 : ENDIF 718 : 719 : C 20 - Elements de zones cohesives 720 : C ================================== 721 : ELSE IF (MFR.EQ.77) THEN 722 : C ===== 723 : C 20.1 - Tridimensionnel 724 : C ===== 725 : IF (IFOU.GT.0) THEN 726 : NBROBL=3 727 : SEGINI,NOMID 728 : LESOBL(1)='FRS1' 729 : LESOBL(2)='FRS2' 730 : LESOBL(3)='FRN ' 731 : C ===== 732 : C 20.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie 733 : C ===== 734 : ELSE IF (IFOU.LE.0) THEN 735 : NBROBL=2 736 : SEGINI,NOMID 737 : LESOBL(1)='FRSN' 738 : LESOBL(2)='FRN ' 739 : ENDIF 740 : 741 : ENDIF 742 : 743 : C Par DEFAUT : segment VIDE 744 : C =========================== 745 : 100 continue 746 : IF (NOMID.EQ.0) SEGINI,NOMID 747 : 748 : SEGDES,NOMID 749 : IPNOMC=NOMID 750 : 751 : RETURN 752 : END 753 : 754 : 755 : 756 : 757 :
© Cast3M 2003 - All rights reserved.
Disclaimer