Télécharger la plaquette Cast3M


fr en

Télécharger uo2ce2.eso

Retour à la liste

Numérotation des lignes : 

uo2ce2

  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.  
  49. -INC PPARAM
  50. -INC CCOPTIO
  51. C
  52.       PARAMETER (XZER=0.D0)
  53. C
  54.       DIMENSION MM(*),SIGMA(*),DXV1(*)
  55.       DIMENSION W(*),XLTR(*),XINVL(*)
  56.       DIMENSION FCRIT(*),VF(3,*),JRESU(*)
  57.       DIMENSION FCF(3)
  58. C
  59. C -----------------------
  60. C     INITIALISATIONS
  61. C -----------------------
  62.       KERRE=0
  63.       DO IC=1,15
  64.         FCRIT(IC)=-1.D4*PRECIZ
  65.       ENDDO
  66.       CALL IANUL(JRESU,3)
  67. C
  68. C
  69.       LAPPEL=0
  70.       DO IC=1,MC
  71.          JC=MM(IC)
  72.          IF(JC.GE.1.AND.JC.LE.3) THEN
  73.             LAPPEL=1
  74.             GO TO 200
  75.          ENDIF
  76.       ENDDO
  77. C
  78.  200  CONTINUE
  79. C
  80. C --- pas d appel a OTTOFI en cas de pre-fissuration
  81. C
  82.       IF(XINVL(1)*XINVL(2)*XINVL(3).EQ.0.D0.AND.LAPPEL.EQ.1) THEN
  83.           CALL OTTOFI(SIGMA,NFISSU,NVF,XLTR,XINVL,
  84.      &                VF,JRESU,FCF,PRECIZ,KERRE)
  85.           IF(KERRE.NE.0) RETURN
  86.            IF(IIMPI.EQ.42) THEN
  87.               WRITE(IOIMP,76600)
  88. 76600         FORMAT(// 2X, ' UO2CE2 - APRES OTTOFI '  /)
  89.               WRITE(IOIMP,76601) (JRESU(IC),IC=1,3)
  90. 76601         FORMAT(2X, ' JRESU ',2X,3I3/)
  91.             ENDIF
  92.       ENDIF
  93. C
  94. C
  95.       DO IC=1,MC
  96.          JC=MM(IC)
  97. C
  98.          GO TO (1,1,1,4,4,4,7,7,7),JC
  99. C
  100.          PRINT *,'UO2CE2 - CAS IMPOSSIBLE JC = ',JC
  101.          KERRE=989
  102.          RETURN
  103. C
  104.  1       CONTINUE
  105. C
  106. C       nouvelle fissure selon une direction
  107. C       ------------------------------------
  108.            IF(JRESU(JC).NE.0) THEN
  109.               FCRIT(JC)=FCF(JC)
  110.            ENDIF
  111.          GO TO 100
  112. C
  113.  4       CONTINUE
  114.          KC=JC-3
  115. C
  116. C        reouverture selon la direction KC ( 1, 2 OU 3 )
  117. C        -----------------------------------------------
  118.          FCRIT(JC)= SIGMA(KC)
  119.          GO TO 100
  120. C
  121.  7       CONTINUE
  122.          KC=JC-6
  123. C
  124. C        refermeture selon la direction KC ( 1, 2 OU 3 )
  125. C        -----------------------------------------------
  126.          WAUX=W(KC)+DXV1(KC)/XINVL(KC)
  127.          FCRIT(JC)=YOUN*(XZER-WAUX)*XINVL(KC)
  128.          GO TO 100
  129. C
  130. C
  131.  100     CONTINUE
  132. C
  133.       ENDDO
  134. C
  135. C     IF(IIMPI.EQ.42) THEN
  136. C       WRITE(IOIMP,77000)
  137. C77000   FORMAT(// 2X, ' EN SORTIE DE UO2CE2 ' /)
  138. C       DO IC=1,MC
  139. C         JC=MM(IC)
  140. C         WRITE(IOIMP,77001) IC,JC,FCRIT(JC)
  141. C77001     FORMAT( 2X, ' UO2CE2 - IC=',I3,2X,'JC=',I3,2X,
  142. C    &      'FCRIT=',1PE12.5/)
  143. C       ENDDO
  144. C     ENDIF
  145. C
  146.       RETURN
  147.       END
  148.  
  149.  

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