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