1 : C UO2VE SOURCE STRU 06/03/20 21:16:14 5340 2 : SUBROUTINE UO2VE(NCRIT,JCRIT,SIGMA,W,WMAX,WRUPT,SMAX, 3 : & BTR,XLTR,XINVL,BILIN,SBILI,PRECIE,PRECIZ,NN, 4 : & NN0,NC,NCA,FC,FC2,PENTE,PENTE2,LEBIL,ISING,KERRE) 5 : C---------------------------------------------------------------------- 6 : C ECOULEMENT MODELE UO2 (OTTOSEN ET GATT_MONERIE) 7 : C PREPARATION DU TERRAIN POUR LE PAS SUIVANT 8 : C---------------------------------------------------------------------- 9 : C 10 : C ENTREES 11 : C ------- 12 : C NCRIT = NOMBRE DE NOUVELLES NON LINEARITES (FISSURATION) 13 : C JCRIT(20) = TYPE DES NOUVELLES NON LINEARITES RENCONTREES 14 : C POUR LA FISSURATION 15 : C SIGMA(6) = CONTR. A LA FIN DE LA DUREE CALCULEE 16 : C W(3) = OUVERTURES DES FISSURES A LA FIN DE LA DUREE CALCULEE 17 : C WMAX(3) = OUVERTURES MAX. DES FISS. A LA FIN DE LA DUREE CALCULEE 18 : C WRUPT(3) = OUVERTURES CONDITIONNANT LA RUPTURE 19 : C SMAX(3) = CONTR. CORRESPONDANT A WMAX 20 : C BTR = PARAMETRE DE FERMETURE 21 : C XLTR(3) = LIMITES EN TRACTION POUR LA FISSURATION 22 : C XINVL(3) = PARAMETRES DE TAILLE 23 : C BILIN(3) = OUVERTURES DEFINISSANT LE CHANGEMENT DE PENTE EN CAS DE 24 : C RELATION BILINEAIRE ENTRE CONTRAINTE ET OUVERTURE 25 : C SBILI(3) = CONTR. CORRESPONDANT A BILIN 26 : C PRECIE = PRECISION POUR TESTS SUR OUVERTURES DE FISSURES 27 : C PRECIZ = PRECISION POUR TESTS SUR CONTRAINTES 28 : C NN(NCA) = NUMEROS DES DIRECTIONS DE FISS. OU UN CRIT. EST ATTEINT 29 : C POUR L ETAT INITIAL 30 : C NC = (=3) NBR. TOTAL DE DIRECTIONS DE FISS. POSSIBLES 31 : C 32 : C SORTIES 33 : C ------- 34 : C FC(NC) = CRITERE DE FISSURATION 35 : C FC2(NC) = DEUXIEME CRITERE DE FISSURATION SI BIFURCATION POSSIBLE 36 : C PENTE(NC) = PENTE DE LA DROITE DE FISSURATION CORRESPONDANT A FCRIT 37 : C PENTE2(NC) = PENTE DE LA DROITE DE FISS. CORRESPONDANT A FCRIT2 38 : C LEBIL(NC) = COMPRESSION/TRACTION 39 : C ISING(NC) = SINGULARITES CORRESPONDANT A UNE BIFURCATION 40 : C---------------------------------------------------------------------- 41 : C 42 : IMPLICIT INTEGER(I-N) 43 : IMPLICIT REAL*8(A-H,O-Z) 44 : -INC CCOPTIO 45 : C 46 : PARAMETER (XZER=0.D0) 47 : C 48 : DIMENSION JCRIT(*),SIGMA(*),W(*),WMAX(*),WRUPT(*),SMAX(*),XLTR(*), 49 : & XINVL(*),BILIN(*),SBILI(*),NN(*),NN0(*),LEBIL(*), 50 : & ISING(*) 51 : DIMENSION FC(*),FC2(*),PENTE(*),PENTE2(*) 52 : C 53 : KERRE=0 54 : C 55 : C 56 : DO I=1,3 57 : IF(ABS(WMAX(I)-WRUPT(I)).LT.WRUPT(I)*PRECIE) THEN 58 : IF(W(I).EQ.WMAX(I)) THEN 59 : W(I)=WRUPT(I) 60 : ENDIF 61 : WMAX(I)=WRUPT(I) 62 : SMAX(I)=XZER 63 : ENDIF 64 : ENDDO 65 : C 66 : C 67 : C --------------------------------------------- 68 : C INCREMENT DES DEF. DE FISS. (OUV. + GLISS.) 69 : C --------------------------------------------- 70 : * DO I=1,6 71 : * IF(ABS(DEFPLA(I)).LT.PRECIE*1.D-2) DEFPLA(I)=XZER 72 : * ENDDO 73 : C 74 : C 75 : C************************* 76 : IF (NCRIT.NE.0) THEN 77 : C************************* 78 : DO IC=1,NCRIT 79 : JC=JCRIT(IC) 80 : C 81 : C CAS DE REOUVERTURE 82 : C ------------------ 83 : IF(JC.GE.10.AND.JC.LE.12) THEN 84 : KC=JC-9 85 : SIGMA(KC)=XZER 86 : W(KC)=BTR*MIN(WMAX(KC),WRUPT(KC)) 87 : C 88 : C CAS DE FERMETURE 89 : C ---------------- 90 : ELSE IF(JC.GE.13.AND.JC.LE.15) THEN 91 : KC=JC-12 92 : SIGMA(KC)=XZER 93 : W(KC)=BTR*MIN(WMAX(KC),WRUPT(KC)) 94 : C 95 : C CAS DE RUPTURE 96 : C -------------- 97 : ELSE IF(JC.GE.7.AND.JC.LE.9) THEN 98 : KC=JC-6 99 : SIGMA(KC)=XZER 100 : W(KC)=WRUPT(KC) 101 : WMAX(KC)=WRUPT(KC) 102 : SMAX(KC)=XZER 103 : C 104 : C CAS DE FISSURATION 105 : C ------------------ 106 : ELSE IF(JC.GE.1.AND.JC.LE.3) THEN 107 : IF(ISING(JC).EQ.2) THEN 108 : SIGMA(JC)=SMAX(JC) 109 : W(JC)=WMAX(JC) 110 : ENDIF 111 : C 112 : C CAS DE POURSUITE DE LA FISSURATION 113 : C ---------------------------------- 114 : ELSE IF(JC.GE.4.AND.JC.LE.6) THEN 115 : KC=JC-3 116 : SIGMA(KC)=SMAX(KC) 117 : W(KC)=WMAX(KC) 118 : ENDIF 119 : ENDDO 120 : C********** 121 : ENDIF 122 : C********** 123 : C 124 : C 125 : C 126 : IF (NCA.NE.0) THEN 127 : DO IC=1,NCA 128 : JC=NN(IC) 129 : C 130 : C FISSURATION AVEC ISING=2 131 : C ------------------------ 132 : IF((JC.GE.1.AND.JC.LE.3).AND.ISING(JC).EQ.2) THEN 133 : SIGMA(JC)=SMAX(JC) 134 : W(JC)=WMAX(JC) 135 : ENDIF 136 : C 137 : ENDDO 138 : ENDIF 139 : C 140 : C 141 : C 142 : DO IC=1,NC 143 : NN0(IC)=IC 144 : ENDDO 145 : C 145 : CALL UO2EC(NC,NC,NN0,SIGMA,W,WMAX,SMAX,BILIN,WRUPT,BTR, 146 : & XLTR,XINVL,SBILI,FC,FC2,PENTE,PENTE2,LEBIL,ISING, 147 : & PRECIE,PRECIZ,KERRE) 148 : IF(KERRE.NE.0) THEN 149 : PRINT *, ' UO2VE - APRES UO2EC KERRE=',KERRE 150 : RETURN 151 : ENDIF 152 : C 153 : DO IC=1,NC 154 : IF(FC(IC).GT.PRECIZ.OR.FC2(IC).GT.PRECIZ) THEN 155 : PRINT *,'CRITERE DEPASSE NUMERO ',IC 156 : PRINT *,'CRITERE FC = ',FC(IC) 157 : PRINT *,'CRITERE FC2 = ',FC2(IC) 158 : KERRE=989 159 : RETURN 160 : ENDIF 161 : ENDDO 162 : C 163 : C 164 : IF(IIMPI.EQ.42) THEN 165 : WRITE(IOIMP,78000) (FC(IC),IC=1,NC) 166 : 78000 FORMAT( 2X, ' UO2VE - FC '/(4(1X,1PE12.5)/)/) 167 : ENDIF 168 : C 169 : RETURN 170 : END 171 : 172 : 173 : 174 :
© Cast3M 2003 - All rights reserved.
Disclaimer