Télécharger ecouv.eso

Retour à la liste

Numérotation des lignes :

  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