Télécharger la plaquette Cast3M


fr en

Télécharger uo2fl.eso

Retour à la liste

Numérotation des lignes : 

uo2fl

  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.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. C
  35.       PARAMETER (XZER=0.D0)
  36. C
  37.       DIMENSION NN(*),DFDS(6,*),DGDS(6,*),HDFDQ(*)
  38.       DIMENSION XINVL(*),PENTE(*),SIGMA(*),GS(*)
  39.       DIMENSION IV1(3),IV2(3)
  40. C
  41.       DATA IV1/4,4,5/
  42.       DATA IV2/5,6,6/
  43. C
  44.       KERRE=0
  45. C
  46.       IF (NCA.EQ.0) THEN
  47.         RETURN
  48.       ENDIF
  49. C
  50. C
  51.       DO IC=1,NCA
  52.         JC=NN(IC)
  53.         DO I=1,6
  54.            DFDS(I,JC)=XZER
  55.            DGDS(I,JC)=XZER
  56.         ENDDO
  57.         HDFDQ(JC)=XZER
  58. C
  59.         GO TO (1,1,1),JC
  60. C
  61.         KERRE=99
  62.         PRINT *,'&&&&&&&& UO2FL  JC = ',JC
  63.         RETURN
  64. C
  65.  1      CONTINUE
  66. C
  67.         IC1=IV1(JC)
  68.         IC2=IV2(JC)
  69. C
  70.         IF (XINVL(JC).NE.XZER) THEN
  71.            IF(IIMPI.EQ.42) THEN
  72.              WRITE(IOIMP,77005) XINVL(JC),PENTE(JC)
  73. 77005        FORMAT(2X,'UO2FL - XINVL=',1PE12.5,2X,'PENTE=',1PE12.5/)
  74.            ENDIF
  75. C
  76. C
  77.            H =PENTE(JC)/XINVL(JC)
  78.            IF(IIMPI.EQ.42) THEN
  79.               WRITE(IOIMP,77009)  H
  80. 77009         FORMAT(2X,'UO2FL - MODULE D ECROUISSAGE H= ',1PE12.5/)
  81.            ENDIF
  82.            HDFDQ(JC)=-H
  83.            DFDS(JC,JC)=1.D0
  84.            DGDS(JC,JC)=1.D0
  85.            DGDS(IC1,JC)=SIGMA(IC1)/GS(JC)/2.D0
  86.            DGDS(IC2,JC)=SIGMA(IC2)/GS(JC)/2.D0
  87.         ENDIF
  88.         GO TO 100
  89. C
  90. 100     CONTINUE
  91.       ENDDO
  92. C
  93.       RETURN
  94.       END
  95.  
  96.  
  97.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales