Esope source files

   1 : C UO2MAJ    SOURCE    STRU      07/05/31    21:15:44     5744           
   2 :        SUBROUTINE UO2MAJ(IFOUR,NC,NCA,NDIM,NN,NGAT,XINVL,
   3 :      &                   TAU,SIG,EPSV,VAR,DX,SIG1,EPSV1,VAR1,DX1,W1,
   4 :      &                   SIGPT,EPSVPT,VARPT,DDX)
   5 : C-----------------------------------------------------------------------
   6 : C                MODELE UO2 (OTTOSEN+GATT_MONERIE)
   7 : C              MISE A JOUR DES GRANDEURS SUIVANTES:
   8 : C  - CONTRAINTES
   9 : C  - DEFORMATIONS VISCOPLASTIQUES
  10 : C  - VARIABLES INTERNES SCALAIRES DE GATT_MONERIE (SAUF BURNUP)
  11 : C  - DEFORMATIONS DUES A LA FISSURATION (OUVERTURES)
  12 : C-----------------------------------------------------------------------
  13 : C                                                                               
  14 : C  ENTREES
  15 : C  -------
  16 : C IFOUR        = INDIQUE L HYPOTHESE DU CALCUL
  17 : C NC           = 3 NBR. TOTAL DE DIRECTIONS DE FISS. POSSIBLES
  18 : C NCA          = NBR. DE DIRECTIONS DE FISS. OU UN CRITERE EST ATTEINT
  19 : C NDIM         = NCA+1 SI CP, NCA SINON
  20 : C NN(NC)       = NUMEROS DES DIRECTIONS DE FISS. OU UN CRIT. EST ATTEINT
  21 : C XINVL(3)     = PARAMETRES DE TAILLE
  22 : C TAU          = SOUS PAS D INTEGRATION
  23 : C SIG(6)       = CONTRAINTES                     AU DEB. DU SOUS PAS
  24 : C EPSV(6)      = DEFORMATIONS VISCOPLASTIQUES    AU DEB. DU SOUS PAS
  25 : C VAR(NGAT)    = VAR. INT. SCAL. DE GATT_MONERIE AU DEB. DU SOUS PAS
  26 : C DX(NC)       = DEF. DE FISSURATION (OUV.)      AU DEB. DU SOUS PAS
  27 : C SIGPT(6)     = VIT. DES CONTRAINTES
  28 : C EPSVPT(6)    = VIT. DES DEFORMATIONS VISCOPLASTIQUES
  29 : C VARPT(NGAT)  = VIT. DES VAR. INT. SCAL. DE GATT_MONERIE
  30 : C DDX(NC1)     = VIT. DES DEF. DE FISSURATION (OUV.)
  31 : C                (+ PREMIERE COMPOSANTE DE LA VIT. DE DEF. TOT. SI CP)
  32 : C
  33 : C  SORTIES
  34 : C  -------
  35 : C SIG1(6)      = CONTRAINTES                     A LA FIN DU SOUS PAS
  36 : C EPSV1(6)     = DEFORMATIONS VISCOPLASTIQUES    A LA FIN DU SOUS PAS
  37 : C VAR1(NGAT)   = VAR. INT. SCAL. DE GATT_MONERIE A LA FIN DU SOUS PAS
  38 : C                (SAUF TAUX DE COMBUSTION)
  39 : C DX1(NC)      = DEF. DE FISSURATION (OUV.)      A LA FIN DU SOUS PAS
  40 : C W1(3)        = OUVERTURES DE FISSURATION       A LA FIN DU SOUS PAS
  41 : C-----------------------------------------------------------------------
  42 :       IMPLICIT INTEGER(I-N)
  43 :        IMPLICIT REAL*8(A-H,O-Z)
  44 : C
  45 :        PARAMETER (XZER=0.D0)                                                 
  46 : C
  47 :        DIMENSION NN(*),XINVL(*),SIG(*),EPSV(*),VAR(*),DX(*)
  48 :        DIMENSION SIG1(*),EPSV1(*),VAR1(*),DX1(*),W1(*)
  49 :        DIMENSION SIGPT(*),EPSVPT(*),VARPT(*),DDX(*)
  50 : C
  50 :        CALL SHIFTD(DX,DX1,NC)
  50 :        CALL ZDANUL(W1,3)
  51 : C
  52 :        DO 10 I=1,6
  53 :          SIG1(I)=SIG(I)+TAU*SIGPT(I)
  54 :          EPSV1(I)=EPSV(I)+TAU*EPSVPT(I)
  55 :   10   CONTINUE
  56 :        IF(IFOUR.EQ.-2) SIG1(1)=XZER
  57 : C
  58 :        VAR1(1)=VAR(1)+TAU*VARPT(1)
  59 :        VAR1(3)=VAR(3)+TAU*VARPT(3)
  60 :        VAR1(4)=VAR(4)+TAU*VARPT(4)
  61 : 
  62 : C ------ Porosite
  63 :        IF (VAR1(1).LT.0.D0)  VAR1(1) = 0.D0
  64 :        IF (VAR1(1).GT.0.3D0) VAR1(1) = 0.3D0
  65 : C
  66 :       IF (NCA.NE.0) THEN
  67 :         DO 30 IJ=1,NCA
  68 :            JJ=NN(IJ)
  69 :            DX1(JJ)=DX(JJ)+TAU*DDX(IJ)
  70 :   30    CONTINUE
  71 :       ENDIF
  72 : C
  73 :       DO 40 I=1,3
  74 :         IF (XINVL(I).NE.XZER) W1(I)=DX1(I)/XINVL(I)
  75 :   40   CONTINUE
  76 : C
  77 :       RETURN
  78 :       END
  79 :  
  80 :  
  81 :