1 : C UO2VE2 SOURCE STRU 06/03/20 21:16:10 5340 2 : SUBROUTINE UO2VE2(NCRIT,JCRIT,SIGMA,W,WMAX, 3 : & XLTR,XINVL,PRECIE,PRECIZ,NN, 4 : & NN0,NC,NCA,FC,FC2,PENTE,LEBIL,KERRE) 5 : C---------------------------------------------------------------------- 6 : C ECOULEMENT MODELE UO2 (OTTOSEN ET GATT_MONERIE) SIMPLIFIE 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 XLTR(3) = LIMITES EN TRACTION POUR LA FISSURATION 19 : C XINVL(3) = PARAMETRES DE TAILLE 20 : C PRECIE = PRECISION POUR TESTS SUR OUVERTURES DE FISSURES 21 : C PRECIZ = PRECISION POUR TESTS SUR CONTRAINTES 22 : C NN(NCA) = NUMEROS DES DIRECTIONS DE FISS. OU UN CRIT. EST ATTEINT 23 : C POUR L ETAT INITIAL 24 : C NC = (=3) NBR. TOTAL DE DIRECTIONS DE FISS. POSSIBLES 25 : C 26 : C SORTIES 27 : C ------- 28 : C FC(NC) = CRITERE DE FISSURATION 29 : C FC2(NC) = DEUXIEME CRITERE DE FISSURATION SI BIFURCATION POSSIBLE 30 : C PENTE(NC) = PENTE DE LA DROITE DE FISSURATION CORRESPONDANT A FCRIT 31 : C LEBIL(NC) = COMPRESSION/TRACTION 32 : C 33 : C REMARQUE 34 : C -------- 35 : C NN0(NC) = TABLEAU DE TRAVAIL CONTENAT TOUTES LES DIRECTIONS 36 : C DE FISS. SOIT 1 2 ET 3 37 : C---------------------------------------------------------------------- 38 : C 39 : IMPLICIT INTEGER(I-N) 40 : IMPLICIT REAL*8(A-H,O-Z) 41 : -INC CCOPTIO 42 : C 43 : PARAMETER (XZER=0.D0) 44 : C 45 : DIMENSION JCRIT(*),SIGMA(*),W(*),WMAX(*),XLTR(*), 46 : & XINVL(*),NN(*),NN0(*),LEBIL(*) 47 : DIMENSION FC(*),FC2(*),PENTE(*) 48 : C 49 : KERRE=0 50 : C 51 : C 52 : DO I=1,3 53 : IF (XINVL(I).NE.XZER) THEN 54 : PRECIW=PRECIE/XINVL(I) 55 : IF(ABS(WMAX(I)).LT.PRECIW) THEN 56 : IF(W(I).EQ.WMAX(I)) W(I)=XZER 57 : WMAX(I)=XZER 58 : ENDIF 59 : ENDIF 60 : ENDDO 61 : C 62 : C 63 : C 64 : C************************* 65 : IF (NCRIT.NE.0) THEN 66 : C************************* 67 : DO IC=1,NCRIT 68 : JC=JCRIT(IC) 69 : C 70 : C CAS DE REOUVERTURE 71 : C ------------------ 72 : IF(JC.GE.4.AND.JC.LE.6) THEN 73 : KC=JC-3 74 : SIGMA(KC)=XZER 75 : W(KC)=XZER 76 : C 77 : C CAS DE FERMETURE 78 : C ---------------- 79 : ELSE IF(JC.GE.7.AND.JC.LE.9) THEN 80 : KC=JC-6 81 : SIGMA(KC)=XZER 82 : W(KC)=XZER 83 : ENDIF 84 : ENDDO 85 : C********** 86 : ENDIF 87 : C********** 88 : C 89 : C 90 : C 91 : C 92 : DO IC=1,NC 93 : NN0(IC)=IC 94 : ENDDO 95 : C 95 : CALL UO2EC2(NC,NC,NN0,SIGMA,W,WMAX, 96 : & XLTR,XINVL,FC,FC2,PENTE,LEBIL, 97 : & PRECIE,PRECIZ,KERRE) 98 : IF(KERRE.NE.0) THEN 99 : PRINT *, ' UO2VE2 - APRES UO2EC2 KERRE=',KERRE 100 : RETURN 101 : ENDIF 102 : C 103 : DO IC=1,NC 104 : IF(FC2(IC).GT.PRECIZ) THEN 105 : PRINT *,' UO2EV2 - CRITERE FISSURATION VIOLE ',IC 106 : PRINT *,' UO2EV2 - CRITERE FC2 = ',FC2(IC) 107 : KERRE=989 108 : RETURN 109 : ENDIF 110 : IF( (FC2(IC).LT.-PRECIZ).AND.(FC(IC).GT.PRECIZ) ) THEN 111 : PRINT *,' UO2EV2 - CRITERE FRACTURATION VIOLE ',IC 112 : PRINT *,' UO2EV2 - CRITERE FC = ',FC(IC) 113 : KERRE=989 114 : RETURN 115 : ENDIF 116 : ENDDO 117 : C 118 : C 119 : IF(IIMPI.EQ.42) THEN 120 : WRITE(IOIMP,78000) (FC(IC),IC=1,NC) 121 : 78000 FORMAT( 2X, ' UO2VE2 - FC '/(4(1X,1PE12.5)/)/) 122 : WRITE(IOIMP,78001) (FC2(IC),IC=1,NC) 123 : 78001 FORMAT( 2X, ' UO2VE2 - FC2 '/(4(1X,1PE12.5)/)/) 124 : ENDIF 125 : C 126 : RETURN 127 : END 128 :
© Cast3M 2003 - All rights reserved.
Disclaimer