Télécharger la plaquette Cast3M


fr en

Télécharger ecouv.eso

Retour à la liste

Numérotation des lignes : 

ecouv

  1. C ECOUV     SOURCE    STRU      08/12/22    21:15:13     6246
  2. c     Sous-Programme du module MISTRAL0
  3. C     --------------------------------------------------------------------------
  4.       SUBROUTINE ECOUV (SEQ,RR,EPSEQ,TT,FI,PSI, DSESDS, IDP, PECOU,
  5.      &                  NPECOU, VEPSP,VEPSE)
  6. C     --------------------------------------------------------------------------
  7. C     Calcul, pour les deformations viscoplastiques de nature IDP :
  8. C     - de la vitesse de deformation equivalente VEPSE,
  9. C     - des vitesses de deformation VEPSP(IJ,IDP),
  10. C     en fonction :
  11. C     - de la contrainte effective equivalente SEQ,
  12. C     - de la contrainte seuil RR,
  13. C     - de la deformation equivalente EPSEQ,
  14. C     - de la temperature TT,
  15. C     - du flux de neutrons rapides FI,
  16. C     - de la variable de durcissement du a l'irradiation PSI,
  17. C     - des derivees partielles DSESDS(IJ) de SEQ par rapport aux contraintes
  18. C     "effectives" S(IJ) = SIG(IJ)-XX(IJ,IDP),
  19. C     - des parametres de la loi d'ecoulement transmis dans PECOU,
  20. C     NPECOU (incremente ici) etant en entree le nombre de ceux
  21. C     correspondant aux deformations de natures IDPP < IDP.
  22. C     --------------------------------------------------------------------------
  23.       IMPLICIT INTEGER (I-N)
  24.       IMPLICIT REAL*8 (A-H, O-Z)
  25.       PARAMETER ( NPTTM = 5 , NPPSIM = 6 )
  26.       PARAMETER ( NPEDIM = 13+6*NPTTM+NPPSIM )
  27. c     (NPEDIM = 49)
  28.       PARAMETER ( NVM=2 , NPATDM = 14 , NPFIM = 9 )
  29.       PARAMETER ( NPERVM = 11+NVM*(NPEDIM+2*NPATDM+NPFIM+11) )
  30. c     (NPERVM = 205)
  31.       DIMENSION DSESDS(1:*), PECOU(1:*)
  32.       DIMENSION PERV(1:NPERVM)
  33.         DIMENSION VEPSP(1:6,0:*)
  34.       NPECOU = NPECOU+1
  35.       NPERV = NINT(PECOU(NPECOU))
  36.         IF (SEQ.LE.RR) THEN
  37.          DO IJ = 1,6
  38.           VEPSP(IJ,IDP) = 0.
  39.          END DO
  40.        VEPSE = 0.
  41.         ELSE
  42.        DO IP = 1,NPERV
  43.         PERV(IP) = PECOU(NPECOU+IP)
  44.        END DO
  45.        CALL ECOREV (SEQ,RR,EPSEQ,TT,FI,PSI, PERV, VEPSE)
  46.        DO IJ = 1,6
  47.           VEPSP(IJ,IDP) = DSESDS(IJ)*VEPSE
  48.          END DO
  49.         END IF
  50.       NPECOU = NPECOU+NPERV
  51.         RETURN
  52.         END
  53.  
  54.  

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