1 : C ELQUOI SOURCE AM 12/02/06 21:15:27 7281 2 : 3 : C======================================================================= 4 : C= Sous-programme ELQUOI = 5 : C= --------------------- = 6 : C= = 7 : C= Ce sous programme fournit les caracteristiques des elements finis. = 8 : C= = 9 : C= Entrees : = 10 : C= --------- = 11 : C= NELE : numero de l'element fini = 12 : C= NBPOIN : nombre de points d'integration contraintes = 13 : C= multicouches ( COQ6 - COQ8 ) = 14 : C= INTTYP : type d'integration utilisee par l'operateur = 15 : C= Si l'on desire un segment d'integration pour un champ aux : = 16 : C= 1 = NOEUDS ( appel a RENOEU ) = 17 : C= 2 = points de GAUSS pour centre de GRAVITE et CHAMP CONSTANT = 18 : C= 3 = points de GAUSS pour la RIGIDITE = 19 : C= 4 = points de GAUSS pour la MASSE = 20 : C= 5 = points de GAUSS pour les CONTRAINTES = 21 : C= 6 = on ne calcule pas de segment d'integration = 22 : C= 7 = FACE = 23 : C= 8 = CENTREP1 = 24 : C= 9 = MSOMMET = 25 : C -10= en mecanique on veut calculer les 5 premiers et on les met = 26 : C dans infmod(3 ...) = 27 : C= Dans un autre cas, infell(11) renverra un pointeur sur un segment = 28 : C= d'integration de type RIGIDITE. = 29 : C= IMODEL : pointeur sur une zone elementaire d'un modele = 30 : C= SEGMENT ACTIF en entree et en sortie = 31 : C= Sorties : = 32 : C= --------- = 33 : C= retourne IPTR pointeur sur un segment INFO ACTIF contenant : = 34 : C= La dimension de infell est fixee par NINFOS = 35 : C= infell( 1) numero de l'element fini = 36 : C= infell( 2) nombre de points d'integration en contraintes = 37 : C= multicouche = 38 : C= infell( 3) nombre de points d'integration pour la masse = 39 : C= infell( 4) nombre de points d'integration pour SIGMA = 40 : C= BSIGMA et KSIGMA = 41 : C= infell( 5) nombre de caracteristiques = 42 : C= infell( 6) nombre de points d'integration pour la RIGIDITE = 43 : C= infell( 7) longueur d'un tableau de travail pour l'element = 44 : C= infell( 8) nombre de fonctions de forme = 45 : C= infell( 9) nbre de d.d.l. dans la matrice de RIGIDITE = 46 : C= infell(10) taille de la matrice de Hooke = 47 : C= infell(11) pointeur sur le segment d'integration = 48 : C= infell(12) pointeur sur le 2nd segment d'integration des = 49 : C= elements homogeneises ou fluides (COQ6 ou COQ8) = 50 : C= infell(13) numero de la formulation de l'element fini = 51 : C= = 1 MASSIF = 52 : C= voir aussi = 3 COQUE = 53 : C= NUMMFR.eso = 5 COQUE EPAISSE ( COQ6 , COQ8 ) = 54 : C= = 7 POUTRE = 55 : C= = 9 COQUE AVEC CISAILLEMENT TRANSVERSE = 56 : C= = 11 LIQUIDE = 57 : C= = 13 TUYAU = 58 : C= = 15 LINESPRING = 59 : C= = 17 TUYAU FISSURE = 60 : C= = 19 RACCORD MASSIF = 61 : C= = 21 RACCORD COQUE = 62 : C= = 23 SURFACE LIBRE = 63 : C= = 25 MEMBRANE = 64 : C= = 26 point support MODAL = 65 : C= = 27 UNIAXIALE = 66 : C= = 28 point support STATIQUE = 67 : C= = 29 THERMIQUE = 68 : C= = 31 INCOMPRESSIBLES = 69 : C= = 33 POREUX = 70 : C= = 35 JOINT = 71 : C= = 37 HOMOGENEISE = 72 : C= = 39 TUYO = 73 : C= = 41 TUYAU ACOUSTIQUE PURE = 74 : C= = 43 RACCORD TUYAU FLUIDE = 75 : C= = 45 HYBRIDE = 76 : C= = 47 FIBRE = 77 : C= = 49 BARRE 3D EXCENTREE = 78 : C= = 51 LIASON A 2 NOEUDS = 79 : C= = 53 JOINT CISAILLEMENT (2D) = 80 : C= = 55 JOINT GENERALISE = 81 : C= = 57 POREUX TYPE Q = 82 : C= = 59 POREUX TYPE R = 83 : C= = 61 MACRO ELEMENT = 84 : C= = 63 XFEM meca rupt (2D/3D) = 85 : C= = 65 THERMOHYDRIQUE = 86 : C= = 69 MAGNETODYNAMIQUE = 87 : C= = 71 ELECTROSTATIQUE = 88 : C= = 73 DIFFUSION = 89 : C= = 75 JOINT UNIDIMENSIONNEL = 90 : C= = 77 ZONE COHESIVE = 91 : C= infell(14) numero de l'element geometrique associe (NUMGEO) = 92 : C= infell(15) nombre maximal de d.d.l. par noeud = 93 : C= infell(16) nombre de composantes de contraintes ou de deform. = 94 : C======================================================================= 95 : 96 : SUBROUTINE ELQUOI(NELE,NBPOIN,INTTYP,IPTR,IMODEL) 97 : 98 : IMPLICIT INTEGER(I-N) 99 : IMPLICIT REAL*8(A-H,O-Z) 100 : 101 : -INC CCOPTIO 102 : -INC CCGEOME 103 : -INC SMMODEL 104 : -INC SMELEME 105 : -INC SMLMOTS 106 : 107 : C Segment (type LISTENTI) contenant les informations sur un element 108 : SEGMENT INFO 109 : INTEGER INFELL(JG) 110 : ENDSEGMENT 111 : POINTEUR INFO1.INFO,INFO2.INFO 112 : 113 : C======================================================================= 114 : C= INTEGR : SAVE local des pointeurs sur les segments d'integration = 115 : C= CAS GENERAL = 116 : C= NELFI : Nombre de types d'element FINI (= LNOMTP de bdata.eso) = 117 : C= NBMODE : Nombre de MODEs de calcul (= nb valeurs possibles IFOMOD) = 118 : C= TYINTE : Type de l'INTEGRATION utilise par le calcul = 119 : C= (correspond a la valeur de INTTYP=infell(11) en entree) = 120 : C= NTYNTE : Nombre de TYPEs de segment D'INTEGRATION possibles = 121 : C= (voir la description de INTTYP ci-dessus) = 122 : C= NINFOS : Dimension (JG) du tableau infell (segment INFO) = 16 = 123 : C= On sauvegarde dans INTEGR(I,J) le pointeur sur un segment d'inte- = 124 : C= gration correspondant a l'element fini de type I et une integration = 125 : C= de type J. = 126 : C======================================================================= 127 : C Valeurs de IFOMOD : -1 0 1 2 3 4 5 128 : C PLAN AXIS FOUR TRID UNID UNID UNID 129 : C PLAN AXIS SPHE 130 : PARAMETER ( NTYNTE = 6 , NELFI = 265 , NBMODE = 7 ) 131 : PARAMETER ( NINTEG = NTYNTE * NELFI * NBMODE , NINFOS = 16 ) 132 : INTEGER INTEGR(NELFI,NTYNTE,NBMODE),TYINTE 133 : logical dcmate,dcmat2 134 : SAVE INTEGR 135 : 136 : DATA INTEGR / NINTEG*0 / 137 : 138 : dcmate = .false. 139 : dcmat2 = .false. 140 : C On se place sur la bonne tranche du tableau INTEGR 141 : IBMODE = IFOMOD + 2 142 : C Recuperation des donnees de infell en entree 143 : MELE = NELE 144 : NBGS = NBPOIN 145 : TYINTE = INTTYP 146 : if(inttyp.eq.-10) tyinte=0 147 : IELE = NUMGEO(MELE) 148 : MFR = NUMMFR(MELE) 149 : do im = 1,matmod(/2) 150 : if (matmod(im).eq.'IMPEDANCE') then 151 : dcmate = .true. 152 : meleme = imamod 153 : segact meleme*nomod 154 : ityp1 = itypel 155 : if (ityp1.eq.1) mele = 45 156 : iele = numgeo(mele) 157 : mfr = nummfr(mele) 158 : endif 159 : enddo 160 : C Nombre de points d'integration dans l'epaisseur (stocke dans IMODEL) 161 : NPINT=0 162 : IF (INFMOD(/1).NE.0) NPINT=INFMOD(1) 163 : c IF(MELE.EQ.263) NPINT=0 164 : C Initialisation 165 : IPT1 = 0 166 : IPT2 = 0 167 : JG=NINFOS 168 : SEGINI,INFO 169 : 170 : C=---------------------------------------------------------------------= 171 : C= REMPLISSAGE DU TABLEAU infell 172 : C=---------------------------------------------------------------------= 173 : C Remplissage de infell(1) : numero de l'element fini 174 : C -------------------------- 175 : INFELL(1)=NELE 176 : 177 : C Remplissage de infell(8) : 178 : C -------------------------- 179 : C Sauf cas particuliers, le nombre de fonctions de forme de l'element 180 : C fini est egal au nombre de noeuds de l'element geometrique support. 181 : INFELL(8) = NBNNE(IELE) 182 : C Massifs : TRIx-QUAx 4 a 11, CUxx a PRxx 14 a 17, TExx-PYxx 23 a 26 183 : IF (MELE.GE.4.AND.MELE.LE.11) THEN 184 : INFELL(8) = 3 185 : ELSE IF (MELE.GE.14.AND.MELE.LE.17) THEN 186 : infell(8) = 3 187 : ELSE IF (MELE.GE.23.AND.MELE.LE.26) THEN 188 : infell(8) = 3 189 : C Polygones : POxx 111 a 122 190 : ELSE IF (MELE.GE.111.AND.MELE.LE.122) THEN 191 : IF (MELE.EQ.111) infell(8) = 3 192 : IF (MELE.EQ.112) infell(8) = 4 193 : IF (MELE.EQ.113) infell(8) = 5 194 : IF (MELE.EQ.114) infell(8) = 6 195 : IF (MELE.EQ.115) infell(8) = 7 196 : IF (MELE.EQ.116) infell(8) = 8 197 : IF (MELE.EQ.117) infell(8) = 9 198 : IF (MELE.EQ.118) infell(8) = 10 199 : IF (MELE.EQ.119) infell(8) = 11 200 : IF (MELE.EQ.120) infell(8) = 12 201 : IF (MELE.EQ.121) infell(8) = 13 202 : IF (MELE.EQ.122) infell(8) = 14 203 : C Elements QxRI 183 a 184 204 : ELSE IF (MELE.EQ.183.OR.MELE.EQ.184) THEN 205 : infell(8) = 3 206 : CTY Element XQ4R (xfem meca rupture en 2D) 207 : ELSEIF(MELE.eq.263) THEN 208 : infell(8)=40 209 : CTY Element XC8R (xfem meca rupture en 3D) 210 : ELSE IF(MELE.eq.264) THEN 211 : infell(8)=80 212 : ENDIF 213 : 214 : C Elements Navier Stokes : cas P1CENTRE 215 : IF(TYINTE.EQ.8) THEN 216 : C cas MACRO 2D QUAF 2D --> TRI3 217 : IF(MELE.EQ.217.OR.MELE.EQ.224.OR.MELE.EQ.231 218 : & .OR.MELE.EQ.238.OR.MELE.EQ.245.OR.MELE.EQ.252) 219 : & infell(8) = 3 220 : C cas MACRO 2D QUAF 2D --> QUA4 221 : IF(MELE.EQ.218.OR.MELE.EQ.225.OR.MELE.EQ.232 222 : & .OR.MELE.EQ.239.OR.MELE.EQ.246.OR.MELE.EQ.253) 223 : & infell(8) = 3 224 : C cas MACRO 3D QUAF 3D 225 : C IF((MELE.GE.226.AND.MELE.LE.229).OR. 226 : C & (MELE.GE.247.AND.MELE.LE.250)) infell(8) = 4 227 : ENDIF 228 : C Elements Navier Stokes : cas MSOMMET 229 : IF(TYINTE.EQ.9) THEN 230 : C cas MACRO QUAF 1D --> SEG2 231 : IF(MELE.EQ.216.OR.MELE.EQ.223.OR.MELE.EQ.230 232 : & .OR.MELE.EQ.237.OR.MELE.EQ.244.OR.MELE.EQ.251) IELE = 2 233 : C cas MACRO QUAF 2D --> TRI3 234 : IF(MELE.EQ.217.OR.MELE.EQ.224.OR.MELE.EQ.231 235 : & .OR.MELE.EQ.238.OR.MELE.EQ.245.OR.MELE.EQ.252) IELE = 4 236 : C cas MACRO QUAF 2D --> QUA4 237 : IF(MELE.EQ.218.OR.MELE.EQ.225.OR.MELE.EQ.232 238 : & .OR.MELE.EQ.239.OR.MELE.EQ.246.OR.MELE.EQ.253) IELE = 8 239 : C cas MACRO QUAF 3D --> CUB8 240 : IF(MELE.EQ.233.OR.MELE.EQ.254) IELE = 14 241 : C cas MACRO QUAF 3D --> PRI6 242 : IF(MELE.EQ.234.OR.MELE.EQ.255) IELE = 16 243 : C cas MACRO QUAF 3D --> TET4 244 : IF(MELE.EQ.235.OR.MELE.EQ.256) IELE = 23 245 : C cas MACRO QUAF 3D --> PYR5 246 : IF(MELE.EQ.236.OR.MELE.EQ.257) IELE = 25 247 : infell(8) = NBNNE(IELE) 248 : infell(14)= IELE 249 : ENDIF 250 : 251 : C Pour les elements xxxP (79 a 83), TUYO (96), JOPx (108 a 110), xxxQ 252 : C (173 a 177), xxxR (178 a 182), JORx (185 a 187) et JOQx(188 a 190), 253 : C infell(8) est affecte lors du remplissage complet de infell. 254 : 255 : C Remplissage de infell(13) et infell(14) : 256 : C ----------------------------------------- 257 : C infell(13) : numero de la FORMULATION - MFR=NUMMFR(NELE) 258 : infell(13) = MFR 259 : do im = 1,matmod(/2) 260 : if (MFR.eq.27.and.matmod(im).eq.'MODAL') infell(13) = 26 261 : if (MFR.eq.27.and.matmod(im).eq.'STATIQUE') infell(13) = 28 262 : enddo 263 : MFR = infell(13) 264 : mfr3 = nummfr(nele) 265 : C mfr2 non nul en cas de formulation massive non mecanique 266 : mfr2 = 0 267 : if (MFR.EQ.26.OR.MFR.EQ.28) mfr3 = mfr 268 : C Cas particuliers : elements HYxx (99 a 103) 269 : C MFR = NUMMFR(mele) = 45 mais infell(13)=1 (MASSIF) 270 : IF (MELE.GE.99.AND.MELE.LE.103) infell(13)=1 271 : C ELECTROSTATIQUE : MFR=71=infell(13) et NUMMFR(mele)=1 (MASSIF) 272 : IF (FORMOD(1).EQ.'ELECTROSTATIQUE') mfr2=71 273 : C DIFFUSION : MFR=73=infell(13) et NUMMFR(mele)=1 (MASSIF) 274 : IF (FORMOD(1).EQ.'DIFFUSION') mfr2=73 275 : IF (mfr2.NE.0) infell(13) = mfr2 276 : C 277 : C infell(14) : numero de l'element GEOMETRIQUE - IELE=NUMGEO(NELE) 278 : infell(14) = IELE 279 : 280 : C Remplissage de infell(10), infell(15) et infell(16) : 281 : C ----------------------------------------------------- 282 : C infell(10) : dimension de la matrice de Hooke 283 : C infell(15) : nombre maximal de ddl par noeud 284 : C infell(16) : nombre de composantes de contraintes et deformations 285 : 286 : C MASSIFs (2D) : TRIx-QUAx 4 a 11 + POLYgones : POLY 111 a 122 + 287 : C Elements QxRI 183 a 184 288 : IF ((MELE.GE.4.AND.MELE.LE.11).OR.(MELE.GE.111.AND.MELE.LE.122) 289 : . .OR.MELE.EQ.183.OR.MELE.EQ.184) THEN 290 : IF (IFOUR.EQ.1) THEN 291 : infell(10)=6 292 : infell(15)=3 293 : infell(16)=6 294 : ELSE 295 : infell(10)=4 296 : infell(15)=2 297 : infell(16)=4 298 : ENDIF 299 : C MASSIFs (3D) : CUxx-PRxx 14 a 17, TExx-PYxx 23 a 26 300 : ELSE IF ((MELE.GE.14.AND.MELE.LE.17).OR. 301 : . (MELE.GE.23.AND.MELE.LE.26)) THEN 302 : infell(10)=6 303 : infell(15)=3 304 : infell(16)=6 305 : C Elements LIQUIDE : Lxxx 35 a 40 306 : ELSE IF (MELE.GE.35.AND.MELE.LE.40) THEN 307 : infell(10)=5 308 : infell(15)=2 309 : infell(16)=0 310 : C Elements de RACCORD FLUIDE-STRUCTURE : 311 : C RAC2 12, LIAx 18 a 19, RACO 47, LICO 55, LIC4 94 312 : ELSE IF (MELE.EQ.12.OR.MELE.EQ.18.OR.MELE.EQ.19.OR.MELE.EQ.47.OR. 313 : . MELE.EQ.55.OR.MELE.EQ.94) THEN 314 : infell(10)=6 315 : C Elements de SURFACE LIBRE : LSU2 48, LSU3 53, LSU4 54 316 : ELSE IF (MELE.EQ.48.OR.MELE.EQ.53.OR.MELE.EQ.54) THEN 317 : infell(10)=6 318 : infell(15)=3 319 : C Elements THERMIQUEs : CVxx 57 a 62, THxx 63 a 68 320 : ELSE IF (MELE.GE.57.AND.MELE.LE.68) THEN 321 : IF (MELE.LE.60) THEN 322 : infell(10)=2 323 : infell(15)=1 324 : infell(16)=2 325 : ELSE IF (MELE.GE.61) THEN 326 : infell(10)=3 327 : infell(15)=1 328 : infell(16)=3 329 : ENDIF 330 : C Elements MASSIFs IMCOMPRESSIBLES : ICxx 69 a 78 331 : ELSE IF (MELE.GE.69.AND.MELE.LE.78) THEN 332 : IF (MELE.LE.72) THEN 333 : infell(10)=4 334 : C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(10)=4 335 : infell(15)=2 336 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(15)=3 337 : infell(16)=4 338 : C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(16)=4 339 : ELSE IF (MELE.GE.73) THEN 340 : infell(10)=6 341 : infell(15)=3 342 : infell(16)=6 343 : ENDIF 344 : C Elements MASSIFs POREUX : xxxP 79 a 83 345 : ELSE IF (MELE.GE.79.AND.MELE.LE.83) THEN 346 : IF (MELE.LE.80) THEN 347 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN 348 : infell(10)=6 349 : infell(15)=4 350 : infell(16)=7 351 : ELSE 352 : infell(10)=4 353 : infell(15)=3 354 : infell(16)=5 355 : ENDIF 356 : ELSE IF (MELE.GE.81) THEN 357 : infell(10)=6 358 : infell(15)=4 359 : infell(16)=7 360 : ENDIF 361 : C Elements MASSIFs POREUX de type Q : xxxQ 173 a 177 362 : ELSE IF (MELE.GE.173.AND.MELE.LE.177) THEN 363 : IF (MELE.LE.174) THEN 364 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN 365 : infell(10)=8 366 : infell(15)=5 367 : infell(16)=8 368 : ELSE 369 : infell(10)=6 370 : infell(15)=4 371 : infell(16)=6 372 : ENDIF 373 : ELSE IF (MELE.GE.175) THEN 374 : infell(10)=6 375 : infell(15)=5 376 : infell(16)=8 377 : ENDIF 378 : C Elements MASSIFs POREUX de type R : xxxR 178 a 182 379 : ELSE IF(MELE.GE.178.AND.MELE.LE.182) THEN 380 : IF (MELE.LE.179) THEN 381 : infell(10)=7 382 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(10)=9 383 : infell(15)=5 384 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(15)=6 385 : infell(16)=7 386 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(16)=9 387 : ELSE IF (MELE.GE.180) THEN 388 : infell(10)=6 389 : infell(15)=6 390 : infell(16)=9 391 : ENDIF 392 : C Elements JOINTS : JOxx 85 a 90 393 : ELSE IF (MELE.GE.85.AND.MELE.LE.90) THEN 394 : IF (MELE.LE.86) THEN 395 : infell(10)=2 396 : infell(16)=2 397 : ELSE IF (MELE.GE.87) THEN 398 : infell(10)=3 399 : infell(16)=3 400 : ENDIF 401 : C 402 : C MACRO ELEMENT CISAILLEMENT FLEXION 403 : C 404 : ELSE IF(MELE.EQ.191) THEN 405 : infell(5)=3 406 : infell(7)=300 407 : GO TO 258 408 : 409 : C Elements JOINTs CISAILLEMENT et GENERALISES : JCxx-JGxx 168 a 172 410 : ELSE IF (MELE.GE.168.AND.MELE.LE.172) THEN 411 : infell(7)=300 412 : infell(9)=infell(8)*IDIM 413 : infell(15)=IDIM 414 : IF (MELE.EQ.168.OR.MELE.EQ.169) THEN 415 : infell(10)=2 416 : infell(16)=2 417 : ELSE IF (MELE.EQ.170) THEN 418 : infell(10)=4 419 : infell(16)=4 420 : ELSE 421 : infell(10)=6 422 : infell(16)=6 423 : ENDIF 424 : C Elements HOMOGENEISES : TRIH 92, QUAH 126, CUBH 127, TRH6 157 425 : ELSE IF (MELE.EQ.92.OR.MELE.EQ.126.OR.MELE.EQ.127.OR. 426 : . MELE.EQ.157) THEN 427 : infell(10)=10 428 : infell(16)=0 429 : C Elements JOINTS POREUX : JOPx 108 a 110 430 : ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN 431 : infell(7)=100 432 : C Elements JOINTS POREUX de type Q : JOQx 185 a 187 433 : ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN 434 : infell(7)=100 435 : C Elements JOINTS POREUX de type R : JORx 188 a 190 436 : ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN 437 : infell(7)=100 438 : C Elements MASSIFS THERMIQUE/MECANIQUE en 1D : T1Dx-M1Dx 191 a 194 439 : ELSE IF (MELE.GE.191.AND.MELE.LE.194) THEN 440 : infell(10)=3 441 : infell(15)=1 442 : infell(16)=3 443 : ENDIF 444 : C Cas particuliers des autres formulations massives (mfr2 > 0) 445 : IF (mfr2.NE.0) THEN 446 : infell(10)=IDIM 447 : infell(15)=1 448 : infell(16)=IDIM 449 : ENDIF 450 : C Pour les autres elements, les valeurs de infell(10,15,16) sont 451 : C eventuellement affectees (0 par defaut) lors du branchement (GOTO) 452 : C ci-dessous qui effectue le remplissage complet du tableau infell. 453 : 454 : C Remplissage complet du tableau infell suivant l'element fini 455 : C -------------------------------------------------------------- 456 : IF(MELE.GE.1.AND.MELE.LE.194) THEN 457 : C CABL SEG2 SEG3 TRI3 TRI4 TRI6 TRI7 QUA4 QUA5 QUA8 QUA9,RAC2 458 : GOTO ( 300, 2, 3, 4, 300, 6, 300, 8, 300, 10, 300, 12 459 : C RAC3 CUB8 CU20 PRI6 PR15 LIA3 LIA4 LIA6 LIA8 MULT TET4 TE10 460 : . , 300, 14, 15, 16, 17, 18, 19, 300, 300, 300, 23, 24 461 : C PYR5 PY13 COQ3 DKT POUT LISP FAC3 FAC4 FAC6 FAC8 LTR3 LQU4 462 : . , 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36 463 : C LCU8 LPR6 LTE4 LPY5 COQ8 TUYA TUFI COQ2 POI1 BARR RACO LSU2 464 : . , 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48 465 : C COQ4 LISM COF3 RES2 LSU3 LSU4 LICO COQ6 CVS2 CVS3 CVT3 CVT6 466 : . , 49, 50, 51, 300, 53, 54, 55, 56, 57, 58, 59, 60 467 : C CVQ4 CVQ8 THP5 TH13 THP6 TH15 THC8 TH20 ICT3 ICQ4 ICT6 ICQ8 468 : . , 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72 469 : C ICC8 ICT4 ICP6 IC20 IC10 IC15 TRIP QUAP CUBP TETP PRIP TIMO 470 : . , 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84 471 : C JOI2 JOI3 JOT3 JOI4 JOI6 JOI8 LISC TRIH DST LIC4 CERC TUYO 472 : . , 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96 473 : C LSE2 LITU HYT3 HYQ4 HYT4 HYP6 HYC8 TRIS QUAS POIS FOR3 JOP3 474 : . , 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 300, 108 475 : C JOP6 JOP8 POL3 POL4 POL5 POL6 POL7 POL8 POL9 PO10 PO11 PO12 476 : . , 109, 110, 111, 111, 111, 111, 111, 111, 111, 111, 111, 111 477 : C PO13 PO14 BAR3 BAEX LIA2 QUAH CUBH ROT3 SEF2 TRF3 QUF4 CUF8 478 : . , 111, 111, 123, 124, 125, 126, 127, 128, 300, 300, 300, 300 479 : C PRF6 TEF4 PYF5 MSE3 MTR6 MQU9 MC27 MP18 MT10 MP14 SEF3 TRF7 480 : . , 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300 481 : C QUF9 CF27 PF21 TF15 PF19 SEG6 TR21 QU36 C216 P126 TE56 PY91 482 : . , 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300 483 : C TRH6 BSE2 BTR4 BQU5 BCU9 BPR7 BTE5 BPY6 FRO4 SEGS POJS JCT3 484 : . , 157, 300, 300, 300, 300, 300, 300, 300, 300, 166, 167, 168 485 : C JCI4 JGI2 JGT3 JGI4 TRIQ QUAQ CUBQ TETQ PRIQ TRIR QUAR CUBR 486 : . , 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180 487 : C TETR PRIR Q4RI Q8RI JOQ3 JOQ6 JOQ8 JOR3 JOR6 JOR8 T1D2 T1D3 488 : . , 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192 489 : C M1D2 M1D3 490 : . , 191, 192 ), MELE 491 : ELSEIF(MELE.GE.195.AND.MELE.LE.268) THEN 492 : C LC03 LC07 LC09 LC27 LC21 LC15 LC19 493 : GOTO ( 195, 196, 197, 300, 300, 300, 300 494 : C LS03 LS07 LS09 LS27 LS21 LS15 LS19 495 : . , 195, 196, 197, 300, 300, 300, 300 496 : C BS03 BS07 BS09 BS27 BS21 BS15 BS19 497 : . , 195, 196, 197, 300, 300, 300, 300 498 : C MC03 MC07 MC09 MC27 MC21 MC15 MC19 499 : . , 216, 217, 218, 300, 300, 300, 300 500 : C M103 M107 M109 M127 M121 M115 M119 501 : . , 223, 224, 225, 300, 300, 300, 300 502 : C MS03 MS07 MS09 MS27 MS21 MS15 MS19 503 : . , 223, 224, 225, 300, 300, 300, 300 504 : C QC03 QC07 QC09 QC27 QC21 QC15 QC19 505 : . , 237, 238, 239, 300, 300, 300, 300 506 : C Q103 Q107 Q109 Q127 Q121 Q115 Q119 507 : . , 237, 238, 239, 300, 300, 300, 300 508 : C QS03 QS07 QS09 QS27 QS21 QS15 QS19 CIFL SURE 509 : . , 237, 238, 239, 300, 300, 300, 300 , 300 ,300 510 : C SHB8 caf2 caf3 xq4r xc8r joi1 ZCO2 ZCO3 ZCO4 511 : . ,260 , 300, 300, 263, 264, 265 , 266, 267, 268), (MELE-194) 512 : ENDIF 513 : C 514 : C Remplissage de infell pour chaque element 515 : C ------------------------------------------- 516 : C Element SEG2 517 : 2 infell(2)=2 518 : infell(4)=2 519 : infell(6)=2 520 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=3 521 : infell(7)=0 522 : infell(9)=4 523 : IF (IFOUR.EQ.1) infell(9)=6 524 : infell(10)=1 525 : infell(15)=2 526 : IF (IFOUR.EQ.1) infell(15)=3 527 : infell(16)=1 528 : GOTO 300 529 : C Element SEG3 530 : 3 infell(2)=3 531 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(2)=4 532 : infell(4)=3 533 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(4)=4 534 : infell(6)=3 535 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4 536 : infell(7)=0 537 : infell(9)=6 538 : IF (IFOUR.EQ.1) infell(9)=9 539 : infell(10)=1 540 : infell(15)=2 541 : IF (IFOUR.EQ.1) infell(15)=3 542 : infell(16)=1 543 : GOTO 300 544 : C Element TRI3 545 : 4 infell(2)=1 546 : infell(3)=4 547 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=7 548 : infell(4)=1 549 : infell(6)=1 550 : C IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4 551 : infell(9)=6 552 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=9 553 : GOTO 300 554 : C Element TRI6 555 : 6 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 556 : infell(2)=7 557 : infell(4)=7 558 : infell(6)=7 559 : ELSE 560 : infell(2)=4 561 : infell(4)=4 562 : infell(6)=4 563 : ENDIF 564 : infell(3)=7 565 : infell(9)=12 566 : IF (IFOUR.EQ.-3) infell(9)=15 567 : IF (IFOUR.EQ.1) infell(9)=18 568 : GOTO 300 569 : C Element QUA4 570 : 8 infell(2)=4 571 : infell(3)=4 572 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 573 : infell(4)=4 574 : C** infell(4)=5 CORRECTION MILL LE 8/8/90 575 : infell(6)=4 576 : C** IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9 577 : infell(9)=8 578 : IF (IFOUR.EQ.-3) infell(9)=11 579 : IF (IFOUR.EQ.1) infell(9)=12 580 : GOTO 300 581 : C Element QUA8 582 : 10 infell(2)=9 583 : infell(3)=9 584 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=16 585 : infell(4)=9 586 : infell(6)=9 587 : C IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 588 : infell(9)=16 589 : IF (IFOUR.EQ.-3) infell(9)=19 590 : IF (IFOUR.EQ.1) infell(9)=24 591 : GOTO 300 592 : C Element RAC2 (raccord a 4 noeuds en 2 dimensions) 593 : 12 infell(2)=2 594 : infell(4)=2 595 : infell(5)=2 596 : infell(6)=2 597 : infell(9)=8 598 : IF (IFOUR.EQ.1) infell(9)=10 599 : infell(15)=3 600 : GOTO 300 601 : C Element CUB8 602 : 14 infell(2)=8 603 : infell(3)=8 604 : infell(4)=8 605 : infell(6)=8 606 : infell(9)=24 607 : GOTO 300 608 : C Element CU20 609 : 15 infell(2)=27 610 : infell(3)=27 611 : infell(4)=27 612 : infell(6)=27 613 : infell(9)=60 614 : GOTO 300 615 : C Element PRI6 616 : 16 infell(2)=6 617 : infell(3)=6 618 : infell(4)=6 619 : infell(6)=6 620 : infell(9)=18 621 : GOTO 300 622 : C Element PR15 623 : 17 infell(2)=12 624 : infell(3)=12 625 : infell(4)=12 626 : infell(6)=12 627 : infell(9)=45 628 : GOTO 300 629 : C Element LIA3 (raccord a 6 noeudes en 3 dimensions) 630 : 18 infell(2)=4 631 : infell(4)=4 632 : infell(5)=3 633 : infell(6)=4 634 : infell(9)=15 635 : infell(13)=19 636 : infell(15)=3 637 : GOTO 300 638 : C Element LIA4 (raccord a 8 noeuds en 3 dimensions) 639 : 19 infell(2)=4 640 : infell(4)=4 641 : infell(5)=3 642 : infell(6)=4 643 : infell(9)=20 644 : infell(15)=3 645 : GOTO 300 646 : C Element TET4 647 : 23 infell(2)=1 648 : infell(3)=4 649 : infell(4)=1 650 : infell(6)=1 651 : infell(9)=12 652 : GOTO 300 653 : C Element TE10 654 : 24 infell(2)=5 655 : infell(3)=5 656 : infell(4)=5 657 : infell(6)=5 658 : infell(9)=30 659 : GOTO 300 660 : C Element PYR5 661 : 25 infell(2)=5 662 : infell(3)=5 663 : infell(4)=5 664 : infell(6)=5 665 : infell(9)=15 666 : GOTO 300 667 : C Element PY13 668 : 26 infell(2)=27 669 : infell(3)=27 670 : infell(4)=27 671 : infell(6)=27 672 : infell(9)=39 673 : GOTO 300 674 : C Element COQ3 675 : 27 infell(2)=3 676 : infell(4)=1 677 : infell(5)=3 678 : infell(6)=1 679 : infell(7)=600 680 : infell(9)=18 681 : infell(10)=6 682 : infell(15)=6 683 : infell(16)=6 684 : GOTO 300 685 : C Element DKT 686 : 28 infell(2)=3 687 : infell(3)=3 688 : infell(5)=3 689 : infell(7)=600 690 : infell(9)=18 691 : infell(15)=6 692 : IF (NPINT.EQ.0) THEN 693 : C Element DKT COQUE GLOBALE 694 : infell(4)=3 695 : infell(6)=3 696 : infell(10)=6 697 : infell(16)=6 698 : ELSE 699 : C Element DKT COQUE INTEGREE dans l'epaisseur 700 : infell(4)=3*NPINT 701 : infell(6)=3*NPINT 702 : infell(10)=6 703 : infell(16)=4 704 : ENDIF 705 : GOTO 300 706 : C Element POUTRE 707 : 29 infell(2)=2 708 : infell(4)=2 709 : infell(6)=2 710 : infell(7)=200 711 : IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN 712 : C Element POUTRE 2D plan 713 : infell(5)=6 714 : infell(9)=6 715 : infell(10)=3 716 : infell(15)=3 717 : infell(16)=3 718 : ELSE 719 : C Element POUTRE 3D 720 : infell(5)=12 721 : infell(9)=12 722 : infell(10)=6 723 : infell(15)=6 724 : infell(16)=6 725 : ENDIF 726 : GOTO 300 727 : C Element LINE SPRING LISP 728 : 30 infell(2)=3 729 : infell(13)=15 730 : infell(6)=3 731 : infell(4)=3 732 : infell(16)=6 733 : infell(9)=24 734 : infell(7)=100 735 : infell(10)=6 736 : infell(15)=6 737 : infell(5)=5 738 : GOTO 300 739 : C Element FAC3 pour les PRESSIONS uniquement 740 : 31 infell(2)=1 741 : infell(4)=1 742 : infell(6)=1 743 : infell(9)=9 744 : infell(15)=3 745 : GOTO 300 746 : C Element FAC4 pour les PRESSIONS uniquement 747 : 32 infell(2)=4 748 : infell(4)=4 749 : infell(6)=4 750 : infell(9)=12 751 : infell(15)=3 752 : GOTO 300 753 : C Element FAC6 pour les PRESSIONS uniquement 754 : 33 infell(2)=4 755 : infell(4)=4 756 : infell(6)=4 757 : infell(9)=18 758 : infell(15)=3 759 : GOTO 300 760 : C Element FAC8 pour les PRESSIONS uniquement 761 : 34 infell(2)=9 762 : infell(4)=9 763 : infell(6)=9 764 : infell(9)=24 765 : infell(15)=3 766 : GOTO 300 767 : C Element LTR3 : element LIQUIDE triangle a 3 noeuds 768 : 35 infell(2)=4 769 : infell(3)=4 770 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=7 771 : infell(4)=4 772 : infell(6)=4 773 : infell(9)=6 774 : GOTO 300 775 : C Element LQU4 : element LIQUIDE quadrangle a 4 noeuds 776 : 36 infell(2)=4 777 : infell(3)=4 778 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 779 : infell(4)=4 780 : infell(6)=4 781 : infell(9)=8 782 : GOTO 300 783 : C Element LCU8 : Element LIQUIDE CUBE a 8 noeuds 784 : 37 infell(2)=8 785 : infell(3)=8 786 : infell(4)=8 787 : infell(6)=8 788 : infell(9)=16 789 : GOTO 300 790 : C Element LPR6 : Element LIQUIDE PRISME a 6 noeuds 791 : 38 infell(2)=6 792 : infell(3)=6 793 : infell(4)=6 794 : infell(6)=6 795 : infell(9)=12 796 : GOTO 300 797 : C Element LTE4 : Element LIQUIDE TETRAEDRE a 4 noeuds 798 : 39 infell(2)=4 799 : infell(3)=4 800 : infell(4)=4 801 : infell(6)=4 802 : infell(9)=8 803 : GOTO 300 804 : C Element LPY5 : Element LIQUIDE PYRAMIDE a 5 noeuds 805 : 40 infell(2)=5 806 : infell(3)=5 807 : infell(4)=5 808 : infell(6)=5 809 : infell(9)=10 810 : GOTO 300 811 : C Element COQ8 : Element COQUE EPAISSE a 8 noeuds 812 : 41 infell(2)=NBGS 813 : IF (NBGS.EQ.0) infell(2)=8 814 : infell(3)=8 815 : infell(4)=NBGS 816 : IF (NBGS.EQ.0) infell(4)=8 817 : infell(5)=2 818 : infell(6)=8 819 : infell(7)=1700 820 : infell(9)=48 821 : infell(10)=5 822 : infell(15)=6 823 : infell(16)=5 824 : GOTO 300 825 : C Element TUYA : Element TUYAU 826 : 42 infell(2)=2 827 : infell(4)=2 828 : infell(5)=7 829 : infell(6)=2 830 : infell(7)=200 831 : infell(9)=12 832 : infell(10)=6 833 : infell(15)=6 834 : infell(16)=6 835 : GOTO 300 836 : C Element TUFI : Element TUYAU FISSURE 837 : 43 infell(2)=1 838 : infell(3)=0 839 : infell(4)=1 840 : infell(5)=9 841 : infell(6)=1 842 : infell(7)=200 843 : infell(9)=12 844 : infell(10)=6 845 : infell(15)=6 846 : infell(16)=8 847 : GOTO 300 848 : C Element COQ2 849 : 44 infell(2)=3 850 : infell(3)=3 851 : infell(4)=3 852 : infell(5)=3 853 : infell(6)=3 854 : infell(7)=368 855 : infell(9)=8 856 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(9)=6 857 : IF (IFOUR.EQ.-3) infell(9)=9 858 : infell(10)=6 859 : IF (IFOUR.LE.0) infell(10)=4 860 : infell(15)=4 861 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(15)=3 862 : infell(16)=6 863 : IF (IFOUR.LE.0) infell(16)=4 864 : GOTO 300 865 : C Element POI1 (MODEs en DEFOrmations PLANes GENEralisees) 866 : C Valeurs a adapter pour certains MODE UNID (1D) 867 : 45 infell(2)=1 868 : infell(3)=1 869 : infell(4)=1 870 : infell(5)=1 871 : infell(6)=1 872 : infell(7)=50 873 : infell(15)=2 874 : IF (IFOUR.EQ.11) THEN 875 : infell(9)=3 876 : infell(10)=2 877 : infell(16)=2 878 : ELSE 879 : infell(9)=5 880 : IF (IFOUR.GE.3.AND.IFOUR.LE.15) infell(9)=2 881 : infell(10)=1 882 : infell(16)=1 883 : ENDIF 884 : if (mfr.eq.26.or.mfr.eq.28) infell(9) = 1 885 : GOTO 300 886 : C Element BARRE 887 : 46 infell(2)=2 888 : infell(4)=2 889 : infell(5)=1 890 : infell(6)=2 891 : infell(7)=200 892 : infell(9)=6 893 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(9)=4 894 : IF (IFOUR.EQ.-3) infell(9)=7 895 : infell(10)=1 896 : infell(15)=3 897 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(15)=2 898 : infell(16)=1 899 : GOTO 300 900 : C Element RACO : Element de RACCORD LIQUIDE-COQUE a 4 noeuds en 2D 901 : 47 infell(2)=1 902 : infell(4)=1 903 : infell(5)=2 904 : infell(6)=1 905 : infell(9)=10 906 : IF (IFOUR.EQ.1) infell(9)=12 907 : infell(15)=4 908 : GOTO 300 909 : C Element LSU2 : Element de SURFACE LIBRE a 2 noeuds en 2D 910 : 48 infell(2)=2 911 : infell(4)=2 912 : infell(6)=2 913 : infell(9)=6 914 : GOTO 300 915 : C Element COQ4 916 : 49 infell(2)=5 917 : infell(3)=5 918 : infell(4)=5 919 : infell(5)=3 920 : infell(6)=5 921 : infell(7)=100 922 : infell(9)=24 923 : infell(10)=8 924 : infell(15)=6 925 : infell(16)=8 926 : GOTO 300 927 : C Element LINE SPRING LISM 928 : 50 infell(2)=1 929 : infell(4)=1 930 : infell(5)=5 931 : infell(6)=3 932 : infell(7)=100 933 : infell(9)=24 934 : infell(10)=6 935 : infell(15)=6 936 : infell(16)=6 937 : GOTO 300 938 : C Element COF3 939 : 51 infell(2)=3 940 : infell(3)=3 941 : infell(4)=3 942 : infell(5)=3 943 : infell(6)=3 944 : infell(9)=12 945 : infell(7)=368 946 : infell(10)=6 947 : infell(15)=4 948 : infell(16)=6 949 : GOTO 300 950 : C Element LSU3 : Element de SURFACE LIBRE a 3 noeuds en 3D 951 : 53 infell(2)=4 952 : infell(4)=4 953 : infell(6)=4 954 : infell(9)=9 955 : GOTO 300 956 : C Element LSU4 : Element de SURFACE LIBRE a 4 noeuds en 3D 957 : 54 infell(2)=4 958 : infell(4)=4 959 : infell(6)=4 960 : infell(9)=12 961 : GOTO 300 962 : C Element LICO : Element de RACCORD LIQUIDE-COQUE a 6 noeuds en 3D 963 : 55 infell(2)=3 964 : infell(4)=3 965 : infell(5)=3 966 : infell(6)=3 967 : infell(9)=24 968 : infell(15)=6 969 : GOTO 300 970 : C Element COQ6 : Element COQUE EPAISSE a 6 noeuds 971 : 56 infell(2)=NBGS 972 : IF(NBGS.EQ.0) infell(2)=6 973 : infell(3)=6 974 : infell(4)=NBGS 975 : IF (NBGS.EQ.0) infell(4)=6 976 : infell(5)=2 977 : infell(6)=6 978 : infell(9)=36 979 : infell(7)=1000 980 : infell(10)=5 981 : infell(15)=6 982 : infell(16)=5 983 : GOTO 300 984 : C Element THT3 (thermique) 985 : 57 infell(2)=1 986 : infell(6)=1 987 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4 988 : infell(9)=3 989 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=9 990 : GOTO 300 991 : C Element THT6 992 : 58 infell(2)=4 993 : infell(4)=4 994 : infell(6)=4 995 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=7 996 : infell(9)=6 997 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=18 998 : GOTO 300 999 : C Element THQ4 1000 : 59 infell(2)=4 1001 : infell(6)=4 1002 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9 1003 : infell(9)=4 1004 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=12 1005 : GOTO 300 1006 : C Element THQ8 1007 : 60 infell(2)=9 1008 : infell(6)=9 1009 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 1010 : infell(9)=8 1011 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=24 1012 : GOTO 300 1013 : C Element THT4 1014 : 61 infell(2)=1 1015 : infell(6)=1 1016 : infell(9)=4 1017 : GOTO 300 1018 : C Element TH10 1019 : 62 infell(2)=5 1020 : infell(6)=5 1021 : infell(9)=10 1022 : GOTO 300 1023 : C Element THP5 1024 : 63 infell(2)=5 1025 : infell(6)=5 1026 : infell(9)=5 1027 : GOTO 300 1028 : C Element TH13 1029 : 64 infell(2)=27 1030 : infell(6)=27 1031 : infell(9)=13 1032 : GOTO 300 1033 : C Element THP6 1034 : 65 infell(2)=2 1035 : infell(6)=2 1036 : infell(9)=6 1037 : GOTO 300 1038 : C Element TH15 1039 : 66 infell(2)=12 1040 : infell(6)=12 1041 : infell(9)=15 1042 : GOTO 300 1043 : C Element THC8 1044 : 67 infell(2)=8 1045 : infell(6)=8 1046 : infell(9)=8 1047 : GOTO 300 1048 : C Element TH20 1049 : 68 infell(2)=27 1050 : infell(6)=27 1051 : infell(9)=20 1052 : GOTO 300 1053 : C Element ICT3 1054 : 69 infell(2)=1 1055 : infell(3)=4 1056 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=7 1057 : infell(4)=1 1058 : infell(6)=1 1059 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4 1060 : infell(9)=6 1061 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=9 1062 : GOTO 300 1063 : C Element ICQ4 1064 : 70 infell(2)=4 1065 : infell(3)=4 1066 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 1067 : infell(4)=4 1068 : C** infell(4)=5 CORRECTION MILL LE 8/8/90 1069 : infell(6)=4 1070 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9 1071 : infell(9)=8 1072 : IF (IFOUR.EQ.-3) infell(9)=11 1073 : IF (IFOUR.EQ.1) infell(9)=12 1074 : GOTO 300 1075 : C Element ICT6 1076 : 71 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 1077 : infell(2)=7 1078 : infell(4)=7 1079 : infell(6)=7 1080 : ELSE 1081 : infell(2)=4 1082 : infell(4)=4 1083 : infell(6)=4 1084 : ENDIF 1085 : infell(3)=7 1086 : infell(9)=12 1087 : IF (IFOUR.EQ.-3) infell(9)=15 1088 : IF (IFOUR.EQ.1) infell(9)=18 1089 : GOTO 300 1090 : C Element ICQ8 1091 : 72 infell(2)=9 1092 : infell(3)=9 1093 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=16 1094 : infell(4)=9 1095 : infell(6)=9 1096 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 1097 : infell(9)=16 1098 : IF (IFOUR.EQ.-3) infell(9)=19 1099 : IF (IFOUR.EQ.1) infell(9)=24 1100 : GOTO 300 1101 : C Element ICC8 1102 : 73 infell(2)=8 1103 : infell(3)=8 1104 : infell(4)=8 1105 : infell(6)=8 1106 : infell(9)=24 1107 : GOTO 300 1108 : C Element ICT4 1109 : 74 infell(2)=4 1110 : infell(3)=5 1111 : infell(4)=4 1112 : infell(6)=4 1113 : infell(9)=12 1114 : GOTO 300 1115 : C Element ICP6 1116 : 75 infell(2)=2 1117 : infell(3)=6 1118 : infell(4)=2 1119 : infell(6)=2 1120 : infell(9)=18 1121 : GOTO 300 1122 : C Element IC20 1123 : 76 infell(2)=27 1124 : infell(3)=27 1125 : infell(4)=27 1126 : infell(6)=27 1127 : infell(9)=60 1128 : GOTO 300 1129 : C Element IC10 1130 : 77 infell(2)=5 1131 : infell(3)=5 1132 : infell(4)=5 1133 : infell(6)=5 1134 : infell(9)=30 1135 : GOTO 300 1136 : C Element IC15 1137 : 78 infell(2)=12 1138 : infell(3)=12 1139 : infell(4)=12 1140 : infell(6)=12 1141 : infell(9)=45 1142 : GOTO 300 1143 : C Element TRIP 1144 : 79 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 1145 : infell(2)=7 1146 : infell(4)=7 1147 : infell(6)=7 1148 : ELSE 1149 : infell(2)=4 1150 : infell(4)=4 1151 : infell(6)=4 1152 : ENDIF 1153 : infell(3)=7 1154 : infell(8)=9 1155 : infell(9)=15 1156 : IF (IFOUR.EQ.-3.OR.IFOUR.EQ.1) infell(9)=21 1157 : GOTO 300 1158 : C Element QUAP 1159 : 80 infell(2)=9 1160 : infell(3)=9 1161 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 1162 : infell(4)=9 1163 : infell(6)=9 1164 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 1165 : infell(8)=12 1166 : infell(9)=20 1167 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=28 1168 : GOTO 300 1169 : C Element CUBP 1170 : 81 infell(2)=27 1171 : infell(3)=27 1172 : infell(4)=27 1173 : infell(6)=27 1174 : infell(8)=28 1175 : infell(9)=68 1176 : GOTO 300 1177 : C Element TETP 1178 : 82 infell(2)=5 1179 : infell(3)=5 1180 : infell(4)=5 1181 : infell(6)=5 1182 : infell(8)=14 1183 : infell(9)=34 1184 : GOTO 300 1185 : C Element PRIP 1186 : 83 infell(2)=12 1187 : infell(3)=12 1188 : infell(4)=12 1189 : infell(6)=12 1190 : infell(8)=21 1191 : infell(9)=51 1192 : GOTO 300 1193 : C Element TIMO : Element POUTRE de TIMOSHENCKO 1194 : 84 infell(2)=1 1195 : infell(4)=1 1196 : infell(6)=1 1197 : infell(7)=200 1198 : IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN 1199 : C Element TIMO 2D plan 1200 : infell(5)=6 1201 : infell(9)=6 1202 : infell(10)=3 1203 : infell(15)=3 1204 : infell(16)=3 1205 : ELSE 1206 : C Element TIMO 3D 1207 : infell(5)=12 1208 : infell(9)=12 1209 : infell(10)=6 1210 : infell(15)=6 1211 : infell(16)=6 1212 : ENDIF 1213 : GOTO 300 1214 : C Element JOI2 1215 : 85 infell(2)=4 1216 : infell(3)=2 1217 : infell(4)=2 1218 : infell(5)=0 1219 : infell(6)=2 1220 : infell(7)=300 1221 : infell(9)=8 1222 : GOTO 300 1223 : C Element JOI3 1224 : 86 infell(2)=3 1225 : infell(3)=2 1226 : infell(4)=3 1227 : infell(5)=0 1228 : infell(6)=3 1229 : infell(7)=100 1230 : infell(9)=12 1231 : GOTO 300 1232 : C Element JOT3 1233 : 87 infell(2)=3 1234 : infell(3)=3 1235 : infell(4)=3 1236 : infell(5)=0 1237 : infell(6)=3 1238 : infell(7)=100 1239 : infell(9)=18 1240 : GOTO 300 1241 : C Element JOI4 1242 : 88 infell(2)=4 1243 : infell(3)=4 1244 : infell(4)=4 1245 : infell(5)=0 1246 : infell(6)=4 1247 : infell(7)=100 1248 : infell(9)=24 1249 : GOTO 300 1250 : C Element JOI6 1251 : 89 infell(2)=5 1252 : infell(4)=5 1253 : infell(6)=4 1254 : infell(9)=36 1255 : GOTO 300 1256 : C Element JOI8 1257 : 90 infell(2)=9 1258 : infell(4)=9 1259 : infell(6)=9 1260 : infell(9)=48 1261 : GOTO 300 1262 : C Element LINE SPRING LISC 1263 : 91 infell(2)=3 1264 : infell(4)=3 1265 : infell(5)=5 1266 : infell(6)=5 1267 : infell(7)=100 1268 : infell(9)=36 1269 : infell(10)=2 1270 : infell(15)=6 1271 : infell(16)=6 1272 : GOTO 300 1273 : C Element TRIH 1274 : 92 infell(2)=4 1275 : infell(6)=4 1276 : infell(7)= 0 1277 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN 1278 : infell(3)=7 1279 : infell(5)=5 1280 : IF (IFOUR.EQ.0) infell(9)=12 1281 : IF (IFOUR.EQ.1) infell(9)=18 1282 : infell(15)=6 1283 : ELSE 1284 : infell(3)=4 1285 : infell(5)=5 1286 : infell(9)=12 1287 : infell(15)=4 1288 : ENDIF 1289 : GOTO 300 1290 : C Element DST 1291 : 93 infell(2)=3 1292 : infell(3)=3 1293 : infell(4)=3 1294 : infell(5)=3 1295 : infell(6)=3 1296 : infell(7)=600 1297 : infell(9)=18 1298 : infell(10)=8 1299 : infell(15)=6 1300 : infell(16)=8 1301 : GOTO 300 1302 : C Element LIC4 : Element de RACCORD LIQUIDE-COQUE a 8 noeuds en 3D 1303 : 94 infell(2)=4 1304 : infell(3)=4 1305 : infell(4)=4 1306 : infell(5)=3 1307 : infell(6)=4 1308 : infell(9)=32 1309 : infell(15)=6 1310 : GOTO 300 1311 : C Element CERCE 1312 : 95 infell(2)=1 1313 : infell(3)=1 1314 : infell(4)=1 1315 : infell(5)=1 1316 : infell(6)=1 1317 : infell(7)=50 1318 : infell(9)=2 1319 : IF (IFOUR.EQ.1) infell(9)=3 1320 : infell(10)=1 1321 : infell(15)=2 1322 : IF (IFOUR.EQ.1) infell(15)=3 1323 : infell(16)=1 1324 : GOTO 300 1325 : C Element TUYO 1326 : 96 NTETA=40 1327 : infell(2)=2*NTETA 1328 : infell(3)=2*NTETA 1329 : infell(4)=2*NTETA 1330 : infell(5)=7 1331 : infell(6)=2*NTETA 1332 : infell(7)=500 1333 : infell(8)=2+4 1334 : infell(9)=12+2*9 1335 : infell(10)=6 1336 : infell(15)=6+9 1337 : infell(16)=6 1338 : GOTO 300 1339 : C Element LSE2 : TUYAU ACOUSTIQUE PURE 1340 : 97 infell(2)=2 1341 : infell(3)=2 1342 : infell(4)=2 1343 : infell(5)=7 1344 : infell(6)=2 1345 : infell(7)=200 1346 : infell(9)=4 1347 : infell(10)=2 1348 : infell(15)=2 1349 : infell(16)=2 1350 : GOTO 300 1351 : C Element LITU : COUPLAGE LIQUIDE TUYAU 1352 : 98 infell(2)=2 1353 : infell(3)=2 1354 : infell(4)=2 1355 : infell(5)=6 1356 : infell(6)=2 1357 : infell(7)=200 1358 : infell(9)=10 1359 : infell(10)=2 1360 : infell(15)=5 1361 : infell(16)=5 1362 : GOTO 300 1363 : C Element HYT3 : Element HYBRIDE TRI3 1364 : 99 infell(3)=4 1365 : infell(6)=4 1366 : GOTO 300 1367 : C Element HYQ4 : Element HYBRIDE QUA4 1368 : 100 infell(3)=4 1369 : infell(6)=4 1370 : GOTO 300 1371 : C Element HYBRIDE TET4 1372 : 101 infell(3)=4 1373 : infell(6)=4 1374 : GOTO 300 1375 : C Element HYBRIDE PRI6 1376 : 102 infell(3)=6 1377 : infell(6)=6 1378 : GOTO 300 1379 : C Element HYBRIDE CUB8 1380 : 103 infell(3)=8 1381 : infell(6)=8 1382 : GOTO 300 1383 : C Element TRIS (TRI3-SECTION) 1384 : 104 infell(2)=4 1385 : infell(3)=4 1386 : infell(4)=4 1387 : infell(5)=2 1388 : infell(6)=4 1389 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 1390 : infell(10)=2 1391 : infell(16)=2 1392 : ELSE IF (IFOUR.EQ.2) THEN 1393 : infell(10)=3 1394 : infell(16)=3 1395 : ENDIF 1396 : GOTO 300 1397 : C Element QUAS (QUA4-SECTION) 1398 : 105 infell(2)=4 1399 : infell(3)=4 1400 : infell(4)=4 1401 : infell(5)=2 1402 : infell(6)=4 1403 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 1404 : infell(10)=2 1405 : infell(16)=2 1406 : ELSE IF (IFOUR.EQ.2) THEN 1407 : infell(10)=3 1408 : infell(16)=3 1409 : ENDIF 1410 : GOTO 300 1411 : C Element POIS (POI1-SECTION) 1412 : 106 infell(2)=1 1413 : infell(3)=1 1414 : infell(4)=1 1415 : infell(5)=1 1416 : CMILL???? infell(5)=3 1417 : infell(6)=1 1418 : infell(10)=1 1419 : infell(16)=1 1420 : GOTO 300 1421 : C Element JOP3 : JOINTS POREUX 1422 : 108 infell(2)=3 1423 : infell(3)=3 1424 : infell(4)=3 1425 : infell(6)=3 1426 : infell(8)=12 1427 : infell(9)=18 1428 : infell(10)=2 1429 : infell(15)=3 1430 : infell(16)=3 1431 : GOTO 300 1432 : C Element JOP6 : JOINTS POREUX 1433 : 109 infell(2)=7 1434 : infell(3)=7 1435 : infell(4)=7 1436 : infell(6)=7 1437 : infell(8)=21 1438 : infell(9)=45 1439 : infell(10)=3 1440 : infell(15)=4 1441 : infell(16)=4 1442 : GOTO 300 1443 : C Element JOP8 : JOINTS POREUX 1444 : 110 infell(2)=9 1445 : infell(3)=9 1446 : infell(4)=9 1447 : infell(6)=9 1448 : infell(8)=28 1449 : infell(9)=60 1450 : infell(10)=3 1451 : infell(15)=4 1452 : infell(16)=4 1453 : GOTO 300 1454 : C Elements POLYGONE a N cotes 1455 : 111 NBSH=infell(8) 1456 : infell(2)=NBSH 1457 : infell(3)=NBSH 1458 : infell(4)=NBSH 1459 : infell(6)=NBSH 1460 : infell(9)=2*NBSH 1461 : IF (IFOUR.EQ.1) infell(9)=3*NBSH 1462 : IF (IFOUR.EQ.-3) infell(9)=2*NBSH+3 1463 : GOTO 300 1464 : C Element BAR3 (barre a 3 noeuds) 1465 : 123 NPOINT=3 1466 : infell(2)=NPOINT 1467 : infell(4)=NPOINT 1468 : infell(5)=1 1469 : infell(6)=NPOINT 1470 : infell(7)=200 1471 : infell(9)=9 1472 : IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) infell(9)=6 1473 : infell(10)=1 1474 : infell(15)=3 1475 : IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) infell(15)=2 1476 : infell(16)=1 1477 : GOTO 300 1478 : C Element BAEX (barre excentree 3D a 2 noeuds) 1479 : 124 infell(2)=2 1480 : infell(3)=2 1481 : infell(4)=2 1482 : infell(5)=6 1483 : infell(6)=2 1484 : infell(7)=700 1485 : infell(9)=12 1486 : infell(10)=1 1487 : infell(15)=6 1488 : infell(16)=1 1489 : GOTO 300 1490 : C Element LIA2 : element de liaison a 2 noeuds (6 ddl par noeuds) 1491 : 125 infell(2)=2 1492 : infell(4)=2 1493 : infell(5)=9 1494 : infell(6)=2 1495 : infell(7)=200 1496 : infell(9)=12 1497 : infell(10)=6 1498 : infell(15)=6 1499 : infell(16)=6 1500 : GOTO 300 1501 : C Element QUAH 1502 : 126 infell(2)=4 1503 : infell(3)=4 1504 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 1505 : infell(5)=5 1506 : infell(6)=4 1507 : infell(7)=0 1508 : IF (IFOUR.EQ.0) infell(9)=16 1509 : IF (IFOUR.EQ.1) infell(9)=24 1510 : infell(15)=6 1511 : GOTO 300 1512 : C Element CUBH 1513 : 127 infell(2)=8 1514 : infell(3)=27 1515 : infell(5)=5 1516 : infell(6)=27 1517 : infell(7)=0 1518 : infell(9)=48 1519 : infell(15)=6 1520 : GOTO 300 1521 : C Element ROT3 (1 pt GAUSS pr MASSE, 7 pts GAUSS pr RIGIDITE) 1522 : 128 infell(3)=1 1523 : infell(6)=7 1524 : infell(9)=3 1525 : infell(15)=1 1526 : GOTO 300 1527 : C Element TR6H 1528 : 157 infell(2)=4 1529 : infell(6)=6 1530 : infell(7)=0 1531 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN 1532 : infell(3)=7 1533 : infell(5)=5 1534 : IF (IFOUR.EQ.0) infell(9)=12 1535 : IF (IFOUR.EQ.1) infell(9)=18 1536 : infell(15)=6 1537 : ELSE 1538 : infell(3)=6 1539 : infell(5)=5 1540 : infell(9)=24 1541 : infell(15)=4 1542 : ENDIF 1543 : GOTO 300 1544 : C Element SEGS (SEG2-SECTION 2D) 1545 : 166 infell(2)=2 1546 : infell(3)=2 1547 : infell(4)=2 1548 : infell(5)=3 1549 : infell(6)=2 1550 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 1551 : infell(10)=2 1552 : infell(16)=2 1553 : ELSE IF (IFOUR.EQ.2) THEN 1554 : infell(10)=3 1555 : infell(16)=3 1556 : ENDIF 1557 : GOTO 300 1558 : C Element POJS (POI1-SECTION) 1559 : 167 infell(2)=1 1560 : infell(3)=1 1561 : infell(4)=1 1562 : infell(5)=3 1563 : infell(6)=1 1564 : IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN 1565 : infell(10)=2 1566 : infell(16)=2 1567 : ELSE IF (IFOUR.EQ.2) THEN 1568 : infell(10)=3 1569 : infell(16)=3 1570 : ENDIF 1571 : GOTO 300 1572 : C Element JCT3 1573 : 168 infell(2)=3 1574 : infell(3)=3 1575 : infell(4)=3 1576 : infell(6)=3 1577 : GOTO 300 1578 : C Element JCI4 1579 : 169 infell(2)=4 1580 : infell(3)=4 1581 : infell(4)=4 1582 : infell(6)=4 1583 : GOTO 300 1584 : C Element JGI2 1585 : 170 infell(2)=4 1586 : infell(3)=2 1587 : infell(4)=2 1588 : infell(6)=2 1589 : IF (IFOUR.EQ.-3) infell(9)=infell(9)+3 1590 : infell(15)=2 1591 : GOTO 300 1592 : C Element JGT3 1593 : 171 infell(2)=3 1594 : infell(3)=3 1595 : infell(4)=3 1596 : infell(6)=3 1597 : GOTO 300 1598 : C Element JGI4 1599 : 172 infell(2)=4 1600 : infell(3)=4 1601 : infell(4)=4 1602 : infell(6)=4 1603 : GOTO 300 1604 : C Element TRIQ 1605 : 173 infell(3)=7 1606 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 1607 : infell(2)=7 1608 : infell(4)=7 1609 : infell(6)=7 1610 : ELSE 1611 : infell(2)=4 1612 : infell(4)=4 1613 : infell(6)=4 1614 : ENDIF 1615 : infell(8)=9 1616 : infell(9)=18 1617 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=24 1618 : GOTO 300 1619 : C Element QUAQ 1620 : 174 infell(2)=9 1621 : infell(3)=9 1622 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 1623 : infell(4)=9 1624 : infell(6)=9 1625 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 1626 : infell(8)=12 1627 : infell(9)=24 1628 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=32 1629 : GOTO 300 1630 : C Element CUBQ 1631 : 175 infell(2)=27 1632 : infell(3)=27 1633 : infell(4)=27 1634 : infell(6)=27 1635 : infell(8)=28 1636 : infell(9)=76 1637 : GOTO 300 1638 : C Element TETQ 1639 : 176 infell(2)=5 1640 : infell(3)=5 1641 : infell(4)=5 1642 : infell(6)=5 1643 : infell(8)=14 1644 : infell(9)=38 1645 : GOTO 300 1646 : C Element PRIQ 1647 : 177 infell(2)=12 1648 : infell(3)=12 1649 : infell(4)=12 1650 : infell(6)=12 1651 : infell(8)=21 1652 : infell(9)=57 1653 : GOTO 300 1654 : C Element TRIR 1655 : 178 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 1656 : infell(2)=7 1657 : infell(6)=7 1658 : infell(4)=7 1659 : ELSE 1660 : infell(2)=4 1661 : infell(4)=4 1662 : infell(6)=4 1663 : ENDIF 1664 : infell(3)=7 1665 : infell(8)=9 1666 : infell(9)=21 1667 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=27 1668 : GOTO 300 1669 : C Element QUAR 1670 : 179 infell(2)=9 1671 : infell(3)=9 1672 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 1673 : infell(4)=9 1674 : infell(6)=9 1675 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 1676 : infell(8)=12 1677 : infell(9)=28 1678 : IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=36 1679 : GOTO 300 1680 : C Element CUBR 1681 : 180 infell(2)=27 1682 : infell(3)=27 1683 : infell(4)=27 1684 : infell(6)=27 1685 : infell(8)=28 1686 : infell(9)=84 1687 : GOTO 300 1688 : C Element TETR 1689 : 181 infell(2)=5 1690 : infell(3)=5 1691 : infell(4)=5 1692 : infell(6)=5 1693 : infell(8)=14 1694 : infell(9)=42 1695 : GOTO 300 1696 : C Element PRIR 1697 : 182 infell(2)=12 1698 : infell(3)=12 1699 : infell(4)=12 1700 : infell(6)=12 1701 : infell(8)=21 1702 : infell(9)=61 1703 : GOTO 300 1704 : C Element Q4RI (QUA4 with 1x1 Gauss points) 1705 : 183 infell(2)=4 1706 : infell(3)=1 1707 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9 1708 : infell(4)=1 1709 : C** infell(4)=5 CORRECTION MILL LE 8/8/90 1710 : infell(6)=1 1711 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9 1712 : infell(9)=8 1713 : IF (IFOUR.EQ.-3) infell(9)=11 1714 : IF (IFOUR.EQ.1) infell(9)=12 1715 : GOTO 300 1716 : C Element Q8RI (QUA8 with 2x2 Gauss points) 1717 : 184 infell(2)=9 1718 : infell(3)=4 1719 : IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=16 1720 : infell(4)=4 1721 : infell(6)=4 1722 : C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16 1723 : infell(9)=16 1724 : IF (IFOUR.EQ.-3) infell(9)=19 1725 : IF (IFOUR.EQ.1) infell(9)=24 1726 : GOTO 300 1727 : C Element JOQ3 1728 : 185 infell(2)=3 1729 : infell(3)=3 1730 : infell(4)=3 1731 : infell(6)=3 1732 : infell(8)=12 1733 : infell(9)=18 1734 : infell(10)=2 1735 : infell(15)=3 1736 : infell(16)=3 1737 : GOTO 300 1738 : C Element JOQ6 1739 : 186 infell(2)=7 1740 : infell(3)=7 1741 : infell(4)=7 1742 : infell(6)=7 1743 : infell(8)=21 1744 : infell(9)=45 1745 : infell(10)=3 1746 : infell(15)=4 1747 : infell(16)=4 1748 : GOTO 300 1749 : C Element JOQ8 1750 : 187 infell(2)=9 1751 : infell(3)=9 1752 : infell(4)=9 1753 : infell(6)=9 1754 : infell(8)=28 1755 : infell(9)=60 1756 : infell(10)=3 1757 : infell(15)=4 1758 : infell(16)=4 1759 : GOTO 300 1760 : C Element JOR3 1761 : 188 infell(2)=3 1762 : infell(3)=3 1763 : infell(4)=3 1764 : infell(6)=3 1765 : infell(8)=12 1766 : infell(9)=18 1767 : infell(10)=2 1768 : infell(15)=3 1769 : infell(16)=3 1770 : GOTO 300 1771 : C Element JOR6 1772 : 189 infell(2)=7 1773 : infell(3)=7 1774 : infell(4)=7 1775 : infell(6)=7 1776 : infell(8)=21 1777 : infell(9)=45 1778 : infell(10)=3 1779 : infell(15)=4 1780 : infell(16)=4 1781 : GOTO 300 1782 : C Element JOR8 1783 : 190 infell(2)=9 1784 : infell(3)=9 1785 : infell(4)=9 1786 : infell(6)=9 1787 : infell(8)=28 1788 : infell(9)=60 1789 : infell(10)=3 1790 : infell(15)=4 1791 : infell(16)=4 1792 : GOTO 300 1793 : C Element T1D2/M1D2 : Element THERMIQUE/MECANIQUE MASSIF 1D a 2 noeuds 1794 : C On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1) 1795 : 191 infell(2)=2 1796 : infell(3)=2 1797 : infell(4)=2 1798 : infell(6)=2 1799 : infell(9)=2 1800 : IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) infell(9)=3 1801 : IF (IFOUR.EQ.11) infell(9)=4 1802 : GOTO 300 1803 : C Element T1D3/M1D3 : Element THERMIQUE/MECANIQUE MASSIF 1D a 3 noeuds 1804 : C On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1) 1805 : 192 infell(2)=3 1806 : infell(3)=3 1807 : infell(4)=3 1808 : infell(6)=3 1809 : infell(9)=3 1810 : IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) infell(9)=4 1811 : IF (IFOUR.EQ.11) infell(9)=5 1812 : GOTO 300 1813 : C Element fluide NAVIER_STOKES 1814 : C LC03 LS03 BS03 1815 : 195 infell(6)=2 1816 : GOTO 300 1817 : C LC07 LS07 BS07 1818 : 196 infell(6)=7 1819 : GOTO 300 1820 : C LC09 LS09 BS09 1821 : 197 infell(6)=4 1822 : GOTO 300 1823 : C MC03 1824 : 216 infell(6)=2 1825 : GOTO 300 1826 : C MC07 1827 : 217 infell(6)=7 1828 : GOTO 300 1829 : C MC09 1830 : 218 infell(6)=4 1831 : GOTO 300 1832 : C M103 MS03 1833 : 223 infell(6)=3 1834 : GOTO 300 1835 : C M107 MS07 1836 : 224 infell(6)=7 1837 : GOTO 300 1838 : C M109 MS09 1839 : 225 infell(6)=16 1840 : GOTO 300 1841 : C QC03 Q103 QS03 1842 : 237 infell(6)=3 1843 : GOTO 300 1844 : C QC07 Q107 QS07 1845 : 238 infell(6)=7 1846 : GOTO 300 1847 : C QC09 Q109 QS09 1848 : 239 infell(6)=16 1849 : GOTO 300 1850 : C 1851 : C CIFL 1852 : C 1853 : 258 CONTINUE 1854 : infell(2)=1 1855 : infell(3)=1 1856 : infell(4)=1 1857 : infell(6)=1 1858 : infell(8)=0 1859 : infell(9)=8 1860 : infell(10)=5 1861 : infell(15)=3 1862 : infell(16)=5 1863 : GO TO 300 1864 : C 1865 : C SHB8 1866 : C 1867 : 260 CONTINUE 1868 : infell(2)=5 1869 : infell(3)=5 1870 : infell(4)=5 1871 : infell(6)=5 1872 : infell(8)=0 1873 : infell(9)=24 1874 : infell(10)=6 1875 : infell(13)= 1 1876 : infell(15)=3 1877 : infell(16)=6 1878 : GO TO 300 1879 : C 1880 : C XQ4R 1881 : C 1882 : 263 INFELl(2)=0 1883 : INFELl(3)=64 1884 : INFELl(4)=64 1885 : INFELl(6)=64 1886 : INFELl(15)=20 1887 : INFELl(16)=4 1888 : INFELl(9) = 4*INFELl(15) 1889 : INFELl(10)=4 1890 : GOTO 300 1891 : C 1892 : C XC8R 1893 : C 1894 : 264 INFELl(2)=0 1895 : INFELl(3)=64 1896 : INFELl(4)=64 1897 : INFELl(6)=64 1898 : INFELl(15)=30 1899 : INFELl(16)=6 1900 : INFELl(9) = 8*INFELl(15) 1901 : INFELl(10)=6 1902 : GOTO 300 1903 : C 1904 : C Element JOI1 : element de liaison a 2 noeuds (6 ddl par noeuds) 1905 : C 1906 : 265 infell(2)=1 1907 : infell(3)=1 1908 : infell(4)=1 1909 : infell(5)=0 1910 : infell(6)=1 1911 : infell(7)=700 1912 : infell(8)=0 1913 : infell(13)=75 1914 : infell(14)=2 1915 : IF(IFOUR.EQ.2) THEN 1916 : infell(9)=12 1917 : infell(10)=6 1918 : infell(15)=6 1919 : infell(16)=6 1920 : ELSE IF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN 1921 : infell(9)=6 1922 : infell(10)=3 1923 : infell(15)=3 1924 : infell(16)=3 1925 : ENDIF 1926 : GOTO 300 1927 : C 1928 : C Element ZCO2 : zone cohesive a 2 noeuds (2 ddls par noeud) 1929 : C 1930 : 266 infell(2)=2 1931 : infell(4)=2 1932 : infell(6)=2 1933 : infell(7)=100 1934 : infell(9)=4 1935 : infell(10)=2 1936 : infell(15)=2 1937 : infell(16)=2 1938 : GOTO 300 1939 : C 1940 : C Element ZCO3 : zone cohesive a 3 noeuds (2 ddls par noeud) 1941 : C 1942 : 267 infell(2)=1 1943 : infell(4)=1 1944 : infell(6)=1 1945 : infell(7)=100 1946 : infell(9)=9 1947 : infell(10)=3 1948 : infell(15)=3 1949 : infell(16)=3 1950 : GOTO 300 1951 : C 1952 : C Element ZCO4 : zone cohesive a 4 noeuds (2 ddls par noeud) 1953 : C 1954 : 268 infell(2)=4 1955 : infell(4)=4 1956 : infell(6)=4 1957 : infell(7)=100 1958 : infell(9)=12 1959 : infell(10)=3 1960 : infell(15)=3 1961 : infell(16)=3 1962 : GOTO 300 1963 : C 1964 : C 1965 : 300 IF (IELE.EQ.0) IELE=infell(14) 1966 : NBGST = infell(4) 1967 : NBG = infell(6) 1968 : NBGM = infell(3) 1969 : IF (NBGM.EQ.0) THEN 1970 : NBGM=NBG 1971 : infell(3)=NBGM 1972 : ENDIF 1973 : NBSH = NBNNE(IELE) 1974 : IF (MELE.GE.79.AND.MELE.LE.83 ) NBSH = infell(8) 1975 : IF (MELE.EQ.96) NBSH = infell(8) 1976 : IF (MELE.GE.108.AND.MELE.LE.110 ) NBSH = infell(8) 1977 : IF (MELE.GE.111.AND.MELE.LE.122 ) NBSH = infell(8) 1978 : IF (MELE.GE.173.AND.MELE.LE.182 ) NBSH = infell(8) 1979 : IF (MELE.GE.185.AND.MELE.LE.190 ) NBSH = infell(8) 1980 : IF (MELE.GE.195.AND.MELE.LE.257 ) NBSH = infell(8) 1981 : IF(MELE.EQ.263.OR.MELE.EQ.264) NBSH = infell(8) 1982 : 1983 : C Cas particulier des formulations massives (mfr2 > 0) 1984 : C Ici : infell(9) = 1 * nb.noeuds.EF 1985 : IF (mfr2.NE.0) THEN 1986 : infell(9) = NBSH 1987 : ENDIF 1988 : * 1989 : * nombre composante contraintes/defo, materiau IMPEDANCE 1990 : 1991 : if (dcmate) then 1991 : call iddepl(mfr3,IFOUR,IPNOMC,NBROBL,NBRFAC) 1992 : if(tymode(/2).gt.0) then 1993 : if (tymode(1).eq.'LISTMOTS') then 1994 : dcmat2 = .true. 1995 : mlmot5 = ivamod(1) 1996 : segact mlmot5 1997 : nbrobl = mlmot5.mots(/2)*2 1998 : endif 1999 : endif 2000 : if (dcmat2.and.mele.eq.2) then 2001 : infell(9) = nbrobl 2002 : infell(10) = nbrobl/2 2003 : infell(15) = nbrobl/2 2004 : infell(16) = nbrobl/2 2005 : else 2006 : infell(15) = nbrobl 2007 : if (nbrfac.gt.0) infell(15) = infell(15) + nbrfac 2008 : infell(16) = infell(15) 2009 : if (mele.eq.45) infell(9) = nbrobl 2010 : endif 2011 : endif 2012 : 2013 : C WRITE(6,*)'TYINTE',TYINTE 2014 : C WRITE(6,*)'MELE:',MELE,'IELE:',IELE,'NBSH:',NBSH,'NBG:',NBG 2015 : 2016 : C Remplissage de infell(12) si MFR=5 cad MELE=41 ou 56 (COQ8 ou COQ6) 2017 : C ------------------------- 2018 : IF (MFR.EQ.5) THEN 2019 : * write(6,*) ' elquoi appel a renoeu ' 2019 : CALL RENOEU(IELE,MELE,NBSH,IPT2,IRT2) 2020 : * write(6,*) ' elquoi sortie renoeu ipt2 irt2' ,ipt2,irt2 2021 : if(inttyp.eq.-10) then 2022 : infmod(8)=IPT2 2023 : else 2024 : infell(12)=ipt2 2025 : endif 2026 : ENDIF 2027 : c 2028 : c cas XFEM : le 2nd segment d integration contient le MINTE de l'EF 2029 : c correspondant non enrichi (avec un nombre reduit de point de Gauss) 2030 : IF (MFR.EQ.63) THEN 2031 : if(MELE.eq.263) then 2032 : MEL2 = 8 2033 : NBG2 = 4 2034 : NBSH2 = 4 2035 : else if(MELE.eq.264) then 2036 : MEL2 = 14 2037 : NBG2 = 8 2038 : NBSH2 = 8 2039 : endif 2040 : c pour linstant on ne distingue pas les cas selon TYINTE 2040 : CALL RESHPT(NBG2,NBSH2,IELE,MEL2,0,IPT2,IRT2) 2041 : if(inttyp.eq.-10) then 2042 : infmod(8)=IPT2 2043 : else 2044 : infell(12)=IPT2 2045 : endif 2046 : ENDIF 2047 : 2048 : C Remplissage de infell(11) : INTEGRATION DE L'ELEMENT FINI 2049 : C --------------------------- 2050 : C Si le segment n'a pas deja ete rempli ou si c'est un DKT : 2051 : C == === ====== === 2052 : IDETYI=0 2053 : 400 continue 2054 : IF(INTTYP.eq.-10) then 2055 : tyinte=tyinte+1 2056 : iin=0 2057 : else 2058 : iin=INTEGR(MELE,TYINTE,IBMODE) 2059 : endif 2060 : IF ( iin.EQ.0 .OR. MELE.EQ. 28) THEN 2061 : C Cas des elements de MECANIQUE : 2062 : C ------------------------------- 2063 : C SI FAUX ? Ne devrait-on pas avoir MELE < 57 OU MELE > 68 ????? 2064 : IF (MELE.LT.57.OR.MELE.GT.64) THEN 2065 : C 1 = Champ aux noeuds 2066 : IF (TYINTE.EQ.1) THEN 2067 : IF(MELE.GE.195.and.mele.ne.260.and.mele.ne.263.and. 2068 : & mele.ne.264) THEN 2068 : CALL RESHPT(NBG,NBSH,IELE,MELE,0,IPT1,IRT1) 2069 : ELSE 2069 : CALL RENOEU(IELE,MELE,NBSH,IPT1,IRT1) 2070 : ENDIF 2071 : C 2 = Point de Gauss, centre de gravite et champ CONSTANT 2072 : ELSE IF (TYINTE.EQ.2) THEN 2072 : CALL RESHPT(1,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2073 : C 3 = Point de Gauss pour la rigidite 2074 : ELSE IF (TYINTE.EQ.3) THEN 2075 : if(mele.eq.263.or.mele.eq.264) then 2075 : CALL RESHPX(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2076 : else 2076 : CALL RESHPT(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2077 : endif 2078 : C 4 = Point de Gauss pour la masse 2079 : ELSE IF (TYINTE.EQ.4) THEN 2080 : if(mele.eq.263.or.mele.eq.264) then 2080 : CALL RESHPX(NBGM,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2081 : else 2081 : CALL RESHPT(NBGM,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2082 : endif 2083 : C 5 = Point de Gauss - calcul des contraintes 2084 : ELSE IF (TYINTE.EQ.5) THEN 2085 : if(mele.eq.263.or.mele.eq.264) then 2085 : CALL RESHPX(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2086 : else 2086 : CALL RESHPT(NBGST,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2087 : endif 2088 : C 6 = Pas de segment d'integration 2089 : ELSE IF (TYINTE.EQ.6) THEN 2090 : IPT1 = 0 2091 : C 7 = Champ aux faces 2092 : C 8 = Champ aux CENTREP1 2093 : ELSE IF (TYINTE.EQ.8) THEN 2093 : CALL RESHPT(NBG,NBSH,IELE,MELE,1,IPT1,IRT1) 2094 : C 9 = Champ aux MSOMMET 2095 : ELSE IF (TYINTE.EQ.9) THEN 2095 : CALL RESHPT(NBG,NBSH,IELE,MELE,2,IPT1,IRT1) 2096 : C Autre cas = Point de Gauss pour la rigidite 2097 : ELSE 2097 : CALL RESHPT(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1) 2098 : ENDIF 2099 : C Cas des elements de THERMIQUE : 2100 : C ------------------------------- 2101 : C ELSE 2101 : C CALL TINTR(MELE,IPT1) 2102 : ENDIF 2103 : if(inttyp.eq.-10) then 2104 : infmod(2+tyinte+idetyi)=ipt1 2105 : if( tyinte.lt.5) goto 400 2106 : if(tyinte.eq.5.and.idetyi.eq.0) then 2107 : if(mele.eq.263) then 2108 : * pour le XQ4R on recommence pour avoir les fonctions de formes 2109 : classique du : qua4 2110 : tyinte=0 2111 : nbsh=4 2112 : nbg=4 2113 : nbgm=4 2114 : nbgst=4 2115 : idetyi=5 2116 : mele=8 2117 : go to 400 2118 : else if(mele.eq.264) then 2119 : * pour le XC8R on recommence pour avoir les fonctions de formes 2120 : classique du : cub8 2121 : tyinte=0 2122 : nbsh=8 2123 : nbg=8 2124 : nbgm=8 2125 : nbgst=8 2126 : idetyi=5 2127 : mele=14 2128 : go to 400 2129 : endif 2130 : endif 2131 : else 2132 : INTEGR(MELE,TYINTE,IBMODE) = IPT1 2133 : infell(11) = IPT1 2133 : IF(IPT1.NE.0)CALL SAVSEG(IPT1) 2134 : endif 2135 : C Si le segment d'integration a deja ete rempli : 2136 : C == ==== ====== 2137 : ELSE 2138 : infell(11) = INTEGR(MELE,TYINTE,IBMODE) 2139 : ENDIF 2140 : 301 CONTINUE 2141 : IPTR=INFO 2142 : c write(6,*)'-> ELQUOI infell',(infell(iou),iou=1,16) 2143 : RETURN 2144 : END 2145 : 2146 : 2147 : 2148 : 2149 : 2150 : 2151 :
© Cast3M 2003 - All rights reserved.
Disclaimer