Télécharger ecdirr.eso

Retour à la liste

Numérotation des lignes :

  1. C ECDIRR SOURCE STRU 08/12/22 21:15:10 6246
  2. c Sous-Programme du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE ECDIRR (TT,EPSEQ,PSI, ID, PEDIR, NPEDIR, RR,
  5. & DRRDTT,DRDEPS,DRDPSI)
  6. C --------------------------------------------------------------------------
  7. C Calcul, pour les deformations plastiques de nature IDP,
  8. C de la contrainte seuil RR
  9. C et de ses derivees partielles DRRDTT, DRDEPS et DRDPSI si ID = 1,
  10. C en fonction :
  11. C - de la temperature TT,
  12. C - de la deformation equivalente EPSEQ,
  13. C - de la variable de durcissement du a l'irradiation PSI,
  14. C et des parametres transmis dans PEDIR,
  15. C NPEDIR (incremente ici) etant en entree le nombre de ceux
  16. C correspondant aux deformations de natures IDPP < IDP.
  17. C --------------------------------------------------------------------------
  18. IMPLICIT INTEGER (I-N)
  19. IMPLICIT REAL*8 (A-H, O-Z)
  20. PARAMETER ( NPTTM = 5 , NPPSIM = 6 )
  21. PARAMETER ( NPEDIM = 13+6*NPTTM+NPPSIM )
  22. c (NPEDIM = 49)
  23. DIMENSION PEDIR(1:*)
  24. DIMENSION PECDI(1:NPEDIM)
  25. NPEDIR = NPEDIR+1
  26. KRR = NINT(PEDIR(NPEDIR+1))
  27. IF (KRR.EQ.0) THEN
  28. NPEDIR = NPEDIR+1
  29. RR = 0.D0
  30. IF (ID.EQ.1) THEN
  31. DRRDTT = 0.D0
  32. DRDEPS = 0.D0
  33. DRDPSI = 0.D0
  34. END IF
  35. ELSE IF (KRR.EQ.-1) THEN
  36. NPEDIR = NPEDIR+2
  37. RR = PEDIR(NPEDIR)
  38. IF (ID.EQ.1) THEN
  39. DRRDTT = 0.D0
  40. DRDEPS = 0.D0
  41. DRDPSI = 0.D0
  42. END IF
  43. ELSE
  44. NPECDI = NINT(PEDIR(NPEDIR))
  45. DO IP = 1,NPECDI
  46. PECDI(IP) = PEDIR(NPEDIR+IP)
  47. END DO
  48. NPEDIR = NPEDIR+NPECDI
  49. CALL FECRDI (TT,EPSEQ,PSI, ID, PECDI, RR,DRRDTT,DRDEPS,DRDPSI)
  50. END IF
  51. RETURN
  52. END
  53.  
  54.  

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