1 : C MISTR1 SOURCE STRU 06/12/14 21:15:23 5611 2 : c Sous-Programme principal du module MISTRAL1 3 : C -------------------------------------------------------------------------- 4 : SUBROUTINE MISTR1 (T, TT,FI,FIT, SIG, 5 : & EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI, 6 : & DTA,SEQA0,RRA0,VEEA0, KPLAST, 7 : & DTI, TTP,FIP, CSIG,CEPS,CSIEP0, 8 : & PDILT,PCOEL, NDPI,NDVP,NXX,NPSI, 9 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI, 10 : & NPINCR,PINCR, DTMI,DTMAX,TMAX, 11 : & NPAS,NCHPL, EPSMAX, IES,FICH,IUM, ISTOP, 12 : & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH, 13 : & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI) 14 : C -------------------------------------------------------------------------- 15 : C Determination, pour la periode d'integration de T a T+DTI : 16 : C - des valeurs a T+DTI : 17 : C - de la temperature TT, 18 : C - du flux de neutrons rapides FI, 19 : C - de la fluence de neutrons rapides FIT, 20 : C - des contraintes SIG(IJ), 21 : C - des deformations totales EPS(IJ), 22 : C - des dilatations thermiques EPSTH(IJ), 23 : C - des deformations elastiques EPSEL(IJ), 24 : C - des deformations plastiques EPSP(IJ,IDP), 25 : C - des deformations de croissance EPSCR(IJ), 26 : C - des deformations plastiques equivalentes EPSE(IDP) et EPSEQ(IDP), 27 : C - des contraintes internes XX(IJ,IXX), 28 : C - de l'eventuelle variable PSI de durcissement du a l'irradiation, 29 : C - de KPLAST (0 : elastique, 1 : plastique), si NDPI=1, 30 : C - de leurs variations entre T et T+DTI, respectivement : 31 : C DELTT, DELFI, DELFIT, DELSIG, DELEPS, DELETH, 32 : C DELEEL, DELEP, DELECR, DELEE, DELEEQ, DELXX ,DELPSI 33 : C en fonction : 34 : C - de leurs valeurs a l'instant T (memes variables qu'a T+DTI), 35 : C - des derivees, supposees constantes durant toute la periode : 36 : C - TTP de la temperature, 37 : C - FIP du flux de neutrons rapides, 38 : C - des conditions mecaniques imposees, 39 : C supposees constantes durant toute la periode : 40 : C CSIG*SIGP + CEPS*VEPS = CSIEP0 , 41 : C - des parametres du modele, pour le materiau au point considere, fournis 42 : C dans PDILT, PCOEL, PCOHI, PECOU, PEDIR, PRVCE, PECRX, PDVDI et PCROI. 43 : C -------------------------------------------------------------------------- 44 : C IJ dans l'ordre : 11, 22, 33, 12, 13, 23 45 : C IDP compris entre 1-NDPI et NDVP 46 : C NXX = nombre de natures de contraintes internes 47 : C -------------------------------------------------------------------------- 48 : C Parametres et variable pour la determination automatique des pas 49 : C de calcul : DTMI, DTMA (=MIN(DTMAX,DTI)), TMAX, PINCR(1:NPINCR), DTA. 50 : C 51 : C Variables pour la determination des changements 52 : C d'etat plastique (KPLAST) : SEQA0, RRA0, VEEA0, DTA. 53 : C 54 : C Incrementation de : 55 : C - NPAS du nombre de pas de calcul durant la periode, 56 : C - NCHPL du nombre de changements d'etat plastique durant la periode. 57 : C 58 : C EPSMAX est le maximum autorise pour la valeur absolue 59 : C de toute composante de deformation totale. 60 : C 61 : C Si IES different de 0, stockage de valeurs par STOCM sur des fichiers 62 : C associes au canal IUM et aux NDPI+NDVP+NXX canaux suivants : 63 : C au debut de la periode d'integration, tous les ABS(IES) pas, 64 : C et en fin de periode si IES<0 (reperage de la derniere periode). 65 : C -------------------------------------------------------------------------- 66 : IMPLICIT INTEGER (I-N) 67 : IMPLICIT REAL*8 (A-H, O-Z) 68 : c std : 69 : PARAMETER ( NDVPMA = 4 , NXXMAX = 3 ) 70 : c std. 71 : PARAMETER ( NPINCM = 16) 72 : DIMENSION SIG(1:*),EPS(1:*),EPSTH(1:*),EPSEL(1:*),EPSCR(1:*), 73 : & EPSP(1:6,0:*),EPSE(0:*),EPSEQ(0:*),XX(1:6,1:*) 74 : DIMENSION PDILT(1:*),PCOEL(1:*),PCOHI(1:*),PECOU(1:*), 75 : & PEDIR(1:*),PRVCE(1:*),PECRX(1:*),PDVDI(1:*),PCROI(1:*), 76 : & PINCR(1:*) 77 : DIMENSION CSIG(1:6,1:*),CEPS(1:6,1:*),CSIEP0(1:*) 78 : c std : 79 : DIMENSION SEQ(0:NDVPMA),RR(0:NDVPMA) 80 : DIMENSION SIGP(1:6),VEPS(1:6),VEPSTH(1:6),VEPSEL(1:6),VEPSCR(1:6), 81 : & VEPSP(1:6,0:NDVPMA),VEPSE(0:NDVPMA), 82 : & EPSEQP(0:NDVPMA),XXP(1:6,1:NXXMAX) 83 : DIMENSION SIGM(1:6),EPSM(1:6),EPSTHM(1:6),EPSELM(1:6),EPSCRM(1:6), 84 : & EPSPM(1:6,0:NDVPMA),EPSEM(0:NDVPMA),EPSEQM(0:NDVPMA), 85 : & XXM(1:6,1:NXXMAX) 86 : DIMENSION SIGPM(1:6),VEPSM(1:6),VEPTHM(1:6),VEPELM(1:6), 87 : & VEPCRM(1:6),VEPSPM(1:6,0:NDVPMA),VEPSEM(0:NDVPMA), 88 : & EPEQPM(0:NDVPMA),XXPM(1:6,1:NXXMAX) 89 : DIMENSION VEPTHA(1:6),VEPELA(1:6),VEPCRA(1:6), 90 : & VEPEA(0:NDVPMA),EPSEPA(0:NDVPMA),XXPA(1:6,1:NXXMAX) 91 : c std. 92 : DIMENSION FPINCR(1:NPINCM) 93 : DIMENSION DELSIG(1:*), 94 : & DELEPS(1:*),DELETH(1:*),DELEEL(1:*),DELECR(1:*), 95 : & DELEP(1:6,0:*),DELEE(0:*),DELEEQ(0:*), 96 : & DELXX(1:6,1:*) 97 : CHARACTER *(*) FICH 98 : C ------------------------------------------------------------------ 99 : IUI = IUM-2 100 : IF (NPSI.EQ.0) THEN 101 : II = 19+2*(NDPI+NDVP)+6*NXX 102 : ELSE 103 : II = 20+2*(NDPI+NDVP)+6*NXX 104 : END IF 105 : C ------------------------------------------------------------------ 106 : C Initialisations au debut de la periode d'integration 107 : C 107 : CALL INCREM (TT,FI,FIT, 108 : & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI, 109 : & TTP,FIP, 110 : & SIGP,VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,VEPSE, 111 : & EPSEQP,XXP,PSIP, 112 : & NDPI,NDVP,NXX, 113 : & 0.D0, 0, 114 : & TT,FI,FIT, 115 : & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI, 116 : & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH, 117 : & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI) 118 : NP = 0 119 : DTIR = DTI 120 : IF (NDPI.EQ.0) THEN 121 : KPLAST = 0 122 : ELSE 123 : FDTI = 1.D-2*DTI 124 : DO IPINCR = 1,NPINCR 125 : FPINCR(IPINCR) = 1.D-1*PINCR(IPINCR) 126 : END DO 127 : ICHPL = 0 128 : IF (T.GT.0.) THEN 129 : ICHPLT = 0 130 : ELSE 131 : ICHPLT = 1 132 : END IF 133 : END IF 134 : C ------------------------------------------------------------------ 135 : C Debut du pas de temps, calcul des derivees 136 : C 137 : 1 CONTINUE 137 : CALL MISTR0 (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST, 138 : & CSIG,CEPS,CSIEP0, IUI, 139 : & PDILT,PCOEL, NDPI,NDVP,NXX, 140 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI, 141 : & SEQ,RR, 142 : & SIGP, VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR, 143 : & VEPSE,EPSEQP,XXP,PSIP, 144 : & ISTOP) 145 : IF (ISTOP.NE.0) THEN 146 : RETURN 147 : END IF 148 : C ------------------------------------------------------------------ 149 : C Determination de la duree DT du pas de temps 150 : C 151 : IF (NDPI.EQ.0) THEN 152 : DTMA = MIN(DTMAX,DTIR) 152 : CALL DTPREC (DTMI,DTMA,TMAX, 153 : & VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP, 154 : & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA, 155 : & DTA,NP, NDPI,NDVP,NXX,NPSI, II, PINCR, DT) 156 : ELSE 157 : IF (ICHPLT.EQ.0) THEN 158 : DTMA = MIN(DTMAX,DTIR) 158 : CALL DTPREC (DTMI,DTMA,TMAX, 159 : & VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP, 160 : & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA, 161 : & DTA,NP, NDPI,NDVP,NXX,NPSI, II, PINCR, DT) 161 : CALL DTCHPL (KPLAST, SEQA0,SEQ(0),RRA0,RR(0), 162 : & VEEA0,VEPSE(0),DTA, DT, ICHPL) 163 : IF ((ICHPL.EQ.1).AND.(DT.LE.0.)) THEN 164 : ICHPLT = KPLAST 165 : KPLAST = 1-KPLAST 166 : NCHPL = NCHPL+1 167 : ICHPL = 0 168 : GOTO 1 169 : END IF 170 : ELSE 171 : DTMA = MIN(DTMAX,DTIR,FDTI) 171 : CALL DTPREC (0.D0,DTMA,TMAX, 172 : & VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP, 173 : & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA, 174 : & DTA,0, NDPI,NDVP,NXX,NPSI, II, FPINCR, DT) 175 : END IF 176 : END IF 177 : C ------------------------------------------------------------------ 178 : C Stockage de valeurs au debut du pas de temps 179 : C Arret eventuel sur deformation excessive 180 : C 181 : IF (IES.NE.0) THEN 182 : IF ((MOD(NPAS,ABS(IES)).EQ.0).OR.(DTIR.EQ.DTI)) THEN 182 : CALL STOCM (T, TT,FI,FIT, PSI, 183 : & SIG, VEPS, EPS,EPSTH,EPSEL,EPSCR, 184 : & SEQ,RR, VEPSE,EPSEQP,VEPSP, EPSE,EPSEQ,EPSP, 185 : & XX, NDPI,NDVP,NXX, FICH,IUM) 186 : END IF 187 : END IF 188 : DO IJ = 1,6 189 : IF (ABS(EPS(IJ)).GT.EPSMAX) THEN 190 : WRITE (IUI,601) T,IJ,EPSMAX 191 : 601 FORMAT (' Arret de MISTRAL a l''instant : ',G22.15,/, 192 : & ' car la composante ',I1,' de deformation totale',/, 193 : & ' depasse ',G22.15,' en valeur absolue') 194 : ISTOP = 1 195 : END IF 196 : END DO 197 : IF (ISTOP.NE.0) THEN 198 : RETURN 199 : END IF 199 : CALL VAM (VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP, 200 : & SEQ,RR, NDPI,NDVP,NXX,NPSI, 201 : & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA, 202 : & SEQA0,RRA0,VEEA0) 203 : C ------------------------------------------------------------------ 204 : C Valeurs au milieu du pas de temps 205 : C 206 : 2 DTM = DT/2. 206 : CALL INCREM (TT,FI,FIT, 207 : & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI, 208 : & TTP,FIP, 209 : & SIGP,VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,VEPSE, 210 : & EPSEQP,XXP,PSIP, 211 : & NDPI,NDVP,NXX, 212 : & DTM, 1, 213 : & TTM,FIM,FITM, 214 : & SIGM,EPSM,EPSTHM,EPSELM,EPSPM,EPSCRM,EPSEM, 215 : & EPSEQM,XXM,PSIM, 216 : & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH, 217 : & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI) 217 : CALL MISTR0 (TTM,TTP,FIM,FITM, SIGM, XXM,EPSEM,EPSEQM,PSIM,KPLAST, 218 : & CSIG,CEPS,CSIEP0, IUI, 219 : & PDILT,PCOEL, NDPI,NDVP,NXX, 220 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI, 221 : & SEQ,RR, 222 : & SIGPM, VEPSM,VEPTHM,VEPELM,VEPSPM,VEPCRM, 223 : & VEPSEM,EPEQPM,XXPM,PSIPM, 224 : & ISTOP) 225 : IF (ISTOP.NE.0) THEN 226 : RETURN 227 : END IF 228 : C ------------------------------------------------------------------ 229 : C Verification et diminution eventuelle de DT 230 : C 231 : DTMA = MIN(DTMAX,DTIR) 231 : CALL DTPREC (DTMI,DTMA,TMAX, 232 : & VEPTHM,VEPELM,VEPCRM,VEPSEM,EPEQPM,XXPM,FIM,PSIPM, 233 : & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA, 234 : & DTM,1, NDPI,NDVP,NXX,NPSI, II, PINCR, DTC) 235 : IF (DTC.LT.DTM) THEN 236 : DT = DTC 237 : IF (NDPI.NE.0) ICHPL = 0 238 : GO TO 2 239 : END IF 240 : C ------------------------------------------------------------------ 241 : C Valeurs a la fin du pas de temps 242 : C 242 : CALL INCREM (TT,FI,FIT, 243 : & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI, 244 : & TTP,FIP, 245 : & SIGPM,VEPSM,VEPTHM,VEPELM,VEPSPM,VEPCRM,VEPSEM, 246 : & EPEQPM,XXPM,PSIPM, 247 : & NDPI,NDVP,NXX, 248 : & DT, 2, 249 : & TT,FI,FIT, 250 : & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI, 251 : & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH, 252 : & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI) 253 : T = T+DT 254 : DTA = DT 255 : NPAS = NPAS+1 256 : NP = NPAS 257 : IF (NDPI.NE.0) THEN 258 : IF (ICHPL.EQ.0) THEN 259 : ICHPLT = 0 260 : ELSE 261 : ICHPLT = KPLAST 262 : KPLAST = 1-KPLAST 263 : NCHPL = NCHPL+1 264 : ICHPL = 0 265 : NP = 0 266 : END IF 267 : END IF 268 : DTIR = DTIR-DT 269 : IF (DTIR.GT.0.) GOTO 1 270 : C ------------------------------------------------------------------ 271 : C Stockage eventuel de resultats 272 : C a la fin de la periode d'integration 273 : C 274 : IF (IES.LT.0) THEN 274 : CALL MISTR0 (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST, 275 : & CSIG,CEPS,CSIEP0, IUI, 276 : & PDILT,PCOEL, NDPI,NDVP,NXX, 277 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI, 278 : & SEQ,RR, 279 : & SIGP, VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR, 280 : & VEPSE,EPSEQP,XXP,PSIP, 281 : & ISTOP) 282 : IF (ISTOP.NE.0) THEN 283 : RETURN 284 : END IF 284 : CALL STOCM (T, TT,FI,FIT, PSI, 285 : & SIG, VEPS, EPS,EPSTH,EPSEL,EPSCR, 286 : & SEQ,RR, VEPSE,EPSEQP,VEPSP, EPSE,EPSEQ,EPSP, 287 : & XX, NDPI,NDVP,NXX, FICH,IUM) 288 : END IF 289 : RETURN 290 : END 291 : 292 :
© Cast3M 2003 - All rights reserved.
Disclaimer