C IDVAR4 SOURCE YURI 11/10/28 22:15:22 7181 C======================================================================= * RECHERCHE DES NOMS DE COMPOSANTES DES VARIABLES INTERNES * -------------------------------------------------------- * * ENTREES : * --------- * MATEPL NUMERO DU MATERIAU PLASTIQUE * MFR NUMERO DE FORMULATION * IFOUR VALEUR DE IFOU DANS CCOPTIO * NPINT NOMBRE DE POINTS D'INTEGRATION DANS L'EPAISSEUR * POUR LES ELEMENTS DE COQUE INTEGRES * * SORTIES : * --------- * IPCOMP POINTEUR SUR LA LISTE DES NOMS DE COMPOSANTES * OBLIGATOIRES ET FACULTATIVES * NBROBL NOMBRE DE COMPOSANTES OBLIGATOIRES * NBRFAC NOMBRE DE COMPOSANTES FACTULTATIVES C======================================================================= SUBROUTINE IDVAR4 (MATEPL,MFR,IFOUR,NPINT,IPCOMP,NBROBL,NBRFAC) IMPLICIT INTEGER(I-N) SEGMENT NOMID CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC) ENDSEGMENT NBROBL=0 NBRFAC=0 NOMID=0 C 1 - Modele de CHABOCHE et de OHNO C =================================== IF (MATEPL.EQ.25) THEN C ===== C 1.1 - Elements MASSIFs C ===== IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.33) THEN C= 1.1.1 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou C= des DEFORMATIONS PLANES GENERALISEES IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN NBROBL=20 SEGINI,NOMID LESOBL( 1)='X1XX' LESOBL( 2)='X1YY' LESOBL( 3)='X1ZZ' LESOBL( 4)='X1XY' LESOBL( 5)='X2XX' LESOBL( 6)='X2YY' LESOBL( 7)='X2ZZ' LESOBL( 8)='X2XY' LESOBL( 9)='GPXX' LESOBL(10)='GPYY' LESOBL(11)='GPZZ' LESOBL(12)='GPXY' LESOBL(13)='VIXX' LESOBL(14)='VIYY' LESOBL(15)='VIZZ' LESOBL(16)='VIXY' LESOBL(17)='EPSE' LESOBL(18)='RR ' LESOBL(19)='QQQ ' LESOBL(20)='QQ ' C= 1.1.2 - Cas AXISYMETRIQUE ELSE IF (IFOUR.EQ.0) THEN NBROBL=20 SEGINI,NOMID LESOBL( 1)='X1RR' LESOBL( 2)='X1ZZ' LESOBL( 3)='X1TT' LESOBL( 4)='X1RZ' LESOBL( 5)='X2RR' LESOBL( 6)='X2ZZ' LESOBL( 7)='X2TT' LESOBL( 8)='X2RZ' LESOBL( 9)='GPRR' LESOBL(10)='GPZZ' LESOBL(11)='GPTT' LESOBL(12)='GPRZ' LESOBL(13)='VIRR' LESOBL(14)='VIZZ' LESOBL(15)='VITT' LESOBL(16)='VIRZ' LESOBL(17)='EPSE' LESOBL(18)='RR ' LESOBL(19)='QQQ ' LESOBL(20)='QQ ' C= 1.1.3 - Cas des series de FOURIER ELSE IF (IFOUR.EQ.1) THEN NBROBL=28 SEGINI,NOMID LESOBL( 1)='X1RR' LESOBL( 2)='X1ZZ' LESOBL( 3)='X1TT' LESOBL( 4)='X1RZ' LESOBL( 5)='X1RT' LESOBL( 6)='X1ZT' LESOBL( 7)='X2RR' LESOBL( 8)='X2ZZ' LESOBL( 9)='X2TT' LESOBL(10)='X2RZ' LESOBL(11)='X2RT' LESOBL(12)='X2ZT' LESOBL(13)='GPRR' LESOBL(14)='GPZZ' LESOBL(15)='GPTT' LESOBL(16)='GPRZ' LESOBL(17)='GPRT' LESOBL(18)='GPZT' LESOBL(19)='VIRR' LESOBL(20)='VIZZ' LESOBL(21)='VITT' LESOBL(22)='VIRZ' LESOBL(23)='VIRT' LESOBL(24)='VIZT' LESOBL(25)='EPSE' LESOBL(26)='RR ' LESOBL(27)='QQQ ' LESOBL(28)='QQ ' C= 1.1.4 - Cas TRIDIMENSIONNEL ELSE IF (IFOUR.EQ.2) THEN NBROBL=28 SEGINI,NOMID LESOBL( 1)='X1XX' LESOBL( 2)='X1YY' LESOBL( 3)='X1ZZ' LESOBL( 4)='X1XY' LESOBL( 5)='X1YZ' LESOBL( 6)='X1ZX' LESOBL( 7)='X2XX' LESOBL( 8)='X2YY' LESOBL( 9)='X2ZZ' LESOBL(10)='X2XY' LESOBL(11)='X2YZ' LESOBL(12)='X2ZX' LESOBL(13)='GPXX' LESOBL(14)='GPYY' LESOBL(15)='GPZZ' LESOBL(16)='GPXY' LESOBL(17)='GPYZ' LESOBL(18)='GPZX' LESOBL(19)='VIXX' LESOBL(20)='VIYY' LESOBL(21)='VIZZ' LESOBL(22)='VIXY' LESOBL(23)='VIYZ' LESOBL(24)='VIZX' LESOBL(25)='EPSE' LESOBL(26)='RR ' LESOBL(27)='QQQ ' LESOBL(28)='QQ ' C= 1.1.5 - Modes de calcul UNIDIMENSIONNELS (1D) PLAN ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN NBROBL=16 SEGINI,NOMID LESOBL( 1)='X1XX' LESOBL( 2)='X1YY' LESOBL( 3)='X1ZZ' LESOBL( 4)='X2XX' LESOBL( 5)='X2YY' LESOBL( 6)='X2ZZ' LESOBL( 7)='GPXX' LESOBL( 8)='GPYY' LESOBL( 9)='GPZZ' LESOBL(10)='VIXX' LESOBL(11)='VIYY' LESOBL(12)='VIZZ' LESOBL(13)='EPSE' LESOBL(14)='RR ' LESOBL(15)='QQQ ' LESOBL(16)='QQ ' C= 1.1.6 - Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN NBROBL=16 SEGINI,NOMID LESOBL( 1)='X1RR' LESOBL( 2)='X1ZZ' LESOBL( 3)='X1TT' LESOBL( 4)='X2RR' LESOBL( 5)='X2ZZ' LESOBL( 6)='X2TT' LESOBL( 7)='GPRR' LESOBL( 8)='GPZZ' LESOBL( 9)='GPTT' LESOBL(10)='VIRR' LESOBL(11)='VIZZ' LESOBL(12)='VITT' LESOBL(13)='EPSE' LESOBL(14)='RR ' LESOBL(15)='QQQ ' LESOBL(16)='QQ ' ENDIF C ===== C 1.2 - Elements COQUEs C ===== ELSE IF (MFR.EQ.3) THEN IF (NPINT.EQ.0) THEN C= 1.2.1 - Modes de calcul TRIDIMENSIONNEL et FOURIER IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN NBROBL=36 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1N4' LESOBL( 5)='X1M1' LESOBL( 6)='X1M2' LESOBL( 7)='X1M3' LESOBL( 8)='X1M4' LESOBL( 9)='X2N1' LESOBL(10)='X2N2' LESOBL(11)='X2N3' LESOBL(12)='X2N4' LESOBL(13)='X2M1' LESOBL(14)='X2M2' LESOBL(15)='X2M3' LESOBL(16)='X2M4' LESOBL(17)='GPN1' LESOBL(18)='GPN2' LESOBL(19)='GPN3' LESOBL(20)='GPN4' LESOBL(21)='GPM1' LESOBL(22)='GPM2' LESOBL(23)='GPM3' LESOBL(24)='GPM4' LESOBL(25)='EPN1' LESOBL(26)='EPN2' LESOBL(27)='EPN3' LESOBL(28)='EPN4' LESOBL(29)='EPM1' LESOBL(30)='EPM2' LESOBL(31)='EPM3' LESOBL(32)='EPM4' LESOBL(33)='EPSE' LESOBL(34)='RR ' LESOBL(35)='QQQ ' LESOBL(36)='QQ ' C= 1.2.2 - Cas AXISYMETRIQUE ELSE IF (IFOUR.EQ.0) THEN NBROBL=28 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1M1' LESOBL( 5)='X1M2' LESOBL( 6)='X1M3' LESOBL( 7)='X2N1' LESOBL( 8)='X2N2' LESOBL( 9)='X2N3' LESOBL(10)='X2M1' LESOBL(11)='X2M2' LESOBL(12)='X2M3' LESOBL(13)='GPN1' LESOBL(14)='GPN2' LESOBL(15)='GPN3' LESOBL(16)='GPM1' LESOBL(17)='GPM2' LESOBL(18)='GPM3' LESOBL(19)='EPN1' LESOBL(20)='EPN2' LESOBL(21)='EPN3' LESOBL(22)='EPM1' LESOBL(23)='EPM2' LESOBL(24)='EPM3' LESOBL(25)='EPSE' LESOBL(26)='RR ' LESOBL(27)='QQQ ' LESOBL(28)='QQ ' C= 1.2.3 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou C= des DEFORMATIONS PLANES GENERALISEES ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-3) THEN NBROBL=28 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1M1' LESOBL( 5)='X1M2' LESOBL( 6)='X1M3' LESOBL( 7)='X2N1' LESOBL( 8)='X2N2' LESOBL( 9)='X2N3' LESOBL(10)='X2M1' LESOBL(11)='X2M2' LESOBL(12)='X2M3' LESOBL(13)='GPN1' LESOBL(14)='GPN2' LESOBL(15)='GPN3' LESOBL(16)='GPM1' LESOBL(17)='GPM2' LESOBL(18)='GPM3' LESOBL(19)='EPN1' LESOBL(20)='EPN2' LESOBL(21)='EPN3' LESOBL(22)='EPM1' LESOBL(23)='EPM2' LESOBL(24)='EPM3' LESOBL(25)='EPSE' LESOBL(26)='RR ' LESOBL(27)='QQQ ' LESOBL(28)='QQ ' ENDIF ELSE C= 1.2.4 - Cas des COQUES INTEGREES en TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=20 SEGINI,NOMID LESOBL( 1)='X1SS' LESOBL( 2)='X1TT' LESOBL( 3)='X1NN' LESOBL( 4)='X1ST' LESOBL( 5)='X2SS' LESOBL( 6)='X2TT' LESOBL( 7)='X2NN' LESOBL( 8)='X2ST' LESOBL( 9)='GPSS' LESOBL(10)='GPTT' LESOBL(11)='GPNN' LESOBL(12)='GPST' LESOBL(13)='VISS' LESOBL(14)='VITT' LESOBL(15)='VINN' LESOBL(16)='VIST' LESOBL(17)='EPSE' LESOBL(18)='RR ' LESOBL(19)='QQQ ' LESOBL(20)='QQ ' ENDIF ENDIF C ===== C 1.3 - Element COQUE EPAISSE C ===== ELSE IF (MFR.EQ.5) THEN C= Mode de calcul TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=28 SEGINI,NOMID LESOBL( 1)='X1SS' LESOBL( 2)='X1TT' LESOBL( 3)='X1NN' LESOBL( 4)='X1ST' LESOBL( 5)='X1SN' LESOBL( 6)='X1TN' LESOBL( 7)='X2SS' LESOBL( 8)='X2TT' LESOBL( 9)='X2NN' LESOBL(10)='X2ST' LESOBL(11)='X2SN' LESOBL(12)='X2TN' LESOBL(13)='GPSS' LESOBL(14)='GPTT' LESOBL(15)='GPNN' LESOBL(16)='GPST' LESOBL(17)='GPSN' LESOBL(18)='GPTN' LESOBL(19)='VISS' LESOBL(20)='VITT' LESOBL(21)='VINN' LESOBL(22)='VIST' LESOBL(23)='VISN' LESOBL(24)='VITN' LESOBL(25)='EPSE' LESOBL(26)='RR ' LESOBL(27)='QQQ ' LESOBL(28)='QQ ' ENDIF C ===== C 1.4 - Elements POUTRE et TUYAU C ===== ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN C= Mode de calcul TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=28 SEGINI,NOMID LESOBL( 1)='X1FX' LESOBL( 2)='X1FY' LESOBL( 3)='X1FZ' LESOBL( 4)='X1MX' LESOBL( 5)='X1MY' LESOBL( 6)='X1MZ' LESOBL( 7)='X2FX' LESOBL( 8)='X2FY' LESOBL( 9)='X2FZ' LESOBL(10)='X2MX' LESOBL(11)='X2MY' LESOBL(12)='X2MZ' LESOBL(13)='GPFX' LESOBL(14)='GPFY' LESOBL(15)='GPFZ' LESOBL(16)='GPMX' LESOBL(17)='GPMY' LESOBL(18)='GPMZ' LESOBL(19)='EPFX' LESOBL(20)='EPFY' LESOBL(21)='EPFZ' LESOBL(22)='EPMX' LESOBL(23)='EPMY' LESOBL(24)='EPMZ' LESOBL(25)='EPSE' LESOBL(26)='RR ' LESOBL(27)='QQQ ' LESOBL(28)='QQ ' ENDIF C ===== C 1.5 - Element COQUE CISAILLEMENT TRANSVERSE C ===== ELSE IF (MFR.EQ.9) THEN C= Mode de calcul TRIDIMENSIONNEL IF (NPINT.EQ.0) THEN IF (IFOUR.EQ.2) THEN NBROBL=36 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1M1' LESOBL( 5)='X1M2' LESOBL( 6)='X1M3' LESOBL( 7)='X1V1' LESOBL( 8)='X1V2' LESOBL( 9)='X2N1' LESOBL(10)='X2N2' LESOBL(11)='X2N3' LESOBL(12)='X2M1' LESOBL(13)='X2M2' LESOBL(14)='X2M3' LESOBL(15)='X2V1' LESOBL(16)='X2V2' LESOBL(17)='GPN1' LESOBL(18)='GPN2' LESOBL(19)='GPN3' LESOBL(20)='GPM1' LESOBL(21)='GPM2' LESOBL(22)='GPM3' LESOBL(23)='GPV1' LESOBL(24)='GPV2' LESOBL(25)='EPN1' LESOBL(26)='EPN2' LESOBL(27)='EPN3' LESOBL(28)='EPM1' LESOBL(29)='EPM2' LESOBL(30)='EPM3' LESOBL(31)='EPV1' LESOBL(32)='EPV2' LESOBL(33)='EPSE' LESOBL(34)='RR ' LESOBL(35)='QQQ ' LESOBL(36)='QQ ' ENDIF ENDIF C ===== C 1.6 - Formulation MEMBRANE C ===== ELSE IF (MFR.EQ.25) THEN NBROBL=16 SEGINI,NOMID LESOBL( 1)='X1SS' LESOBL( 2)='X1TT' LESOBL( 3)='X1ST' LESOBL( 4)='X2SS' LESOBL( 5)='X2TT' LESOBL( 6)='X2ST' LESOBL( 7)='GPSS' LESOBL( 8)='GPTT' LESOBL( 9)='GPST' LESOBL(10)='VISS' LESOBL(11)='VITT' LESOBL(12)='VIST' LESOBL(13)='EPSE' LESOBL(14)='RR ' LESOBL(15)='QQQ ' LESOBL(16)='QQ ' C ===== C 1.7 - Formulation UNIAXIALE C ===== ELSE IF (MFR.EQ.27) THEN NBROBL=8 SEGINI,NOMID LESOBL(1)='X1FX' LESOBL(2)='X2FX' LESOBL(3)='GPFX' LESOBL(4)='EPFX' LESOBL(5)='EPSE' LESOBL(6)='RR ' LESOBL(7)='QQQ ' LESOBL(8)='QQ ' ENDIF C 2 - Modele VISCOPLASTIQUE PARFAIT C =================================== ELSE IF (MATEPL.EQ.43) THEN C= Formulation MASSIF : C= Modele non disponible en PLAN CONT et en FOURIER IF ( (MFR.EQ.1.OR.MFR.EQ.31) .AND. & (IFOUR.NE.-2.AND.IFOUR.NE.1) ) THEN NBROBL=1 SEGINI,NOMID LESOBL(1)='EPSE' ENDIF C 3 - Modele a deux deformations inelastiques (DDI) C =================================================== ELSE IF (MATEPL.EQ.63) THEN C ===== C 3.1 - Elements MASSIFs C ===== IF (MFR.EQ.1.OR.MFR.EQ.31.OR.MFR.EQ.33) THEN C= 3.1.1 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou C= des DEFORMATIONS PLANES GENERALISEES IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN NBROBL=45 SEGINI,NOMID LESOBL( 1)='X1XX' LESOBL( 2)='X1YY' LESOBL( 3)='X1ZZ' LESOBL( 4)='X1XY' LESOBL( 5)='X2XX' LESOBL( 6)='X2YY' LESOBL( 7)='X2ZZ' LESOBL( 8)='X2XY' LESOBL( 9)='A1XX' LESOBL(10)='A1YY' LESOBL(11)='A1ZZ' LESOBL(12)='A1XY' LESOBL(13)='A2XX' LESOBL(14)='A2YY' LESOBL(15)='A2ZZ' LESOBL(16)='A2XY' LESOBL(17)='VIXX' LESOBL(18)='VIYY' LESOBL(19)='VIZZ' LESOBL(20)='VIXY' LESOBL(21)='Y1XX' LESOBL(22)='Y1YY' LESOBL(23)='Y1ZZ' LESOBL(24)='Y1XY' LESOBL(25)='Y2XX' LESOBL(26)='Y2YY' LESOBL(27)='Y2ZZ' LESOBL(28)='Y2XY' LESOBL(29)='B1XX' LESOBL(30)='B1YY' LESOBL(31)='B1ZZ' LESOBL(32)='B1XY' LESOBL(33)='B2XX' LESOBL(34)='B2YY' LESOBL(35)='B2ZZ' LESOBL(36)='B2XY' LESOBL(37)='EVXX' LESOBL(38)='EVYY' LESOBL(39)='EVZZ' LESOBL(40)='EVXY' LESOBL(41)='VIPW' LESOBL(42)='RP ' LESOBL(43)='V ' LESOBL(44)='RV ' LESOBL(45)='EPSE' C= 3.1.2 - Cas AXISYMETRIQUE ELSE IF (IFOUR.EQ.0) THEN NBROBL=45 SEGINI,NOMID LESOBL( 1)='X1RR' LESOBL( 2)='X1ZZ' LESOBL( 3)='X1TT' LESOBL( 4)='X1RZ' LESOBL( 5)='X2RR' LESOBL( 6)='X2ZZ' LESOBL( 7)='X2TT' LESOBL( 8)='X2RZ' LESOBL( 9)='A1RR' LESOBL(10)='A1ZZ' LESOBL(11)='A1TT' LESOBL(12)='A1RZ' LESOBL(13)='A2RR' LESOBL(14)='A2ZZ' LESOBL(15)='A2TT' LESOBL(16)='A2RZ' LESOBL(17)='VIRR' LESOBL(18)='VIZZ' LESOBL(19)='VITT' LESOBL(20)='VIRZ' LESOBL(21)='Y1RR' LESOBL(22)='Y1ZZ' LESOBL(23)='Y1TT' LESOBL(24)='Y1RZ' LESOBL(25)='Y2RR' LESOBL(26)='Y2ZZ' LESOBL(27)='Y2TT' LESOBL(28)='Y2RZ' LESOBL(29)='B1RR' LESOBL(30)='B1ZZ' LESOBL(31)='B1TT' LESOBL(32)='B1RZ' LESOBL(33)='B2RR' LESOBL(34)='B2ZZ' LESOBL(35)='B2TT' LESOBL(36)='B2RZ' LESOBL(37)='EVRR' LESOBL(38)='EVZZ' LESOBL(39)='EVTT' LESOBL(40)='EVRZ' LESOBL(41)='VIPW' LESOBL(42)='RP ' LESOBL(43)='V ' LESOBL(44)='RV ' LESOBL(45)='EPSE' C= 3.1.3 - Cas des series de FOURIER ELSE IF (IFOUR.EQ.1) THEN NBROBL=65 SEGINI,NOMID LESOBL( 1)='X1RR' LESOBL( 2)='X1ZZ' LESOBL( 3)='X1TT' LESOBL( 4)='X1RZ' LESOBL( 5)='X1RT' LESOBL( 6)='X1ZT' LESOBL( 7)='X2RR' LESOBL( 8)='X2ZZ' LESOBL( 9)='X2TT' LESOBL(10)='X2RZ' LESOBL(11)='X2RT' LESOBL(12)='X2ZT' LESOBL(13)='A1RR' LESOBL(14)='A1ZZ' LESOBL(15)='A1TT' LESOBL(16)='A1RZ' LESOBL(17)='A1RT' LESOBL(18)='A1ZT' LESOBL(19)='A2RR' LESOBL(20)='A2ZZ' LESOBL(21)='A2TT' LESOBL(22)='A2RZ' LESOBL(23)='A2RT' LESOBL(24)='A2ZT' LESOBL(25)='VIRR' LESOBL(26)='VIZZ' LESOBL(27)='VITT' LESOBL(28)='VIRZ' LESOBL(29)='VIRT' LESOBL(30)='VIZT' LESOBL(31)='Y1RR' LESOBL(32)='Y1ZZ' LESOBL(33)='Y1TT' LESOBL(34)='Y1RZ' LESOBL(35)='Y1RT' LESOBL(36)='Y1ZT' LESOBL(37)='Y2RR' LESOBL(38)='Y2ZZ' LESOBL(39)='Y2TT' LESOBL(40)='Y2RZ' LESOBL(41)='Y2RT' LESOBL(42)='Y2ZT' LESOBL(43)='B1RR' LESOBL(44)='B1ZZ' LESOBL(45)='B1TT' LESOBL(46)='B1RZ' LESOBL(47)='B1RT' LESOBL(48)='B1ZT' LESOBL(49)='B2RR' LESOBL(50)='B2ZZ' LESOBL(51)='B2TT' LESOBL(52)='B2RZ' LESOBL(53)='B2RT' LESOBL(54)='B2ZT' LESOBL(55)='EVRR' LESOBL(56)='EVZZ' LESOBL(57)='EVTT' LESOBL(58)='EVRZ' LESOBL(59)='EVRT' LESOBL(60)='EVZT' LESOBL(61)='VIPW' LESOBL(62)='RP ' LESOBL(63)='V ' LESOBL(64)='RV ' LESOBL(65)='EPSE' C= 3.1.4 - Cas TRIDIMENSIONNEL ELSE IF (IFOUR.EQ.2) THEN NBROBL=65 SEGINI,NOMID LESOBL( 1)='X1XX' lESOBL( 2)='X1YY' LESOBL( 3)='X1ZZ' LESOBL( 4)='X1XY' LESOBL( 5)='X1YZ' LESOBL( 6)='X1ZX' LESOBL( 7)='X2XX' LESOBL( 8)='X2YY' LESOBL( 9)='X2ZZ' LESOBL(10)='X2XY' LESOBL(11)='X2YZ' LESOBL(12)='X2ZX' LESOBL(13)='A1XX' LESOBL(14)='A1YY' LESOBL(15)='A1ZZ' LESOBL(16)='A1XY' LESOBL(17)='A1YZ' LESOBL(18)='A1ZX' LESOBL(19)='A2XX' LESOBL(20)='A2YY' LESOBL(21)='A2ZZ' LESOBL(22)='A2XY' LESOBL(23)='A2YZ' LESOBL(24)='A2ZX' LESOBL(25)='VIXX' LESOBL(26)='VIYY' LESOBL(27)='VIZZ' LESOBL(28)='VIXY' LESOBL(29)='VIYZ' LESOBL(30)='VIZX' LESOBL(31)='Y1XX' LESOBL(32)='Y1YY' LESOBL(33)='Y1ZZ' LESOBL(34)='Y1XY' LESOBL(35)='Y1YZ' LESOBL(36)='Y1ZX' LESOBL(37)='Y2XX' LESOBL(38)='Y2YY' LESOBL(39)='Y2ZZ' LESOBL(40)='Y2XY' LESOBL(41)='Y2YZ' LESOBL(42)='Y2ZX' LESOBL(43)='B1XX' LESOBL(44)='B1YY' LESOBL(45)='B1ZZ' LESOBL(46)='B1XY' LESOBL(47)='B1YZ' LESOBL(48)='B1ZX' LESOBL(49)='B2XX' LESOBL(50)='B2YY' LESOBL(51)='B2ZZ' LESOBL(52)='B2XY' LESOBL(53)='B2YZ' LESOBL(54)='B2ZX' LESOBL(55)='EVXX' LESOBL(56)='EVYY' LESOBL(57)='EVZZ' LESOBL(58)='EVXY' LESOBL(59)='EVYZ' LESOBL(60)='EVZX' LESOBL(61)='VIPW' LESOBL(62)='RP ' LESOBL(63)='V ' LESOBL(64)='RV ' LESOBL(65)='EPSE' C= 3.1.5 - Modes de calcul UNIDIMENSIONNELS (1D) PLAN ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN NBROBL=35 SEGINI,NOMID LESOBL( 1)='X1XX' LESOBL( 2)='X1YY' LESOBL( 3)='X1ZZ' LESOBL( 4)='X2XX' LESOBL( 5)='X2YY' LESOBL( 6)='X2ZZ' LESOBL( 7)='A1XX' LESOBL( 8)='A1YY' LESOBL( 9)='A1ZZ' LESOBL(10)='A2XX' LESOBL(11)='A2YY' LESOBL(12)='A2ZZ' LESOBL(13)='VIXX' LESOBL(14)='VIYY' LESOBL(15)='VIZZ' LESOBL(16)='Y1XX' LESOBL(17)='Y1YY' LESOBL(18)='Y1ZZ' LESOBL(19)='Y2XX' LESOBL(20)='Y2YY' LESOBL(21)='Y2ZZ' LESOBL(22)='B1XX' LESOBL(23)='B1YY' LESOBL(24)='B1ZZ' LESOBL(25)='B2XX' LESOBL(26)='B2YY' LESOBL(27)='B2ZZ' LESOBL(28)='EVXX' LESOBL(29)='EVYY' LESOBL(30)='EVZZ' LESOBL(31)='VIPW' LESOBL(32)='RP ' LESOBL(33)='V ' LESOBL(34)='RV ' LESOBL(35)='EPSE' C= 3.1.6 - Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN NBROBL=35 SEGINI,NOMID LESOBL( 1)='X1RR' LESOBL( 2)='X1ZZ' LESOBL( 3)='X1TT' LESOBL( 4)='X2RR' LESOBL( 5)='X2ZZ' LESOBL( 6)='X2TT' LESOBL( 7)='A1RR' LESOBL( 8)='A1ZZ' LESOBL( 9)='A1TT' LESOBL(10)='A2RR' LESOBL(11)='A2ZZ' LESOBL(12)='A2TT' LESOBL(13)='VIRR' LESOBL(14)='VIZZ' LESOBL(15)='VITT' LESOBL(16)='Y1RR' LESOBL(17)='Y1ZZ' LESOBL(18)='Y1TT' LESOBL(19)='Y2RR' LESOBL(20)='Y2ZZ' LESOBL(21)='Y2TT' LESOBL(22)='B1RR' LESOBL(23)='B1ZZ' LESOBL(24)='B1TT' LESOBL(25)='B2RR' LESOBL(26)='B2ZZ' LESOBL(27)='B2TT' LESOBL(28)='EVRR' LESOBL(29)='EVZZ' LESOBL(30)='EVTT' LESOBL(31)='VIPW' LESOBL(32)='RP ' LESOBL(33)='V ' LESOBL(34)='RV ' LESOBL(35)='EPSE' ENDIF C ===== C 3.2 - Elements COQUEs C ===== ELSE IF (MFR.EQ.3) THEN IF (NPINT.EQ.0) THEN C= 3.2.1 - Modes de calcul TRIDIMENSIONNEL et FOURIER IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN NBROBL=85 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1N4' LESOBL( 5)='X1M1' LESOBL( 6)='X1M2' LESOBL( 7)='X1M3' LESOBL( 8)='X1M4' LESOBL( 9)='X2N1' LESOBL(10)='X2N2' LESOBL(11)='X2N3' LESOBL(12)='X2N4' LESOBL(13)='X2M1' LESOBL(14)='X2M2' LESOBL(15)='X2M3' LESOBL(16)='X2M4' LESOBL(17)='A1N1' LESOBL(18)='A1N2' LESOBL(19)='A1N3' LESOBL(20)='A1N4' LESOBL(21)='A1M1' LESOBL(22)='A1M2' LESOBL(23)='A1M3' LESOBL(24)='A1M4' LESOBL(25)='A2N1' LESOBL(26)='A2N2' LESOBL(27)='A2N3' LESOBL(28)='A2N4' LESOBL(29)='A2M1' LESOBL(30)='A2M2' LESOBL(31)='A2M3' LESOBL(32)='A2M4' LESOBL(33)='EPN1' LESOBL(34)='EPN2' LESOBL(35)='EPN3' LESOBL(36)='EPN4' LESOBL(37)='EPM1' LESOBL(38)='EPM2' LESOBL(39)='EPM3' LESOBL(40)='EPM4' LESOBL(41)='Y1N1' LESOBL(42)='Y1N2' LESOBL(43)='Y1N3' LESOBL(44)='Y1N4' LESOBL(45)='Y1M1' LESOBL(46)='Y1M2' LESOBL(47)='Y1M3' LESOBL(48)='Y1M4' LESOBL(49)='Y2N1' LESOBL(50)='Y2N2' LESOBL(51)='Y2N3' LESOBL(52)='Y2N4' LESOBL(53)='Y2M1' LESOBL(54)='Y2M2' LESOBL(55)='Y2M3' LESOBL(56)='Y2M4' LESOBL(57)='B1N1' LESOBL(58)='B1N2' LESOBL(59)='B1N3' LESOBL(60)='B1N4' LESOBL(61)='B1M1' LESOBL(62)='B1M2' LESOBL(63)='B1M3' LESOBL(64)='B1M4' LESOBL(65)='B2N1' LESOBL(66)='B2N2' LESOBL(67)='B2N3' LESOBL(68)='B2N4' LESOBL(69)='B2M1' LESOBL(70)='B2M2' LESOBL(71)='B2M3' LESOBL(72)='A2M4' LESOBL(73)='EVN1' LESOBL(74)='EVN2' LESOBL(75)='EVN3' LESOBL(76)='EVN4' LESOBL(77)='EVM1' LESOBL(78)='EVM2' LESOBL(79)='EVM3' LESOBL(80)='EVM4' LESOBL(81)='VIPW' LESOBL(82)='RP ' LESOBL(83)='V ' LESOBL(84)='RV ' LESOBL(85)='EPSE' C= 3.2.2 - Cas AXISYMETRIQUE ELSE IF (IFOUR.EQ.0) THEN NBROBL=65 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1M1' LESOBL( 5)='X1M2' LESOBL( 6)='X1M3' LESOBL( 7)='X2N1' LESOBL( 8)='X2N2' LESOBL( 9)='X2N3' LESOBL(10)='X2M1' LESOBL(11)='X2M2' LESOBL(12)='X2M3' LESOBL(13)='A1N1' LESOBL(14)='A1N2' LESOBL(15)='A1N3' LESOBL(16)='A1M1' LESOBL(17)='A1M2' LESOBL(18)='A1M3' LESOBL(19)='A2N1' LESOBL(20)='A2N2' LESOBL(21)='A2N3' LESOBL(22)='A2M1' LESOBL(23)='A2M2' LESOBL(24)='A2M3' LESOBL(25)='EPN1' LESOBL(26)='EPN2' LESOBL(27)='EPN3' LESOBL(28)='EPM1' LESOBL(29)='EPM2' LESOBL(30)='EPM3' LESOBL(31)='Y1N1' LESOBL(32)='Y1N2' LESOBL(33)='Y1N3' LESOBL(34)='Y1M1' LESOBL(35)='Y1M2' LESOBL(36)='Y1M3' LESOBL(37)='Y2N1' LESOBL(38)='Y2N2' LESOBL(39)='Y2N3' LESOBL(40)='Y2M1' LESOBL(41)='Y2M2' LESOBL(42)='Y2M3' LESOBL(43)='B1N1' LESOBL(44)='B1N2' LESOBL(45)='B1N3' LESOBL(46)='B1M1' LESOBL(47)='B1M2' LESOBL(48)='B1M3' LESOBL(49)='B2N1' LESOBL(50)='B2N2' LESOBL(51)='B2N3' LESOBL(52)='B2M1' LESOBL(53)='B2M2' LESOBL(54)='B2M3' LESOBL(55)='EVN1' LESOBL(56)='EVN2' LESOBL(57)='EVN3' LESOBL(58)='EVM1' LESOBL(59)='EVM2' LESOBL(60)='EVM3' LESOBL(61)='VIPW' LESOBL(62)='RP ' LESOBL(63)='V ' LESOBL(64)='RV ' LESOBL(65)='EPSE' C= 3.2.3 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou C= des DEFORMATIONS PLANES GENERALISEES ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-3) THEN NBROBL=65 SEGINI,NOMID LESOBL( 1)='X1N1' LESOBL( 2)='X1N2' LESOBL( 3)='X1N3' LESOBL( 4)='X1M1' LESOBL( 5)='X1M2' LESOBL( 6)='X1M3' LESOBL( 7)='X2N1' LESOBL( 8)='X2N2' LESOBL( 9)='X2N3' LESOBL(10)='X2M1' LESOBL(11)='X2M2' LESOBL(12)='X2M3' LESOBL(13)='A1N1' LESOBL(14)='A1N2' LESOBL(15)='A1N3' LESOBL(16)='A1M1' LESOBL(17)='A1M2' LESOBL(18)='A1M3' LESOBL(19)='A2N1' LESOBL(20)='A2N2' LESOBL(21)='A2N3' LESOBL(22)='A2M1' LESOBL(23)='A2M2' LESOBL(24)='A2M3' LESOBL(25)='EPN1' LESOBL(26)='EPN2' LESOBL(27)='EPN3' LESOBL(28)='EPM1' LESOBL(29)='EPM2' LESOBL(30)='EPM3' LESOBL(31)='Y1N1' LESOBL(32)='Y1N2' LESOBL(33)='Y1N3' LESOBL(34)='Y1M1' LESOBL(35)='Y1M2' LESOBL(36)='Y1M3' LESOBL(37)='Y2N1' LESOBL(38)='Y2N2' LESOBL(39)='Y2N3' LESOBL(40)='Y2M1' LESOBL(41)='Y2M2' LESOBL(42)='Y2M3' LESOBL(43)='B1N1' LESOBL(44)='B1N2' LESOBL(45)='B1N3' LESOBL(46)='B1M1' LESOBL(47)='B1M2' LESOBL(48)='B1M3' LESOBL(49)='B2N1' LESOBL(50)='B2N2' LESOBL(51)='B2N3' LESOBL(52)='B2M1' LESOBL(53)='B2M2' LESOBL(54)='B2M3' LESOBL(55)='EVN1' LESOBL(56)='EVN2' LESOBL(57)='EVN3' LESOBL(58)='EVM1' LESOBL(59)='EVM2' LESOBL(60)='EVM3' LESOBL(61)='VIPW' LESOBL(62)='RP ' LESOBL(63)='V ' LESOBL(64)='RV ' LESOBL(65)='EPSE' ENDIF ELSE C= 3.2.4 - Cas des COQUES INTEGREES en TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=45 SEGINI,NOMID LESOBL( 1)='X1SS' LESOBL( 2)='X1TT' LESOBL( 3)='X1NN' LESOBL( 4)='X1ST' LESOBL( 5)='X2SS' LESOBL( 6)='X2TT' LESOBL( 7)='X2NN' LESOBL( 8)='X2ST' LESOBL( 9)='A1SS' LESOBL(10)='A1TT' LESOBL(11)='A1NN' LESOBL(12)='A1ST' LESOBL(13)='A2SS' LESOBL(14)='A2TT' LESOBL(15)='A2NN' LESOBL(16)='A2ST' LESOBL(17)='VISS' LESOBL(18)='VITT' LESOBL(19)='VINN' LESOBL(20)='VIST' LESOBL(21)='Y1SS' LESOBL(22)='Y1TT' LESOBL(23)='Y1NN' LESOBL(24)='Y1ST' LESOBL(25)='Y2SS' LESOBL(26)='Y2TT' LESOBL(27)='Y2NN' LESOBL(28)='Y2ST' LESOBL(29)='B1SS' LESOBL(30)='B1TT' LESOBL(31)='B1NN' LESOBL(32)='B1ST' LESOBL(33)='B2SS' LESOBL(34)='B2TT' LESOBL(35)='B2NN' LESOBL(36)='B2ST' LESOBL(37)='EVSS' LESOBL(38)='EVTT' LESOBL(39)='EVNN' LESOBL(40)='EVST' LESOBL(41)='VIPW ' LESOBL(42)='RP ' LESOBL(43)='V ' LESOBL(44)='RV ' LESOBL(45)='EPSE' ENDIF ENDIF C ===== C 3.3 - Element COQUE EPAISSE C ===== ELSE IF (MFR.EQ.5) THEN C= Mode de calcul TRIDIMENSIONNEL NBROBL=65 SEGINI,NOMID LESOBL( 1)='X1SS' LESOBL( 2)='X1TT' LESOBL( 3)='X1NN' LESOBL( 4)='X1ST' LESOBL( 5)='X1SN' LESOBL( 6)='X1TN' LESOBL( 7)='X2SS' LESOBL( 8)='X2TT' LESOBL( 9)='X2NN' LESOBL(10)='X2ST' LESOBL(11)='X2SN' LESOBL(12)='X2TN' LESOBL(13)='A1SS' LESOBL(14)='A1TT' LESOBL(15)='A1NN' LESOBL(16)='A1ST' LESOBL(17)='A1SN' LESOBL(18)='A1TN' LESOBL(19)='A2SS' LESOBL(20)='A2TT' LESOBL(21)='A2NN' LESOBL(22)='A2ST' LESOBL(23)='A2SN' LESOBL(24)='A2TN' LESOBL(25)='VISS' LESOBL(26)='VITT' LESOBL(27)='VINN' LESOBL(28)='VIST' LESOBL(29)='VISN' LESOBL(30)='VITN' LESOBL(31)='Y1SS' LESOBL(32)='Y1TT' LESOBL(33)='Y1NN' LESOBL(34)='Y1ST' LESOBL(35)='Y1SN' LESOBL(36)='Y1TN' LESOBL(37)='Y2SS' LESOBL(38)='Y2TT' LESOBL(39)='Y2NN' LESOBL(40)='Y2ST' LESOBL(41)='Y2SN' LESOBL(42)='Y2TN' LESOBL(43)='B1SS' LESOBL(44)='B1TT' LESOBL(45)='B1NN' LESOBL(46)='B1ST' LESOBL(47)='B1SN' LESOBL(48)='B1TN' LESOBL(49)='B2SS' LESOBL(50)='B2TT' LESOBL(51)='B2NN' LESOBL(52)='B2ST' LESOBL(53)='B2SN' LESOBL(54)='B2TN' LESOBL(55)='EVSS' LESOBL(56)='EVTT' LESOBL(57)='EVNN' LESOBL(58)='EVST' LESOBL(59)='EVSN' LESOBL(60)='EVTN' LESOBL(61)='VIPW' LESOBL(62)='RP ' LESOBL(63)='V ' LESOBL(64)='RV ' LESOBL(65)='EPSE' ENDIF C 4 - Modele ELASTO-VISCOPLASTIQUE DE KOCKS C =========================================== ELSE IF (MATEPL.EQ.70) THEN C= Formulation MASSIF : C= Modele non disponible en PLAN CONT et en FOURIER IF (MFR.EQ.1.AND.(IFOUR.NE.-2.AND.IFOUR.NE.1)) THEN NBROBL=2 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='S ' ENDIF C 5 - Modele MISTRAL C ==================== ELSE IF (MATEPL.EQ.94) THEN C ===== C 5.1 - Elements MASSIFs (seule formulation disponible actuellement) C ===== C= Note : Les indices finaux 1, 2, 3, 4, 5, 6 C= correspondent respectivement aux composantes : C= 11, 22, 33, 12, 13, 23 IF (MFR.EQ.1.OR.MFR.EQ.31) THEN IF (IFOUR.NE.1) THEN NBROBL=101 SEGINI,NOMID LESOBL( 1)='ETH1' LESOBL( 2)='ETH2' LESOBL( 3)='ETH3' LESOBL( 4)='ETH4' LESOBL( 5)='ETH5' LESOBL( 6)='ETH6' LESOBL( 7)='EEL1' LESOBL( 8)='EEL2' LESOBL( 9)='EEL3' LESOBL(10)='EEL4' LESOBL(11)='EEL5' LESOBL(12)='EEL6' LESOBL(13)='EP01' LESOBL(14)='EP02' LESOBL(15)='EP03' LESOBL(16)='EP04' LESOBL(17)='EP05' LESOBL(18)='EP06' LESOBL(19)='EP11' LESOBL(20)='EP12' LESOBL(21)='EP13' LESOBL(22)='EP14' LESOBL(23)='EP15' LESOBL(24)='EP16' LESOBL(25)='EP21' LESOBL(26)='EP22' LESOBL(27)='EP23' LESOBL(28)='EP24' LESOBL(29)='EP25' LESOBL(30)='EP26' LESOBL(31)='EP31' LESOBL(32)='EP32' LESOBL(33)='EP33' LESOBL(34)='EP34' LESOBL(35)='EP35' LESOBL(36)='EP36' LESOBL(37)='EP41' LESOBL(38)='EP42' LESOBL(39)='EP43' LESOBL(40)='EP44' LESOBL(41)='EP45' LESOBL(42)='EP46' LESOBL(43)='ECR1' LESOBL(44)='ECR2' LESOBL(45)='ECR3' LESOBL(46)='ECR4' LESOBL(47)='ECR5' LESOBL(48)='ECR6' LESOBL(49)='EP0E' LESOBL(50)='EP1E' LESOBL(51)='EP2E' LESOBL(52)='EP3E' LESOBL(53)='EP4E' LESOBL(54)='EP0Q' LESOBL(55)='EP1Q' LESOBL(56)='EP2Q' LESOBL(57)='EP3Q' LESOBL(58)='EP4Q' LESOBL(59)='X11 ' LESOBL(60)='X12 ' LESOBL(61)='X13 ' LESOBL(62)='X14 ' LESOBL(63)='X15 ' LESOBL(64)='X16 ' LESOBL(65)='X21 ' LESOBL(66)='X22 ' LESOBL(67)='X23 ' LESOBL(68)='X24 ' LESOBL(69)='X25 ' LESOBL(70)='X26 ' LESOBL(71)='X31 ' LESOBL(72)='X32 ' LESOBL(73)='X33 ' LESOBL(74)='X34 ' LESOBL(75)='X35 ' LESOBL(76)='X36 ' LESOBL(77)='FIT ' LESOBL(78)='PSI ' LESOBL(79)='DTA ' LESOBL(80)='SEA0' LESOBL(81)='RRA0' LESOBL(82)='VEA0' LESOBL(83)='KPLA' LESOBL(84)='EPSE' LESOBL(85)='VET1' LESOBL(86)='VET2' LESOBL(87)='VET3' LESOBL(88)='VET4' LESOBL(89)='VET5' LESOBL(90)='VET6' LESOBL(91)='VEI1' LESOBL(92)='VEI2' LESOBL(93)='VEI3' LESOBL(94)='VEI4' LESOBL(95)='VEI5' LESOBL(96)='VEI6' LESOBL(97)='VP0E' LESOBL(98)='VP1E' LESOBL(99)='VP2E' LESOBL(100)='VP3E' LESOBL(101)='VP4E' ENDIF ENDIF C 6 - Modele de GATT_MONERIE C ============================ C= Modele disponible uniquement en formulation MASSIF ELSE IF (MATEPL.EQ.107) THEN IF (MFR.EQ.1.OR.MFR.EQ.31) THEN IF (IFOUR.NE.1) THEN NBROBL=6 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='PORO' LESOBL(3)='BU ' LESOBL(4)='TETA' LESOBL(5)='EPSD' LESOBL(6)='EPSG' ENDIF ENDIF C 7 - Modele UO2 (= OTTOSEN + GATT_MONERIE) C =========================================== ELSE IF (MATEPL.EQ.108) THEN C ===== C 7.1 - Elements MASSIFs C ===== IF (MFR.EQ.1.OR.MFR.EQ.31) THEN C= 7.1.1 - Cas TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=30 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='W3MA' LESOBL( 5)='W1 ' LESOBL( 6)='W2 ' LESOBL( 7)='W3 ' LESOBL( 8)='VF1X' LESOBL( 9)='VF1Y' LESOBL(10)='VF1Z' LESOBL(11)='VF2X' LESOBL(12)='VF2Y' LESOBL(13)='VF2Z' LESOBL(14)='VF3X' LESOBL(15)='VF3Y' LESOBL(16)='VF3Z' LESOBL(17)='IVL1' LESOBL(18)='IVL2' LESOBL(19)='IVL3' LESOBL(20)='PREC' LESOBL(21)='EVXX' LESOBL(22)='EVYY' LESOBL(23)='EVZZ' LESOBL(24)='GVXY' LESOBL(25)='GVXZ' LESOBL(26)='GVYZ' LESOBL(27)='PORO' LESOBL(28)='BU ' LESOBL(29)='EPSF' LESOBL(30)='TETA' C= 7.1.2 - Cas AXISYMETRIQUE ELSE IF (IFOUR.EQ.0) THEN NBROBL=23 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='EPSM' LESOBL( 5)='W1 ' LESOBL( 6)='W2 ' LESOBL( 7)='VEPS' LESOBL( 8)='VF1X' LESOBL( 9)='VF1Y' LESOBL(10)='VF2X' LESOBL(11)='VF2Y' LESOBL(12)='IVL1' LESOBL(13)='IVL2' LESOBL(14)='IVL3' LESOBL(15)='PREC' LESOBL(16)='EVRR' LESOBL(17)='EVZZ' LESOBL(18)='EVTT' LESOBL(19)='GVRZ' LESOBL(20)='PORO' LESOBL(21)='BU ' LESOBL(22)='EPSF' LESOBL(23)='TETA' C= 7.1.3 - Cas des series de FOURIER ELSE IF (IFOUR.EQ.1) THEN NBROBL=25 SEGINI NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='EPSM' LESOBL( 5)='W1 ' LESOBL( 6)='W2 ' LESOBL( 7)='VEPS' LESOBL( 8)='VF1X' LESOBL( 9)='VF1Y' LESOBL(10)='VF2X' LESOBL(11)='VF2Y' LESOBL(12)='IVL1' LESOBL(13)='IVL2' LESOBL(14)='IVL3' LESOBL(15)='PREC' LESOBL(16)='EVRR' LESOBL(17)='EVZZ' LESOBL(18)='EVTT' LESOBL(19)='GVRZ' LESOBL(20)='GVRT' LESOBL(21)='GVZT' LESOBL(22)='PORO' LESOBL(23)='BU ' LESOBL(24)='EPSF' LESOBL(25)='TETA' C= 7.1.4 - Cas DEFOrmations PLANES (GENEralisees ou non) ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN NBROBL=23 SEGINI NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='EPSM' LESOBL( 5)='W1 ' LESOBL( 6)='W2 ' LESOBL( 7)='VEPS' LESOBL( 8)='VF1X' LESOBL( 9)='VF1Y' LESOBL(10)='VF2X' LESOBL(11)='VF2Y' LESOBL(12)='IVL1' LESOBL(13)='IVL2' LESOBL(14)='IVL3' LESOBL(15)='PREC' LESOBL(16)='EVXX' LESOBL(17)='EVYY' LESOBL(18)='EVZZ' LESOBL(19)='GVXY' LESOBL(20)='PORO' LESOBL(21)='BU ' LESOBL(22)='EPSF' LESOBL(23)='TETA' C= 7.1.5 - Cas CONTraintes PLANes ELSE IF (IFOUR.EQ.-2) THEN NBROBL=20 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='W1 ' LESOBL( 5)='W2 ' LESOBL( 6)='VF1X' LESOBL( 7)='VF1Y' LESOBL( 8)='VF2X' LESOBL( 9)='VF2Y' LESOBL(10)='IVL1' LESOBL(11)='IVL2' LESOBL(12)='PREC' LESOBL(13)='EVXX' LESOBL(14)='EVYY' LESOBL(15)='EVZZ' LESOBL(16)='GVXY' LESOBL(17)='PORO' LESOBL(18)='BU ' LESOBL(19)='EPSF' LESOBL(20)='TETA' C= 7.1.6 - Modes de calcul UNIDIMENSIONNELs (1D) PLAN ELSE IF (IFOUR.GE.3.AND.IFOUR.LE.11) THEN C* NBROBL= C* SEGINI,NOMID C* LESOBL( 1)='EPSE' C* LESOBL( 2)='W1MA' C* LESOBL( 3)='W2MA' C* LESOBL( 4)='W1 ' C* LESOBL( 5)='W2 ' C* LESOBL( 6)='VF1X' C* LESOBL( 7)='VF1Y' C* LESOBL( 8)='VF2X' C* LESOBL( 9)='VF2Y' C* LESOBL(10)='IVL1' C* LESOBL(11)='IVL2' C* LESOBL(12)='PREC' C* LESOBL(13)='EVXX' C* LESOBL(14)='EVYY' C* LESOBL(15)='EVZZ' C* LESOBL(16)='PORO' C* LESOBL(17)='BU ' C* LESOBL(18)='EPSF' C* LESOBL(19)='TETA' C= 7.1.7 - Modes de calcul UNIDIMENSIONNELs (1D) AXIS et SPHE ELSE IF (IFOUR.GE.12.AND.IFOUR.LE.15) THEN C* NBROBL= C* SEGINI,NOMID C* LESOBL( 1)='EPSE' C* LESOBL( 2)='W1MA' C* LESOBL( 3)='W2MA' C* LESOBL( 4)='W1 ' C* LESOBL( 5)='W2 ' C* LESOBL( 6)='VF1X' C* LESOBL( 7)='VF1Y' C* LESOBL( 8)='VF2X' C* LESOBL( 9)='VF2Y' C* LESOBL(10)='IVL1' C* LESOBL(11)='IVL2' C* LESOBL(12)='PREC' C* LESOBL(13)='EVRR' C* LESOBL(14)='EVZZ' C* LESOBL(15)='EVTT' C* LESOBL(16)='PORO' C* LESOBL(17)='BU ' C* LESOBL(18)='EPSF' C* LESOBL(19)='TETA' ENDIF C ===== C 7.2 - Elements COQUEs MINCEs C ===== ELSE IF (MFR.EQ.3) THEN C= 7.2.1 - Mode de calcul TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=22 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='W1 ' LESOBL( 5)='W2 ' LESOBL( 6)='VF1X' LESOBL( 7)='VF1Y' LESOBL( 8)='VF2X' LESOBL( 9)='VF2Y' LESOBL(10)='IVL1' LESOBL(11)='IVL2' LESOBL(12)='PREC' LESOBL(13)='EVSS' LESOBL(14)='EVTT' LESOBL(15)='GVST' LESOBL(16)='RVSS' LESOBL(17)='RVTT' LESOBL(18)='RVST' LESOBL(19)='PORO' LESOBL(20)='BU ' LESOBL(21)='EPSF' LESOBL(22)='TETA' ENDIF C ===== C 7.3 - Elements COQUEs en CISAILLEMENT TRANSVERSE$ C ===== ELSE IF (MFR.EQ.9) THEN C= 7.3.1 - Mode de calcul TRIDIMENSIONNEL IF (IFOUR.EQ.2) THEN NBROBL=24 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='W1MA' LESOBL( 3)='W2MA' LESOBL( 4)='W1 ' LESOBL( 5)='W2 ' LESOBL( 6)='VF1X' LESOBL( 7)='VF1Y' LESOBL( 8)='VF2X' LESOBL( 9)='VF2Y' LESOBL(10)='IVL1' LESOBL(11)='IVL2' LESOBL(12)='PREC' LESOBL(13)='EVSS' LESOBL(14)='EVTT' LESOBL(15)='GVST' LESOBL(16)='RVSS' LESOBL(17)='RVTT' LESOBL(18)='RVST' LESOBL(19)='GVSN' LESOBL(20)='GVTN' LESOBL(21)='PORO' LESOBL(22)='BU ' LESOBL(23)='EPSF' LESOBL(24)='TETA' ENDIF ENDIF C 8 - Modele visqueux avec double endommagement C ============================ C= Modele disponible uniquement en formulation MASSIF ELSE IF (MATEPL.EQ.130) THEN IF (MFR.EQ.1.OR.MFR.EQ.31) THEN IF (IFOUR.EQ.0) THEN NBROBL=4 SEGINI NOMID LESOBL(1)='R' LESOBL(2)='EPSE' LESOBL(3)='DD' LESOBL(4)='DC' IPCOMP=NOMID SEGDES NOMID ELSE IF (IFOUR.EQ.2) THEN NBROBL=4 SEGINI NOMID LESOBL(1)='R' LESOBL(2)='EPSE' LESOBL(3)='DD' LESOBL(4)='DC' IPCOMP=NOMID SEGDES NOMID ENDIF ENDIF C 9 - Modele visqueux CHAB_SINH_R et CHAB_NOR_R C ============================================== C= Modele disponible uniquement en formulation MASSIF ELSE IF (MATEPL.EQ.136.OR.MATEPL.EQ.138) THEN IF (MFR.EQ.1.OR.MFR.EQ.31) THEN C--------AXIS IF (IFOUR.EQ.0) THEN NBROBL=2 SEGINI NOMID LESOBL(1)='R' LESOBL(2)='EPSE' IPCOMP=NOMID SEGDES NOMID C--------3D ELSE IF (IFOUR.EQ.2) THEN NBROBL=2 SEGINI NOMID LESOBL(1)='R' LESOBL(2)='EPSE' IPCOMP=NOMID SEGDES NOMID ENDIF ENDIF C 9 - Modele visqueux CHAB_SINH_X et CHAB_NOR_X C ============================================== C= Modele disponible uniquement en formulation MASSIF ELSE IF (MATEPL.EQ.137.OR.MATEPL.EQ.139) THEN IF (MFR.EQ.1.OR.MFR.EQ.31) THEN C--------AXIS IF (IFOUR.EQ.0) THEN NBROBL=10 SEGINI NOMID LESOBL( 1)='A1RR' LESOBL( 2)='A1ZZ' LESOBL( 3)='A1TT' LESOBL( 4)='A1RZ' LESOBL( 5)='A2RR' LESOBL( 6)='A2ZZ' LESOBL( 7)='A2TT' LESOBL( 8)='A2RZ' LESOBL(9)='r' LESOBL(10)='EPSE' IPCOMP=NOMID SEGDES NOMID C--------3D ELSE IF (IFOUR.EQ.2) THEN C write(6,*) 'on est en 3D',IFOUR NBROBL=14 SEGINI NOMID LESOBL( 1)='A1XX' LESOBL( 2)='A1YY' LESOBL( 3)='A1ZZ' LESOBL( 4)='A1XY' LESOBL( 5)='A1YZ' LESOBL( 6)='A1ZX' LESOBL( 7)='A2XX' LESOBL( 8)='A2YY' LESOBL( 9)='A2ZZ' LESOBL(10)='A2XY' LESOBL(11)='A2YZ' LESOBL(12)='A2ZX' LESOBL(13)='r' LESOBL(14)='EPSE' IPCOMP=NOMID SEGDES NOMID ENDIF ENDIF C 10 - Modeles de SYMONDS et COWPER (SYCO1 et SYCO2) C =================================== ELSE IF (MATEPL.EQ.153.OR.MATEPL.EQ.154) THEN C element massif ou xfem IF (MFR.EQ.1.OR.MFR.EQ.63) THEN C Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou C= tridimensionnel IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.2) THEN NBROBL=2 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='VP' ENDIF ENDIF ENDIF IPCOMP=NOMID IF (IPCOMP.NE.0) THEN SEGDES,NOMID ELSE CALL ERREUR(19) ENDIF RETURN END