Télécharger la plaquette Cast3M


fr en

Télécharger ric2nl.eso

Retour à la liste

Numérotation des lignes : 

ric2nl

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

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