1 : C IDDEFO SOURCE AM 12/02/06 21:15:40 7281 2 : 3 : C======================================================================= 4 : C= RECHERCHE DES NOMS DE COMPOSANTES DE DEFORMATIONS = 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 IDDEFO (IPMODL,IFOU,IPNOMC,NBROBL,NBRFAC) 19 : 20 : IMPLICIT INTEGER(I-N) 21 : 22 : -INC CCOPTIO 23 : -INC SMMODEL 24 : logical dcmate 25 : 26 : NBROBL=0 27 : NBRFAC=0 28 : NOMID=0 29 : 30 : IMODEL=IPMODL 31 : MELE=NEFMOD 32 : MFR=NUMMFR(MELE) 33 : dcmate = .false. 34 : 35 : IF (FORMOD(1).eq.'MECANIQUE') THEN 36 : DO imat = 1 , matmod(/2) 37 : if (matmod(imat).eq.'IMPEDANCE') then 38 : MFR2 = infele(13) 39 : if (mfr2.eq.26.or.mfr2.eq.28) mfr = mfr2 39 : call iddepl(MFR,IFOU,IPNOMC,NBROBL,NBRFAC) 40 : dcmate = .true. 41 : nomid = ipnomc 42 : segact nomid*mod 43 : do io = 1,nbrobl 44 : lesobl(io)(2:4)=lesobl(io)(1:3) 45 : lesobl(io)(1:1)='E' 46 : enddo 47 : do io = 1,nbrfac 48 : lesfac(io)(2:4)=lesfac(io)(1:3) 49 : lesfac(io)(1:1)='E' 50 : enddo 51 : segdes nomid 52 : return 53 : endif 54 : ENDDO 55 : ENDIF 56 : IF (FORMOD(1).EQ.'ELECTROSTATIQUE') MFR = 71 57 : IF (FORMOD(1).EQ.'DIFFUSION') MFR = 73 58 : 59 : mfr2=0 59 : CALL PLACE(FORMOD,FORMOD(/2),ITHEHY,'THERMOHYDRIQUE') 59 : CALL PLACE(FORMOD,FORMOD(/2),ITTHE,'CONVECTION') 59 : CALL PLACE(FORMOD,FORMOD(/2),ITHER,'THERMIQUE') 60 : ITHER=ITHER+ITTHE 60 : CALL PLACE(FORMOD,FORMOD(/2),IMAGN,'MAGNETODYNAMIQUE') 61 : if(ithehy.eq.1) mfr2=65 62 : if(ITHER.eq.1 ) mfr2=29 63 : if(IMAGN.eq.1 ) mfr2=69 64 : IF(mfr2.ne.0) go to 100 65 : 66 : C Coque integree ou pas ? 67 : NPINT=0 68 : IF (INFMOD(/1).NE.0) NPINT=INFMOD(1) 69 : C 70 : * 71 : * -1/ MACRO ELEMENT 72 : * 73 : IF (MFR.EQ.61)THEN 74 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN 75 : * 76 : * 0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES 77 : * 78 : NBROBL=5 79 : SEGINI NOMID 80 : LESOBL(1)='EPS1' 81 : LESOBL(2)='CZ1 ' 82 : LESOBL(3)='GXY ' 83 : LESOBL(4)='EPS2' 84 : LESOBL(5)='CZ2 ' 85 : ENDIF 86 : 87 : C 1 - Element JOINT 88 : C =================== 89 : ELSE IF (MFR.EQ.35) THEN 90 : C ===== 91 : C 1.1 - Tridimensionnel 92 : C ===== 93 : IF (IFOU.EQ.2) THEN 94 : NBROBL=3 95 : SEGINI,NOMID 96 : LESOBL(1)='DRS1' 97 : LESOBL(2)='DRS2' 98 : LESOBL(3)='DRN ' 99 : C ===== 100 : C 1.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie 101 : C ===== 102 : ELSE IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR. 103 : . (IFOU.EQ.0)) THEN 104 : NBROBL=2 105 : SEGINI,NOMID 106 : LESOBL(1)='DRSN' 107 : LESOBL(2)='DRN ' 108 : ENDIF 109 : 110 : C 2 - Element JOINT CISAILLEMENT (2D) 111 : C Bidimensionnel PLAN (CP/DP/DPGE) 112 : C ====================================== 113 : ELSE IF (MFR.EQ.53) THEN 114 : IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1)) THEN 115 : NBROBL=2 116 : SEGINI,NOMID 117 : LESOBL(1)='DRS1' 118 : LESOBL(2)='DRS2' 119 : ENDIF 120 : 121 : C 3 - Element JOINT GENERALISE 122 : C ============================== 123 : ELSE IF (MFR.EQ.55) THEN 124 : C ===== 125 : C 3.1 - Bidimensionnel PLAN (CP/DP/DPGE) 126 : C ===== 127 : IF ((IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.(IFOU.EQ.-3)) THEN 128 : NBROBL=4 129 : SEGINI,NOMID 130 : LESOBL(1)='DRSN' 131 : LESOBL(2)='DRN ' 132 : LESOBL(3)='EPTT' 133 : LESOBL(4)='EPZZ' 134 : C ===== 135 : C 3.2 - Tridimensionnel 136 : C ===== 137 : ELSE IF (IFOU.EQ.2) THEN 138 : NBROBL=6 139 : SEGINI,NOMID 140 : LESOBL(1)='DRS1' 141 : LESOBL(2)='DRS2' 142 : LESOBL(3)='DRN ' 143 : LESOBL(4)='EPTT' 144 : LESOBL(5)='EPUU' 145 : LESOBL(6)='GATU' 146 : ENDIF 147 : 148 : C 4 - Elements MASSIFS 149 : C ====================== 150 : ELSE IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.63) THEN 151 : C ===== 152 : C 4.1 - Bidimensionnel PLAN (CP/DP/DPGE) 153 : C ===== 154 : IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 155 : NBROBL=4 156 : SEGINI,NOMID 157 : LESOBL(1)='EPXX' 158 : LESOBL(2)='EPYY' 159 : LESOBL(3)='EPZZ' 160 : LESOBL(4)='GAXY' 161 : C ===== 162 : C 4.2 - Axisymetrie 163 : C ===== 164 : ELSE IF (IFOU.EQ.0) THEN 165 : NBROBL=4 166 : SEGINI,NOMID 167 : LESOBL(1)='EPRR' 168 : LESOBL(2)='EPZZ' 169 : LESOBL(3)='EPTT' 170 : LESOBL(4)='GARZ' 171 : C ===== 172 : C 4.3 - Fourier 173 : C ===== 174 : ELSE IF (IFOU.EQ.1) THEN 175 : NBROBL=6 176 : SEGINI,NOMID 177 : LESOBL(1)='EPRR' 178 : LESOBL(2)='EPZZ' 179 : LESOBL(3)='EPTT' 180 : LESOBL(4)='GARZ' 181 : LESOBL(5)='GART' 182 : LESOBL(6)='GAZT' 183 : C ===== 184 : C 4.4 - Tridimensionnel 185 : C ===== 186 : ELSE IF (IFOU.EQ.2) THEN 187 : NBROBL=6 188 : SEGINI,NOMID 189 : LESOBL(1)='EPXX' 190 : LESOBL(2)='EPYY' 191 : LESOBL(3)='EPZZ' 192 : LESOBL(4)='GAXY' 193 : LESOBL(5)='GAXZ' 194 : LESOBL(6)='GAYZ' 195 : C ===== 196 : C 4.5 - Unidimensionnel (1D) 197 : C ===== 198 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN 199 : NBROBL=3 200 : SEGINI,NOMID 201 : IF (IFOU.LE.11) THEN 202 : LESOBL(1)='EPXX' 203 : LESOBL(2)='EPYY' 204 : LESOBL(3)='EPZZ' 205 : ELSE IF (IFOU.GE.12.AND.IFOU.LE.14) THEN 206 : LESOBL(1)='EPRR' 207 : LESOBL(2)='EPZZ' 208 : LESOBL(3)='EPTT' 209 : ELSE IF (IFOU.EQ.15) THEN 210 : LESOBL(1)='EPRR' 211 : C Voir TC pour LESOBL(2)='EPFF' ? 212 : LESOBL(2)='EPZZ' 213 : LESOBL(3)='EPTT' 214 : ENDIF 215 : ENDIF 216 : 217 : C 5 - Elements COQUES 218 : C ===================== 219 : ELSE IF (MFR.EQ.3.OR.MFR.EQ.39) THEN 220 : IF (NPINT.EQ.0) THEN 221 : C ===== 222 : C 5.1 - Tridimensionnel et Fourier 223 : C ===== 224 : IF (IFOU.EQ.2.OR.IFOU.EQ.1) THEN 225 : NBROBL=6 226 : SEGINI,NOMID 227 : LESOBL(1)='EPSS' 228 : LESOBL(2)='EPTT' 229 : LESOBL(3)='GAST' 230 : LESOBL(4)='RTSS' 231 : LESOBL(5)='RTTT' 232 : LESOBL(6)='RTST' 233 : C ===== 234 : C 5.2 - Axisymetrie 235 : C ===== 236 : ELSE IF (IFOU.EQ.0) THEN 237 : NBROBL=4 238 : SEGINI,NOMID 239 : LESOBL(1)='EPSS' 240 : LESOBL(2)='EPTT' 241 : LESOBL(3)='RTSS' 242 : LESOBL(4)='RTTT' 243 : C ===== 244 : C 5.3 - Bidimensionnel PLAN (CP/DP/DPGE) 245 : C ===== 246 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 247 : NBROBL=4 248 : SEGINI,NOMID 249 : LESOBL(1)='EPSS' 250 : LESOBL(2)='EPZZ' 251 : LESOBL(3)='RTSS' 252 : LESOBL(4)='RTZZ' 253 : ENDIF 254 : ELSE 255 : C ===== 256 : C 5.4 - Tridimensionnel - Cas particulier 257 : C ===== 258 : IF (MFR.EQ.3.AND.IFOU.EQ.2) THEN 259 : NBROBL=4 260 : SEGINI,NOMID 261 : LESOBL(1)='EPSS' 262 : LESOBL(2)='EPTT' 263 : LESOBL(3)='EPNN' 264 : LESOBL(4)='GAST' 265 : ENDIF 266 : ENDIF 267 : 268 : C 6 - Elements de COQUE EPAISSE 269 : C Tridimensionnel, axisymetrie et Fourier 270 : C ============================================= 271 : ELSE IF (MFR.EQ.5) THEN 272 : NBROBL=5 273 : SEGINI,NOMID 274 : LESOBL(1)='EPSS' 275 : LESOBL(2)='EPTT' 276 : LESOBL(3)='GAST' 277 : LESOBL(4)='GASN' 278 : LESOBL(5)='GATN' 279 : 280 : C 7 - LIA2: Element 3D de liaison a 2 noeuds 281 : C ============================================= 282 : ELSE IF (MFR.EQ.51) THEN 283 : NBROBL=6 284 : SEGINI,NOMID 285 : LESOBL(1)='EPS ' 286 : LESOBL(2)='GXY ' 287 : LESOBL(3)='GXZ ' 288 : LESOBL(4)='CX ' 289 : LESOBL(5)='CY ' 290 : LESOBL(6)='CZ ' 291 : C 7bis - JOI1: Element 3D de liaison a 2 noeuds 292 : C ============================================= 293 : ELSE IF (MFR.EQ.75) THEN 294 : IF (IFOU.EQ.2) THEN 295 : NBROBL=6 296 : SEGINI,NOMID 297 : LESOBL(1)='EXX ' 298 : LESOBL(2)='GXY ' 299 : LESOBL(3)='GXZ ' 300 : LESOBL(4)='CXX ' 301 : LESOBL(5)='CXY ' 302 : LESOBL(6)='CXZ ' 303 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 304 : NBROBL=3 305 : SEGINI,NOMID 306 : LESOBL(1)='EXX ' 307 : LESOBL(2)='GXY ' 308 : LESOBL(3)='CXZ ' 309 : ENDIF 310 : 311 : C 8 - BAEX : BARRE EXCENTREE 312 : C ============================ 313 : ELSE IF (MFR.EQ.49) THEN 314 : NBROBL=1 315 : SEGINI,NOMID 316 : LESOBL(1)='EPSS' 317 : 318 : C 9 - Elements POUTRE et TUYAU 319 : C ============================== 320 : ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN 321 : C ===== 322 : C 9.1 - Tridimensionnel 323 : C ===== 324 : IF (IFOU.EQ.2) THEN 325 : NBROBL=6 326 : SEGINI,NOMID 327 : LESOBL(1)='EPS ' 328 : LESOBL(2)='GXY ' 329 : LESOBL(3)='GXZ ' 330 : LESOBL(4)='CX ' 331 : LESOBL(5)='CY ' 332 : LESOBL(6)='CZ ' 333 : C ===== 334 : C 9.2 - Bidimensionnel PLAN (CP/DP/DPGE) 335 : C ===== 336 : ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 337 : NBROBL=3 338 : SEGINI,NOMID 339 : LESOBL(1)='EPS ' 340 : LESOBL(2)='GXY ' 341 : LESOBL(3)='CZ ' 342 : ENDIF 343 : 344 : C 10 - Element de COQUE en CISAILLEMENT TRANSVERSE 345 : C ================================================== 346 : ELSE IF (MFR.EQ.9) THEN 347 : IF (NPINT.EQ.0) THEN 348 : C ====== 349 : C 10.1 - Tridimensionnel 350 : C ====== 351 : IF (IFOU.EQ.2) THEN 352 : NBROBL=8 353 : SEGINI,NOMID 354 : LESOBL(1)='EPSS' 355 : LESOBL(2)='EPTT' 356 : LESOBL(3)='GAST' 357 : LESOBL(4)='RTSS' 358 : LESOBL(5)='RTTT' 359 : LESOBL(6)='RTST' 360 : LESOBL(7)='GASN' 361 : LESOBL(8)='GATN' 362 : ENDIF 363 : ENDIF 364 : 365 : C 11 - Elements LINESPRING (tridimensionnel) 366 : C ============================================ 367 : ELSE IF (MFR.EQ.15) THEN 368 : IF (IFOU.EQ.2) THEN 369 : NBROBL=6 370 : SEGINI,NOMID 371 : LESOBL(1)='EPZZ' 372 : LESOBL(2)='GAXZ' 373 : LESOBL(3)='GAYZ' 374 : LESOBL(4)='RTXX' 375 : LESOBL(5)='RTZZ' 376 : LESOBL(6)='DJP ' 377 : ENDIF 378 : 379 : C 12 - Elements de MEMBRANE 380 : C =========================== 381 : ELSE IF (MFR.EQ.25) THEN 382 : NBROBL=3 383 : SEGINI,NOMID 384 : LESOBL(1)='EPSS' 385 : LESOBL(2)='EPTT' 386 : LESOBL(3)='GAST' 387 : 388 : C 13 - Formulation UNIAXIALE 389 : C ============================= 390 : ELSE IF (MFR.EQ.27) THEN 391 : NBROBL=1 392 : SEGINI,NOMID 393 : LESOBL(1)='EPSS' 394 : 395 : C 14 - Formulation THERMIQUE 396 : C ============================ 397 : ELSE IF (MFR.EQ.29) THEN 398 : C ====== 399 : C 14.1 - Tridimensionnel 400 : C ====== 401 : IF (IFOU.EQ.2) THEN 402 : NBROBL=3 403 : SEGINI,NOMID 404 : LESOBL(1)='GRDX' 405 : LESOBL(2)='GRDY' 406 : LESOBL(3)='GRDZ' 407 : C ====== 408 : C 14.2 - Bidimensionnel PLAN (CP/DP/DPGE) 409 : C ====== 410 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 411 : NBROBL=2 412 : SEGINI,NOMID 413 : LESOBL(1)='GRDX' 414 : LESOBL(2)='GRDY' 415 : C ====== 416 : C 14.3 - Axisymetrie et Fourier 417 : C ====== 418 : ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN 419 : NBROBL=2 420 : SEGINI,NOMID 421 : LESOBL(1)='GRDX' 422 : LESOBL(2)='GRDY' 423 : C*OF LESOBL(1)='GRDR' 424 : C*OF LESOBL(2)='GRDZ' 425 : C ====== 426 : C 14.4 - Unidimensionnel (1D) 427 : C ====== 428 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN 429 : NBROBL=1 430 : SEGINI,NOMID 431 : C*OF IF (IFOU.LE.11) THEN 432 : LESOBL(1)='GRDX' 433 : C*OF ELSE 434 : C*OF LESOBL(1)='GRDR' 435 : C*OF ENDIF 436 : ENDIF 437 : 438 : C 15 - Formulation POREUX 439 : C ========================= 440 : ELSE IF (MFR.EQ.33) THEN 441 : C ===== 442 : C 15.A - Elements POREUX 443 : C ====== 444 : IF (MELE.GE.79.AND.MELE.LE.83) THEN 445 : C ======== 446 : C 15.A.1 - Bidimensionnel PLAN (CP/DP) 447 : C ======== 448 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 449 : NBROBL=5 450 : SEGINI,NOMID 451 : LESOBL(1)='EPXX' 452 : LESOBL(2)='EPYY' 453 : LESOBL(3)='EPZZ' 454 : LESOBL(4)='GAXY' 455 : LESOBL(5)='PRES' 456 : C ======== 457 : C 15.A.2 - Axisymetrie 458 : C ======== 459 : ELSE IF (IFOU.EQ.0) THEN 460 : NBROBL=5 461 : SEGINI,NOMID 462 : LESOBL(1)='EPRR' 463 : LESOBL(2)='EPZZ' 464 : LESOBL(3)='EPTT' 465 : LESOBL(4)='GARZ' 466 : LESOBL(5)='PRES' 467 : C ======== 468 : C 15.A.3 - Fourier 469 : C ======== 470 : ELSE IF (IFOU.EQ.1) THEN 471 : NBROBL=7 472 : SEGINI,NOMID 473 : LESOBL(1)='EPRR' 474 : LESOBL(2)='EPZZ' 475 : LESOBL(3)='EPTT' 476 : LESOBL(4)='GARZ' 477 : LESOBL(5)='GART' 478 : LESOBL(6)='GAZT' 479 : LESOBL(7)='PRES' 480 : C ======== 481 : C 15.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 482 : C ======== 483 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 484 : NBROBL=7 485 : SEGINI,NOMID 486 : LESOBL(1)='EPXX' 487 : LESOBL(2)='EPYY' 488 : LESOBL(3)='EPZZ' 489 : LESOBL(4)='GAXY' 490 : LESOBL(5)='GAXZ' 491 : LESOBL(6)='GAYZ' 492 : LESOBL(7)='PRES' 493 : ENDIF 494 : C ====== 495 : C 15.B - Elements JOINT POREUX 496 : C ====== 497 : ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN 498 : C ======== 499 : C 15.B.1 - Tridimensionnel 500 : C ======== 501 : IF (IFOU.EQ.2) THEN 502 : NBROBL=4 503 : SEGINI,NOMID 504 : LESOBL(1)='DRS1' 505 : LESOBL(2)='DRS2' 506 : LESOBL(3)='DRN ' 507 : LESOBL(4)='PRES' 508 : C ======== 509 : C 15.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie 510 : C ======== 511 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN 512 : NBROBL=3 513 : SEGINI,NOMID 514 : LESOBL(1)='DRSN' 515 : LESOBL(2)='DRN ' 516 : LESOBL(3)='PRES' 517 : ENDIF 518 : ENDIF 519 : 520 : C 16 - Formulation POREUX type Q 521 : C ================================ 522 : ELSE IF (MFR.EQ.57) THEN 523 : C ===== 524 : C 16.A - Elements POREUX Q 525 : C ====== 526 : IF (MELE.GE.173.AND.MELE.LE.177) THEN 527 : C ======== 528 : C 16.A.1 - Bidimensionnel PLAN (CP/DP) 529 : C ======== 530 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 531 : NBROBL=6 532 : SEGINI,NOMID 533 : LESOBL(1)='EPXX' 534 : LESOBL(2)='EPYY' 535 : LESOBL(3)='EPZZ' 536 : LESOBL(4)='GAXY' 537 : LESOBL(5)='PRES' 538 : LESOBL(6)='PREQ' 539 : C ======== 540 : C 16.A.2 - Axisymetrie 541 : C ======== 542 : ELSE IF (IFOU.EQ.0) THEN 543 : NBROBL=6 544 : SEGINI,NOMID 545 : LESOBL(1)='EPRR' 546 : LESOBL(2)='EPZZ' 547 : LESOBL(3)='EPTT' 548 : LESOBL(4)='GARZ' 549 : LESOBL(5)='PRES' 550 : LESOBL(6)='PREQ' 551 : C ======== 552 : C 16.A.3 - Fourier 553 : C ======== 554 : ELSE IF (IFOU.EQ.1) THEN 555 : NBROBL=8 556 : SEGINI,NOMID 557 : LESOBL(1)='EPRR' 558 : LESOBL(2)='EPZZ' 559 : LESOBL(3)='EPTT' 560 : LESOBL(4)='GARZ' 561 : LESOBL(5)='GART' 562 : LESOBL(6)='GAZT' 563 : LESOBL(7)='PRES' 564 : LESOBL(8)='PREQ' 565 : C ======== 566 : C 16.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 567 : C ======== 568 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 569 : NBROBL=8 570 : SEGINI,NOMID 571 : LESOBL(1)='EPXX' 572 : LESOBL(2)='EPYY' 573 : LESOBL(3)='EPZZ' 574 : LESOBL(4)='GAXY' 575 : LESOBL(5)='GAXZ' 576 : LESOBL(6)='GAYZ' 577 : LESOBL(7)='PRES' 578 : LESOBL(8)='PREQ' 579 : ENDIF 580 : C ====== 581 : C 16.B - Elements JOINT POREUX Q 582 : C ====== 583 : ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN 584 : C ======== 585 : C 16.B.1 - Tridimensionnel 586 : C ======== 587 : IF (IFOU.EQ.2) THEN 588 : NBROBL=5 589 : SEGINI,NOMID 590 : LESOBL(1)='DRS1' 591 : LESOBL(2)='DRS2' 592 : LESOBL(3)='DRN ' 593 : LESOBL(4)='PRES' 594 : LESOBL(5)='PREQ' 595 : C ======== 596 : C 16.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie 597 : C ======== 598 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN 599 : NBROBL=4 600 : SEGINI,NOMID 601 : LESOBL(1)='DRSN' 602 : LESOBL(2)='DRN ' 603 : LESOBL(3)='PRES' 604 : LESOBL(4)='PREQ' 605 : ENDIF 606 : ENDIF 607 : 608 : C 17 - Formulation POREUX type R 609 : C ================================ 610 : ELSE IF (MFR.EQ.59) THEN 611 : C ===== 612 : C 17.A - Elements POREUX R 613 : C ====== 614 : IF (MELE.GE.178.AND.MELE.LE.182) THEN 615 : C ======== 616 : C 17.A.1 - Bidimensionnel PLAN (CP/DP) 617 : C ======== 618 : IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 619 : NBROBL=7 620 : SEGINI NOMID 621 : LESOBL(1)='EPXX' 622 : LESOBL(2)='EPYY' 623 : LESOBL(3)='EPZZ' 624 : LESOBL(4)='GAXY' 625 : LESOBL(5)='PRES' 626 : LESOBL(6)='PREQ' 627 : LESOBL(7)='TPT ' 628 : C ======== 629 : C 17.A.2 - Axisymetrie 630 : C ======== 631 : ELSE IF (IFOU.EQ.0) THEN 632 : NBROBL=7 633 : SEGINI NOMID 634 : LESOBL(1)='EPRR' 635 : LESOBL(2)='EPZZ' 636 : LESOBL(3)='EPTT' 637 : LESOBL(4)='GARZ' 638 : LESOBL(5)='PRES' 639 : LESOBL(6)='PREQ' 640 : LESOBL(7)='TPT ' 641 : C ======== 642 : C 17.A.3 - Fourier 643 : C ======== 644 : ELSE IF (IFOU.EQ.1) THEN 645 : NBROBL=9 646 : SEGINI,NOMID 647 : LESOBL(1)='EPRR' 648 : LESOBL(2)='EPZZ' 649 : LESOBL(3)='EPTT' 650 : LESOBL(4)='GARZ' 651 : LESOBL(5)='GART' 652 : LESOBL(6)='GAZT' 653 : LESOBL(7)='PRES' 654 : LESOBL(8)='PREQ' 655 : LESOBL(9)='TPT ' 656 : C ======== 657 : C 17.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE 658 : C ======== 659 : ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN 660 : NBROBL=9 661 : SEGINI,NOMID 662 : LESOBL(1)='EPXX' 663 : LESOBL(2)='EPYY' 664 : LESOBL(3)='EPZZ' 665 : LESOBL(4)='GAXY' 666 : LESOBL(5)='GAXZ' 667 : LESOBL(6)='GAYZ' 668 : LESOBL(7)='PRES' 669 : LESOBL(8)='PREQ' 670 : LESOBL(9)='TPT ' 671 : ENDIF 672 : C ====== 673 : C 17.B - Elements JOINT POREUX R 674 : C ====== 675 : ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN 676 : C ======== 677 : C 17.B.1 - Tridimensionnel 678 : C ======== 679 : IF (IFOU.EQ.2) THEN 680 : NBROBL=6 681 : SEGINI,NOMID 682 : LESOBL(1)='DRS1' 683 : LESOBL(2)='DRS2' 684 : LESOBL(3)='DRN ' 685 : LESOBL(4)='PRES' 686 : LESOBL(5)='PREQ' 687 : LESOBL(6)='TPT ' 688 : C ======== 689 : C 17.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie 690 : C ======== 691 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN 692 : NBROBL=5 693 : SEGINI,NOMID 694 : LESOBL(1)='DRSN' 695 : LESOBL(2)='DRN ' 696 : LESOBL(3)='PRES' 697 : LESOBL(4)='PREQ' 698 : LESOBL(5)='TPT ' 699 : ENDIF 700 : ENDIF 701 : 702 : C 18 - Elements TUYAU FISSURE (tridimensionnel) 703 : C =============================================== 704 : ELSE IF (MFR.EQ.17) THEN 705 : IF (IFOU.EQ.2) THEN 706 : NBROBL=8 707 : SEGINI,NOMID 708 : LESOBL(1)='EPS ' 709 : LESOBL(2)='GXY ' 710 : LESOBL(3)='GXZ ' 711 : LESOBL(4)='CX ' 712 : LESOBL(5)='CY ' 713 : LESOBL(6)='CZ ' 714 : LESOBL(7)='EPS7' 715 : LESOBL(8)='EPS8' 716 : ENDIF 717 : 718 : C 19 - Elements de SECTION/FIBRE 719 : C ================================ 720 : ELSE IF (MFR.EQ.47) THEN 721 : C ====== 722 : C 19.1 - Bidimensionnel PLAN (CP/DP/DPGE) 723 : C ====== 724 : IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 725 : NBROBL=2 726 : SEGINI,NOMID 727 : LESOBL(1)='EPXX' 728 : LESOBL(2)='GAXY' 729 : C ====== 730 : C 19.2 - Tridimensionnel 731 : C ====== 732 : ELSE IF (IFOU.EQ.2) THEN 733 : NBROBL=3 734 : SEGINI,NOMID 735 : LESOBL(1)='EPXX' 736 : LESOBL(2)='GAXY' 737 : LESOBL(3)='GAXZ' 738 : ENDIF 739 : 740 : C 20 - Elements de zones cohesives 741 : C ================================== 742 : ELSE IF (MFR.EQ.77) THEN 743 : C ===== 744 : C 20.1 - Tridimensionnel 745 : C ===== 746 : IF (IFOU.GT.0) THEN 747 : NBROBL=3 748 : SEGINI,NOMID 749 : LESOBL(1)='DRS1' 750 : LESOBL(2)='DRS2' 751 : LESOBL(3)='DRN ' 752 : C ===== 753 : C 20.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie 754 : C ===== 755 : ELSE IF (IFOU.LE.0) THEN 756 : NBROBL=2 757 : SEGINI,NOMID 758 : LESOBL(1)='DRSN' 759 : LESOBL(2)='DRN ' 760 : ENDIF 761 : 762 : C 21 - Formulation ELECTROSTATIQUE (base MASSIF) 763 : C ================================== 764 : ELSE IF (MFR.EQ.71) THEN 765 : C ====== 766 : C 21.1 - Tridimensionnel 767 : C ====== 768 : IF (IFOU.EQ.2) THEN 769 : NBROBL=3 770 : SEGINI,NOMID 771 : LESOBL(1)='EELX ' 772 : LESOBL(2)='EELY ' 773 : LESOBL(3)='EELZ ' 774 : C ====== 775 : C 21.2 - Bidimensionnel PLAN (CP/DP/DPGE) 776 : C ====== 777 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 778 : NBROBL=2 779 : SEGINI,NOMID 780 : LESOBL(1)='EELX ' 781 : LESOBL(2)='EELY ' 782 : C ====== 783 : C 21.3 - Axisymetrie et Fourier 784 : C ====== 785 : ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN 786 : NBROBL=2 787 : SEGINI,NOMID 788 : LESOBL(1)='EELR ' 789 : LESOBL(2)='EELZ ' 790 : C ====== 791 : C 21.4 - Unidimensionnel (1D) 792 : C ====== 793 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN 794 : NBROBL=1 795 : SEGINI,NOMID 796 : IF (IFOU.LE.11) THEN 797 : LESOBL(1)='EELX ' 798 : ELSE 799 : LESOBL(1)='EELR ' 800 : ENDIF 801 : ENDIF 802 : 803 : C 22 - Formulation DIFFUSION (base MASSIF) 804 : C ============================ 805 : ELSE IF (MFR.EQ.73) THEN 806 : C ====== 807 : C 22.1 - Tridimensionnel 808 : C ====== 809 : IF (IFOU.EQ.2) THEN 810 : NBROBL=3 811 : SEGINI,NOMID 812 : LESOBL(1)='GX ' 813 : LESOBL(2)='GY ' 814 : LESOBL(3)='GZ ' 815 : C ====== 816 : C 22.2 - Bidimensionnel PLAN (CP/DP/DPGE) 817 : C ====== 818 : ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN 819 : NBROBL=2 820 : SEGINI,NOMID 821 : LESOBL(1)='GX ' 822 : LESOBL(2)='GY ' 823 : C ====== 824 : C 22.3 - Axisymetrie et Fourier 825 : C ====== 826 : ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN 827 : NBROBL=2 828 : SEGINI,NOMID 829 : LESOBL(1)='GR ' 830 : LESOBL(2)='GZ ' 831 : C ====== 832 : C 22.4 - Unidimensionnel (1D) 833 : C ====== 834 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN 835 : NBROBL=1 836 : SEGINI,NOMID 837 : IF (IFOU.LE.11) THEN 838 : LESOBL(1)='GX ' 839 : ELSE 840 : LESOBL(1)='GR ' 841 : ENDIF 842 : ENDIF 843 : C* Ajout du nom de l'INCOnnue (DDL) contenu dans TYMODE(1) 844 : DO i=1,NBROBL 845 : C* LESOBL(i)(3:8)=TYMODE(1)(1:6) 846 : LESOBL(i)(3:4)=TYMODE(1)(1:2) 847 : ENDDO 848 : 849 : ENDIF 850 : 851 : C Par DEFAUT : segment VIDE 852 : C =========================== 853 : 100 continue 854 : IF (NOMID.EQ.0) THEN 855 : SEGINI,NOMID 856 : ELSE 857 : if (ifomod.eq.6.and.(.not.dcmate)) then 858 : nbrfa0 = nbrfac 859 : NBRFAC = NBROBL + (nbrfa0*2) 860 : segadj nomid 861 : do imo = 1,nbrobl 862 : lesfac(nbrfa0 + imo)(2:4) = lesobl(imo)(1:3) 863 : lesfac(nbrfa0 + imo)(1:1) = 'I' 864 : enddo 865 : do imo = 1,nbrfa0 866 : lesfac(nbrfa0+nbrobl+imo)(2:4) = lesfac(imo)(1:3) 867 : lesfac(nbrfa0+nbrobl+imo)(1:1) = 'I' 868 : enddo 869 : endif 870 : ENDIF 871 : 872 : SEGDES,NOMID 873 : IPNOMC=NOMID 874 : 875 : RETURN 876 : END 877 : 878 : 879 : 880 : 881 : 882 :
© Cast3M 2003 - All rights reserved.
Disclaimer