1 : C POUDRB SOURCE PV 09/12/22 21:15:17 6589 2 : C 3 : C PLGDebut : implantation des equations relatives au modele 4 : C viscoplastique poudre. 5 : C 6 : C Modification ChD division par des rĹels et non pas des entiers 7 : C COEC doit etre plus grand que 1.0 et non pas 0.0 8 : C 9 : SUBROUTINE POUDRB(SIG,EPSV,VAR,EPSVPT,VARPT,XMAT,NSTRS, 10 : + NVARI,NCOMAT,KERRE) 11 : C 12 : IMPLICIT INTEGER(I-N) 13 : IMPLICIT REAL*8(A-H,O-Z) 14 : -INC CCOPTIO 15 : DIMENSION SIG(*),EPSV(*),VAR(*),EPSVPT(*),VARPT(*) 16 : DIMENSION XMAT(*) 17 : DIMENSION XX(6) 18 : C---------------------------------------------------------------- 19 : C ====> Connaissant l'etat, c'est a dire SIG, EPSV et VAR on 20 : C calcule les derivees EPSVPT et VARPT(1) en utilisant les 21 : C constantes XMAT 22 : C---------------------------------------------------------------- 23 : C 24 : C-------calcul des invariants S1=Tr(sigma) 25 : C S2C=S2*S2=Tr(Dev(sigma):Dev(sigma)) 26 : C avec Dev(sigma)=sigma-1/3Tr(sigma) 27 : S1 = SIG(1)+SIG(2)+SIG(3) 28 : C Modif Ch.D du 20/10/94 : Cf DEA Canel page A18 assure tr(XX) = 0 29 : DO 30 I=1,2 30 : XX(I) = SIG(I)-S1/3.0D0 31 : 30 CONTINUE 32 : XX(3) = -XX(1)-XX(2) 33 : C Fin de modification Ch.D du 20/10/94 34 : DO 40 I=4,NSTRS,1 35 : XX(I) = SIG(I) 36 : 40 CONTINUE 37 : S2C = PROCON(XX,XX,NSTRS) 38 : C 39 : RHOR = VAR(NVARI) 40 : A = XMAT(6) 41 : RN = XMAT(5) 42 : QSRT = XMAT(7) 43 : F0 = XMAT(8) 44 : F1 = XMAT(9) 45 : F2 = XMAT(10) 46 : F3 = XMAT(11) 47 : F4 = XMAT(12) 48 : F5 = XMAT(13) 49 : C0 = XMAT(14) 50 : C1 = XMAT(15) 51 : C2 = XMAT(16) 52 : C3 = XMAT(17) 53 : C4 = XMAT(18) 54 : C5 = XMAT(19) 55 : G0 = XMAT(20) 56 : G1 = XMAT(21) 57 : G2 = XMAT(22) 58 : G3 = XMAT(23) 59 : G4 = XMAT(24) 60 : G5 = XMAT(25) 61 : G6 = XMAT(26) 62 : C 63 : RMIN = F5+1.0D-03 64 : IF (RHOR.LT.RMIN) THEN 65 : RHOR1 = RMIN 66 : ELSE 67 : RHOR1 = RHOR 68 : ENDIF 69 : C 70 : C 71 : C-------calcul de EPSVPT en fonction du parametre F4 : 72 : C------- F4 = 0 : loi ABOUAF avec taille de grain 73 : C------- F4 = 1 : loi de Stutz 74 : C 75 : *pv pour ętre sur d'avoir la bonne valeur 76 : NF4=INT(F4+0.1)+1 77 : GOTO(10,20),NF4 78 : C 79 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 80 : C 81 : C Loi Abouaf avec calcul de la taille de grain 82 : C 83 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 84 : C 85 : 10 CONTINUE 86 : C 87 : C-------calcul de F et C : fonctions de rho 88 : C 89 : IF ((1.0D0-RHOR1).LT.(1.D-06)) THEN 90 : COEF = 0.0 91 : COEC = 1.0 92 : ELSE 93 : COEF = F0*(((1-RHOR1)/(RHOR1-F5))**(F1*RHOR1+F2)) 94 : COEC = 1.0D0+C0*(((1-RHOR1)/(RHOR1-F5))**(C1*RHOR1+C2)) 95 : ENDIF 96 : IF (COEF.GT.(C5)) COEF = C5 97 : IF (COEC.GT.(C5)) COEC = C5 98 : C 99 : C-------calcul de sigma equivalent : 100 : C 101 : SIGEQ = SQRT(COEF*S1*S1+1.5D0*COEC*S2C) 102 : DO 500 I=1,NSTRS,1 103 : B = 0.D0 104 : IF (I.LE.3) B=1.D0 105 : EPSVPT(I) = 1.D0*(10**A)*(SIGEQ**(RN-1.D0))*RHOR 106 : & *(COEF*S1*B+1.5D0*COEC*XX(I)) 107 : 500 CONTINUE 108 : GOTO 1000 109 : C 110 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 111 : C 112 : C Loi de Stutz 113 : C 114 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 115 : C 116 : 20 CONTINUE 117 : COEH = F0*(((1.-RHOR)/(RHOR-F5))**F1) 118 : COEE = C0*(RHOR-RHOR**2) 119 : COES = G0-(G1*(1.-RHOR))+(G2*((1.-RHOR)**2)) 120 : COF1 = (2./3.)**((RN+1)/2)*COEH/3. 121 : COF2 = (1.-COEE)*(COES**(-RN))+2.*((3./2.)**((RN+1)/2))*COF1 122 : COF3 = COEE*(COES**(-RN))*((2./3.)**((RN+5)/4)) 123 : & -2.*COF1*((3./2.)**((RN+1)/4)) 124 : C 125 : C-------calcul de EPSVPT : 126 : C 127 : C 128 : DO 600 I=1,NSTRS,1 129 : B = 0.0 130 : IF (I.LE.3) B=1.0 131 : SIGEQ = (S1**2*S2C)**((RN-3)/4) 132 : EPSVPT(I) = (3./2.)**((RN+1)/2)*A*EXP(-QSRT)*RHOR 133 : & *(3.*S1*(COF1*S1**(RN-1)+COF3/2.*S2C*SIGEQ)*B 134 : & +(COF2*S2C**((RN-1)/2)+COF3/2.*S1**2*SIGEQ)*XX(I)) 135 : 600 CONTINUE 136 : C 137 : C-------calcul de VARPT(1) derivee de EPSE C 138 : 1000 CONTINUE 139 : PPT = SQRT(2.0/3.0*PROCON(EPSVPT,EPSVPT,NSTRS)) 140 : VARPT(1) = PPT 141 : C 142 : C--------on rentre 143 : C 144 : RETURN 145 : END 146 : C 147 : 148 : 149 : 150 :
© Cast3M 2003 - All rights reserved.
Disclaimer