Télécharger la plaquette Cast3M


fr en

Télécharger uo2ca.eso

Retour à la liste

Numérotation des lignes : 

uo2ca

  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.  
  75. -INC PPARAM
  76. -INC CCOPTIO
  77. C
  78.       PARAMETER (XZER=0.D0,NC=3,NGATT=4)
  79. C
  80.       DIMENSION SIG0(*),DEPST(*),EPSV(*),VARGA0(*)
  81.       DIMENSION W(*),WMAX(*),SMAX(*),WRUPT(*),BILIN(*),SBILI(*),XLTR(*)
  82.       DIMENSION XINVL(*),WREOUV(*)
  83.       DIMENSION GS(*),XMAT(*)
  84.       DIMENSION VF(3,*),IRESU(*)
  85.       DIMENSION SIGF(*),DEFP(*),EPSVF(*),VARGAF(*)
  86. C
  87.       DIMENSION EPSPT(6),SIGG0(6),EPSV0(6),VAR0(NGATT)
  88.       DIMENSION W0(3),WMAX0(3),SMAX0(3),DX0(NC),DX00(NC),WREOU0(3)
  89.       DIMENSION SIGGF(6),EPSVVF(6),VARF(NGATT),WF(3),DXF(NC),EPSFF(6)
  90.       DIMENSION DGDS0(6,NC),DGDSF(6,NC),DFDS(6,NC),HDFDQ(NC)
  91.       DIMENSION WMAXF(3),SMAXF(3),WREOUF(3),DXV1(NC)
  92. C
  93.       DIMENSION NN(NC),FC(NC),FC2(NC),PENTE(NC),PENTE2(NC),LEBIL(NC)
  94.       DIMENSION FCV(20),FCV1(20),MM(20)
  95. C
  96.       DIMENSION ISING(NC),IFERM(NC),IBRUP(NC),NN0(NC)
  97. C
  98.       DIMENSION JFIS(3)
  99. C
  100.       DATA ITMAX/25/
  101.       DATA IRMAX/5/
  102. C
  103. C--------------------
  104. C   INITIALISATIONS
  105. C--------------------
  106.       KERRE=0
  107.       NC1=NC+1
  108.       IRED=0
  109.       IFLAG1=1
  110.       YOUN=XMAT(1)
  111. C
  112. C --- PAS DE TEMPS CALCULE INITIAL
  113. C
  114.       DTCALC=DT
  115. C
  116. C ------------ LA VITESSE DE LA DEFORMATION TOTALE EST CTE ------------
  117. C
  118.       DO 10 I=1,6
  119.         EPSPT(I)=DEPST(I)/DT
  120. 10    CONTINUE
  121.       IF (IFOUR.EQ.-2) EPSPT(1)=XZER
  122. C
  123. C
  124.       IF(IIMPI.EQ.42) THEN
  125.          WRITE(IOIMP,74433) (EPSPT(I),I=1,6)
  126. 74433    FORMAT(1X,' EPSPT ',6(1X,1PE12.5)/)
  127.       ENDIF
  128. C
  129. C ------------ DEFORMATION DE FISSURATION EN OUVERTURE ------------
  130. C
  131.       DO I=1,3
  132.          DX00(I)=W(I)*XINVL(I)
  133.       ENDDO
  134. C
  135. C
  136. 888   CONTINUE
  137. C
  138.       IRED=IRED+1
  139.       IF(IRED.GT.IRMAX) THEN
  140.          PRINT *,'&&&&&&&& UO2CA  IRED = ',IRED
  141.          KERRE=989
  142.          RETURN
  143.       ENDIF
  144. C
  145.       IF(IRED.GT.1) DTCALC=0.5*DTCALC
  146. C
  147. C
  148. C
  149.       DTLEFT = DTCALC
  150.       TAU    = DTLEFT
  151.       XX     = XZER
  152. C
  153.       T0=TETA1
  154.       FI0=FI1
  155. C
  156. C ----- CONTRAINTES - DEF. VISCOPLAST. -----
  157. C
  158.       DO I=1,6
  159.         SIGG0(I)=SIG0(I)
  160.         EPSV0(I)=EPSV(I)
  161.       ENDDO
  162. C
  163. C ------- VARIABLES INTERNES DE GATT_MONERIE -------
  164. C
  165.       DO I=1,NGAT
  166.         VAR0(I)=VARGA0(I)
  167.       ENDDO
  168. C
  169. C ------------ VARIABLES DE FISSURATION ------------
  170. C
  171.       DO I=1,3
  172.          W0(I)=W(I)
  173.          WMAX0(I)=WMAX(I)
  174.          WREOU0(I) = BTR*MIN(WMAX(I),WRUPT(I))
  175.          SMAX0(I)=SMAX(I)
  176.          DX0(I)=W0(I)*XINVL(I)
  177.       ENDDO
  178. C
  179. C
  180.       NSSINC=0
  181. C
  182. C --------------------------------------------------------------
  183. C       DEBUT DES ITERATIONS EN SSINCREMENTS / FIN SI DTLEFT = 0
  184. C --------------------------------------------------------------
  185. 1400  CONTINUE
  186.       ITERU=0
  187. 1500  CONTINUE
  188. C
  189.       NSSINC = NSSINC + 1
  190.       IF (NSSINC.GT.MSOUPA) THEN
  191.         KERRE=268
  192.         PRINT *,' UO2CA - PAS DE CONVERGENCE EN VISCOPLASTICITE'
  193.         RETURN
  194.       ENDIF
  195. C
  196. C-----------------------------------------------------------------------
  197. C  Etat de fissuration au temps XX et crit. susceptibles d etre atteints
  198. C-----------------------------------------------------------------------
  199.       CALL ZDANUL(DXV1,3)
  200. C
  201.       CALL UO2ET(NC,NN,SIGG0,W0,WMAX0,SMAX0,BILIN,WRUPT,BTR,
  202.      &           XLTR,XINVL,SBILI,PRECIE,PRECIZ,FC,FC2,
  203.      &           LEBIL,NFISSU,NVF,VF,FCV,PENTE,PENTE2,
  204.      &           DXV1,YOUN,NCA,MC,MM,ISING,IFERM,IBRUP,KERRE)
  205.       IF(KERRE.NE.0) RETURN
  206. C
  207.       DO I=1,NC
  208.          ISING(I)=0
  209.          IF(LEBIL(I).EQ.2) ISING(I)=1
  210.       ENDDO
  211. C
  212. C-------------------------------------------------------------
  213. C   En cas de bifurcation, il faut choisir son chemin
  214. C   si la cond. ci-dessous est vraie NCA est forcement non nul
  215. C-------------------------------------------------------------
  216.       IF(ISING(1)+ISING(2)+ISING(3).NE.0) THEN
  217.          CALL UO2IN(ISING,IFERM,IBRUP,LEBIL,PENTE,PENTE2,NCA,NN,MC,
  218.      &              MM,SIGG0,GS,FC,XINVL,XLTR,WMAX0,WRUPT,PRECIZ,T0,
  219.      &              TPOINT,FI0,FPOINT,PRECIS,XMAT,NCOMAT,NSIMP,AAD,BTR,
  220.      &              EPSPT,EPSV0,VAR0,W0,DX0,NGAT,TAU,KERRE)
  221.          IF(KERRE.NE.0) RETURN
  222.       ENDIF
  223. C
  224. C
  225. C
  226.       ITERI=0
  227. 444   CONTINUE
  228. C
  229.       IF (ITERI.NE.0) THEN
  230. C     --------------------
  231.         IF(IIMPI.EQ.42) THEN
  232.           WRITE(IOIMP,74078)
  233. 74078     FORMAT(/2X,' ON REINITIALISE LE PROCESSUS'//)
  234.         ENDIF
  235.         IF(ITERI.GT.NC) THEN
  236.           KERRE=989
  237.           PRINT *,' UO2CA - PAS DE CONVERGENCE - ITERI ',ITERI
  238.           RETURN
  239.         ENDIF
  240. C
  241. C       calcul de PENTE pour tableau NN reduit NN(ITERI+1)<NN(ITERI)
  242. C
  243.         CALL UO2CR(NC,NCA,NN,SIGG0,W0,WMAX0,SMAX0,BILIN,WRUPT,BTR,
  244.      &             XLTR,XINVL,SBILI,FC,PENTE,LEBIL,PRECIE,PRECIZ,KERRE)
  245.         IF(KERRE.NE.0) RETURN
  246. C
  247.       ENDIF
  248. C     -----
  249. C
  250.       NDIM=NCA
  251.       IF(IFOUR.EQ.-2) NDIM=NCA+1
  252. C
  253. C --- determination d un etat converge pour TAU inf. ou egal a DTCALC
  254. C
  255.       CALL UO2DCN(IFLAG1,T0,TPOINT,FI0,FPOINT,PRECIS,PRECIZ,
  256.      &            XMAT,NCOMAT,NSIMP,AAD,BTR,GS,WRUPT,LEBIL,XINVL,PENTE,
  257.      &            EPSPT,SIGG0,EPSV0,VAR0,W0,WMAX0,WREOU0,DX0,
  258.      &            NGAT,NC1,NCA,NDIM,NN,TAU,TAUNEX,SIGGF,EPSVVF,
  259.      &            VARF,WF,DXF,WMAXF,WREOUF,TF,FIF,KERRE)
  260.       IF (KERRE.NE.0) THEN
  261.         RETURN
  262.       ENDIF
  263. C
  264. C --- Viole-t-on un critere de fissuration pour l etat converge ?
  265. C
  266.       CALL ZDANUL(DXV1,3)
  267. C
  268.       CALL UO2CE(MC,MM,SIGG0,DXV1,W0,WMAX0,SMAX0,WRUPT,XLTR,XINVL,BTR,
  269.      &           NFISSU,NVF,FCV1,VF,YOUN,PRECIZ,JFIS,KERRE)
  270.       IF(KERRE.NE.0) RETURN
  271. C
  272.       CALL UO2XX(T0,TPOINT,FI0,FPOINT,PRECIS,PRECIZ,
  273.      &            XMAT,NCOMAT,NSIMP,AAD,GS,WRUPT,LEBIL,XINVL,PENTE,BTR,
  274.      &            YOUN,EPSPT,SIGG0,EPSV0,VAR0,W0,WMAX0,WREOU0,DX0,
  275.      &            SMAX0,XLTR,NFISSU,NVF,VF,MC,MM,FCV1,NGAT,NC1,NCA,NDIM,
  276.      &            NN,TAU,SIGGF,EPSVVF,VARF,WF,DXF,WMAXF,WREOUF,
  277.      &            TF,FIF,MCA,XXX,IRESU,NRESU,LERED,KERRE)
  278.       IF(KERRE.NE.0) RETURN
  279.       IF(LERED.EQ.1) GO TO 888
  280. C
  281.       IF(XXX.LE.0.D0) THEN
  282.          KERRE=989
  283.          PRINT *,' UO2CA - XXX EST NEGATIF OU NUL '
  284.          RETURN
  285.       ENDIF
  286. C
  287. C
  288. C --- Controle de la coherence des resultats pour XXX*TAU
  289. C     ------------------
  290.       IF (NCA.NE.0) THEN
  291. C     ------------------
  292.       NCA0=NCA
  293.       NCAA=0
  294.       DO IJ=1,NCA
  295.          NN0(IJ)=NN(IJ)
  296.          JJ=NN(IJ)
  297.          DXSOM=DXF(JJ)-DX0(JJ)
  298. C
  299.          IF(DXSOM.GE.0.D0) THEN
  300.             NCAA=NCAA+1
  301.             NN(NCAA)=JJ
  302.          ELSE IF(DXSOM.LT.0.D0) THEN
  303. C
  304.             IF(LEBIL(JJ).EQ.1) THEN
  305. C
  306.                IF(IFERM(JJ).NE.1) THEN
  307.                   NCAA=NCAA+1
  308.                   NN(NCAA)=JJ
  309.                ELSE IF(IFERM(JJ).EQ.1) THEN
  310.                   IF(WF(JJ).GT.WREOUF(JJ)) THEN
  311.                      NCAA=NCAA+1
  312.                      NN(NCAA)=JJ
  313.                   ENDIF
  314.                ENDIF
  315. C
  316.             ENDIF
  317.          ENDIF
  318.       ENDDO
  319. C
  320.       IF(IIMPI.EQ.42) THEN
  321.          WRITE(IOIMP,77008) NCA,NCAA,(NN(IJ),IJ=1,NCAA)
  322. 77008    FORMAT(5X,'NCA=',I3,2X,'NCAA=',I3/2X,'NN ' ,8I3/)
  323.       ENDIF
  324. C
  325.       IF(NCAA.EQ.0) THEN
  326.         IF(IIMPI.EQ.42) THEN
  327.            PRINT *,' UO2CA - NCAA EST NUL '
  328.         ENDIF
  329.            NCA=NCAA
  330.            ITERI=ITERI+1
  331.            GO TO 444
  332.       ENDIF
  333. C
  334. C
  335.       IF(NCA.NE.NCAA) THEN
  336.          IF(NCAA.GT.NCA) THEN
  337.             PRINT *,' UO2CA - NCAA EST SUPERIEUR A NCA '
  338.             KERRE=989
  339.             RETURN
  340.          ENDIF
  341.          NCA=NCAA
  342.          ITERI=ITERI+1
  343.          GO TO 444
  344.       ENDIF
  345. C     -----
  346.       ENDIF
  347. C     -----
  348. C
  349. C     -------------------------------------
  350.       IF ((NCA.EQ.0).AND.(ITERI.NE.0)) THEN
  351. C     -------------------------------------
  352.       DO IJ=1,NCA0
  353.          JJ=NN0(IJ)
  354.          IF (IFERM(JJ).EQ.1) THEN
  355.             SXDIF=SIGGF(JJ)-SIGG0(JJ)
  356.             IF (SXDIF.GT.XZER) THEN
  357.                PRINT *,'ATTENTION FOIREUX IB IGAU',IB,IGAU
  358.                ITERU=ITERU+1
  359.                IF(ITERU.GT.IRMAX) THEN
  360.                  PRINT *,'&&&&&&&& UO2CA  ITERU = ',ITERU
  361.                  KERRE=989
  362.                  RETURN
  363.                ENDIF
  364.                TAU=0.5*TAU
  365.                GOTO 1500
  366.             ENDIF
  367.          ENDIF
  368.       ENDDO
  369. C     -----
  370.       ENDIF
  371. C     -----
  372. C --- Fin du controle de la coherence des resultats pour XXX*TAU
  373. C
  374. C
  375.       DO I=1,3
  376.          IF(WF(I).GE.WMAXF(I)) THEN
  377.             WMAXF(I)=WF(I)
  378.             CALL OTTOFU(XINVL,XLTR,WRUPT,BTR,BILIN,
  379.      &                  SBILI,WF,WMAXF,SMAXF,PENTE,I)
  380.          ELSE
  381.             SMAXF(I)=SMAX0(I)
  382.          ENDIF
  383.          WREOUF(I)=BTR*MIN(WMAXF(I),WRUPT(I))
  384.       ENDDO
  385. C
  386. C
  387. C ---------------------
  388. C ENCORE QUELQUES TESTS
  389. C ---------------------
  390.       CALL UO2VE(NRESU,IRESU,SIGGF,WF,WMAXF,WRUPT,SMAXF,
  391.      &           BTR,XLTR,XINVL,BILIN,SBILI,PRECIE,PRECIZ,
  392.      &           NN,NN0,NC,NCA,FC,FC2,PENTE,PENTE2,LEBIL,ISING,KERRE)
  393.       IF(KERRE.NE.0) RETURN
  394. C
  395. C
  396.       IF (MCA.NE.0) THEN
  397.         XX=XX+(XXX*TAU)
  398.         GOTO 777
  399.       ELSE
  400.         XX=XX+TAU
  401.       ENDIF
  402. C
  403. C ----- MAJ CONTRAINTES - DEF. VISCOPLAST. A t0+XX -----
  404. C
  405.       DO I=1,6
  406.         SIGG0(I)=SIGGF(I)
  407.         EPSV0(I)=EPSVVF(I)
  408.       ENDDO
  409. C
  410. C ------- MAJ VARIABLES INTERNES DE GATT_MONERIE A t0+XX -------
  411. C
  412.       DO I=1,NGAT
  413.         VAR0(I)=VARF(I)
  414.       ENDDO
  415. C
  416. C ------------ MAJ VARIABLES DE FISSURATION A t0+XX ------------
  417. C
  418.       DO I=1,3
  419.          W0(I)=WF(I)
  420.          DX0(I)=DXF(I)
  421.          WMAX0(I)=WMAXF(I)
  422.          SMAX0(I)=SMAXF(I)
  423.       ENDDO
  424. C
  425.       T0=TF
  426.       FI0=FIF
  427. C
  428. C ----------------------------------------------------
  429. C  TEST DE FIN SS INCREMENTS /  MISE A JOUR DE TAU
  430. C    fin des boucles en ss increments  si tau = dtleft
  431. C ----------------------------------------------------
  432. C
  433.       IF (TAU.LT.DTLEFT) THEN
  434.         DTLEFT = DTLEFT - TAU
  435.         IF (TAUNEX.GT.DTLEFT) TAUNEX = DTLEFT
  436.         TAU=TAUNEX
  437.         GOTO 1400
  438.       ENDIF
  439. C
  440.       IF (ABS(TAU-DTLEFT).GT.(TAU/1000.)) THEN
  441.          PRINT *, ' UO2CA - PROBLEME  TAU > DTLEFT '
  442.          KERRE=989
  443.          RETURN
  444.       ENDIF
  445. C
  446.       XX=DTCALC
  447. C
  448.  777  CONTINUE
  449. C
  450. C -------------------------
  451. C MAJ DES VARIABLES A t0+XX
  452. C -------------------------
  453. C ------------ MAJ VARIABLES DE FISSURATION ------------
  454. C
  455.       DO I=1,3
  456.          W(I)=WF(I)
  457.          WMAX(I)=WMAXF(I)
  458.          SMAX(I)=SMAXF(I)
  459.          WREOUV(I)=WREOUF(I)
  460.       ENDDO
  461. C
  462. C ----- MAJ DEF. VISCOPLAST. -----
  463. C
  464.       DO I=1,6
  465.         EPSVF(I)=EPSVVF(I)
  466.       ENDDO
  467. C
  468. C ----- MAJ INCREMENT DES DEF. DE FISS. (OUV. + GLISS.)
  469. C
  470.       CALL ZDANUL(EPSFF,6)
  471.       IF (NCA.NE.0) THEN
  472. C
  473. C      Calcul de DGDS0/DGDSF pour SIG0/SIGGF (etats initial/duree calculee)
  474.         CALL UO2FL(NCA,NN,XINVL,PENTE,SIG0,GS,DFDS,DGDS0,HDFDQ,KERRE)
  475.         CALL UO2FL(NCA,NN,XINVL,PENTE,SIGGF,GS,DFDS,DGDSF,HDFDQ,KERRE)
  476. C
  477.         DO IJ=1,NCA
  478.            JJ=NN(IJ)
  479.            DO I=1,6
  480.              DXF(IJ)=W(IJ)*XINVL(IJ)
  481.              EPSFF(I)=EPSFF(I)+
  482.      &                     (DGDSF(I,JJ)*DXF(IJ))-(DGDS0(I,JJ)*DX00(IJ))
  483.            ENDDO
  484.         ENDDO
  485.       ENDIF
  486. C
  487. C ----- MAJ INCREMENT DES DEF. INELAST. (FISS.+VISCOPLAST.)
  488. C
  489.       DO I=1,6
  490.         EPSFF(I)=EPSFF(I)+(EPSVF(I)-EPSV(I))
  491.       ENDDO
  492. C
  493. C ------- MAJ VARIABLES INTERNES DE GATT_MONERIE -------
  494. C
  495.       DO I=1,NGAT
  496.         VARGAF(I)=VARF(I)
  497.       ENDDO
  498. C
  499. C --------------------------------------------------
  500. C MAJ CONTR. - INCR. DES DEF. INELAST. - SORTIE
  501. C --------------------------------------------------
  502.       DO I=1,6
  503.         SIGF(I)=SIGGF(I)
  504.         DEFP(I)=EPSFF(I)
  505.       ENDDO
  506. C
  507. C
  508.       IF(IIMPI.EQ.42) THEN
  509.          WRITE(IOIMP,77894) XX
  510. 77894    FORMAT(1X,'  UO2CA - DUREE CALCULEE XX '/1X,1PE12.5/)
  511.          WRITE(IOIMP,77895) (SIGF(IJ),IJ=1,6)
  512. 77895    FORMAT(1X,'  UO2CA - SIGF    '/(6(1X,1PE12.5))/)
  513.          WRITE(IOIMP,77896) (DEFP(IJ),IJ=1,6)
  514. 77896    FORMAT(1X,'  UO2CA - DEFP '/(6(1X,1PE12.5))/)
  515.       ENDIF
  516. C
  517.       RETURN
  518.       END
  519.  
  520.  
  521.  
  522.  
  523.  
  524.  
  525.  
  526.  

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