C IDENDO SOURCE BR232186 12/08/01 22:28:49 7456 SUBROUTINE IDENDO (MFR,MOOBL,MOFAC,IPLAC,IRET) *--------------------------------------------------------------------* * Noms de composantes de materiaux endommageables * *--------------------------------------------------------------------* * * * ENTREES: * * ------- * * MFR num{ro de formulation * * MOOBL pointeur sur les mots obligatoires * * MOFAC POINTEUR SUR LES MOTS FACULTATIFS * * IPLAC numero du materiau dans MODEND * * * * SORTIE : * * ------- * * * * IRET : 0 SI PROBLEME * * * *--------------------------------------------------------------------* * IMPLICIT INTEGER(I-N) -INC SMLMOTS -INC CCOPTIO * JGN=8 IRET=1 * *--------------------------------------------------------------------* * * * T Y P E E N D O M M A G E A B L E * * * *--------------------------------------------------------------------* IF (IPLAC.EQ.1) THEN * * MODELE ENDOMMAGEMENT SCALAIRE DE MAZARS * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+6 SEGADJ MLMOTS MOOBL=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' ELSE IF (IPLAC.EQ.2) THEN * * MODELE ENDOMMAGEMENT UNILATERAL LA BORDERIE * 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' ELSE IF (IPLAC.EQ.3) THEN * * MODELE ROTATING_CRACK * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+2 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'EPCR' MOTS(JGA+2)= 'MUP ' ELSE IF (IPLAC.EQ.4) THEN * * MODELE SIC_SIC * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+15 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'G1DC' MOTS(JGA+2)= 'G1Y0' MOTS(JGA+3)= 'G1YC' MOTS(JGA+4)= 'G1P' MOTS(JGA+5)= 'G2DC' MOTS(JGA+6)= 'G2Y0' MOTS(JGA+7)= 'G2YC' MOTS(JGA+8)= 'G2P' MOTS(JGA+9)= 'G3DC' MOTS(JGA+10)= 'G3Y0' MOTS(JGA+11)= 'G3YC' MOTS(JGA+12)= 'G3P' MOTS(JGA+13)= 'H23' MOTS(JGA+14)= 'H13' MOTS(JGA+15)= 'H12' MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+3 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1)= 'A12' MOTS(JGA+2)= 'A23' MOTS(JGA+3)= 'A13' * ELSE IF (IPLAC.EQ.5 .OR. IPLAC .EQ. 7) THEN * * MODELE ENDOMMAGEMENT MICROPLANS * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'KTR0' MOTS(JGA+2)= 'BTRA' MOTS(JGA+3)= 'FPC ' MOTS(JGA+4)= 'MP ' MOTS(JGA+5)= 'ALF1' * ELSE IF (IPLAC.EQ.6) THEN * * MODELE ENDOMMAGEMENT UNILATERAL LA BORDERIE * + VISCOSITE J.F. DUBE * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 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' MOTS(JGA+10)= 'M1 ' MOTS(JGA+11)= 'M2 ' MOTS(JGA+12)= 'ALP1' MOTS(JGA+13)= 'ALP2' * BCN ELSE IF (IPLAC.EQ.8) THEN * * Modified Von mises model * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'Y0 ' MOTS(JGA+2)= 'B1 ' MOTS(JGA+3)= 'B2 ' MOTS(JGA+4)= 'RATI' MOTS(JGA+5)= 'LOI ' * BCN * RDj ELSE IF (IPLAC.EQ.9) THEN * * MODELE ENDOMMAGEMENT SCALAIRE AVEC DEFORMATION ANELASTIQUE * CRUSH_CRACK MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+9 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'KCO0' MOTS(JGA+2)= 'KTR0' MOTS(JGA+3)= 'ACRU' MOTS(JGA+4)= 'BCRU' MOTS(JGA+5)= 'ACRA' MOTS(JGA+6)= 'BCRA' MOTS(JGA+7)= 'ECOM' MOTS(JGA+8)= 'ECOB' MOTS(JGA+9)= 'ETRA' * RDj * CCd ELSE IF (IPLAC.EQ.10) THEN * *+++++++++++++++++++++++++++++++++++ * MODELE DAMAGE_TENSION&COMPRESSION = DAMAGE_TC *+++++++++++++++++++++++++++++++++++ * * COMPULSORY CARACTERISTICS * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+15 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='HLEN' MOTS(JGA+2)='GVAL' MOTS(JGA+3)='FTUL' MOTS(JGA+4)='REDC' MOTS(JGA+5)='FC01' MOTS(JGA+6)='RT45' MOTS(JGA+7)='FCU1' MOTS(JGA+8)='STRU' MOTS(JGA+9)='EXTP' MOTS(JGA+10)='STRP' MOTS(JGA+11)='EXT1' MOTS(JGA+12)='STR1' MOTS(JGA+13)='EXT2' MOTS(JGA+14)='STR2' MOTS(JGA+15)='NCRI' * * FACULTATIVE CARACTERISTICS * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+10 SEGADJ MLMOTS MOFAC=MLMOTS C WARNING !!!!!!!!!! RHO et ALPH doivent etre a la fin! voir defmat!!! DO ie1=1,JGA MOTS(JGM+1-ie1)=MOTS(JGA+1-ie1) ENDDO C MOTS(1)='TANB' MOTS(2)='BPLC' MOTS(3)='TR0T' MOTS(4)='TR0C' MOTS(5)='TRUC' MOTS(6)='AVAT' MOTS(7)='BVAT' MOTS(8)='AVAC' MOTS(9)='BVAC' MOTS(10)='DVOL' * * CCd * ELSE IF (IPLAC.EQ.11) THEN * * MODELE SICSCAL * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+44 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1)= 'H1N' MOTS(JGA+2)= 'H1HP' MOTS(JGA+3)= 'H1P' MOTS(JGA+4)= 'H2N' MOTS(JGA+5)= 'H2HP' MOTS(JGA+6)= 'H2P' MOTS(JGA+7)= 'H3N' MOTS(JGA+8)= 'H3P' MOTS(JGA+9)= 'Y0N1' MOTS(JGA+10)='Y0N2' MOTS(JGA+11)='Y0N3' MOTS(JGA+12)= 'YCN1' MOTS(JGA+13)= 'YCN2' MOTS(JGA+14)= 'YCN3' MOTS(JGA+15)='Y0T1' MOTS(JGA+16)='Y0T2' MOTS(JGA+17)='Y0T3' MOTS(JGA+18)='YCT1' MOTS(JGA+19)='YCT2' MOTS(JGA+20)='YCT3' MOTS(JGA+21)='DCT1' MOTS(JGA+22)='DCT2' MOTS(JGA+23)='DCT3' MOTS(JGA+24)='DCN1' MOTS(JGA+25)='DCN2' MOTS(JGA+26)='DCN3' MOTS(JGA+27)='PN1' MOTS(JGA+28)='PN2' MOTS(JGA+29)='PN3' MOTS(JGA+30)='PT1' MOTS(JGA+31)='PT2' MOTS(JGA+32)='PT3' MOTS(JGA+33)='B' MOTS(JGA+34)='DTAL' MOTS(JGA+35)='TER0' MOTS(JGA+36)='SIF1' MOTS(JGA+37)='SIF2' MOTS(JGA+38)='SIF3' MOTS(JGA+39)='AIF1' MOTS(JGA+40)='AIF2' MOTS(JGA+41)='AIF3' MOTS(JGA+42)='ETA1' MOTS(JGA+43)='ETA2' MOTS(JGA+44)='ETA3' * ELSE IF (IPLAC.EQ.12) THEN * * MODELE SICTENS * MLMOTS=MOFAC JGA=MOTS(/2) JGM=JGA+51 SEGADJ MLMOTS MOFAC=MLMOTS MOTS(JGA+1)= 'H1N' MOTS(JGA+2)= 'H1HP' MOTS(JGA+3)= 'H1P' MOTS(JGA+4)= 'H2N' MOTS(JGA+5)= 'H2HP' MOTS(JGA+6)= 'H2P' MOTS(JGA+7)= 'H3N' MOTS(JGA+8)= 'H3P' MOTS(JGA+9)= 'H4N' MOTS(JGA+10)= 'H4HP' MOTS(JGA+11)= 'H4P' MOTS(JGA+12)= 'Y01' MOTS(JGA+13)='Y02' MOTS(JGA+14)='Y03' MOTS(JGA+15)='Y04' MOTS(JGA+16)='Y05' MOTS(JGA+17)= 'YC1' MOTS(JGA+18)= 'YC2' MOTS(JGA+19)= 'YC3' MOTS(JGA+20)= 'YC4' MOTS(JGA+21)= 'YC5' MOTS(JGA+22)='DC1' MOTS(JGA+23)='DC2' MOTS(JGA+24)='DC3' MOTS(JGA+25)='DC4' MOTS(JGA+26)='DC5' MOTS(JGA+27)='PY1' MOTS(JGA+28)='PY2' MOTS(JGA+29)='PY3' MOTS(JGA+30)='PY4' MOTS(JGA+31)='PY5' MOTS(JGA+32)='B1' MOTS(JGA+33)='B2' MOTS(JGA+34)='B3' MOTS(JGA+35)='DTAL' MOTS(JGA+36)='TER0' MOTS(JGA+37)='SIF1' MOTS(JGA+38)='SIF2' MOTS(JGA+39)='SIF3' MOTS(JGA+40)='SIF4' MOTS(JGA+41)='SIF5' MOTS(JGA+42)='AIF1' MOTS(JGA+43)='AIF2' MOTS(JGA+44)='AIF3' MOTS(JGA+45)='AIF4' MOTS(JGA+46)='AIF5' MOTS(JGA+47)='ETA1' MOTS(JGA+48)='ETA2' MOTS(JGA+49)='ETA3' MOTS(JGA+50)='ETA4' MOTS(JGA+51)='ETA5' * ELSE IF (IPLAC.EQ.13) THEN * * MODELE DESMORAT * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='K0' MOTS(JGA+2)='A' MOTS(JGA+3)='a' MOTS(JGA+4)='etaC' MOTS(JGA+5)='etaT' MOTS(JGA+6)='k' MOTS(JGA+7)='Dc' * ELSE IF (IPLAC.EQ.14) THEN * * MODELE ENDOMMAGEMENT FATIGUE 'FATSIN' * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+8 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)= 'KTR0' MOTS(JGA+2)= 'BETA' MOTS(JGA+3)= 'LOI ' * pour la loi "L2R" MOTS(JGA+4)= 'C ' MOTS(JGA+5)= 'ALFA' * pour la loi "L3R" MOTS(JGA+6)= 'ALF1' MOTS(JGA+7)= 'ALF2' MOTS(JGA+8)= 'ALF3' ELSE IF (IPLAC.EQ.15) THEN * * MODELE ENDOMMAGEMENT 'RICRAG' * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+5 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='FT ' MOTS(JGA+2)='ALIN' MOTS(JGA+3)='GAM1' MOTS(JGA+4)='A1 ' MOTS(JGA+5)='ALDI' C +BR ELSE IF (IPLAC.EQ.16) THEN * * MODELE ENDOMMAGEMENT 'GLRC_DM' * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+7 SEGADJMLMOTS MOOBL=MLMOTS MOTS(JGA+1)='YOUF' MOTS(JGA+2)='NUF ' MOTS(JGA+3)='GAMT' MOTS(JGA+4)='GAMC' MOTS(JGA+5)='GAMF' MOTS(JGA+6)='SEUI' MOTS(JGA+7)='ALF ' ELSE IF (IPLAC.EQ.17) THEN * * MODELE ENDOMMAGEMENT 'RICBET' * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+13 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+ 1)= 'FT ' MOTS(JGA+ 2)= 'ALDI' MOTS(JGA+ 3)= 'GAM1' MOTS(JGA+ 4)= 'A1 ' MOTS(JGA+ 5)= 'SREF' MOTS(JGA+ 6)= 'AF ' MOTS(JGA+ 7)= 'BF ' MOTS(JGA+ 8)= 'AG ' MOTS(JGA+ 9)= 'BG ' MOTS(JGA+10)= 'AC ' MOTS(JGA+11)= 'BC ' MOTS(JGA+12)= 'SIGU' MOTS(JGA+13)= 'FC ' ELSE IF (IPLAC.EQ.18) THEN * * MODELE ENDOMMAGEMENT 'RICCOQ' * MLMOTS=MOOBL JGA=MOTS(/2) JGM=JGA+4 SEGADJ MLMOTS MOOBL=MLMOTS MOTS(JGA+1)='FT ' MOTS(JGA+2)='EPUT' MOTS(JGA+3)='FC ' MOTS(JGA+4)='EPUC' C -BR ELSE * * SI PROBLEME , ON SORT AVEC IRET A 0 * IRET=0 CALL ERREUR(261) ENDIF RETURN END