Télécharger la plaquette Cast3M


fr en

Télécharger plast.eso

Retour à la liste

Numérotation des lignes : 

plast

  1. C PLAST     SOURCE    STRU      08/12/22    21:15:26     6246
  2. c     Sous-Programme du module MISTRAL0
  3. C     --------------------------------------------------------------------------
  4.       SUBROUTINE PLAST (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST,
  5.      &                  NDPI,NDVP, NXX,
  6.      &                  PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI,
  7.      &                  SEQ,RR, VEPSP,VEPSE,EPSEQP,EPEQP0, XXP,XXP0,
  8.      &                  PSIP,PSIP0, DSESDS,A77,B7)
  9. C     --------------------------------------------------------------------------
  10. C       Calcul :
  11. C       - de la contrainte equivalente SEQ(IDP),
  12. C       - de la contrainte seuil RR(IDP),
  13. C       - des vitesses de deformation VEPSP(IJ,IDP),
  14. C       sauf pour IDP=0 si KPLAST=1 (VEPSP(IJ,0) contient alors DSESDS(IJ)),
  15. C       - de la vitesse de deformation equivalente VEPSE(IDP),
  16. C       sauf pour IDP=0 si KPLAST=1 (VEPSE(0) contient alors 1.),
  17. C     - des termes EPSEQP(IDP) et EPEQP0(IDP) de l'expression suivante des
  18. C       derivees par rapport au temps des deformations equivalentes EPSEQ(IDP)
  19. C               dEPSEQ(IDP)/dt = EPSEQP(IDP) + EPEQP0(IDP)*VEPSE(0)
  20. C       pour les deformations (visco)plastiques de chaque nature IDP
  21. C       (IDP allant de 1-NDPI a NDVP),
  22. C     - des termes XXP(IJ,IXX) et XXP0(IJ,IXX) de l'expression suivante
  23. C       des derivees par rapport au temps des contraintes internes :
  24. C            dXX(IJ,IXX)/dt = XXP(IJ,IXX) + XXP0(IJ,IXX)*VEPSE(0)
  25. C       pour les contraintes internes de chaque nature IXX
  26. C       (IXX allant de 1 a NXX),
  27. C     - de PSIP et PSIP0 definis par : dPSI/dt = PSIP + PSIP0*VEPSE(0)
  28. C       PSI etant la variable de durcissement du a l'irradiation,
  29. C         si celle-ci existe,
  30. C     et, si KPLAST=1 :
  31. C     - des derivees partielles DSESDS(IJ) de SEQ(0) par rapport aux S(IJ)
  32. C       ou SIG(IJ),
  33. C       - du facteur A77 de VEPSE(0) et du coefficient B7 du second membre
  34. C       dans la relation de consistance,
  35. C       en fonction :
  36. C       - de la temperature TT,
  37. C       - de sa derivee par rapport au temps TTP (si KPLAST=1),
  38. C       - du flux de neutrons rapides FI,
  39. C       - de la fluence de neutrons rapides FIT,
  40. C       - des contraintes SIG(IJ),
  41. C       - des contraintes internes XX(IJ,IXX),
  42. C       - des deformations equivalentes EPSE(IDP) et EPSEQ(IDP),
  43. C     - de l'eventuelle variable PSI,
  44. C     - de la variable KPLAST (0 : elastique, 1 : plastique),
  45. C     - des parametres du modele de plasticite transmis dans
  46. C       PCOHI, PECOU, PEDIR, PRVCE, PECRX et PDVDI.
  47. C     --------------------------------------------------------------------------
  48.       IMPLICIT INTEGER (I-N)
  49.       IMPLICIT REAL*8 (A-H, O-Z)
  50. c     std :
  51.       PARAMETER ( NDVPMA = 4 , NXXMAX = 3 )
  52. c     std.
  53.       DIMENSION SIG(1:*),XX(1:6,1:*),EPSE(0:*),EPSEQ(0:*)
  54.       DIMENSION PCOHI(1:*),PECOU(1:*),PEDIR(1:*),PRVCE(1:*),
  55.      &          PECRX(1:*),PDVDI(1:*)
  56. c     std :
  57.       DIMENSION S(1:6),XXIDP(1:6),FIX0(NXXMAX)
  58.       DIMENSION PS(0:NDVPMA),XX0P(1:6),XX0P0(1:6)
  59. c     std.
  60.       DIMENSION SEQ(0:*),DSESDS(1:*),RR(0:*),
  61.      &          VEPSP(1:6,0:*),VEPSE(0:*),EPSEQP(0:*),EPEQP0(0:*)
  62.       DIMENSION XXP(1:6,1:*),XXP0(1:6,1:*)
  63.       NPCOHI = 0
  64.       NPECOU = 0
  65.       NPEDIR = 0
  66.       NPRVCE = 0
  67.       NPECRX = 0
  68.       NPDVDI = 0
  69.       NPSI = 0
  70.         DO IDP = 1-NDPI,NDVP
  71.        DO IJ = 1,6
  72.         XXIDP(IJ) = 0.
  73.        END DO
  74.        DO IXX = 1,NXX
  75.         NPECRX = NPECRX+1
  76.         FIXIP = PECRX(NPECRX)
  77.         IF (IDP.EQ.0) THEN
  78.          FIX0(IXX) = FIXIP
  79.         END IF
  80.         DO IJ = 1,6
  81.          XXIDP(IJ) = XXIDP(IJ)+FIXIP*XX(IJ,IXX)
  82.         END DO
  83.        END DO
  84.          DO IJ = 1,6
  85.           S(IJ) = SIG(IJ)-XXIDP(IJ)
  86.          END DO
  87.        IF ((IDP.EQ.0).AND.(KPLAST.EQ.1)) THEN
  88.         ID = 1
  89.        ELSE
  90.         ID = 0
  91.        END IF
  92.        NPDVDI = NPDVDI+1
  93.        KVDI = NINT(PDVDI(NPDVDI))
  94.        IF (KVDI.EQ.0) THEN
  95.         PSIIDP = FIT
  96.        ELSE
  97.         PSIIDP = PSI
  98.         NPSI = NPSI+1
  99.        END IF
  100.        PS(IDP) = PSIIDP
  101.        CALL HILL (TT,PSIIDP,S, ID, PCOHI, NPCOHI, SEQ(IDP),DSESDS,
  102.      &            DSEDTT,DSDPSI)
  103.        CALL ECDIRR (TT,EPSEQ(IDP),PSIIDP, ID, PEDIR, NPEDIR, RR(IDP),
  104.      &              DRRDTT,DRDEPS,DRDPSI)
  105.        IF (IDP.EQ.0) THEN
  106.           IF (KPLAST.EQ.0) THEN
  107.            DO IJ = 1,6
  108.             VEPSP(IJ,IDP) = 0.
  109.            END DO
  110.          VEPSE(IDP) = 0.
  111.           ELSE
  112.            DO IJ = 1,6
  113.             VEPSP(IJ,IDP) = DSESDS(IJ)
  114.            END DO
  115.          VEPSE(IDP) = 1.
  116.           END IF
  117.          ELSE
  118.         CALL ECOUV (SEQ(IDP),RR(IDP),EPSEQ(IDP),TT,FI,PSIIDP,
  119.      &              DSESDS, IDP, PECOU, NPECOU, VEPSP,VEPSE(IDP))
  120.          END IF
  121.         END DO
  122.       NPEDIR = 0
  123.         DO IDP = 1-NDPI,NDVP
  124.        CALL ECDIRR (TT,0.D0,PS(IDP), 0, PEDIR, NPEDIR, RR0,
  125.      &              DRRDTT,DRDEPS,DRDPSI)
  126.        CALL RVCEPE (VEPSE, TT,FI,PS(IDP), SEQ(IDP),RR(IDP),RR0,
  127.      &              EPSE(IDP),EPSEQ(IDP),IDP, NDPI,NDVP,
  128.      &              PRVCE, NPRVCE, EPSEQP(IDP),EPEQP0(IDP))
  129.         END DO
  130.       DO IXX = 1,NXX
  131.        CALL ECRXX (VEPSP,VEPSE,XX, TT,EPSEQ,FI,FIT,PSI, IXX,NXX,
  132.      &             NDPI,NDVP, NPECRX,PECRX,PDVDI, XXP,XXP0)
  133.         END DO
  134.       IF (NPSI.EQ.0) THEN
  135.        PSIP = 0.
  136.        PSIP0 = 0.
  137.       ELSE
  138.        CALL DERVDI (FI,TT,PSI,SEQ,EPSEQ,VEPSE, NDPI,NDVP,NXX,
  139.      &              PDVDI, PSIP,PSIP0)
  140.       END IF
  141.         IF (KPLAST.EQ.1) THEN
  142.          DO IJ = 1,6
  143.           DSESDS(IJ) = VEPSP(IJ,0)
  144.          END DO
  145.        KVDI = NINT(PDVDI(1))
  146.        IF (KVDI.EQ.0) THEN
  147.         PSI0P0 = FI
  148.         PSI0P1 = 0.
  149.        ELSE
  150.         PSI0P0 = PSIP
  151.         PSI0P1 = PSIP0
  152.        END IF
  153.        A77 = -DRDEPS-(DRDPSI-DSDPSI)*PSI0P1
  154.        B7 =
  155.      &  (DRRDTT-DSEDTT)*TTP+DRDEPS*EPSEQP(0)+(DRDPSI-DSDPSI)*PSI0P0
  156.        IF (NXX.NE.0) THEN
  157.         DO IJ = 1,6
  158.          XX0P(IJ) = 0.
  159.          XX0P0(IJ) = 0.
  160.         END DO
  161.         DO IXX = 1,NXX
  162.          DO IJ = 1,6
  163.           XX0P(IJ) = XX0P(IJ)+FIX0(IXX)*XXP(IJ,IXX)
  164.           XX0P0(IJ) = XX0P0(IJ)+FIX0(IXX)*XXP0(IJ,IXX)
  165.          END DO
  166.         END DO
  167.         DO IJ = 1,3
  168.          A77 = A77-DSESDS(IJ)*XX0P0(IJ)
  169.          B7 = B7+DSESDS(IJ)*XX0P(IJ)
  170.         END DO
  171.         DO IJ = 4,6
  172.          A77 = A77-2.*DSESDS(IJ)*XX0P0(IJ)
  173.          B7 = B7+2.*DSESDS(IJ)*XX0P(IJ)
  174.         END DO
  175.          END IF
  176.         END IF
  177.         RETURN
  178.         END
  179.  
  180.  

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