1 : C MISTR0 SOURCE STRU 06/12/14 21:15:21 5611 2 : c Sous-Programme principal du module MISTRAL0 3 : C -------------------------------------------------------------------------- 4 : SUBROUTINE MISTR0 (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST, 5 : & CSIG,CEPS,CSIEP0, IUI, 6 : & PDILT,PCOEL, NDPI,NDVP,NXX, 7 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI, 8 : & SEQ,RR, 9 : & SIGP, VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR, 10 : & VEPSE,EPSEQP,XXP,PSIP, 11 : & ISTOP) 12 : C -------------------------------------------------------------------------- 13 : C Determination, a un instant quelconque : 14 : C - des contraintes equivalentes SEQ(IDP), 15 : C - des contraintes seuil RR(IDP), 16 : C - des derivees des contraintes par rapport au temps SIGP(IJ), 17 : C - des vitesses de deformation totale VEPS(IJ), 18 : C - des vitesses de dilatation thermique VEPSTH(IJ), 19 : C - des vitesses de deformation elastique VEPSEL(IJ), 20 : C - des vitesses de deformations plastiques VEPSP(IJ,IDP), 21 : C - des vitesses de croissance VEPSCR(IJ), 22 : C - des vitesses de deformations plastiques equivalentes VEPSE(IDP), 23 : C - des derivees des deformations plastiques equivalentes EPSEQP(IDP), 24 : C - des derivees des contraintes internes XXP(IJ,IXX), 25 : C - de la derivee PSIP de la variable de durcissement du a l'irradiation 26 : C PSI, si elle existe, 27 : C en fonction : 28 : C - de la temperature TT, 29 : C - de sa derivee par rapport au temps TTP, 30 : C - du flux de neutrons rapides FI, 31 : C - de la fluence de neutrons rapides FIT, 32 : C - des contraintes SIG(IJ), 33 : C - des contraintes internes XX(IJ,IXX), 34 : C - des deformations equivalentes EPSE(IDP) et EPSEQ(IDP), 35 : C - de l'eventuelle variable PSI de durcissement du a l'irradiation, 36 : C - de la variable KPLAST (0 : elastique, 1 : plastique), 37 : C - des conditions mecaniques imposees : CSIG*SIGP + CEPS*VEPS = CSIEP0 38 : C - des parametres du modele, pour le materiau au point considere, transmis 39 : C dans PDILT, PCOEL, PCOHI, PECOU, PEDIR, PRVCE, PECRX, PDVDI et PCROI. 40 : C -------------------------------------------------------------------------- 41 : C IJ dans l'ordre : 11, 22, 33, 12, 13, 23 42 : C IDP compris entre 1-NDPI et NDVP 43 : C IXX compris entre 1 et NXX 44 : C -------------------------------------------------------------------------- 45 : C Resolution du systeme d'equations suivant : 46 : C - conditions mecaniques imposees : 47 : C CSIG*SIGP + CEPS*VEPS = CSIEP0 48 : C - equations issues des lois de deformation du materiau : 49 : C NDVP 50 : C VEPS(IJ) = VEPSTH(IJ) + VEPSEL(IJ) + SOMME ( VEPSP(IJ,IDP) ) + VEPSCR(IJ) 51 : C IDP=1-KPLAST 52 : C avec : VEPSEL = VEPEL0 + SOEL*SIGP 53 : C et, si KPLAST=1 : VEPSP(IJ,0) = DSESDS(IJ)*VEPSE(0) 54 : C soit : VEPS = VEPS0 + SOEL*SIGP (+VEPSE(0)*DSESDS) 55 : C - relation de consistance, si KPLAST=1 : 56 : C SOMME ( DSESDS(IJ)*SIGP(IJ) ) + A77*VEPSE(0) = B7 57 : C IJ 58 : C -------------------------------------------------------------------------- 59 : IMPLICIT INTEGER (I-N) 60 : IMPLICIT REAL*8 (A-H, O-Z) 61 : c std : 62 : PARAMETER ( NDVPMA = 4 , NXXMAX = 3 ) 63 : c std. 64 : DIMENSION SIG(1:*),XX(1:6,1:*),EPSE(0:*),EPSEQ(0:*) 65 : DIMENSION CSIG(1:6,1:*),CEPS(1:6,1:*),CSIEP0(1:*) 66 : DIMENSION PDILT(1:*),PCOEL(1:*),PCOHI(1:*),PECOU(1:*), 67 : & PEDIR(1:*),PRVCE(1:*),PECRX(1:*),PDVDI(1:*),PCROI(1:*) 68 : DIMENSION SOEL(1:6,1:6),VEPEL0(1:6) 69 : c std : 70 : DIMENSION VEPS0(1:6),DSESDS(1:6),EPEQP0(0:NDVPMA), 71 : & XXP0(1:6,1:NXXMAX) 72 : c std. 73 : DIMENSION A0(1:6,1:6),B0(1:6),B00(1:6,1:1), 74 : & A1(1:7,1:7),B1(1:7,1:1),ISOL(1:1) 75 : DIMENSION SIGP0(1:6,1:1),SIPVE(1:7,1:1) 76 : DIMENSION SEQ(0:*),RR(0:*), 77 : & SIGP(1:*),VEPS(1:*),VEPSTH(1:*),VEPSEL(1:*), 78 : & VEPSP(1:6,0:*),VEPSCR(1:*), 79 : & VEPSE(0:*),EPSEQP(0:*),XXP(1:6,1:*) 79 : CALL DILTH (TT,TTP, PDILT, VEPSTH) 79 : CALL ELAST (TT,TTP,SIG, PCOEL, SOEL,VEPEL0) 79 : CALL PLAST (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST, 80 : & NDPI,NDVP, NXX, 81 : & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, 82 : & SEQ,RR, VEPSP,VEPSE,EPSEQP,EPEQP0, XXP,XXP0, 83 : & PSIP,PSIP0, DSESDS,A77,B7) 83 : CALL CROIS (TT,FI,FIT, PCROI, VEPSCR) 84 : DO L = 1,6 85 : VEPS0(L) = VEPSTH(L)+VEPEL0(L)+VEPSCR(L) 86 : DO IDP = 1,NDVP 87 : VEPS0(L) = VEPS0(L)+VEPSP(L,IDP) 88 : END DO 89 : END DO 89 : CALL PRODMA(CEPS,SOEL,6,6,6, A0) 89 : CALL PRODMA(CEPS,VEPS0,6,6,1, B0) 90 : IF (KPLAST.EQ.0) THEN 91 : DO L = 1,6 92 : DO K = 1,6 93 : A0(K,L) = CSIG(K,L)+A0(K,L) 94 : END DO 95 : B00(L,1) = CSIEP0(L)-B0(L) 96 : END DO 96 : CALL GAUSS (A0,B00,6,1, SIGP0, ISOL) 97 : DO IJ = 1,6 98 : SIGP(IJ) = SIGP0(IJ,1) 99 : END DO 100 : ELSE 101 : DO L = 1,6 102 : DO K = 1,6 103 : A1(K,L) = CSIG(K,L)+A0(K,L) 104 : END DO 105 : B1(L,1) = CSIEP0(L)-B0(L) 106 : END DO 106 : CALL PRODMA(CEPS,DSESDS,6,6,1, B0) 107 : DO IJ = 1,3 108 : A1(7,IJ) = B0(IJ) 109 : A1(IJ,7) = DSESDS(IJ) 110 : END DO 111 : DO IJ = 4,6 112 : A1(7,IJ) = B0(IJ) 113 : A1(IJ,7) = 2.*DSESDS(IJ) 114 : END DO 115 : A1(7,7) = A77 116 : B1(7,1) = B7 116 : CALL GAUSS (A1,B1,7,1, SIPVE, ISOL) 117 : VEPSE(0) = SIPVE(7,1) 118 : DO IJ = 1,6 119 : SIGP(IJ) = SIPVE(IJ,1) 120 : VEPSP(IJ,0) = DSESDS(IJ)*VEPSE(0) 121 : END DO 122 : END IF 123 : IF (ISOL(1).NE.1) THEN 124 : WRITE (IUI,600) 125 : WRITE (IUI,601) ISOL(1) 126 : 601 FORMAT (' Arret au sous-programme MISTR0 :',/, 127 : & ' 0 : aucune solution',/, 128 : & ' 2 : infinite de solutions',/,2X,I1) 129 : ISTOP = 1 130 : RETURN 131 : END IF 131 : CALL PRODMA(SOEL,SIGP,6,6,1, VEPSEL) 132 : DO IJ = 1,6 133 : VEPS(IJ) = VEPS0(IJ)+VEPSEL(IJ) 134 : IF (KPLAST.EQ.1) VEPS(IJ) = VEPS(IJ)+VEPSP(IJ,0) 135 : VEPSEL(IJ) = VEPEL0(IJ)+VEPSEL(IJ) 136 : END DO 137 : IF (KPLAST.EQ.1) THEN 138 : DO IDP = 0,NDVP 139 : EPSEQP(IDP) = EPSEQP(IDP)+EPEQP0(IDP)*VEPSE(0) 140 : END DO 141 : DO IXX = 1,NXX 142 : DO IJ = 1,6 143 : XXP(IJ,IXX) = XXP(IJ,IXX)+XXP0(IJ,IXX)*VEPSE(0) 144 : END DO 145 : END DO 146 : PSIP = PSIP+PSIP0*VEPSE(0) 147 : END IF 148 : RETURN 149 : 600 FORMAT (/,'iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii',/) 150 : END 151 : 152 :
© Cast3M 2003 - All rights reserved.
Disclaimer