Télécharger la plaquette Cast3M


fr en

Télécharger ared3m.eso

Retour à la liste

Numérotation des lignes : 

ared3m

  1. C ARED3M    SOURCE    CHAT      05/01/12    21:23:28     5004
  2.        SUBROUTINE ARED3M(SIGT,EPST,ENDT,XLIMT,
  3.      *                   SIGD,EPSD,ENDD,XLIMD,
  4.      *                   DEPS,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1,
  5.      *                   VAL2,VAL3,VP2,VP3,EPSTT,FDDD,D,SIFF,VDEII,
  6.      *                   VDEP2,EPSC,PPLUS,U,PROD,ENDOPL,VDEI,KM,AAA)
  7. C======================================================================
  8. C     TEST D ENDOMMAGEMENT INITIAL POUR LES ELEMENTS TRIDIMENSIONNELS
  9. C======================================================================
  10. C
  11. C
  12. C     MODIF.   : F.CORMERY
  13. C                E.N.S.M.A - LMPM
  14. C                JUIN 1993
  15. C
  16. C======================================================================
  17. C     ENTREES : SIGT,EPST,ENDT,    VALEURS AU PAS PRECEDENT
  18. C               DSIG,DEPS          ACCROISSEMENTS INITIAUX
  19. C     SORTIES   SIGD,EPSD,ENDD,    VALEURS A LA FIN DU PAS
  20. C               DEPS               ACCROISSEMENTS
  21. C               VDEP               MATRICE TANGENTE
  22. C======================================================================
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25.       REAL*8 LAMB,MUT
  26. C----------------------------------------------------------------------
  27. C                                                                COMMON
  28. C----------------------------------------------------------------------
  29.       COMMON /APENTE/ LAMB,MUT,G,ALP1,BETA
  30. C----------------------------------------------------------------------
  31. C                                                    DATA ET DIMENSIONS
  32. C----------------------------------------------------------------------
  33. C                                        N9
  34. CC      DIMENSION SIGD(6),EPSD(6),ENDD(6),D(3,3),VAL1(3,3),VP1(3),
  35. C                                        N15        N51
  36. CC     *          SIGT(6),EPST(6),ENDT(6),SIFF(6),VDEII(6,6),
  37. C
  38. CC     *          DEPS(6),VDEP(6,6),VDEP1(6,6),VAL(3,3),VP(3),
  39. C                                                   N87
  40. CC     *          VAL2(3,3),VAL3(3,3),VP3(3),VP2(3),VDEP2(6,6)
  41. C                 N93              N174           N255
  42. CC      DIMENSION EPSC(6),EPSPLS(6),PPLUS(3,3,3,3),U(3,3,3,3)
  43. C                 N261   N267      N303      N309
  44. CC      DIMENSION PROD(6),ENDOPL(6),VDEI(6,6),KM(6),EPSTT(6)
  45. C
  46. C
  47.       DIMENSION AAA(*)
  48.       DIMENSION SIGD(*),EPSD(*),ENDD(*),D(3,*),VAL1(3,*),VP1(*),
  49. C                                        N15        N51
  50.      *          SIGT(*),EPST(*),ENDT(*),SIFF(*),VDEII(6,*),
  51. C
  52.      *          DEPS(*),VDEP(6,*),VDEP1(6,*),VAL(3,*),VP(*),
  53. C                                                   N87
  54.      *          VAL2(3,*),VAL3(3,*),VP3(*),VP2(*),VDEP2(6,*)
  55. C                 N93              N174           N255
  56.       DIMENSION EPSC(*),EPSPLS(*),PPLUS(3,3,3,*),U(3,3,3,*)
  57. C                 N261   N267      N303      N309
  58.       DIMENSION PROD(*),ENDOPL(*),VDEI(6,*),KM(*),EPSTT(*)
  59.       DIMENSION NII(6),MII(6),NIII(6),MIII(6)
  60.       DATA NII/1,2,3,3,3,2/
  61.       DATA MII/1,2,3,2,1,1/
  62.       DATA NIII/1,2,3,2,3,3/
  63.       DATA MIII/1,2,3,1,2,1/
  64.       DATA ZERO/0.D0/
  65.       DATA IMA/6/,DEUX/2.D0/,UD/0.5D0/
  66. C----------------------------------------------------------------------
  67.       MT=6
  68.       IKK=1
  69. C----------------------------------------------------------------------
  70. C                   TENSEUR DE DEFORMATION ACTUALISEE
  71. C----------------------------------------------------------------------
  72.       MMI=1
  73.       GOTO 121
  74. C------------
  75. 3     CONTINUE
  76. C------------
  77.        DO 144 I=1,6
  78.           M= MII(I)
  79.           N=NII(I)
  80.           DO 144 J=1,6
  81.             K=MII(J)
  82.             L=NII(J)
  83.             VDEI(I,J)=U(M,N,K,L)
  84.             VDEII(I,J)=U(M,N,K,L)
  85. 144       CONTINUE
  86. C---------------------------Calcul de la deformation initiale----------
  87.       IF(NPA.EQ.1)THEN
  88. C------------
  89.           CALL INVERE(VDEI,6,6,KM,DET)
  90. C------------
  91.       DO 4 I=1,6
  92.           EPST(I)=0.D0
  93.       DO 4 J=1,6
  94.           EPST(I)=EPST(I)+VDEI(I,J)*(-G*ENDT(J))
  95. 4     CONTINUE
  96. C------------
  97.           MMI=0
  98.       ENDIF
  99. C----------------------------------------------------------------------
  100. 234   CONTINUE
  101.       DO 10 I=1,IMA
  102.          EPSD(I)=EPST(I)+DEPS(I)
  103.          EPSC(I)=EPSD(I)
  104. 10    CONTINUE
  105. C----------------------------------------------------------------------
  106. C               CALCUL DU TENSEUR DE DEFORMATION POSITIVE
  107. C----------------------------------------------------------------------
  108. 7        CALL EPPLUM(EPSC,PPLUS,1,VAL1,VP1,AAA(442),
  109.      .               AAA(478),AAA(514),AAA(517),AAA(526),
  110.      .               AAA(532))
  111. C----------------
  112.       EPSD(4)=EPSD(4)/2
  113.       EPSD(5)=EPSD(5)/2
  114.       EPSD(6)=EPSD(6)/2
  115. C----------------
  116.       CALL T4CT2(PPLUS,EPSD,EPSPLS,6)
  117. C----------------
  118.       DO 324 I=1,6
  119.       IF(ABS(EPSPLS(I)).LE.1E-15)EPSPLS(I)=0.D0
  120. 324   CONTINUE
  121. C----------------
  122.       EPSD(4)=EPSD(4)*2
  123.       EPSD(5)=EPSD(5)*2
  124.       EPSD(6)=EPSD(6)*2
  125.       EPSPLS(4)=EPSPLS(4)*2
  126.       EPSPLS(5)=EPSPLS(5)*2
  127.       EPSPLS(6)=EPSPLS(6)*2
  128. C----------------------------------------------------------------------
  129. C                          TEST D ENDOMMAGEMENT
  130.  
  131. C----------------------------------------------------------------------
  132.       CALL AVENDO(EPSPLS,ENDT,DEQ)
  133.       CRIT=DEQ-XLIMT
  134.       IF(CRIT.LE.ZERO)GO TO 100
  135. C----------------------------------------------------------------------
  136. C                        CALCUL DU MULTIPLICATEUR
  137. C----------------------------------------------------------------------
  138.       XLMD=ZERO
  139.       CALL AMULTM(EPSPLS,ENDT,XLMD,XLIMT)
  140. C----------------------------------------------------------------------
  141. C                      ACTUALISATION DES GRANDEURS
  142. C----------------------------------------------------------------------
  143.       CALL AACTUM(SIGT,EPST,ENDT,XLIMT,SIGD,EPSD,ENDD,XLIMD,
  144.      *              VDEP,EPSC,EPSPLS,XLMD,PPLUS,
  145.      *              VDEP1,SIFF,VDEP2,FDDD,AAA(442),AAA(448),
  146.      *              AAA(454),AAA(460),AAA(466),AAA(472),
  147.      *              AAA(478),AAA(484),AAA(493),AAA(502),
  148.      *              AAA(511),AAA(520),AAA(529),AAA(538),
  149.      *              AAA(547),AAA(628),AAA(709),AAA(790),
  150.      *              AAA(871))
  151.       ISS=2
  152.       GO TO 899
  153. C**********************************************************************
  154. C                            CAS ELASTIQUE
  155. C**********************************************************************
  156. 100   CONTINUE
  157.          ISS=1
  158. C----------------------------------------------------------------------
  159.       DO 110 I=1,IMA
  160.          ENDD(I)=ENDT(I)
  161. 110   CONTINUE
  162.          XLIMD=XLIMT
  163. C--------------------------------CALCUL   PPLUS:D=D+  ------------
  164. *      MODEE1=1
  165. *      IF(MODEE1.EQ.1)THEN
  166.          DO 123 I=1,6
  167. 123      ENDOPL(I)=ENDD(I)
  168. *        GOTO 108
  169. *        ENDIF
  170. *         CALL T4CT2(PPLUS,ENDD,ENDOPL,6)
  171. C--------------------------------CALCUL DES CONTRAINTES----------------
  172. *108   TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)
  173.       TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)
  174.      *     +EPSD(4)*ENDD(4)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  175.       TRAC3=EPSD(1)+EPSD(2)+EPSD(3)
  176.       PROD(1)=DEUX*EPSD(1)*ENDD(1)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  177.       PROD(2)=DEUX*EPSD(2)*ENDD(2)+EPSD(4)*ENDD(4)+EPSD(6)*ENDD(6)
  178.       PROD(3)=DEUX*EPSD(3)*ENDD(3)+EPSD(5)*ENDD(5)+EPSD(4)*ENDD(4)
  179.       PROD(4)=UD*EPSD(6)*ENDD(5)+EPSD(2)*ENDD(4)+UD*EPSD(4)*ENDD(3)+
  180.      *       UD*EPSD(5)*ENDD(6)+UD*EPSD(4)*ENDD(2)+EPSD(3)*ENDD(4)
  181.       PROD(5)=UD*EPSD(5)*ENDD(1)+UD*EPSD(4)*ENDD(6)+EPSD(3)*ENDD(5)+
  182.      *       EPSD(1)*ENDD(5)+UD*EPSD(6)*ENDD(4)+UD*EPSD(5)*ENDD(3)
  183.       PROD(6)=UD*EPSD(6)*ENDD(1)+EPSD(2)*ENDD(6)+UD*EPSD(4)*ENDD(5)+
  184.      *       EPSD(1)*ENDD(6)+UD*EPSD(6)*ENDD(2)+UD*EPSD(5)*ENDD(4)
  185. C---------------------
  186.       DO 115 I=1,3
  187.          SIGD(I)=G*ENDOPL(I)+LAMB*TRAC3+ALP1*(TRAC2+TRAC3*ENDD(I))
  188.      *           +DEUX*BETA*PROD(I)+DEUX*MUT*EPSD(I)
  189. 115   CONTINUE
  190. C---------------------
  191.       DO 116 I=4,6
  192.          SIGD(I)=G*ENDOPL(I)+ALP1*TRAC3*ENDD(I)
  193.      *           +DEUX*BETA*PROD(I)+MUT*EPSD(I)
  194. 116   CONTINUE
  195. C------------------
  196.       DO 17 I=1,6
  197.          SIFF(I)=0
  198.       DO 17 J=1,6
  199.       SIFF(I)=SIFF(I)+VDEII(I,J)*EPSC(J)
  200. 17    CONTINUE
  201.       DO 18 I=1,6
  202.       SIFF(I)=SIFF(I)+G*ENDOPL(I)
  203. 18    CONTINUE
  204. C----------------------------------------------------------------------
  205.       MMI=0
  206. C----------------------------------------------------------------------
  207. C            CALCUL MATRICE TANGENTE ELASTIQUE ANISOTROPE
  208. C----------------------------------------------------------------------
  209. C
  210. C---------------------------------Passage en notation indicielle-------
  211. 121    CONTINUE
  212.           DO 42 I=1,3
  213.           D(I,I)=ENDT(I)
  214.           DO 43 J=I+1,3
  215.              K = 9 - I -J
  216. *             IF ((I.EQ.1).AND.(J.EQ.2)) K=6
  217. *             IF ((I.EQ.1).AND.(J.EQ.3)) K=5
  218. *             IF ((I.EQ.2).AND.(J.EQ.3)) K=4
  219.              D(I,J)=ENDT(K)
  220.              D(J,I)=ENDT(K)
  221. 43        CONTINUE
  222. 42     CONTINUE
  223.  
  224. C----------------------------Calcul de la matrice tangente-------------
  225.        DO 51 I=1,3
  226.           DO 51 J=1,3
  227.              DO 51 K=1,3
  228.                 DO 51 L=1,3
  229. C-------------------------
  230.          U(I,J,K,L)=LAMB*DELT(I,J)*DELT(K,L)+MUT*(DELT(I,K)*DELT(J,L)
  231.      * +DELT(I,L)*DELT(J,K))+ALP1*(DELT(I,J)*D(K,L)+D(I,J)*DELT(K,L))
  232.      *    +BETA*(D(I,K)*DELT(J,L)+D(L,J)*DELT(I,K)+DELT(I,L)
  233.      * *D(J,K)+D(I,L)*DELT(J,K))
  234. C-------------------------
  235. 51     CONTINUE
  236. C-------------------------------Notation de VOIGT--------------------
  237.        IF(MMI.EQ.1)GOTO 3
  238.        DO 134 I=1,6
  239.           M=MIII(I)
  240.           N=NIII(I)
  241.           DO 134 J=1,6
  242.             K= MIII(J)
  243.             L= NIII(J)
  244.             VDEP(I,J)=U(M,N,K,L)
  245. 134       CONTINUE
  246. C-------------------------------------------------------------------
  247.         DO 142 I=1,6
  248.         DO 142 J=1,6
  249.              VDEP1(I,J)=VDEP(I,J)
  250.              VDEP2(I,J)=VDEP(I,J)
  251. 142     CONTINUE
  252. C********************************************************************
  253. C
  254. C********************************************************************
  255. 899   DO 321 I=1,6
  256.         IF(ABS(SIGD(I)).LE.1E-15)SIGD(I)=0.D0
  257. 321   CONTINUE
  258. C--------------------------------------------------------------------
  259.       CALL VECT(SIGD,VAL,VP,AAA(442),AAA(451),AAA(457))
  260.       DO 332 I=4,6
  261. 332   EPSPLS(I)=EPSPLS(I)/2
  262.       CALL VECT(EPSPLS,VAL2,VP2,AAA(442),AAA(451),AAA(457))
  263.       DO 333 I=4,6
  264. 333   EPSPLS(I)=EPSPLS(I)*2
  265.       CALL VECT(ENDD,VAL3,VP3,AAA(442),AAA(451),AAA(457))
  266. C--------------------------------------------------------------------
  267.           CALL INVERE(VDEP2,6,6,KM,DET)
  268. C------------
  269.       DO 401 I=1,6
  270.           EPSTT(I)=0.D0
  271.       DO 401 J=1,6
  272.           EPSTT(I)=EPSTT(I)+VDEP2(I,J)*(-G*ENDD(J))
  273. 401     CONTINUE
  274.        RETURN
  275.        END
  276.  
  277.  
  278.  
  279.  

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