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    CB215821  26/08/24    21:15:10     12622          
  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 900 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.   900  CONTINUE
  87. C---------------------------Calcul de la deformation initiale----------
  88.       IF(NPA.EQ.1)THEN
  89. C------------
  90.           CALL INVERE(VDEI,6,6,KM,DET)
  91. C------------
  92.       DO 901 I=1,6
  93.           EPST(I)=0.D0
  94.       DO 4 J=1,6
  95.           EPST(I)=EPST(I)+VDEI(I,J)*(-G*ENDT(J))
  96. 4     CONTINUE
  97.   901 CONTINUE
  98. C------------
  99.           MMI=0
  100.       ENDIF
  101. C----------------------------------------------------------------------
  102. 234   CONTINUE
  103.       DO 10 I=1,IMA
  104.          EPSD(I)=EPST(I)+DEPS(I)
  105.          EPSC(I)=EPSD(I)
  106. 10    CONTINUE
  107. C----------------------------------------------------------------------
  108. C               CALCUL DU TENSEUR DE DEFORMATION POSITIVE
  109. C----------------------------------------------------------------------
  110. 7        CALL EPPLUM(EPSC,PPLUS,1,VAL1,VP1,AAA(442),
  111.      .               AAA(478),AAA(514),AAA(517),AAA(526),
  112.      .               AAA(532))
  113. C----------------
  114.       EPSD(4)=EPSD(4)/2
  115.       EPSD(5)=EPSD(5)/2
  116.       EPSD(6)=EPSD(6)/2
  117. C----------------
  118.       CALL T4CT2(PPLUS,EPSD,EPSPLS,6)
  119. C----------------
  120.       DO 324 I=1,6
  121.       IF(ABS(EPSPLS(I)).LE.1E-15)EPSPLS(I)=0.D0
  122. 324   CONTINUE
  123. C----------------
  124.       EPSD(4)=EPSD(4)*2
  125.       EPSD(5)=EPSD(5)*2
  126.       EPSD(6)=EPSD(6)*2
  127.       EPSPLS(4)=EPSPLS(4)*2
  128.       EPSPLS(5)=EPSPLS(5)*2
  129.       EPSPLS(6)=EPSPLS(6)*2
  130. C----------------------------------------------------------------------
  131. C                          TEST D ENDOMMAGEMENT
  132.  
  133. C----------------------------------------------------------------------
  134.       CALL AVENDO(EPSPLS,ENDT,DEQ)
  135.       CRIT=DEQ-XLIMT
  136.       IF(CRIT.LE.ZERO)GO TO 100
  137. C----------------------------------------------------------------------
  138. C                        CALCUL DU MULTIPLICATEUR
  139. C----------------------------------------------------------------------
  140.       XLMD=ZERO
  141.       CALL AMULTM(EPSPLS,ENDT,XLMD,XLIMT)
  142. C----------------------------------------------------------------------
  143. C                      ACTUALISATION DES GRANDEURS
  144. C----------------------------------------------------------------------
  145.       CALL AACTUM(SIGT,EPST,ENDT,XLIMT,SIGD,EPSD,ENDD,XLIMD,
  146.      *              VDEP,EPSC,EPSPLS,XLMD,PPLUS,
  147.      *              VDEP1,SIFF,VDEP2,FDDD,AAA(442),AAA(448),
  148.      *              AAA(454),AAA(460),AAA(466),AAA(472),
  149.      *              AAA(478),AAA(484),AAA(493),AAA(502),
  150.      *              AAA(511),AAA(520),AAA(529),AAA(538),
  151.      *              AAA(547),AAA(628),AAA(709),AAA(790),
  152.      *              AAA(871))
  153.       ISS=2
  154.       GO TO 899
  155. C**********************************************************************
  156. C                            CAS ELASTIQUE
  157. C**********************************************************************
  158. 100   CONTINUE
  159.          ISS=1
  160. C----------------------------------------------------------------------
  161.       DO 110 I=1,IMA
  162.          ENDD(I)=ENDT(I)
  163. 110   CONTINUE
  164.          XLIMD=XLIMT
  165. C--------------------------------CALCUL   PPLUS:D=D+  ------------
  166. *      MODEE1=1
  167. *      IF(MODEE1.EQ.1)THEN
  168.          DO 123 I=1,6
  169.          ENDOPL(I)=ENDD(I)
  170.   123    CONTINUE
  171. *        GOTO 108
  172. *        ENDIF
  173. *         CALL T4CT2(PPLUS,ENDD,ENDOPL,6)
  174. C--------------------------------CALCUL DES CONTRAINTES----------------
  175. *108   TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)
  176.       TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)
  177.      *     +EPSD(4)*ENDD(4)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  178.       TRAC3=EPSD(1)+EPSD(2)+EPSD(3)
  179.       PROD(1)=DEUX*EPSD(1)*ENDD(1)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  180.       PROD(2)=DEUX*EPSD(2)*ENDD(2)+EPSD(4)*ENDD(4)+EPSD(6)*ENDD(6)
  181.       PROD(3)=DEUX*EPSD(3)*ENDD(3)+EPSD(5)*ENDD(5)+EPSD(4)*ENDD(4)
  182.       PROD(4)=UD*EPSD(6)*ENDD(5)+EPSD(2)*ENDD(4)+UD*EPSD(4)*ENDD(3)+
  183.      *       UD*EPSD(5)*ENDD(6)+UD*EPSD(4)*ENDD(2)+EPSD(3)*ENDD(4)
  184.       PROD(5)=UD*EPSD(5)*ENDD(1)+UD*EPSD(4)*ENDD(6)+EPSD(3)*ENDD(5)+
  185.      *       EPSD(1)*ENDD(5)+UD*EPSD(6)*ENDD(4)+UD*EPSD(5)*ENDD(3)
  186.       PROD(6)=UD*EPSD(6)*ENDD(1)+EPSD(2)*ENDD(6)+UD*EPSD(4)*ENDD(5)+
  187.      *       EPSD(1)*ENDD(6)+UD*EPSD(6)*ENDD(2)+UD*EPSD(5)*ENDD(4)
  188. C---------------------
  189.       DO 115 I=1,3
  190.          SIGD(I)=G*ENDOPL(I)+LAMB*TRAC3+ALP1*(TRAC2+TRAC3*ENDD(I))
  191.      *           +DEUX*BETA*PROD(I)+DEUX*MUT*EPSD(I)
  192. 115   CONTINUE
  193. C---------------------
  194.       DO 116 I=4,6
  195.          SIGD(I)=G*ENDOPL(I)+ALP1*TRAC3*ENDD(I)
  196.      *           +DEUX*BETA*PROD(I)+MUT*EPSD(I)
  197. 116   CONTINUE
  198. C------------------
  199.       DO 902 I=1,6
  200.          SIFF(I)=0
  201.       DO 17 J=1,6
  202.       SIFF(I)=SIFF(I)+VDEII(I,J)*EPSC(J)
  203. 17    CONTINUE
  204.   902 CONTINUE
  205.       DO 18 I=1,6
  206.       SIFF(I)=SIFF(I)+G*ENDOPL(I)
  207. 18    CONTINUE
  208. C----------------------------------------------------------------------
  209.       MMI=0
  210. C----------------------------------------------------------------------
  211. C            CALCUL MATRICE TANGENTE ELASTIQUE ANISOTROPE
  212. C----------------------------------------------------------------------
  213. C
  214. C---------------------------------Passage en notation indicielle-------
  215. 121    CONTINUE
  216.           DO 42 I=1,3
  217.           D(I,I)=ENDT(I)
  218.           DO 43 J=I+1,3
  219.              K = 9 - I -J
  220. *             IF ((I.EQ.1).AND.(J.EQ.2)) K=6
  221. *             IF ((I.EQ.1).AND.(J.EQ.3)) K=5
  222. *             IF ((I.EQ.2).AND.(J.EQ.3)) K=4
  223.              D(I,J)=ENDT(K)
  224.              D(J,I)=ENDT(K)
  225. 43        CONTINUE
  226. 42     CONTINUE
  227.  
  228. C----------------------------Calcul de la matrice tangente-------------
  229.        DO 905 I=1,3
  230.           DO 904 J=1,3
  231.              DO 903 K=1,3
  232.                 DO 51 L=1,3
  233. C-------------------------
  234.          U(I,J,K,L)=LAMB*DELT(I,J)*DELT(K,L)+MUT*(DELT(I,K)*DELT(J,L)
  235.      * +DELT(I,L)*DELT(J,K))+ALP1*(DELT(I,J)*D(K,L)+D(I,J)*DELT(K,L))
  236.      *    +BETA*(D(I,K)*DELT(J,L)+D(L,J)*DELT(I,K)+DELT(I,L)
  237.      * *D(J,K)+D(I,L)*DELT(J,K))
  238. C-------------------------
  239. 51     CONTINUE
  240.   903        CONTINUE
  241.   904     CONTINUE
  242.   905  CONTINUE
  243. C-------------------------------Notation de VOIGT--------------------
  244.        IF(MMI.EQ.1)GOTO 3
  245.        DO 906 I=1,6
  246.           M=MIII(I)
  247.           N=NIII(I)
  248.           DO 134 J=1,6
  249.             K= MIII(J)
  250.             L= NIII(J)
  251.             VDEP(I,J)=U(M,N,K,L)
  252. 134       CONTINUE
  253.   906  CONTINUE
  254. C-------------------------------------------------------------------
  255.         DO 907 I=1,6
  256.         DO 142 J=1,6
  257.              VDEP1(I,J)=VDEP(I,J)
  258.              VDEP2(I,J)=VDEP(I,J)
  259. 142     CONTINUE
  260.   907   CONTINUE
  261. C********************************************************************
  262. C
  263. C********************************************************************
  264. 899   DO 321 I=1,6
  265.         IF(ABS(SIGD(I)).LE.1E-15)SIGD(I)=0.D0
  266. 321   CONTINUE
  267. C--------------------------------------------------------------------
  268.       CALL VECT(SIGD,VAL,VP,AAA(442),AAA(451),AAA(457))
  269.       DO 332 I=4,6
  270.       EPSPLS(I)=EPSPLS(I)/2
  271.   332 CONTINUE
  272.       CALL VECT(EPSPLS,VAL2,VP2,AAA(442),AAA(451),AAA(457))
  273.       DO 333 I=4,6
  274.       EPSPLS(I)=EPSPLS(I)*2
  275.   333 CONTINUE
  276.       CALL VECT(ENDD,VAL3,VP3,AAA(442),AAA(451),AAA(457))
  277. C--------------------------------------------------------------------
  278.           CALL INVERE(VDEP2,6,6,KM,DET)
  279. C------------
  280.       DO 908 I=1,6
  281.           EPSTT(I)=0.D0
  282.       DO 401 J=1,6
  283.           EPSTT(I)=EPSTT(I)+VDEP2(I,J)*(-G*ENDD(J))
  284. 401     CONTINUE
  285.   908 CONTINUE
  286.        RETURN
  287.        END
  288.  
  289.  
  290.  
  291.  
  292.  

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