Esope source files

   1 : C UO2FL     SOURCE    CHAT      05/01/13    03:59:06     5004           
   2 :       SUBROUTINE UO2FL(NCA,NN,XINVL,PENTE,SIGMA,GS,
   3 :      &                 DFDS,DGDS,HDFDQ,KERRE)
   4 : C-----------------------------------------------------------------------
   5 : C                MODELE UO2 (OTTOSEN+GATT_MONERIE)
   6 : C-----------------------------------------------------------------------
   7 : C                                                                               
   8 : C  ENTREES
   9 : C  -------
  10 : C NCA          = NBR. DE CRITERES ATTEINTS
  11 : C NN(NC)       = DIRECTIONS CORRESPONDANTES AUX CRITERES ATTEINTS
  12 : C XINVL(3)     = INVERSES DES PARAMETRES DE TAILLE
  13 : C PENTE(NC)    = PENTES DES CRITERES
  14 : C                (DERIVEE DES CRIT. / DEF. DE FISS. EN OUV.)
  15 : C SIGMA(6)     = CONTRAINTES AU DEB. DU SOUS PAS
  16 : C GS(3)        = RESISTANCES AU CISAILLEMENT
  17 : C
  18 : C  SORTIES
  19 : C  -------
  20 : C DFDS(6,NC)   = DERIVEES DES CRITERES DE FISSURATION / CONTRAINTES
  21 : C DGDS(6,NC)   = MATRICE DE PASSAGE DU VECTEUR DES DEF. DE FISSURATION
  22 : C                EN OUVERTURE AU VECTEUR DE FISSURATION COMPLET
  23 : C                (OUVERTURES + GLISSEMENTS)
  24 : C HDFDQ(NC)    = DERIVEES DES CRITERES DE FISSURATION /
  25 : C                VECTEUR DES DEF. DE FISSURATION EN OUVERTURE
  26 : C KERRE        = GESTION DES ERREURS
  27 : C-----------------------------------------------------------------------
  28 : C
  29 :       IMPLICIT INTEGER(I-N)
  30 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  31 : -INC CCOPTIO                                                                    
  32 : C
  33 :       PARAMETER (XZER=0.D0)
  34 : C
  35 :       DIMENSION NN(*),DFDS(6,*),DGDS(6,*),HDFDQ(*)
  36 :       DIMENSION XINVL(*),PENTE(*),SIGMA(*),GS(*)
  37 :       DIMENSION IV1(3),IV2(3)
  38 : C
  39 :       DATA IV1/4,4,5/
  40 :       DATA IV2/5,6,6/
  41 : C
  42 :       KERRE=0
  43 : C
  44 :       IF (NCA.EQ.0) THEN
  45 :         RETURN
  46 :       ENDIF
  47 : C
  48 : C
  49 :       DO IC=1,NCA
  50 :         JC=NN(IC)
  51 :         DO I=1,6
  52 :            DFDS(I,JC)=XZER
  53 :            DGDS(I,JC)=XZER
  54 :         ENDDO
  55 :         HDFDQ(JC)=XZER
  56 : C
  57 :         GO TO (1,1,1),JC
  58 : C
  59 :         KERRE=99
  60 :         PRINT *,'&&&&&&&& UO2FL  JC = ',JC
  61 :         RETURN
  62 : C
  63 :  1      CONTINUE
  64 : C
  65 :         IC1=IV1(JC)
  66 :         IC2=IV2(JC)
  67 : C
  68 :         IF (XINVL(JC).NE.XZER) THEN
  69 :            IF(IIMPI.EQ.42) THEN
  70 :              WRITE(IOIMP,77005) XINVL(JC),PENTE(JC)
  71 : 77005        FORMAT(2X,'UO2FL - XINVL=',1PE12.5,2X,'PENTE=',1PE12.5/)
  72 :            ENDIF
  73 : C
  74 : C
  75 :            H =PENTE(JC)/XINVL(JC)
  76 :            IF(IIMPI.EQ.42) THEN
  77 :               WRITE(IOIMP,77009)  H
  78 : 77009         FORMAT(2X,'UO2FL - MODULE D ECROUISSAGE H= ',1PE12.5/)
  79 :            ENDIF
  80 :            HDFDQ(JC)=-H
  81 :            DFDS(JC,JC)=1.D0
  82 :            DGDS(JC,JC)=1.D0
  83 :            DGDS(IC1,JC)=SIGMA(IC1)/GS(JC)/2.D0
  84 :            DGDS(IC2,JC)=SIGMA(IC2)/GS(JC)/2.D0
  85 :         ENDIF
  86 :         GO TO 100
  87 : C
  88 : 100     CONTINUE
  89 :       ENDDO
  90 : C
  91 :       RETURN
  92 :       END
  93 :  
  94 :