1 : C RVCEPE SOURCE STRU 08/12/22 21:15:31 6246 2 : c Sous-Programme du module MISTRAL0 3 : C -------------------------------------------------------------------------- 4 : SUBROUTINE RVCEPE (VEPSE, TT,FI,PSI, SEQ,RR,RR0,EPSE,EPSEQ, 5 : & IDP, NDPI,NDVP, PRVCE, NPRVCE, EPSEQP,EPEQP0) 6 : C -------------------------------------------------------------------------- 7 : C Calcul des termes EPSEQP et EPEQP0 de l'expression suivante de la 8 : C derivee par rapport au temps de la deformation equivalente EPSEQ 9 : C relative aux deformations plastiques de nature IDP : 10 : C dEPSEQ/dt = EPSEQP + ( EPEQP0*VEPSE(0) ) 11 : C en fonction : 12 : C - de la vitesse de deformation equivalente VEPSE(IDP), 13 : C - de la temperature TT, 14 : C - du flux de neutrons rapides FI, 15 : C - de la variable de durcissement du a l'irradiation PSI *, 16 : C - de la contrainte equivalente SEQ *, 17 : C - des contraintes seuil eventuelles RR et RR0(EPSEQ=0) *, 18 : C - des deformations equivalentes EPSE et EPSEQ *, 19 : C (* relatives aux deformations plastiques de nature IDP) 20 : C - des vitesses de deformation equivalente VEPSE(IDPP) pour les autres 21 : C natures de deformations (allant de 1-NDPI a NDVP), 22 : C - des parametres du materiau transmis dans PRVCE, 23 : C NPRVCE (incremente ici) etant en entree le nombre de ceux 24 : C correspondant aux deformations de natures IDPP < IDP. 25 : C -------------------------------------------------------------------------- 26 : IMPLICIT INTEGER (I-N) 27 : IMPLICIT REAL*8 (A-H, O-Z) 28 : PARAMETER ( NPTTM = 5 , NPPSIM = 6 ) 29 : PARAMETER ( NPEDIM = 13+6*NPTTM+NPPSIM ) 30 : c (NPEDIM = 49) 31 : PARAMETER ( NVM=2 , NPATDM = 14 , NPFIM = 9 ) 32 : PARAMETER ( NPERVM = 11+NVM*(NPEDIM+2*NPATDM+NPFIM+11) ) 33 : c (NPERVM = 205) 34 : DIMENSION VEPSE(0:*) 35 : DIMENSION PRVCE(1:*) 36 : DIMENSION PERV(1:NPERVM) 37 : C ------------------------------------------------------------------ 38 : IF (IDP.EQ.0) THEN 39 : EPSEQP = 0. 40 : EPEQP0 = 1. 41 : ELSE 42 : EPSEQP = VEPSE(IDP) 43 : EPEQP0 = 0. 44 : END IF 45 : C ----------------------- 46 : C Restauration eventuelle 47 : C 48 : NPRVCE = NPRVCE+1 49 : KR = NINT(PRVCE(NPRVCE)) 50 : IF (KR.NE.0) THEN 51 : NPRVCE = NPRVCE+1 52 : NPERV = NINT(PRVCE(NPRVCE)) 53 : DO IP = 1,NPERV 54 : PERV(IP) = PRVCE(NPRVCE+IP) 55 : END DO 56 : NPRVCE = NPRVCE+NPERV 57 : IF (KR.EQ.1) THEN 58 : S = EPSEQ 59 : R = 0. 60 : ELSE 61 : S = RR 62 : R = RR0 63 : END IF 64 : IF (S.GT.R) THEN 64 : CALL ECOREV (S,R,EPSE,TT,FI,PSI, PERV, EPEQPR) 65 : EPSEQP = EPSEQP-EPEQPR 66 : END IF 67 : END IF 68 : C ----------------------- 69 : C Vieillissement eventuel 70 : C 71 : NPRVCE = NPRVCE+1 72 : KV = NINT(PRVCE(NPRVCE)) 73 : IF (KV.NE.0) THEN 74 : NPRVCE = NPRVCE+1 75 : NPERV = NINT(PRVCE(NPRVCE)) 76 : DO IP = 1,NPERV 77 : PERV(IP) = PRVCE(NPRVCE+IP) 78 : END DO 79 : NPRVCE = NPRVCE+NPERV 79 : CALL ECOREV (SEQ,0.D0,EPSEQ,TT,FI,PSI, PERV, EPEQPV) 80 : EPSEQP = EPSEQP+EPEQPV 81 : END IF 82 : C ------------------- 83 : C Couplages eventuels 84 : C 85 : NPRVCE = NPRVCE+1 86 : KC = NINT(PRVCE(NPRVCE)) 87 : IF (KC.NE.0) THEN 88 : DO IDPP = 1-NDPI,NDVP 89 : IF (IDPP.NE.IDP) THEN 90 : NPRVCE = NPRVCE+1 91 : IF (IDPP.EQ.0) THEN 92 : EPEQP0 = PRVCE(NPRVCE) 93 : ELSE 94 : EPSEQP = EPSEQP+PRVCE(NPRVCE)*VEPSE(IDPP) 95 : END IF 96 : END IF 97 : END DO 98 : END IF 99 : C ------------------------------------------------------------------ 100 : RETURN 101 : END 102 :
© Cast3M 2003 - All rights reserved.
Disclaimer