Télécharger mistr0.eso

Retour à la liste

Numérotation des lignes :

  1. C MISTR0 SOURCE STRU 06/12/14 21:15:21 5611
  2. c Sous-Programme principal du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE MISTR0 (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST,
  5. & CSIG,CEPS,CSIEP0, IUI,
  6. & PDILT,PCOEL, NDPI,NDVP,NXX,
  7. & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI,
  8. & SEQ,RR,
  9. & SIGP, VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,
  10. & VEPSE,EPSEQP,XXP,PSIP,
  11. & ISTOP)
  12. C --------------------------------------------------------------------------
  13. C Determination, a un instant quelconque :
  14. C - des contraintes equivalentes SEQ(IDP),
  15. C - des contraintes seuil RR(IDP),
  16. C - des derivees des contraintes par rapport au temps SIGP(IJ),
  17. C - des vitesses de deformation totale VEPS(IJ),
  18. C - des vitesses de dilatation thermique VEPSTH(IJ),
  19. C - des vitesses de deformation elastique VEPSEL(IJ),
  20. C - des vitesses de deformations plastiques VEPSP(IJ,IDP),
  21. C - des vitesses de croissance VEPSCR(IJ),
  22. C - des vitesses de deformations plastiques equivalentes VEPSE(IDP),
  23. C - des derivees des deformations plastiques equivalentes EPSEQP(IDP),
  24. C - des derivees des contraintes internes XXP(IJ,IXX),
  25. C - de la derivee PSIP de la variable de durcissement du a l'irradiation
  26. C PSI, si elle existe,
  27. C en fonction :
  28. C - de la temperature TT,
  29. C - de sa derivee par rapport au temps TTP,
  30. C - du flux de neutrons rapides FI,
  31. C - de la fluence de neutrons rapides FIT,
  32. C - des contraintes SIG(IJ),
  33. C - des contraintes internes XX(IJ,IXX),
  34. C - des deformations equivalentes EPSE(IDP) et EPSEQ(IDP),
  35. C - de l'eventuelle variable PSI de durcissement du a l'irradiation,
  36. C - de la variable KPLAST (0 : elastique, 1 : plastique),
  37. C - des conditions mecaniques imposees : CSIG*SIGP + CEPS*VEPS = CSIEP0
  38. C - des parametres du modele, pour le materiau au point considere, transmis
  39. C dans PDILT, PCOEL, PCOHI, PECOU, PEDIR, PRVCE, PECRX, PDVDI et PCROI.
  40. C --------------------------------------------------------------------------
  41. C IJ dans l'ordre : 11, 22, 33, 12, 13, 23
  42. C IDP compris entre 1-NDPI et NDVP
  43. C IXX compris entre 1 et NXX
  44. C --------------------------------------------------------------------------
  45. C Resolution du systeme d'equations suivant :
  46. C - conditions mecaniques imposees :
  47. C CSIG*SIGP + CEPS*VEPS = CSIEP0
  48. C - equations issues des lois de deformation du materiau :
  49. C NDVP
  50. C VEPS(IJ) = VEPSTH(IJ) + VEPSEL(IJ) + SOMME ( VEPSP(IJ,IDP) ) + VEPSCR(IJ)
  51. C IDP=1-KPLAST
  52. C avec : VEPSEL = VEPEL0 + SOEL*SIGP
  53. C et, si KPLAST=1 : VEPSP(IJ,0) = DSESDS(IJ)*VEPSE(0)
  54. C soit : VEPS = VEPS0 + SOEL*SIGP (+VEPSE(0)*DSESDS)
  55. C - relation de consistance, si KPLAST=1 :
  56. C SOMME ( DSESDS(IJ)*SIGP(IJ) ) + A77*VEPSE(0) = B7
  57. C IJ
  58. C --------------------------------------------------------------------------
  59. IMPLICIT INTEGER (I-N)
  60. IMPLICIT REAL*8 (A-H, O-Z)
  61. c std :
  62. PARAMETER ( NDVPMA = 4 , NXXMAX = 3 )
  63. c std.
  64. DIMENSION SIG(1:*),XX(1:6,1:*),EPSE(0:*),EPSEQ(0:*)
  65. DIMENSION CSIG(1:6,1:*),CEPS(1:6,1:*),CSIEP0(1:*)
  66. DIMENSION PDILT(1:*),PCOEL(1:*),PCOHI(1:*),PECOU(1:*),
  67. & PEDIR(1:*),PRVCE(1:*),PECRX(1:*),PDVDI(1:*),PCROI(1:*)
  68. DIMENSION SOEL(1:6,1:6),VEPEL0(1:6)
  69. c std :
  70. DIMENSION VEPS0(1:6),DSESDS(1:6),EPEQP0(0:NDVPMA),
  71. & XXP0(1:6,1:NXXMAX)
  72. c std.
  73. DIMENSION A0(1:6,1:6),B0(1:6),B00(1:6,1:1),
  74. & A1(1:7,1:7),B1(1:7,1:1),ISOL(1:1)
  75. DIMENSION SIGP0(1:6,1:1),SIPVE(1:7,1:1)
  76. DIMENSION SEQ(0:*),RR(0:*),
  77. & SIGP(1:*),VEPS(1:*),VEPSTH(1:*),VEPSEL(1:*),
  78. & VEPSP(1:6,0:*),VEPSCR(1:*),
  79. & VEPSE(0:*),EPSEQP(0:*),XXP(1:6,1:*)
  80. CALL DILTH (TT,TTP, PDILT, VEPSTH)
  81. CALL ELAST (TT,TTP,SIG, PCOEL, SOEL,VEPEL0)
  82. CALL PLAST (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST,
  83. & NDPI,NDVP, NXX,
  84. & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI,
  85. & SEQ,RR, VEPSP,VEPSE,EPSEQP,EPEQP0, XXP,XXP0,
  86. & PSIP,PSIP0, DSESDS,A77,B7)
  87. CALL CROIS (TT,FI,FIT, PCROI, VEPSCR)
  88. DO L = 1,6
  89. VEPS0(L) = VEPSTH(L)+VEPEL0(L)+VEPSCR(L)
  90. DO IDP = 1,NDVP
  91. VEPS0(L) = VEPS0(L)+VEPSP(L,IDP)
  92. END DO
  93. END DO
  94. CALL PRODMA(CEPS,SOEL,6,6,6, A0)
  95. CALL PRODMA(CEPS,VEPS0,6,6,1, B0)
  96. IF (KPLAST.EQ.0) THEN
  97. DO L = 1,6
  98. DO K = 1,6
  99. A0(K,L) = CSIG(K,L)+A0(K,L)
  100. END DO
  101. B00(L,1) = CSIEP0(L)-B0(L)
  102. END DO
  103. CALL GAUSS (A0,B00,6,1, SIGP0, ISOL)
  104. DO IJ = 1,6
  105. SIGP(IJ) = SIGP0(IJ,1)
  106. END DO
  107. ELSE
  108. DO L = 1,6
  109. DO K = 1,6
  110. A1(K,L) = CSIG(K,L)+A0(K,L)
  111. END DO
  112. B1(L,1) = CSIEP0(L)-B0(L)
  113. END DO
  114. CALL PRODMA(CEPS,DSESDS,6,6,1, B0)
  115. DO IJ = 1,3
  116. A1(7,IJ) = B0(IJ)
  117. A1(IJ,7) = DSESDS(IJ)
  118. END DO
  119. DO IJ = 4,6
  120. A1(7,IJ) = B0(IJ)
  121. A1(IJ,7) = 2.*DSESDS(IJ)
  122. END DO
  123. A1(7,7) = A77
  124. B1(7,1) = B7
  125. CALL GAUSS (A1,B1,7,1, SIPVE, ISOL)
  126. VEPSE(0) = SIPVE(7,1)
  127. DO IJ = 1,6
  128. SIGP(IJ) = SIPVE(IJ,1)
  129. VEPSP(IJ,0) = DSESDS(IJ)*VEPSE(0)
  130. END DO
  131. END IF
  132. IF (ISOL(1).NE.1) THEN
  133. WRITE (IUI,600)
  134. WRITE (IUI,601) ISOL(1)
  135. 601 FORMAT (' Arret au sous-programme MISTR0 :',/,
  136. & ' 0 : aucune solution',/,
  137. & ' 2 : infinite de solutions',/,2X,I1)
  138. ISTOP = 1
  139. RETURN
  140. END IF
  141. CALL PRODMA(SOEL,SIGP,6,6,1, VEPSEL)
  142. DO IJ = 1,6
  143. VEPS(IJ) = VEPS0(IJ)+VEPSEL(IJ)
  144. IF (KPLAST.EQ.1) VEPS(IJ) = VEPS(IJ)+VEPSP(IJ,0)
  145. VEPSEL(IJ) = VEPEL0(IJ)+VEPSEL(IJ)
  146. END DO
  147. IF (KPLAST.EQ.1) THEN
  148. DO IDP = 0,NDVP
  149. EPSEQP(IDP) = EPSEQP(IDP)+EPEQP0(IDP)*VEPSE(0)
  150. END DO
  151. DO IXX = 1,NXX
  152. DO IJ = 1,6
  153. XXP(IJ,IXX) = XXP(IJ,IXX)+XXP0(IJ,IXX)*VEPSE(0)
  154. END DO
  155. END DO
  156. PSIP = PSIP+PSIP0*VEPSE(0)
  157. END IF
  158. RETURN
  159. 600 FORMAT (/,'iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii',/)
  160. END
  161.  
  162.  
  163.  

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