Télécharger la plaquette Cast3M


fr en

Télécharger uo2cr2.eso

Retour à la liste

Numérotation des lignes : 

uo2cr2

  1. C UO2CR2    SOURCE    STRU      06/03/20    21:15:35     5340
  2. C   responsable STRUB
  3.       SUBROUTINE UO2CR2(NC,NCA,NN,SIGMA,W,WMAX,
  4.      &      XINVL,FCRIT,PENTE,LEBIL,PRECIE,PRECIZ,KERRE)
  5. C----------------------------------------------------------------------
  6. C      ECOULEMENT MODELE UO2 (OTTOSEN ET GATT_MONERIE) SIMPLIFIE
  7. C----------------------------------------------------------------------
  8. C
  9. C  ENTREES
  10. C  -------
  11. C NC           = 3 NBR. TOTAL DE DIRECTIONS DE FISS. POSSIBLES
  12. C NCA          = NBR. DE DIRECTIONS DE FISS. OU UN CRITERE EST ATTEINT
  13. C NN(NC)       = NUMEROS DES DIRECTIONS DE FISS. OU UN CRIT. EST ATTEINT
  14. C SIGMA(6)     = CONTRAINTES
  15. C W(3)         = OUVERTURES DES FISSURES
  16. C WMAX(3)      = OUVERTURES MAXIMALES DES FISSURES
  17. C XINVL(3)     = PARAMETRES DE TAILLE
  18. C LEBIL(NC)    = COMPRESSION/TRACTION
  19. C PRECIE       = PRECISION POUR TESTS SUR OUVERTURES DE FISSURES
  20. C PRECIZ       = PRECISION POUR TESTS SUR CONTRAINTES
  21. C
  22. C  SORTIES
  23. C  -------
  24. C FCRIT(NC)    = CRITERE DE FISSURATION
  25. C PENTE(NC)    = PENTE DE LA DROITE DE FISSURATION CORRESPONDANT A FCRIT
  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,UN=1.D0)
  36. C
  37.       DIMENSION NN(*),SIGMA(*),W(*),WMAX(*)
  38.       DIMENSION XINVL(*),FCRIT(*),PENTE(*),LEBIL(*)
  39. C
  40. C     INITIALISATIONS
  41. C
  42.       KERRE=0
  43.       DO IC=1,NC
  44.         FCRIT(IC)=-1.D4*PRECIZ
  45.       ENDDO
  46. C
  47.       DO IC=1,NCA
  48.          JC=NN(IC)
  49. C
  50.          GO TO (1,1,1),JC
  51. C
  52.          PRINT *,'UO2CR2 - CAS IMPOSSIBLE JC = ',JC
  53.          KERRE=989
  54.          RETURN
  55. C
  56.   1      CONTINUE
  57. C
  58.          IF (XINVL(JC).NE.XZER) THEN
  59. C
  60.            PRECIW=PRECIE/XINVL(JC)
  61.  
  62.            IF(IIMPI.EQ.42) THEN
  63.               PRINT *,'  '
  64.               PRINT *,'UO2CR2 - JC     =',JC
  65.               PRINT *,'UO2CR2 - W      =',W(JC)
  66.               PRINT *,'UO2CR2 - LEBIL  =',LEBIL(JC)
  67.            ENDIF
  68. C
  69. C
  70. C          fissuration ou fracturation sans ouverture prealable
  71. C          ----------------------------------------------------
  72. C
  73.            IF (WMAX(JC).EQ.XZER) THEN
  74.  
  75.              IF (LEBIL(JC).EQ.1) THEN
  76.                FCRIT(JC)=SIGMA(JC)
  77.                PENTE(JC)=XZER
  78.              ENDIF
  79. C
  80.              IF (LEBIL(JC).EQ.0) THEN
  81.                FCRIT(JC)=SIGMA(JC)
  82.                IF (FCRIT(JC).GE.-PRECIZ) PENTE(JC)=XZER
  83.              ENDIF
  84. C
  85.            ELSE
  86. C
  87. C          fracturation
  88. C          ------------
  89. C
  90.              FCRIT(JC)=SIGMA(JC)
  91.              IF (ABS(W(JC)).LT.PRECIW) THEN
  92.                IF (FCRIT(JC).GE.-PRECIZ) PENTE(JC)=XZER
  93.              ELSE IF(W(JC).GT.XZER) THEN
  94.                PENTE(JC)=XZER
  95.              ENDIF
  96. C
  97.            ENDIF
  98. C
  99. C          fin tests sur ouvertures
  100. C          ------------------------
  101. C
  102.          ENDIF
  103. C
  104.       ENDDO
  105. C
  106.       IF(IIMPI.EQ.42) THEN
  107.         WRITE(IOIMP,77000) (FCRIT(IC),IC=1,3)
  108. 77000   FORMAT( 2X, ' UO2CR2 - FCRIT '/(3(1X,1PE12.5)/)/)
  109.         WRITE(IOIMP,77003) (PENTE(IC),IC=1,3)
  110. 77003   FORMAT( 2X, ' UO2CR2 - PENTE '/(3(1X,1PE12.5)/)/)
  111.       ENDIF
  112. C
  113.       RETURN
  114.       END
  115.  
  116.  

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