Télécharger la plaquette Cast3M


fr en

Télécharger ricbet3d.eso

Retour à la liste

Numérotation des lignes : 

ricbet3d

  1. C RICBET3D  SOURCE    CB215821  16/04/21    21:18:15     8920
  2.           SUBROUTINE RICBET3D(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 RICBET (B. Richard) pour representer
  7. C                le comportement 3D d'un béton ordinaire sous
  8. C                sollicitations cycliques
  9. C
  10. C Traits       : - Endommagement scalaire
  11. C                - Glissement frottant
  12. C                - 3D
  13. C                - Monotone et cyclique
  14. C                - Effet unilatteral
  15. C                - Comportement lineaire en compression
  16. C                - Non localite portant sur le seuil d'endommagement (Fd)
  17. C
  18. C Auteurs      : B. Richard (Dr - Ing.) & F. Ragueneau (Pr) - LMT/ENS
  19. C Date         : Mars 2011
  20. C====&===1=========2=========3=========4=========5=========6=========7==
  21. C
  22. C-----DECLARATION GENERALE----------------------------------------------
  23. C
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25. C
  26. C-----DECLARATION GENERALE----------------------------------------------
  27. C
  28.       REAL*8 X(3,3),UNIT(3,3),EPSPI(3,3),EPSI(3,3),SIGMA(3,3)
  29.       REAL*8 EPSDEV(3,3),SIGMA1(3,3),EPSPIDEV(3,3)
  30.       REAL*8 EPSIPOS(3,3),EPSINEG(3,3),SIGPI(3,3),DEPSN(3,3)
  31.       REAL*8 TERME1(3,3),TERME2(3,3),TERME3(3,3),TERME4(3,3)
  32.       REAL*8 SIGMOX(3,3),TRAV(3,3),SIG1DEV(3,3),TERME5(3,3)
  33.       REAL*8 XDT,XTRAV2(3,3),XTRAV3(3,3),DEFN(3,3),EPSIT(3,3)
  34.       REAL*8 XMAT(*),SIG0(*),SIGF(*),VAR0(*),VARF(*),DEPST(*)
  35.       REAL*8 XDFP(3,3),XDGP(3,3),XDJ2P(3,3),DEVGP(3,3)
  36. C
  37. C-----DEFINITION DE PARAMETRES------------------------------------------
  38. C
  39.       PARAMETER (IMAX = 1000)
  40. C
  41. C-----MISE EN DONNEE----------------------------------------------------
  42. C
  43.       XE = XMAT(1)
  44.       XNU = XMAT(2)
  45.       XFT = XMAT(5)
  46.       ALDIR = XMAT(6)
  47.       GAMMA1 = XMAT(7)
  48.       A1 = XMAT(8)
  49.       SIGRF = XMAT(9)
  50.  
  51.       XAF = XMAT(10)
  52.       XBF = XMAT(11)
  53.       XAG = XMAT(12)
  54.       XBG = XMAT(13)
  55.       XAC = XMAT(14)
  56.       XBC = XMAT(15)
  57.       XSIGU = XMAT(16)
  58.       XFC = XMAT(17)
  59.  
  60.       XK = XE/((1.0D0-2.D0*XNU)*3.0D0)
  61.       XG = XE/(2.0D0*(1.D0+XNU))
  62. C
  63. C-----DEFINITION DE LA MATRICE UNITE D ORDRE 2------------------------------
  64. C
  65.           DO I=1,3
  66.          DO J=1,3
  67.             IF (I.EQ.J) THEN
  68.                UNIT(I,J) = 1.0D0
  69.             ELSE
  70.                UNIT(I,J) = 0.0D0
  71.             ENDIF
  72.          ENDDO
  73.       ENDDO
  74. C
  75. C-----VARIABLES INTERNES----------------------------------------------------
  76. C
  77.  
  78. C     Tenseur de glissement interne
  79.       EPSPI(1,1) = VAR0(2)
  80.       EPSPI(2,2) = VAR0(3)
  81.       EPSPI(3,3) = VAR0(4)
  82.       EPSPI(1,2) = VAR0(5)
  83.       EPSPI(1,3) = VAR0(6)
  84.       EPSPI(2,3) = VAR0(7)
  85.       EPSPI(2,1) = VAR0(5)
  86.       EPSPI(3,1) = VAR0(6)
  87.       EPSPI(3,2) = VAR0(7)
  88.  
  89. C     Tenseur de la contrainte de rappel
  90.       X(1,1) = VAR0(8)
  91.       X(2,2) = VAR0(9)
  92.       X(3,3) = VAR0(10)
  93.       X(1,2) = VAR0(11)
  94.       X(1,3) = VAR0(12)
  95.       X(2,3) = VAR0(13)
  96.       X(2,1) = VAR0(11)
  97.       X(3,1) = VAR0(12)
  98.       X(3,2) = VAR0(13)
  99.  
  100. C     Endommagement
  101.       XDT = VAR0(14)
  102.  
  103. C     Force associee a l ecrouissage isotrope
  104.       XFFT = VAR0(15)
  105.  
  106. C     Deformations totales
  107.       EPSI(1,1) = VAR0(16)
  108.       EPSI(2,2) = VAR0(17)
  109.       EPSI(3,3) = VAR0(18)
  110.       EPSI(1,2) = VAR0(19)
  111.       EPSI(1,3) = VAR0(20)
  112.       EPSI(2,3) = VAR0(21)
  113.       EPSI(2,1) = VAR0(19)
  114.       EPSI(3,1) = VAR0(20)
  115.       EPSI(3,2) = VAR0(21)
  116.  
  117. C     Tenseur des deformations plastiques
  118.       DEFN(1,1) = VAR0(35)
  119.       DEFN(2,2) = VAR0(36)
  120.       DEFN(3,3) = VAR0(37)
  121.       DEFN(1,2) = VAR0(38)
  122.       DEFN(1,3) = VAR0(39)
  123.       DEFN(2,3) = VAR0(40)
  124.       DEFN(2,1) = VAR0(38)
  125.       DEFN(3,1) = VAR0(39)
  126.       DEFN(3,2) = VAR0(40)
  127.  
  128. C     Ecrouissage plastique
  129.       XECO = VAR0(41)
  130.  
  131. C     Deformation plastique cumulee
  132.       XDPP = VAR0(42)
  133.  
  134. C     Contraintes tests
  135.       SIGMA1(1,1) = VAR0(22)
  136.       SIGMA1(2,2) = VAR0(23)
  137.       SIGMA1(3,3) = VAR0(24)
  138.       SIGMA1(1,2) = VAR0(25)
  139.       SIGMA1(1,3) = VAR0(26)
  140.       SIGMA1(2,3) = VAR0(27)
  141.       SIGMA1(2,1) = VAR0(25)
  142.       SIGMA1(3,1) = VAR0(26)
  143.       SIGMA1(3,2) = VAR0(27)
  144.  
  145. C     Indicateur
  146.       XBEN = VAR0(22)+VAR0(23)+VAR0(24)
  147.  
  148. C     Contrainte initiales
  149.       SIGMA(1,1) = SIG0(1)
  150.       SIGMA(2,2) = SIG0(2)
  151.       SIGMA(3,3) = SIG0(3)
  152.       SIGMA(1,2) = SIG0(4)
  153.       SIGMA(1,3) = SIG0(5)
  154.       SIGMA(2,3) = SIG0(6)
  155.       SIGMA(2,1) = SIG0(4)
  156.       SIGMA(3,1) = SIG0(5)
  157.       SIGMA(3,2) = SIG0(6)
  158.  
  159. C
  160. C-----ACTUALISATION DES DEFORMATIONS------------------------------------
  161. C
  162.       EPSI(1,1) = EPSI(1,1) + DEPST(1)
  163.       EPSI(2,2) = EPSI(2,2) + DEPST(2)
  164.       EPSI(3,3) = EPSI(3,3) + DEPST(3)
  165.       EPSI(1,2) = EPSI(1,2) + (DEPST(4)/2.0D0)
  166.       EPSI(1,3) = EPSI(1,3) + (DEPST(5)/2.0D0)
  167.       EPSI(2,3) = EPSI(2,3) + (DEPST(6)/2.0D0)
  168.       EPSI(2,1) = EPSI(2,1) + (DEPST(4)/2.0D0)
  169.       EPSI(3,1) = EPSI(3,1) + (DEPST(5)/2.0D0)
  170.       EPSI(3,2) = EPSI(3,2) + (DEPST(6)/2.0D0)
  171.  
  172. C
  173. C-----CALCUL DE QUELQUES QUANTITES--------------------------------------
  174. C
  175.  
  176. C     Calcul de la trace de EPSI
  177.       TRA = EPSI(1,1)+EPSI(2,2)+EPSI(3,3)
  178.  
  179. C     Calcul du deviateur de EPSI
  180.       CALL DEVIAT(EPSI,EPSDEV)
  181.  
  182. C     Calcul de la trace de EPSPI
  183.       TRAPI = EPSPI(1,1)+EPSPI(2,2)+EPSPI(3,3)
  184.  
  185. C     Calcul du deviateur de EPSPI
  186.       CALL DEVIAT(EPSPI,EPSPIDEV)
  187.  
  188. C     Calcul du deviateur de DEFN
  189.       CALL DEVIAT(DEFN,DEPSN)
  190.  
  191. C     Calcul de la trace de DEFN
  192.       TRADN = DEFN(1,1)+DEFN(2,2)+DEFN(3,3)
  193.  
  194. C
  195. C-----PARTIE POSITIVES ET NEGATIVES DE (EPS)-----------------------------
  196. C
  197.       DO I=1,3
  198.          DO J=1,3
  199.             EPSIT(I,J) = EPSI(I,J)-DEFN(I,J)
  200.          ENDDO
  201.       ENDDO
  202.  
  203.       CALL POSNEG(EPSIT,EPSIPOS,EPSINEG)
  204. C
  205. C-----SEUIL ENDOMMAGEMENT------------------------------------------------
  206. C
  207.  
  208. C     Calcul de l enerigie Y+
  209.       YP = 0.0D0
  210.  
  211.       DO I=1,3
  212.          DO J=1,3
  213.             YP = YP + EPSIPOS(I,J)*EPSIPOS(I,J)
  214.          ENDDO
  215.       ENDDO
  216.  
  217.       YP = 0.5D0*YP*XE
  218.  
  219. C
  220. C-----NON LOCAL----------------------------------------------------------
  221. C
  222.       IF (ISTEP.EQ.0) THEN
  223.          SEUIL1=YP
  224.          VARF(1)=YP
  225.       ELSE IF (ISTEP.EQ.1) THEN
  226.          VARF(1)=YP
  227.          DO I=2,42
  228.             VARF(I)=VAR0(I)
  229.          ENDDO
  230.          GOTO 2000
  231.       ELSE IF (ISTEP.EQ.2) THEN
  232.          SEUIL1=VAR0(1)
  233.          VARF(1)=SEUIL1
  234.       ENDIF
  235. C
  236. C-----EVOLUTION DE L ENDOMMAGEMENT---------------------------------------
  237. C
  238.  
  239. C     Calcul du seuil
  240.       XFD = SEUIL1 - ((XFT**2)/(2.0D0*XE) + XFFT)
  241.  
  242. C     Evolution explicite de l'endommagement
  243.       IF ((XFD.GT.0.0D0).AND.(XBEN.GE.0.0D0)) THEN
  244.          XDT  = 1.0D0-1.0D0/(1.0D0+ALDIR*(SEUIL1-((XFT**2)/(2.0D0*XE))))
  245.          XFFT = SEUIL1 - ((XFT**2)/(2.0D0*XE))
  246.       ENDIF
  247. C
  248. C-----SEUIL GLISSEMENT INTERNE-------------------------------------------
  249. C
  250.  
  251. C     Calcul de sigma_pi
  252.       DO I=1,3
  253.          DO J=1,3
  254.             SIGPI(I,J) = XDT*XK*(TRA-TRAPI-TRADN)*UNIT(I,J)+
  255.      &      2.0D0*XG*XDT*(EPSDEV(I,J)-EPSPIDEV(I,J)-DEPSN(I,J))
  256.          ENDDO
  257.       ENDDO
  258.  
  259. C     Calcul de la quantite (sigma_pi - X)
  260.       DO I=1,3
  261.          DO J=1,3
  262.             SIGMOX(I,J) = SIGPI(I,J)-X(I,J)
  263.          ENDDO
  264.       ENDDO
  265.  
  266. C     Calcul de I2(sigma_pi-X)
  267.       CALL I2(SIGMOX,XFPI)
  268.  
  269. C     Sigma_test = Sigma_finale
  270.       DO I=1,3
  271.          DO J=1,3
  272.             SIGMA(I,J)=SIGMA1(I,J)
  273.          ENDDO
  274.       ENDDO
  275.  
  276. C     Mise a 1 du multiplicateur plastique
  277.       XLAM = 1.0D0
  278.  
  279. C     CAlcul de la trace de Sigma_finale
  280.       TRASIG1 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  281.  
  282. C     Test sur la positivite du seuil de glissement interne
  283.       IF ((XFPI.GT.0.0D0).AND.(XBEN.GE.0.D0)) THEN
  284.  
  285. C         On retient la valeur initiale du seuil
  286.           XFPI0 = XFPI
  287.  
  288. C         On itere en bloquant le nombre d iteration max a IMAX
  289.           DO NI=1,IMAX
  290.  
  291. C             On realise un test de sortie
  292.               IF ((ABS(XFPI/XFPI0 ).LE.1.0D-8).OR.
  293.      &            (ABS(XLAM).LE.(1.0D-10))) THEN
  294.  
  295.                  GOTO 400
  296.  
  297.               ELSE
  298.  
  299. C                Calcul des derivees dfpi/dsigpi et dphipi/dsigpi
  300.                  CALL DI2(SIGMOX,TERME1)
  301.                  CALL DI2(SIGMOX,TERME2)
  302.  
  303. C                Calcul des derivees dfpi/dX et dphipi/dX
  304.                  DO I=1,3
  305.                     DO J=1,3
  306.                        TERME3(I,J) = 0.0D0
  307.                        TERME4(I,J) = 0.0D0
  308.  
  309.                        TERME3(I,J) = -1.0D0*TERME2(I,J)
  310.                        TERME4(I,J) = -1.0D0*TERME2(I,J)+A1*X(I,J)
  311.                     ENDDO
  312.                  ENDDO
  313.  
  314. C                On calcul le terme C:dfpi/dsigpi
  315.                  TRA_SE = TERME2(1,1)+TERME2(2,2)+TERME2(3,3)
  316.  
  317.                  CALL DEVIAT(TERME2,TRAV)
  318.  
  319.                  DO I=1,3
  320.                     DO J=1,3
  321.                        TERME5(I,J) = 0.0D0
  322.                        TERME5(I,J) = XK*TRA_SE*UNIT(I,J)+
  323.      &                 2.0D0*XG*TRAV(I,J)
  324.                     ENDDO
  325.                  ENDDO
  326.  
  327. C                On calcul le denominateur du multiplicateur plastique
  328.                  XT1 = 0.0D0
  329.                  DO I=1,3
  330.                     DO J=1,3
  331.                        XT1 = XT1 + TERME1(I,J)*TERME5(I,J)
  332.                     ENDDO
  333.                  ENDDO
  334.  
  335.                  XT2 = 0.0D0
  336.                  DO I=1,3
  337.                     DO J=1,3
  338.                        XT2 = XT2 + TERME3(I,J)*TERME4(I,J)
  339.                     ENDDO
  340.                  ENDDO
  341.  
  342. C                Calcul du multiplicateur plastique
  343.                  XLAM = XFPI/(XDT*XT1 + GAMMA1*XT2)
  344.  
  345. C                On actualisatise X et sigpi
  346.                  DO I=1,3
  347.                     DO J=1,3
  348.                        SIGPI(I,J)=SIGPI(I,J)-XDT*XLAM*TERME5(I,J)
  349.                        X(I,J)=X(I,J)-GAMMA1*XLAM*TERME4(I,J)
  350.                     ENDDO
  351.                  ENDDO
  352.  
  353. C                On recalcule le nouveau seuil
  354.                  DO I=1,3
  355.                     DO J=1,3
  356.                        SIGMOX(I,J) = 0.0D0
  357.                        SIGMOX(I,J) = SIGPI(I,J)-X(I,J)
  358.                     ENDDO
  359.                  ENDDO
  360.  
  361.                  CALL I2(SIGMOX,XFPI)
  362.  
  363.               ENDIF
  364.  
  365.           ENDDO
  366.  
  367.  400      CONTINUE
  368.  
  369.  
  370. C         On calcul epspi
  371.           TRA_SE=SIGPI(1,1)+SIGPI(2,2)+SIGPI(3,3)
  372.           CALL DEVIAT(SIGPI,TRAV)
  373.  
  374.           TRAPI = (TRA-TRADN) - (TRA_SE/(XDT*3.0D0*XK))
  375.  
  376.           DO I=1,3
  377.              DO J=1,3
  378.                 EPSPIDEV(I,J) = (EPSDEV(I,J)-DEPSN(I,J))-
  379.      &          (1.0D0/(2.0D0*XG*XDT))*TRAV(I,J)
  380.              ENDDO
  381.           ENDDO
  382.  
  383.           DO I=1,3
  384.              DO J=1,3
  385.                 EPSPI(I,J) = (TRAPI/3.0D0)*UNIT(I,J)+EPSPIDEV(I,J)
  386.              ENDDO
  387.           ENDDO
  388.  
  389.       ENDIF
  390. C
  391. C-----CONTRAINTES ACTUALISEE
  392. C
  393.       TRA1 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  394.  
  395.       IF (XBEN.GE.0.0D0) THEN
  396.          DO I=1,3
  397.             DO J=1,3
  398.                SIGMA1(I,J)=XK*(TRA-XDT*TRAPI-TRADN)*UNIT(I,J)+
  399.      &          2.0D0*XG*
  400.      &          (EPSDEV(I,J)-XDT*EPSPIDEV(I,J)-DEPSN(I,J))
  401.             ENDDO
  402.          ENDDO
  403.  
  404.       ELSEIF (XBEN.LE.SIGRF) THEN
  405.  
  406.  
  407. C        La contrainte effective
  408.          DO I=1,3
  409.             DO J=1,3
  410.                SIGMA1(I,J)=XK*(TRA-TRADN)*UNIT(I,J)+
  411.      &         2.0D0*XG*(EPSDEV(I,J)-DEPSN(I,J))
  412.             ENDDO
  413.          ENDDO
  414.  
  415. C        On verifie le seuile de plasticite
  416.          CALL J2(SIGMA1,TEMP1)
  417.  
  418.          TEMP2=SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  419.  
  420.          TEMP2 = TEMP2/3.00
  421.  
  422.          XFP = XAF*TEMP1+XBF*TEMP2 - (XFC + XECO)
  423.  
  424. C        On verifie la positivite du seuil
  425.          IF (XFP.GT.0.0D0) THEN
  426.  
  427. C           On retient la valeur initiale pour le critere relatif
  428.             XFP0 = XFP
  429.  
  430. C           On commence les choses serieuses
  431.             DO NP=1,IMAX
  432.  
  433. C           La derivee du critere
  434.             CALL DI2(SIGMA1,XDJ2P)
  435.  
  436.             DO I=1,3
  437.             DO J=1,3
  438.             XDFP(I,J) = XAF*XDJ2P(I,J)+XBF*UNIT(I,J)/3.0D0
  439.             XDGP(I,J) = XAG*XDJ2P(I,J)+XBG*UNIT(I,J)/3.0D0
  440.             ENDDO
  441.             ENDDO
  442.  
  443. C           Le hessien elasique
  444.             CALL DEVIAT(XDGP,DEVGP)
  445.  
  446.             XTGP = XDGP(1,1)+XDGP(2,2)+XDGP(3,3)
  447.  
  448.             DO I=1,3
  449.             DO J=1,3
  450.             TRAV(I,J)=0.0D0
  451.             TRAV(I,J)=XK*XTGP*UNIT(I,J)+2.0D0*XG*DEVGP(I,J)
  452.             ENDDO
  453.             ENDDO
  454.  
  455. C           La derivee de l ecrouissage
  456.             XTEMP= EXP(-1.0D0*XBC*XDPP)*(XAC-(XAC*XDPP+XFC+XSIGU)*XBC)
  457.  
  458. C           Le critere actuel
  459.  
  460.             CALL J2(SIGMA1,TEMP1)
  461.  
  462.             TEMP2=SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  463.  
  464.             TEMP2 = TEMP2/3.00
  465.  
  466.             XFP = XAF*TEMP1+XBF*TEMP2 - (XFC + XECO)
  467.  
  468. C           Le multiplicateur plastique
  469.             CALL DBLECONT(XDFP,TRAV,TEMP9)
  470.  
  471.             XLAMP = XFP/(TEMP9+XTEMP)
  472.  
  473. C           Actualisation des variables internes
  474.             XECO = XECO + XLAMP*XTEMP
  475.  
  476.             XDPP = XDPP + XLAMP
  477.  
  478.             DO I=1,3
  479.             DO J=1,3
  480.             DEFN(I,J)=DEFN(I,J) + XLAMP*XDGP(I,J)
  481.             ENDDO
  482.             ENDDO
  483.  
  484.             DO I=1,3
  485.             DO J=1,3
  486.             SIGMA1(I,J)=SIGMA1(I,J) - XLAMP*TRAV(I,J)
  487.             ENDDO
  488.             ENDDO
  489.  
  490. C           On calcul le seuil pour un eventuel arret
  491.             CALL J2(SIGMA1,TEMP1)
  492.  
  493.             TEMP2=SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  494.  
  495.             TEMP2 = TEMP2/3.00
  496.  
  497.             XFP = XAF*TEMP1+XBF*TEMP2 - (XFC + XECO)
  498.  
  499. C           Le critere d arret des iterations internes
  500.             IF ((ABS(XFP/XFP0).LE.(1.0D-8)).OR.(NP.GE.IMAX)) THEN
  501.  
  502.                GOTO 294
  503.  
  504.             ENDIF
  505.  
  506.  
  507.             ENDDO
  508.  
  509.  294     CONTINUE
  510.  
  511.          ENDIF
  512.  
  513.       ELSE
  514.  
  515. C        Calcul du deviateur de DEFN
  516.          CALL DEVIAT(DEFN,DEPSN)
  517.  
  518. C        Calcul de la trace de DEFN
  519.          TRADN = DEFN(1,1)+DEFN(2,2)+DEFN(3,3)
  520.  
  521.          TRASIG1 = 3.0D0*XK*(TRA-XDT*TRAPI-TRADN)/
  522.      &             (1.0D0-(3.0D0*XK/SIGRF)*XDT*TRAPI)
  523.  
  524.          DO I=1,3
  525.             DO J=1,3
  526.                SIG1DEV(I,J)=2.00D0*XG*(EPSDEV(I,J)-XDT*EPSPIDEV(I,J)*
  527.      &         (1.0D0-(TRASIG1/SIGRF))-DEPSN(I,J))
  528.             ENDDO
  529.          ENDDO
  530.  
  531.          DO I=1,3
  532.             DO J=1,3
  533.                SIGMA1(I,J)=TRASIG1*UNIT(I,J)/3.0D0+SIG1DEV(I,J)
  534.             ENDDO
  535.          ENDDO
  536.  
  537.       ENDIF
  538.  
  539. C
  540. C-----CALCUL DE L ENDOMMAGEMENT EFFECTIF (GRANDEUR INDICATIVE)
  541. C-----PAS DE SENS THERMO
  542. C
  543.       XFSIG = 1.0D0-(SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3))/SIGRF
  544.  
  545.       IF ((XFSIG.LE.1.0D0).AND.(XFSIG.GE.0.0D0)) THEN
  546.           XDEFF = XDT*XFSIG
  547.       ELSE
  548.           XDEFF = XDT
  549.       ENDIF
  550.  
  551. C
  552. C-----STOCKAGE DES CONTRAINTES AVEC SYMMETRISATION POUR ROBUSTESSE
  553. C
  554.       SIGF(1)= SIGMA1(1,1)
  555.       SIGF(2)= SIGMA1(2,2)
  556.       SIGF(3)= SIGMA1(3,3)
  557.       SIGF(4)= (SIGMA1(1,2)+SIGMA1(2,1))/2.0D0
  558.       SIGF(5)= (SIGMA1(1,3)+SIGMA1(3,1))/2.0D0
  559.       SIGF(6)= (SIGMA1(2,3)+SIGMA1(3,2))/2.0D0
  560. C
  561. C-----STOCKAGE DES DEFORMATIONS DE GLISSEMENT
  562. C
  563.       VARF(2)=EPSPI(1,1)
  564.       VARF(3)=EPSPI(2,2)
  565.       VARF(4)=EPSPI(3,3)
  566.       VARF(5)=EPSPI(1,2)
  567.       VARF(6)=EPSPI(1,3)
  568.       VARF(7)=EPSPI(2,3)
  569. C
  570. C-----STOCKAGE DE L ECROUISSAGE CINEMATIQUE
  571. C
  572.       VARF(8)=X(1,1)
  573.       VARF(9)=X(2,2)
  574.       VARF(10)=X(3,3)
  575.       VARF(11)=X(1,2)
  576.       VARF(12)=X(1,3)
  577.       VARF(13)=X(2,3)
  578. C
  579. C-----STOCKAGE DES CONTRAINTES TESTS
  580. C
  581.       VARF(22)=SIGMA1(1,1)
  582.       VARF(23)=SIGMA1(2,2)
  583.       VARF(24)=SIGMA1(3,3)
  584.       VARF(25)=SIGMA1(1,2)
  585.       VARF(26)=SIGMA1(1,3)
  586.       VARF(27)=SIGMA1(2,3)
  587. C
  588. C-----TENSEUR DES DEFORMATIONS PLASTIQUES
  589. C
  590.       VARF(35) = DEFN(1,1)
  591.       VARF(36) = DEFN(2,2)
  592.       VARF(37) = DEFN(3,3)
  593.       VARF(38) = DEFN(1,2)
  594.       VARF(39) = DEFN(1,3)
  595.       VARF(40) = DEFN(2,3)
  596. C
  597. C-----Ecrouissage plastique
  598. C
  599.       VARF(41) = XECO
  600. C
  601. C-----Deformation plastique cumulee
  602. C
  603.       VARF(42) = XDPP
  604. C
  605. C-----STOCKAGE DES CONTRAINTES TESTS
  606. C
  607.       VARF(29)=SIGPI(1,1)
  608.       VARF(30)=SIGPI(2,2)
  609.       VARF(31)=SIGPI(3,3)
  610.       VARF(32)=SIGPI(1,2)
  611.       VARF(33)=SIGPI(1,3)
  612.       VARF(34)=SIGPI(2,3)
  613. C
  614. C-----ENDOMMAGEMENT EFFECTIF
  615. C
  616.       VARF(28) = XDEFF
  617. C
  618. C-----STOCKAGE DE L ENDOMMAGEMENT
  619. C
  620.       VARF(14) = XDT
  621. C
  622. C-----STOCKAGE DES FORCES THERMODYNAMIQUES ASSOCIEES A L ECROUISSAGE ISOTROPE
  623. C
  624.       VARF(15) = XFFT
  625. C
  626. C-----STOCKAGE DES DEFORMATIONS TOTALES
  627. C
  628.       VARF(16)=EPSI(1,1)
  629.       VARF(17)=EPSI(2,2)
  630.       VARF(18)=EPSI(3,3)
  631.       VARF(19)=EPSI(1,2)
  632.       VARF(20)=EPSI(1,3)
  633.       VARF(21)=EPSI(2,3)
  634. C
  635. C-----BALISE DE SORTIE POUR LE NON LOCAL
  636. C
  637.  2000 CONTINUE
  638.  
  639.       RETURN
  640.           END
  641.  
  642.  
  643.  
  644.  
  645.  
  646.  
  647.  
  648.  
  649.  
  650.  
  651.  
  652.  
  653.  
  654.  
  655.  
  656.  
  657.  
  658.  
  659.  
  660.  
  661.  
  662.  
  663.  
  664.  

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