Télécharger plast.eso

Retour à la liste

Numérotation des lignes :

  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