idvar6
C IDVAR6 SOURCE JB251061 22/12/09 21:15:01 11515 *_______________________________________________________________________ * 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 *_______________________________________________________________________ . NBRFAC) IMPLICIT INTEGER(I-N) SEGMENT NOMID CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC) ENDSEGMENT CHARACTER*4 NOMBET(17),NOMTUP(5) CHARACTER*4 NOMIN1(32),NOMIN2(6),NOMIN4(4),NOMI11(14) DATA NOMBET/'EPSE','DAM1','DAM2','DAM3', . 'OUV1','OUV2','OUV3','LAMD', . 'VI1X','VI1Y','VI1Z', . 'VI2X','VI2Y','VI2Z', . 'VI3X','VI3Y','VI3Z'/ DATA NOMTUP/'VAM1','VAM2','VHTA','VHJP','EPSE'/ DATA NOMIN1/'EPSE','TAN1','TAN2','IPL1','IFI2','RTM1', . 'RTM2','IFI1','JFI1','VHGL','JFI2','EQC1', . 'EPE1','EPT1','EPT2','EQC2','EPE2','EPC1', . 'EPC2','IPL2','OUV1','EDC1','EDC2','ETS1', . 'ETS2','EDT1','EDT2','OUV2','DFR1','DFR2', . 'VHAS','VHMP'/ DATA NOMI11/'EQC3','EPE3','EPT3','EPC3','DFR3','RTM3', . 'EDC3','ETS3','EDT3','TAN3','OUV3','IPL3', . 'IFI3','JFI3'/ DATA NOMIN2/'EP01','EP02','EP03','EP04','EP05','EP06'/ DATA NOMIN4/'EPSE','EPT1','EPT2','EPT3'/ IPCOMP=0 NBROBL=0 NBRFAC=0 NOMID=0 C write(6,*) ' matepl ', matepl C Materiau ELASTIQUE : EPSE facultative IF (MATEPL.EQ.0) THEN NBRFAC = 1 SEGINI,NOMID LESFAC(1)='EPSE' GOTO 666 ENDIF C= Materiaux 'VISCO_EXTERNE' : le GOTO indexe ne marche pas car C= =========================== valeur negative de MATEPL IF (MATEPL.EQ.-101) THEN NBROBL=4 SEGINI,NOMID LESOBL(1)='EC0 ' LESOBL(2)='ESW0' LESOBL(3)='P ' LESOBL(4)='QTLD' GOTO 666 ENDIF C on ne traite ici que les 50 premiers matepl la suite C est dans idva6b, idva6c ... GOTO (1, 2, 1, 4, 1,97, 4,97, 9,97,97, 4,97,14, 1,16,17,97,97,97, . 97,97,97,97,97, 3,97,28,29,97,97,32,33,34,35,36,97,38,39,40, . 41,42,97,44,45,46,47,48,49,97,98,98,97,97,97,97,97,97,97,97, . 97,98,97,98,98,98,98,98,98,97,98,97,97,97,98,98,97,98,98,98, . 97,97,98,97,97,97,98,97,97,97,97,98,98,97,97,97,97,97,97,98, . 98,97,29,98,98,98,97,97,97,97,97,97,97,97,97,97,97,97,97,98, . 98,98,97,97,97,97,98,98, 1, 1, 1, 1,98, 1,98, 1, 1, 1, 1,98, . 98,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98,98, . 98,98,98,98,98,98,98,98,98,98,98,98,98,98,98), . MATEPL 98 CONTINUE . NBRFAC) RETURN 1 CONTINUE NBROBL=1 SEGINI,NOMID LESOBL(1)='EPSE' GOTO 666 2 CONTINUE NBROBL=2 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)=NOMTUP(4) GOTO 666 C 3 - MATERIAUX ENDOMMAGEMENT ET ECROUISSAGE ISOTROPES SI MATEPL=26 C =================================================================== 3 CONTINUE C ===== C 3.1 - Element MASSIF C ===== IF (MFR.NE.1.AND.MFR.NE.31) GOTO 1201 C= 3.1.1 - Cas des CONTRAINTES PLANES ou des DEFORMATIONS PLANES ou C= des DEFORMATIONS PLANES GENERALISEES IF (IFOUR.NE.-2.AND.IFOUR.NE.-1.AND.IFOUR.NE.-3) GOTO 1101 NBROBL=7 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VIXX' LESOBL(5)='VIYY' LESOBL(6)='VIZZ' LESOBL(7)='VIXY' GOTO 666 C= 3.1.2 - Cas AXISYMETRIQUE 1101 IF (IFOUR.NE.0) GOTO 1111 NBROBL=7 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VIRR' LESOBL(5)='VIZZ' LESOBL(6)='VITT' LESOBL(7)='VIRZ' GOTO 666 C= 3.1.3 - Cas des series de FOURIER 1111 IF (IFOUR.NE.1) GOTO 1121 NBROBL=9 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VIRR' LESOBL(5)='VIZZ' LESOBL(6)='VITT' LESOBL(7)='VIRZ' LESOBL(8)='VIRT' LESOBL(9)='VIZT' GOTO 666 C= 3.1.4 - Cas TRIDIMENSIONNEL 1121 IF (IFOUR.NE.2) GOTO 1131 NBROBL=9 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VIXX' LESOBL(5)='VIYY' LESOBL(6)='VIZZ' LESOBL(7)='VIXY' LESOBL(8)='VIXZ' LESOBL(9)='VIYZ' GOTO 666 C= 3.1.5 - Modes de calcul UNIDIMENSIONNELS (1D) PLAN 1131 IF (IFOUR.LT.3.OR.IFOUR.GT.11) GOTO 1141 NBROBL=6 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VIXX' LESOBL(5)='VIYY' LESOBL(6)='VIZZ' GOTO 666 C= 3.1.6 - Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE 1141 IF (IFOUR.LT.12.OR.IFOUR.GT.15) GOTO 999 NBROBL=6 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VIRR' LESOBL(5)='VIZZ' LESOBL(6)='VITT' GOTO 666 C COQUE C 1201 CONTINUE IF (NPINT.NE.0) THEN IF (MFR.EQ.3.AND.IFOUR.EQ.2) THEN NBROBL=7 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' LESOBL(4)='VISS' LESOBL(5)='VITT' LESOBL(6)='VINN' LESOBL(7)='VIST' GOTO 666 ENDIF ENDIF GOTO 999 C C MATERIAU GAUVAIN POUR LES POUTRES UNIQUEMENT C 16 CONTINUE IF (MFR.NE.7) GOTO 999 NBROBL=13 SEGINI NOMID LESOBL( 1)='EPSE' LESOBL( 2)='EPY0' LESOBL( 3)='SYIN' LESOBL( 4)='EYIN' LESOBL( 5)='SYSU' LESOBL( 6)='EYSU' LESOBL( 7)='NUMY' LESOBL( 8)='EPZ0' LESOBL( 9)='SZIN' LESOBL(10)='EZIN' LESOBL(11)='SZSU' LESOBL(12)='EZSU' LESOBL(13)='NUMZ' GOTO 666 C 4 CONTINUE C 17 CONTINUE IF (MFR.NE.1.AND.MFR.NE.31.AND.MFR.NE.33.AND.MFR.NE.63) GOTO 1001 C C ELEMENT MASSIF C IF (IFOUR.NE.-2.AND.IFOUR.NE.-1.AND.IFOUR.NE.-3) GOTO 1040 C C CAS CONTRAINTES PLANES OU DEFORMATIONS PLANES C OU DEFORMATIONS PLANES C NBROBL=5 IF (MATEPL.EQ.12) NBROBL=9 IF (MATEPL.EQ.17) NBROBL=22 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMXX' LESOBL(3)='VMYY' LESOBL(4)='VMZZ' LESOBL(5)='VMXY' IF (MATEPL.EQ.12) THEN LESOBL(6)='XMXX' LESOBL(7)='XMYY' LESOBL(8)='XMZZ' LESOBL(9)='XMXY' ELSE IF (MATEPL.EQ.17) THEN LESOBL( 2)='VHP1' LESOBL( 3)='VLP2' LESOBL( 4)='VLP4' LESOBL( 5)='VMXX' LESOBL( 6)='VMYY' LESOBL( 7)='VMZZ' LESOBL( 8)='VMXY' LESOBL( 9)='VP1M' LESOBL(10)='VHWK' LESOBL(11)='VIXX' LESOBL(12)='VIYY' LESOBL(13)='VIZZ' LESOBL(14)='VIXY' LESOBL(15)='KVXX' LESOBL(16)='KVYY' LESOBL(17)='KVZZ' LESOBL(18)='KVXY' LESOBL(19)='KDXX' LESOBL(20)='KDYY' LESOBL(21)='KDZZ' LESOBL(22)='KDXY' ENDIF GOTO 666 C 1040 IF (IFOUR.NE.0) GOTO 1041 C C CAS AXISYMETRIQUE C NBROBL=5 IF(MATEPL.EQ.12) NBROBL=9 IF(MATEPL.EQ.17) NBROBL=22 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMRR' LESOBL(3)='VMZZ' LESOBL(4)='VMTT' LESOBL(5)='VMRZ' IF (MATEPL.EQ.12) THEN LESOBL(6)='XMRR' LESOBL(7)='XMZZ' LESOBL(8)='XMTT' LESOBL(9)='XMRZ' ELSE IF (MATEPL.EQ.17) THEN LESOBL( 2)='VHP1' LESOBL( 3)='VLP2' LESOBL( 4)='VLP4' LESOBL( 5)='VMRR' LESOBL( 6)='VMZZ' LESOBL( 7)='VMTT' LESOBL( 8)='VMRZ' LESOBL( 9)='VP1M' LESOBL(10)='VHWK' LESOBL(11)='VIRR' LESOBL(12)='VIZZ' LESOBL(13)='VITT' LESOBL(14)='VIRZ' LESOBL(15)='KVRR' LESOBL(16)='KVZZ' LESOBL(17)='KVTT' LESOBL(18)='KVRZ' LESOBL(19)='KDRR' LESOBL(20)='KDZZ' LESOBL(21)='KDTT' LESOBL(22)='KDRZ' ENDIF GOTO 666 C 1041 IF (IFOUR.NE.1) GOTO 1042 C C CAS SERIE DE FOURIER C NBROBL=7 IF (MATEPL.EQ.12) NBROBL=13 IF (MATEPL.EQ.17) NBROBL=30 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMRR' LESOBL(3)='VMZZ' LESOBL(4)='VMTT' LESOBL(5)='VMRZ' LESOBL(6)='VMRT' LESOBL(7)='VMZT' IF (MATEPL.EQ.12) THEN LESOBL( 8)='XMRR' LESOBL( 9)='XMZZ' LESOBL(10)='XMTT' LESOBL(11)='XMRZ' LESOBL(12)='XMRT' LESOBL(13)='XMZT' ELSE IF (MATEPL.EQ.17) THEN LESOBL( 2)='VHP1' LESOBL( 3)='VLP2' LESOBL( 4)='VLP4' LESOBL( 5)='VMRR' LESOBL( 6)='VMZZ' LESOBL( 7)='VMTT' LESOBL( 8)='VMRZ' LESOBL( 9)='VMRT' LESOBL(10)='VMZT' LESOBL(11)='VP1M' LESOBL(12)='VHWK' LESOBL(13)='VIRR' LESOBL(14)='VIZZ' LESOBL(15)='VITT' LESOBL(16)='VIRZ' LESOBL(17)='VIRT' LESOBL(18)='VIZT' LESOBL(19)='KVRR' LESOBL(20)='KVZZ' LESOBL(21)='KVTT' LESOBL(22)='KVRZ' LESOBL(23)='KVRT' LESOBL(24)='KVZT' LESOBL(25)='KDRR' LESOBL(26)='KDZZ' LESOBL(27)='KDTT' LESOBL(28)='KDRZ' LESOBL(29)='KDRT' LESOBL(30)='KDZT' ENDIF GOTO 666 C 1042 IF(IFOUR.NE.2) GOTO 1043 C C CAS TRIDIMENSIONNEL C NBROBL=7 IF(MATEPL.EQ.12) NBROBL=13 IF(MATEPL.EQ.17) NBROBL=30 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMXX' LESOBL(3)='VMYY' LESOBL(4)='VMZZ' LESOBL(5)='VMXY' LESOBL(6)='VMXZ' LESOBL(7)='VMYZ' IF (MATEPL.EQ.12) THEN LESOBL( 8)='XMXX' LESOBL( 9)='XMYY' LESOBL(10)='XMZZ' LESOBL(11)='XMXY' LESOBL(12)='XMXZ' LESOBL(13)='XMYZ' ELSE IF (MATEPL.EQ.17) THEN LESOBL( 2)='VHP1' LESOBL( 3)='VLP2' LESOBL( 4)='VLP4' LESOBL( 5)='VMXX' LESOBL( 6)='VMYY' LESOBL( 7)='VMZZ' LESOBL( 8)='VMXY' LESOBL( 9)='VMXZ' LESOBL(10)='VMYZ' LESOBL(11)='VP1M' LESOBL(12)='VHWK' LESOBL(13)='VIXX' LESOBL(14)='VIYY' LESOBL(15)='VIZZ' LESOBL(16)='VIXY' LESOBL(17)='VIXZ' LESOBL(18)='VIYZ' LESOBL(19)='KVXX' LESOBL(20)='KVYY' LESOBL(21)='KVZZ' LESOBL(22)='KVXY' LESOBL(23)='KVXZ' LESOBL(24)='KVYZ' LESOBL(25)='KDXX' LESOBL(26)='KDYY' LESOBL(27)='KDZZ' LESOBL(28)='KDXY' LESOBL(29)='KDXZ' LESOBL(30)='KDYZ' ENDIF GOTO 666 C= Modes de calcul UNIDIMENSIONNELS (1D) PLAN 1043 IF (IFOUR.LT.3.OR.IFOUR.GT.11) GOTO 1044 IF (MATEPL.EQ.12) THEN NBROBL=7 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VMXX' LESOBL(3)='VMYY' LESOBL(4)='VMZZ' LESOBL(5)='XMXX' LESOBL(6)='XMYY' LESOBL(7)='XMZZ' ELSE IF (MATEPL.EQ.17) THEN NBROBL=18 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='VHP1' LESOBL( 3)='VLP2' LESOBL( 4)='VLP4' LESOBL( 5)='VMXX' LESOBL( 6)='VMYY' LESOBL( 7)='VMZZ' LESOBL( 8)='VP1M' LESOBL( 9)='VHWK' LESOBL(10)='VIXX' LESOBL(11)='VIYY' LESOBL(12)='VIZZ' LESOBL(13)='KVXX' LESOBL(14)='KVYY' LESOBL(15)='KVZZ' LESOBL(16)='KDXX' LESOBL(17)='KDYY' LESOBL(18)='KDZZ' ELSE NBROBL=4 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VMXX' LESOBL(3)='VMYY' LESOBL(4)='VMZZ' ENDIF GOTO 666 C= Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE 1044 IF (IFOUR.LT.12.OR.IFOUR.GT.15) GOTO 999 IF (MATEPL.EQ.12) THEN NBROBL=7 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VMRR' LESOBL(3)='VMZZ' LESOBL(4)='VMTT' LESOBL(5)='XMRR' LESOBL(6)='XMZZ' LESOBL(7)='XMTT' ELSE IF (MATEPL.EQ.17) THEN NBROBL=18 SEGINI,NOMID LESOBL( 1)='EPSE' LESOBL( 2)='VHP1' LESOBL( 3)='VLP2' LESOBL( 4)='VLP4' LESOBL( 5)='VMRR' LESOBL( 6)='VMZZ' LESOBL( 7)='VMTT' LESOBL( 8)='VP1M' LESOBL( 9)='VHWK' LESOBL(10)='VIRR' LESOBL(11)='VIZZ' LESOBL(12)='VITT' LESOBL(13)='KVRR' LESOBL(14)='KVZZ' LESOBL(15)='KVTT' LESOBL(16)='KDRR' LESOBL(17)='KDZZ' LESOBL(18)='KDTT' ELSE NBROBL=4 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VMRR' LESOBL(3)='VMZZ' LESOBL(4)='VMTT' ENDIF GOTO 666 C 44 CONTINUE IF (MFR.NE.1.AND.MFR.NE.31.AND.MFR.NE.33) GOTO 999 IF (IFOUR.NE.-2.AND.IFOUR.NE.-1.AND.IFOUR.NE.-3) & GOTO 1330 C C CAS CONTRAINTES PLANES OU DEFORMATIONS PLANES C OU DEFORMATIONS PLANES GENERALISEES C NBROBL=6 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VIXX' LESOBL(3)='VIYY' LESOBL(4)='VIZZ' LESOBL(5)='VIXY' LESOBL(6)='VHOR' GOTO 666 C 1330 IF (IFOUR.NE.0) GOTO 1331 C C CAS AXISYMETRIQUE C NBROBL=6 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VIRR' LESOBL(3)='VIZZ' LESOBL(4)='VITT' LESOBL(5)='VIRZ' LESOBL(6)='VHOR' GOTO 666 C 1331 IF(IFOUR.NE.2) GOTO 1332 C C CAS TRIDIMENSIONNEL C NBROBL=8 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VIXX' LESOBL(3)='VIYY' LESOBL(4)='VIZZ' LESOBL(5)='VIXY' LESOBL(6)='VIXZ' LESOBL(7)='VIYZ' LESOBL(8)='VHOR' GOTO 666 C= Modes de calcul UNIDIMENSIONNELS (1D) PLAN 1332 IF (IFOUR.LT.3.OR.IFOUR.GT.11) GOTO 1333 NBROBL=5 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VIXX' LESOBL(3)='VIYY' LESOBL(4)='VIZZ' LESOBL(5)='VHOR' GOTO 666 C= Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE 1333 IF (IFOUR.LT.12.OR.IFOUR.GT.15) GOTO 999 NBROBL=5 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VIRR' LESOBL(3)='VIZZ' LESOBL(4)='VITT' LESOBL(5)='VHOR' GOTO 666 C 45 CONTINUE IF (MFR.NE.1.AND.MFR.NE.31.AND.MFR.NE.33) GOTO 999 IF (IFOUR.NE.-2.AND.IFOUR.NE.-1.AND.IFOUR.NE.-3) & GOTO 1430 C C CAS CONTRAINTES PLANES OU DEFORMATIONS PLANES C OU DEFORMATIONS PLANES GENERALISEES C NBROBL=8 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VIXX' LESOBL(3)='VIYY' LESOBL(4)='VIZZ' LESOBL(5)='VIXY' LESOBL(6)='GRAI' LESOBL(7)='DENE' LESOBL(8)='VHOR' GOTO 666 C 1430 IF (IFOUR.NE.0) GOTO 1431 C C CAS AXISYMETRIQUE C NBROBL=8 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VIRR' LESOBL(3)='VIZZ' LESOBL(4)='VITT' LESOBL(5)='VIRZ' LESOBL(6)='GRAI' LESOBL(7)='DENE' LESOBL(8)='VHOR' GOTO 666 C 1431 IF(IFOUR.NE.2) GOTO 1432 C C CAS TRIDIMENSIONNEL C NBROBL=1 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VIXX' LESOBL(3)='VIYY' LESOBL(4)='VIZZ' LESOBL(5)='VIXY' LESOBL(6)='VIXZ' LESOBL(7)='VIYZ' LESOBL(8)='GRAI' LESOBL(9)='DENE' LESOBL(10)='VHOR' GOTO 666 C= Modes de calcul UNIDIMENSIONNELS (1D) PLAN 1432 IF (IFOUR.LT.3.OR.IFOUR.GT.11) GOTO 1433 NBROBL=7 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VIXX' LESOBL(3)='VIYY' LESOBL(4)='VIZZ' LESOBL(5)='GRAI' LESOBL(6)='DENE' LESOBL(7)='VHOR' GOTO 666 C= Modes de calcul UNIDIMENSIONNELS (1D) AXIS et SPHE 1433 IF (IFOUR.LT.12.OR.IFOUR.GT.15) GOTO 999 NBROBL=7 SEGINI,NOMID LESOBL(1)='EPSE' LESOBL(2)='VIRR' LESOBL(3)='VIZZ' LESOBL(4)='VITT' LESOBL(5)='GRAI' LESOBL(6)='DENE' LESOBL(7)='VHOR' GOTO 666 C 1001 IF(MFR.NE.3) GOTO 1003 C C COQUE C IF(NPINT.EQ.0)THEN C COQUE GLOBALE IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN C C CAS TRIDIMENSIONNEL FOURIER C NBROBL=7 IF(MATEPL.EQ.12) NBROBL=13 IF(MATEPL.EQ.17) NBROBL=19 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VN11' LESOBL(3)='VN22' LESOBL(4)='VN12' LESOBL(5)='VM11' LESOBL(6)='VM22' LESOBL(7)='VM12' IF(MATEPL.EQ.12) THEN LESOBL(8)='XN11' LESOBL(9)='XN22' LESOBL(10)='XN12' LESOBL(11)='XM11' LESOBL(12)='XM22' LESOBL(13)='XM12' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VN11' LESOBL(6)='VN22' LESOBL(7)='VN12' LESOBL(8)='VM11' LESOBL(9)='VM22' LESOBL(10)='VM12' LESOBL(11)='VP1M' LESOBL(12)='VHWK' LESOBL(13)='VEP1' LESOBL(14)='VISS' LESOBL(15)='VITT' LESOBL(16)='VIST' LESOBL(17)='VRSS' LESOBL(18)='VRTT' LESOBL(19)='VRST' ENDIF C ELSE IF(IFOUR.EQ.0) THEN C C CAS AXISYMETRIQUE C NBROBL=5 IF(MATEPL.EQ.12) NBROBL=9 IF(MATEPL.EQ.17) NBROBL=15 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VN11' LESOBL(3)='VN22' LESOBL(4)='VM11' LESOBL(5)='VM22' IF(MATEPL.EQ.12) THEN LESOBL(6)='XN11' LESOBL(7)='XN22' LESOBL(8)='XM11' LESOBL(9)='XM22' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VN11' LESOBL(6)='VN22' LESOBL(7)='VM11' LESOBL(8)='VM22' LESOBL(9)='VP1M' LESOBL(10)='VHWK' LESOBL(11)='VEP1' LESOBL(12)='VISS' LESOBL(13)='VITT' LESOBL(14)='VRSS' LESOBL(15)='VRTT' ENDIF C ELSE IF(IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-3) THEN C C CAS CONTRAINTES OU DEFORMATIONS PLANES C OU DEFORMATIONS PLANES GENERALISEES C NBROBL=5 IF(MATEPL.EQ.12) NBROBL=9 IF(MATEPL.EQ.17) NBROBL=15 SEGINI NOMID LESOBL(1)='VISE' LESOBL(2)='VN11' LESOBL(3)='VNZZ' LESOBL(4)='VM11' LESOBL(5)='VMZZ' IF(MATEPL.EQ.12) THEN LESOBL(6)='XN11' LESOBL(7)='NZZX' LESOBL(8)='XM11' LESOBL(9)='MZZX' ELSE IF (MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VN11' LESOBL(6)='VNZZ' LESOBL(7)='VM11' LESOBL(8)='VMZZ' LESOBL(9)='VP1M' LESOBL(10)='VHWK' LESOBL(11)='VEP1' LESOBL(12)='VISS' LESOBL(13)='VIZZ' LESOBL(14)='VRSS' LESOBL(15)='VRZZ' ENDIF ENDIF GOTO 666 ELSE C COQUE INTEGREE IF(IFOUR.EQ.2)THEN NBROBL=5 IF(MATEPL.EQ.12) NBROBL=9 IF(MATEPL.EQ.17) NBROBL=15 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMSS' LESOBL(3)='VMTT' LESOBL(4)='VMNN' LESOBL(5)='VMST' IF(MATEPL.EQ.12) THEN LESOBL(6)='XMSS' LESOBL(7)='XMTT' LESOBL(8)='XMNN' LESOBL(9)='XMST' ELSE IF (MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VMSS' LESOBL(6)='VMTT' LESOBL(7)='VMNN' LESOBL(8)='VMST' LESOBL(9)='VP1M' LESOBL(10)='VHWK' LESOBL(11)='VEP1' LESOBL(12)='VISS' LESOBL(13)='VITT' LESOBL(14)='VINN' LESOBL(15)='VIST' ENDIF GOTO 666 ELSE GOTO 999 ENDIF ENDIF C 1003 IF(MFR.NE.5) GOTO 1005 C C COQUE EPAISSE C IF (IFOUR.NE.2) GOTO 999 C C CAS TRIDIMENSIONNEL C NBROBL=6 IF(MATEPL.EQ.12) NBROBL=11 IF(MATEPL.EQ.17) NBROBL=30 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMSS' LESOBL(3)='VMTT' LESOBL(4)='VMST' LESOBL(5)='VMSN' LESOBL(6)='VMTN' IF(MATEPL.EQ.12) THEN LESOBL(7)='XMSS' LESOBL(8)='XMTT' LESOBL(9)='XMST' LESOBL(10)='XMSN' LESOBL(11)='XMTN' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VMSS' LESOBL(6)='VMTT' LESOBL(7)='VMNN' LESOBL(8)='VMST' LESOBL(9)='VMSN' LESOBL(10)='VMTN' LESOBL(11)='VP1M' LESOBL(12)='VHWK' LESOBL(13)='VISS' LESOBL(14)='VITT' LESOBL(15)='VINN' LESOBL(16)='VIST' LESOBL(17)='VISN' LESOBL(18)='VITN' LESOBL(19)='KVSS' LESOBL(20)='KVTT' LESOBL(21)='KVNN' LESOBL(22)='KVST' LESOBL(23)='KVSN' LESOBL(24)='KVTN' LESOBL(25)='KDSS' LESOBL(26)='KDTT' LESOBL(27)='KDNN' LESOBL(28)='KDST' LESOBL(29)='KDSN' LESOBL(30)='KDTN' ENDIF GOTO 666 C 1005 IF(MFR.NE.7.AND.MFR.NE.13) GOTO 1007 C C POUTRE ET TUYAU C IF (IFOUR.NE.2) GOTO 999 C C CAS TRIDIMENSIONNEL C NBROBL=7 IF(MATEPL.EQ.12) NBROBL=13 IF(MATEPL.EQ.17) NBROBL=19 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VFFX' LESOBL(3)='VFFY' LESOBL(4)='VFFZ' LESOBL(5)='VOMX' LESOBL(6)='VOMY' LESOBL(7)='VOMZ' IF(MATEPL.EQ.12) THEN LESOBL(8)='XEFX' LESOBL(9)='XEFY' LESOBL(10)='XEFZ' LESOBL(11)='XMOX' LESOBL(12)='XMOY' LESOBL(13)='XMOZ' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VFFX' LESOBL(6)='VFFY' LESOBL(7)='VFFZ' LESOBL(8)='VOMX' LESOBL(9)='VOMY' LESOBL(10)='VOMZ' LESOBL(11)='VP1M' LESOBL(12)='VHWK' LESOBL(13)='VEP1' LESOBL(14)='VEPS' LESOBL(15)='VGXY' LESOBL(16)='VGXZ' LESOBL(17)='VHCX' LESOBL(18)='VHCY' LESOBL(19)='VHCZ' ENDIF GOTO 666 C 1007 IF(MFR.NE.9) GOTO 1009 C C COQUE CISAILLEMENT TRANSVERSE C IF(NPINT.EQ.0)THEN IF (IFOUR.NE.2) GOTO 999 C C CAS TRIDIMENSIONNEL C NBROBL=9 IF(MATEPL.EQ.12) NBROBL=17 IF(MATEPL.EQ.17) NBROBL=23 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VN11' LESOBL(3)='VN22' LESOBL(4)='VN12' LESOBL(5)='VM11' LESOBL(6)='VM22' LESOBL(7)='VM12' LESOBL(8)='VHV1' LESOBL(9)='VHV2' IF(MATEPL.EQ.12) THEN LESOBL(10)='XN11' LESOBL(11)='XN22' LESOBL(12)='XN12' LESOBL(13)='XM11' LESOBL(14)='XM22' LESOBL(15)='XM12' LESOBL(16)='VXV1' LESOBL(17)='VXV2' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VN11' LESOBL(6)='VN22' LESOBL(7)='VN12' LESOBL(8)='VM11' LESOBL(9)='VM22' LESOBL(10)='VM12' LESOBL(11)='VHV1' LESOBL(12)='VHV2' LESOBL(13)='VP1M' LESOBL(14)='VHWK' LESOBL(15)='VEP1' LESOBL(16)='VISS' LESOBL(17)='VITT' LESOBL(18)='VIST' LESOBL(19)='VRSS' LESOBL(20)='VRTT' LESOBL(21)='VRST' LESOBL(22)='VISN' LESOBL(23)='VITN' ENDIF GOTO 666 ELSE GOTO 999 ENDIF C 1009 CONTINUE C C FORMULATION MEMBRANE C IF(MFR.NE.25) GOTO 1025 NBROBL=4 IF(MATEPL.EQ.12) NBROBL=7 IF(MATEPL.EQ.17) NBROBL=13 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VMSS' LESOBL(3)='VMTT' LESOBL(4)='VMST' IF(MATEPL.EQ.12) THEN LESOBL(5)='XMSS' LESOBL(6)='XMTT' LESOBL(7)='XMST' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VMSS' LESOBL(6)='VMTT' LESOBL(7)='VMST' LESOBL(8)='VP1M' LESOBL(9)='VHWK' LESOBL(10)='VEP1' LESOBL(11)='VISS' LESOBL(12)='VITT' LESOBL(13)='VIST' ENDIF GOTO 666 C C FORMULATION UNIAXIALE C 1025 CONTINUE IF(MFR.NE.27) GOTO 999 NBROBL=2 IF(MATEPL.EQ.12) NBROBL=3 IF(MATEPL.EQ.17) NBROBL=9 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VFFX' IF(MATEPL.EQ.12) THEN LESOBL(3)='XEFX' ELSE IF(MATEPL.EQ.17) THEN LESOBL(2)='VHP1' LESOBL(3)='VLP2' LESOBL(4)='VLP4' LESOBL(5)='VFFX' LESOBL(6)='VP1M' LESOBL(7)='VHWK' LESOBL(8)='VEP1' LESOBL(9)='VISS' ENDIF GOTO 666 C 9 CONTINUE NBRFAC=17 SEGINI NOMID DO 991 I=1,NBRFAC LESFAC(I)=NOMBET(I) 991 CONTINUE GOTO 666 C 14 CONTINUE NBROBL=5 SEGINI NOMID LESOBL(1)=NOMTUP(1) LESOBL(2)=NOMTUP(2) LESOBL(3)=NOMTUP(3) LESOBL(4)=NOMTUP(4) LESOBL(5)=NOMTUP(5) GOTO 666 C C MATERIAU UBIQUITOUS C 28 CONTINUE NBROBL=7 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='EPN1' LESOBL(3)='GAP1' LESOBL(4)='LAM1' LESOBL(5)='EPN2' LESOBL(6)='GAP2' LESOBL(7)='LAM2' GOTO 666 C C MATERIAU VISCOPLASTIQUE ENDOMMAGEABLE C 29 CONTINUE NBROBL=3 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VVIR' LESOBL(3)='VHWD' GOTO 666 C C MATERIAU GLOBAL POUR LES POUTRES UNIQUEMENT C 32 CONTINUE IF(MFR.NE.7) GOTO 999 NBROBL=31 SEGINI NOMID C LESOBL(1)='EPSE' C LESOBL(2)='FLGC' LESOBL(3)='RCIN' LESOBL(4)='FMXC' LESOBL(5)='DMXC' LESOBL(6)='FMNC' LESOBL(7)='DMNC' LESOBL(8)='DEPC' LESOBL(9)='XINT' LESOBL(10)='PLAC' C LESOBL(11)='FLGF' LESOBL(12)='RFIN' LESOBL(13)='FMXF' LESOBL(14)='DMXF' LESOBL(15)='FMNF' LESOBL(16)='DMNF' LESOBL(17)='X00F' LESOBL(18)='DEPF' LESOBL(19)='PLAF' C LESOBL(20)='FLGT' LESOBL(21)='RTIN' LESOBL(22)='FMXT' LESOBL(23)='DMXT' LESOBL(24)='FMNT' LESOBL(25)='DMNT' LESOBL(26)='X01T' LESOBL(27)='X02T' LESOBL(28)='FINT' LESOBL(29)='FLIMT' LESOBL(30)='DEPT' LESOBL(31)='PLAT' GOTO 666 C C CAM-CLAY C 33 CONTINUE NBROBL=2 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='TRAP' GOTO 666 C C MATERIAU MOHR COULOMB C 34 CONTINUE IF (MFR.EQ.35) THEN NBROBL=4 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='EPOU' LESOBL(3)='STAT' LESOBL(4)='LAM1' GOTO 666 ELSE IF (MFR.EQ.75) THEN NBROBL=9 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='EPOU' LESOBL(3)='STAT' LESOBL(4)='LAM1' LESOBL(5)='PLA1' LESOBL(6)='PLA2' LESOBL(7)='PLA3' LESOBL(8)='EPCO' LESOBL(9)='LAM3' GOTO 666 ELSE GOTO 999 ENDIF C C MATERIAU JOINT-DILATANT (FORMULATION JOINT 2D-PLAN) C 35 CONTINUE IF(MFR.NE.35) GOTO 999 C IF (IFOUR.NE.-2.AND.IFOUR.NE.-1) GOTO 999 C NBROBL=10 SEGINI NOMID LESOBL(1) ='VDS0' LESOBL(2) ='VDN0' LESOBL(3) ='VRSN' LESOBL(4) ='VDRN' LESOBL(5) ='DPSN' LESOBL(6) ='VDPN' LESOBL(7) ='HC11' LESOBL(8) ='HC12' LESOBL(9) ='HC21' LESOBL(10)='HC22' GOTO 666 C C MATERIAU BETON_AXI C 36 CONTINUE NBROBL=10 SEGINI NOMID LESOBL(1) ='EPSE' LESOBL(2) ='DAM1' LESOBL(3) ='DAM2' LESOBL(4) ='DAM3' LESOBL(5) ='VHGL' LESOBL(6) ='LAM1' LESOBL(7) ='LAM2' LESOBL(8) ='LAM3' LESOBL(9) ='LAM4' LESOBL(10)='LAM5' GOTO 666 C C MATERIAU DE GURSON C 38 CONTINUE IF ((IFOUR.NE.0).AND.(IFOUR.NE.2) & .AND.(IFOUR.NE.-1)) GOTO 999 NBROBL=3 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='VPOR' LESOBL(3)='VHNS' GOTO 666 C C MATERIAU BETON_UNI C 39 CONTINUE NBROBL=6 SEGINI NOMID LESOBL( 1) ='EPSO' LESOBL( 2) ='EMAX' LESOBL( 3) ='EMIN' LESOBL( 4) ='STMN' LESOBL( 5) ='VHNG' LESOBL( 6) ='VHND' GOTO 666 C C MATERIAU ACIER_UNI C 40 CONTINUE NBROBL=12 SEGINI NOMID LESOBL( 1) ='FLAG' LESOBL( 2) ='EPSO' LESOBL( 3) ='RFAC' LESOBL( 4) ='VHNG' LESOBL( 5) ='EPLR' LESOBL( 6) ='EPIP' LESOBL( 7) ='EPMP' LESOBL( 8) ='EPMN' LESOBL( 9) ='STLR' LESOBL(10) ='STIP' LESOBL(11) ='STMP' LESOBL(12) ='STMN' GO TO 666 C C FORMULATION SECTION C 41 CONTINUE NBROBL=2 SEGINI NOMID LESOBL(1) ='VONS' LESOBL(2) ='VAIS' GOTO 666 C----------------------------- C FORMULATION OTTOSEN C----------------------------- 42 CONTINUE C C-------CAS TRIDIMENSIONNEL MASSIF------ C IF (IFOUR.EQ.2.AND.MFR.EQ.1) THEN NBROBL=20 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)='LBDC' GOTO 666 ENDIF C C-------CAS DEFO PLANES/AXISYMETRIQUE-------- C IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-3 & .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN NBROBL=15 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)='LBDC' GOTO 666 ENDIF C C-------CAS CONT PLANES-------- C -----OU CAS TRIDIMENSIONNEL COQUES------ C IF (IFOUR.EQ.-2.OR. . (IFOUR.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9))) THEN NBROBL=12 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)='LBDC' GOTO 666 ENDIF C*OF --> A voir pour le 1D C C MATERIAU MAZARS_FIB C 46 CONTINUE NBROBL=5 SEGINI NOMID LESOBL( 1) ='EPTI' LESOBL( 2) ='VHNI' LESOBL( 3) ='VISX' LESOBL( 4) ='VISY' LESOBL( 5) ='VISZ' GOTO 666 C C MATERIAU AMADEI_SAEB C 47 CONTINUE IF (MFR.NE.35) GOTO 999 NBROBL=7 SEGINI NOMID LESOBL( 1) ='EPSE' LESOBL( 2) ='EPOU' LESOBL( 3) ='STAT' LESOBL( 4) ='UEQU' LESOBL( 5) ='EPS1' LESOBL( 6) ='EPS2' LESOBL( 7) ='EPS3' GOTO 666 C C HUJEUX C 48 CONTINUE NBROBL=2 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='TRAP' GOTO 666 C 49 CONTINUE NBROBL=14 SEGINI NOMID LESOBL(1)='EPSE' LESOBL(2)='EPS1' LESOBL(3)='EPS2' LESOBL(4)='EPS3' LESOBL(5)='EPS4' LESOBL(6)='EPS5' LESOBL(7)='EPS6' LESOBL(8)='GTIN' LESOBL(9)='EPP1' LESOBL(10)='EPP2' LESOBL(11)='EPP3' LESOBL(12)='EPP4' LESOBL(13)='EPP5' LESOBL(14)='EPP6' GOTO 666 C C****************************************************** 97 continue 999 CONTINUE SEGINI,NOMID 666 CONTINUE IPCOMP=NOMID SEGDES,NOMID RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales