Télécharger la plaquette Cast3M


fr en

Télécharger elocrak1.eso

Retour à la liste

Numérotation des lignes : 

elocrak1

  1. C ELOCRAK1  SOURCE    OF166741  25/11/04    21:15:53     12349          
  2.       SUBROUTINE ELOCRAK1(WRK52,WRK53,WRK54,NVARI,Iecou)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires :   Subroutine permettant de mettre en oeuvre le
  6. C                          modele ELOCRACK (CEA) pour representer
  7. C                  le comportement cyclique du béton sous
  8. C                  sollicitations sismiques
  9. C
  10. C Traits       : -
  11. C                -
  12. C                -
  13. C                -
  14. C                -
  15. C
  16. C Auteur       : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
  17. C====&===1=========2=========3=========4=========5=========6=========7==
  18. C
  19. C----DECLARATION GENERALES----------------------------------------------
  20. C
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8(A-H,O-Z)
  23. C
  24. C----LES INCLUDES-------------------------------------------------------
  25. C
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29. -INC DECHE
  30. C
  31. C----DECLARATION PARAMETRES---------------------------------------------
  32. C
  33. -INC TECOU
  34.  
  35.       SEGMENT WRKK2
  36.          REAL*8 EPSILI(NSTRSV)
  37.       END SEGMENT
  38.  
  39.       SEGMENT WRK6
  40.         REAL*8 SIG0S(NSTRS),DEPSTS(NSTRS)
  41.       END SEGMENT
  42.  
  43.       INTEGER NVARI,KCAS,IRTD,ISTRS
  44.  
  45.       REAL*8 PREC,EPAI,FAC,AUX,AUX1,AUX2,YOUN,XNU,UN
  46.  
  47.       PARAMETER (UN=1.D0)
  48.  
  49.       KERRE=0
  50. C
  51. C----MISE EN DONNEES----------------------------------------------------
  52. C
  53.       YOUN = XMAT(1)
  54.       XNU  = XMAT(2)
  55.  
  56.       NSTRS1 = iecou.NSTRSS
  57. C
  58. C----CAS DE LA FORMULATION COQUES MINCES--------------------------------
  59. C
  60.       IF (MFR1 .EQ. 9) THEN
  61.  
  62.          EPAI=XCARB(1)
  63.          SEGINI WRK6
  64.  
  65.          DO 101 ISTRS=1,NSTRS1
  66.             SIG0S(ISTRS)=SIG0(ISTRS)
  67.             DEPSTS(ISTRS)=DEPST(ISTRS)
  68.   101    CONTINUE
  69.  
  70.          NSTRSV=4
  71.          IFOUR2=-2
  72.          JDIM =2
  73.  
  74.          SIG0(1)=SIG0S(1)/EPAI
  75.          SIG0(2)=SIG0S(2)/EPAI
  76.          SIG0(3)=0.D0
  77.          SIG0(4)=SIG0S(3)/EPAI
  78.  
  79.          DEPST(3)=0.D0
  80.          DEPST(4)=DEPSTS(3)
  81.  
  82. C----CAS DE LA FORMULATION MASSIVE--------------------------------------
  83. C
  84.       ELSE IF (MFR1 .EQ. 1) THEN
  85.  
  86.          NSTRSV=NSTRS1
  87.          IFOUR2=IFOUR
  88.          JDIM = IDIM
  89.  
  90.       ELSE
  91.  
  92.          PRINT*,'MFR1=',MFR1
  93.          KERRE=57
  94.          RETURN
  95.  
  96.       END IF
  97.  
  98.       SEGINI WRKK2
  99.  
  100.       CMATE = 'ISOTROPE'
  101.       KCAS=1
  102.  
  103.       CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,KCAS,DDHOOK,IRTD)
  104.  
  105.       IF ( IRTD .EQ. 1) THEN
  106.  
  107.          CALL MATVE1 (DDHOOK,DEPST,NSTRSV,NSTRSV,DSIGT,1)
  108.  
  109.          PREC=1.D-08
  110.  
  111.          CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,2,DDHOOK,IRTD)
  112.  
  113.          CALL INVALM(DDHOOK,NSTRSV,NSTRSV,IRTD,PREC)
  114.  
  115.          IF (IRTD.EQ.0)THEN
  116.  
  117.          CALL MATVE1 (DDHOOK,SIG0,NSTRSV,NSTRSV,EPSILI,1)
  118.  
  119.  
  120. C----CALCUL DES CONTRAINTES DE CAUCHY DANS LE CAS 3D--------------------
  121. C
  122.          icarbi=iecou.icara
  123.  
  124.          CALL  ELOCRAK2(WRK52,WRK53,WRK54,WRKK2,NSTRSV,NVARI,
  125.      1                  ICARbi,JDIM,IFOUR2)
  126.  
  127. C----CALCUL DES CONTRAINTES VRAIES DANS LE CAS DES COQUES MINCES--------
  128. C
  129.          IF (MFR1 .EQ. 9) THEN
  130.  
  131.             FAC=(EPAI**3)/12.D0
  132.             AUX =FAC*YOUN/(1.D0-XNU*XNU)
  133.             AUX1=FAC*YOUN*.5D0/(1.D0+XNU)
  134.             AUX2=EPAI*YOUN*.5D0/(1.d0+XNU)/1.2d0
  135.  
  136.             DO 103 ISTRS=1,NSTRS1
  137.                SIG0(ISTRS)=SIG0S(ISTRS)
  138.                DEPST(ISTRS)=DEPSTS(ISTRS)
  139.  103        CONTINUE
  140.  
  141.             SIGF(1)=SIGF(1)*EPAI
  142.             SIGF(2)=SIGF(2)*EPAI
  143.             SIGF(3)=SIGF(4)*EPAI
  144.             SIGF(4)=SIG0(4)+AUX*(DEPST(4)+XNU*DEPST(5))
  145.             SIGF(5)=SIG0(5)+AUX*(DEPST(5)+XNU*DEPST(4))
  146.             SIGF(6)=SIG0(6)+AUX1*DEPST(6)
  147.             SIGF(7)=SIG0(7)+AUX2*DEPST(7)
  148.             SIGF(8)=SIG0(8)+AUX2*DEPST(8)
  149.             SEGSUP WRK6
  150.          END IF
  151.  
  152.          ELSE
  153.             PRINT*,'ERREUR DANS INVALM'
  154.             KERRE=56
  155.          END IF
  156.       ELSE
  157.          PRINT*,'ERREUR DANS DOHMAS'
  158.          KERRE=56
  159.       END IF
  160.       SEGSUP WRKK2
  161.  
  162.       RETURN
  163.       END
  164.  
  165.  
  166.  

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