1 : C PLAST SOURCE STRU 08/12/22 21:15:26 6246 2 : c Sous-Programme du module MISTRAL0 3 : C -------------------------------------------------------------------------- 4 : SUBROUTINE PLAST (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST, 5 : & NDPI,NDVP, NXX, 6 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, 7 : & SEQ,RR, VEPSP,VEPSE,EPSEQP,EPEQP0, XXP,XXP0, 8 : & PSIP,PSIP0, DSESDS,A77,B7) 9 : C -------------------------------------------------------------------------- 10 : C Calcul : 11 : C - de la contrainte equivalente SEQ(IDP), 12 : C - de la contrainte seuil RR(IDP), 13 : C - des vitesses de deformation VEPSP(IJ,IDP), 14 : C sauf pour IDP=0 si KPLAST=1 (VEPSP(IJ,0) contient alors DSESDS(IJ)), 15 : C - de la vitesse de deformation equivalente VEPSE(IDP), 16 : C sauf pour IDP=0 si KPLAST=1 (VEPSE(0) contient alors 1.), 17 : C - des termes EPSEQP(IDP) et EPEQP0(IDP) de l'expression suivante des 18 : C derivees par rapport au temps des deformations equivalentes EPSEQ(IDP) 19 : C dEPSEQ(IDP)/dt = EPSEQP(IDP) + EPEQP0(IDP)*VEPSE(0) 20 : C pour les deformations (visco)plastiques de chaque nature IDP 21 : C (IDP allant de 1-NDPI a NDVP), 22 : C - des termes XXP(IJ,IXX) et XXP0(IJ,IXX) de l'expression suivante 23 : C des derivees par rapport au temps des contraintes internes : 24 : C dXX(IJ,IXX)/dt = XXP(IJ,IXX) + XXP0(IJ,IXX)*VEPSE(0) 25 : C pour les contraintes internes de chaque nature IXX 26 : C (IXX allant de 1 a NXX), 27 : C - de PSIP et PSIP0 definis par : dPSI/dt = PSIP + PSIP0*VEPSE(0) 28 : C PSI etant la variable de durcissement du a l'irradiation, 29 : C si celle-ci existe, 30 : C et, si KPLAST=1 : 31 : C - des derivees partielles DSESDS(IJ) de SEQ(0) par rapport aux S(IJ) 32 : C ou SIG(IJ), 33 : C - du facteur A77 de VEPSE(0) et du coefficient B7 du second membre 34 : C dans la relation de consistance, 35 : C en fonction : 36 : C - de la temperature TT, 37 : C - de sa derivee par rapport au temps TTP (si KPLAST=1), 38 : C - du flux de neutrons rapides FI, 39 : C - de la fluence de neutrons rapides FIT, 40 : C - des contraintes SIG(IJ), 41 : C - des contraintes internes XX(IJ,IXX), 42 : C - des deformations equivalentes EPSE(IDP) et EPSEQ(IDP), 43 : C - de l'eventuelle variable PSI, 44 : C - de la variable KPLAST (0 : elastique, 1 : plastique), 45 : C - des parametres du modele de plasticite transmis dans 46 : C PCOHI, PECOU, PEDIR, PRVCE, PECRX et PDVDI. 47 : C -------------------------------------------------------------------------- 48 : IMPLICIT INTEGER (I-N) 49 : IMPLICIT REAL*8 (A-H, O-Z) 50 : c std : 51 : PARAMETER ( NDVPMA = 4 , NXXMAX = 3 ) 52 : c std. 53 : DIMENSION SIG(1:*),XX(1:6,1:*),EPSE(0:*),EPSEQ(0:*) 54 : DIMENSION PCOHI(1:*),PECOU(1:*),PEDIR(1:*),PRVCE(1:*), 55 : & PECRX(1:*),PDVDI(1:*) 56 : c std : 57 : DIMENSION S(1:6),XXIDP(1:6),FIX0(NXXMAX) 58 : DIMENSION PS(0:NDVPMA),XX0P(1:6),XX0P0(1:6) 59 : c std. 60 : DIMENSION SEQ(0:*),DSESDS(1:*),RR(0:*), 61 : & VEPSP(1:6,0:*),VEPSE(0:*),EPSEQP(0:*),EPEQP0(0:*) 62 : DIMENSION XXP(1:6,1:*),XXP0(1:6,1:*) 63 : NPCOHI = 0 64 : NPECOU = 0 65 : NPEDIR = 0 66 : NPRVCE = 0 67 : NPECRX = 0 68 : NPDVDI = 0 69 : NPSI = 0 70 : DO IDP = 1-NDPI,NDVP 71 : DO IJ = 1,6 72 : XXIDP(IJ) = 0. 73 : END DO 74 : DO IXX = 1,NXX 75 : NPECRX = NPECRX+1 76 : FIXIP = PECRX(NPECRX) 77 : IF (IDP.EQ.0) THEN 78 : FIX0(IXX) = FIXIP 79 : END IF 80 : DO IJ = 1,6 81 : XXIDP(IJ) = XXIDP(IJ)+FIXIP*XX(IJ,IXX) 82 : END DO 83 : END DO 84 : DO IJ = 1,6 85 : S(IJ) = SIG(IJ)-XXIDP(IJ) 86 : END DO 87 : IF ((IDP.EQ.0).AND.(KPLAST.EQ.1)) THEN 88 : ID = 1 89 : ELSE 90 : ID = 0 91 : END IF 92 : NPDVDI = NPDVDI+1 93 : KVDI = NINT(PDVDI(NPDVDI)) 94 : IF (KVDI.EQ.0) THEN 95 : PSIIDP = FIT 96 : ELSE 97 : PSIIDP = PSI 98 : NPSI = NPSI+1 99 : END IF 100 : PS(IDP) = PSIIDP 100 : CALL HILL (TT,PSIIDP,S, ID, PCOHI, NPCOHI, SEQ(IDP),DSESDS, 101 : & DSEDTT,DSDPSI) 101 : CALL ECDIRR (TT,EPSEQ(IDP),PSIIDP, ID, PEDIR, NPEDIR, RR(IDP), 102 : & DRRDTT,DRDEPS,DRDPSI) 103 : IF (IDP.EQ.0) THEN 104 : IF (KPLAST.EQ.0) THEN 105 : DO IJ = 1,6 106 : VEPSP(IJ,IDP) = 0. 107 : END DO 108 : VEPSE(IDP) = 0. 109 : ELSE 110 : DO IJ = 1,6 111 : VEPSP(IJ,IDP) = DSESDS(IJ) 112 : END DO 113 : VEPSE(IDP) = 1. 114 : END IF 115 : ELSE 115 : CALL ECOUV (SEQ(IDP),RR(IDP),EPSEQ(IDP),TT,FI,PSIIDP, 116 : & DSESDS, IDP, PECOU, NPECOU, VEPSP,VEPSE(IDP)) 117 : END IF 118 : END DO 119 : NPEDIR = 0 120 : DO IDP = 1-NDPI,NDVP 120 : CALL ECDIRR (TT,0.D0,PS(IDP), 0, PEDIR, NPEDIR, RR0, 121 : & DRRDTT,DRDEPS,DRDPSI) 121 : CALL RVCEPE (VEPSE, TT,FI,PS(IDP), SEQ(IDP),RR(IDP),RR0, 122 : & EPSE(IDP),EPSEQ(IDP),IDP, NDPI,NDVP, 123 : & PRVCE, NPRVCE, EPSEQP(IDP),EPEQP0(IDP)) 124 : END DO 125 : DO IXX = 1,NXX 125 : CALL ECRXX (VEPSP,VEPSE,XX, TT,EPSEQ,FI,FIT,PSI, IXX,NXX, 126 : & NDPI,NDVP, NPECRX,PECRX,PDVDI, XXP,XXP0) 127 : END DO 128 : IF (NPSI.EQ.0) THEN 129 : PSIP = 0. 130 : PSIP0 = 0. 131 : ELSE 131 : CALL DERVDI (FI,TT,PSI,SEQ,EPSEQ,VEPSE, NDPI,NDVP,NXX, 132 : & PDVDI, PSIP,PSIP0) 133 : END IF 134 : IF (KPLAST.EQ.1) THEN 135 : DO IJ = 1,6 136 : DSESDS(IJ) = VEPSP(IJ,0) 137 : END DO 138 : KVDI = NINT(PDVDI(1)) 139 : IF (KVDI.EQ.0) THEN 140 : PSI0P0 = FI 141 : PSI0P1 = 0. 142 : ELSE 143 : PSI0P0 = PSIP 144 : PSI0P1 = PSIP0 145 : END IF 146 : A77 = -DRDEPS-(DRDPSI-DSDPSI)*PSI0P1 147 : B7 = 148 : & (DRRDTT-DSEDTT)*TTP+DRDEPS*EPSEQP(0)+(DRDPSI-DSDPSI)*PSI0P0 149 : IF (NXX.NE.0) THEN 150 : DO IJ = 1,6 151 : XX0P(IJ) = 0. 152 : XX0P0(IJ) = 0. 153 : END DO 154 : DO IXX = 1,NXX 155 : DO IJ = 1,6 156 : XX0P(IJ) = XX0P(IJ)+FIX0(IXX)*XXP(IJ,IXX) 157 : XX0P0(IJ) = XX0P0(IJ)+FIX0(IXX)*XXP0(IJ,IXX) 158 : END DO 159 : END DO 160 : DO IJ = 1,3 161 : A77 = A77-DSESDS(IJ)*XX0P0(IJ) 162 : B7 = B7+DSESDS(IJ)*XX0P(IJ) 163 : END DO 164 : DO IJ = 4,6 165 : A77 = A77-2.*DSESDS(IJ)*XX0P0(IJ) 166 : B7 = B7+2.*DSESDS(IJ)*XX0P(IJ) 167 : END DO 168 : END IF 169 : END IF 170 : RETURN 171 : END 172 :
© Cast3M 2003 - All rights reserved.
Disclaimer