1 : C CRIDAM SOURCE CHAT 05/01/12 22:30:15 5004 2 : SUBROUTINE CRIDAM(SIGG,NSTRS,VARR,XMAT,NCOMAT,XCAR,ICARA, 3 : 1INPLAS,NTRAC,TRACC,NC,NPNTER,NNVARI,MFR,XF,KERRE,SIGGD, 4 : 2YY,T,TRAC,NCOURB) 5 : C 6 : C =================================================================== 7 : C CE SOUS-PROGRAMME EST APPELE PAR "ENDOM". 8 : C IL CALCULE LA VALEUR DU CRITERE POUR CHAQUE MODELE DE MATERIAU 9 : C ENDOMMAGEABLE. 10 : C 11 : C ENTREES: 12 : C ------- 13 : C NSTRS = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM. 14 : C SIGG(NSTRS) = CONTR. 15 : C NNVARI = NBR. DE VARIABLES INTERNES PILOTANT LES EQUATIONS 16 : C DU MODELE 17 : C VARR(2)... = VARIABLES INTERNES PILOTANT LES EQUATIONS DU MODELE 18 : C ...VARR(1+NNVARI) 19 : C NCOMAT = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU 20 : C XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU 21 : C MFR = INDICE DE LA FORMULATION MECANIQUE; SEULEMENT 22 : C MASSIF OU COQUE POUR LES MATERIAUX ENDOMMAGEABLES 23 : C ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS 24 : C XCAR(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS 25 : C INPLAS = INDICE DU MODELE DECRIVANT LE MATERIAU ENDOMMAGEABLE 26 : C NTRAC = DIMENSION DE TRACC 27 : C TRACC(NTRAC) = COURBES DE TRACTION 28 : C NC = NBR. DE COURBES DE TRACTION 29 : C NPNTER(NC)= TABLEAU CONTENANT LE NOMBRE DE POINTS POUR CHAQUE 30 : C COURBE DE TRACTION (A T FIXEE) 31 : C T = TEMPERATURE A LAQUELLE ON CALCULE LE CRITERE 32 : C 33 : C SORTIES: 34 : C ------- 35 : C XF= VALEUR DU CRITERE 36 : C KERRE= INDICE QUI REGIT LES ERREURS 37 : C = 77 SI LA DEFORM. PLAST. CUMULEE ENDOMMAGEE (VARIABLE VIR) 38 : C EST EN DEHORS DE LA COURBE DE TRACTION, DS. LE CAS DE 39 : C L'ECROUISSAGE ET DE L'ENDOMMAGEMENT ISOTROPES 40 : C = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE POUR 41 : C LE MODELE CONSIDERE 42 : C 43 : C ================================================================== 44 : C 45 : IMPLICIT INTEGER(I-N) 46 : IMPLICIT REAL*8 (A-H,O-Z) 47 : DIMENSION SIGG(*),VARR(*),XMAT(*),XCAR(*),TRACC(*),SIGGD(*), 48 : 1 TRAC(*),NPNTER(*) 49 : C 50 : KERRE=0 51 : IPLAS=INPLAS-25 52 : C 53 : C ------------------------------------------------------------- 54 : C LE CALCUL DU CRITERE DEPEND DU MODELE DONC DE INPLAS OU IPLAS 55 : C ------------------------------------------------------------- 56 : GOTO (10,1000,1000) IPLAS 57 : C 58 : C 59 : C ================================================== 60 : C ECROUISSAGE ET ENDOMMAGEMENT ISOTROPES - VON MISES 61 : C ================================================== 62 : 10 VIR=VARR(2) 63 : XD=VARR(3) 64 : C 65 : IF (MFR.NE.1) GOTO 50 66 : C 67 : C ------- 68 : C MASSIFS 69 : C ------- 70 : C 70 : CALL TRACTT(TRACC,NPNTER,TRAC,NTRAC,NC,NCOURB,T) 71 : N2COUR=2*NCOURB 71 : CALL DERTRA(N2COUR,TRAC,VIR,YY,YPRIM,XINF,XSUP) 72 : C 73 : TRACE=SIGG(1)+SIGG(2)+SIGG(3) 74 : C 75 : C .............................. 76 : C CALCUL DU DEVIATEUR DES CONTR. 77 : C .............................. 77 : CALL ZDANUL(SIGGD,NSTRS) 78 : DO 30 I=1,NSTRS 79 : A=0.D0 80 : IF (I.LE.3) A=1.D0/3.D0 81 : SIGGD(I)=SIGG(I)-A*TRACE 82 : 30 CONTINUE 83 : C 84 : C ............................................ 85 : C CALCUL DE LA CONTR. EQUIVALENTE DE VON MISES 86 : C ............................................ 87 : SIGEQI=PROCON(SIGGD,SIGGD,NSTRS) 88 : SIGEQI=DSQRT((3.D0/2.D0)*SIGEQI) 89 : C 90 : C ................. 91 : C CALCUL DU CRITERE 92 : C ................. 93 : XF=SIGEQI/(1.D0-XD)-YY 94 : GOTO 1000 95 : C 96 : 50 IF (MFR.NE.5) GOTO 998 97 : C 98 : C --------------------------------------- 99 : C COQUES EPAISSES - NON ENCORE DISPONIBLE 100 : C --------------------------------------- 101 : C 102 : 998 KERRE=99 103 : 1000 CONTINUE 104 : RETURN 105 : END 106 : 107 :
© Cast3M 2003 - All rights reserved.
Disclaimer