1 : C POUDRA SOURCE PV 09/12/22 21:15:02 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 : C Modification ChD nouvelles equations de c et f selon 10 : C les equations de DA pour sa these. 11 : C 12 : C Modification LS nouvelle loi viscoplastique a trois parametres 13 : C Lois viscoplastiques d Abouaf et 3S sont disponibles ... 14 : C 15 : SUBROUTINE POUDRA(SIG,EPSV,VAR,EPSVPT,VARPT,XMAT,NSTRS, 16 : & NVARI,NCOMAT,KERRE) 17 : C 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : -INC CCOPTIO 21 : DIMENSION SIG(*),EPSV(*),VAR(*),EPSVPT(*),VARPT(*) 22 : DIMENSION XMAT(*) 23 : DIMENSION XX(6) 24 : C 25 : C--------------------------------------------------------------------- 26 : C ====> Connaissant l'etat, c'est a dire SIG, EPSV et VAR on 27 : C calcule les derivees EPSVPT et VARPT(1) en utilisant les 28 : C constantes XMAT 29 : C----- calcul des invariants S1=Tr(sigma) 30 : C S2C=S2*S2=Tr(Dev(sigma):Dev(sigma)) 31 : C avec Dev(sigma)=sigma-1/3Tr(sigma) 32 : C 33 : S1 = SIG(1)+SIG(2)+SIG(3) 34 : C 35 : C Modif Ch.D du 20/10/94 : Cf DEA Canel page A18 36 : C assure tr(XX) = 0 37 : DO I=1,2 38 : XX(I) = SIG(I)-S1/3.0D0 39 : END DO 40 : XX(3) = -XX(1)-XX(2) 41 : C Fin de modification Ch.D du 20/10/94 42 : C 43 : DO I=4,NSTRS,1 44 : XX(I) = SIG(I) 45 : END DO 46 : S2C = PROCON(XX,XX,NSTRS) 47 : C 48 : C----- calcul de F et C : fonctions de rho 49 : RHOR = VAR(NVARI) 50 : A = XMAT(6) 51 : RN = XMAT(5) 52 : QSRT = XMAT(7) 53 : F0 = XMAT(8) 54 : F1 = XMAT(9) 55 : F2 = XMAT(10) 56 : F3 = XMAT(11) 57 : F4 = XMAT(12) 58 : F5 = XMAT(13) 59 : C0 = XMAT(14) 60 : C1 = XMAT(15) 61 : C2 = XMAT(16) 62 : C3 = XMAT(17) 63 : C4 = XMAT(18) 64 : C5 = XMAT(19) 65 : RMIN = F5+1.0D-03 66 : IF (RHOR.LT.RMIN) THEN 67 : RHOR1 = RMIN 68 : ELSE 69 : RHOR1 = RHOR 70 : ENDIF 71 : C 72 : C----- Modification L3S du PROGRAMME : 73 : C Integration de la loi a 3 parametres 74 : C Le 15/10/1997 par Laurent SANCHEZ ... 75 : C COEF : Coefficient F(Rhor) 76 : C COEC : Coefficient C(Rhor) 77 : C COEE : Coefficient E(Rhor) 78 : C COSS : Coefficient S(Rhor) 79 : C 80 : C COF1 : Coefficient f1(Rhor) 81 : C COF2 : Coefficient f2(Rhor) 82 : C COF3 : Coefficient f3(Rhor) 83 : C 84 : IF ((1.0D0-RHOR1).LT.1.0D-06) THEN 85 : COEF = 0.0D0 86 : COEC = 1.0D0 87 : ELSE 88 : COEF = F0*(((1.D0-RHOR1)/(RHOR1-F5))**(F1*RHOR1+F2)) 89 : COEC = 1.0D0+C0*(((1.D0-RHOR1)/(RHOR1-F5))**(C1*RHOR1+C2)) 90 : ENDIF 91 : IF (COEF.GT.(C5)) COEF = C5 92 : IF (COEC.GT.(C5)) COEC = C5 93 : C 94 : C----- calcul de sigma equivalent : 95 : SIGEQ = SQRT(COEF*S1*S1+1.5D0*COEC*S2C) 96 : C 97 : C----- calcul de EPSVPT en fonction du parametre A : 98 : C A : INT(F4) = 0 99 : C Log10(A) : INT(F4) = 1 100 : C loi a 3 parametres : INT(F4) = 2 101 : C 102 : *pv pour ętre sur d'avoir la bonne valeur, on augmente un peu f4 103 : NF4=INT(F4+0.1) 104 : C 105 : C----- Loi d'Abouaf : 106 : IF (NF4.EQ.1) THEN 107 : DO I=1,NSTRS,1 108 : B = 0.D0 109 : IF (I.LE.3) B=1.D0 110 : EPSVPT(I) = 1.D0*(10.D0**A)*(SIGEQ**(RN-1.D0))*RHOR 111 : & *(COEF*S1*B+1.5D0*COEC*XX(I)) 112 : END DO 113 : C 114 : C----- Loi viscoplastique 3S a trois parametres : 115 : ELSE IF (NF4.EQ.2) THEN 116 : IF ((1.D0-RHOR1).LT.1.D-6) THEN 117 : COEE = 0.0D0 118 : COES = 1.0D0 119 : COF1 = 0.D0 120 : COF2 = 1.D0 121 : COF3 = 0.D0 122 : ELSE 123 : COEE = (3.D0*COEF)/(COEF+COEC) 124 : COSS = (COEF+COEC)**(-(RN+1.D0)/(2.D0*RN)) 125 : COF1 = (2.D0*COEF/3.D0)**((RN+1.D0)/2.D0) 126 : COF2 = (1.D0-COEE)*(COSS**(-RN)) 127 : & + 2D0*( (1.5D0)**((RN+1.D0)/2.D0) )*COF1 128 : COF3 = COEE*(COSS**(-RN))*((2.D0/3.D0)**((RN+5.D0)/4.D0)) 129 : & - 2.D0*( (1.5D0)**((RN+1.D0)/4.D0) )*COF1 130 : ENDIF 131 : DO I=1,NSTRS,1 132 : B = 0.D0 133 : IF (I.LE.3) B=1.D0 134 : EPSVPT(I) = (3.D0/2.D0)**((RN+1.D0)/2.D0)*EXP(-QSRT) 135 : & *(10.D0**A)*RHOR *((COF1*S1*(S1**2.D0)**((RN-1.D0)/2.D0) 136 : & +0.5D0*COF3*S1*S2C*(S1**2.D0*S2C)**((RN-3.D0)/4.D0)) 137 : & *B+(COF2*S2C**(RN/2.D0)+0.5D0*COF3*S1**2.D0*S2C**0.5D0 138 : & *(S1**2.D0*S2C)**((RN-3.D0)/4.D0))*XX(I)) 139 : END DO 140 : C 141 : C----- Loi viscoplastique par defaut 142 : ELSE 143 : DO I=1,NSTRS,1 144 : B = 0.D0 145 : IF (I.LE.3) B=1.D0 146 : EPSVPT(I) = 1.D0*A*(SIGEQ**(RN-1.D0))*EXP(-QSRT)*RHOR 147 : & *(COEF*S1*B+1.5D0*COEC*XX(I)) 148 : END DO 149 : ENDIF 150 : C 151 : C----- calcul de VARPT(1) derivee de EPSE 152 : PPT = SQRT(2.D0/3.D0*PROCON(EPSVPT,EPSVPT,NSTRS)) 153 : VARPT(1) = PPT 154 : C--------------------------------------------------------------------- 155 : C 156 : RETURN 157 : END 158 : 159 : 160 : 161 :
© Cast3M 2003 - All rights reserved.
Disclaimer