Télécharger la plaquette Cast3M


fr en

Télécharger aactum.eso

Retour à la liste

Numérotation des lignes : 

aactum

  1. C AACTUM    SOURCE    CB215821  16/04/21    21:15:00     8920
  2.       SUBROUTINE AACTUM(SIGT,EPST,ENDT,XLIMT,SIGD,EPSD,ENDD,XLIMD,
  3.      *                    VDEP,EPSC,EPSPLS,XLMD,PPLUS,
  4.      *                    VDEP1,SIFF,VDEP2,FDDD,DEPS,DSIG,EDED,FD,
  5.      *                    PROD,ENDOPL,EPL,E,EP,D,CONT3,CONT4,EPLUS,
  6.      *                    EDPLUS,U,V,R,W1,W2)
  7. C===================================================================
  8. C
  9. C     CALCUL DES GRANDEURS A LA FIN DU PAS
  10. C
  11. C====================================================================
  12. C
  13. C     CREATION : F.CORMERY
  14. C                E.N.S.M.A - LMPM
  15. C                JUIN 1993
  16. C
  17. C====================================================================
  18.       IMPLICIT INTEGER(I-N)
  19.       IMPLICIT REAL*8 (A-H,O-Z)
  20.       REAL*8 LAMB,MUT
  21. C-----------------------------------------------------------------------
  22. C                           COMMON SPECIFIQUE
  23. C-----------------------------------------------------------------------
  24.       COMMON/APENTE/ LAMB,MUT,G,ALP1,BETA
  25.       COMMON/MAED/  B,C0,CPRIME
  26.       COMMON/DSIGEP/ H
  27. C-----------------------------------------------------------------------
  28. C                           DATA ET DIMENSION
  29. C-----------------------------------------------------------------------
  30. C                                                  N6
  31. CC      DIMENSION  SIGT(6),EPST(6),ENDT(6),SIFF(6),DEPS(6),
  32. C                                          N12       N18
  33. CC     *           SIGD(6),EPSD(6),ENDD(6),DSIG(6),EDED(6),
  34. C                                       N24
  35. CC     *           VDEP(6,6),VDEP1(6,6),FD(6),
  36. C                                    N30        N36        N42
  37. CC     *           EPSC(6),EPSPLS(6),PROD(6),ENDOPL(6),EPL(6)
  38. C                  N51     N60   N69       N78    N87          N96
  39. CC      DIMENSION  E(3,3),EP(3,3),D(3,3),CONT3(3,3),CONT4(3,3),EPLUS(3,3)
  40. C                   N105
  41. CC     *           ,EDPLUS(3,3),VDEP2(6,6)
  42. C                  N186       N267       N348
  43. CC      DIMENSION  U(3,3,3,3),V(3,3,3,3),R(3,3,3,3),PPLUS(3,3,3,3),
  44. C                  N429          N510
  45. CC     *           W1(3,3,3,3),W2(3,3,3,3)
  46. C
  47.       DIMENSION  SIGT(*),EPST(*),ENDT(*),SIFF(*),DEPS(*),
  48.      *           SIGD(*),EPSD(*),ENDD(*),DSIG(*),EDED(*),
  49.      *           VDEP(6,*),VDEP1(6,*),FD(*),
  50.      *           EPSC(*),EPSPLS(*),PROD(*),ENDOPL(*),EPL(*)
  51.       DIMENSION  E(3,*),EP(3,*),D(3,*),CONT3(3,*),CONT4(3,*),EPLUS(3,*)
  52.      *           ,EDPLUS(3,*),VDEP2(6,*)
  53.       DIMENSION  U(3,3,3,*),V(3,3,3,*),R(3,3,3,*),PPLUS(3,3,3,*),
  54.      *           W1(3,3,3,*),W2(3,3,3,*)
  55. C--------------
  56.       DATA ZERO/0.D0/,UN/1.D0/,DEUX/2.D0/
  57.       DATA UD/0.5D0/
  58. C-----------------------------------------------------------------------
  59. C                            CALCUL DES TRACES
  60. C-----------------------------------------------------------------------
  61.        TRAC1=EPSPLS(1)**2+EPSPLS(2)**2+EPSPLS(3)**2
  62.      *      +UD*(EPSPLS(4)**2+EPSPLS(5)**2+EPSPLS(6)**2)
  63. C----------------
  64.        TRAC3=EPSD(1)+EPSD(2)+EPSD(3)
  65. C----------------
  66.        TRAC4=EPSPLS(1)*EPSD(1)+EPSPLS(2)*EPSD(2)+EPSPLS(3)*EPSD(3)
  67.      *      +UD*(EPSPLS(4)*EPSD(4)+EPSPLS(5)*EPSD(5)
  68.      *                      +EPSPLS(6)*EPSD(6))
  69. C----------------
  70.        TRAC5=EPSPLS(1)+EPSPLS(2)+EPSPLS(3)
  71. C************************************************************************
  72. C                   ACTUALISATION DES GRANDEURS
  73. C************************************************************************
  74. C
  75. C------------------------------------------------------------------------
  76. C                     CALCUL DE D n+1
  77. C------------------------------------------------------------------------
  78.        DO 10 I=1,6
  79.           IF(I.LE.3) ENDD(I)=ENDT(I)+(XLMD/(UN-B*XLMD))*
  80.      *           (EPSPLS(I)/SQRT(DEUX*TRAC1)+B*ENDT(I))
  81.           IF(I.GE.4) ENDD(I)=ENDT(I)+(XLMD/(UN-B*XLMD))*
  82.      *           ((EPSPLS(I)/DEUX)/SQRT(DEUX*TRAC1)+B*ENDT(I))
  83.           IF(ABS(ENDD(I)).LE.1E-15) ENDD(I)=0.D0
  84. 10     CONTINUE
  85. C------------------------------------------------------------------------
  86. C                    CALCUL DES FORCES THERMO
  87. C------------------------------------------------------------------------
  88. C---------CALCUL DU PRODUIT EPS.EPS
  89.        EDED(1)=EPSD(1)**2+EPSD(6)**2+EPSD(5)**2
  90.        EDED(2)=EPSD(2)**2+EPSD(6)**2+EPSD(4)**2
  91.        EDED(1)=EPSD(3)**2+EPSD(4)**2+EPSD(5)**2
  92.        EDED(4)=EPSD(6)*EPSD(5)+EPSD(2)*EPSD(4)+EPSD(4)*EPSD(3)
  93.        EDED(5)=EPSD(1)*EPSD(5)+EPSD(6)*EPSD(4)+EPSD(5)*EPSD(3)
  94.        EDED(6)=EPSD(6)*EPSD(1)+EPSD(2)*EPSD(6)+EPSD(4)*EPSD(5)
  95.        DO 11 I=1,6
  96.          FD(I)=-G*EPSD(I)-ALP1*TRAC3*EPSD(I)-2*BETA*EDED(I)
  97. 11     CONTINUE
  98. C------------------------------------------------------------------------
  99. C                    CALCUL DE FD:DD
  100. C------------------------------------------------------------------------
  101.         FDDD=ZERO
  102.         DO 12 I=1,6
  103.           FDDD=FDDD+FD(I)*(ENDD(I)-ENDT(I))
  104. 12      CONTINUE
  105. C------------------------------------------------------------------------
  106. C                    CALCUL DES TRACES
  107. C------------------------------------------------------------------------
  108.        TRAC2=EPSPLS(1)*ENDD(1)+EPSPLS(2)*ENDD(2)+EPSPLS(3)*ENDD(3)
  109.      *      +(EPSPLS(4)*ENDD(4)+EPSPLS(5)*ENDD(5)
  110.      *                      +EPSPLS(6)*ENDD(6))
  111.        TRAC6=ENDD(1)+ENDD(2)+ENDD(3)
  112. C------------------------------------------------------------------------
  113. C                   CALCUL DE   C0+C1 tr(Dn+1)
  114. C------------------------------------------------------------------------
  115.        XLIMD=C0+CPRIME*(ENDD(1)+ENDD(2)+ENDD(3))
  116. C------------------------------------------------------------------------
  117. C                   CALCUL DE D+ n+1
  118. C------------------------------------------------------------------------
  119.        MOD1=0
  120.        IF(MOD1.EQ.1)THEN
  121.        CALL T4CT2(PPLUS,ENDD,ENDOPL,6)
  122.        GOTO 61
  123.        ENDIF
  124.        DO 62 I=1,6
  125.        ENDOPL(I)=ENDD(I)
  126. 62     CONTINUE
  127. 61     CONTINUE
  128. C------------------------------------------------------------------------
  129. C
  130. C-------------------------------------------------------------------------
  131.        DO 13 I=4,6
  132.        EPSPLS(I)=EPSPLS(I)/2
  133. 13     CONTINUE
  134.        CALL T4CT2(PPLUS,EPSPLS,EPL,6)
  135.        DO 14 I=4,6
  136.        EPSPLS(I)=EPSPLS(I)*2
  137. 14     CONTINUE
  138. C-------------------------------------------------------------------------
  139. C                    CALCUL DE SIG n+1
  140. C-------------------------------------------------------------------------
  141.        TRACE2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)+
  142.      *        EPSD(4)*ENDD(4)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  143. C-------
  144.        PROD(1)=DEUX*(EPSD(1)*ENDD(1))+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  145.        PROD(2)=DEUX*(EPSD(2)*ENDD(2))+EPSD(4)*ENDD(4)+EPSD(6)*ENDD(6)
  146.        PROD(3)=DEUX*(EPSD(3)*ENDD(3))+EPSD(5)*ENDD(5)+EPSD(4)*ENDD(4)
  147.        PROD(4)=UD*EPSD(6)*ENDD(5)+EPSD(2)*ENDD(4)+UD*EPSD(4)*ENDD(3)+
  148.      *      UD*EPSD(5)*ENDD(6)+UD*EPSD(4)*ENDD(2)+EPSD(3)*ENDD(4)
  149.        PROD(5)=UD*EPSD(5)*ENDD(1)+UD*EPSD(4)*ENDD(6)+EPSD(3)*ENDD(5)+
  150.      *         EPSD(1)*ENDD(5)+UD*EPSD(6)*ENDD(4)+UD*EPSD(5)*ENDD(3)
  151.        PROD(6)=UD*EPSD(6)*ENDD(1)+EPSD(2)*ENDD(6)+UD*EPSD(4)*ENDD(5)+
  152.      *       EPSD(1)*ENDD(6)+UD*EPSD(6)*ENDD(2)+UD*EPSD(5)*ENDD(4)
  153. C-------
  154.         NP=1
  155.         IF(NP.EQ.1)THEN
  156.         DO 30 I=1,3
  157.           SIGD(I)=G*ENDOPL(I)+LAMB*TRAC3+ALP1*(TRACE2+TRAC3*ENDD(I))
  158.      *            +DEUX*BETA*PROD(I)+DEUX*MUT*EPSD(I)
  159. 30      CONTINUE
  160.         DO 35 I=4,6
  161.           SIGD(I)=G*ENDOPL(I)+ALP1*TRAC3*ENDD(I)+MUT*EPSD(I)
  162.      *            +DEUX*BETA*PROD(I)
  163. 35      CONTINUE
  164.         GOTO 21
  165.         ENDIF
  166. C------------------------------------------------------------------
  167. C                CALCUL MATRICE TANGENTE continue
  168. C                    ET MATRICE TANGENTE consistante
  169. C------------------------------------------------------------------
  170. C
  171. C----------------------------------Calcul des coef.----------------
  172. 21     XA=DEUX*TRAC1
  173.        XM=B*ABS(G)*SQRT(TRAC1/DEUX)+(B**2)*ABS(G)*TRAC2
  174.      *     -CPRIME*TRAC5/SQRT(DEUX*TRAC1)-CPRIME*B*TRAC6
  175.        V2=-XLMD/(SQRT(XA)*TRAC5)
  176.        V3=XLMD/SQRT(XA)
  177.        V5=ABS(G)/(SQRT(XA)*CPRIME*TRAC5)
  178. C---------------------------------Passage en notation indicielle---
  179.        DO 42 I=1,3
  180.           D(I,I)=ENDD(I)
  181.           EP(I,I)=EPSPLS(I)
  182.           E(I,I)=EPSD(I)
  183.           EDPLUS(I,I)=ENDOPL(I)
  184.           EPLUS(I,I)=EPL(I)
  185.           DO 43 J=I+1,3
  186.              IF ((I.EQ.1).AND.(J.EQ.2)) K=6
  187.              IF ((I.EQ.1).AND.(J.EQ.3)) K=5
  188.              IF ((I.EQ.2).AND.(J.EQ.3)) K=4
  189.              D(I,J)=ENDD(K)
  190.              D(J,I)=ENDD(K)
  191.              E(I,J)=UD*EPSD(K)
  192.              E(J,I)=UD*EPSD(K)
  193.              EP(I,J)=UD*EPSPLS(K)
  194.              EP(J,I)=UD*EPSPLS(K)
  195.              EDPLUS(I,J)=ENDOPL(K)
  196.              EDPLUS(J,I)=ENDOPL(K)
  197.              EPLUS(I,J)=EPL(K)
  198.              EPLUS(J,I)=EPL(K)
  199. 43        CONTINUE
  200. 42     CONTINUE
  201.  
  202. C--------------------- Calcul des produit contractes -----------------
  203.        CONT1=ZERO
  204.        CONT2=ZERO
  205.        DO 44 I=1,3
  206.           DO 45 J=1,3
  207.              CONT1=CONT1+E(I,J)*EP(I,J)
  208.              CONT2=CONT2+E(I,J)*D(I,J)
  209.              CONT3(I,J)=ZERO
  210.              CONT4(I,J)=ZERO
  211.              DO 46 K=1,3
  212.              CONT3(I,J)=CONT3(I,J)+E(I,K)*EP(K,J)+EP(I,K)*E(K,J)
  213.              CONT4(I,J)=CONT4(I,J)+E(I,K)*D(K,J)+D(I,K)*E(K,J)
  214. 46           CONTINUE
  215. 45        CONTINUE
  216. 44     CONTINUE
  217.  
  218. C----------------------------Calcul de la matrice tangente---------
  219.        DO 50 I=1,3
  220.           DO 60 J=1,3
  221.              DO 70 K=1,3
  222.                 DO 80 L=1,3
  223. C-------------------------
  224.          U(I,J,K,L)=LAMB*DELT(I,J)*DELT(K,L)+MUT*(DELT(I,K)*DELT(J,L)
  225.      * +DELT(I,L)*DELT(J,K))+ALP1*(DELT(I,J)*D(K,L)+D(I,J)*DELT(K,L))
  226.      *    +BETA*(D(I,K)*DELT(J,L)+D(L,J)*DELT(I,K)+DELT(I,L)
  227.      * *D(J,K)+D(I,L)*DELT(J,K))
  228. C-------------------------
  229.          V(I,J,K,L)=(1/XM)*(ALP1*G*CONT1*DELT(I,J)/XA+ALP1*B*G
  230.      *     *CONT2*DELT(I,J)/SQRT(XA)+(ALP1*G/XA)*EP(I,J)
  231.      *     *TRAC3+(ALP1*B*G/SQRT(XA))*D(I,J)*TRAC3+(DEUX*BETA*G/XA)
  232.      *    *CONT3(I,J)+DEUX*BETA*(B*G/SQRT(XA))*CONT4(I,J)
  233.      *     +(G**2/XA)*EPLUS(I,J)+(B*(G**2)/SQRT(XA))
  234.      *     *EDPLUS(I,J))*EPLUS(K,L)
  235. C-------------------------
  236.          R(I,J,K,L)=(1/XM)*((ALP1*B*G/SQRT(XA))*CONT1*DELT(I,J)
  237.      *     +ALP1*(B**2)*G*CONT2*DELT(I,J)+(ALP1*B*G/SQRT(XA))
  238.      *     *EP(I,J)*TRAC3+ALP1*(B**2)*G*D(I,J)*TRAC3+DEUX*BETA
  239.      *     *(B*G/SQRT(XA))*CONT3(I,J)+DEUX*BETA*(B**2)*G*
  240.      *     CONT4(I,J)+(B*(G**2)/SQRT(XA))*EPLUS(I,J)
  241.      *     +B**2*G**2*EDPLUS(I,J))*EDPLUS(K,L)
  242. C-------------------------
  243.          W1(I,J,K,L)=0.D0
  244.          W2(I,J,K,L)=0.D0
  245.       DO 81 M=1,3
  246.       DO 81 N=1,3
  247.          W1(I,J,K,L)=(G*V5*EPLUS(I,J)*E(M,N)+G*V2*EPLUS(I,J)*DELT(M,N)
  248.      *     +G*V3*PPLUS(I,J,M,N)+DEUX*BETA*V5*CONT3(I,J)*E(M,N)
  249.      *     +DEUX*BETA*V2*CONT3(I,J)*DELT(M,N)+BETA*V3*
  250.      *     (E(I,M)*DELT(J,N)+DELT(J,M)*E(I,N)+DELT(I,M)*E(J,N)+
  251.      *     E(J,M)*DELT(I,N))+ALP1*V5*TRAC4*DELT(I,J)*EP(M,N)
  252.      *     +ALP1*V2*TRAC4*DELT(I,J)*DELT(M,N)+ALP1*V3*
  253.      *     DELT(I,J)*E(M,N)+ALP1*V5*TRAC3*EP(I,J)*EP(M,N)
  254.      *     )*PPLUS(M,N,K,L)+W1(I,J,K,L)
  255.          W2(I,J,K,L)=(ALP1*V2*TRAC3*EP(I,J)*DELT(M,N)
  256.      *   +ALP1*V3*TRAC3*(DELT(I,M)*DELT(N,J)+DELT(J,M)*DELT(N,I))
  257.      *    /2.D0)*PPLUS(M,N,K,L)+W2(I,J,K,L)
  258. 81    CONTINUE
  259. C-------------------------
  260. 80              CONTINUE
  261. 70           CONTINUE
  262. 60        CONTINUE
  263. 50     CONTINUE
  264. C--------------------------
  265.        DO 110 I=1,6
  266.           DO 120 J=1,6
  267.           IF (I.EQ.1) THEN
  268.                        M=1
  269.                        N=1
  270.           ENDIF
  271.           IF (I.EQ.2) THEN
  272.                        M=2
  273.                        N=2
  274.           ENDIF
  275.           IF (I.EQ.3) THEN
  276.                        M=3
  277.                        N=3
  278.           ENDIF
  279.           IF (I.EQ.4) THEN
  280.                        M=1
  281.                        N=2
  282.           ENDIF
  283.           IF (I.EQ.5) THEN
  284.                        M=2
  285.                        N=3
  286.           ENDIF
  287.           IF (I.EQ.6) THEN
  288.                        M=1
  289.                        N=3
  290.           ENDIF
  291.           IF (J.EQ.1) THEN
  292.                        K=1
  293.                        L=1
  294.           ENDIF
  295.           IF (J.EQ.2) THEN
  296.                        K=2
  297.                        L=2
  298.           ENDIF
  299.           IF (J.EQ.3) THEN
  300.                        K=3
  301.                        L=3
  302.           ENDIF
  303.           IF (J.EQ.4) THEN
  304.                        K=1
  305.                        L=2
  306.           ENDIF
  307.           IF (J.EQ.5) THEN
  308.                        K=2
  309.                        L=3
  310.           ENDIF
  311.           IF (J.EQ.6) THEN
  312.                        K=1
  313.                        L=3
  314.           ENDIF
  315.           VDEP2(I,J)=U(M,N,K,L)
  316.           VDEP(I,J)=U(M,N,K,L)+V(M,N,K,L)+R(M,N,K,L)
  317.           VDEP1(I,J)=U(M,N,K,L)+W1(M,N,K,L)+W2(M,N,K,L)
  318.           IF(ABS(VDEP2(I,J)).LE.1E-15)VDEP2(I,J)=0.D0
  319.           IF(ABS(VDEP1(I,J)).LE.1E-15)VDEP1(I,J)=0.D0
  320.           IF(ABS(VDEP(I,J)).LE.1E-15)VDEP(I,J)=0.D0
  321. 120       CONTINUE
  322. 110   CONTINUE
  323. C---------------------------------------------------------------------
  324. C                     CALCUL DE SIG n+1
  325. C---------------------------------------------------------------------
  326.        DO 111 I=1,6
  327.           DO 121 J=1,6
  328.           IF (I.EQ.1) THEN
  329.                        M=1
  330.                        N=1
  331.           ENDIF
  332.           IF (I.EQ.2) THEN
  333.                        M=2
  334.                        N=2
  335.           ENDIF
  336.           IF (I.EQ.3) THEN
  337.                        M=3
  338.                        N=3
  339.           ENDIF
  340.           IF (I.EQ.4) THEN
  341.                        M=2
  342.                        N=3
  343.           ENDIF
  344.           IF (I.EQ.5) THEN
  345.                        M=1
  346.                        N=3
  347.           ENDIF
  348.           IF (I.EQ.6) THEN
  349.                        M=1
  350.                        N=2
  351.           ENDIF
  352.           IF (J.EQ.1) THEN
  353.                        K=1
  354.                        L=1
  355.           ENDIF
  356.           IF (J.EQ.2) THEN
  357.                        K=2
  358.                        L=2
  359.           ENDIF
  360.           IF (J.EQ.3) THEN
  361.                        K=3
  362.                        L=3
  363.           ENDIF
  364.           IF (J.EQ.4) THEN
  365.                        K=2
  366.                        L=3
  367.           ENDIF
  368.           IF (J.EQ.5) THEN
  369.                        K=1
  370.                        L=3
  371.           ENDIF
  372.           IF (J.EQ.6) THEN
  373.                        K=1
  374.                        L=2
  375.           ENDIF
  376.           VDEP2(I,J)=U(M,N,K,L)
  377. 121       CONTINUE
  378. 111    CONTINUE
  379. C-------------------------------------------------------------
  380.        NP=2
  381.        IF(NP.EQ.2)THEN
  382.        DO 302 I=1,6
  383.        SIFF(I)=0.D0
  384.        DO 301 J=1,6
  385.         SIFF(I)=VDEP2(I,J)*EPSD(J)+SIFF(I)
  386. 301    CONTINUE
  387.         SIFF(I)=G*ENDOPL(I)+SIFF(I)
  388. 302    CONTINUE
  389.        ENDIF
  390. C-------------------------------------------------------------
  391.        H=0
  392.        DO 307 I=1,6
  393.        DEPS(I)=EPSD(I)-EPST(I)
  394.        DSIG(I)=SIGD(I)-SIGT(I)
  395.        H=H+DEPS(I)*DSIG(I)
  396. 307    CONTINUE
  397. C       WRITE(10,1203)H
  398. C1203   FORMAT(/5X,'DSIG : DEPS = ',E12.5)
  399.        RETURN
  400.        END
  401.  
  402.  
  403.  
  404.  
  405.  
  406.  

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