1 : C IDGRAD SOURCE FANDEUR 12/07/05 21:15:56 7422 2 : 3 : C======================================================================= 4 : C= RECHERCHE DES NOMS DE COMPOSANTES DE GRADIENTS = 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 IDGRAD (MFR,IFOU,IPNOMC,NBROBL,NBRFAC) 19 : 20 : IMPLICIT INTEGER(I-N) 21 : 22 : -INC SMMODEL 23 : 24 : NBROBL=0 25 : NBRFAC=0 26 : NOMID=0 27 : 28 : C FORMULATION THERMOHYDRIQUE 29 : C ========================== 30 : IF(MFR.EQ.65) THEN 31 : C ===== 32 : C Tridimensionnel 33 : C ===== 34 : IF (IFOU.EQ.2) THEN 35 : NBROBL=9 36 : SEGINI,NOMID 37 : LESOBL(1)='PG,X ' 38 : LESOBL(2)='PG,Y ' 39 : LESOBL(3)='PG,Z ' 40 : LESOBL(4)='PC,X ' 41 : LESOBL(5)='PC,Y ' 42 : LESOBL(6)='PC,Z ' 43 : LESOBL(7)='T,X ' 44 : LESOBL(8)='T,Y ' 45 : LESOBL(9)='T,Z ' 46 : C ===== 47 : C Bidimensionnel PLAN (DP/CP/DPGE) 48 : C ===== 49 : ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 50 : NBROBL=6 51 : SEGINI,NOMID 52 : LESOBL(1)='PG,X ' 53 : LESOBL(2)='PG,Y ' 54 : LESOBL(3)='PC,X ' 55 : LESOBL(4)='PC,Y ' 56 : LESOBL(5)='T,X ' 57 : LESOBL(6)='T,Y ' 58 : C ===== 59 : C Axisymetrie et Fourier 60 : C ===== 61 : ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN 62 : NBROBL=6 63 : SEGINI,NOMID 64 : LESOBL(1)='PG,R ' 65 : LESOBL(2)='PG,Z ' 66 : LESOBL(3)='PC,R ' 67 : LESOBL(4)='PC,Z ' 68 : LESOBL(5)='T,R ' 69 : LESOBL(6)='T,Z ' 70 : ENDIF 71 : 72 : 73 : C 1 - Element JOINT JOI1 JOI2 74 : C ============================ 75 : 76 : ELSEIF (MFR.EQ.35.OR.MFR.EQ.75) THEN 77 : C ===== 78 : C 1.1 - Tridimensionnel 79 : C ===== 80 : IF (IFOU.EQ.2) THEN 81 : NBROBL=9 82 : SEGINI,NOMID 83 : LESOBL(1)='UX,X ' 84 : LESOBL(2)='UX,Y ' 85 : LESOBL(3)='UX,Z ' 86 : LESOBL(4)='UY,X ' 87 : LESOBL(5)='UY,Y ' 88 : LESOBL(6)='UY,Z ' 89 : LESOBL(7)='UZ,X ' 90 : LESOBL(8)='UZ,Y ' 91 : LESOBL(9)='UZ,Z ' 92 : C ===== 93 : C 1.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie 94 : C ===== 95 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 96 : NBROBL=4 97 : SEGINI,NOMID 98 : LESOBL(1)='UX,X ' 99 : LESOBL(2)='UX,Y ' 100 : LESOBL(3)='UY,X ' 101 : LESOBL(4)='UY,Y ' 102 : C ===== 103 : C 1.3 - Axisymetrie 104 : C ===== 105 : ELSE IF (IFOU.EQ.0) THEN 106 : NBROBL=4 107 : SEGINI,NOMID 108 : LESOBL(1)='UR,R ' 109 : LESOBL(2)='UR,Z ' 110 : LESOBL(3)='UZ,R ' 111 : LESOBL(4)='UZ,Z ' 112 : ENDIF 113 : 114 : C 2 - Elements MASSIFS (standard, incompressibles et XFEM) 115 : C ====================== 116 : C* ELSE IF (MFR.EQ.1) THEN 117 : ELSE IF (MFR.EQ.1 .OR. MFR.EQ.31 .OR. MFR.EQ.63) THEN 118 : C ===== 119 : C 2.1 - Bidimensionnel PLAN (CP/DP/DPGE) et Tridimensionnel 120 : C ===== 121 : IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.2) THEN 122 : NBROBL=9 123 : SEGINI,NOMID 124 : LESOBL(1)='UX,X ' 125 : LESOBL(2)='UX,Y ' 126 : LESOBL(3)='UX,Z ' 127 : LESOBL(4)='UY,X ' 128 : LESOBL(5)='UY,Y ' 129 : LESOBL(6)='UY,Z ' 130 : LESOBL(7)='UZ,X ' 131 : LESOBL(8)='UZ,Y ' 132 : LESOBL(9)='UZ,Z ' 133 : C ===== 134 : C 2.2 - Axisymetrie et Fourier et Unidimensionnel AXIS et SPHE (1D) 135 : C ===== 136 : ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1.OR.IFOU.EQ.4.OR.IFOU.EQ.5) THEN 137 : NBROBL=9 138 : SEGINI,NOMID 139 : LESOBL(1)='UR,R ' 140 : LESOBL(2)='UR,Z ' 141 : LESOBL(3)='UR,T ' 142 : LESOBL(4)='UZ,R ' 143 : LESOBL(5)='UZ,Z ' 144 : LESOBL(6)='UZ,T ' 145 : LESOBL(7)='UT,R ' 146 : LESOBL(8)='UT,Z ' 147 : LESOBL(9)='UT,T ' 148 : C ===== 149 : C 2.3 - Massif Unidimensionnel (1D) PLAN et AXIS/SPHE 150 : C ===== 151 : ELSE IF (MFR.EQ.1 .AND. (IFOU.GE.3.AND.IFOU.LE.15)) THEN 152 : NBROBL=9 153 : SEGINI,NOMID 154 : IF (IFOU.GE.3.AND.IFOU.LE.11) THEN 155 : LESOBL(1)='UX,X ' 156 : LESOBL(2)='UX,Y ' 157 : LESOBL(3)='UX,Z ' 158 : LESOBL(4)='UY,X ' 159 : LESOBL(5)='UY,Y ' 160 : LESOBL(6)='UY,Z ' 161 : LESOBL(7)='UZ,X ' 162 : LESOBL(8)='UZ,Y ' 163 : LESOBL(9)='UZ,Z ' 164 : ELSE 165 : LESOBL(1)='UR,R ' 166 : LESOBL(2)='UR,Z ' 167 : LESOBL(3)='UR,T ' 168 : LESOBL(4)='UZ,R ' 169 : LESOBL(5)='UZ,Z ' 170 : LESOBL(6)='UZ,T ' 171 : LESOBL(7)='UT,R ' 172 : LESOBL(8)='UT,Z ' 173 : LESOBL(9)='UT,T ' 174 : ENDIF 175 : ENDIF 176 : 177 : C 3 - Elements COQUES 178 : C ===================== 179 : ELSE IF (MFR.EQ.3) THEN 180 : C ===== 181 : C 3.1 - Tridimensionnel 182 : C ===== 183 : IF (IFOU.EQ.2) THEN 184 : NBROBL=9 185 : SEGINI,NOMID 186 : LESOBL(1)='UX,X' 187 : LESOBL(2)='UX,Y' 188 : LESOBL(3)='UX,Z' 189 : LESOBL(4)='UY,X' 190 : LESOBL(5)='UY,Y' 191 : LESOBL(6)='UY,Z' 192 : LESOBL(7)='UZ,X' 193 : LESOBL(8)='UZ,Y' 194 : LESOBL(9)='UZ,Z' 195 : C ===== 196 : C 3.2 - Axisymetrie, Fourier, Bidimensionnel PLAN (CP/DP) 197 : C ===== 198 : ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0.OR. 199 : . IFOU.EQ.-1.OR.IFOU.EQ.-2.or.ifou.eq.-3) THEN 200 : NBROBL=9 201 : SEGINI,NOMID 202 : LESOBL(1)='UU,S' 203 : LESOBL(2)='UU,T' 204 : LESOBL(3)='UU,N' 205 : LESOBL(4)='UV,S' 206 : LESOBL(5)='UV,T' 207 : LESOBL(6)='UV,N' 208 : LESOBL(7)='UW,S' 209 : LESOBL(8)='UW,T' 210 : LESOBL(9)='UW,N' 211 : ENDIF 212 : 213 : C 4 - Elements COQUES EPAISSES (COQ6,COQ8) 214 : C ============================== 215 : ELSE IF (MFR.EQ.5) THEN 216 : C ===== 217 : C 4.1 - Tridimensionnel 218 : C ===== 219 : IF (IFOU.EQ.2) THEN 220 : NBROBL=9 221 : SEGINI,NOMID 222 : LESOBL(1)='UX,X' 223 : LESOBL(2)='UX,Y' 224 : LESOBL(3)='UX,Z' 225 : LESOBL(4)='UY,X' 226 : LESOBL(5)='UY,Y' 227 : LESOBL(6)='UY,Z' 228 : LESOBL(7)='UZ,X' 229 : LESOBL(8)='UZ,Y' 230 : LESOBL(9)='UZ,Z' 231 : C ===== 232 : C 4.2 - Axisymetrie et Fourier 233 : C ===== 234 : ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN 235 : NBROBL=9 236 : SEGINI,NOMID 237 : LESOBL(1)='UR,R' 238 : LESOBL(2)='UR,Z' 239 : LESOBL(3)='UR,T' 240 : LESOBL(4)='UZ,R' 241 : LESOBL(5)='UZ,Z' 242 : LESOBL(6)='UZ,T' 243 : LESOBL(7)='UT,R' 244 : LESOBL(8)='UT,Z' 245 : LESOBL(9)='UT,T' 246 : ENDIF 247 : 248 : C 5 - Element de COQUE en CISAILLEMENT TRANSVERSE (tridimensionnel) 249 : C ================================================= 250 : ELSE IF (MFR.EQ.9) THEN 251 : IF (IFOU.EQ.2) THEN 252 : NBROBL=9 253 : SEGINI,NOMID 254 : LESOBL(1)='UX,X' 255 : LESOBL(2)='UX,Y' 256 : LESOBL(3)='UX,Z' 257 : LESOBL(4)='UY,X' 258 : LESOBL(5)='UY,Y' 259 : LESOBL(6)='UY,Z' 260 : LESOBL(7)='UZ,X' 261 : LESOBL(8)='UZ,Y' 262 : LESOBL(9)='UZ,Z' 263 : ENDIF 264 : 265 : C 6 - Formulation THERMIQUE MASSIF 266 : C ================================== 267 : ELSE IF (MFR.EQ.29) THEN 268 : C ===== 269 : C 6.1 - Tridimensionnel 270 : C ===== 271 : IF (IFOU.EQ.2) THEN 272 : NBROBL=3 273 : SEGINI,NOMID 274 : LESOBL(1)='T,X ' 275 : LESOBL(2)='T,Y ' 276 : LESOBL(3)='T,Z ' 277 : C ===== 278 : C 6.2 - Bidimensionnel PLAN (DP/CP/DPGE) 279 : C ===== 280 : ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 281 : NBROBL=2 282 : SEGINI,NOMID 283 : LESOBL(1)='T,X ' 284 : LESOBL(2)='T,Y ' 285 : C ===== 286 : C 6.3 - Axisymetrie et Fourier 287 : C ===== 288 : ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN 289 : NBROBL=2 290 : SEGINI,NOMID 291 : LESOBL(1)='T,R ' 292 : LESOBL(2)='T,Z ' 293 : C ===== 294 : C 6.4 - Unidimensionnel (1D) 295 : C ===== 296 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.11) THEN 297 : NBROBL=1 298 : SEGINI,NOMID 299 : LESOBL(1)='T,X ' 300 : ELSE IF (IFOU.GE.12.AND.IFOU.LE.15) THEN 301 : NBROBL=1 302 : SEGINI,NOMID 303 : LESOBL(1)='T,R ' 304 : ENDIF 305 : 306 : C 7 - Formulation THERMIQUE COQUE 307 : C ================================= 308 : ELSE IF (MFR.EQ.39) THEN 309 : C ===== 310 : C 7.1 - Tridimensionnel 311 : C ===== 312 : IF (IFOU.EQ.2) THEN 313 : NBROBL=2 314 : SEGINI,NOMID 315 : LESOBL(1)='T,X ' 316 : LESOBL(2)='T,Y ' 317 : C ===== 318 : C 7.2 - Axisymetrie et Fourier 319 : C ===== 320 : ELSE IF (IFOU.EQ.0) THEN 321 : NBROBL=1 322 : SEGINI,NOMID 323 : LESOBL(1)='T,Z ' 324 : ENDIF 325 : 326 : C 8 - Formulation MAGNETODYNAMIQUE COQUE (tridimensionnel) 327 : C ======================================== 328 : ELSE IF (MFR.EQ.69) THEN 329 : NBROBL=3 330 : SEGINI,NOMID 331 : LESOBL(1)='FC,X ' 332 : LESOBL(2)='FC,Y ' 333 : LESOBL(3)='FC,Z ' 334 : 335 : C 9 - formulation XFEM meca rupt 336 : C ===================================== 337 : c ELSE IF (MFR.EQ.63) THEN 338 : CTY - Element XQ4R (xfem meca rupture en 2D) 339 : c IF(IFOU.EQ.-2.OR.ifou.EQ.-1) then 340 : c NBROBL=4 341 : c NBRFAC=36 342 : c NBRFAC=20 343 : c SEGINI,NOMID 344 : c LESOBL(1)='UX,X' 345 : c LESOBL(2)='UX,Y' 346 : c LESOBL(3)='UY,X' 347 : c LESOBL(4)='UY,Y' 348 : c LESFAC(1)='AX,X' 349 : c LESFAC(2)='AX,Y' 350 : c LESFAC(3)='AY,X' 351 : c LESFAC(4)='AY,Y' 352 : c LESFAC(5)='B1XX' 353 : c LESFAC(6)='B1XY' 354 : c LESFAC(7)='B1YX' 355 : c LESFAC(8)='B1YY' 356 : c LESFAC(9)='C1XX' 357 : c LESFAC(10)='C1XY' 358 : c LESFAC(11)='C1YX' 359 : c LESFAC(12)='C1YY' 360 : c LESFAC(13)='D1XX' 361 : c LESFAC(14)='D1XY' 362 : c LESFAC(15)='D1YX' 363 : c LESFAC(16)='D1YY' 364 : c LESFAC(17)='E1XX' 365 : c LESFAC(18)='E1XY' 366 : c LESFAC(19)='E1YX' 367 : c LESFAC(20)='E1YY' 368 : c LESFAC(21)='B2XX' 369 : c LESFAC(22)='B2XY' 370 : c LESFAC(23)='B2YX' 371 : c LESFAC(24)='B2YY' 372 : c LESFAC(25)='C2XX' 373 : c LESFAC(26)='C2XY' 374 : c LESFAC(27)='C2YX' 375 : c LESFAC(28)='C2YY' 376 : c LESFAC(29)='D2XX' 377 : c LESFAC(30)='D2XY' 378 : c LESFAC(31)='D2YX' 379 : c LESFAC(32)='D2YY' 380 : c LESFAC(33)='E2XX' 381 : c LESFAC(34)='E2XY' 382 : c LESFAC(35)='E2YX' 383 : c LESFAC(36)='E2YY' 384 : c ELSE IF (IFOU.EQ.2) THEN 385 : CTY - Element XC8R (xfem meca rupture en 3D) 386 : c NBROBL=9 387 : c NBRFAC=45 388 : c NBRFAC=81 389 : c SEGINI,NOMID 390 : c LESOBL(1)='UX,X' 391 : c LESOBL(2)='UX,Y' 392 : c LESOBL(3)='UX,Z' 393 : c LESOBL(4)='UY,X' 394 : c LESOBL(5)='UY,Y' 395 : c LESOBL(6)='UY,Z' 396 : c LESOBL(7)='UZ,X' 397 : c LESOBL(8)='UZ,Y' 398 : c LESOBL(9)='UZ,Z' 399 : c LESFAC(1)='AX,X' 400 : c LESFAC(2)='AX,Y' 401 : c LESFAC(3)='AX,Z' 402 : c LESFAC(4)='AY,X' 403 : c LESFAC(5)='AY,Y' 404 : c LESFAC(6)='AY,Z' 405 : c LESFAC(7)='AZ,X' 406 : c LESFAC(8)='AZ,Y' 407 : c LESFAC(9)='AZ,Z' 408 : c LESFAC(10)='B1XX' 409 : c LESFAC(11)='B1XY' 410 : c LESFAC(12)='B1XZ' 411 : c LESFAC(13)='B1YX' 412 : c LESFAC(14)='B1YY' 413 : c LESFAC(15)='B1YZ' 414 : c LESFAC(16)='B1ZX' 415 : c LESFAC(17)='B1ZY' 416 : c LESFAC(18)='B1ZZ' 417 : c LESFAC(19)='C1XX' 418 : c LESFAC(20)='C1XY' 419 : c LESFAC(21)='C1XZ' 420 : c LESFAC(22)='C1YX' 421 : c LESFAC(23)='C1YY' 422 : c LESFAC(24)='C1YZ' 423 : c LESFAC(25)='C1ZX' 424 : c LESFAC(26)='C1ZY' 425 : c LESFAC(27)='C1ZZ' 426 : c LESFAC(28)='D1XX' 427 : c LESFAC(29)='D1XY' 428 : c LESFAC(30)='D1XZ' 429 : c LESFAC(31)='D1YX' 430 : c LESFAC(32)='D1YY' 431 : c LESFAC(33)='D1YZ' 432 : c LESFAC(34)='D1ZX' 433 : c LESFAC(35)='D1ZY' 434 : c LESFAC(36)='D1ZZ' 435 : c LESFAC(37)='E1XX' 436 : c LESFAC(38)='E1XY' 437 : c LESFAC(39)='E1XZ' 438 : c LESFAC(40)='E1YX' 439 : c LESFAC(41)='E1YY' 440 : c LESFAC(42)='E1YZ' 441 : c LESFAC(43)='E1ZX' 442 : c LESFAC(44)='E1ZY' 443 : c LESFAC(45)='E1ZZ' 444 : c LESFAC(46)='B2XX' 445 : c LESFAC(47)='B2XY' 446 : c LESFAC(48)='B2XZ' 447 : c LESFAC(49)='B2YX' 448 : c LESFAC(50)='B2YY' 449 : c LESFAC(51)='B2YZ' 450 : c LESFAC(52)='B2ZX' 451 : c LESFAC(53)='B2ZY' 452 : c LESFAC(54)='B2ZZ' 453 : c LESFAC(55)='C2XX' 454 : c LESFAC(56)='C2XY' 455 : c LESFAC(57)='C2XZ' 456 : c LESFAC(58)='C2YX' 457 : c LESFAC(59)='C2YY' 458 : c LESFAC(60)='C2YZ' 459 : c LESFAC(61)='C2ZX' 460 : c LESFAC(62)='C2ZY' 461 : c LESFAC(63)='C2ZZ' 462 : c LESFAC(64)='D2XX' 463 : c LESFAC(65)='D2XY' 464 : c LESFAC(66)='D2XZ' 465 : c LESFAC(67)='D2YX' 466 : c LESFAC(68)='D2YY' 467 : c LESFAC(69)='D2YZ' 468 : c LESFAC(70)='D2ZX' 469 : c LESFAC(71)='D2ZY' 470 : c LESFAC(72)='D2ZZ' 471 : c LESFAC(73)='E2XX' 472 : c LESFAC(74)='E2XY' 473 : c LESFAC(75)='E2XZ' 474 : c LESFAC(76)='E2YX' 475 : c LESFAC(77)='E2YY' 476 : c LESFAC(78)='E2YZ' 477 : c LESFAC(79)='E2ZX' 478 : c LESFAC(80)='E2ZY' 479 : c LESFAC(81)='E2ZZ' 480 : c ENDIF 481 : 482 : C 10- Formulation MECANIQUE POUTRES ET TUYAUX 483 : C ============================================ 484 : ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN 485 : C ====== 486 : C 10.1 - Tridimensionnel 487 : C ====== 488 : IF (IFOU.EQ.2) THEN 489 : NBROBL=3 490 : SEGINI,NOMID 491 : LESOBL(1)='U,S ' 492 : LESOBL(2)='V,S ' 493 : LESOBL(3)='W,S ' 494 : C ====== 495 : C 10.2 - Autres cas 496 : C ====== 497 : ELSE 498 : NBROBL=2 499 : SEGINI,NOMID 500 : LESOBL(1)='U,S ' 501 : LESOBL(2)='V,S ' 502 : ENDIF 503 : 504 : C 11 - Formulation ELECTROSTATIQUE (base MASSIF) 505 : C ================================== 506 : ELSE IF (MFR.EQ.71) THEN 507 : C ====== 508 : C 11.1 - Tridimensionnel 509 : C ====== 510 : IF (IFOU.EQ.2) THEN 511 : NBROBL=3 512 : SEGINI,NOMID 513 : LESOBL(1)='VE,X ' 514 : LESOBL(2)='VE,Y ' 515 : LESOBL(3)='VE,Z ' 516 : C ====== 517 : C 11.2 - Bidimensionnel PLAN (DP/CP/DPGE) 518 : C ====== 519 : ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 520 : NBROBL=2 521 : SEGINI,NOMID 522 : LESOBL(1)='VE,X ' 523 : LESOBL(2)='VE,Y ' 524 : C ====== 525 : C 11.3 - Axisymetrie et Fourier 526 : C ====== 527 : ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN 528 : NBROBL=2 529 : SEGINI,NOMID 530 : LESOBL(1)='VE,R ' 531 : LESOBL(2)='VE,Z ' 532 : C ====== 533 : C 11.4 - Unidimensionnel (1D) 534 : C ====== 535 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.11) THEN 536 : NBROBL=1 537 : SEGINI,NOMID 538 : LESOBL(1)='VE,X ' 539 : ELSE IF (IFOU.GE.12.AND.IFOU.LE.15) THEN 540 : NBROBL=1 541 : SEGINI,NOMID 542 : LESOBL(1)='VE,R ' 543 : ENDIF 544 : 545 : C 12 - Formulation DIFFUSION (base MASSIF) 546 : C ============================ 547 : ELSE IF (MFR.EQ.73) THEN 548 : C* L'INCOnnue (DDL) contenue dans TYMODE(1) est ajoutee ailleurs 549 : C ====== 550 : C 12.1 - Tridimensionnel 551 : C ====== 552 : IF (IFOU.EQ.2) THEN 553 : NBROBL=3 554 : SEGINI,NOMID 555 : C*8 LESOBL(1)=' ,X' 556 : C*8 LESOBL(2)=' ,Y' 557 : C*8 LESOBL(3)=' ,Z' 558 : LESOBL(1)=' ,X ' 559 : LESOBL(2)=' ,Y ' 560 : LESOBL(3)=' ,Z ' 561 : C ====== 562 : C 12.2 - Bidimensionnel PLAN (DP/CP/DPGE) 563 : C ====== 564 : ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 565 : NBROBL=2 566 : SEGINI,NOMID 567 : C*8 LESOBL(1)=' ,X' 568 : C*8 LESOBL(2)=' ,Y' 569 : LESOBL(1)=' ,X ' 570 : LESOBL(2)=' ,Y ' 571 : C ====== 572 : C 12.3 - Axisymetrie et Fourier 573 : C ====== 574 : ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN 575 : NBROBL=2 576 : SEGINI,NOMID 577 : C*8 LESOBL(1)=' ,R' 578 : C*8 LESOBL(2)=' ,Z' 579 : LESOBL(1)=' ,R ' 580 : LESOBL(2)=' ,Z ' 581 : C ====== 582 : C 12.4 - Unidimensionnel (1D) 583 : C ====== 584 : ELSE IF (IFOU.GE.3.AND.IFOU.LE.11) THEN 585 : NBROBL=1 586 : SEGINI,NOMID 587 : C*8 LESOBL(1)=' ,X' 588 : LESOBL(1)=' ,X ' 589 : ELSE IF (IFOU.GE.12.AND.IFOU.LE.15) THEN 590 : NBROBL=1 591 : SEGINI,NOMID 592 : C*8 LESOBL(1)=' ,R' 593 : LESOBL(1)=' ,R ' 594 : ENDIF 595 : 596 : C 13 - Formulation ZONE COHESIVE 597 : C ================================ 598 : ELSE IF (MFR.EQ.77) THEN 599 : C ===== 600 : C 13.1 - Tridimensionnel 601 : C ===== 602 : IF (IFOU.EQ.2) THEN 603 : NBROBL=9 604 : SEGINI,NOMID 605 : LESOBL(1)='AX,X ' 606 : LESOBL(2)='AX,Y ' 607 : LESOBL(3)='AX,Z ' 608 : LESOBL(4)='AY,X ' 609 : LESOBL(5)='AY,Y ' 610 : LESOBL(6)='AY,Z ' 611 : LESOBL(7)='AZ,X ' 612 : LESOBL(8)='AZ,Y ' 613 : LESOBL(9)='AZ,Z ' 614 : C ===== 615 : C 13.2 - Bidimensionnel PLAN (CP/DP/DPGE) 616 : C ===== 617 : ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN 618 : NBROBL=4 619 : SEGINI,NOMID 620 : LESOBL(1)='AX,X ' 621 : LESOBL(2)='AX,Y ' 622 : LESOBL(3)='AY,X ' 623 : LESOBL(4)='AY,Y ' 624 : ENDIF 625 : 626 : ENDIF 627 : 628 : C Par DEFAUT : segment VIDE 629 : C =========================== 630 : IF (NOMID.EQ.0) SEGINI,NOMID 631 : 632 : SEGDES,NOMID 633 : IPNOMC=NOMID 634 : 635 : RETURN 636 : END 637 : 638 :
© Cast3M 2003 - All rights reserved.
Disclaimer