Télécharger la plaquette Cast3M


fr en

Télécharger concyc2.eso

Retour à la liste

Numérotation des lignes : 

concyc2

  1. C CONCYC2   SOURCE    CB215821  16/04/21    21:15:55     8920
  2.       SUBROUTINE CONCYC2(XMAT,SIG0,SIGF,VAR0,VARF,DEPST,ISTEP)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de mettre en oeuvre le
  6. C                modele CONCYC (M. Vassaux & B. Richard) pour representer
  7. C                le comportement 2D/3D d'un beton ordinaire sous
  8. C                sollicitations cycliques
  9. C
  10. C Traits       : - Endommagement scalaire
  11. C                - Boucles hysteretiques
  12. C                - Deformations permanentes
  13. C                - Refermeture lineaire des fissures
  14. C                - Effet unilatteral complet
  15. C                - Non localite portant sur le seuil d'endommagement (Fd)
  16. C
  17. C Auteur       : M. Vassaux
  18. C
  19. C Co-auteur    : B. Richard
  20. C
  21. C Date         : 2010 - 2011
  22. C====&===1=========2=========3=========4=========5=========6=========7==
  23. C
  24. C----DECLARATION GENERALES----------------------------------------------
  25. C
  26.       IMPLICIT REAL*8(A-H,O-Z)
  27. C
  28. C-----DECLARATION GENERALE----------------------------------------------
  29. C
  30.       REAL*8 UNIT(3,3)
  31.       REAL*8 SIGTOT(3,3),SIGMAT(3,3),SIGFIS(3,3),QF(3,3),SIGFISEFF(3,3)
  32.       REAL*8 EPSTOT(3,3),EPSTOTPREC(3,3)
  33.       REAL*8 EPSFIS(3,3),DEPSFIS(3,3),EPSFISMAX(3,3),EPSFISP(3,3)
  34.       REAL*8 J11(3,3),J22(3,3),J12(3,3),J21(3,3),DSFIEFIJ12(3,3)
  35.       REAL*8 SIGFISEFFDEV(3,3),J12DEV(3,3)
  36.       REAL*8 EPSTOTDEV(3,3),EPSFISDEV(3,3),DEPSFISDEV(3,3)
  37.       REAL*8 EPSTOTPOS(3,3),EPSTOTNEG(3,3)
  38.       REAL*8 DI2SIGFISEFFDEV(3,3),EPSIPP(3),VALP33(3,3)
  39.       REAL*8 SIGELA(3,3), SIGELADEV(3,3),SIGIPP(3),VALS33(3,3)
  40.       REAL*8 XMAT(*),SIG0(*),SIGF(*),VAR0(*),VARF(*),DEPST(*)
  41. C
  42. C-----DEFINITION DE LA MATRICE UNITE D ORDRE 3------------------------------
  43. C
  44.       DO I=1,3
  45.          DO J=1,3
  46.             IF (I.EQ.J) THEN
  47.                UNIT(I,J) = 1.0D0
  48.             ELSE
  49.                UNIT(I,J) = 0.0D0
  50.             ENDIF
  51.          ENDDO
  52.       ENDDO
  53. C
  54. C-----DEFINITION DE PARAMETRES------------------------------------------
  55. C
  56.       IMAX = 100
  57. C
  58. C-----CHARGEMENT DES PARAMETRES MATERIAUX----------------------------------------------------
  59. C
  60.       XE    = XMAT(1)
  61.       XNU   = XMAT(2)
  62.       XK    = XE/((1.0D0-2.D0*XNU)*3.0D0)
  63.       XG    = XE/(2.0D0*(1.D0+XNU))
  64.  
  65.       XEND  = XMAT(5)
  66.  
  67.       XSIGT = XMAT(6)
  68.       ATRA  = XMAT(7)
  69.       BTRA  = XMAT(8)
  70.  
  71.       XQP   = XMAT(9)
  72.       XCF   = XMAT(10)
  73.  
  74. C     Ecrouissage cinematique
  75.       XHF   = 0.0
  76. C     Non-linearite de l'ecrouissage cinematique
  77.       XAF   = 0.0
  78.  
  79. C
  80. C-----CHARGEMENT DES VARIABLES INTERNES----------------------------------------------------
  81. C
  82. C     Deformations totales
  83.       EPSTOT(1,1) = VAR0(1)
  84.       EPSTOT(2,2) = VAR0(2)
  85.       EPSTOT(3,3) = VAR0(3)
  86.       EPSTOT(1,2) = VAR0(4)
  87.       EPSTOT(1,3) = VAR0(5)
  88.       EPSTOT(2,3) = VAR0(6)
  89.       EPSTOT(2,1) = VAR0(4)
  90.       EPSTOT(3,1) = VAR0(5)
  91.       EPSTOT(3,2) = VAR0(6)
  92.  
  93. C     Contrainte initiales
  94.       SIGTOT(1,1) = SIG0(1)
  95.       SIGTOT(2,2) = SIG0(2)
  96.       SIGTOT(3,3) = SIG0(3)
  97.       SIGTOT(1,2) = SIG0(4)
  98.       SIGTOT(1,3) = SIG0(5)
  99.       SIGTOT(2,3) = SIG0(6)
  100.       SIGTOT(2,1) = SIG0(4)
  101.       SIGTOT(3,1) = SIG0(5)
  102.       SIGTOT(3,2) = SIG0(6)
  103.  
  104. C     Contrainte fissures initiales
  105.       SIGFIS(1,1) = VAR0(7)
  106.       SIGFIS(2,2) = VAR0(8)
  107.       SIGFIS(3,3) = VAR0(9)
  108.       SIGFIS(1,2) = VAR0(10)
  109.       SIGFIS(1,3) = VAR0(11)
  110.       SIGFIS(2,3) = VAR0(12)
  111.       SIGFIS(2,1) = VAR0(10)
  112.       SIGFIS(3,1) = VAR0(11)
  113.       SIGFIS(3,2) = VAR0(12)
  114.  
  115. C     Tenseur d'ecrouissage cinematique
  116.       QF(1,1) = VAR0(13)
  117.       QF(2,2) = VAR0(14)
  118.       QF(3,3) = VAR0(15)
  119.       QF(1,2) = VAR0(16)
  120.       QF(1,3) = VAR0(17)
  121.       QF(2,3) = VAR0(18)
  122.       QF(2,1) = VAR0(16)
  123.       QF(3,1) = VAR0(17)
  124.       QF(3,2) = VAR0(18)
  125.  
  126. C     Endommagement
  127.       XDT = VAR0(19)
  128.       XDT0= VAR0(19)
  129.       XZT = VAR0(20)
  130.  
  131. C     Fonctions de refermeture associees a la trace et au deviateur
  132.       XPHI = VAR0(21)
  133.  
  134. C     Deformations associées a la fissuration maximale enregistrees
  135.       EPSFISMAX(1,1) = VAR0(22)
  136.       EPSFISMAX(2,2) = VAR0(23)
  137.       EPSFISMAX(3,3) = VAR0(24)
  138.       EPSFISMAX(1,2) = VAR0(25)
  139.       EPSFISMAX(1,3) = VAR0(26)
  140.       EPSFISMAX(2,3) = VAR0(27)
  141.       EPSFISMAX(2,1) = VAR0(25)
  142.       EPSFISMAX(3,1) = VAR0(26)
  143.       EPSFISMAX(3,2) = VAR0(27)
  144.  
  145. C     Deformations plastiques associées au glissement
  146.       EPSFISP(1,1) = VAR0(28)
  147.       EPSFISP(2,2) = VAR0(29)
  148.       EPSFISP(3,3) = VAR0(30)
  149.       EPSFISP(1,2) = VAR0(31)
  150.       EPSFISP(1,3) = VAR0(32)
  151.       EPSFISP(2,3) = VAR0(33)
  152.       EPSFISP(2,1) = VAR0(31)
  153.       EPSFISP(3,1) = VAR0(32)
  154.       EPSFISP(3,2) = VAR0(33)
  155.  
  156. C     Stockage des variables au pas precedent
  157.       XDTPREC = XDT
  158.       DO I=1,3
  159.          DO J=1,3
  160.             EPSTOTPREC(I,J) = EPSTOT(I,J)
  161.          ENDDO
  162.       ENDDO
  163. C     Fonctions de refermeture associees a la trace et au deviateur
  164.       XEDISS = VAR0(34)
  165. C
  166. C-----ACTUALISATION DES DEFORMATIONS------------------------------------
  167. C
  168.       EPSTOT(1,1) = EPSTOT(1,1) + DEPST(1)
  169.       EPSTOT(2,2) = EPSTOT(2,2) + DEPST(2)
  170.       EPSTOT(3,3) = EPSTOT(3,3) + DEPST(3)
  171.       EPSTOT(1,2) = EPSTOT(1,2) + (DEPST(4)/2.0D0)
  172.       EPSTOT(1,3) = EPSTOT(1,3) + (DEPST(5)/2.0D0)
  173.       EPSTOT(2,3) = EPSTOT(2,3) + (DEPST(6)/2.0D0)
  174.       EPSTOT(2,1) = EPSTOT(2,1) + (DEPST(4)/2.0D0)
  175.       EPSTOT(3,1) = EPSTOT(3,1) + (DEPST(5)/2.0D0)
  176.       EPSTOT(3,2) = EPSTOT(3,2) + (DEPST(6)/2.0D0)
  177.  
  178.       CALL JACOB3(EPSTOT,3,EPSIPP,VALP33)
  179.  
  180.       EPSTIL=MAX(0.0D0,EPSIPP(1))**2.0D0 +
  181.      1       MAX(0.0D0,EPSIPP(2))**2.0D0 +
  182.      2       MAX(0.0D0,EPSIPP(3))**2.0D0
  183.  
  184.       EPSTIL=SQRT(EPSTIL)
  185.  
  186.       XEPSTOTTRAC = EPSTOT(1,1)+EPSTOT(2,2)+EPSTOT(3,3)
  187.  
  188.       CALL DEVIAT(EPSTOT,EPSTOTDEV)
  189.  
  190.       DO I=1,3
  191.          DO J=1,3
  192.             SIGELA(I,J) = (XK*XEPSTOTTRAC*UNIT(I,J)
  193.      &                          + 2.0D0*XG*EPSTOTDEV(I,J))
  194.          ENDDO
  195.       ENDDO
  196.       CALL JACOB3(SIGELA,3,SIGIPP,VALP33)
  197.  
  198. C-----Modele endommagement RICRAG classique ou avec consolidation modifiee
  199.       IF (XEND.EQ.1.OR.XEND.EQ.2) THEN
  200.  
  201.          YP = 0.5D0*XSIGT*EPSTIL
  202.  
  203. C-----Modele endommagement RICRAG avec critere modifie
  204.       ELSEIF (XEND.EQ.3) THEN
  205.  
  206.          XSI1 = MIN(0.0D0,SIGIPP(1)) + MIN(0.0D0,SIGIPP(2))
  207.      &                 + MIN(0.0D0,SIGIPP(3))
  208.  
  209.          EPSTIL=EPSTIL + BTRA*MIN(0.0D0,XSI1)/(3.0*XK)
  210.  
  211.          YP = 0.5D0*XSIGT*EPSTIL
  212.       ENDIF
  213. C
  214. C-----NON LOCAL----------------------------------------------------------
  215. C
  216.       IF (ISTEP.EQ.0) THEN
  217.          SEUIL1=YP
  218.          VARF(35)=YP
  219.       ELSE IF (ISTEP.EQ.1) THEN
  220.          VARF(35)=YP
  221.          DO I=1,34
  222.             VARF(I)=VAR0(I)
  223.          ENDDO
  224.          GOTO 2000
  225.       ELSE IF (ISTEP.EQ.2) THEN
  226.          SEUIL1=VAR0(35)
  227.          VARF(35)=SEUIL1
  228.       ENDIF
  229. C
  230. C-----EVOLUTION DE L ENDOMMAGEMENT---------------------------------------
  231. C
  232. C     Modele endommagement RICRAG classique ou avec critere modifie
  233.       IF (XEND.EQ.1.OR.XEND.EQ.3) THEN
  234.  
  235.          XFD = SEUIL1 - ((XSIGT**2)/(2.0D0*XE) + XZT)
  236.  
  237.          XYY0 = (XSIGT**2)/(2.0D0*XE*SEUIL1)
  238.  
  239.          IF (XFD.GT.0.0D0) THEN
  240.             XDT  = 1.0D0
  241.      1        - XYY0*EXP(-ATRA*(SEUIL1-((XSIGT**2)/(2.0D0*XE))))
  242.             XZT = SEUIL1 - ((XSIGT**2)/(2.0D0*XE))
  243.          ENDIF
  244.  
  245. C     Modele endommagement RICRAG avec consolidation modifiee
  246.       ELSEIF (XEND.EQ.2) THEN
  247.          XSIGEELA = SIGIPP(1)**2.0D0
  248.      1               + SIGIPP(2)**2.0D0
  249.      2                  + SIGIPP(3)**2.0D0
  250.          XSIGEELANEG = MIN(0.0D0,SIGIPP(1))**2.0D0
  251.      1                  + MIN(0.0D0,SIGIPP(2))**2.0D0
  252.      2                     + MIN(0.0D0,SIGIPP(3))**2.0D0
  253.  
  254.          XCONF  = (1 + BTRA*(XSIGEELANEG/XSIGEELA)**0.5)
  255.  
  256.          XYY0 = (XSIGT**2)/(2.0D0*XE*SEUIL1)
  257.  
  258.          XFD = SEUIL1 - ((XSIGT**2)/(2.0D0*XE) + XZT)
  259.  
  260.          IF (XFD.GT.0.0D0) THEN
  261.             XDT  = 1.0D0
  262.      &      - XYY0*EXP(-(ATRA/XCONF)*(SEUIL1-((XSIGT**2)/(2.0D0*XE))))
  263.             XZT = SEUIL1 - ((XSIGT**2)/(2.0D0*XE))
  264.          ENDIF
  265.  
  266.       ENDIF
  267. C
  268. C-----ACTUALISATION DE LA CONTRAINTE DANS LA MATRICE FISSUREE------------
  269. C
  270. C     Calcul de la trace de EPSTOT
  271.       XEPSTOTTRAC = EPSTOT(1,1)+EPSTOT(2,2)+EPSTOT(3,3)
  272.  
  273. C     Calcul du deviateur de EPSTOT
  274.       CALL DEVIAT(EPSTOT,EPSTOTDEV)
  275.  
  276.       DO I=1,3
  277.          DO J=1,3
  278.             SIGMAT(I,J) = (1.0D0-XDT)*(XK*XEPSTOTTRAC*UNIT(I,J)
  279.      &                          + 2.0D0*XG*EPSTOTDEV(I,J))
  280.          ENDDO
  281.       ENDDO
  282.       XI1SIGMAT = (SIGMAT(1,1)+SIGMAT(2,2)+SIGMAT(3,3))/3.0D0
  283. C
  284. C-----ACTUALISATION DE LA DEFORMATION ASSOCIEE AUX FISSURES-------------
  285. C
  286.       DO I=1,3
  287.          DO J=1,3
  288.             EPSFIS(I,J) = EPSTOT(I,J)*XDT
  289.          ENDDO
  290.       ENDDO
  291.       DO I=1,3
  292.          DO J=1,3
  293.             DEPSFIS(I,J) = EPSTOT(I,J)*XDT
  294.      &          - EPSTOTPREC(I,J)*XDTPREC
  295.          ENDDO
  296.       ENDDO
  297.  
  298. C     Calcul de la trace de EPSFIS
  299.       XEPSFISTRAC = EPSFIS(1,1) + EPSFIS(2,2) + EPSFIS(3,3)
  300.       IF (ABS(XEPSFISTRAC).LE.3.0D-12) THEN
  301.           XEPSFISTRAC = 3.0D-12
  302.       ENDIF
  303. C     XEPSFISTRAC = EPSFIS(3,3)
  304. C     Calcul du deviateur de DEPSFIS
  305.       CALL DEVIAT(EPSFIS,EPSFISDEV)
  306.  
  307. C     Calcul de la trace de EPSFIS
  308.       XEPSFISMAXTRAC = EPSFISMAX(1,1) + EPSFISMAX(2,2) + EPSFISMAX(3,3)
  309.  
  310. C     Securite...
  311.       IF (XEPSFISMAXTRAC.LE.0.0D0) THEN
  312.           XEPSFISMAXTRAC = 3.0D-12
  313.       ENDIF
  314.  
  315. C     Mise a jour de EPSFISMAX
  316.       IF (XEPSFISTRAC.GT.XEPSFISMAXTRAC) THEN
  317.           DO I=1,3
  318.               DO J=1,3
  319.                   EPSFISMAX(I,J) = EPSFIS(I,J)
  320.               ENDDO
  321.           ENDDO
  322.       ENDIF
  323. C
  324. C-----CALCUL DES FONCTIONS DE REFERMETURE DE FISSURE--------------------
  325. C
  326. C     Calcul de la trace de EPSFISMAX
  327.       XEPSFISMAXTRAC = EPSFISMAX(1,1) + EPSFISMAX(2,2) + EPSFISMAX(3,3)
  328.  
  329. C     Securite...
  330.       IF (XEPSFISMAXTRAC.LE.0.0D0) THEN
  331.           XEPSFISMAXTRAC = 3.0D-12
  332.       ENDIF
  333.  
  334. C     Calcul de la trace de EPSFISP
  335.       XEPSFISPTRAC = EPSFISP(1,1) + EPSFISP(2,2) + EPSFISP(3,3)
  336.  
  337. C     Calcul de la fonction de refermeture
  338.       XPHI = 1.0D0 - 1.0D0/(1.0D0 + EXP(-1.0D0*XQP*
  339.      &    XEPSFISTRAC/XEPSFISMAXTRAC))
  340. C
  341. C-----CALCUL DE LA CONTRAINTE DANS LES FISSURES: TEST--------------------
  342. C
  343. C     Calcul de la trace de EPSFIS
  344.       XEPSFISTRAC = EPSFIS(1,1) + EPSFIS(2,2) + EPSFIS(3,3)
  345.  
  346. C     Calcul de la trace de DEPSFIS
  347.       XDEPSFISTRAC = DEPSFIS(1,1)+DEPSFIS(2,2)+DEPSFIS(3,3)
  348.  
  349. C     Calcul du deviateur de DEPSFIS
  350.       CALL DEVIAT(DEPSFIS,DEPSFISDEV)
  351.  
  352.       DO I=1,3
  353.          DO J=1,3
  354.             SIGFIS(I,J) = SIGFIS(I,J)
  355.      &          + XPHI*(XK*XDEPSFISTRAC)*UNIT(I,J)
  356.      &          + XPHI*(2.0D0*XG*DEPSFISDEV(I,J))
  357.          ENDDO
  358.       ENDDO
  359. C
  360. C-----SEUIL GLISSEMENT INTERNE FISSURES--------------------------------
  361. C
  362. C     Calcul du I1 de SIGFIS
  363.       XI1SIGFIS = (SIGFIS(1,1)+SIGFIS(2,2)+SIGFIS(3,3))/3.0D0
  364.       XI1SIGTOT = XI1SIGMAT + XI1SIGFIS
  365.  
  366. C     Calcul de la quantite (SIGFIS - QF)
  367.       DO I=1,3
  368.          DO J=1,3
  369.             SIGFISEFF(I,J) = SIGFIS(I,J) - QF(I,J)
  370.          ENDDO
  371.       ENDDO
  372.  
  373. C     Calcul de I2(DEV(SIGFIS - QF))
  374.       CALL DEVIAT(SIGFISEFF,SIGFISEFFDEV)
  375.       CALL I2(SIGFISEFFDEV,XI2SIGFISEFFDEV)
  376.  
  377. C     Calcul de la non-linearite de la dependance au confinement
  378.       XFPI = XI2SIGFISEFFDEV + XCF*(XI1SIGFIS)
  379.  
  380.       XTEDISS = 0.0D0
  381. C     Test sur la positivite du seuil de glissement interne
  382.       IF (XFPI.GT.0.0D0) THEN
  383.           XFPI0 = XFPI
  384.           XLAM = 1.0D0
  385. C         Debut du "Return-Mapping"
  386.           DO NI=1,IMAX
  387. C             On realise un test de sortie
  388.               IF ((ABS(XFPI/XFPI0).LE.(1.0D-8)).OR.
  389.      &            (ABS(XLAM).LE.(1.0D-10))) THEN
  390.                  GOTO 400
  391. C             Si le test n est pas concluant...
  392.               ELSE
  393. C                Calcul de dI2(DEV(SIGFIS - QF))
  394.                  IF (XI2SIGFISEFFDEV.EQ.0.0D0) THEN
  395.                     DO I=1,3
  396.                        DO J=1,3
  397.                           DI2SIGFISEFFDEV(I,J) = 0.0D0
  398.                        ENDDO
  399.                     ENDDO
  400.                  ELSE
  401.                     DO I=1,3
  402.                        DO J=1,3
  403.                           DI2SIGFISEFFDEV(I,J) =
  404.      &                         1.5D0*SIGFISEFFDEV(I,J)/XI2SIGFISEFFDEV
  405.                        ENDDO
  406.                     ENDDO
  407.                  ENDIF
  408.  
  409. C                Calcul des jacobiens
  410.                  DO I=1,3
  411.                     DO J=1,3
  412.                        J11(I,J) = DI2SIGFISEFFDEV(I,J)
  413.      &                             + XCF*(UNIT(I,J)/3.0D0)
  414.                        J12(I,J) = DI2SIGFISEFFDEV(I,J)
  415.                        J21(I,J) = -1.0D0*DI2SIGFISEFFDEV(I,J)
  416.                        J22(I,J) = -1.0D0*DI2SIGFISEFFDEV(I,J)
  417.      &                             + XAF*QF(I,J)
  418.                     ENDDO
  419.                  ENDDO
  420.  
  421. C                Calcul de la trace de DEPSFIS
  422.                  XJ12TRAC = J12(1,1)+J12(2,2)+J12(3,3)
  423.  
  424. C                Calcul du deviateur de DEPSFIS
  425.                  CALL DEVIAT(J12,J12DEV)
  426.  
  427. C                Calcul de (dsigfis/depsfis):J12
  428.                  DO I=1,3
  429.                     DO J=1,3
  430.                         DSFIEFIJ12(I,J)=
  431.      &                          +1.0D0*((XPHI)*XJ12TRAC*UNIT(I,J)
  432.      &                               +  2.0D0*XG*(XPHI)*J12DEV(I,J))
  433.                     ENDDO
  434.                  ENDDO
  435.  
  436. C                Calcul de produits contractes de matrices
  437.                  XJ11JXY = 0.0D0
  438.                  XJ21J22 = 0.0D0
  439.                  DO I=1,3
  440.                     DO J=1,3
  441.                        XJ11JXY = XJ11JXY + J11(I,J)*DSFIEFIJ12(I,J)
  442.                        XJ21J22 = XJ21J22 + J21(I,J)*J22(I,J)
  443.                     ENDDO
  444.                  ENDDO
  445.  
  446.                  XDFPI = (XJ11JXY + XHF*XJ21J22)
  447.  
  448. C                Calcul du multiplicateur plastique
  449.                  IF(XDFPI.NE.0.0D0) THEN
  450.                         XLAM = XFPI/XDFPI
  451.                  ELSE
  452.                         GOTO 400
  453.                  ENDIF
  454.  
  455. C                Actualisation de SIGFIS et QF
  456.                  DO I=1,3
  457.                     DO J=1,3
  458.                        SIGFIS(I,J)= SIGFIS(I,J) - XLAM*DSFIEFIJ12(I,J)
  459.                        QF(I,J) = QF(I,J) - XLAM*(XHF*J22(I,J))
  460.                        EPSFISP(I,J) = EPSFISP(I,J) + XLAM*J12(I,J)
  461.                     ENDDO
  462.                  ENDDO
  463.  
  464. C                Calcul du I1 de SIGFIS
  465.                  XI1SIGFIS =
  466.      &                  (SIGFIS(1,1)+SIGFIS(2,2)+SIGFIS(3,3))/3.0D0
  467.                  XI1SIGTOT = XI1SIGMAT + XI1SIGFIS
  468.  
  469. C                Calcul de la quantite (SIGFIS - QF)
  470.                  DO I=1,3
  471.                     DO J=1,3
  472.                        SIGFISEFF(I,J) = SIGFIS(I,J) - QF(I,J)
  473.                     ENDDO
  474.                  ENDDO
  475.  
  476. C                Calcul de I2(DEV(SIGFIS - QF))
  477.                  CALL DEVIAT(SIGFISEFF,SIGFISEFFDEV)
  478.                  CALL I2(SIGFISEFFDEV,XI2SIGFISEFFDEV)
  479.  
  480. C                Calcul du critere de glissement
  481.                  XFPI = XI2SIGFISEFFDEV + XCF*(XI1SIGFIS)
  482.  
  483. C                Calcul de l'increment d'energie dissipee
  484.                  XTEDISS = 0.0D0
  485.                  DO I=1,3
  486.                    DO J=1,3
  487.                      XTEDISS = XTEDISS
  488.      &                    + ABS(SIGFISEFFDEV(I,J)*XLAM*J12(I,J))
  489.                    ENDDO
  490.                  ENDDO
  491.               ENDIF
  492.  
  493.           ENDDO
  494.  
  495.  400      CONTINUE
  496.       ENDIF
  497. C
  498. C-----CALCUL DE LA CONTRAINTE TOTALE
  499. C
  500.       DO I=1,3
  501.          DO J=1,3
  502.            SIGTOT(I,J) = SIGMAT(I,J) + SIGFIS(I,J)
  503.          ENDDO
  504.       ENDDO
  505.  
  506. C     PRINT *,'----------------FIN   ITERATION------------------'
  507. C
  508. C-----CALCUL DE L ENERGIE DISSIPEE
  509. C
  510.       XEDISS = XEDISS + XTEDISS
  511. C
  512. C-----STOCKAGE DES DEFORMATIONS TOTALES
  513. C
  514.       VARF(1)=EPSTOT(1,1)
  515.       VARF(2)=EPSTOT(2,2)
  516.       VARF(3)=EPSTOT(3,3)
  517.       VARF(4)=EPSTOT(1,2)
  518.       VARF(5)=EPSTOT(1,3)
  519.       VARF(6)=EPSTOT(2,3)
  520. C
  521. C-----STOCKAGE DES CONTRAINTES AVEC SYMMETRISATION POUR ROBUSTESSE
  522. C
  523.       SIGF(1)= SIGTOT(1,1)
  524.       SIGF(2)= SIGTOT(2,2)
  525.       SIGF(3)= SIGTOT(3,3)
  526.       SIGF(4)= (SIGTOT(1,2)+SIGTOT(2,1))/2.0D0
  527.       SIGF(5)= (SIGTOT(1,3)+SIGTOT(3,1))/2.0D0
  528.       SIGF(6)= (SIGTOT(2,3)+SIGTOT(3,2))/2.0D0
  529.  
  530. C
  531. C-----STOCKAGE DES CONTRAINTES DANS LES FISSURES
  532. C
  533.       VARF(7)=SIGFIS(1,1)
  534.       VARF(8)=SIGFIS(2,2)
  535.       VARF(9)=SIGFIS(3,3)
  536.       VARF(10)=(SIGFIS(1,2)+SIGFIS(2,1))/2.0D0
  537.       VARF(11)=(SIGFIS(1,3)+SIGFIS(3,1))/2.0D0
  538.       VARF(12)=(SIGFIS(2,3)+SIGFIS(3,2))/2.0D0
  539. C
  540. C-----STOCKAGE DE L ECROUISSAGE CINEMATIQUE DE GLISSEMENT
  541. C
  542.       VARF(13) =QF(1,1)
  543.       VARF(14) =QF(2,2)
  544.       VARF(15) =QF(3,3)
  545.       VARF(16)=(QF(1,2)+QF(2,1))/2.0D0
  546.       VARF(17)=(QF(1,3)+QF(3,1))/2.0D0
  547.       VARF(18)=(QF(2,3)+QF(3,2))/2.0D0
  548. C
  549. C-----STOCKAGE DE L ENDOMMAGEMENT
  550. C
  551.       VARF(19) = XDT
  552.       VARF(20) = XZT
  553. C
  554. C-----STOCKAGE DES FONCTIONS DE REFERMETURE
  555. C
  556.       VARF(21) = XPHI
  557. C
  558. C-----STOCKAGE DEFORMATIONS MAXIMALES ASSOCIEES AU FISSURES
  559. C
  560.       VARF(22) = EPSFISMAX(1,1)
  561.       VARF(23) = EPSFISMAX(2,2)
  562.       VARF(24) = EPSFISMAX(3,3)
  563.       VARF(25) = (EPSFISMAX(1,2)+EPSFISMAX(2,1))/2.0D0
  564.       VARF(26) = (EPSFISMAX(1,3)+EPSFISMAX(3,1))/2.0D0
  565.       VARF(27) = (EPSFISMAX(2,3)+EPSFISMAX(3,2))/2.0D0
  566. C-----STOCKAGE DEFORMATIONS PLASTIQUES ASSOCIEES AU GLISSEMENT
  567. C
  568.       VARF(28) = EPSFISP(1,1)
  569.       VARF(29) = EPSFISP(2,2)
  570.       VARF(30) = EPSFISP(3,3)
  571.       VARF(31) = (EPSFISP(1,2)+EPSFISP(2,1))/2.0D0
  572.       VARF(32) = (EPSFISP(1,3)+EPSFISP(3,1))/2.0D0
  573.       VARF(33) = (EPSFISP(2,3)+EPSFISP(3,2))/2.0D0
  574. C
  575. C-----STOCKAGE DE L ENERGIE DISSIPEE PAR FRICTION
  576. C
  577.       VARF(34) = XEDISS
  578. C
  579. C-----BALISE DE SORTIE POUR LE NON LOCAL
  580. C
  581.  2000 CONTINUE
  582.  
  583.       RETURN
  584.       END
  585.  
  586.  
  587.  
  588.  
  589.  
  590.  
  591.  
  592.  
  593.  
  594.  
  595.  
  596.  
  597.  
  598.  
  599.  
  600.  
  601.  
  602.  
  603.  
  604.  
  605.  
  606.  
  607.  
  608.  
  609.  

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