Download uo2ca.eso
Back to the list
1 : C UO2CA SOURCE STRU 07/05/31 21:15:34 5744
2 : C responsable STRUB
3 : SUBROUTINE UO2CA(DT,MSOUPA,PRECIE,PRECIZ,SIG0,DEPST,EPSV,VARGA0,
4 : & NGAT,W,WMAX,SMAX,WRUPT,BILIN,SBILI,XLTR,
5 : & XINVL,WREOUV,GS,BTR,XMAT,NCOMAT,NSIMP,AAD,
6 : & TETA1,TPOINT,FI1,FPOINT,PRECIS,IB,IGAU,
7 : & NFISSU,NVF,VF,XX,IRESU,NRESU,
8 : & SIGF,DEFP,EPSVF,VARGAF,KERRE)
9 : C----------------------------------------------------------------------
10 : C ECOULEMENT MODELE UO2 (OTTOSEN ET GATT_MONERIE)
11 : C----------------------------------------------------------------------
12 : C
13 : C ENTREES
14 : C -------
15 : C DT = PAS D INTEGRATION
16 : C MSOUPA = NBR. MAX. DE SOUS PAS AUTORISES POUR LA VISCOPLASTICITE
17 : C PRECIE = PRECISION POUR TESTS SUR OUVERTURES DE FISSURES
18 : C PRECIZ = PRECISION POUR TESTS SUR CONTRAINTES
19 : C SIG0(6) = CONTR. AU DEBUT DU SOUS PAS D'INTEGRATION
20 : C DEPST(6) = INCR. DES DEFORM. TOTALES SUR LE PAS D INTEGRATION
21 : C EPSV(6) = DEFORM. VISCOPLAST. AU DEBUT DU SOUS PAS D'INTEGRATION
22 : C VARGA0(NGAT) = VAR. INT. SCAL. DE GATT_MONERIE AU DEB. DU SS PAS
23 : C NGAT = NOMBRE DE VAR. INT. SCALAIRES DU MODELE GATT_MONERIE
24 : C W(3) = OUVERTURES DES FISSURES AU DEBUT DU SOUS PAS
25 : C WMAX(3) = OUVERTURES MAXIMALES DES FISSURES AU DEBUT DU SOUS PAS
26 : C SMAX(3) = CONTR. CORRESPONDANT A WMAX
27 : C WRUPT(3) = OUVERTURES CONDITIONNANT LA RUPTURE
28 : C BILIN(3) = OUVERTURES DEFINISSANT LE CHANGEMENT DE PENTE EN CAS DE
29 : C RELATION BILINEAIRE ENTRE CONTRAINTE ET OUVERTURE
30 : C SBILI(3) = CONTR. CORRESPONDANT A BILIN
31 : C XLTR(3) = LIMITES EN TRACTION POUR LA FISSURATION
32 : C XINVL(3) = PARAMETRES DE TAILLE
33 : C GS(3) = RESISTANCES AU CISAILLEMENT
34 : C BTR = PARAMETRE DE FERMETURE
35 : C XMAT(NCOMAT) = CARACTERISTIQUES THERMOMECANIQUES DU MATERIAU
36 : C NSIMP = POINTE SUR LA CARACTERISTIQUE FACULTATIVE 'SIMP' DE XMAT
37 : C AAD = COEFFICIENT INTERVENANT DANS LE CALCUL DE L INCREMENT
38 : C DE LA DEFORMATION DE DENSIFICATION
39 : C TETA1 = TEMPERATURE AU DEBUT DU PAS D INTEGRATION
40 : C TPOINT = VITESSE DE TEMPERATURE CTE SUR LE PAS D INTEGRATION
41 : C FI1 = DENSITE DE FISSION AU DEBUT DU PAS D INTEGRATION
42 : C FPOINT = VIT. DE DENSITE DE FISSION CTE SUR LE PAS D INTEGRATION
43 : C PRECIS = PRECISION POUR LA VISCOPLASTICITE
44 : C IB = NUMERO DE L ELEMENT
45 : C IGAU = NUMERO DU POINT D INTEGRATION
46 : C NFISSU = NOMBRE DE FISSURES
47 : C NVF = NBR DE DIRECTIONS IMPOSEES POUVANT DEVENIR
48 : C DES DIRECTIONS DE FISSURATION
49 : C
50 : C SORTIES
51 : C -------
52 : C W(3) = OUVERTURES DES FISSURES A LA FIN DE LA DUREE CALCULEE
53 : C WMAX(3) = OUVERTURES MAX. DES FISS. A LA FIN DE LA DUREE CALCULEE
54 : C WREOUV(3) = LIMITES DE FERMETURE A LA FIN DE LA DUREE CALCULEE
55 : C SMAX(3) = CONTR. CORRESPONDANT A WMAX
56 : C VF(3,3) = VECTEURS DES DIRECTIONS DE FISSURATION
57 : C XX = DUREE QUI VIENT D ETRE CALCULEE
58 : C IRESU(20) = TYPE DES NOUVELLES NON LINEARITES RENCONTREES
59 : C POUR LA FISSURATION
60 : C NRESU = NOMBRE DE NOUVELLES NON LINEARITES (FISSURATION)
61 : C SIGF(6) = CONTR. A LA FIN DE LA DUREE CALCULEE
62 : C DEFP(6) = INC. DES DEFORM. INELASTIQUES (FISSURATION+VISCOPLAST.)
63 : C A LA FIN DE LA DUREE CALCULEE
64 : C EPSVF(6) = DEFORM. VISCOPLAST. A LA FIN DE LA DUREE CALCULEE
65 : C VARGAF(NGAT) = VAR. INT. SCAL. A LA FIN DE LA DUREE CALCULEE
66 : C KERRE = GESTION DES ERREURS
67 : C
68 : C RAPPEL : EN DEFO PLANES,CONTRAINTES PLANES OU AXI,
69 : C LA DIRECTION 1 EST CELLE PERPENDICULAIRE AU MAILLAGE
70 : C-----------------------------------------------------------------------
71 : C
72 : IMPLICIT INTEGER(I-N)
73 : IMPLICIT REAL*8(A-H,O-Z)
74 : -INC CCOPTIO
75 : C
76 : PARAMETER (XZER=0.D0,NC=3,NGATT=4)
77 : C
78 : DIMENSION SIG0(*),DEPST(*),EPSV(*),VARGA0(*)
79 : DIMENSION W(*),WMAX(*),SMAX(*),WRUPT(*),BILIN(*),SBILI(*),XLTR(*)
80 : DIMENSION XINVL(*),WREOUV(*)
81 : DIMENSION GS(*),XMAT(*)
82 : DIMENSION VF(3,*),IRESU(*)
83 : DIMENSION SIGF(*),DEFP(*),EPSVF(*),VARGAF(*)
84 : C
85 : DIMENSION EPSPT(6),SIGG0(6),EPSV0(6),VAR0(NGATT)
86 : DIMENSION W0(3),WMAX0(3),SMAX0(3),DX0(NC),DX00(NC),WREOU0(3)
87 : DIMENSION SIGGF(6),EPSVVF(6),VARF(NGATT),WF(3),DXF(NC),EPSFF(6)
88 : DIMENSION DGDS0(6,NC),DGDSF(6,NC),DFDS(6,NC),HDFDQ(NC)
89 : DIMENSION WMAXF(3),SMAXF(3),WREOUF(3),DXV1(NC)
90 : C
91 : DIMENSION NN(NC),FC(NC),FC2(NC),PENTE(NC),PENTE2(NC),LEBIL(NC)
92 : DIMENSION FCV(20),FCV1(20),MM(20)
93 : C
94 : DIMENSION ISING(NC),IFERM(NC),IBRUP(NC),NN0(NC)
95 : C
96 : DIMENSION JFIS(3)
97 : C
98 : DATA ITMAX/25/
99 : DATA IRMAX/5/
100 : C
101 : C--------------------
102 : C INITIALISATIONS
103 : C--------------------
104 : KERRE=0
105 : NC1=NC+1
106 : IRED=0
107 : IFLAG1=1
108 : YOUN=XMAT(1)
109 : C
110 : C --- PAS DE TEMPS CALCULE INITIAL
111 : C
112 : DTCALC=DT
113 : C
114 : C ------------ LA VITESSE DE LA DEFORMATION TOTALE EST CTE ------------
115 : C
116 : DO 10 I=1,6
117 : EPSPT(I)=DEPST(I)/DT
118 : 10 CONTINUE
119 : IF (IFOUR.EQ.-2) EPSPT(1)=XZER
120 : C
121 : C
122 : IF(IIMPI.EQ.42) THEN
123 : WRITE(IOIMP,74433) (EPSPT(I),I=1,6)
124 : 74433 FORMAT(1X,' EPSPT ',6(1X,1PE12.5)/)
125 : ENDIF
126 : C
127 : C ------------ DEFORMATION DE FISSURATION EN OUVERTURE ------------
128 : C
129 : DO I=1,3
130 : DX00(I)=W(I)*XINVL(I)
131 : ENDDO
132 : C
133 : C
134 : 888 CONTINUE
135 : C
136 : IRED=IRED+1
137 : IF(IRED.GT.IRMAX) THEN
138 : PRINT *,'&&&&&&&& UO2CA IRED = ',IRED
139 : KERRE=989
140 : RETURN
141 : ENDIF
142 : C
143 : IF(IRED.GT.1) DTCALC=0.5*DTCALC
144 : C
145 : C
146 : C
147 : DTLEFT = DTCALC
148 : TAU = DTLEFT
149 : XX = XZER
150 : C
151 : T0=TETA1
152 : FI0=FI1
153 : C
154 : C ----- CONTRAINTES - DEF. VISCOPLAST. -----
155 : C
156 : DO I=1,6
157 : SIGG0(I)=SIG0(I)
158 : EPSV0(I)=EPSV(I)
159 : ENDDO
160 : C
161 : C ------- VARIABLES INTERNES DE GATT_MONERIE -------
162 : C
163 : DO I=1,NGAT
164 : VAR0(I)=VARGA0(I)
165 : ENDDO
166 : C
167 : C ------------ VARIABLES DE FISSURATION ------------
168 : C
169 : DO I=1,3
170 : W0(I)=W(I)
171 : WMAX0(I)=WMAX(I)
172 : WREOU0(I) = BTR*MIN(WMAX(I),WRUPT(I))
173 : SMAX0(I)=SMAX(I)
174 : DX0(I)=W0(I)*XINVL(I)
175 : ENDDO
176 : C
177 : C
178 : NSSINC=0
179 : C
180 : C --------------------------------------------------------------
181 : C DEBUT DES ITERATIONS EN SSINCREMENTS / FIN SI DTLEFT = 0
182 : C --------------------------------------------------------------
183 : 1400 CONTINUE
184 : ITERU=0
185 : 1500 CONTINUE
186 : C
187 : NSSINC = NSSINC + 1
188 : IF (NSSINC.GT.MSOUPA) THEN
189 : KERRE=268
190 : PRINT *,' UO2CA - PAS DE CONVERGENCE EN VISCOPLASTICITE'
191 : RETURN
192 : ENDIF
193 : C
194 : C-----------------------------------------------------------------------
195 : C Etat de fissuration au temps XX et crit. susceptibles d etre atteints
196 : C-----------------------------------------------------------------------
196 : CALL ZDANUL(DXV1,3)
197 : C
197 : CALL UO2ET(NC,NN,SIGG0,W0,WMAX0,SMAX0,BILIN,WRUPT,BTR,
198 : & XLTR,XINVL,SBILI,PRECIE,PRECIZ,FC,FC2,
199 : & LEBIL,NFISSU,NVF,VF,FCV,PENTE,PENTE2,
200 : & DXV1,YOUN,NCA,MC,MM,ISING,IFERM,IBRUP,KERRE)
201 : IF(KERRE.NE.0) RETURN
202 : C
203 : DO I=1,NC
204 : ISING(I)=0
205 : IF(LEBIL(I).EQ.2) ISING(I)=1
206 : ENDDO
207 : C
208 : C-------------------------------------------------------------
209 : C En cas de bifurcation, il faut choisir son chemin
210 : C si la cond. ci-dessous est vraie NCA est forcement non nul
211 : C-------------------------------------------------------------
212 : IF(ISING(1)+ISING(2)+ISING(3).NE.0) THEN
212 : CALL UO2IN(ISING,IFERM,IBRUP,LEBIL,PENTE,PENTE2,NCA,NN,MC,
213 : & MM,SIGG0,GS,FC,XINVL,XLTR,WMAX0,WRUPT,PRECIZ,T0,
214 : & TPOINT,FI0,FPOINT,PRECIS,XMAT,NCOMAT,NSIMP,AAD,BTR,
215 : & EPSPT,EPSV0,VAR0,W0,DX0,NGAT,TAU,KERRE)
216 : IF(KERRE.NE.0) RETURN
217 : ENDIF
218 : C
219 : C
220 : C
221 : ITERI=0
222 : 444 CONTINUE
223 : C
224 : IF (ITERI.NE.0) THEN
225 : C --------------------
226 : IF(IIMPI.EQ.42) THEN
227 : WRITE(IOIMP,74078)
228 : 74078 FORMAT(/2X,' ON REINITIALISE LE PROCESSUS'//)
229 : ENDIF
230 : IF(ITERI.GT.NC) THEN
231 : KERRE=989
232 : PRINT *,' UO2CA - PAS DE CONVERGENCE - ITERI ',ITERI
233 : RETURN
234 : ENDIF
235 : C
236 : C calcul de PENTE pour tableau NN reduit NN(ITERI+1)UO2CR(NC,NCA,NN,SIGG0,W0,WMAX0,SMAX0,BILIN,WRUPT,BTR,
238 : & XLTR,XINVL,SBILI,FC,PENTE,LEBIL,PRECIE,PRECIZ,KERRE)
239 : IF(KERRE.NE.0) RETURN
240 : C
241 : ENDIF
242 : C -----
243 : C
244 : NDIM=NCA
245 : IF(IFOUR.EQ.-2) NDIM=NCA+1
246 : C
247 : C --- determination d un etat converge pour TAU inf. ou egal a DTCALC
248 : C
248 : CALL UO2DCN(IFLAG1,T0,TPOINT,FI0,FPOINT,PRECIS,PRECIZ,
249 : & XMAT,NCOMAT,NSIMP,AAD,BTR,GS,WRUPT,LEBIL,XINVL,PENTE,
250 : & EPSPT,SIGG0,EPSV0,VAR0,W0,WMAX0,WREOU0,DX0,
251 : & NGAT,NC1,NCA,NDIM,NN,TAU,TAUNEX,SIGGF,EPSVVF,
252 : & VARF,WF,DXF,WMAXF,WREOUF,TF,FIF,KERRE)
253 : IF (KERRE.NE.0) THEN
254 : RETURN
255 : ENDIF
256 : C
257 : C --- Viole-t-on un critere de fissuration pour l etat converge ?
258 : C
258 : CALL ZDANUL(DXV1,3)
259 : C
259 : CALL UO2CE(MC,MM,SIGG0,DXV1,W0,WMAX0,SMAX0,WRUPT,XLTR,XINVL,BTR,
260 : & NFISSU,NVF,FCV1,VF,YOUN,PRECIZ,JFIS,KERRE)
261 : IF(KERRE.NE.0) RETURN
262 : C
262 : CALL UO2XX(T0,TPOINT,FI0,FPOINT,PRECIS,PRECIZ,
263 : & XMAT,NCOMAT,NSIMP,AAD,GS,WRUPT,LEBIL,XINVL,PENTE,BTR,
264 : & YOUN,EPSPT,SIGG0,EPSV0,VAR0,W0,WMAX0,WREOU0,DX0,
265 : & SMAX0,XLTR,NFISSU,NVF,VF,MC,MM,FCV1,NGAT,NC1,NCA,NDIM,
266 : & NN,TAU,SIGGF,EPSVVF,VARF,WF,DXF,WMAXF,WREOUF,
267 : & TF,FIF,MCA,XXX,IRESU,NRESU,LERED,KERRE)
268 : IF(KERRE.NE.0) RETURN
269 : IF(LERED.EQ.1) GO TO 888
270 : C
271 : IF(XXX.LE.0.D0) THEN
272 : KERRE=989
273 : PRINT *,' UO2CA - XXX EST NEGATIF OU NUL '
274 : RETURN
275 : ENDIF
276 : C
277 : C
278 : C --- Controle de la coherence des resultats pour XXX*TAU
279 : C ------------------
280 : IF (NCA.NE.0) THEN
281 : C ------------------
282 : NCA0=NCA
283 : NCAA=0
284 : DO IJ=1,NCA
285 : NN0(IJ)=NN(IJ)
286 : JJ=NN(IJ)
287 : DXSOM=DXF(JJ)-DX0(JJ)
288 : C
289 : IF(DXSOM.GE.0.D0) THEN
290 : NCAA=NCAA+1
291 : NN(NCAA)=JJ
292 : ELSE IF(DXSOM.LT.0.D0) THEN
293 : C
294 : IF(LEBIL(JJ).EQ.1) THEN
295 : C
296 : IF(IFERM(JJ).NE.1) THEN
297 : NCAA=NCAA+1
298 : NN(NCAA)=JJ
299 : ELSE IF(IFERM(JJ).EQ.1) THEN
300 : IF(WF(JJ).GT.WREOUF(JJ)) THEN
301 : NCAA=NCAA+1
302 : NN(NCAA)=JJ
303 : ENDIF
304 : ENDIF
305 : C
306 : ENDIF
307 : ENDIF
308 : ENDDO
309 : C
310 : IF(IIMPI.EQ.42) THEN
311 : WRITE(IOIMP,77008) NCA,NCAA,(NN(IJ),IJ=1,NCAA)
312 : 77008 FORMAT(5X,'NCA=',I3,2X,'NCAA=',I3/2X,'NN ' ,8I3/)
313 : ENDIF
314 : C
315 : IF(NCAA.EQ.0) THEN
316 : IF(IIMPI.EQ.42) THEN
317 : PRINT *,' UO2CA - NCAA EST NUL '
318 : ENDIF
319 : NCA=NCAA
320 : ITERI=ITERI+1
321 : GO TO 444
322 : ENDIF
323 : C
324 : C
325 : IF(NCA.NE.NCAA) THEN
326 : IF(NCAA.GT.NCA) THEN
327 : PRINT *,' UO2CA - NCAA EST SUPERIEUR A NCA '
328 : KERRE=989
329 : RETURN
330 : ENDIF
331 : NCA=NCAA
332 : ITERI=ITERI+1
333 : GO TO 444
334 : ENDIF
335 : C -----
336 : ENDIF
337 : C -----
338 : C
339 : C -------------------------------------
340 : IF ((NCA.EQ.0).AND.(ITERI.NE.0)) THEN
341 : C -------------------------------------
342 : DO IJ=1,NCA0
343 : JJ=NN0(IJ)
344 : IF (IFERM(JJ).EQ.1) THEN
345 : SXDIF=SIGGF(JJ)-SIGG0(JJ)
346 : IF (SXDIF.GT.XZER) THEN
347 : PRINT *,'ATTENTION FOIREUX IB IGAU',IB,IGAU
348 : ITERU=ITERU+1
349 : IF(ITERU.GT.IRMAX) THEN
350 : PRINT *,'&&&&&&&& UO2CA ITERU = ',ITERU
351 : KERRE=989
352 : RETURN
353 : ENDIF
354 : TAU=0.5*TAU
355 : GOTO 1500
356 : ENDIF
357 : ENDIF
358 : ENDDO
359 : C -----
360 : ENDIF
361 : C -----
362 : C --- Fin du controle de la coherence des resultats pour XXX*TAU
363 : C
364 : C
365 : DO I=1,3
366 : IF(WF(I).GE.WMAXF(I)) THEN
367 : WMAXF(I)=WF(I)
367 : CALL OTTOFU(XINVL,XLTR,WRUPT,BTR,BILIN,
368 : & SBILI,WF,WMAXF,SMAXF,PENTE,I)
369 : ELSE
370 : SMAXF(I)=SMAX0(I)
371 : ENDIF
372 : WREOUF(I)=BTR*MIN(WMAXF(I),WRUPT(I))
373 : ENDDO
374 : C
375 : C
376 : C ---------------------
377 : C ENCORE QUELQUES TESTS
378 : C ---------------------
378 : CALL UO2VE(NRESU,IRESU,SIGGF,WF,WMAXF,WRUPT,SMAXF,
379 : & BTR,XLTR,XINVL,BILIN,SBILI,PRECIE,PRECIZ,
380 : & NN,NN0,NC,NCA,FC,FC2,PENTE,PENTE2,LEBIL,ISING,KERRE)
381 : IF(KERRE.NE.0) RETURN
382 : C
383 : C
384 : IF (MCA.NE.0) THEN
385 : XX=XX+(XXX*TAU)
386 : GOTO 777
387 : ELSE
388 : XX=XX+TAU
389 : ENDIF
390 : C
391 : C ----- MAJ CONTRAINTES - DEF. VISCOPLAST. A t0+XX -----
392 : C
393 : DO I=1,6
394 : SIGG0(I)=SIGGF(I)
395 : EPSV0(I)=EPSVVF(I)
396 : ENDDO
397 : C
398 : C ------- MAJ VARIABLES INTERNES DE GATT_MONERIE A t0+XX -------
399 : C
400 : DO I=1,NGAT
401 : VAR0(I)=VARF(I)
402 : ENDDO
403 : C
404 : C ------------ MAJ VARIABLES DE FISSURATION A t0+XX ------------
405 : C
406 : DO I=1,3
407 : W0(I)=WF(I)
408 : DX0(I)=DXF(I)
409 : WMAX0(I)=WMAXF(I)
410 : SMAX0(I)=SMAXF(I)
411 : ENDDO
412 : C
413 : T0=TF
414 : FI0=FIF
415 : C
416 : C ----------------------------------------------------
417 : C TEST DE FIN SS INCREMENTS / MISE A JOUR DE TAU
418 : C fin des boucles en ss increments si tau = dtleft
419 : C ----------------------------------------------------
420 : C
421 : IF (TAU.LT.DTLEFT) THEN
422 : DTLEFT = DTLEFT - TAU
423 : IF (TAUNEX.GT.DTLEFT) TAUNEX = DTLEFT
424 : TAU=TAUNEX
425 : GOTO 1400
426 : ENDIF
427 : C
428 : IF (ABS(TAU-DTLEFT).GT.(TAU/1000.)) THEN
429 : PRINT *, ' UO2CA - PROBLEME TAU > DTLEFT '
430 : KERRE=989
431 : RETURN
432 : ENDIF
433 : C
434 : XX=DTCALC
435 : C
436 : 777 CONTINUE
437 : C
438 : C -------------------------
439 : C MAJ DES VARIABLES A t0+XX
440 : C -------------------------
441 : C ------------ MAJ VARIABLES DE FISSURATION ------------
442 : C
443 : DO I=1,3
444 : W(I)=WF(I)
445 : WMAX(I)=WMAXF(I)
446 : SMAX(I)=SMAXF(I)
447 : WREOUV(I)=WREOUF(I)
448 : ENDDO
449 : C
450 : C ----- MAJ DEF. VISCOPLAST. -----
451 : C
452 : DO I=1,6
453 : EPSVF(I)=EPSVVF(I)
454 : ENDDO
455 : C
456 : C ----- MAJ INCREMENT DES DEF. DE FISS. (OUV. + GLISS.)
457 : C
457 : CALL ZDANUL(EPSFF,6)
458 : IF (NCA.NE.0) THEN
459 : C
460 : C Calcul de DGDS0/DGDSF pour SIG0/SIGGF (etats initial/duree calculee)
460 : CALL UO2FL(NCA,NN,XINVL,PENTE,SIG0,GS,DFDS,DGDS0,HDFDQ,KERRE)
460 : CALL UO2FL(NCA,NN,XINVL,PENTE,SIGGF,GS,DFDS,DGDSF,HDFDQ,KERRE)
461 : C
462 : DO IJ=1,NCA
463 : JJ=NN(IJ)
464 : DO I=1,6
465 : DXF(IJ)=W(IJ)*XINVL(IJ)
466 : EPSFF(I)=EPSFF(I)+
467 : & (DGDSF(I,JJ)*DXF(IJ))-(DGDS0(I,JJ)*DX00(IJ))
468 : ENDDO
469 : ENDDO
470 : ENDIF
471 : C
472 : C ----- MAJ INCREMENT DES DEF. INELAST. (FISS.+VISCOPLAST.)
473 : C
474 : DO I=1,6
475 : EPSFF(I)=EPSFF(I)+(EPSVF(I)-EPSV(I))
476 : ENDDO
477 : C
478 : C ------- MAJ VARIABLES INTERNES DE GATT_MONERIE -------
479 : C
480 : DO I=1,NGAT
481 : VARGAF(I)=VARF(I)
482 : ENDDO
483 : C
484 : C --------------------------------------------------
485 : C MAJ CONTR. - INCR. DES DEF. INELAST. - SORTIE
486 : C --------------------------------------------------
487 : DO I=1,6
488 : SIGF(I)=SIGGF(I)
489 : DEFP(I)=EPSFF(I)
490 : ENDDO
491 : C
492 : C
493 : IF(IIMPI.EQ.42) THEN
494 : WRITE(IOIMP,77894) XX
495 : 77894 FORMAT(1X,' UO2CA - DUREE CALCULEE XX '/1X,1PE12.5/)
496 : WRITE(IOIMP,77895) (SIGF(IJ),IJ=1,6)
497 : 77895 FORMAT(1X,' UO2CA - SIGF '/(6(1X,1PE12.5))/)
498 : WRITE(IOIMP,77896) (DEFP(IJ),IJ=1,6)
499 : 77896 FORMAT(1X,' UO2CA - DEFP '/(6(1X,1PE12.5))/)
500 : ENDIF
501 : C
502 : RETURN
503 : END
504 :
505 :
506 :
507 :
508 :
509 :
510 :