Télécharger la plaquette Cast3M


fr en

Télécharger ric3nl.eso

Retour à la liste

Numérotation des lignes : 

ric3nl

  1. C RIC3NL    SOURCE    OF166741  25/11/04    21:16:04     12349          
  2. C sub ricnl3d
  3.        SUBROUTINE RIC3NL(wrk52,wrk53,wrk54,nvari,iecou)
  4. C
  5. C====&===1=========2=========3=========4=========5=========6=========7==
  6. C Commentaires : Subroutine permettant de mettre en oeuvre le
  7. C                    modele RICRAG (F. Ragueneau & B. Richard) pour représenter
  8. C                le comportement 3D d'un béton ordinaire sous
  9. C                sollicitations complexes en non local
  10. C Traits       : - Endommagement scalaire
  11. C                - Glissement frottant
  12. C                - 3D
  13. C                - Monotone et cyclique
  14. C                - Non localité portant sur le seuil d'endommagement (Fd)
  15. C                - Bon comportement unilatteral
  16. C Auteurs      : B. Richard (doctorant) & F. Ragueneau (LCPC, FDOA & BCC)
  17. C Date         : Mai 2008
  18. C====&===1=========2=========3=========4=========5=========6=========7==
  19. C
  20. C-----DECLARATION GENERALE----------------------------------------------
  21. C
  22.       IMPLICIT REAL*8(A-H,O-Z)
  23. C
  24. C-----LISTE DES INCLUDES A CHARGER--------------------------------------
  25. C
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29. -INC DECHE
  30. -INC CCREEL
  31. C
  32. C-----DECLARATION GENERALE----------------------------------------------
  33. C
  34. C       REAL*8 XMAT(9),VAR0(30),VARF(30),SIGF(6),DEPST(6)
  35.       REAL*8 XZER,UN,DEUX,EPSI(3,3)
  36.       REAL*8 E,NU,FT,ALIND,GAMMA1,A1,ALDIR
  37.       REAL*8 EPS0,K,G,Y0,DEPSDEV(3,3),DTRA
  38.       REAL*8 UNIT(3,3),DEPSI(3,3),EPSTFBIS(6)
  39.       REAL*8 EPSIPI(3,3),ALPHA(3,3),Z,D,DONMDNL,DONMINL,FFFZ,DOLD,FORCED
  40.       REAL*8 FORCEI,NRUPT,SEUILT,DEVPI(3,3)
  41.       REAL*8 TRA,EPSDEV(3,3),SIGMA(3,3),VPEPS(3),SEUIL1NL
  42.       REAL*8 YD,YN,EPSDIR(3,3),EPSIND(3,3),Dn
  43.       REAL*8 SIGMPI(3,3),SIGMOD(3,3),SIGMOX(3,3),X(3,3),TERME6(3,3)
  44.       REAL*8 TERME7(3,3),TERME8(3,3),TERME9(3,3),TROGO,DEVTRO(3,3)
  45.       REAL*8 DEVDGS(3,3)
  46. C       CHARACTER*8 CMATE
  47. C
  48. C-----DEFINITION DE PARAMETRES------------------------------------------
  49. C
  50.       PARAMETER (XZER=0.0D0 , UN=1.0D0 , DEUX=2.0D0, XCOEF=0.00D0)
  51. C
  52. C-----OUVERTURE DE SEGMENTS DE TRAVAIL----------------------------------
  53. C
  54. C-----MISE EN DONNEE----------------------------------------------------
  55. C
  56.       CMATE='ISOTROPE'
  57.       E      = XMAT(1)
  58.       NU     = XMAT(2)
  59.       FT     = XMAT(5)
  60.       ALIND  = XMAT(6)
  61.       GAMMA1 = XMAT(7)
  62.       A1     = XMAT(8)
  63.       ALDIR  = XMAT(9)
  64.  
  65.       EPS0   = FT/E
  66.       K      = E/(1.D0-2.D0*NU)
  67.       G      = E/(2.D0*(1.D0+NU))
  68.       Y0     = (K/6.0D0*EPS0**2.D0)
  69.  
  70.         DO J=1,3
  71.          DO I=1,3
  72.             IF (I.EQ.J) THEN
  73.                UNIT(I,J)=UN
  74.             ELSE
  75.                UNIT(I,J)=XZER
  76.             ENDIF
  77.          ENDDO
  78.       ENDDO
  79. C
  80. C-----PASSAGE DES VARIABLES INTERNES------------------------------------
  81. C
  82. C     GLISSEMENT (EPSILON_PI)
  83. C
  84.       EPSIPI(1,1)=VAR0(2)
  85.       EPSIPI(2,2)=VAR0(3)
  86.       EPSIPI(3,3)=VAR0(4)
  87.  
  88.       EPSIPI(1,2)=VAR0(5)
  89.       EPSIPI(1,3)=VAR0(6)
  90.       EPSIPI(2,3)=VAR0(7)
  91.  
  92.       EPSIPI(2,1)=VAR0(5)
  93.       EPSIPI(3,1)=VAR0(6)
  94.       EPSIPI(3,2)=VAR0(7)
  95. C
  96. C     ECROUISSAGE CINEMATIQUE (ALPHA)
  97. C
  98.       ALPHA(1,1)=VAR0(8)
  99.       ALPHA(2,2)=VAR0(9)
  100.       ALPHA(3,3)=VAR0(10)
  101. C
  102.       ALPHA(1,2)=VAR0(11)
  103.       ALPHA(1,3)=VAR0(12)
  104.       ALPHA(2,3)=VAR0(13)
  105. C
  106.       ALPHA(2,1)=VAR0(11)
  107.       ALPHA(3,1)=VAR0(12)
  108.       ALPHA(3,2)=VAR0(13)
  109. C
  110. C     ENDOMMAGEMENT (D)
  111. C
  112.       D=VAR0(14)
  113. C
  114. C     FORCE ECROUISSAGE ISOTROPE
  115. C
  116.       FFFZ=VAR0(17)
  117. C
  118. C     DEFORMATION TOTALE MISE DANS EPSI
  119. C
  120.       DO I=1,6
  121.            EPSTFBIS(I)=VAR0(18+I)+DEPST(I)
  122.       ENDDO
  123.  
  124.       DO I=1,3
  125.            EPSI (I,I)=VAR0(18+I)+DEPST(I)
  126.            DEPSI(I,I)=DEPST(I)
  127.       ENDDO
  128.  
  129.       EPSI(1,2)=VAR0(22)+DEPST(4)*0.5D0
  130.       EPSI(1,3)=VAR0(23)+DEPST(5)*0.5D0
  131.       EPSI(2,3)=VAR0(24)+DEPST(6)*0.5D0
  132.       EPSI(2,1)=VAR0(22)+DEPST(4)*0.5D0
  133.       EPSI(3,1)=VAR0(23)+DEPST(5)*0.5D0
  134.       EPSI(3,2)=VAR0(24)+DEPST(6)*0.5D0
  135.  
  136.       DEPSI(1,2)=DEPST(4)*0.5D0
  137.       DEPSI(1,3)=DEPST(5)*0.5D0
  138.       DEPSI(2,3)=DEPST(6)*0.5D0
  139.       DEPSI(2,1)=DEPST(4)*0.5D0
  140.       DEPSI(3,1)=DEPST(5)*0.5D0
  141.       DEPSI(3,2)=DEPST(6)*0.5D0
  142.  
  143.       SIGMA(1,1)=SIG0(1)
  144.       SIGMA(2,2)=SIG0(2)
  145.       SIGMA(3,3)=SIG0(3)
  146.       SIGMA(1,2)=SIG0(4)
  147.       SIGMA(1,3)=SIG0(5)
  148.       SIGMA(2,3)=SIG0(6)
  149.       SIGMA(2,1)=SIG0(4)
  150.       SIGMA(3,1)=SIG0(5)
  151.       SIGMA(3,2)=SIG0(6)
  152. C
  153. C     ECROUISSAGE
  154. C
  155.       FORCED = VAR0(25)
  156.       FORCEI = VAR0(26)
  157. C
  158. C     MATRICE DE HOOKE
  159. C
  160.       AUX0=E/((UN+NU)*(UN-NU-NU))
  161.       AUX=AUX0*(UN-NU)
  162.       AUX1=AUX0*NU
  163.       XGEGE=0.5D0*E/(UN+NU)
  164.  
  165.       DDHOOK(1,1)=AUX
  166.       DDHOOK(2,1)=AUX1
  167.       DDHOOK(3,1)=AUX1
  168.  
  169.       DDHOOK(1,2)=AUX1
  170.       DDHOOK(2,2)=AUX
  171.       DDHOOK(3,2)=AUX1
  172.  
  173.       DDHOOK(1,3)=AUX1
  174.       DDHOOK(2,3)=AUX1
  175.       DDHOOK(3,3)=AUX
  176.  
  177.       DDHOOK(4,4)=XGEGE
  178.       DDHOOK(5,5)=XGEGE
  179.       DDHOOK(6,6)=XGEGE
  180. C-----------------------------------------------------------------------
  181. C-----DEBUT INTEGRATION NUMERIQUE---------------------------------------
  182. C-----------------------------------------------------------------------
  183.       TRA     = (EPSI(1,1)+EPSI(2,2)+EPSI(3,3))/3.D0
  184.       DTRA    = (DEPSI(1,1)+DEPSI(2,2)+DEPSI(3,3))/3.D0
  185.  
  186.         DO J=1,3
  187.          DO I=1,3
  188.             EPSDEV(I,J)  = EPSI(I,J)-TRA*UNIT(I,J)
  189.             DEPSDEV(I,J) = DEPSI(I,J)-DTRA*UNIT(I,J)
  190.          ENDDO
  191.       ENDDO
  192.  
  193.         DO J=1,3
  194.          DO I=1,3
  195.                 SIGMA(I,J)=SIGMA(I,J)+K*DTRA+2.0D0*G*DEPSDEV(I,J)
  196.          ENDDO
  197.       ENDDO
  198. C
  199. C-----ENDOMMAGEMENT------------------------------------------------------
  200. C
  201. C       CALCUL DES VALEURS PRINCIPALES
  202. C
  203.         CALL JACOD3(EPSI ,3,VPEPS)
  204. C
  205. C       CALCUL DE EPSILON DIRECT
  206. C
  207.         IF (TRA.GE.XZER) THEN
  208.          XPDC=1.0D0
  209. c          GAMMA1 = 7.0D8
  210. c          A1     = 7.0D-5
  211.       ELSE
  212.          XPDC=XZER
  213.       ENDIF
  214.  
  215.         DO J=1,3
  216.          DO I=1,3
  217.             IF (I.EQ.J) THEN
  218.                EPSDIR(I,J)=0.5D0*(ABS(VPEPS(I))+VPEPS(I))*XPDC
  219.             ELSE
  220.                EPSDIR(I,J)=XZER
  221.             ENDIF
  222.          ENDDO
  223.       ENDDO
  224. C
  225. C       CALCUL DE EPSILON INDUIT
  226. C
  227.         DO J=1,3
  228.          DO I=1,3
  229.             IF (I.EQ.J) THEN
  230.                EPSIND(I,J)=0.5D0*(ABS(VPEPS(I))+VPEPS(I))-EPSDIR(I,J)
  231.             ELSE
  232.                EPSIND(I,J)=XZER
  233.             ENDIF
  234.          ENDDO
  235.       ENDDO
  236. C
  237. C       TAUX D ENERGIES CORRESPONDANTES
  238. C
  239.       DONMDNL=XZER
  240.       DONMINL=XZER
  241.  
  242.       DO J=1,3
  243.            DO I=1,3
  244.               DONMDNL=DONMDNL+EPSDIR(I,J)*EPSDIR(I,J)
  245.               DONMINL=DONMINL+EPSIND(I,J)*EPSIND(I,J)
  246.            ENDDO
  247.       ENDDO
  248.  
  249.         SEUIL1NL=(K/6.0D0) * (DONMDNL*XPDC+DONMINL*(1.0D0-XPDC))
  250.  
  251.       IF (ISTEP.EQ.0) THEN
  252.          SEUIL1=SEUIL1NL
  253.          VARF(1)=SEUIL1
  254.       ELSE IF (ISTEP.EQ.1) THEN
  255.          VARF(1)=SEUIL1NL
  256.          DO I=2,26
  257.           VARF(I)=VAR0(I)
  258.          ENDDO
  259.          GOTO 2000
  260.       ELSE IF (ISTEP.EQ.2) THEN
  261.          SEUIL1=VAR0(1)
  262.          VARF(1)=SEUIL1
  263.       ENDIF
  264. C
  265. C       CALCUL NON LOCAL - ON MOYENNE LE SEUIL
  266. C
  267.       FFFZ=FORCED*XPDC+FORCEI*(1.0D0-XPDC)
  268.         SEUILT=SEUIL1-(FFFZ+Y0)
  269. C
  270. C       EVOLUTION ENDOMMAGEMENT
  271. C
  272.         IF (SEUILT.GT.XZER) THEN
  273.            IF (TRA.GE.XZER) THEN
  274.               Dn      = 1.0D0-1.0D0/(1.0D0+ALDIR*(SEUIL1-Y0))
  275.               FORCED = SEUIL1-Y0
  276.            ELSE
  277.               Dn      = 1.0D0-1.0D0/(1.0D0+ALIND*(SEUIL1-Y0))
  278.               FORCEI = SEUIL1-Y0
  279.            ENDIF
  280. C
  281. C     ON VERIFIE D COISSANT
  282. C
  283.         IF (Dn.GE.D) THEN
  284.            D=Dn
  285.         ENDIF
  286. C
  287. C     ON LIMITE LA VALEUR DE D
  288. C
  289.       IF (D.GE.0.99999D0) THEN
  290.          D=0.99999D0
  291.       ENDIF
  292. C
  293. C       ON BLOQUE LES ECROUISSAGES
  294. C
  295.            FFFZ=FORCED*XPDC+FORCEI*(1.0D0-XPDC)
  296.         ENDIF
  297.  
  298.       DO J=1,3
  299.            DO I=1,3
  300.               X(I,J)=GAMMA1*ALPHA(I,J)
  301.            ENDDO
  302.       ENDDO
  303. C
  304. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  305. C
  306. C                               SUBSTEPPING
  307. C
  308. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  309. C
  310.         NENTIER=100
  311.         DO KIN=1,NENTIER
  312.  
  313.         DO I=1,6
  314.            EPSTFBIS(I)=VAR0(18+I)+KIN*DEPST(I)/NENTIER
  315.       ENDDO
  316.  
  317.         EPSI(1,1)=EPSTFBIS(1)
  318.         EPSI(2,2)=EPSTFBIS(2)
  319.         EPSI(3,3)=EPSTFBIS(3)
  320.         EPSI(1,2)=EPSTFBIS(4)/2.0D0
  321.         EPSI(1,3)=EPSTFBIS(5)/2.0D0
  322.         EPSI(2,3)=EPSTFBIS(6)/2.0D0
  323.         EPSI(2,1)=EPSTFBIS(4)/2.0D0
  324.         EPSI(3,1)=EPSTFBIS(5)/2.0D0
  325.         EPSI(3,2)=EPSTFBIS(6)/2.0D0
  326.  
  327.         TRA   = (EPSI(1,1)+EPSI(2,2)+EPSI(3,3))/3.D0
  328.         TRAPI = (EPSIPI(1,1)+EPSIPI(2,2)+EPSIPI(3,3))/3.D0
  329.  
  330.         DO J=1,3
  331.           DO I=1,3
  332.             EPSDEV(I,J)=EPSI(I,J)-TRA*UNIT(I,J)
  333.              DEVPI(I,J)=EPSIPI(I,J)-TRAPI*UNIT(I,J)
  334.           ENDDO
  335.       ENDDO
  336.  
  337.         DO J=1,3
  338.          DO I=1,3
  339.                 SIGMA(I,J)  = (1.0D0-D)*K*TRA*UNIT(I,J)+
  340.      &                    (1.0D0-D)*2.0D0*G*EPSDEV(I,J)+
  341.      &                    K*D*(TRA-TRAPI)*UNIT(I,J)+
  342.      &                2.0D0*G*D*(EPSDEV(I,J)-DEVPI(I,J))
  343.          ENDDO
  344.       ENDDO
  345.  
  346.          DO J=1,3
  347.           DO I=1,3
  348.             SIGMPI(I,J)=K*D*(TRA-TRAPI)*UNIT(I,J)+
  349.      &              2.0D0*G*D*(EPSDEV(I,J)-DEVPI(I,J))
  350.           ENDDO
  351.        ENDDO
  352. C
  353. C====&===1=========2=========3=========4=========5=========6=========7==
  354. C-----DEBUT INTEGRATION GLISSEMENT--------------------------------------
  355. C====&===1=========2=========3=========4=========5=========6=========7==
  356. C
  357.         DO 500 IREMP=1,500
  358. C
  359. C     1) EVALUATION DU SEUIL
  360. C     ----------------------
  361.         DO 101 J=1,3
  362.           DO 201 I=1,3
  363.             SIGMOX(I,J)=XZER
  364.             SIGMOX(I,J)=SIGMPI(I,J)-X(I,J)
  365.  201      CONTINUE
  366.  101    CONTINUE
  367.  
  368.  
  369.         TRASX= (SIGMOX(1,1)+SIGMOX(2,2)+SIGMOX(3,3))/3.D0
  370.  
  371.         DO 102 J=1,3
  372.           DO 202 I=1,3
  373.             SIGMOD(I,J)=XZER
  374.             SIGMOD(I,J)=SIGMOX(I,J)-TRASX*UNIT(I,J)
  375.  202      CONTINUE
  376.  102    CONTINUE
  377.  
  378.         TERME4=XZER
  379.         DO 103 J=1,3
  380.           DO 203 I=1,3
  381.             TERME4=TERME4+SIGMOD(I,J)*SIGMOD(I,J)
  382.  203      CONTINUE
  383.  103    CONTINUE
  384. C
  385. C         CALCUL DU J2(SIGMPI-X)
  386. C
  387.         TRASIN = SIGMA(1,1)+SIGMA(2,2)+SIGMA(3,3)
  388.         TRASIN = 0.5D0*(-ABS(TRASIN)+TRASIN)*0.0D0
  389.  
  390.         TERME4=SQRT(1.5D0*TERME4)+(TRASIN/3.0D0)*0.80D0
  391.  
  392.         SEUIL2=TERME4+XCOEF*TRASX
  393.         IF (IREMP.EQ.1) THEN
  394.           SEUINI=SEUIL2
  395.           DLAM2=XZER
  396.           COCO=2
  397.         ENDIF
  398.  
  399.         IF (SEUINI.GT.1.0D0) THEN
  400.           CRIT=SEUIL2/SEUINI
  401.         ELSE
  402.           CRIT=XZER
  403.         ENDIF
  404. C
  405. C       2) ON VERIFIE LA VALEUR RELATIVE DU SEUIL
  406. C       -----------------------------------------
  407.         IF ((CRIT.LE.1.0D-5).OR.(DLAM2.LE.1.0D-10.AND.IREMP.GT.1).OR.
  408.      &       (D.EQ.XZER))THEN
  409.           GOTO 666
  410.         ELSE
  411.         IF (IREMP.EQ.1) COCO=1
  412. C
  413. C       3) CALCUL DES DERIVÉES POUR LA REDESCENTE
  414. C       -----------------------------------------
  415.         DO J=1,3
  416.           DO I=1,3
  417.             TERME6(I,J)=XZER
  418.             TERME7(I,J)=XZER
  419.             TERME8(I,J)=XZER
  420.             TERME9(I,J)=XZER
  421.             ENDDO
  422.           ENDDO
  423. C
  424. C       CALCUL DE DF/DS....................................................
  425. C
  426.         DO J=1,3
  427.           DO I=1,3
  428.             TERME6(I,J)=1.5D0*SIGMOX(I,J)/TERME4+XCOEF*1.0D0/3.D0*
  429.      &                  UNIT(I,J)
  430.             ENDDO
  431.           ENDDO
  432. C
  433. C       CALCUL DE DF/DX....................................................
  434. C
  435.         DO J=1,3
  436.           DO I=1,3
  437.             TERME7(I,J)=-TERME6(I,J)
  438.             ENDDO
  439.           ENDDO
  440. C
  441. C       CALCUL DE DG/DS....................................................
  442. C
  443.         DO J=1,3
  444.           DO I=1,3
  445.             TERME8(I,J)=1.5D0*SIGMOX(I,J)/TERME4+XCOEF*1.0D0/3.D0*
  446.      &                  UNIT(I,J)
  447.             ENDDO
  448.           ENDDO
  449. C
  450. C       CALCUL DE DG/DX....................................................
  451. C
  452.         DO J=1,3
  453.           DO I=1,3
  454.             TERME9(I,J)=-1.5D0*SIGMOX(I,J)/TERME4+A1*X(I,J)-
  455.      &                  XCOEF*1.0D0/3.D0*UNIT(I,J)
  456.             ENDDO
  457.           ENDDO
  458. C
  459. C       4) CALCUL DU MULTIPLICATEUR DE GLISSEMENT
  460. C       -----------------------------------------
  461.           TRDGZ=(TERME8(1,1)+TERME8(2,2)+TERME8(3,3))/3.D0
  462.  
  463.           DO J=1,3
  464.           DO I=1,3
  465.             DEVDGS(I,J)=TERME8(I,J)-TRDGZ*UNIT(I,J)
  466.             ENDDO
  467.           ENDDO
  468.  
  469.         TERM10=XZER
  470.         DO J=1,3
  471.           DO I=1,3
  472.             TERM10=TERM10+TERME6(I,J)*(K*TRDGZ*UNIT(I,J)+
  473.      &             2.0D0*G*DEVDGS(I,J))
  474.             ENDDO
  475.           ENDDO
  476.  
  477.         TERM11=XZER
  478.         DO J=1,3
  479.           DO I=1,3
  480.             TERM11=TERM11+TERME7(I,J)*TERME9(I,J)
  481.             ENDDO
  482.           ENDDO
  483.  
  484.         DLAM2=SEUIL2/(TERM10+GAMMA1*TERM11)
  485. C
  486. C       5) ACTUALISATION DES VARIABLES FORCES
  487. C       -------------------------------------
  488.         DO J=1,3
  489.           DO I=1,3
  490.             SIGMPI(I,J)=SIGMPI(I,J)-DLAM2*D*(K*TRDGZ*UNIT(I,J)+
  491.      &             2.0D0*G*DEVDGS(I,J))
  492.             ENDDO
  493.           ENDDO
  494.  
  495.         DO J=1,3
  496.           DO I=1,3
  497.             X(I,J)=X(I,J)-GAMMA1*DLAM2*TERME9(I,J)
  498.             ENDDO
  499.           ENDDO
  500.  
  501.         ENDIF
  502.  500  CONTINUE
  503. C
  504. C-----ACTUALISATION DES VARIABLES FLUX (S'IL Y A LIEU)------------------
  505. C
  506. C     1) ALPHA :
  507. C       -------
  508.  666  CONTINUE
  509.       IF (COCO.EQ.1) THEN
  510.         DO 303 J=1,3
  511.           DO 403 I=1,3
  512.             ALPHA(I,J)=X(I,J)/GAMMA1
  513.  403      CONTINUE
  514.  303    CONTINUE
  515. C
  516. C     2) EPSILON_PI :
  517. C        ------------
  518.         dplus=max(d,xpetit)
  519.         TRA_Z = (SIGMPI(1,1)+SIGMPI(2,2)+SIGMPI(3,3))/3.D0
  520.         TROGO = TRA - TRA_Z/(Dplus*K)
  521.         DO J=1,3
  522.           DO I=1,3
  523.             DEVTRO(I,J)=EPSDEV(I,J)-(SIGMPI(I,J)-TRA_Z*UNIT(I,J))
  524.      &                                /(2.0D0*G*Dplus)
  525.             ENDDO
  526.         ENDDO
  527.  
  528.         DO J=1,3
  529.           DO I=1,3
  530.             EPSIPI(I,J)=TROGO*UNIT(I,J)+DEVTRO(I,J)
  531.             ENDDO
  532.         ENDDO
  533.  
  534.       ENDIF
  535. C
  536. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  537. C
  538.         ENDDO
  539. C
  540. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  541. C
  542. C
  543. C-----CORRECTION DES CONTRAINTES
  544. C
  545.       TRAPOS = ( ABS(TRA)+TRA)*0.5D0
  546.       TRANEG = (-ABS(TRA)+TRA)*0.5D0
  547.  
  548.        DO J=1,3
  549.          DO I=1,3
  550.               SIGMA(I,J)  = (1.0D0-D)*K*TRAPOS*UNIT(I,J)+
  551.      &                              K*TRANEG*UNIT(I,J)+
  552.      &                          (1.0D0-D)*2.0D0*G*EPSDEV(I,J)+
  553.      &                           SIGMPI(I,J)
  554.          ENDDO
  555.       ENDDO
  556. C------------------------------------------------------------------------
  557. C-----FIN INTEGRATION NUMERIQUE------------------------------------------
  558. C------------------------------------------------------------------------
  559.       SIGF(1)= SIGMA(1,1)
  560.       SIGF(2)= SIGMA(2,2)
  561.       SIGF(3)= SIGMA(3,3)
  562.       SIGF(4)= (SIGMA(1,2)+SIGMA(2,1))*0.5D0
  563.       SIGF(5)= (SIGMA(1,3)+SIGMA(3,1))*0.5D0
  564.       SIGF(6)= (SIGMA(2,3)+SIGMA(3,2))*0.5D0
  565. C
  566. C     GLISSEMENT (EPSILON_PI)
  567. C
  568.       VARF(2)=EPSIPI(1,1)
  569.       VARF(3)=EPSIPI(2,2)
  570.       VARF(4)=EPSIPI(3,3)
  571.       VARF(5)=EPSIPI(1,2)
  572.       VARF(6)=EPSIPI(1,3)
  573.       VARF(7)=EPSIPI(2,3)
  574. C
  575. C     ECROUISSAGE CINEMATIQUE (ALPHA)
  576. C
  577.       VARF(8 )=ALPHA(1,1)
  578.       VARF(9 )=ALPHA(2,2)
  579.       VARF(10)=ALPHA(3,3)
  580.       VARF(11)=ALPHA(1,2)
  581.       VARF(12)=ALPHA(1,3)
  582.       VARF(13)=ALPHA(2,3)
  583. C
  584. C     ENDOMMAGEMENT (D)
  585. C
  586.       VARF(14)=D
  587. C
  588. C     FORCE LIEE A L ECROUISSAGE ISOTROPE
  589. C
  590.       VARF(17)=FFFZ
  591. C
  592. C     DEFORMATION TOTALE
  593. C
  594.       DO I=1,6
  595.          VARF(18+I)=EPSTFBIS(I)
  596.         ENDDO
  597. C
  598. C     ECROUISSAGES
  599. C
  600.         VARF(25)=FORCED
  601.         VARF(26)=FORCEI
  602. C
  603. C     CONTRAINTES DE FROTTEMENT
  604. C
  605. *      VARF(27)=SIGMPI(1,1)
  606. *      VARF(28)=SIGMPI(2,2)
  607. *      VARF(29)=SIGMPI(3,3)
  608. *      VARF(30)=SIGMPI(1,3)
  609.  
  610.  2000 CONTINUE
  611.  
  612.         RETURN
  613.         END
  614.  
  615.  
  616.  

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