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    PV        17/12/08    21:17:23     9660           
  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.        SEGMENT IECOU
  34.        INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7,
  35.      1 icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16,
  36.      2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  37.      3 icow25,icow26,icow27,icow28,icow29,icow30,ICARA,
  38.      4 icow32,icow33,NSTRS1,MFR1,icow36,icow37,icow38,
  39.      5 icow39,icow40,icow41,icow42,icow43,icow44
  40.  
  41.        INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  42.      1 icow51,icow52,icow53,icow54,icow55,icow56
  43.      2 icow57,icow58
  44.        ENDSEGMENT
  45.  
  46.       SEGMENT WRKK2
  47.          REAL*8 EPSILI(NSTRSV)
  48.       END SEGMENT
  49.  
  50.       SEGMENT WRK6
  51.         REAL*8 SIG0S(NSTRS),DEPSTS(NSTRS)
  52.       END SEGMENT
  53.  
  54.       INTEGER NVARI,KCAS,IRTD,ISTRS
  55.  
  56.       REAL*8 PREC,EPAI,FAC,AUX,AUX1,AUX2,YOUN,XNU,UN
  57.  
  58.       PARAMETER (UN=1.D0)
  59.  
  60.       KERRE=0
  61. C
  62. C----MISE EN DONNEES----------------------------------------------------
  63. C
  64.       YOUN = XMAT(1)
  65.       XNU  = XMAT(2)
  66. C
  67. C----CAS DE LA FORMULATION COQUES MINCES--------------------------------
  68. C
  69.       IF (MFR1 .EQ. 9) THEN
  70.  
  71.          EPAI=XCARB(1)
  72.          SEGINI WRK6
  73.  
  74.          DO 101 ISTRS=1,NSTRS1
  75.             SIG0S(ISTRS)=SIG0(ISTRS)
  76.             DEPSTS(ISTRS)=DEPST(ISTRS)
  77.   101    CONTINUE
  78.  
  79.          NSTRSV=4
  80.          IFOUR2=-2
  81.          JDIM =2
  82.  
  83.          DO 102 ISTRS=1,2
  84.             SIG0(ISTRS)=SIG0S(ISTRS)/EPAI
  85.   102    CONTINUE
  86.  
  87.          DEPST(3)=0.D0
  88.          DEPST(4)=DEPSTS(3)
  89.  
  90.          SIG0(3)=0.D0
  91.          SIG0(4)=SIG0S(3)/EPAI
  92. C
  93. C----CAS DE LA FORMULATION MASSIVE--------------------------------------
  94. C
  95.       ELSE IF (MFR1 .EQ. 1) THEN
  96.  
  97.          NSTRSV=NSTRS1
  98.          IFOUR2=IFOUR
  99.          JDIM = IDIM
  100.  
  101.       ELSE
  102.  
  103.          PRINT*,'MFR1=',MFR1
  104.          KERRE=57
  105.          RETURN
  106.  
  107.       END IF
  108.  
  109.       SEGINI WRKK2
  110.  
  111.       CMATE = 'ISOTROPE'
  112.       KCAS=1
  113.  
  114.       CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,KCAS,DDHOOK,IRTD)
  115.  
  116.       IF ( IRTD .EQ. 1) THEN
  117.  
  118.          CALL MATVE1 (DDHOOK,DEPST,NSTRSV,NSTRSV,DSIGT,1)
  119.  
  120.          PREC=1.D-08
  121.  
  122.          CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,2,DDHOOK,IRTD)
  123.  
  124.          CALL INVALM(DDHOOK,NSTRSV,NSTRSV,IRTD,PREC)
  125.  
  126.          IF (IRTD.EQ.0)THEN
  127.  
  128.          CALL MATVE1 (DDHOOK,SIG0,NSTRSV,NSTRSV,EPSILI,1)
  129.  
  130.          icarbi=icara
  131. C
  132. C----CALCUL DES CONTRAINTES DE CAUCHY DANS LE CAS 3D--------------------
  133. C
  134.          CALL  ELOCRAK2(WRK52,WRK53,WRK54,WRKK2,NSTRSV,NVARI,
  135.      1                  ICARbi,JDIM,IFOUR2)
  136.  
  137.          icara=icarbi
  138. C
  139. C----CALCUL DES CONTRAINTES VRAIES DANS LE CAS DES COQUES MINCES--------
  140. C
  141.          IF (MFR1 .EQ. 9) THEN
  142.  
  143.             FAC=(EPAI**3)/12.D0
  144.             AUX =FAC*YOUN/(1.D0-XNU*XNU)
  145.             AUX1=FAC*YOUN*.5D0/(1.D0+XNU)
  146.             AUX2=EPAI*YOUN*.5D0/(1.d0+XNU)/1.2d0
  147.  
  148.             DO 103 ISTRS=1,NSTRS1
  149.                SIG0(ISTRS)=SIG0S(ISTRS)
  150.                DEPST(ISTRS)=DEPSTS(ISTRS)
  151.  103        CONTINUE
  152.  
  153.             DO 104 ISTRS=1,2
  154.                SIGF(ISTRS)=SIGF(ISTRS)*EPAI
  155.  104        CONTINUE
  156.  
  157.             SIGF(3)=SIGF(4)*EPAI
  158.             SIGF(4)=SIG0(4)+AUX*(DEPST(4)+XNU*DEPST(5))
  159.             SIGF(5)=SIG0(5)+AUX*(DEPST(5)+XNU*DEPST(4))
  160.             SIGF(6)=SIG0(6)+AUX1*DEPST(6)
  161.             SIGF(7)=SIG0(7)+AUX2*DEPST(7)
  162.             SIGF(8)=SIG0(8)+AUX2*DEPST(8)
  163.             SEGSUP WRK6
  164.          END IF
  165.  
  166.          ELSE
  167.             PRINT*,'ERREUR DANS INVALM'
  168.             KERRE=56
  169.          END IF
  170.       ELSE
  171.          PRINT*,'ERREUR DANS DOHMAS'
  172.          KERRE=56
  173.       END IF
  174.       SEGSUP WRKK2
  175.       RETURN
  176.       END
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  

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