Télécharger mistr1.eso

Retour à la liste

Numérotation des lignes :

  1. C MISTR1 SOURCE STRU 06/12/14 21:15:23 5611
  2. c Sous-Programme principal du module MISTRAL1
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE MISTR1 (T, TT,FI,FIT, SIG,
  5. & EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI,
  6. & DTA,SEQA0,RRA0,VEEA0, KPLAST,
  7. & DTI, TTP,FIP, CSIG,CEPS,CSIEP0,
  8. & PDILT,PCOEL, NDPI,NDVP,NXX,NPSI,
  9. & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI,
  10. & NPINCR,PINCR, DTMI,DTMAX,TMAX,
  11. & NPAS,NCHPL, EPSMAX, IES,FICH,IUM, ISTOP,
  12. & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH,
  13. & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI)
  14. C --------------------------------------------------------------------------
  15. C Determination, pour la periode d'integration de T a T+DTI :
  16. C - des valeurs a T+DTI :
  17. C - de la temperature TT,
  18. C - du flux de neutrons rapides FI,
  19. C - de la fluence de neutrons rapides FIT,
  20. C - des contraintes SIG(IJ),
  21. C - des deformations totales EPS(IJ),
  22. C - des dilatations thermiques EPSTH(IJ),
  23. C - des deformations elastiques EPSEL(IJ),
  24. C - des deformations plastiques EPSP(IJ,IDP),
  25. C - des deformations de croissance EPSCR(IJ),
  26. C - des deformations plastiques equivalentes EPSE(IDP) et EPSEQ(IDP),
  27. C - des contraintes internes XX(IJ,IXX),
  28. C - de l'eventuelle variable PSI de durcissement du a l'irradiation,
  29. C - de KPLAST (0 : elastique, 1 : plastique), si NDPI=1,
  30. C - de leurs variations entre T et T+DTI, respectivement :
  31. C DELTT, DELFI, DELFIT, DELSIG, DELEPS, DELETH,
  32. C DELEEL, DELEP, DELECR, DELEE, DELEEQ, DELXX ,DELPSI
  33. C en fonction :
  34. C - de leurs valeurs a l'instant T (memes variables qu'a T+DTI),
  35. C - des derivees, supposees constantes durant toute la periode :
  36. C - TTP de la temperature,
  37. C - FIP du flux de neutrons rapides,
  38. C - des conditions mecaniques imposees,
  39. C supposees constantes durant toute la periode :
  40. C CSIG*SIGP + CEPS*VEPS = CSIEP0 ,
  41. C - des parametres du modele, pour le materiau au point considere, fournis
  42. C dans PDILT, PCOEL, PCOHI, PECOU, PEDIR, PRVCE, PECRX, PDVDI et PCROI.
  43. C --------------------------------------------------------------------------
  44. C IJ dans l'ordre : 11, 22, 33, 12, 13, 23
  45. C IDP compris entre 1-NDPI et NDVP
  46. C NXX = nombre de natures de contraintes internes
  47. C --------------------------------------------------------------------------
  48. C Parametres et variable pour la determination automatique des pas
  49. C de calcul : DTMI, DTMA (=MIN(DTMAX,DTI)), TMAX, PINCR(1:NPINCR), DTA.
  50. C
  51. C Variables pour la determination des changements
  52. C d'etat plastique (KPLAST) : SEQA0, RRA0, VEEA0, DTA.
  53. C
  54. C Incrementation de :
  55. C - NPAS du nombre de pas de calcul durant la periode,
  56. C - NCHPL du nombre de changements d'etat plastique durant la periode.
  57. C
  58. C EPSMAX est le maximum autorise pour la valeur absolue
  59. C de toute composante de deformation totale.
  60. C
  61. C Si IES different de 0, stockage de valeurs par STOCM sur des fichiers
  62. C associes au canal IUM et aux NDPI+NDVP+NXX canaux suivants :
  63. C au debut de la periode d'integration, tous les ABS(IES) pas,
  64. C et en fin de periode si IES<0 (reperage de la derniere periode).
  65. C --------------------------------------------------------------------------
  66. IMPLICIT INTEGER (I-N)
  67. IMPLICIT REAL*8 (A-H, O-Z)
  68. c std :
  69. PARAMETER ( NDVPMA = 4 , NXXMAX = 3 )
  70. c std.
  71. PARAMETER ( NPINCM = 16)
  72. DIMENSION SIG(1:*),EPS(1:*),EPSTH(1:*),EPSEL(1:*),EPSCR(1:*),
  73. & EPSP(1:6,0:*),EPSE(0:*),EPSEQ(0:*),XX(1:6,1:*)
  74. DIMENSION PDILT(1:*),PCOEL(1:*),PCOHI(1:*),PECOU(1:*),
  75. & PEDIR(1:*),PRVCE(1:*),PECRX(1:*),PDVDI(1:*),PCROI(1:*),
  76. & PINCR(1:*)
  77. DIMENSION CSIG(1:6,1:*),CEPS(1:6,1:*),CSIEP0(1:*)
  78. c std :
  79. DIMENSION SEQ(0:NDVPMA),RR(0:NDVPMA)
  80. DIMENSION SIGP(1:6),VEPS(1:6),VEPSTH(1:6),VEPSEL(1:6),VEPSCR(1:6),
  81. & VEPSP(1:6,0:NDVPMA),VEPSE(0:NDVPMA),
  82. & EPSEQP(0:NDVPMA),XXP(1:6,1:NXXMAX)
  83. DIMENSION SIGM(1:6),EPSM(1:6),EPSTHM(1:6),EPSELM(1:6),EPSCRM(1:6),
  84. & EPSPM(1:6,0:NDVPMA),EPSEM(0:NDVPMA),EPSEQM(0:NDVPMA),
  85. & XXM(1:6,1:NXXMAX)
  86. DIMENSION SIGPM(1:6),VEPSM(1:6),VEPTHM(1:6),VEPELM(1:6),
  87. & VEPCRM(1:6),VEPSPM(1:6,0:NDVPMA),VEPSEM(0:NDVPMA),
  88. & EPEQPM(0:NDVPMA),XXPM(1:6,1:NXXMAX)
  89. DIMENSION VEPTHA(1:6),VEPELA(1:6),VEPCRA(1:6),
  90. & VEPEA(0:NDVPMA),EPSEPA(0:NDVPMA),XXPA(1:6,1:NXXMAX)
  91. c std.
  92. DIMENSION FPINCR(1:NPINCM)
  93. DIMENSION DELSIG(1:*),
  94. & DELEPS(1:*),DELETH(1:*),DELEEL(1:*),DELECR(1:*),
  95. & DELEP(1:6,0:*),DELEE(0:*),DELEEQ(0:*),
  96. & DELXX(1:6,1:*)
  97. CHARACTER *(*) FICH
  98. C ------------------------------------------------------------------
  99. IUI = IUM-2
  100. IF (NPSI.EQ.0) THEN
  101. II = 19+2*(NDPI+NDVP)+6*NXX
  102. ELSE
  103. II = 20+2*(NDPI+NDVP)+6*NXX
  104. END IF
  105. C ------------------------------------------------------------------
  106. C Initialisations au debut de la periode d'integration
  107. C
  108. CALL INCREM (TT,FI,FIT,
  109. & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI,
  110. & TTP,FIP,
  111. & SIGP,VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,VEPSE,
  112. & EPSEQP,XXP,PSIP,
  113. & NDPI,NDVP,NXX,
  114. & 0.D0, 0,
  115. & TT,FI,FIT,
  116. & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI,
  117. & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH,
  118. & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI)
  119. NP = 0
  120. DTIR = DTI
  121. IF (NDPI.EQ.0) THEN
  122. KPLAST = 0
  123. ELSE
  124. FDTI = 1.D-2*DTI
  125. DO IPINCR = 1,NPINCR
  126. FPINCR(IPINCR) = 1.D-1*PINCR(IPINCR)
  127. END DO
  128. ICHPL = 0
  129. IF (T.GT.0.) THEN
  130. ICHPLT = 0
  131. ELSE
  132. ICHPLT = 1
  133. END IF
  134. END IF
  135. C ------------------------------------------------------------------
  136. C Debut du pas de temps, calcul des derivees
  137. C
  138. 1 CONTINUE
  139. CALL MISTR0 (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST,
  140. & CSIG,CEPS,CSIEP0, IUI,
  141. & PDILT,PCOEL, NDPI,NDVP,NXX,
  142. & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI,
  143. & SEQ,RR,
  144. & SIGP, VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,
  145. & VEPSE,EPSEQP,XXP,PSIP,
  146. & ISTOP)
  147. IF (ISTOP.NE.0) THEN
  148. RETURN
  149. END IF
  150. C ------------------------------------------------------------------
  151. C Determination de la duree DT du pas de temps
  152. C
  153. IF (NDPI.EQ.0) THEN
  154. DTMA = MIN(DTMAX,DTIR)
  155. CALL DTPREC (DTMI,DTMA,TMAX,
  156. & VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP,
  157. & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA,
  158. & DTA,NP, NDPI,NDVP,NXX,NPSI, II, PINCR, DT)
  159. ELSE
  160. IF (ICHPLT.EQ.0) THEN
  161. DTMA = MIN(DTMAX,DTIR)
  162. CALL DTPREC (DTMI,DTMA,TMAX,
  163. & VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP,
  164. & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA,
  165. & DTA,NP, NDPI,NDVP,NXX,NPSI, II, PINCR, DT)
  166. CALL DTCHPL (KPLAST, SEQA0,SEQ(0),RRA0,RR(0),
  167. & VEEA0,VEPSE(0),DTA, DT, ICHPL)
  168. IF ((ICHPL.EQ.1).AND.(DT.LE.0.)) THEN
  169. ICHPLT = KPLAST
  170. KPLAST = 1-KPLAST
  171. NCHPL = NCHPL+1
  172. ICHPL = 0
  173. GOTO 1
  174. END IF
  175. ELSE
  176. DTMA = MIN(DTMAX,DTIR,FDTI)
  177. CALL DTPREC (0.D0,DTMA,TMAX,
  178. & VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP,
  179. & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA,
  180. & DTA,0, NDPI,NDVP,NXX,NPSI, II, FPINCR, DT)
  181. END IF
  182. END IF
  183. C ------------------------------------------------------------------
  184. C Stockage de valeurs au debut du pas de temps
  185. C Arret eventuel sur deformation excessive
  186. C
  187. IF (IES.NE.0) THEN
  188. IF ((MOD(NPAS,ABS(IES)).EQ.0).OR.(DTIR.EQ.DTI)) THEN
  189. CALL STOCM (T, TT,FI,FIT, PSI,
  190. & SIG, VEPS, EPS,EPSTH,EPSEL,EPSCR,
  191. & SEQ,RR, VEPSE,EPSEQP,VEPSP, EPSE,EPSEQ,EPSP,
  192. & XX, NDPI,NDVP,NXX, FICH,IUM)
  193. END IF
  194. END IF
  195. DO IJ = 1,6
  196. IF (ABS(EPS(IJ)).GT.EPSMAX) THEN
  197. WRITE (IUI,601) T,IJ,EPSMAX
  198. 601 FORMAT (' Arret de MISTRAL a l''instant : ',G22.15,/,
  199. & ' car la composante ',I1,' de deformation totale',/,
  200. & ' depasse ',G22.15,' en valeur absolue')
  201. ISTOP = 1
  202. END IF
  203. END DO
  204. IF (ISTOP.NE.0) THEN
  205. RETURN
  206. END IF
  207. CALL VAM (VEPSTH,VEPSEL,VEPSCR,VEPSE,EPSEQP,XXP,FI,PSIP,
  208. & SEQ,RR, NDPI,NDVP,NXX,NPSI,
  209. & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA,
  210. & SEQA0,RRA0,VEEA0)
  211. C ------------------------------------------------------------------
  212. C Valeurs au milieu du pas de temps
  213. C
  214. 2 DTM = DT/2.
  215. CALL INCREM (TT,FI,FIT,
  216. & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI,
  217. & TTP,FIP,
  218. & SIGP,VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,VEPSE,
  219. & EPSEQP,XXP,PSIP,
  220. & NDPI,NDVP,NXX,
  221. & DTM, 1,
  222. & TTM,FIM,FITM,
  223. & SIGM,EPSM,EPSTHM,EPSELM,EPSPM,EPSCRM,EPSEM,
  224. & EPSEQM,XXM,PSIM,
  225. & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH,
  226. & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI)
  227. CALL MISTR0 (TTM,TTP,FIM,FITM, SIGM, XXM,EPSEM,EPSEQM,PSIM,KPLAST,
  228. & CSIG,CEPS,CSIEP0, IUI,
  229. & PDILT,PCOEL, NDPI,NDVP,NXX,
  230. & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI,
  231. & SEQ,RR,
  232. & SIGPM, VEPSM,VEPTHM,VEPELM,VEPSPM,VEPCRM,
  233. & VEPSEM,EPEQPM,XXPM,PSIPM,
  234. & ISTOP)
  235. IF (ISTOP.NE.0) THEN
  236. RETURN
  237. END IF
  238. C ------------------------------------------------------------------
  239. C Verification et diminution eventuelle de DT
  240. C
  241. DTMA = MIN(DTMAX,DTIR)
  242. CALL DTPREC (DTMI,DTMA,TMAX,
  243. & VEPTHM,VEPELM,VEPCRM,VEPSEM,EPEQPM,XXPM,FIM,PSIPM,
  244. & VEPTHA,VEPELA,VEPCRA,VEPEA,EPSEPA,XXPA,FIA,PSIPA,
  245. & DTM,1, NDPI,NDVP,NXX,NPSI, II, PINCR, DTC)
  246. IF (DTC.LT.DTM) THEN
  247. DT = DTC
  248. IF (NDPI.NE.0) ICHPL = 0
  249. GO TO 2
  250. END IF
  251. C ------------------------------------------------------------------
  252. C Valeurs a la fin du pas de temps
  253. C
  254. CALL INCREM (TT,FI,FIT,
  255. & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI,
  256. & TTP,FIP,
  257. & SIGPM,VEPSM,VEPTHM,VEPELM,VEPSPM,VEPCRM,VEPSEM,
  258. & EPEQPM,XXPM,PSIPM,
  259. & NDPI,NDVP,NXX,
  260. & DT, 2,
  261. & TT,FI,FIT,
  262. & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,EPSE,EPSEQ,XX,PSI,
  263. & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH,
  264. & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI)
  265. T = T+DT
  266. DTA = DT
  267. NPAS = NPAS+1
  268. NP = NPAS
  269. IF (NDPI.NE.0) THEN
  270. IF (ICHPL.EQ.0) THEN
  271. ICHPLT = 0
  272. ELSE
  273. ICHPLT = KPLAST
  274. KPLAST = 1-KPLAST
  275. NCHPL = NCHPL+1
  276. ICHPL = 0
  277. NP = 0
  278. END IF
  279. END IF
  280. DTIR = DTIR-DT
  281. IF (DTIR.GT.0.) GOTO 1
  282. C ------------------------------------------------------------------
  283. C Stockage eventuel de resultats
  284. C a la fin de la periode d'integration
  285. C
  286. IF (IES.LT.0) THEN
  287. CALL MISTR0 (TT,TTP,FI,FIT, SIG, XX,EPSE,EPSEQ,PSI,KPLAST,
  288. & CSIG,CEPS,CSIEP0, IUI,
  289. & PDILT,PCOEL, NDPI,NDVP,NXX,
  290. & PCOHI,PECOU,PEDIR,PRVCE,PECRX,PDVDI, PCROI,
  291. & SEQ,RR,
  292. & SIGP, VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,
  293. & VEPSE,EPSEQP,XXP,PSIP,
  294. & ISTOP)
  295. IF (ISTOP.NE.0) THEN
  296. RETURN
  297. END IF
  298. CALL STOCM (T, TT,FI,FIT, PSI,
  299. & SIG, VEPS, EPS,EPSTH,EPSEL,EPSCR,
  300. & SEQ,RR, VEPSE,EPSEQP,VEPSP, EPSE,EPSEQ,EPSP,
  301. & XX, NDPI,NDVP,NXX, FICH,IUM)
  302. END IF
  303. RETURN
  304. END
  305.  
  306.  
  307.  

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