Esope source files

   1 : C UO2EC2    SOURCE    STRU      06/03/20    21:15:44     5340           
   2 :       SUBROUTINE UO2EC2(NC,NCA,NN,SIGMA,W,WMAX,
   3 :      &    XLTR,XINVL,FCRIT,FCRIT2,PENTE,LEBIL,
   4 :      &    PRECIE,PRECIZ,KERRE)
   5 : C----------------------------------------------------------------------
   6 : C         MODELE UO2 (OTTOSEN ET GATT_MONERIE) SIMPLIFIE
   7 : C                  ENERGIE DE FISSURATION NULLE 
   8 : C             CARACTERISATION DE L ETAT DE FISSURATION
   9 : C----------------------------------------------------------------------
  10 : C
  11 : C  ENTREES
  12 : C  -------
  13 : C NC           = (3) NBR. TOTAL DE DIRECTIONS DE FISS. POSSIBLES
  14 : C NCA          = NBR. DE DIRECTIONS DE FISS. OU UN CRITERE EST ATTEINT
  15 : C                (NCA=NC=3 A L APPEL PAR UO2ET)
  16 : C NN(NC)       = NUMEROS DES DIRECTIONS DE FISS. OU UN CRIT. EST ATTEINT
  17 : C                (CONTIENT NN(I)=I I=1,3 A L APPEL PAR UO2ET)
  18 : C SIGMA(6)     = CONTRAINTES INITIALES
  19 : C W(3)         = OUVERTURES INITIALES DES FISSURES
  20 : C WMAX(3)      = OUVERTURES MAXIMALES INITIALES DES FISSURES
  21 : C XLTR(3)      = LIMITES EN TRACTION POUR LA FISSURATION
  22 : C XINVL(3)     = PARAMETRES DE TAILLE
  23 : C PRECIE       = PRECISION POUR TESTS SUR OUVERTURES DE FISSURES
  24 : C PRECIZ       = PRECISION POUR TESTS SUR CONTRAINTES
  25 : C 
  26 : C
  27 : C
  28 : C  SORTIES
  29 : C  -------
  30 : C FCRIT(NC)    = CRITERE DE FRACTURATION
  31 : C FCRIT2(NC)   = CRITERE DE FISSURATION
  32 : C PENTE(NC)    = PENTE DE LA DROITE DE FRACTURATION CORRESPONDANT A FCRIT
  33 : C LEBIL(NC)    = COMPRESSION/TRACTION
  34 : C----------------------------------------------------------------------
  35 : C
  36 :       IMPLICIT INTEGER(I-N)
  37 :       IMPLICIT REAL*8(A-H,O-Z)      
  38 : -INC CCOPTIO
  39 : C
  40 :       PARAMETER (XZER=0.D0,UNDEMI=.5D0,UN=1.D0,DEUX=2.D0,TROIS=3.D0)
  41 : C
  42 :       DIMENSION NN(*),SIGMA(*),W(*),WMAX(*)
  43 :       DIMENSION XLTR(*),XINVL(*)
  44 :       DIMENSION FCRIT(*),FCRIT2(*),PENTE(*),LEBIL(*)
  45 : C
  46 : C -----------------------
  47 : C     INITIALISATIONS
  48 : C -----------------------
  49 : C
  50 :       KERRE=0
  51 :       DO IC=1,NC
  52 :         FCRIT(IC)=-1.D4*PRECIZ
  53 :         FCRIT2(IC)=-1.D4*PRECIZ
  54 :         LEBIL(IC)=0
  55 :       ENDDO
  56 : C
  57 : C
  58 :       DO IC=1,NCA
  59 :          JC=NN(IC)
  60 : C
  61 :          GO TO (1,1,1),JC
  62 : C
  63 :          PRINT *,'UO2EC2 - CAS IMPOSSIBLE JC = ',JC
  64 :          KERRE=989
  65 :          RETURN
  66 : C
  67 :   1      CONTINUE
  68 : C-----------------
  69 : C
  70 :          IF (XINVL(JC).NE.XZER) THEN
  71 : C
  72 :            PRECIW=PRECIE/XINVL(JC)
  73 : C
  74 :            FCRIT(JC)=SIGMA(JC)
  75 :            PENTE(JC)=XZER
  76 : C
  77 :            IF(IIMPI.EQ.42) THEN
  78 :                   PRINT *,'UO2EC2 - JC     =',JC
  79 :                   PRINT *,'UO2EC2 - W      =',W(JC)
  80 :                   PRINT *,'UO2EC2 - WMAX   =',WMAX(JC)
  81 :             ENDIF
  82 : C
  83 : C
  84 : C          fissuration ou fracturation sans ouverture prealable
  85 : C          ----------------------------------------------------          
  86 : C
  87 :            IF (WMAX(JC).EQ.XZER) THEN
  88 : 
  89 :              IF (FCRIT(JC).GT.PRECIZ) THEN
  90 :                FCRIT2(JC) = SIGMA(JC) - XLTR(JC)
  91 :                LEBIL(JC)=1
  92 :              ELSE
  93 :                IF (FCRIT(JC).GE.-PRECIZ) LEBIL(JC)=1
  94 :              ENDIF
  95 : 
  96 :            ELSE
  97 : C
  98 : C          fracturation
  99 : C          ------------
 100 : C
 101 :              IF (ABS(W(JC)).LT.PRECIW) THEN
 102 :                IF (FCRIT(JC).GE.-PRECIZ) LEBIL(JC)=1
 103 :              ELSE IF(W(JC).GT.XZER) THEN
 104 :                LEBIL(JC)=1
 105 :              ENDIF             
 106 : C
 107 :            ENDIF
 108 : C
 109 : C          fin tests sur ouvertures
 110 : C          ------------------------
 111 : C
 112 :          ENDIF
 113 : C
 114 :       ENDDO
 115 : C
 116 :       IF(IIMPI.EQ.42) THEN
 117 :         WRITE(IOIMP,77000) (FCRIT(IC),IC=1,NC)
 118 : 77000   FORMAT( 2X, ' UO2EC2 - FCRIT '/(3(1X,1PE12.5)/)/)
 119 :         WRITE(IOIMP,77001) (FCRIT2(IC),IC=1,NC)
 120 : 77001   FORMAT( 2X, ' UO2EC2 - FCRIT '/(3(1X,1PE12.5)/)/)
 121 :         WRITE(IOIMP,77002) (PENTE(IC),IC=1,NC)
 122 : 77002   FORMAT( 2X, ' UO2EC2 - PENTE '/(3(1X,1PE12.5)/)/)
 123 :         WRITE(IOIMP,77004) (LEBIL(IC),IC=1,NC)
 124 : 77004   FORMAT( 2X, ' UO2EC2 - LEBIL '/(3I4)/)
 125 :       ENDIF
 126 : C
 127 :       RETURN
 128 :       END
 129 :