C PRQUOI SOURCE OF166741 26/02/23 21:15:26 12480 C======================================================================= C= Sous-programme PRQUOI = C= --------------------- = C= = C= Ce sous programme fournit les caracteristiques des elements finis. = C= = C= Entrees : = C= --------- = C= IMODEL : pointeur sur un IMODEL, zone elementaire d'un MMODEL = C= = C= Sorties : = C= --------- = C= L'element INFELE du segment IMODEL a ete renseigne et contient : = C= INFELE( 1) : numero de l'element fini = C= INFELE( 2) : nb de points d'integration en contrainte multicouche = C= INFELE( 3) : nb de points d'integration pour MASSE = C= INFELE( 4) : nb de points d'integration pour SIGMA BSIGMA KSIGMA = C= INFELE( 5) : nb de caracteristiques = C= INFELE( 6) : nb de points d'integration pour RIGI/COND = C= INFELE( 7) : longueur d'un tableau de travail pour l'element = C= INFELE( 8) : nombre de fonctions de forme = C= INFELE( 9) : nb de d.d.l. dans la matrice de RIGIDITE = C= INFELE(10) : taille de la matrice de Hooke = C= INFELE(11) : pointeur sur le segment d'integration = C= INFELE(12) : inutilise = C= INFELE(13) : numero de la formulation de l'element fini (NUMMFR) = C= INFELE(14) : numero de l'element geometrique associe (NUMGEO) = C= INFELE(15) : nb maximal de d.d.l. par noeud = C= INFELE(16) : nb de composantes de contraintes ou de deformations = C======================================================================= SUBROUTINE PRQUOI(IMODEL) C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) C -INC PPARAM -INC CCOPTIO -INC CCGEOME -INC CCHAMP C= DEB === Formulation HHO ===== Include specifique ==================== -INC CCHHOPA C= FIN === Formulation HHO ============================================= -INC SMMODEL -INC SMELEME -INC SMLMOTS C logical dcmate,dcmat2,bmodi C IF (IMODEL.INFMOD(/1).LT.1) THEN WRITE(IOIMP,*) WRITE(IOIMP,*) 'PRQUOI : INFMOD(/1) incorrect',INFMOD(/1),' < 1' CALL ERREUR(5) ENDIF C C (0) petite verrue pour la formulation CONTACT IF (IMODEL.FORMOD(1) .EQ. 'CONTACT ') RETURN C Recuperation des donnees de INFELE en entree MELE = IMODEL.NEFMOD IELE = NUMGEO(MELE) MFR = NUMMFR(MELE) C C (1) petite verrue pour les MULT IF (MELE .EQ. 22) RETURN C (2) petite verrue pour les SURE C Travail fait par ailleurs : IF (MELE .EQ. 259) RETURN C==== Formulation HHO ===== Appel specifique a HHOELQ ================== IF (MELE .EQ. HHO_NUM_ELEMENT) THEN C test equivalent IF (MFR .EQ. HHO_MFR_ELEMENT) THEN CALL HHOELQ(IMODEL) RETURN END IF C==== Formulation HHO ================================================== dcmate = .false. dcmat2 = .false. bmodi=.false. do im = 1, imodel.matmod(/2) if (imodel.matmod(im).eq.'IMPEDANCE') then dcmate = .true. meleme = imodel.imamod ityp1 = meleme.itypel mfr = nummfr(mele) if (ityp1.eq.1) then bmodi=.true. mele = 45 melep=imodel.nefmod imodel.nefmod=mele endif iele = numgeo(mele) endif enddo C C Nb de points d'integration dans l'epaisseur NPINT = IMODEL.INFMOD(1) C C ================================================================= C INFELE(1) : Numero de l'element fini C ================================================================= INFELE(1) = MELE C ================================================================= C INFELE(8) : Nombre de fonctions de forme de l'element C ================================================================= INFELE(8) = NBNNE(IELE) C ================================================================= C INFELE(13) : Numero de la formulation de l'element fini C ================================================================= INFELE(13) = MFR IF (FORMOD(1).EQ.'LIAISON'.AND.MFR.EQ.27) INFELE(13) = 24 DO IM = 1,MATMOD(/2) IF (MFR.EQ.27.AND.MATMOD(IM).EQ.'MODAL') INFELE(13) = 26 IF (MFR.EQ.27.AND.MATMOD(IM).EQ.'STATIQUE') INFELE(13) = 28 ENDDO IF (FORMOD(1).EQ.'CHARGEMENT') THEN IF (MFR.EQ.1) INFELE(13) = 72 IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.7.OR.MFR.EQ.9.OR. & MFR.EQ.13) INFELE(13) = 74 ENDIF MFR = INFELE(13) IF (MELE.GE.99.AND.MELE.LE.103) INFELE(13)=1 C ELECTROSTATIQUE : MFR=71=INFELE(13) et NUMMFR(mele)=1 (MASSIF) IF (FORMOD(1).EQ.'ELECTROSTATIQUE') INFELE(13)=71 C DIFFUSION : MFR=73=INFELE(13) et NUMMFR(mele)=1 (MASSIF) IF (FORMOD(1).EQ.'DIFFUSION') THEN IF (MFR.EQ.1) INFELE(13)=73 ENDIF C ================================================================= C INFELE(14) : Numero de l'element geometrique C ================================================================= INFELE(14) = NUMGEO(MELE) C Redirection suivant l'element fini C ---------------------------------------------------------- IF(MELE.GE.1.AND.MELE.LE.100) THEN C CABL SEG2 SEG3 TRI3 TRI4 TRI6 TRI7 QUA4 QUA5 QUA8 GOTO ( 300, 2, 3, 4, 300, 6, 7, 8, 300, 10 C QUA9 RAC2 RAC3 CUB8 CU20 PRI6 PR15 LIA3 LIA4 LIA6 1 , 11, 12, 3, 14, 15, 16, 17, 18, 19, 6 C LIA8 MULT TET4 TE10 PYR5 PY13 COQ3 DKT POUT LISP 2 , 10, 300, 23, 24, 25, 26, 27, 28, 29, 30 C FAC3 FAC4 FAC6 FAC8 LTR3 LQU4 LCU8 LPR6 LTE4 LPY5 3 , 31, 32, 33, 34, 35, 36, 37, 38, 39, 40 C COQ8 TUYA TUFI COQ2 POI1 BARR RACO LSU2 COQ4 LISM 4 , 41, 42, 43, 44, 45, 46, 47, 48, 49, 50 C COF3 RES2 LSU3 LSU4 LICO COQ6 CVS2 CVS3 CVT3 CVT6 5 , 51, 300, 53, 54, 55, 56, 57, 58, 59, 60 C CVQ4 CVQ8 THP5 TH13 THP6 TH15 THC8 TH20 ICT3 ICQ4 6 , 61, 62, 63, 64, 65, 66, 67, 68, 69, 70 C ICT6 ICQ8 ICC8 ICT4 ICP6 IC20 IC10 IC15 TRIP QUAP 7 , 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 C CUBP TETP PRIP TIMO JOI2 JOI3 JOT3 JOI4 JOI6 JOI8 8 , 81, 82, 83, 84, 85, 86, 87, 88, 89, 90 C LISC TRIH DST LIC4 CERC TUYO LSE2 LITU HYT3 HYQ4 9 , 91, 92, 93, 94, 95, 96, 97, 98, 4, 8) C & ,MELE ELSEIF(MELE.GE.101.AND.MELE.LE.200) THEN C HYT4 HYP6 HYC8 TRIS QUAS POIS FOR3 JOP3 JOP6 JOP8 GOTO ( 23, 16, 14, 104, 105, 106, 300, 108, 109, 110 C POL3 POL4 POL5 POL6 POL7 POL8 POL9 PO10 PO11 PO12 1 , 111, 111, 111, 111, 111, 111, 111, 111, 111, 111 C PO13 PO14 BAR3 BAEX LIA2 QUAH CUBH ROT3 SEF2 TRF3 2 , 111, 111, 123, 124, 125, 126, 127, 128, 2, 4 C QUF4 CUF8 PRF6 TEF4 PYF5 MSE3 MTR6 MQU9 MC27 MP18 3 , 8, 14, 16, 23, 25, 3, 6, 11, 275, 300 C MT10 MP14 SEF3 TRF7 QUF9 CF27 PF21 TF15 PF19 SEG6 4 , 24, 300, 3, 7, 11, 275, 276, 277, 278, 300 C TR21 QU36 C216 P126 TE56 PY91 TRH6 BSE2 BTR4 BQU5 5 , 300, 300, 300, 300, 300, 300, 157, 300, 300, 300 C BCU9 BPR7 BTE5 BPY6 FRO4 SEGS POJS JCT3 JCI4 JGI2 6 , 300, 300, 300, 300, 300, 166, 167, 168, 169, 170 C JGT3 JGI4 TRIQ QUAQ CUBQ TETQ PRIQ TRIR QUAR CUBR 7 , 171, 172, 173, 174, 175, 176, 177, 178, 179, 180 C TETR PRIR Q4RI Q8RI JOQ3 JOQ6 JOQ8 JOR3 JOR6 JOR8 8 , 181, 182, 183, 184, 185, 186, 187, 188, 189, 190 C T1D2 T1D3 M1D2 M1D3 LC03 LC07 LC09 LC27 LC21 LC15 9 , 191, 192, 191, 192, 195, 196, 197, 300, 300, 300) C & ,MELE-100 ELSEIF(MELE.GE.201.AND.MELE.LE.300) THEN C LC19 LS03 LS07 LS09 LS27 LS21 LS15 LS19 BS03 BS07 GOTO ( 300, 195, 196, 197, 300, 300, 300, 300, 195, 196 C BS09 BS27 BS21 BS15 BS19 MC03 MC07 MC09 MC27 MC21 1 , 197, 300, 300, 300, 300, 216, 217, 218, 300, 300 C MC15 MC19 M103 M107 M109 M127 M121 M115 M119 MS03 2 , 300, 300, 223, 224, 225, 300, 300, 300, 300, 223 C MS07 MS09 MS27 MS21 MS15 MS19 QC03 QC07 QC09 QC27 3 , 224, 225, 300, 300, 300, 300, 237, 238, 239, 300 C QC21 QC15 QC19 Q103 Q107 Q109 Q127 Q121 Q115 Q119 4 , 300, 300, 300, 237, 238, 239, 300, 300, 300, 300 C QS03 QS07 QS09 QS27 QS21 QS15 QS19 CIFL SURE SHB8 5 , 237, 238, 239, 300, 300, 300, 300, 300, 300, 260 C CAF2 CAF3 XQ4R XC8R JOI1 ZCO2 ZCO3 ZCO4 TUY2 TUY3 6 , 300, 300, 263, 264, 265, 266, 267, 268, 269, 270 C COS2 COA2 ICY5 IC13 CU27 PR21 TE15 PY19 C20R P15R 7 , 271, 272, 25, 26, 275, 276, 277, 278, 279, 280 C HHO .... .... .... .... .... .... .... .... .... 8 , 999, 999, 999, 999, 999, 999, 999, 999, 999, 999 C .... .... .... .... .... .... .... .... .... .... 9 , 999, 999, 999, 999, 999, 999, 999, 999, 999, 999) C & ,MELE-200 ENDIF C C ================================================================= C ERREUR 999 CONTINUE CALL ERREUR(5) C C ================================================================= C SEG2 C ================================================================= 2 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(6)=2 INFELE(7)=0 INFELE(9)=4 INFELE(10)=1 INFELE(15)=2 INFELE(16)=1 IF (IFOUR.EQ.1) THEN INFELE(9)=6 INFELE(15)=3 ENDIF GOTO 300 C ================================================================= C SEG3 C ================================================================= 3 CONTINUE INFELE(2)=3 INFELE(4)=3 INFELE(6)=3 INFELE(7)=0 INFELE(9)=6 INFELE(10)=1 INFELE(15)=2 INFELE(16)=1 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 IF (IFOUR.EQ.1) THEN INFELE(9)=9 INFELE(15)=3 ENDIF ENDIF IF (MELE.EQ.99) THEN INFELE(3)=4 INFELE(6)=4 ENDIF GOTO 300 C ================================================================= C TRI3 C ================================================================= 4 CONTINUE INFELE(2)=1 INFELE(3)=4 INFELE(4)=1 INFELE(6)=1 INFELE(9)=6 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=7 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=9 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C TRI6 C ================================================================= 6 CONTINUE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=7 INFELE(4)=7 INFELE(6)=7 ELSE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 ENDIF INFELE(3)=7 INFELE(9)=12 IF (IFOUR.EQ.-3) INFELE(9)=15 IF (IFOUR.EQ.1) INFELE(9)=18 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C TRI7 C ================================================================= 7 CONTINUE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=7 INFELE(4)=7 INFELE(6)=7 ELSE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 ENDIF INFELE(3)=7 INFELE(9)=14 IF (IFOUR.EQ.-3) INFELE(9)=17 IF (IFOUR.EQ.1) INFELE(9)=21 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C QUA4 C ================================================================= 8 CONTINUE INFELE(2)=4 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=4 C** INFELE(4)=5 CORRECTION MILL LE 8/8/90 INFELE(6)=4 C** IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9 INFELE(9)=8 IF (IFOUR.EQ.-3) INFELE(9)=11 IF (IFOUR.EQ.1) INFELE(9)=12 IF (MELE.EQ.100) THEN INFELE(3)=4 INFELE(6)=4 ENDIF IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C QUA8 C ================================================================= 10 CONTINUE INFELE(2)=9 INFELE(3)=9 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16 INFELE(4)=9 INFELE(6)=9 C IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(9)=16 IF (IFOUR.EQ.-3) INFELE(9)=19 IF (IFOUR.EQ.1) INFELE(9)=24 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C QUA9 C ================================================================= 11 CONTINUE INFELE(2)=9 INFELE(3)=9 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16 INFELE(4)=9 INFELE(6)=9 C IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(9)=18 IF (IFOUR.EQ.-3) INFELE(9)=21 IF (IFOUR.EQ.1) INFELE(9)=27 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C RAC2 (raccord a 4 noeuds en 2 dimensions) C ================================================================= 12 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(5)=2 INFELE(6)=2 INFELE(9)=8 IF (IFOUR.EQ.1) INFELE(9)=10 INFELE(10)=6 INFELE(15)=3 GOTO 300 C ================================================================= C CUB8 C ================================================================= 14 CONTINUE INFELE(2)=8 INFELE(3)=8 INFELE(4)=8 INFELE(6)=8 INFELE(9)=24 IF (MELE.EQ.103) THEN INFELE(3)=8 INFELE(6)=8 ENDIF INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C CU20 C ================================================================= 15 CONTINUE INFELE(2)=27 INFELE(3)=27 INFELE(4)=27 INFELE(6)=27 INFELE(9)=60 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C PRI6 C ================================================================= 16 CONTINUE INFELE(2)=6 INFELE(3)=8 INFELE(4)=6 INFELE(6)=6 INFELE(9)=18 IF (MELE.EQ.102) THEN INFELE(3)=6 INFELE(6)=6 ENDIF INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C PR15 C ================================================================= 17 CONTINUE INFELE(2)=12 INFELE(3)=21 INFELE(4)=12 INFELE(6)=12 INFELE(9)=45 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C LIA3 (raccord a 6 noeudes en 3 dimensions) C ================================================================= 18 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(5)=3 INFELE(6)=4 INFELE(9)=15 INFELE(10)=6 INFELE(15)=3 GOTO 300 C ================================================================= C LIA4 (raccord a 8 noeuds en 3 dimensions) C ================================================================= 19 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(5)=3 INFELE(6)=4 INFELE(9)=20 INFELE(10)=6 INFELE(15)=3 GOTO 300 C ================================================================= C TET4 C ================================================================= 23 CONTINUE INFELE(2)=1 INFELE(3)=4 INFELE(4)=1 INFELE(6)=1 INFELE(9)=12 IF (MELE.EQ.101) THEN INFELE(3)=4 INFELE(6)=4 ENDIF INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C TE10 C ================================================================= 24 CONTINUE INFELE(2)=5 cbp INFELE(3)=5 INFELE(3)=15 INFELE(4)=5 INFELE(6)=5 INFELE(9)=30 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C PYR5 C ================================================================= 25 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(9)=15 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C PY13 C ================================================================= 26 CONTINUE INFELE(2)=10 INFELE(3)=17 INFELE(4)=10 INFELE(6)=10 INFELE(9)=39 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 IF (MELE.EQ.274) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ENDIF GOTO 300 C ================================================================= C COQ3 C ================================================================= 27 CONTINUE INFELE(2)=3 INFELE(4)=1 INFELE(5)=3 INFELE(6)=1 INFELE(7)=600 INFELE(9)=18 IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=9 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C DKT C ================================================================= 28 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(5)=3 INFELE(7)=600 INFELE(9)=18 INFELE(15)=6 IF (NPINT.EQ.0) THEN C Element DKT COQUE GLOBALE INFELE(4)=3 INFELE(6)=3 INFELE(10)=6 INFELE(16)=6 ELSE C Element DKT COQUE INTEGREE dans l'epaisseur INFELE(4)=3*NPINT INFELE(6)=3*NPINT INFELE(10)=6 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C POUTRE C ================================================================= 29 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(6)=2 INFELE(7)=200 IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN C Element POUTRE 2D plan INFELE(5)=6 INFELE(9)=6 INFELE(10)=3 INFELE(15)=3 INFELE(16)=3 ELSE C Element POUTRE 3D INFELE(5)=12 INFELE(9)=12 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 ENDIF GOTO 300 C ================================================================= C LINE SPRING LISP C ================================================================= 30 CONTINUE INFELE(2)=3 INFELE(4)=3 INFELE(5)=5 INFELE(6)=3 INFELE(7)=100 INFELE(9)=24 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C FAC3 pour les PRESSIONS uniquement C ================================================================= 31 CONTINUE INFELE(2)=1 INFELE(4)=1 INFELE(6)=1 INFELE(9)=9 INFELE(15)=3 GOTO 300 C ================================================================= C FAC4 pour les PRESSIONS uniquement C ================================================================= 32 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 INFELE(9)=12 INFELE(15)=3 GOTO 300 C ================================================================= C FAC6 pour les PRESSIONS uniquement C ================================================================= 33 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 INFELE(9)=18 INFELE(15)=3 GOTO 300 C ================================================================= C FAC8 pour les PRESSIONS uniquement C ================================================================= 34 CONTINUE INFELE(2)=9 INFELE(4)=9 INFELE(6)=9 INFELE(9)=24 INFELE(15)=3 GOTO 300 C ================================================================= C LTR3 : element LIQUIDE triangle a 3 noeuds C ================================================================= 35 CONTINUE INFELE(2)=4 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=7 INFELE(4)=4 INFELE(6)=4 INFELE(9)=6 INFELE(10)=5 INFELE(15)=2 INFELE(16)=0 GOTO 300 C ================================================================= C LQU4 : element LIQUIDE quadrangle a 4 noeuds C ================================================================= 36 CONTINUE INFELE(2)=4 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=4 INFELE(6)=4 INFELE(9)=8 INFELE(10)=5 INFELE(15)=2 INFELE(16)=0 GOTO 300 C ================================================================= C LCU8 : Element LIQUIDE CUBE a 8 noeuds C ================================================================= 37 CONTINUE INFELE(2)=8 INFELE(3)=8 INFELE(4)=8 INFELE(6)=8 INFELE(9)=16 INFELE(10)=5 INFELE(15)=2 INFELE(16)=0 GOTO 300 C ================================================================= C LPR6 : Element LIQUIDE PRISME a 6 noeuds C ================================================================= 38 CONTINUE INFELE(2)=6 INFELE(3)=6 INFELE(4)=6 INFELE(6)=6 INFELE(9)=12 INFELE(10)=5 INFELE(15)=2 INFELE(16)=0 GOTO 300 C ================================================================= C LTE4 : Element LIQUIDE TETRAEDRE a 4 noeuds C ================================================================= 39 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(6)=4 INFELE(9)=8 INFELE(10)=5 INFELE(15)=2 INFELE(16)=0 GOTO 300 C ================================================================= C LPY5 : Element LIQUIDE PYRAMIDE a 5 noeuds C ================================================================= 40 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(9)=10 INFELE(10)=5 INFELE(15)=2 INFELE(16)=0 GOTO 300 C ================================================================= C COQ8 : Element COQUE EPAISSE a 8 noeuds C ================================================================= 41 CONTINUE INFELE(2)=NPINT IF (NPINT.EQ.0) INFELE(2)=8 INFELE(3)=8 INFELE(4)=NPINT IF (NPINT.EQ.0) INFELE(4)=8 INFELE(5)=2 INFELE(6)=8 INFELE(7)=1700 INFELE(9)=48 IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=24 INFELE(10)=5 INFELE(15)=6 INFELE(16)=5 GOTO 300 C ================================================================= C TUYA : Element TUYAU C ================================================================= 42 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(5)=7 INFELE(6)=2 INFELE(7)=200 INFELE(9)=12 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C TUFI : Element TUYAU FISSURE C ================================================================= 43 CONTINUE INFELE(2)=1 INFELE(3)=0 INFELE(4)=1 INFELE(5)=9 INFELE(6)=1 INFELE(7)=200 INFELE(9)=12 INFELE(10)=6 INFELE(15)=6 INFELE(16)=8 GOTO 300 C ================================================================= C COQ2 C ================================================================= 44 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(5)=3 INFELE(6)=3 INFELE(7)=368 INFELE(9)=8 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(9)=6 IF (IFOUR.EQ.-3) INFELE(9)=9 INFELE(10)=6 IF (IFOUR.LE.0) INFELE(10)=4 INFELE(15)=4 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(15)=3 INFELE(16)=6 IF (IFOUR.LE.0) INFELE(16)=4 GOTO 300 C ================================================================= C POI1 (MODEs en DEFOrmations PLANes GENEralisees) C ================================================================= 45 CONTINUE C Valeurs a adapter pour certains MODE UNID (1D) INFELE(2)=1 INFELE(3)=1 INFELE(4)=1 INFELE(5)=1 INFELE(6)=1 INFELE(7)=50 INFELE(15)=2 IF (IFOUR.EQ.11) THEN INFELE(9)=3 INFELE(10)=2 INFELE(16)=2 ELSE INFELE(9)=5 IF (IFOUR.GE.3.AND.IFOUR.LE.15) INFELE(9)=2 INFELE(10)=1 INFELE(16)=1 ENDIF if (mfr.eq.26.or.mfr.eq.28) INFELE(9) = 1 GOTO 300 C ================================================================= C BARRE C ================================================================= 46 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(5)=1 INFELE(6)=2 INFELE(7)=200 INFELE(9)=6 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(9)=4 IF (IFOUR.EQ.-3) INFELE(9)=7 IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=2 INFELE(10)=1 INFELE(15)=3 IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(15)=1 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) INFELE(15)=2 INFELE(16)=1 GOTO 300 C ================================================================= C RACO : Element de RACCORD LIQUIDE-COQUE a 4 noeuds en 2D C ================================================================= 47 CONTINUE INFELE(2)=1 INFELE(4)=1 INFELE(5)=2 INFELE(6)=1 INFELE(9)=10 IF (IFOUR.EQ.1) INFELE(9)=12 INFELE(15)=4 INFELE(10)=6 GOTO 300 C ================================================================= C LSU2 : Element de SURFACE LIBRE a 2 noeuds en 2D C ================================================================= 48 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(6)=2 INFELE(9)=6 INFELE(10)=6 INFELE(15)=3 GOTO 300 C ================================================================= C COQ4 C ================================================================= 49 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(5)=3 INFELE(6)=5 INFELE(7)=100 INFELE(9)=24 IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=12 INFELE(10)=8 INFELE(15)=6 INFELE(16)=8 GOTO 300 C ================================================================= C LINE SPRING LISM C ================================================================= 50 CONTINUE INFELE(2)=1 INFELE(4)=1 INFELE(5)=5 INFELE(6)=3 INFELE(7)=100 INFELE(9)=24 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C COF3 C ================================================================= 51 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(5)=3 INFELE(6)=3 INFELE(7)=368 INFELE(9)=12 INFELE(10)=6 INFELE(15)=4 INFELE(16)=6 GOTO 300 C ================================================================= C LSU3 : Element de SURFACE LIBRE a 3 noeuds en 3D C ================================================================= 53 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 INFELE(9)=9 INFELE(10)=6 INFELE(15)=3 GOTO 300 C ================================================================= C LSU4 : Element de SURFACE LIBRE a 4 noeuds en 3D C ================================================================= 54 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 INFELE(9)=12 INFELE(10)=6 INFELE(15)=3 GOTO 300 C ================================================================= C LICO : Element de RACCORD LIQUIDE-COQUE a 6 noeuds en 3D C ================================================================= 55 CONTINUE INFELE(2)=3 INFELE(4)=3 INFELE(5)=3 INFELE(6)=3 INFELE(9)=24 INFELE(10)=6 INFELE(15)=6 GOTO 300 C ================================================================= C COQ6 : Element COQUE EPAISSE a 6 noeuds C ================================================================= 56 CONTINUE INFELE(2)=NPINT IF(NPINT.EQ.0) INFELE(2)=6 INFELE(3)=6 INFELE(4)=NPINT IF (NPINT.EQ.0) INFELE(4)=6 INFELE(5)=2 INFELE(6)=6 INFELE(7)=1000 INFELE(9)=36 IF (FORMOD(1) .EQ. 'DIFFUSION') INFELE(9)=18 INFELE(10)=5 INFELE(15)=6 INFELE(16)=5 GOTO 300 C ================================================================= C THT3 (thermique) C ================================================================= 57 CONTINUE INFELE(2)=1 INFELE(6)=1 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=4 INFELE(9)=3 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=9 INFELE(10)=2 INFELE(15)=1 INFELE(16)=2 GOTO 300 C ================================================================= C THT6 C ================================================================= 58 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=7 INFELE(9)=6 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=18 INFELE(10)=2 INFELE(15)=1 INFELE(16)=2 GOTO 300 C ================================================================= C THQ4 C ================================================================= 59 CONTINUE INFELE(2)=4 INFELE(6)=4 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9 INFELE(9)=4 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=12 INFELE(10)=2 INFELE(15)=1 INFELE(16)=2 GOTO 300 C ================================================================= C THQ8 C ================================================================= 60 CONTINUE INFELE(2)=9 INFELE(6)=9 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(9)=8 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=24 INFELE(10)=2 INFELE(15)=1 INFELE(16)=2 GOTO 300 C ================================================================= C THT4 C ================================================================= 61 CONTINUE INFELE(2)=1 INFELE(6)=1 INFELE(9)=4 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C TH10 C ================================================================= 62 CONTINUE INFELE(2)=5 INFELE(6)=5 INFELE(9)=10 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C THP5 C ================================================================= 63 CONTINUE INFELE(2)=5 INFELE(6)=5 INFELE(9)=5 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C TH13 C ================================================================= 64 CONTINUE INFELE(2)=27 INFELE(6)=27 INFELE(9)=13 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C THP6 C ================================================================= 65 CONTINUE INFELE(2)=2 INFELE(6)=2 INFELE(9)=6 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C TH15 C ================================================================= 66 CONTINUE INFELE(2)=12 INFELE(6)=12 INFELE(9)=15 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C THC8 C ================================================================= 67 CONTINUE INFELE(2)=8 INFELE(6)=8 INFELE(9)=8 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C TH20 C ================================================================= 68 CONTINUE INFELE(2)=27 INFELE(6)=27 INFELE(9)=20 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C ICT3 C ================================================================= 69 CONTINUE INFELE(2)=1 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=7 INFELE(4)=1 INFELE(6)=1 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=4 INFELE(9)=6 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=9 INFELE(10)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4 INFELE(15)=2 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3 INFELE(16)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4 GOTO 300 C ================================================================= C ICQ4 C ================================================================= 70 CONTINUE INFELE(2)=4 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=4 C** INFELE(4)=5 CORRECTION MILL LE 8/8/90 INFELE(6)=4 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9 INFELE(9)=8 IF (IFOUR.EQ.-3) INFELE(9)=11 IF (IFOUR.EQ.1) INFELE(9)=12 INFELE(10)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4 INFELE(15)=2 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3 INFELE(16)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4 GOTO 300 C ================================================================= C ICT6 C ================================================================= 71 CONTINUE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=7 INFELE(4)=7 INFELE(6)=7 ELSE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 ENDIF INFELE(3)=7 INFELE(9)=12 IF (IFOUR.EQ.-3) INFELE(9)=15 IF (IFOUR.EQ.1) INFELE(9)=18 INFELE(10)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4 INFELE(15)=2 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3 INFELE(16)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4 GOTO 300 C ================================================================= C ICQ8 C ================================================================= 72 CONTINUE INFELE(2)=9 INFELE(3)=9 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16 INFELE(4)=9 INFELE(6)=9 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(9)=16 IF (IFOUR.EQ.-3) INFELE(9)=19 IF (IFOUR.EQ.1) INFELE(9)=24 INFELE(10)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(10)=4 INFELE(15)=2 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(15)=3 INFELE(16)=4 C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(16)=4 GOTO 300 C ================================================================= C ICC8 C ================================================================= 73 CONTINUE INFELE(2)=8 INFELE(3)=8 INFELE(4)=8 INFELE(6)=8 INFELE(9)=24 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C ICT4 C ================================================================= 74 CONTINUE INFELE(2)=1 INFELE(3)=4 INFELE(4)=1 INFELE(6)=1 INFELE(9)=12 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C ICP6 C ================================================================= 75 CONTINUE INFELE(2)=6 INFELE(3)=8 INFELE(4)=6 INFELE(6)=6 INFELE(9)=18 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C IC20 C ================================================================= 76 CONTINUE INFELE(2)=27 INFELE(3)=27 INFELE(4)=27 INFELE(6)=27 INFELE(9)=60 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C IC10 C ================================================================= 77 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(9)=30 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C IC15 C ================================================================= 78 CONTINUE INFELE(2)=12 INFELE(3)=21 INFELE(4)=12 INFELE(6)=12 INFELE(9)=45 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C TRIP C ================================================================= 79 CONTINUE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=7 INFELE(4)=7 INFELE(6)=7 ELSE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 ENDIF INFELE(3)=7 INFELE(8)=9 INFELE(9)=15 IF (IFOUR.EQ.-3.OR.IFOUR.EQ.1) INFELE(9)=21 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN INFELE(10)=6 INFELE(15)=4 INFELE(16)=7 ELSE INFELE(10)=4 INFELE(15)=3 INFELE(16)=5 ENDIF GOTO 300 C ================================================================= C QUAP C ================================================================= 80 CONTINUE INFELE(2)=9 INFELE(3)=9 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=9 INFELE(6)=9 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(8)=12 INFELE(9)=20 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=28 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN INFELE(10)=6 INFELE(15)=4 INFELE(16)=7 ELSE INFELE(10)=4 INFELE(15)=3 INFELE(16)=5 ENDIF GOTO 300 C ================================================================= C CUBP C ================================================================= 81 CONTINUE INFELE(2)=27 INFELE(3)=27 INFELE(4)=27 INFELE(6)=27 INFELE(8)=28 INFELE(9)=68 INFELE(10)=6 INFELE(15)=4 INFELE(16)=7 GOTO 300 C ================================================================= C TETP C ================================================================= 82 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(8)=14 INFELE(9)=34 INFELE(10)=6 INFELE(15)=4 INFELE(16)=7 GOTO 300 C ================================================================= C PRIP C ================================================================= 83 CONTINUE INFELE(2)=12 INFELE(3)=12 INFELE(4)=12 INFELE(6)=12 INFELE(8)=21 INFELE(9)=51 INFELE(10)=6 INFELE(15)=4 INFELE(16)=7 GOTO 300 C ================================================================= C TIMO : Element POUTRE de TIMOSHENCKO C ================================================================= 84 CONTINUE INFELE(2)=1 INFELE(4)=1 INFELE(6)=1 INFELE(7)=200 IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN C Element TIMO 2D plan INFELE(5)=6 INFELE(9)=6 INFELE(10)=3 INFELE(15)=3 INFELE(16)=3 ELSE C Element TIMO 3D INFELE(5)=12 INFELE(9)=12 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 ENDIF GOTO 300 C ================================================================= C JOI2 C ================================================================= 85 CONTINUE INFELE(2)=4 INFELE(3)=2 INFELE(4)=2 INFELE(5)=0 INFELE(6)=2 INFELE(7)=300 INFELE(9)=8 INFELE(10)=2 INFELE(16)=2 GOTO 300 C ================================================================= C JOI3 C ================================================================= 86 CONTINUE INFELE(2)=3 INFELE(3)=2 INFELE(4)=3 INFELE(5)=0 INFELE(6)=3 INFELE(7)=100 INFELE(9)=12 INFELE(10)=2 INFELE(16)=2 GOTO 300 C ================================================================= C JOT3 C ================================================================= 87 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(5)=0 INFELE(6)=3 INFELE(7)=100 INFELE(9)=18 INFELE(10)=3 INFELE(16)=3 GOTO 300 C ================================================================= C JOI4 C ================================================================= 88 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(5)=0 INFELE(6)=4 INFELE(7)=100 INFELE(9)=24 INFELE(10)=3 INFELE(16)=3 GOTO 300 C ================================================================= C JOI6 C ================================================================= 89 CONTINUE INFELE(2)=5 INFELE(4)=5 INFELE(6)=4 INFELE(9)=36 INFELE(10)=3 INFELE(16)=3 GOTO 300 C ================================================================= C JOI8 C ================================================================= 90 CONTINUE INFELE(2)=9 INFELE(4)=9 INFELE(6)=9 INFELE(9)=48 INFELE(10)=3 INFELE(16)=3 GOTO 300 C ================================================================= C LINE SPRING LISC C ================================================================= 91 CONTINUE INFELE(2)=3 INFELE(4)=3 INFELE(5)=5 INFELE(6)=5 INFELE(7)=100 INFELE(9)=36 INFELE(10)=2 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C TRIH C ================================================================= 92 CONTINUE INFELE(2)=4 INFELE(6)=4 INFELE(7)= 0 IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN INFELE(3)=7 INFELE(5)=5 IF (IFOUR.EQ.0) INFELE(9)=12 IF (IFOUR.EQ.1) INFELE(9)=18 INFELE(15)=6 ELSE INFELE(3)=4 INFELE(5)=5 INFELE(9)=12 INFELE(15)=4 ENDIF INFELE(10)=10 INFELE(16)=0 GOTO 300 C ================================================================= C DST C ================================================================= 93 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(5)=3 INFELE(6)=3 INFELE(7)=600 INFELE(9)=18 INFELE(10)=8 INFELE(15)=6 INFELE(16)=8 GOTO 300 C ================================================================= C LIC4 : Element de RACCORD LIQUIDE-COQUE a 8 noeuds en 3D C ================================================================= 94 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(5)=3 INFELE(6)=4 INFELE(9)=32 INFELE(10)=6 INFELE(15)=6 GOTO 300 C ================================================================= C CERCE C ================================================================= 95 CONTINUE INFELE(2)=1 INFELE(3)=1 INFELE(4)=1 INFELE(5)=1 INFELE(6)=1 INFELE(7)=50 INFELE(9)=2 IF (IFOUR.EQ.1) INFELE(9)=3 INFELE(10)=1 INFELE(15)=2 IF (IFOUR.EQ.1) INFELE(15)=3 INFELE(16)=1 GOTO 300 C ================================================================= C TUYO C ================================================================= 96 CONTINUE NTETA=40 INFELE(2)=2*NTETA INFELE(3)=2*NTETA INFELE(4)=2*NTETA INFELE(5)=7 INFELE(6)=2*NTETA INFELE(7)=500 INFELE(8)=2+4 INFELE(9)=12+2*9 INFELE(10)=6 INFELE(15)=6+9 INFELE(16)=6 GOTO 300 C ================================================================= C LSE2 : TUYAU ACOUSTIQUE PURE C ================================================================= 97 CONTINUE INFELE(2)=2 INFELE(3)=2 INFELE(4)=2 INFELE(5)=7 INFELE(6)=2 INFELE(7)=200 INFELE(9)=4 INFELE(10)=2 INFELE(15)=2 INFELE(16)=2 GOTO 300 C ================================================================= C LITU : COUPLAGE LIQUIDE TUYAU C ================================================================= 98 CONTINUE INFELE(2)=2 INFELE(3)=2 INFELE(4)=2 INFELE(5)=6 INFELE(6)=2 INFELE(7)=200 INFELE(9)=10 INFELE(10)=2 INFELE(15)=5 INFELE(16)=5 GOTO 300 C ================================================================= C HYT3 : Element HYBRIDE TRI3 C ================================================================= 99 CONTINUE INFELE(3)=4 INFELE(6)=4 GOTO 300 C ================================================================= C HYQ4 : Element HYBRIDE QUA4 C ================================================================= 100 CONTINUE INFELE(3)=4 INFELE(6)=4 GOTO 300 C ================================================================= C HYT4 : Element HYBRIDE TET4 C ================================================================= 101 CONTINUE INFELE(3)=4 INFELE(6)=4 GOTO 300 C ================================================================= C HYP6 : Element HYBRIDE PRI6 C ================================================================= 102 CONTINUE INFELE(3)=6 INFELE(6)=6 GOTO 300 C ================================================================= C HYC8 : Element HYBRIDE CUB8 C ================================================================= 103 CONTINUE INFELE(3)=8 INFELE(6)=8 GOTO 300 C ================================================================= C TRIS (TRI3-SECTION) C ================================================================= 104 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(5)=2 INFELE(6)=4 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN INFELE(10)=2 INFELE(16)=2 ELSE IF (IFOUR.EQ.2) THEN INFELE(10)=3 INFELE(16)=3 ENDIF GOTO 300 C ================================================================= C QUAS (QUA4-SECTION) C ================================================================= 105 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(5)=2 INFELE(6)=4 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN INFELE(10)=2 INFELE(16)=2 ELSE IF (IFOUR.EQ.2) THEN INFELE(10)=3 INFELE(16)=3 ENDIF GOTO 300 C ================================================================= C POIS (POI1-SECTION) C ================================================================= 106 CONTINUE INFELE(2)=1 INFELE(3)=1 INFELE(4)=1 INFELE(5)=1 INFELE(6)=1 INFELE(10)=1 INFELE(16)=1 GOTO 300 C ================================================================= C JOP3 : JOINTS POREUX C ================================================================= 108 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(6)=3 INFELE(7)=100 INFELE(8)=12 INFELE(9)=18 INFELE(10)=2 INFELE(15)=3 INFELE(16)=3 GOTO 300 C ================================================================= C JOP6 : JOINTS POREUX C ================================================================= 109 CONTINUE INFELE(2)=7 INFELE(3)=7 INFELE(4)=7 INFELE(6)=7 INFELE(7)=100 INFELE(8)=21 INFELE(9)=45 INFELE(10)=3 INFELE(15)=4 INFELE(16)=4 GOTO 300 C ================================================================= C JOP8 : JOINTS POREUX C ================================================================= 110 CONTINUE INFELE(2)=9 INFELE(3)=9 INFELE(4)=9 INFELE(6)=9 INFELE(7)=100 INFELE(8)=28 INFELE(9)=60 INFELE(10)=3 INFELE(15)=4 INFELE(16)=4 GOTO 300 C ================================================================= C POLYGONE a N cotes C ================================================================= 111 CONTINUE INFELE(8) = MELE - 108 NBSH=INFELE(8) INFELE(2)=NBSH INFELE(3)=NBSH INFELE(4)=NBSH INFELE(6)=NBSH INFELE(9)=2*NBSH IF (IFOUR.EQ.1) INFELE(9)=3*NBSH IF (IFOUR.EQ.-3) INFELE(9)=2*NBSH+3 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C BAR3 (barre a 3 noeuds) C ================================================================= 123 CONTINUE NPOINT=3 INFELE(2)=NPOINT INFELE(4)=NPOINT INFELE(5)=1 INFELE(6)=NPOINT INFELE(7)=200 INFELE(9)=9 IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) INFELE(9)=6 INFELE(10)=1 INFELE(15)=3 IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) INFELE(15)=2 INFELE(16)=1 GOTO 300 C ================================================================= C BAEX (barre excentree 3D a 2 noeuds) C ================================================================= 124 CONTINUE INFELE(2)=2 INFELE(3)=2 INFELE(4)=2 INFELE(5)=6 INFELE(6)=2 INFELE(7)=700 INFELE(9)=12 INFELE(10)=1 INFELE(15)=6 INFELE(16)=1 GOTO 300 C ================================================================= C LIA2 : element de liaison a 2 noeuds (6 ddl par noeuds) C ================================================================= 125 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(5)=9 INFELE(6)=2 INFELE(7)=200 INFELE(9)=12 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C QUAH C ================================================================= 126 CONTINUE INFELE(2)=4 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(5)=5 INFELE(6)=4 INFELE(7)=0 IF (IFOUR.EQ.0) INFELE(9)=16 IF (IFOUR.EQ.1) INFELE(9)=24 INFELE(10)=10 INFELE(15)=6 INFELE(16)=0 GOTO 300 C ================================================================= C CUBH C ================================================================= 127 CONTINUE INFELE(2)=8 INFELE(3)=27 INFELE(5)=5 INFELE(6)=27 INFELE(7)=0 INFELE(9)=48 INFELE(10)=10 INFELE(15)=6 INFELE(16)=0 GOTO 300 C ================================================================= C ROT3 (1 pt GAUSS pr MASSE, 7 pts GAUSS pr RIGIDITE) C ================================================================= 128 CONTINUE INFELE(3)=1 INFELE(4)=INFELE(3) INFELE(6)=7 INFELE(9)=3 INFELE(15)=1 GOTO 300 C ================================================================= C TR6H C ================================================================= 157 CONTINUE INFELE(2)=4 INFELE(6)=6 INFELE(7)=0 IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN INFELE(3)=7 INFELE(5)=5 IF (IFOUR.EQ.0) INFELE(9)=12 IF (IFOUR.EQ.1) INFELE(9)=18 INFELE(15)=6 ELSE INFELE(3)=6 INFELE(5)=5 INFELE(9)=24 INFELE(15)=4 ENDIF INFELE(10)=10 INFELE(16)=0 GOTO 300 C ================================================================= C SEGS (SEG2-SECTION 2D) C ================================================================= 166 CONTINUE INFELE(2)=2 INFELE(3)=2 INFELE(4)=2 INFELE(5)=3 INFELE(6)=2 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN INFELE(10)=2 INFELE(16)=2 ELSE IF (IFOUR.EQ.2) THEN INFELE(10)=3 INFELE(16)=3 ENDIF GOTO 300 C ================================================================= C POJS (POI1-SECTION) C ================================================================= 167 CONTINUE INFELE(2)=1 INFELE(3)=1 INFELE(4)=1 INFELE(5)=3 INFELE(6)=1 IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN INFELE(10)=2 INFELE(16)=2 ELSE IF (IFOUR.EQ.2) THEN INFELE(10)=3 INFELE(16)=3 ENDIF GOTO 300 C ================================================================= C JCT3 C ================================================================= 168 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(6)=3 INFELE(7)=300 INFELE(9)=INFELE(8)*IDIM INFELE(10)=2 INFELE(15)=IDIM INFELE(16)=2 GOTO 300 C ================================================================= C JCI4 C ================================================================= 169 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(6)=4 INFELE(7)=300 INFELE(9)=INFELE(8)*IDIM INFELE(10)=2 INFELE(15)=IDIM INFELE(16)=2 GOTO 300 C ================================================================= C JGI2 C ================================================================= 170 CONTINUE INFELE(2)=4 INFELE(3)=2 INFELE(4)=2 INFELE(6)=2 IF (IFOUR.EQ.-3) INFELE(9)=INFELE(9)+3 INFELE(7)=300 INFELE(9)=INFELE(8)*IDIM INFELE(10)=4 INFELE(15)=IDIM INFELE(16)=4 GOTO 300 C ================================================================= C JGT3 C ================================================================= 171 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(6)=3 INFELE(7)=300 INFELE(9)=INFELE(8)*IDIM INFELE(10)=6 INFELE(15)=IDIM INFELE(16)=6 GOTO 300 C ================================================================= C JGI4 C ================================================================= 172 CONTINUE INFELE(2)=4 INFELE(3)=4 INFELE(4)=4 INFELE(6)=4 INFELE(7)=300 INFELE(9)=INFELE(8)*IDIM INFELE(10)=6 INFELE(15)=IDIM INFELE(16)=6 GOTO 300 C ================================================================= C TRIQ C ================================================================= 173 CONTINUE INFELE(3)=7 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=7 INFELE(4)=7 INFELE(6)=7 ELSE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 ENDIF INFELE(8)=9 INFELE(9)=18 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=24 IF (IFOUR.EQ.1) THEN INFELE(10)=8 INFELE(15)=5 INFELE(16)=8 ELSE INFELE(10)=6 INFELE(15)=4 INFELE(16)=6 ENDIF GOTO 300 C ================================================================= C QUAQ C ================================================================= 174 CONTINUE INFELE(2)=9 INFELE(3)=9 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=9 INFELE(6)=9 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(8)=12 INFELE(9)=24 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=32 IF (IFOUR.EQ.1) THEN INFELE(10)=8 INFELE(15)=5 INFELE(16)=8 ELSE INFELE(10)=6 INFELE(15)=4 INFELE(16)=6 ENDIF GOTO 300 C ================================================================= C CUBQ C ================================================================= 175 CONTINUE INFELE(2)=27 INFELE(3)=27 INFELE(4)=27 INFELE(6)=27 INFELE(8)=28 INFELE(9)=76 INFELE(10)=8 INFELE(15)=5 INFELE(16)=8 GOTO 300 C ================================================================= C TETQ C ================================================================= 176 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(8)=14 INFELE(9)=38 INFELE(10)=8 INFELE(15)=5 INFELE(16)=8 GOTO 300 C ================================================================= C PRIQ C ================================================================= 177 CONTINUE INFELE(2)=12 INFELE(3)=12 INFELE(4)=12 INFELE(6)=12 INFELE(8)=21 INFELE(9)=57 INFELE(10)=8 INFELE(15)=5 INFELE(16)=8 GOTO 300 C ================================================================= C TRIR C ================================================================= 178 CONTINUE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN INFELE(2)=7 INFELE(6)=7 INFELE(4)=7 ELSE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 ENDIF INFELE(3)=7 INFELE(8)=9 INFELE(9)=21 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=27 IF (IFOUR.EQ.1) THEN INFELE(10)=9 INFELE(15)=6 INFELE(16)=9 ELSE INFELE(10)=7 INFELE(15)=5 INFELE(16)=7 ENDIF GOTO 300 C ================================================================= C QUAR C ================================================================= 179 CONTINUE INFELE(2)=9 INFELE(3)=9 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=9 INFELE(6)=9 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(8)=12 INFELE(9)=28 IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) INFELE(9)=36 IF (IFOUR.EQ.1) THEN INFELE(10)=9 INFELE(15)=6 INFELE(16)=9 ELSE INFELE(10)=7 INFELE(15)=5 INFELE(16)=7 ENDIF GOTO 300 C ================================================================= C CUBR C ================================================================= 180 CONTINUE INFELE(2)=27 INFELE(3)=27 INFELE(4)=27 INFELE(6)=27 INFELE(8)=28 INFELE(9)=84 INFELE(10)=9 INFELE(15)=6 INFELE(16)=9 GOTO 300 C ================================================================= C TETR C ================================================================= 181 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(8)=14 INFELE(9)=42 INFELE(10)=9 INFELE(15)=6 INFELE(16)=9 GOTO 300 C ================================================================= C PRIR C ================================================================= 182 CONTINUE INFELE(2)=12 INFELE(3)=12 INFELE(4)=12 INFELE(6)=12 INFELE(8)=21 INFELE(9)=61 INFELE(10)=9 INFELE(15)=6 INFELE(16)=9 GOTO 300 C ================================================================= C Q4RI (QUA4 with 1x1 Gauss points) C ================================================================= 183 CONTINUE INFELE(2)=4 INFELE(3)=1 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=9 INFELE(4)=1 C** INFELE(4)=5 CORRECTION MILL LE 8/8/90 INFELE(6)=1 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=9 INFELE(9)=8 IF (IFOUR.EQ.-3) INFELE(9)=11 IF (IFOUR.EQ.1) INFELE(9)=12 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C Q8RI (QUA8 with 2x2 Gauss points) C ================================================================= 184 CONTINUE INFELE(2)=9 INFELE(3)=4 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(3)=16 INFELE(4)=4 INFELE(6)=4 C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) INFELE(6)=16 INFELE(9)=16 IF (IFOUR.EQ.-3) INFELE(9)=19 IF (IFOUR.EQ.1) INFELE(9)=24 IF (IFOUR.EQ.1) THEN INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 ELSE INFELE(10)=4 INFELE(15)=2 INFELE(16)=4 ENDIF GOTO 300 C ================================================================= C JOQ3 C ================================================================= 185 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(6)=3 INFELE(7)=150 INFELE(8)=12 INFELE(9)=24 INFELE(10)=4 INFELE(15)=4 INFELE(16)=4 GOTO 300 C ================================================================= C JOQ6 C ================================================================= 186 CONTINUE INFELE(2)=7 INFELE(3)=7 INFELE(4)=7 INFELE(6)=7 INFELE(7)=150 INFELE(8)=21 INFELE(9)=54 INFELE(10)=5 INFELE(15)=5 INFELE(16)=5 GOTO 300 C ================================================================= C JOQ8 C ================================================================= 187 CONTINUE INFELE(2)=9 INFELE(3)=9 INFELE(4)=9 INFELE(6)=9 INFELE(7)=150 INFELE(8)=28 INFELE(9)=72 INFELE(10)=5 INFELE(15)=5 INFELE(16)=5 GOTO 300 C ================================================================= C JOR3 C ================================================================= 188 CONTINUE INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(6)=3 INFELE(7)=200 INFELE(8)=12 INFELE(9)=30 INFELE(10)=5 INFELE(15)=5 INFELE(16)=5 GOTO 300 C ================================================================= C JOR6 C ================================================================= 189 CONTINUE INFELE(2)=7 INFELE(3)=7 INFELE(4)=7 INFELE(6)=7 INFELE(7)=200 INFELE(8)=21 INFELE(9)=63 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C JOR8 C ================================================================= 190 CONTINUE INFELE(2)=9 INFELE(3)=9 INFELE(4)=9 INFELE(6)=9 INFELE(7)=200 INFELE(8)=28 INFELE(9)=84 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 GOTO 300 C ================================================================= C T1D2/M1D2 : Element THERMIQUE/MECANIQUE MASSIF 1D a 2 noeuds C ================================================================= 191 CONTINUE C On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1) INFELE(2)=2 INFELE(3)=2 INFELE(4)=2 INFELE(5)=3 INFELE(6)=2 INFELE(7)=300 INFELE(9)=2 IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) INFELE(9)=3 IF (IFOUR.EQ.11) INFELE(9)=4 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C T1D3/M1D3 : Element THERMIQUE/MECANIQUE MASSIF 1D a 3 noeuds C ================================================================= 192 CONTINUE C On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1) INFELE(2)=3 INFELE(3)=3 INFELE(4)=3 INFELE(6)=3 INFELE(9)=3 IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) INFELE(9)=4 IF (IFOUR.EQ.11) INFELE(9)=5 INFELE(10)=3 INFELE(15)=1 INFELE(16)=3 GOTO 300 C ================================================================= C LC03 LS03 BS03 C ================================================================= 195 CONTINUE INFELE(6)=2 GOTO 300 C ================================================================= C LC07 LS07 BS07 C ================================================================= 196 CONTINUE INFELE(6)=7 GOTO 300 C ================================================================= C LC09 LS09 BS09 C ================================================================= 197 CONTINUE INFELE(6)=4 GOTO 300 C ================================================================= C MC03 C ================================================================= 216 CONTINUE INFELE(6)=2 GOTO 300 C ================================================================= C MC07 C ================================================================= 217 CONTINUE INFELE(6)=7 GOTO 300 C ================================================================= C MC09 C ================================================================= 218 CONTINUE INFELE(6)=4 GOTO 300 C ================================================================= C M103 MS03 C ================================================================= 223 CONTINUE INFELE(6)=3 GOTO 300 C ================================================================= C M107 MS07 C ================================================================= 224 CONTINUE INFELE(6)=7 GOTO 300 C ================================================================= C M109 MS09 C ================================================================= 225 CONTINUE INFELE(6)=16 GOTO 300 C ================================================================= C QC03 Q103 QS03 C ================================================================= 237 CONTINUE INFELE(6)=3 GOTO 300 C ================================================================= C QC07 Q107 QS07 C ================================================================= 238 CONTINUE INFELE(6)=7 GOTO 300 C ================================================================= C QC09 Q109 QS09 C ================================================================= 239 CONTINUE INFELE(6)=16 GOTO 300 C ================================================================= C CIFL C ================================================================= 258 CONTINUE INFELE(2)=1 INFELE(3)=1 INFELE(4)=1 INFELE(6)=1 INFELE(8)=0 INFELE(9)=8 INFELE(10)=5 INFELE(15)=3 INFELE(16)=5 GOTO 300 C ================================================================= C SHB8 C ================================================================= 260 CONTINUE INFELE(2)=5 INFELE(3)=5 INFELE(4)=5 INFELE(6)=5 INFELE(9)=24 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C XQ4R C ================================================================= 263 CONTINUE INFELE(2)=0 INFELE(3)=64 INFELE(4)=64 INFELE(6)=64 INFELE(8)=40 INFELE(10)=4 INFELE(15)=20 INFELE(16)=4 INFELE(9)=4*INFELE(15) GOTO 300 C ================================================================= C XC8R C ================================================================= 264 CONTINUE INFELE(2)=0 INFELE(3)=64 INFELE(4)=64 INFELE(6)=64 INFELE(8)=80 INFELE(10)=6 INFELE(15)=30 INFELE(16)=6 INFELE(9)=8*INFELE(15) GOTO 300 C ================================================================= C JOI1 : element de liaison a 2 noeuds (6 ddl par noeuds) C ================================================================= 265 CONTINUE INFELE(2)=1 INFELE(3)=1 INFELE(4)=1 INFELE(5)=0 INFELE(6)=1 INFELE(7)=700 IF(IFOUR.EQ.2) THEN INFELE(9)=12 INFELE(10)=6 INFELE(15)=6 INFELE(16)=6 ELSE IF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN INFELE(9)=6 INFELE(10)=3 INFELE(15)=3 INFELE(16)=3 ENDIF GOTO 300 C ================================================================= C ZCO2 : zone cohesive a 2 noeuds (2 ddls par noeud) C ================================================================= 266 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(6)=2 INFELE(7)=100 INFELE(9)=4 INFELE(10)=2 INFELE(15)=2 INFELE(16)=2 GOTO 300 C ================================================================= C ZCO3 : zone cohesive a 3 noeuds (2 ddls par noeud) C ================================================================= 267 CONTINUE INFELE(2)=1 INFELE(4)=1 INFELE(6)=1 INFELE(7)=100 INFELE(9)=9 INFELE(10)=3 INFELE(15)=3 INFELE(16)=3 GOTO 300 C ================================================================= C ZCO4 : zone cohesive a 4 noeuds (2 ddls par noeud) C ================================================================= 268 CONTINUE INFELE(2)=4 INFELE(4)=4 INFELE(6)=4 INFELE(7)=100 INFELE(9)=12 INFELE(10)=3 INFELE(15)=3 INFELE(16)=3 GOTO 300 C ================================================================= C TUY2 : tuyau pour modele d'advection en thermique C ================================================================= 269 CONTINUE INFELE(2)=2 INFELE(4)=2 INFELE(5)=5 INFELE(6)=2 INFELE(7)=200 INFELE(9)=2 INFELE(10)=1 INFELE(15)=1 INFELE(16)=1 GOTO 300 C ================================================================= C TUY3 : tuyau pour modele d'advection en thermique C ================================================================= 270 CONTINUE INFELE(2)=3 INFELE(4)=3 INFELE(5)=5 INFELE(6)=3 INFELE(7)=200 INFELE(9)=2 INFELE(10)=1 INFELE(15)=1 INFELE(16)=1 GOTO 300 C ================================================================= c COS2 (3D pour liaison acier-beton) C ================================================================= 271 CONTINUE INFELE(2)=1 INFELE(4)=2 INFELE(5)=2 INFELE(6)=2 INFELE(7)=100 INFELE(9)=4*idim INFELE(10)=3 INFELE(15)=IDIM INFELE(16)=3 GOTO 300 C ================================================================= c COA2 (3D pour liaison acier-beton) C ================================================================= 272 CONTINUE INFELE(2)=1 INFELE(4)=2 INFELE(5)=2 INFELE(6)=2 INFELE(7)=100 INFELE(9)=4*idim INFELE(10)=3 INFELE(15)=IDIM INFELE(16)=3 GOTO 300 C ================================================================= c CU27, comme le CU20 C ================================================================= 275 CONTINUE INFELE(2)=27 INFELE(3)=27 INFELE(4)=27 INFELE(6)=27 INFELE(9)=81 GOTO 300 C ================================================================= c PR21 C ================================================================= 276 CONTINUE INFELE(2)=21 cbp INFELE(3)=12 INFELE(3)=21 INFELE(4)=21 INFELE(6)=21 INFELE(9)=63 GOTO 300 C ================================================================= c TE15 C ================================================================= 277 CONTINUE INFELE(2)=15 cbp INFELE(3)=5 INFELE(3)=15 INFELE(4)=15 INFELE(6)=15 INFELE(9)=45 GOTO 300 C ================================================================= c PY19 C ================================================================= 278 CONTINUE INFELE(2)=17 INFELE(3)=17 INFELE(4)=17 INFELE(6)=17 INFELE(9)=57 GOTO 300 C ================================================================= C C20R C ================================================================= 279 CONTINUE INFELE(2)=27 INFELE(3)=8 INFELE(4)=8 INFELE(6)=8 INFELE(9)=60 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= C P15R C ================================================================= 280 CONTINUE INFELE(2)=12 INFELE(3)=8 INFELE(4)=6 INFELE(6)=6 INFELE(9)=45 INFELE(10)=6 INFELE(15)=3 INFELE(16)=6 GOTO 300 C ================================================================= 300 CONTINUE C IF (INFELE(3).EQ.0) INFELE(3)=INFELE(6) C C Cas particuliers des autres formulations massives IF (FORMOD(1).EQ.'ELECTROSTATIQUE') THEN INFELE(9) = NBSH INFELE(10)=IDIM INFELE(15)=1 INFELE(16)=IDIM ENDIF * nombre composante contraintes/defo, materiau IMPEDANCE if (dcmate) then MFR3=INFELE(13) call idprim(IMODEL,MFR3,IPNOMC,NBROBL,NBRFAC) if(tymode(/2).gt.0) then if (tymode(1).eq.'LISTMOTS') then dcmat2 = .true. mlmot5 = ivamod(1) segact mlmot5 nbrobl = mlmot5.mots(/2)*2 endif endif if (dcmat2.and.mele.eq.2) then INFELE(9) = nbrobl INFELE(10) = nbrobl/2 INFELE(15) = nbrobl/2 INFELE(16) = nbrobl/2 else INFELE(15) = nbrobl if (nbrfac.gt.0) INFELE(15) = INFELE(15) + nbrfac INFELE(16) = INFELE(15) if (mele.eq.45) INFELE(9) = nbrobl if (mele.eq.2) INFELE(9) = nbrobl*2 endif endif CJK148537 en attendant les fonctions de forme pr18, pr14, py19 if (MELE.EQ.140.OR.MELE.EQ.142.OR.MELE.EQ.149) GOTO 950 C INTTYP = -5 CALL ELQUOI(IMODEL,INTTYP,IPTR) INFELE(11) = IPTR C if (bmodi) imodel.nefmod=melep C 950 CONTINUE END