Esope source files

   1 : C UO2CE2    SOURCE    STRU      06/03/20    21:15:31     5340           
   2 :       SUBROUTINE UO2CE2(MC,MM,SIGMA,DXV1,W,XLTR,XINVL,
   3 :      &                 NFISSU,NVF,FCRIT,VF,YOUN,PRECIZ,JRESU,KERRE)
   4 : C----------------------------------------------------------------------
   5 : C             MODELE UO2 (OTTOSEN ET GATT_MONERIE) SIMPLIFIE
   6 : C----------------------------------------------------------------------
   7 : C
   8 : C SIGNIFICATION DU CONTENU DU TABLEAU MM
   9 : C --------------------------------------
  10 : C       1   NOUVELLE FISSURATION SELON LA DIRECTION 1
  11 : C       2   NOUVELLE FISSURATION SELON LA DIRECTION 2
  12 : C       3   NOUVELLE FISSURATION SELON LA DIRECTION 3
  13 : C       4   REOUVERTURE SELON LA DIRECTION 1
  14 : C       5   REOUVERTURE SELON LA DIRECTION 2
  15 : C       6   REOUVERTURE SELON LA DIRECTION 3
  16 : C       7   FERMETURE SELON LA DIRECTION 1
  17 : C       8   FERMETURE SELON LA DIRECTION 2
  18 : C       9   FERMETURE SELON LA DIRECTION 3
  19 : C
  20 : C  ENTREES
  21 : C  -------
  22 : C MM(20)       = TYPES DES CRIT. DE FISS. SUSCEPTIBLES D ETRE ATTEINTS
  23 : C MC           = NBR. DE CRITERES DE FISS. SUSCEPTIBLES D ETRE ATTEINTS
  24 : C SIGMA(6)     = CONTRAINTES
  25 : C DXV1(3)      = INCREMENT DES DEFORMATIONS DE FISSURATION (OUV.)
  26 : C W(3)         = OUVERTURES DES FISSURES
  27 : C XLTR(3)      = LIMITES EN TRACTION POUR LA FISSURATION
  28 : C XINVL(3)     = PARAMETRES DE TAILLE
  29 : C NFISSU       = NOMBRE DE FISSURES
  30 : C NVF          = NBR DE DIRECTIONS IMPOSEES POUVANT DEVENIR
  31 : C                DES DIRECTIONS DE FISSURATION
  32 : C YOUN         = MODULE D YOUNG
  33 : C PRECIZ       = PRECISION POUR TESTS SUR CONTRAINTES
  34 : C
  35 : C  SORTIES
  36 : C  -------
  37 : C FCRIT(20)    = CRIT. DE FISS. SUSCEPTIBLES D ETRE ATTEINTS
  38 : C VF(3,3)      = VECTEURS DES DIRECTIONS DE FISSURATION
  39 : C                NON NULS SEULEMENT SI NOUVELLES FISSURATIONS POSSIBLES
  40 : C                ET CRITERES CORRESPONDANTS ATTEINTS
  41 : C JRESU(3)     = 1 2 OU 3 SELON QUE LA DIRECTION POTENTIELLE DE FISS. 
  42 : C                EST CONNUE, INCONNUE EN PLAN OU INCONNUE EN 3D
  43 : C KERRE        = GESTION DES ERREURS
  44 : C----------------------------------------------------------------------
  45 : C
  46 :       IMPLICIT INTEGER(I-N)
  47 :       IMPLICIT REAL*8(A-H,O-Z)      
  48 : -INC CCOPTIO
  49 : C
  50 :       PARAMETER (XZER=0.D0)
  51 : C
  52 :       DIMENSION MM(*),SIGMA(*),DXV1(*)
  53 :       DIMENSION W(*),XLTR(*),XINVL(*)
  54 :       DIMENSION FCRIT(*),VF(3,*),JRESU(*)
  55 :       DIMENSION FCF(3)
  56 : C
  57 : C -----------------------
  58 : C     INITIALISATIONS
  59 : C -----------------------
  60 :       KERRE=0
  61 :       DO IC=1,15
  62 :         FCRIT(IC)=-1.D4*PRECIZ
  63 :       ENDDO
  63 :       CALL IANUL(JRESU,3)
  64 : C
  65 : C
  66 :       LAPPEL=0
  67 :       DO IC=1,MC
  68 :          JC=MM(IC)
  69 :          IF(JC.GE.1.AND.JC.LE.3) THEN
  70 :             LAPPEL=1
  71 :             GO TO 200
  72 :          ENDIF
  73 :       ENDDO
  74 : C
  75 :  200  CONTINUE
  76 : C
  77 : C --- pas d appel a OTTOFI en cas de pre-fissuration
  78 : C
  79 :       IF(XINVL(1)*XINVL(2)*XINVL(3).EQ.0.D0.AND.LAPPEL.EQ.1) THEN
  79 :           CALL OTTOFI(SIGMA,NFISSU,NVF,XLTR,XINVL,
  80 :      &                VF,JRESU,FCF,PRECIZ,KERRE)
  81 :           IF(KERRE.NE.0) RETURN
  82 :            IF(IIMPI.EQ.42) THEN
  83 :               WRITE(IOIMP,76600)
  84 : 76600         FORMAT(// 2X, ' UO2CE2 - APRES OTTOFI '  /)
  85 :               WRITE(IOIMP,76601) (JRESU(IC),IC=1,3)
  86 : 76601         FORMAT(2X, ' JRESU ',2X,3I3/)
  87 :             ENDIF
  88 :       ENDIF
  89 : C
  90 : C
  91 :       DO IC=1,MC
  92 :          JC=MM(IC)
  93 : C
  94 :          GO TO (1,1,1,4,4,4,7,7,7),JC
  95 : C
  96 :          PRINT *,'UO2CE2 - CAS IMPOSSIBLE JC = ',JC
  97 :          KERRE=989
  98 :          RETURN
  99 : C
 100 :  1       CONTINUE
 101 : C
 102 : C       nouvelle fissure selon une direction
 103 : C       ------------------------------------
 104 :            IF(JRESU(JC).NE.0) THEN
 105 :               FCRIT(JC)=FCF(JC)
 106 :            ENDIF
 107 :          GO TO 100
 108 : C
 109 :  4       CONTINUE
 110 :          KC=JC-3
 111 : C
 112 : C        reouverture selon la direction KC ( 1, 2 OU 3 )
 113 : C        -----------------------------------------------
 114 :          FCRIT(JC)= SIGMA(KC)
 115 :          GO TO 100
 116 : C
 117 :  7       CONTINUE
 118 :          KC=JC-6
 119 : C
 120 : C        refermeture selon la direction KC ( 1, 2 OU 3 )
 121 : C        -----------------------------------------------
 122 :          WAUX=W(KC)+DXV1(KC)/XINVL(KC)
 123 :          FCRIT(JC)=YOUN*(XZER-WAUX)*XINVL(KC)
 124 :          GO TO 100
 125 : C
 126 : C
 127 :  100     CONTINUE
 128 : C
 129 :       ENDDO
 130 : C
 131 : C     IF(IIMPI.EQ.42) THEN
 132 : C       WRITE(IOIMP,77000)
 133 : C77000   FORMAT(// 2X, ' EN SORTIE DE UO2CE2 ' /)
 134 : C       DO IC=1,MC
 135 : C         JC=MM(IC)
 136 : C         WRITE(IOIMP,77001) IC,JC,FCRIT(JC)
 137 : C77001     FORMAT( 2X, ' UO2CE2 - IC=',I3,2X,'JC=',I3,2X,
 138 : C    &      'FCRIT=',1PE12.5/)
 139 : C       ENDDO
 140 : C     ENDIF
 141 : C
 142 :       RETURN
 143 :       END
 144 :