C IDVISC SOURCE YURI 11/10/28 22:15:28 7181 C======================================================================= C= Noms de composantes de materiaux viscoplastiques = C= = C= ENTREES : = C= --------- = C= MFR Numero de formulation = C= MOOBL Pointeur sur les mots OBLIGATOIRES = C= MOFAC Pointeur sur les mots FACULTATIFS = C= IPLAC Numero du materiau dans MODVIS = C= = C= SORTIE : = C= -------- = C= IRET =0 si probleme = C======================================================================= SUBROUTINE IDVISC (MFR,MOOBL,MOFAC,IPLAC,IRET) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) -INC CCOPTIO -INC SMLMOTS JGN=8 IRET=1 C === C 1 - Modele viscoplastique de GUIONNET C === IF (IPLAC.EQ.1) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+20 SEGADJ,MLMOTS MOTS(JGA+ 1)='N ' MOTS(JGA+ 2)='K ' MOTS(JGA+ 3)='CK ' MOTS(JGA+ 4)='R0 ' MOTS(JGA+ 5)='CD ' MOTS(JGA+ 6)='M ' MOTS(JGA+ 7)='A ' MOTS(JGA+ 8)='A1 ' MOTS(JGA+ 9)='C ' MOTS(JGA+10)='C1 ' MOTS(JGA+11)='G ' MOTS(JGA+12)='R ' MOTS(JGA+13)='NN ' MOTS(JGA+14)='C2 ' MOTS(JGA+15)='Q ' MOTS(JGA+16)='G1 ' MOTS(JGA+17)='R1 ' MOTS(JGA+18)='BETA' MOTS(JGA+19)='C0 ' MOTS(JGA+20)='P1M0' MOOBL=MLMOTS RETURN C === C 2 - Modele viscoplastique de CHABOCHE C === ELSE IF (IPLAC.EQ.2) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+22 SEGADJ,MLMOTS MOTS(JGA+ 1)='ALF ' MOTS(JGA+ 2)='N ' MOTS(JGA+ 3)='KK ' MOTS(JGA+ 4)='K0 ' MOTS(JGA+ 5)='ALFK' MOTS(JGA+ 6)='ALFR' MOTS(JGA+ 7)='C1 ' MOTS(JGA+ 8)='A1 ' MOTS(JGA+ 9)='PHI ' MOTS(JGA+10)='C2 ' MOTS(JGA+11)='A2 ' MOTS(JGA+12)='B ' MOTS(JGA+13)='BET1' MOTS(JGA+14)='BET2' MOTS(JGA+15)='R1 ' MOTS(JGA+16)='R2 ' MOTS(JGA+17)='GAMA' MOTS(JGA+18)='M ' MOTS(JGA+19)='ETA ' MOTS(JGA+20)='MU ' MOTS(JGA+21)='QMAX' MOTS(JGA+22)='QSTA' MOOBL=MLMOTS RETURN C === C 3 - Modele d'ecrouissage-VISCOsite multiplicative avec enDOMMAGEment C === ELSE IF (IPLAC.EQ.3) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+8 SEGADJ,MLMOTS MOTS(JGA+1)='N ' MOTS(JGA+2)='M ' MOTS(JGA+3)='KK ' MOTS(JGA+4)='ALP1' MOTS(JGA+5)='BLP1' MOTS(JGA+6)='R ' MOTS(JGA+7)='A ' MOTS(JGA+8)='EVOL' MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+1 SEGADJ,MLMOTS MOTS(JGA+1)='SMAX' MOFAC=MLMOTS RETURN C === C 4 - Modele viscoplastique PARFAIT C === ELSE IF (IPLAC.EQ.4) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+3 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='N ' MOTS(JGA+3)='K ' MOOBL=MLMOTS RETURN C === C 5 - Modele viscoplastique POUDRE_A C === ELSE IF (IPLAC.EQ.5) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+16 SEGADJ,MLMOTS MOTS(JGA+ 1)='N ' MOTS(JGA+ 2)='A ' MOTS(JGA+ 3)='QSRT' MOTS(JGA+ 4)='F0 ' MOTS(JGA+ 5)='F1 ' MOTS(JGA+ 6)='F2 ' MOTS(JGA+ 7)='F3 ' MOTS(JGA+ 8)='F4 ' MOTS(JGA+ 9)='F5 ' MOTS(JGA+10)='C0 ' MOTS(JGA+11)='C1 ' MOTS(JGA+12)='C2 ' MOTS(JGA+13)='C3 ' MOTS(JGA+14)='C4 ' MOTS(JGA+15)='C5 ' MOTS(JGA+16)='RHOR' MOOBL=MLMOTS RETURN C === C 6 - Modele viscoplastique POUDRE_B C === ELSE IF (IPLAC.EQ.6) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+23 SEGADJ,MLMOTS MOTS(JGA+ 1)='N ' MOTS(JGA+ 2)='A ' MOTS(JGA+ 3)='QSRT' MOTS(JGA+ 4)='F0 ' MOTS(JGA+ 5)='F1 ' MOTS(JGA+ 6)='F2 ' MOTS(JGA+ 7)='F3 ' MOTS(JGA+ 8)='F4 ' MOTS(JGA+ 9)='F5 ' MOTS(JGA+10)='C0 ' MOTS(JGA+11)='C1 ' MOTS(JGA+12)='C2 ' MOTS(JGA+13)='C3 ' MOTS(JGA+14)='C4 ' MOTS(JGA+15)='C5 ' MOTS(JGA+16)='G0 ' MOTS(JGA+17)='G1 ' MOTS(JGA+18)='G2 ' MOTS(JGA+19)='G3 ' MOTS(JGA+20)='G4 ' MOTS(JGA+21)='G5 ' MOTS(JGA+22)='G6 ' MOTS(JGA+23)='RHOR' MOOBL=MLMOTS RETURN C === C 7 - Modele viscoplastique de OHNO C === ELSE IF (IPLAC.EQ.7) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+24 SEGADJ,MLMOTS MOTS(JGA+ 1)='ALF ' MOTS(JGA+ 2)='N ' MOTS(JGA+ 3)='KK ' MOTS(JGA+ 4)='K0 ' MOTS(JGA+ 5)='ALFK' MOTS(JGA+ 6)='ALFR' MOTS(JGA+ 7)='C1 ' MOTS(JGA+ 8)='A1 ' MOTS(JGA+ 9)='PHI ' MOTS(JGA+10)='C2 ' MOTS(JGA+11)='A2 ' MOTS(JGA+12)='B ' MOTS(JGA+13)='BET1' MOTS(JGA+14)='BET2' MOTS(JGA+15)='R1 ' MOTS(JGA+16)='R2 ' MOTS(JGA+17)='GAMA' MOTS(JGA+18)='M ' MOTS(JGA+19)='ETA ' MOTS(JGA+20)='MU ' MOTS(JGA+21)='QMAX' MOTS(JGA+22)='QSTA' MOTS(JGA+23)='EXP1' MOTS(JGA+24)='EXP2' MOOBL=MLMOTS RETURN C === C 8 - Modele viscoplastique DDI C === ELSE IF (IPLAC.EQ.8) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+18 SEGADJ,MLMOTS MOTS(JGA+ 1)='RP0 ' MOTS(JGA+ 2)='QP ' MOTS(JGA+ 3)='BP ' MOTS(JGA+ 4)='CP1 ' MOTS(JGA+ 5)='CP2 ' MOTS(JGA+ 6)='DP1 ' MOTS(JGA+ 7)='DP2 ' MOTS(JGA+ 8)='KS ' MOTS(JGA+ 9)='N ' MOTS(JGA+10)='RV0 ' MOTS(JGA+11)='QV ' MOTS(JGA+12)='BV ' MOTS(JGA+13)='CV1 ' MOTS(JGA+14)='CV2 ' MOTS(JGA+15)='DV1 ' MOTS(JGA+16)='DV2 ' MOTS(JGA+17)='CVP1' MOTS(JGA+18)='CVP2' MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+10 SEGADJ,MLMOTS MOTS(JGA+ 1)='ALP7' MOTS(JGA+ 2)='BETA' MOTS(JGA+ 3)='M ' MOTS(JGA+ 4)='SU ' MOTS(JGA+ 5)='SL ' MOTS(JGA+ 6)='B ' MOTS(JGA+ 7)='C ' MOTS(JGA+ 8)='A ' MOTS(JGA+ 9)='R ' MOTS(JGA+10)='KK ' MOFAC=MLMOTS RETURN C === C 9 - Modele viscoplastique de KOCK C === ELSE IF (IPLAC.EQ.9) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+10 SEGADJ,MLMOTS MOTS(JGA+ 1)='A ' MOTS(JGA+ 2)='B ' MOTS(JGA+ 3)='M ' MOTS(JGA+ 4)='Q ' MOTS(JGA+ 5)='R ' MOTS(JGA+ 6)='H0 ' MOTS(JGA+ 7)='AP ' MOTS(JGA+ 8)='SB ' MOTS(JGA+ 9)='N ' MOTS(JGA+10)='S0 ' MOOBL=MLMOTS RETURN C ==== C 10 - Modele viscoplastique NOUAILLAS_A C ==== ELSE IF (IPLAC.EQ.10) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+22 SEGADJ,MLMOTS MOTS(JGA+ 1)='ALF ' MOTS(JGA+ 2)='N ' MOTS(JGA+ 3)='KK ' MOTS(JGA+ 4)='K0 ' MOTS(JGA+ 5)='ALFK' MOTS(JGA+ 6)='ALFR' MOTS(JGA+ 7)='CL1 ' MOTS(JGA+ 8)='DNL1' MOTS(JGA+ 9)='PHI ' MOTS(JGA+10)='CL2 ' MOTS(JGA+11)='DNL2' MOTS(JGA+12)='B ' MOTS(JGA+13)='GDM1' MOTS(JGA+14)='GDM2' MOTS(JGA+15)='PTM1' MOTS(JGA+16)='PTM2' MOTS(JGA+17)='GAMA' MOTS(JGA+18)='M ' MOTS(JGA+19)='ETA ' MOTS(JGA+20)='MU ' MOTS(JGA+21)='QMAX' MOTS(JGA+22)='QSTA' MOOBL=MLMOTS RETURN C ==== C 11 - Modele viscoplastique NOUAILLAS_B C ==== ELSE IF (IPLAC.EQ.11) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+14 SEGADJ,MLMOTS MOTS(JGA+ 1)='ALF ' MOTS(JGA+ 2)='N ' MOTS(JGA+ 3)='KK ' MOTS(JGA+ 4)='K0 ' MOTS(JGA+ 5)='CL1 ' MOTS(JGA+ 6)='DNL1' MOTS(JGA+ 7)='PHI ' MOTS(JGA+ 8)='CL2 ' MOTS(JGA+ 9)='DNL2' MOTS(JGA+10)='B ' MOTS(JGA+11)='GDM1' MOTS(JGA+12)='PTM1' MOTS(JGA+13)='RMAX' MOTS(JGA+14)='BR ' MOOBL=MLMOTS RETURN C ==== C 12 - Modele viscoplastique VISK2 C ==== ELSE IF (IPLAC.EQ.12) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='H ' MOTS(JGA+3)='ETA ' MOTS(JGA+4)='HVIS' MOTS(JGA+5)='N ' MOOBL=MLMOTS RETURN C ==== C 13 - Modele VISCOHINTE C ==== ELSE IF (IPLAC.EQ.13) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ,MLMOTS MOTS(JGA+1)='GAM1' MOTS(JGA+2)='AL ' MOTS(JGA+3)='KK ' MOTS(JGA+4)='NN ' MOTS(JGA+5)='YO ' MOTS(JGA+6)='YC ' MOTS(JGA+7)='MM ' MOOBL=MLMOTS RETURN C ==== C 14 - Modele viscoplastique MISTRAL C ==== ELSE IF (IPLAC.EQ.14) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 SEGADJ,MLMOTS MOTS(JGA+ 1)='SIP1' MOTS(JGA+ 2)='SIP2' MOTS(JGA+ 3)='IBAS' MOTS(JGA+ 4)='DILT' MOTS(JGA+ 5)='NBRE' MOTS(JGA+ 6)='COHI' MOTS(JGA+ 7)='ACOU' MOTS(JGA+ 8)='EDIR' MOTS(JGA+ 9)='RVCE' MOTS(JGA+10)='ECRX' MOTS(JGA+11)='DVDI' MOTS(JGA+12)='CROI' MOTS(JGA+13)='INCR' MOOBL=MLMOTS RETURN C ==== C 15 - Modele viscoplastique BETON_URGC C ==== ELSE IF (IPLAC.EQ.15) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+10 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' MOTS(JGA+10)='ETA' MOOBL=MLMOTS RETURN C ==== C 16 - Modele LEMENDO (fluage de type LEMaitre avec ENDOmmagement) C ==== ELSE IF (IPLAC.EQ.16) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='N ' MOTS(JGA+3)='M ' MOTS(JGA+4)='KK ' MOTS(JGA+5)='ALP1' MOTS(JGA+6)='BLP1' MOTS(JGA+7)='R ' MOTS(JGA+8)='A ' MOTS(JGA+9)='EXK ' MOOBL=MLMOTS RETURN C ==== C 17 - Modele FLUNOR2 (fluage de type Norton) C ==== ELSE IF (IPLAC.EQ.17) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+6 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='AF1 ' MOTS(JGA+3)='AF2 ' MOTS(JGA+4)='AF3 ' MOTS(JGA+5)='SFFS' MOTS(JGA+6)='EFFS' MOOBL=MLMOTS RETURN C ==== C 18 - Modele GATT_MONERIE C ==== ELSE IF (IPLAC.EQ.18) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+32 SEGADJ,MLMOTS MOTS(JGA+ 1)='R ' MOTS(JGA+ 2)='DG ' MOTS(JGA+ 3)='K1 ' MOTS(JGA+ 4)='M1 ' MOTS(JGA+ 5)='Q1 ' MOTS(JGA+ 6)='N1 ' MOTS(JGA+ 7)='K2 ' MOTS(JGA+ 8)='M2 ' MOTS(JGA+ 9)='Q2 ' MOTS(JGA+10)='N2 ' MOTS(JGA+11)='OMEG' MOTS(JGA+12)='Q ' MOTS(JGA+13)='H ' MOTS(JGA+14)='BETA' MOTS(JGA+15)='K ' MOTS(JGA+16)='A ' MOTS(JGA+17)='Q3 ' MOTS(JGA+18)='N3 ' MOTS(JGA+19)='ADEN' MOTS(JGA+20)='KGON' MOTS(JGA+21)='POR0' MOTS(JGA+22)='BUMI' MOTS(JGA+23)='EFIS' MOTS(JGA+24)='DG0 ' MOTS(JGA+25)='CR1 ' MOTS(JGA+26)='CR2 ' MOTS(JGA+27)='CR3 ' MOTS(JGA+28)='CR ' MOTS(JGA+29)='KP ' MOTS(JGA+30)='AP ' MOTS(JGA+31)='BP ' MOTS(JGA+32)='QP ' MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+6 SEGADJ,MLMOTS MOTS(JGA+ 1)='TYPE' MOTS(JGA+ 2)='COMP' MOTS(JGA+ 3)='DYN ' MOTS(JGA+ 4)='DYN1' MOTS(JGA+ 5)='DYN2' MOTS(JGA+ 6)='DYN3' MOFAC=MLMOTS RETURN C ==== C 19 - Modele UO2 = OTTOSEN + GATT_MONERIE C ==== ELSE IF (IPLAC.EQ.19) THEN C= 19.1 - Modes de calcul PLAN DEFO, PLAN GENE, Axisymetrie, Fourier IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-1.OR. . IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+39 SEGADJ,MLMOTS MOTS(JGA+ 1)='R ' MOTS(JGA+ 2)='DG ' MOTS(JGA+ 3)='K1 ' MOTS(JGA+ 4)='M1 ' MOTS(JGA+ 5)='Q1 ' MOTS(JGA+ 6)='N1 ' MOTS(JGA+ 7)='K2 ' MOTS(JGA+ 8)='M2 ' MOTS(JGA+ 9)='Q2 ' MOTS(JGA+10)='N2 ' MOTS(JGA+11)='OMEG' MOTS(JGA+12)='Q ' MOTS(JGA+13)='H ' MOTS(JGA+14)='BETA' MOTS(JGA+15)='K ' MOTS(JGA+16)='A ' MOTS(JGA+17)='Q3 ' MOTS(JGA+18)='N3 ' MOTS(JGA+19)='ADEN' MOTS(JGA+20)='KGON' MOTS(JGA+21)='POR0' MOTS(JGA+22)='BUMI' MOTS(JGA+23)='EFIS' MOTS(JGA+24)='DG0 ' MOTS(JGA+25)='CR1 ' MOTS(JGA+26)='CR2 ' MOTS(JGA+27)='CR3 ' MOTS(JGA+28)='CR ' MOTS(JGA+29)='KP ' MOTS(JGA+30)='AP ' MOTS(JGA+31)='BP ' MOTS(JGA+32)='QP ' IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-1) THEN MOTS(JGA+33)='LXX ' MOTS(JGA+34)='LYY ' MOTS(JGA+35)='LXY ' MOTS(JGA+36)='LZZ ' MOTS(JGA+37)='PXX ' MOTS(JGA+38)='PYY ' MOTS(JGA+39)='PXY ' C* ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN ELSE MOTS(JGA+33)='LRR ' MOTS(JGA+34)='LZZ ' MOTS(JGA+35)='LRZ ' MOTS(JGA+36)='LOO ' MOTS(JGA+37)='PRR ' MOTS(JGA+38)='PZZ ' MOTS(JGA+39)='PRZ ' ENDIF MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+28 SEGADJ,MLMOTS MOTS(JGA+ 1)='GFTR' MOTS(JGA+ 2)='LTR ' MOTS(JGA+ 3)='GS ' MOTS(JGA+ 4)='BTR ' MOTS(JGA+ 5)='WRUP' MOTS(JGA+ 6)='BILI' MOTS(JGA+ 7)='EPSR' MOTS(JGA+ 8)='GFT1' MOTS(JGA+ 9)='GFT2' MOTS(JGA+10)='GFT3' MOTS(JGA+11)='GS1 ' MOTS(JGA+12)='GS2 ' MOTS(JGA+13)='GS3 ' MOTS(JGA+14)='LTR1' MOTS(JGA+15)='LTR2' MOTS(JGA+16)='LTR3' MOTS(JGA+17)='WRU1' MOTS(JGA+18)='WRU2' MOTS(JGA+19)='BIL1' MOTS(JGA+20)='BIL2' MOTS(JGA+21)='EPSB' MOTS(JGA+22)='SIMP' MOTS(JGA+23)='TYPE' MOTS(JGA+24)='COMP' MOTS(JGA+25)='DYN ' MOTS(JGA+26)='DYN1' MOTS(JGA+27)='DYN2' MOTS(JGA+28)='DYN3' MOFAC=MLMOTS RETURN C= 19.2 - Modes de calcul PLAN CONT et TRIDimensionnel COQUEs minces ELSE IF ((IFOUR.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9)).OR. . IFOUR.EQ.-2) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+39 SEGADJ,MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)='R ' MOTS(JGA+ 2)='DG ' MOTS(JGA+ 3)='K1 ' MOTS(JGA+ 4)='M1 ' MOTS(JGA+ 5)='Q1 ' MOTS(JGA+ 6)='N1 ' MOTS(JGA+ 7)='K2 ' MOTS(JGA+ 8)='M2 ' MOTS(JGA+ 9)='Q2 ' MOTS(JGA+10)='N2 ' MOTS(JGA+11)='OMEG' MOTS(JGA+12)='Q ' MOTS(JGA+13)='H ' MOTS(JGA+14)='BETA' MOTS(JGA+15)='K ' MOTS(JGA+16)='A ' MOTS(JGA+17)='Q3 ' MOTS(JGA+18)='N3 ' MOTS(JGA+19)='ADEN' MOTS(JGA+20)='KGON' MOTS(JGA+21)='POR0' MOTS(JGA+22)='BUMI' MOTS(JGA+23)='EFIS' MOTS(JGA+24)='DG0 ' MOTS(JGA+25)='CR1 ' MOTS(JGA+26)='CR2 ' MOTS(JGA+27)='CR3 ' MOTS(JGA+28)='CR ' MOTS(JGA+29)='KP ' MOTS(JGA+30)='AP ' MOTS(JGA+31)='BP ' MOTS(JGA+32)='QP ' IF (IFOUR.EQ.-2) THEN MOTS(JGA+33)='LXX ' MOTS(JGA+34)='LYY ' MOTS(JGA+35)='LXY ' MOTS(JGA+36)='LZZ ' MOTS(JGA+37)='PXX ' MOTS(JGA+38)='PYY ' MOTS(JGA+39)='PXY ' C* ELSE IF (IFOUR.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9)) THEN ELSE MOTS(JGA+33)='LSS ' MOTS(JGA+34)='LTT ' MOTS(JGA+35)='LST ' MOTS(JGA+36)='LNN ' MOTS(JGA+37)='PSS ' MOTS(JGA+38)='PTT ' MOTS(JGA+39)='PST ' ENDIF MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+23 SEGADJ,MLMOTS MOFAC=MLMOTS MOTS(JGA+ 1)='GFTR' MOTS(JGA+ 2)='LTR ' MOTS(JGA+ 3)='GS ' MOTS(JGA+ 4)='BTR ' MOTS(JGA+ 5)='WRUP' MOTS(JGA+ 6)='BILI' MOTS(JGA+ 7)='GFT1' MOTS(JGA+ 8)='GFT2' MOTS(JGA+ 9)='GS1 ' MOTS(JGA+10)='GS2 ' MOTS(JGA+11)='LTR1' MOTS(JGA+12)='LTR2' MOTS(JGA+13)='WRU1' MOTS(JGA+14)='WRU2' MOTS(JGA+15)='BIL1' MOTS(JGA+16)='BIL2' MOTS(JGA+17)='SIMP' MOTS(JGA+18)='TYPE' MOTS(JGA+19)='COMP' MOTS(JGA+20)='DYN ' MOTS(JGA+21)='DYN1' MOTS(JGA+22)='DYN2' MOTS(JGA+23)='DYN3' RETURN C= 19.3 - Formulation TRIDIMENSIONNELLE MASSIVE ELSE IF (IFOUR.EQ.2.AND.MFR.EQ.1) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+44 SEGADJ,MLMOTS MOTS(JGA+ 1)='R ' MOTS(JGA+ 2)='DG ' MOTS(JGA+ 3)='K1 ' MOTS(JGA+ 4)='M1 ' MOTS(JGA+ 5)='Q1 ' MOTS(JGA+ 6)='N1 ' MOTS(JGA+ 7)='K2 ' MOTS(JGA+ 8)='M2 ' MOTS(JGA+ 9)='Q2 ' MOTS(JGA+10)='N2 ' MOTS(JGA+11)='OMEG' MOTS(JGA+12)='Q ' MOTS(JGA+13)='H ' MOTS(JGA+14)='BETA' MOTS(JGA+15)='K ' MOTS(JGA+16)='A ' MOTS(JGA+17)='Q3 ' MOTS(JGA+18)='N3 ' MOTS(JGA+19)='ADEN' MOTS(JGA+20)='KGON' MOTS(JGA+21)='POR0' MOTS(JGA+22)='BUMI' MOTS(JGA+23)='EFIS' MOTS(JGA+24)='DG0 ' MOTS(JGA+25)='CR1 ' MOTS(JGA+26)='CR2 ' MOTS(JGA+27)='CR3 ' MOTS(JGA+28)='CR ' MOTS(JGA+29)='KP ' MOTS(JGA+30)='AP ' MOTS(JGA+31)='BP ' MOTS(JGA+32)='QP ' MOTS(JGA+33)='LXX ' MOTS(JGA+34)='LYY ' MOTS(JGA+35)='LZZ ' MOTS(JGA+36)='LXY ' MOTS(JGA+37)='LXZ ' MOTS(JGA+38)='LYZ ' MOTS(JGA+39)='PXX ' MOTS(JGA+40)='PYY ' MOTS(JGA+41)='PZZ ' MOTS(JGA+42)='PXY ' MOTS(JGA+43)='PXZ ' MOTS(JGA+44)='PYZ ' MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+28 SEGADJ,MLMOTS MOTS(JGA+ 1)='GFTR' MOTS(JGA+ 2)='LTR ' MOTS(JGA+ 3)='GS ' MOTS(JGA+ 4)='BTR ' MOTS(JGA+ 5)='WRUP' MOTS(JGA+ 6)='BILI' MOTS(JGA+ 7)='GFT1' MOTS(JGA+ 8)='GFT2' MOTS(JGA+ 9)='GFT3' MOTS(JGA+10)='GS1 ' MOTS(JGA+11)='GS2 ' MOTS(JGA+12)='GS3 ' MOTS(JGA+13)='LTR1' MOTS(JGA+14)='LTR2' MOTS(JGA+15)='LTR3' MOTS(JGA+16)='WRU1' MOTS(JGA+17)='WRU2' MOTS(JGA+18)='WRU3' MOTS(JGA+19)='BIL1' MOTS(JGA+20)='BIL2' MOTS(JGA+21)='BIL3' MOTS(JGA+22)='SIMP' MOTS(JGA+23)='TYPE' MOTS(JGA+24)='COMP' MOTS(JGA+25)='DYN ' MOTS(JGA+26)='DYN1' MOTS(JGA+27)='DYN2' MOTS(JGA+28)='DYN3' MOFAC=MLMOTS RETURN C= 19.4 - Formulation UNIDIMENSIONNELLE (1D) MASSIVE ELSE IF ((IFOUR.GE.3.AND.IFOUR.LE.15).AND.MFR.EQ.1) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+37 SEGADJ,MLMOTS MOTS(JGA+ 1)='R ' MOTS(JGA+ 2)='DG ' MOTS(JGA+ 3)='K1 ' MOTS(JGA+ 4)='M1 ' MOTS(JGA+ 5)='Q1 ' MOTS(JGA+ 6)='N1 ' MOTS(JGA+ 7)='K2 ' MOTS(JGA+ 8)='M2 ' MOTS(JGA+ 9)='Q2 ' MOTS(JGA+10)='N2 ' MOTS(JGA+11)='OMEG' MOTS(JGA+12)='Q ' MOTS(JGA+13)='H ' MOTS(JGA+14)='BETA' MOTS(JGA+15)='K ' MOTS(JGA+16)='A ' MOTS(JGA+17)='Q3 ' MOTS(JGA+18)='N3 ' MOTS(JGA+19)='ADEN' MOTS(JGA+20)='KGON' MOTS(JGA+21)='POR0' MOTS(JGA+22)='BUMI' MOTS(JGA+23)='EFIS' MOTS(JGA+24)='DG0 ' MOTS(JGA+25)='CR1 ' MOTS(JGA+26)='CR2 ' MOTS(JGA+27)='CR3 ' MOTS(JGA+28)='CR ' MOTS(JGA+29)='KP ' MOTS(JGA+30)='AP ' MOTS(JGA+31)='BP ' MOTS(JGA+32)='QP ' C*OF A finir en fonction des definitions de ces grandeurs C*OF et en fonction de IFOUR IF (IFOUR.LE.11) THEN MOTS(JGA+33)='LXX ' MOTS(JGA+34)='LYY ' MOTS(JGA+35)='LZZ ' MOTS(JGA+36)='PXX ' MOTS(JGA+37)='PYY ' ELSE MOTS(JGA+33)='LRR ' MOTS(JGA+34)='LZZ ' MOTS(JGA+35)='LOO ' MOTS(JGA+36)='PRR ' MOTS(JGA+37)='PZZ ' ENDIF MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+28 C* A finir en fonction de la def. des parametres ! SEGADJ,MLMOTS MOTS(JGA+ 1)='GFTR' MOTS(JGA+ 2)='LTR ' MOTS(JGA+ 3)='GS ' MOTS(JGA+ 4)='BTR ' MOTS(JGA+ 5)='WRUP' MOTS(JGA+ 6)='BILI' MOTS(JGA+ 7)='GFT1' MOTS(JGA+ 8)='GFT2' MOTS(JGA+ 9)='GFT3' MOTS(JGA+10)='GS1 ' MOTS(JGA+11)='GS2 ' MOTS(JGA+12)='GS3 ' MOTS(JGA+13)='LTR1' MOTS(JGA+14)='LTR2' MOTS(JGA+15)='LTR3' MOTS(JGA+16)='WRU1' MOTS(JGA+17)='WRU2' MOTS(JGA+18)='WRU3' MOTS(JGA+19)='BIL1' MOTS(JGA+20)='BIL2' MOTS(JGA+21)='BIL3' MOTS(JGA+22)='SIMP' MOTS(JGA+23)='TYPE' MOTS(JGA+24)='COMP' MOTS(JGA+25)='DYN ' MOTS(JGA+26)='DYN1' MOTS(JGA+27)='DYN2' MOTS(JGA+28)='DYN3' MOFAC=MLMOTS RETURN ENDIF C ==== C 20 - Modele viscoplastique BETON_URGC_ENDO C ==== ELSE IF (IPLAC.EQ.20) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+12 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)='DC ' MOTS(JGA+ 9)='DT ' MOTS(JGA+10)='IASO' MOTS(JGA+11)='ICOU' MOTS(JGA+12)='ETA ' MOOBL=MLMOTS RETURN C ==== C 21 - Modele viscoplastique double endommagement C ==== ELSE IF (IPLAC.EQ.21) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+11 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='B' MOTS(JGA+3)='RI' MOTS(JGA+4)='K' MOTS(JGA+5)='N' MOTS(JGA+6)='SD' MOTS(JGA+7)='RD' MOTS(JGA+8)='PD' MOTS(JGA+9)='SC' MOTS(JGA+10)='RC' MOTS(JGA+11)='PC' MOOBL=MLMOTS RETURN C ==== C 22 - Modele viscoplastique CHAB_SINH_R C ==== ELSE IF (IPLAC.EQ.22) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+6 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='B' MOTS(JGA+3)='RI' MOTS(JGA+4)='K' MOTS(JGA+5)='N' MOTS(JGA+6)='epso' MOOBL=MLMOTS RETURN C ==== C 23 - Modele viscoplastique CHAB_SINH_X C ==== ELSE IF (IPLAC.EQ.23) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='B' MOTS(JGA+3)='RI' MOTS(JGA+4)='K' MOTS(JGA+5)='N' MOTS(JGA+6)='epso' MOTS(JGA+7)='C1' MOTS(JGA+8)='gam1' MOTS(JGA+9)='C2' MOOBL=MLMOTS RETURN C ==== C 24 - Modele viscoplastique CHAB_NOR_R C ==== ELSE IF (IPLAC.EQ.24) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='B' MOTS(JGA+3)='RI' MOTS(JGA+4)='K' MOTS(JGA+5)='N' MOOBL=MLMOTS RETURN C ==== C 25 - Modele viscoplastique CHAB_NOR_X C ==== ELSE IF (IPLAC.EQ.25) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+8 SEGADJ,MLMOTS MOTS(JGA+1)='SIGY' MOTS(JGA+2)='B' MOTS(JGA+3)='RI' MOTS(JGA+4)='K' MOTS(JGA+5)='N' MOTS(JGA+6)='C1' MOTS(JGA+7)='gam1' MOTS(JGA+8)='C2' MOOBL=MLMOTS RETURN C === C 26 - Modele d'ecrouissage-VISCOsite multiplicative avec enDOMMAGEment anisotrope C === ELSE IF (IPLAC.EQ.26) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+8 SEGADJ,MLMOTS MOTS(JGA+1)='N ' MOTS(JGA+2)='M ' MOTS(JGA+3)='KK ' MOTS(JGA+4)='ALF2' MOTS(JGA+5)='BET2' MOTS(JGA+6)='R ' MOTS(JGA+7)='A ' MOTS(JGA+8)='Q ' MOOBL=MLMOTS MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+1 SEGADJ,MLMOTS MOTS(JGA+1)='SMAX' MOFAC=MLMOTS RETURN C === C 27 - Modele FLUTRA C === ELSE IF (IPLAC.EQ.27) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+24 SEGADJ,MLMOTS MOTS(JGA+1)= 'KTR0' MOTS(JGA+2)= 'ATRA' MOTS(JGA+3)= 'BTRA' MOTS(JGA+4)= 'ACOM' MOTS(JGA+5)= 'BCOM' MOTS(JGA+6)= 'EKV1' MOTS(JGA+7)= 'TKV1' MOTS(JGA+8)= 'EKV2' MOTS(JGA+9)= 'TKV2' MOTS(JGA+10)= 'EKV3' MOTS(JGA+11)= 'TKV3' MOTS(JGA+12)= 'EKV4' MOTS(JGA+13)= 'TKV4' MOTS(JGA+14)= 'EKV5' MOTS(JGA+15)= 'TKV5' MOTS(JGA+16)= 'EKV6' MOTS(JGA+17)= 'TKV6' MOTS(JGA+18)= 'BETA' MOTS(JGA+19)= 'ETA1' MOTS(JGA+20)= 'NUPO' MOTS(JGA+21)= 'XFDE' MOTS(JGA+22)= 'SLMA' MOTS(JGA+23)= 'VGEA' MOTS(JGA+24)= 'VGEB' MOOBL=MLMOTS RETURN C === C 28 - Modele SYMONDS & COWPER SYCO1 C === ELSE IF (IPLAC.EQ.28) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOTS(JGA+1)='PSYC' MOTS(JGA+2)='DSYC' MOOBL=MLMOTS RETURN C === C 29 - Modele SYMONDS & COWPER SYCO2 C === ELSE IF (IPLAC.EQ.29) THEN MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOTS(JGA+1)='PSYC' MOTS(JGA+2)='ASYC' MOTS(JGA+3)='BSYC' MOTS(JGA+4)='CSYC' MOOBL=MLMOTS RETURN ENDIF C ========== C ERREUR : modele non reference (IRET=0) C ========== IRET=0 CALL ERREUR(261) RETURN END