C IDPLAS SOURCE BR232186 13/12/09 21:15:12 7841 SUBROUTINE IDPLAS (MFR,MOOBL,MOFAC,IPLAC,IRET) *--------------------------------------------------------------------* * Noms de composantes de mat{riaux plastiques * *--------------------------------------------------------------------* * * * ENTREES: * * ------- * * MFR num{ro de formulation * * IPMODL objet mod}le {l{mentaire ( segment actif ) * * MOOBL pointeur sur les mots obligatoires * * MOFAC pointeur sur les mots facultatif * * IPLAC numero du materiau dans MODPLA * * * * SORTIE : * * ------- * * * * IRET : 0 SI PROBLEME * * * *--------------------------------------------------------------------* * IMPLICIT INTEGER(I-N) -INC SMLMOTS -INC CCOPTIO * JGN=8 IRET=1 * *--------------------------------------------------------------------* * * * T Y P E P L A S T I Q U E * * * *--------------------------------------------------------------------* IF (IPLAC.EQ.1) THEN * * MODELE DE VON MISES ECROUISSAGE ISOTROPE * IF (MFR.EQ.15.OR.MFR.EQ.17) GOTO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='TRAC' GOTO 9999 ELSE IF (IPLAC.EQ.2) THEN * * MODELE PLASTIQUE PARFAIT * IF (MFR.EQ.15.OR.MFR.EQ.17) GOTO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIGY' GOTO 9999 ELSE IF (IPLAC.EQ.3) THEN * * MODELE DE VON MISES ECROUISSAGE CINEMATIQUE BILIEAIRE * IF (MFR.EQ.15.OR.MFR.EQ.17) GOTO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='H ' GOTO 9999 ELSE IF (IPLAC.EQ.4) THEN * * MODELE PLASTIQUE DRUCKER GENERAL A ECOULEMENT NON ASSOCIER * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='ALFA' MOTS(JGA+2)='BETA' MOTS(JGA+3)='K ' MOTS(JGA+4)='GAMM' MOTS(JGA+5)='DELT' MOTS(JGA+6)='ETA ' MOTS(JGA+7)='MU ' MOTS(JGA+8)='KL ' MOTS(JGA+9)='H ' GOTO 9999 ELSE IF (IPLAC.EQ.5) THEN * * MODELE PLASTIQUE DRUCKER PRAGER * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'LTR ' MOTS(JGA+2)= 'LCS ' GOTO 9999 ELSE IF (IPLAC.EQ.6) THEN * * MODELE PLASTIQUE BETON * TOUS LES PARAMETRES DU MODELE SONT EN FACULTATIF * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+28 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1 )='LTR1' MOTS(JGA+2 )='ETR1' MOTS(JGA+3 )='LTT1' MOTS(JGA+4 )='ETT1' MOTS(JGA+5 )='ERS1' MOTS(JGA+6 )='LTR2' MOTS(JGA+7 )='ETR2' MOTS(JGA+8 )='LTT2' MOTS(JGA+9 )='ETT2' MOTS(JGA+10)='ERS2' MOTS(JGA+11)='LTR3' MOTS(JGA+12)='ETR3' MOTS(JGA+13)='LTT3' MOTS(JGA+14)='ETT3' MOTS(JGA+15)='ERS3' MOTS(JGA+16)='BETR' MOTS(JGA+17)='LCS ' MOTS(JGA+18)='ECS ' MOTS(JGA+19)='LBIC' MOTS(JGA+20)='VF1X' MOTS(JGA+21)='VF1Y' MOTS(JGA+22)='VF1Z' MOTS(JGA+23)='VF2X' MOTS(JGA+24)='VF2Y' MOTS(JGA+25)='VF2Z' MOTS(JGA+26)='VF3X' MOTS(JGA+27)='VF3Y' MOTS(JGA+28)='VF3Z' GOTO 9999 ELSE IF (IPLAC.EQ.7) THEN * * MODELE PLASTIQUE CHABOCHE A UN CENTRE ET ECROUISSAGE ISOTROPE * IF (MFR.EQ.15.OR.MFR.EQ.17) GOTO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='A ' MOTS(JGA+2)='C ' MOTS(JGA+3)='R0 ' MOTS(JGA+4)='PSI ' MOTS(JGA+5)='OMEG' * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOFAC=MLMOTS JGM=JGA SEGINI MLMOT1 DO 10 IC=1,JGA MLMOT1.MOTS(IC)=MOTS(IC) 10 CONTINUE DO 11 IC=1,JGA MOTS(IC+2)=MLMOT1.MOTS(IC) 11 CONTINUE SEGSUP MLMOT1 MOTS(1)='RM ' MOTS(2)='B ' GOTO 9999 ELSE IF (IPLAC.EQ.8) THEN * * MODELE PLASTIQUE CHABOCHE A DEUX CENTRES (ET ECROUISSAGE ISOTRO * IF (MFR.EQ.15.OR.MFR.EQ.17) GOTO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='A1 ' MOTS(JGA+2)='C1 ' MOTS(JGA+3)='A2 ' MOTS(JGA+4)='C2 ' MOTS(JGA+5)='R0 ' MOTS(JGA+6)='PSI ' MOTS(JGA+7)='OMEG' * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOFAC=MLMOTS JGM=JGA SEGINI MLMOT1 DO 20 IC=1,JGA MLMOT1.MOTS(IC)=MOTS(IC) 20 CONTINUE DO 21 IC=1,JGA MOTS(IC+2)=MLMOT1.MOTS(IC) 21 CONTINUE SEGSUP MLMOT1 MOTS(1)='RM ' MOTS(2)='B ' GOTO 9999 ELSE IF (IPLAC.EQ.9) THEN * * MODELE PLASTIQUE PARFAIT ET ECROUISSAGE POUR L'ELEMENT TUYAU FI * IF (MFR.NE.17) GOTO 99 MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOFAC=MLMOTS JGM=JGA SEGINI MLMOT1 DO 30 IC=1,JGA MLMOT1.MOTS(IC)=MOTS(IC) 30 CONTINUE DO 31 IC=1,JGA MOTS(IC+5)=MLMOT1.MOTS(IC) 31 CONTINUE SEGSUP MLMOT1 MOTS(1)='SIGF' MOTS(2)='J1C ' MOTS(3)='TMAT' MOTS(4)='TRAC' MOTS(5)='JDA ' GOTO 9999 ELSE IF (IPLAC.EQ.10) THEN * * ECROUISSAGE ET ENDOMMAGEMENT ISOTROPES - PLAST. ASS.- VON MISES * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='TRAC' MOTS(JGA+2)='EPSD' MOTS(JGA+3)='DC ' MOTS(JGA+4)='EPSR' GOTO 9999 ELSE IF (IPLAC.EQ.11) THEN * * MODELE DE GAUVAIN * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+3 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='TRAC' MOTS(JGA+2)='STOR' MOTS(JGA+3)='SCOM' GOTO 9999 ELSE IF (IPLAC.EQ.12) THEN * * PLASTIQUE PARFAIT ET ECROUISSAGE POUR L'ELEMENT LINESPRING * IF (MFR.NE.15) GOTO 99 MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOFAC=MLMOTS JGM=JGA SEGINI MLMOT1 DO 40 IC=1,JGA MLMOT1.MOTS(IC)=MOTS(IC) 40 CONTINUE DO 41 IC=1,JGA MOTS(IC+2)=MLMOT1.MOTS(IC) 41 CONTINUE SEGSUP MLMOT1 MOTS(1)='SIGF' MOTS(2)='TRAC' GOTO 9999 ELSE IF (IPLAC.EQ.13) THEN * * MODELE UBIQUITOUS * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='NCRI' MOTS(JGA+2)='ANG1' MOTS(JGA+3)='TRA1' MOTS(JGA+4)='PHI1' MOTS(JGA+5)='PSI1' * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1)='ANG2' MOTS(JGA+2)='TRA2' MOTS(JGA+3)='PHI2' MOTS(JGA+4)='PSI2' GOTO 9999 ELSE IF (IPLAC.EQ.14) THEN * * MODELE GLOBAL * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1)='COMP' MOTS(JGA+2)='FLXY' MOTS(JGA+3)='FLXZ' MOTS(JGA+4)='CISY' MOTS(JGA+5)='CISZ' GOTO 9999 ELSE IF (IPLAC.EQ.15) THEN * * MODELE CAM-CLAY * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='E0 ' MOTS(JGA+2)='P0 ' MOTS(JGA+3)='COHE' MOTS(JGA+4)='M ' MOTS(JGA+5)='LAMD' MOTS(JGA+6)='KAPA' MOTS(JGA+7)='G1 ' GOTO 9999 ELSE IF (IPLAC.EQ.16) THEN * * MODELE MOHR-COULOMB POUR LES JOINTS * IF (MFR.NE.35.AND.MFR.NE.75) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) IF(MFR.EQ.35) THEN JGM=JGA+4 ELSE IF (MFR.EQ.75) THEN JGM=JGA+5 ENDIF SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='EF ' MOTS(JGA+2)='ECN ' MOTS(JGA+3)='COHE' MOTS(JGA+4)='FRIC' IF (MFR.EQ.75) MOTS(JGA+5)='TYPE' MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1) ='FTRC' GOTO 9999 ELSE IF (IPLAC.EQ.17) THEN * * MODELE COULOMB_DILATANT POUR LES JOINTS * IF (MFR.NE.35) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+3 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='PHI ' MOTS(JGA+2)='MU ' MOTS(JGA+3)='FTRC' GOTO 9999 ELSE IF (IPLAC.EQ.18) THEN * * MODELE BETON_AXI * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+13 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1) ='LTR1' MOTS(JGA+2) ='LTR2' MOTS(JGA+3) ='LTR3' MOTS(JGA+4) ='ANGL' MOTS(JGA+5) ='LCS ' MOTS(JGA+6) ='ECS ' MOTS(JGA+7) ='LCT ' MOTS(JGA+8) ='PCT ' MOTS(JGA+9) ='LCD ' MOTS(JGA+10)='PCD ' MOTS(JGA+11)='ECD ' MOTS(JGA+12)='LPH ' MOTS(JGA+13)='PENT' GOTO 9999 ELSE IF (IPLAC.EQ.19) THEN * * MODELE DE GURSON POUR LA RUPTURE DUCTILE * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'SIGY' MOTS(JGA+2)= 'H ' MOTS(JGA+3)= 'SBAR' MOTS(JGA+4)= 'PORO' GOTO 9999 * * MODELE BETON LIBRE ET BETON CONFINE (MODELE A FIBRE ET ELEMENTS BARRE) * ELSE IF (IPLAC.EQ.20) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'STFC' MOTS(JGA+2)= 'EZER' MOTS(JGA+3)= 'STFT' MOTS(JGA+4)= 'ALF1' MOTS(JGA+5)= 'OME1' MOTS(JGA+6)= 'ZETA' MOTS(JGA+7)= 'ST85' MOTS(JGA+8)= 'TRAF' MOTS(JGA+9)= 'FACL' MOTS(JGA+10)= 'FAMX' MOTS(JGA+11)= 'STPT' MOTS(JGA+12)= 'FAM1' MOTS(JGA+13)= 'FAM2' GOTO 9999 ELSE IF (IPLAC.EQ.21) THEN * * MODELE ACIER_UNI (MODELE A FIBRE ET ELEMENTS BARRE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+12 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'STSY' MOTS(JGA+ 2)= 'STSU' MOTS(JGA+ 3)= 'EPSH' MOTS(JGA+ 4)= 'EPSU' MOTS(JGA+ 5)= 'ROFA' MOTS(JGA+ 6)= 'BFAC' MOTS(JGA+ 7)= 'A1FA' MOTS(JGA+ 8)= 'A2FA' MOTS(JGA+ 9)= 'FALD' MOTS(JGA+10)= 'A6FA' MOTS(JGA+11)= 'CFAC' MOTS(JGA+12)= 'AFAC' GOTO 9999 * ELSE IF (IPLAC.EQ.51) THEN * * MODELE ACIER_ANCRAGE AVEC LES LOIS ELIGEHAUSEN ET ACIER_UNI * MLMOTS=MOOBL JGA=MOTS(/2) C IF (MFR.EQ.47) THEN C C Modele à fibre C JGM=JGA+21 ELSE C C Elements de barre C JGM=JGA+20 C ENDIF C SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'STSY' MOTS(JGA+ 2)= 'STSU' MOTS(JGA+ 3)= 'EPSH' MOTS(JGA+ 4)= 'EPSU' MOTS(JGA+ 5)= 'ROFA' MOTS(JGA+ 6)= 'BFAC' MOTS(JGA+ 7)= 'A1FA' MOTS(JGA+ 8)= 'A2FA' MOTS(JGA+ 9)= 'FALD' MOTS(JGA+10)= 'A6FA' MOTS(JGA+11)= 'CFAC' MOTS(JGA+12)= 'AFAC' MOTS(JGA+13)= 'G12 ' MOTS(JGA+14)= 'S1T ' MOTS(JGA+15)= 'S2T ' MOTS(JGA+16)= 'S3T ' MOTS(JGA+17)= 'T1T ' MOTS(JGA+18)= 'T3T ' MOTS(JGA+19)= 'ALFA' MOTS(JGA+20)= 'LANC' C IF (MFR.EQ.47) THEN C C Modele à fibre C MOTS(JGA+21)= 'SECB' ENDIF C GOTO 9999 * ELSE IF (IPLAC.EQ.44) THEN * * MODELE FRAGILE_UNI (MODELE A FIBRE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+6 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'FC ' MOTS(JGA+ 2)= 'FC_R' MOTS(JGA+ 3)= 'STRC' MOTS(JGA+ 4)= 'FT ' MOTS(JGA+ 5)= 'FT_R' MOTS(JGA+ 6)= 'STRT' GOTO 9999 * ELSE IF (IPLAC.EQ.45) THEN * * MODELE BETON_BAEL (MODELE A FIBRE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'FC ' GOTO 9999 * ELSE IF (IPLAC.EQ.46) THEN * * MODELE PARFAIT_UNI (MODELE A FIBRE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'SIGY' MOTS(JGA+ 2)= 'H ' GOTO 9999 * ELSE IF (IPLAC.EQ.50) THEN * * MODELE PARFAIT_ANCRAGE AVEC LES LOIS ELIGEHAUSEN ET PARFAIT_UNI * MLMOTS=MOOBL JGA=MOTS(/2) C IF (MFR.EQ.47) THEN C C Modele à fibre C JGM=JGA+11 ELSE C C Elements de barre C JGM=JGA+10 C ENDIF SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'SIGY' MOTS(JGA+ 2)= 'H ' MOTS(JGA+ 3)= 'G12 ' MOTS(JGA+ 4)= 'S1T ' MOTS(JGA+ 5)= 'S2T ' MOTS(JGA+ 6)= 'S3T ' MOTS(JGA+ 7)= 'T1T ' MOTS(JGA+ 8)= 'T3T ' MOTS(JGA+ 9)= 'ALFA' MOTS(JGA+10)= 'LANC' C IF (MFR.EQ.47) THEN C C Modele à fibre C MOTS(JGA+11)= 'SECB' ENDIF C GOTO 9999 C ELSE IF (IPLAC.EQ.53) THEN C C MODELE UNILATERAL (MODELE A FIBRE) C MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'YS1 ' MOTS(JGA+ 2)= 'YS2 ' MOTS(JGA+ 3)= 'A1 ' MOTS(JGA+ 4)= 'A2 ' MOTS(JGA+ 5)= 'B1 ' MOTS(JGA+ 6)= 'B2 ' MOTS(JGA+ 7)= 'BET1' MOTS(JGA+ 8)= 'BET2' MOTS(JGA+ 9)= 'SIGF' GOTO 9999 C ELSE IF (IPLAC.EQ.47) THEN C C MODELE STRUT_UNI (MODELE A FIBRE) C MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+23 SEGADJ MLMOTS MOOBL=MLMOTS * Beton MOTS(JGA+1)= 'STFC' MOTS(JGA+2)= 'EZER' MOTS(JGA+3)= 'STFT' MOTS(JGA+4)= 'ALF1' MOTS(JGA+5)= 'OME1' MOTS(JGA+6)= 'ZETA' MOTS(JGA+7)= 'ST85' MOTS(JGA+8)= 'TRAF' MOTS(JGA+21)= 'STPT' MOTS(JGA+22)= 'FACL' MOTS(JGA+23)= 'FAMX' * Acier MOTS(JGA+ 9)= 'YOUS' MOTS(JGA+10)= 'STSY' MOTS(JGA+11)= 'STSU' MOTS(JGA+12)= 'EPSH' MOTS(JGA+13)= 'EPSU' MOTS(JGA+14)= 'ROFA' MOTS(JGA+15)= 'BFAC' MOTS(JGA+16)= 'A1FA' MOTS(JGA+17)= 'A2FA' * Diagonale MOTS(JGA+18)= 'ROST' MOTS(JGA+19)= 'THET' MOTS(JGA+20)= 'EULT' * GOTO 9999 * * ELSE IF (IPLAC.EQ.22) THEN * * MODELE SECTION * * RIEN DE SPECIAL A FAIRE GOTO 9999 ELSE IF (IPLAC.EQ.23) THEN C C========================================================= C MODELE BETON OTTOSEN C========================================================= C C Formulation DEFORMATIONS PLANES /AXISYMETRIQUE IF ((IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR. & (IFOMOD.EQ.0.OR.IFOMOD.EQ.1)) THEN * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS IF (IFOMOD.EQ.-1.AND.IFOUR.NE.-2) THEN MOTS(JGA+1)='LXX ' MOTS(JGA+2)='LYY ' MOTS(JGA+3)='LXY ' MOTS(JGA+4)='LZZ ' MOTS(JGA+5)='PXX ' MOTS(JGA+6)='PYY ' MOTS(JGA+7)='PXY ' * ELSE IF (IFOMOD.EQ.0.OR.IFOMOD.EQ.1) THEN MOTS(JGA+1)='LRR ' MOTS(JGA+2)='LZZ ' MOTS(JGA+3)='LRZ ' MOTS(JGA+4)='LOO ' MOTS(JGA+5)='PRR ' MOTS(JGA+6)='PZZ ' MOTS(JGA+7)='PRZ ' ENDIF * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+30 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1) ='GFTR' MOTS(JGA+2) ='LTR ' MOTS(JGA+3) ='GFCS' MOTS(JGA+4) ='LCS ' MOTS(JGA+5) ='GS ' MOTS(JGA+6) ='BTR ' MOTS(JGA+7) ='BCS ' MOTS(JGA+8) ='WRUP' MOTS(JGA+9) ='BILI' MOTS(JGA+10)='EPSR' MOTS(JGA+11)='GFT1' MOTS(JGA+12)='GFT2' MOTS(JGA+13)='GFT3' MOTS(JGA+14)='GS1 ' MOTS(JGA+15)='GS2 ' MOTS(JGA+16)='GS3 ' MOTS(JGA+17)='LTR1' MOTS(JGA+18)='LTR2' MOTS(JGA+19)='LTR3' MOTS(JGA+20)='EPCM' MOTS(JGA+21)='EPCU' MOTS(JGA+22)='LCBI' MOTS(JGA+23)='XK2 ' MOTS(JGA+24)='XGB ' MOTS(JGA+25)='XPA ' MOTS(JGA+26)='WRU1' MOTS(JGA+27)='WRU2' MOTS(JGA+28)='BIL1' MOTS(JGA+29)='BIL2' MOTS(JGA+30)='EPSB' GOTO 9999 ENDIF C C Formulation CONTRAINTES PLANES / 3D COQUES MINCES IF ((IFOMOD.EQ.-1.AND.IFOUR.EQ.-2).OR. & (IFOMOD.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9))) THEN * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS IF (IFOMOD.EQ.-1.AND.IFOUR.EQ.-2) THEN MOTS(JGA+1)='LXX ' MOTS(JGA+2)='LYY ' MOTS(JGA+3)='LXY ' MOTS(JGA+4)='LZZ ' MOTS(JGA+5)='PXX ' MOTS(JGA+6)='PYY ' MOTS(JGA+7)='PXY ' * ELSE IF (IFOMOD.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9)) THEN MOTS(JGA+1)='LSS ' MOTS(JGA+2)='LTT ' MOTS(JGA+3)='LST ' MOTS(JGA+4)='LNN ' MOTS(JGA+5)='PSS ' MOTS(JGA+6)='PTT ' MOTS(JGA+7)='PST ' ENDIF * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+25 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1) ='GFTR' MOTS(JGA+2) ='LTR ' MOTS(JGA+3) ='GFCS' MOTS(JGA+4) ='LCS ' MOTS(JGA+5) ='GS ' MOTS(JGA+6) ='BTR ' MOTS(JGA+7) ='BCS ' MOTS(JGA+8) ='WRUP' MOTS(JGA+9) ='BILI' MOTS(JGA+10)='GFT1' MOTS(JGA+11)='GFT2' MOTS(JGA+12)='GS1 ' MOTS(JGA+13)='GS2 ' MOTS(JGA+14)='LTR1' MOTS(JGA+15)='LTR2' MOTS(JGA+16)='EPCM' MOTS(JGA+17)='EPCU' MOTS(JGA+18)='LCBI' MOTS(JGA+19)='XK2 ' MOTS(JGA+20)='XGB ' MOTS(JGA+21)='XPA ' MOTS(JGA+22)='WRU1' MOTS(JGA+23)='WRU2' MOTS(JGA+24)='BIL1' MOTS(JGA+25)='BIL2' GOTO 9999 ENDIF C C Formulation TRIDIMENSIONNELLE MASSIVE IF (IFOMOD.EQ.2.AND.MFR.EQ.1) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+12 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1) ='LXX ' MOTS(JGA+2) ='LYY ' MOTS(JGA+3) ='LZZ ' MOTS(JGA+4) ='LXY ' MOTS(JGA+5) ='LXZ ' MOTS(JGA+6) ='LYZ ' MOTS(JGA+7) ='PXX ' MOTS(JGA+8) ='PYY ' MOTS(JGA+9) ='PZZ ' MOTS(JGA+10)='PXY ' MOTS(JGA+11)='PXZ ' MOTS(JGA+12)='PYZ ' MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+30 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1)='GFTR' MOTS(JGA+2)='LTR ' MOTS(JGA+3)='GFCS' MOTS(JGA+4)='LCS ' MOTS(JGA+5)='GS ' MOTS(JGA+6)='BTR ' MOTS(JGA+7)='BCS ' MOTS(JGA+8)='WRUP' MOTS(JGA+9)='BILI' MOTS(JGA+10)='GFT1' MOTS(JGA+11)='GFT2' MOTS(JGA+12)='GFT3' MOTS(JGA+13)='GS1 ' MOTS(JGA+14)='GS2 ' MOTS(JGA+15)='GS3 ' MOTS(JGA+16)='LTR1' MOTS(JGA+17)='LTR2' MOTS(JGA+18)='LTR3' MOTS(JGA+19)='EPCM' MOTS(JGA+20)='EPCU' MOTS(JGA+21)='LCBI' MOTS(JGA+22)='XK2 ' MOTS(JGA+23)='XGB ' MOTS(JGA+24)='XPA ' MOTS(JGA+25)='WRU1' MOTS(JGA+26)='WRU2' MOTS(JGA+27)='WRU3' MOTS(JGA+28)='BIL1' MOTS(JGA+29)='BIL2' MOTS(JGA+30)='BIL3' GOTO 9999 ENDIF * ELSE IF (IPLAC.EQ.24) THEN * * MODELE MAZARS_FIB (MODELE A FIBRE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+6 SEGADJ MLMOTS MOTS(JGA+1)= 'KTR0' MOTS(JGA+2)= 'ACOM' MOTS(JGA+3)= 'BCOM' MOTS(JGA+4)= 'ATRA' MOTS(JGA+5)= 'BTRA' MOTS(JGA+6)= 'BETA' GOTO 9999 * ELSE IF (IPLAC.EQ.25) THEN * * MODELE AMADEI-SAEB POUR LES JOINTS * IF (MFR.NE.35) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+10 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='FIMU' MOTS(JGA+2)='SGMT' MOTS(JGA+3)='I0 ' MOTS(JGA+4)='S0 ' MOTS(JGA+5)='B0 ' MOTS(JGA+6)='UR ' MOTS(JGA+7)='UP ' MOTS(JGA+8)='KNI ' MOTS(JGA+9)='FI0 ' MOTS(JGA+10)='VM ' GOTO 9999 * ELSE IF (IPLAC.EQ.26) THEN * * MODELE HUJEUX * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+10 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='P0 ' MOTS(JGA+2)='COHE' MOTS(JGA+3)='M ' MOTS(JGA+4)='BETA' MOTS(JGA+5)='A ' MOTS(JGA+6)='B ' MOTS(JGA+7)='R0 ' MOTS(JGA+8)='P1 ' MOTS(JGA+9)='N ' MOTS(JGA+10)='E1 ' GOTO 9999 * ELSE IF (IPLAC.EQ.27) THEN * * MODELE STEINBERG COCHRAN GUINAN * * IF (MFR.EQ.15.OR.MFR.EQ.17) GOTO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='BETA' MOTS(JGA+3)='N' MOTS(JGA+4)='EPSI' MOTS(JGA+5)='GP' MOTS(JGA+6)='GT' MOTS(JGA+7)='YMAX' MOTS(JGA+8)='TMO' MOTS(JGA+9)='MU' GOTO 9999 * ELSE IF (IPLAC.EQ.28) THEN * * MODELE ZERILLI ARMSTRONG * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+10 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='DYG' MOTS(JGA+2)='C1' MOTS(JGA+3)='C2' MOTS(JGA+4)='C3' MOTS(JGA+5)='C4' MOTS(JGA+6)='C5' MOTS(JGA+7)='N' MOTS(JGA+8)='K' MOTS(JGA+9)='L' MOTS(JGA+10)='TYPE' GOTO 9999 * ELSE IF (IPLAC.EQ.29) THEN * * MODELE PRESTON TONKS WALLACE * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+15 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='TAU' MOTS(JGA+2)='P' MOTS(JGA+3)='S0' MOTS(JGA+4)='SINF' MOTS(JGA+5)='K' MOTS(JGA+6)='G' MOTS(JGA+7)='Y0' MOTS(JGA+8)='YINF' MOTS(JGA+9)='Y1' MOTS(JGA+10)='Y2' MOTS(JGA+11)='BETA' MOTS(JGA+12)='GP' MOTS(JGA+13)='GT' MOTS(JGA+14)='MU' MOTS(JGA+15)='TMO' GOTO 9999 * ELSE IF (IPLAC.EQ.30) THEN * * MODELE BETOCYCL EN CONTRAINTE PLANE * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+11 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='HHH1' MOTS(JGA+2)='FTPE' MOTS(JGA+3)='FCPE' MOTS(JGA+4)='FTGR' MOTS(JGA+5)='FCGR' MOTS(JGA+6)='EPSO' MOTS(JGA+7)='WOR0' MOTS(JGA+8)='TREV' MOTS(JGA+9)='COEV' MOTS(JGA+10)='LCAT' MOTS(JGA+11)='LCAC' GOTO 9999 * ELSE IF (IPLAC.EQ.31) THEN C C MODELE ROTATING CRACK EN CONTRAINTE PLANE C MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'FTRA' MOTS(JGA+2)= 'EPSR' MOTS(JGA+3)= 'FRES' MOTS(JGA+4)= 'BETA' GOTO 9999 * ELSE IF (IPLAC.EQ.32) THEN * * MODELE JOINT_SOFT POUR LES JOINTS * IF (MFR.NE.35) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) * JGM=JGA+6 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='PNOR' MOTS(JGA+2)='CPLG' MOTS(JGA+3)='BETA' MOTS(JGA+4)='SJCB' MOTS(JGA+5)='SJSB' MOTS(JGA+6)='SJTB' * GOTO 9999 * ELSE IF (IPLAC.EQ.49) THEN * * MODELE ANCRAGE_ELIGEHAUSEN POUR LES JOINTS * IF (MFR.NE.35) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='S1T ' MOTS(JGA+2)='S2T ' MOTS(JGA+3)='S3T ' MOTS(JGA+4)='T1T ' MOTS(JGA+5)='T3T ' MOTS(JGA+6)='ALFA' MOTS(JGA+7)='PERI' GOTO 9999 * ELSE IF (IPLAC.EQ.33) THEN * * MODELE BILINEAIRE MOMY * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOTS(JGA+1)='EAYI' MOTS(JGA+2)='YMOM' GOTO 9999 * ELSE IF (IPLAC.EQ.34) THEN * * MODELE BILINEAIRE EFFZ * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOTS(JGA+1)='EAYI' MOTS(JGA+2)='YEFF' GOTO 9999 * ELSE IF ((IPLAC.EQ.35).OR.(IPLAC.EQ.36)) THEN * * MODELE PLASTIQUE TAKEMO_MOMY OU TAKEMO_EFFZ * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOTS(JGA+1)='TRAC' MOTS(JGA+2)='SFDP' MOTS(JGA+3)='SFDN' MOTS(JGA+4)='PINP' MOTS(JGA+5)='PINN' MOTS(JGA+6)='SRDP' MOTS(JGA+7)='SRDN' GOTO 9999 C ELSE IF(IPLAC.EQ.37) THEN * * MODELE BETON_INSA * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+12 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1 )='ALFA' MOTS(JGA+2 )='LCS ' MOTS(JGA+3 )='EMAX' MOTS(JGA+4 )='EPUT' MOTS(JGA+5 )='FTC ' MOTS(JGA+6 )='ICOM' MOTS(JGA+7 )='GFTR' * MOTS(JGA+8 )='IFLU' MOTS(JGA+9 )='RMOY' MOTS(JGA+10)='PHIF' MOTS(JGA+11)='TEM1' MOTS(JGA+12)='POAR' GOTO 9999 * ELSE IF (IPLAC.EQ.38.AND.(MFR.EQ.27.OR.MFR.EQ.49)) THEN * * ECROUIS_INSA (MATERIAU ECROUISSABLE UNIDIRECTIONEL) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='TRAC' GOTO 9999 * ELSE IF (IPLAC.EQ.38.AND.MFR.NE.27.AND.MFR.NE.49) THEN * * ECROUIS_INSA (MATERIAU ORTHOTROPE ECROUISSABLE DECOUPLE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='TRA1' MOTS(JGA+2)='TRA2' GOTO 9999 * ELSE IF (IPLAC.EQ.39.AND.(MFR.EQ.27.OR.MFR.EQ.49)) THEN * * PARFAIT_INSA (MATERIAU PLASTIQUE PARFAIT UNIDIRECTIONEL) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIGY' GOTO 9999 * ELSE IF (IPLAC.EQ.39.AND.MFR.NE.27.AND.MFR.NE.49) THEN * * PARFAIT_INSA (MATERIAU ORTHOTROPE PLASTIQUE PARFAIT DECOUPLE) * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIG1' MOTS(JGA+2)='SIG2' GOTO 9999 * ELSE IF (IPLAC.EQ.40) THEN * * ALONSO * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+12 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='KS' MOTS(JGA+2)='LAMS' MOTS(JGA+3)='M' MOTS(JGA+4)='KK' MOTS(JGA+5)='PC0' MOTS(JGA+6)='P0' MOTS(JGA+7)='LAM1' MOTS(JGA+8)='T0' MOTS(JGA+9)='TAU' MOTS(JGA+10)='GTAN' MOTS(JGA+11)='KTAN' MOTS(JGA+12)='E0' * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1 )='KALF' MOTS(JGA+2 )='KBET' MOTS(JGA+3 )='KGAM' MOTS(JGA+4 )='PREF' GOTO 9999 * ELSE IF (IPLAC.EQ.41) THEN * * PAKZAD * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+15 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='N' MOTS(JGA+2)='KA' MOTS(JGA+3)='GA' MOTS(JGA+4)='PATM' MOTS(JGA+5)='PC0' MOTS(JGA+6)='M' MOTS(JGA+7)='BETA' MOTS(JGA+8)='A0' MOTS(JGA+9)='PREF' MOTS(JGA+10)='S0' MOTS(JGA+11)='M1' MOTS(JGA+12)='M2' MOTS(JGA+13)='M3' MOTS(JGA+14)='BSEC' MOTS(JGA+15)='TAU' GOTO 9999 * ELSE IF (IPLAC.EQ.42) THEN * * MODELE PLASTIQUE INFILL_UNI * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+8 SEGADJ MLMOTS MOTS(JGA+1)='DELA' MOTS(JGA+2)='DMAX' MOTS(JGA+3)='BETA' MOTS(JGA+4)='GAMM' MOTS(JGA+5)='GAMP' MOTS(JGA+6)='ALFA' MOTS(JGA+7)='TETA' * MOTS(JGA+8)='NETA' MOTS(JGA+8)='MONO' GOTO 9999 * ELSE IF (IPLAC.EQ.43) THEN * * MODELE PLASTIQUE CISAIL_NL * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ MLMOTS MOTS(JGA+1)='DELP' MOTS(JGA+2)='DELN' MOTS(JGA+3)='DMAP' MOTS(JGA+4)='DMAN' MOTS(JGA+5)='BETA' MOTS(JGA+6)='ALFA' MOTS(JGA+7)='TETA' MOTS(JGA+8)='MONP' MOTS(JGA+9)='MONN' GOTO 9999 ELSE IF (IPLAC.EQ.48) THEN * *+++++++++++++++++++++++++++++++++ * MODELE HINTE *+++++++++++++++++++++++++++++++++++ * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+6 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='GAM1' MOTS(JGA+2)='AL ' MOTS(JGA+3)='DCRI' MOTS(JGA+4)='NN ' MOTS(JGA+5)='YO ' MOTS(JGA+6)='YC ' GOTO 9999 ELSE IF (IPLAC.EQ.52) THEN * * MODELE PLASTIQUE BETON_URGC * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ MLMOTS MOTS(JGA+1)='ALFA' MOTS(JGA+2)='BETA' MOTS(JGA+3)='LCS ' MOTS(JGA+4)='GFC' MOTS(JGA+5)='GFT' MOTS(JGA+6)='ITON' MOTS(JGA+7)='IMOD' MOTS(JGA+8)='ICT ' MOTS(JGA+9)='ICC' GOTO 9999 * BCN ELSE IF (IPLAC.EQ.54) THEN CCC CCC MRS-Lade model CCC MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+18 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='PC' MOTS(JGA+2)='PA' MOTS(JGA+3)='QA' MOTS(JGA+4)='EXPM' MOTS(JGA+5)='E' MOTS(JGA+6)='K1' MOTS(JGA+7)='K2' MOTS(JGA+8)='ETAB' MOTS(JGA+9)='EXPV' MOTS(JGA+10)='EPSI' MOTS(JGA+11)='N' MOTS(JGA+12)='CCON' MOTS(JGA+13)='EXPL' MOTS(JGA+14)='PCAP' MOTS(JGA+15)='EXPR' MOTS(JGA+16)='CCAP' MOTS(JGA+17)='PHI' MOTS(JGA+18)='ALP' GOTO 9999 ELSE IF (IPLAC.EQ.55) THEN CCC CCC J2 model CCC MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIG0' MOTS(JGA+2)='SIGI' MOTS(JGA+3)='KISO' MOTS(JGA+4)='VELO' GOTO 9999 ELSE IF ((IPLAC.EQ.56).or.(IPLAC.EQ.58)) THEN CCC CCC Rounded Hyperbolic Mohr_coulomb model CCC MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='COHE' MOTS(JGA+2)='PHI' GOTO 9999 ELSE IF (IPLAC.EQ.57) THEN CCC CCC VMT_FEFP model CCC MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIG0' MOTS(JGA+2)='SIGI' MOTS(JGA+3)='KISO' MOTS(JGA+4)='VELO' MOTS(JGA+5)='MSHA' GOTO 9999 ELSE IF (IPLAC.EQ.59) THEN CCC CCC POWDER_FEFP model CCC MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='NNN1' MOTS(JGA+3)='NNN2' MOTS(JGA+4)='ETA0' GOTO 9999 ELSE IF (IPLAC.EQ.60) THEN CCC CCC POWDERCAP_FEFP model CCC MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+8 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='COHE' MOTS(JGA+2)='PHI' MOTS(JGA+3)='ETA0' MOTS(JGA+4)='PHI0' MOTS(JGA+5)='NNNC' MOTS(JGA+6)='SIGY' MOTS(JGA+7)='NNN1' MOTS(JGA+8)='NNN2' GOTO 9999 * BCN ELSE IF (IPLAC.EQ.61) THEN * * MODELE JOINT_COAT POUR LES JOINTS CISAILLEMENT * IF (MFR.NE.53) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='SJSB' MOTS(JGA+2)='BETA' GOTO 9999 * ELSE IF (IPLAC.EQ.62) THEN * * MODELE JOINT_COAT POUR LES JOINTS CISAILLEMENT * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='SIT0' MOTS(JGA+2 )='SIFC' MOTS(JGA+3 )='SIGT' MOTS(JGA+4 )='EPTU' MOTS(JGA+5 )='SIGC' MOTS(JGA+6 )='EPSC' MOTS(JGA+7 )='SICU' MOTS(JGA+8 )='EPCU' MOTS(JGA+9 )='COED' MOTS(JGA+10)='COER' MOTS(JGA+11)='BETA' MOTS(JGA+12)='BET1' MOTS(JGA+13)='COEU' GOTO 9999 C+PPf C+PPm * ELSE IF (IPLAC.EQ.63) THEN * * MODELE MUR_SHEAR POUR LES JOINTS CISAILLEMENT * IF(MFR.NE.7.AND.MFR.NE.61) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='CUFP' MOTS(JGA+2 )='CUKP' MOTS(JGA+3 )='CULP' MOTS(JGA+4 )='CUFM' MOTS(JGA+5 )='CUKM' MOTS(JGA+6 )='CULM' MOTS(JGA+7 )='DTRP' MOTS(JGA+8 )='DTRM' MOTS(JGA+9 )='BETA' MOTS(JGA+10)='NELA' MOTS(JGA+11)='TRFA' MOTS(JGA+12)='DOCP' MOTS(JGA+13)='DOCM' * C+2004 !!!!!!!!!! RHO et ALPH doivent etre a la fin! voir defmat!!! MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+3 SEGADJ MLMOTS MOFAC=MLMOTS DO ie1=1,JGA MOTS(JGM+1-ie1)=MOTS(JGA+1-ie1) ENDDO MOTS(1)='EXPN' MOTS(2)='FNRF' MOTS(3)='FTDI' C+2004 GOTO 9999 * ELSE IF (IPLAC.EQ.64) THEN * * MODELE INCO * IF (MFR.NE.35) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='LTR ' MOTS(JGA+2)='GFTR' MOTS(JGA+3)='HS ' MOTS(JGA+4)='DINT' MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1) ='FTRC' GOTO 9999 * modif DP2 ELSE IF (IPLAC.EQ.65) THEN * * MODELE PLASTIQUE DRUCKER PRAGER A DEUX SURFACES * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='LTR' MOTS(JGA+2)='LCS' MOTS(JGA+3)='LBI' MOTS(JGA+4)='SIGY' MOTS(JGA+5)='EPM' MOTS(JGA+6)='EPU' MOTS(JGA+7)='LCAR' MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+1 SEGADJ MLMOTS MOFAC=MLMOTS JGM=JGA SEGINI MLMOT1 DO IC=1,JGA MLMOT1.MOTS(IC)=MOTS(IC) ENDDO DO IC=1,JGA MOTS(IC+1)=MLMOT1.MOTS(IC) ENDDO SEGSUP MLMOT1 MOTS(1)='GF ' GOTO 9999 ELSE IF (IPLAC.EQ.69) THEN * * MODELE PLASTIQUE BILIN_EFFX * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOTS(JGA+1)='ETAN' MOTS(JGA+2)='YIEF' GOTO 9999 ELSE IF (IPLAC.EQ.67) THEN * * MODELE PLASTIQUE INTIMP * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+14 SEGADJ MLMOTS MOTS(JGA+ 1)= 'SOCT' MOTS(JGA+ 2)= 'DCS ' MOTS(JGA+ 3)= 'SOGS' MOTS(JGA+ 4)= 'TCS ' MOTS(JGA+ 5)= 'MS ' MOTS(JGA+ 6)= 'KS ' MOTS(JGA+ 7)= 'GCEO' MOTS(JGA+ 8)= 'AD ' MOTS(JGA+ 9)= 'ACOE' MOTS(JGA+10)= 'GAMC' MOTS(JGA+11)= 'LCCO' MOTS(JGA+12)= 'EPSC' MOTS(JGA+13)= 'TCI ' MOTS(JGA+14)= 'CALA' GOTO 9999 ELSE IF (IPLAC.EQ.68) THEN * * MODELE PLASTIQUE INTJOI2D * IF (MFR.NE.35) GO TO 99 MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='AD ' MOTS(JGA+2 )='Y0 ' MOTS(JGA+3 )='ALPA' MOTS(JGA+4 )='GAIN' MOTS(JGA+5 )='AAIN' MOTS(JGA+6 )='Q1CO' MOTS(JGA+7 )='Q2CO' MOTS(JGA+8 )='Q3CO' MOTS(JGA+9 )='SYCO' MOTS(JGA+10)='NCOE' MOTS(JGA+11)='KCOE' MOTS(JGA+12)='TC ' MOTS(JGA+13)='GONF' GOTO 9999 * ELSE IF (IPLAC.EQ.70) THEN * * MODELE PLASTIQUE ISS_GRANGE * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+15 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='DIAM ' MOTS(JGA+2 )='LX ' MOTS(JGA+3 )='LY ' C MOTS(JGA+4 )='KELZ ' C MOTS(JGA+5 )='KELH ' C MOTS(JGA+6 )='KERY ' C MOTS(JGA+7 )='KERX ' C MOTS(JGA+8 )='KELT ' MOTS(JGA+4 )='XA ' MOTS(JGA+5 )='XB ' MOTS(JGA+6 )='XC ' MOTS(JGA+7 )='XD ' MOTS(JGA+8 )='XE ' MOTS(JGA+9 )='XF ' MOTS(JGA+10)='QMAX ' MOTS(JGA+11)='A9 ' MOTS(JGA+12)='A6 ' MOTS(JGA+13)='ETA3 ' MOTS(JGA+14)='XTIM ' MOTS(JGA+15)='A8 ' GOTO 9999 * ELSE IF (IPLAC.EQ.71) THEN * * MODELE PLASTIQUE RUP_THER * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+10 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='KB ' MOTS(JGA+2 )='YBO ' MOTS(JGA+3 )='C1 ' MOTS(JGA+4 )='D1 ' MOTS(JGA+5 )='ALPB' MOTS(JGA+6 )='BETB' MOTS(JGA+7 )='KA ' MOTS(JGA+8 )='YAO ' MOTS(JGA+9 )='ALPA' MOTS(JGA+10 )='BETA' GOTO 9999 C +BR ELSE IF (IPLAC.EQ.72) THEN * * MODELE PLASTIQUE RICBET_UNI * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+12 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='FT ' MOTS(JGA+2 )='ALDI' MOTS(JGA+3 )='GAM1' MOTS(JGA+4 )='A1 ' MOTS(JGA+5 )='SIGF' MOTS(JGA+6 )='FC ' MOTS(JGA+7 )='AF ' MOTS(JGA+8 )='AG ' MOTS(JGA+9 )='AC ' MOTS(JGA+10)='BC ' MOTS(JGA+11)='SIGU' MOTS(JGA+12)='HYST' GOTO 9999 ELSE IF (IPLAC.EQ.73) THEN * * MODELE RAGMOUT * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1 )='UELA' MOTS(JGA+2 )='FPLA' MOTS(JGA+3 )='PFIS' MOTS(JGA+4 )='HCIN' MOTS(JGA+5 )='QFRA' MOTS(JGA+6 )='APIH' MOTS(JGA+7 )='BPIH' GOTO 9999 C -BR * * ENDIF GOTO 99 * 9999 CONTINUE IRET=1 RETURN * * SI PROBLEME , ON SORT AVEC IRET A 0 * 99 CONTINUE IRET=0 CALL ERREUR(261) RETURN END