Esope source files

   1 : C BARBET    SOURCE    CHAT      05/01/12    21:33:18     5004           
   2 :       SUBROUTINE BARBET(XMAT,XCAR,DEPST,SIG0,VAR0,SIGF,VARF,DEFP)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 :       DIMENSION XMAT(17),VAR0(6),VARF(6)
   6 : C----------------------------------------------------------------------
   7 : C
   8 : C     COMPORTEMENT UNIAXIAL DE BETON_LIBRE ET DE BETON_CONFINE
   9 : C
  10 : C----------------------------------------------------------------------
  11 : C      Pierre Pegon (ISPRA) Juillet 1994 +- clone de FIBETO
  12 : C----------------------------------------------------------------------
  13 : C
  14 : C      XCAR     SECT
  15 : C
  16 : C      XMAT( 1) YOUNG
  17 : C      XMAT( 2) NU
  18 : C      XMAT( 3) RHO
  19 : C      XMAT( 4) ALPHA
  20 : C      XMAT( 5) MAXIMUM COMPRESSION STRESS
  21 : C      XMAT( 6) STRAIN OF THE MAXIMUM COMPRESSION STRESS
  22 : C      XMAT( 7) MAXIMUM TENSION STRESS
  23 : C      XMAT( 8) ALPHA 1 FACTOR
  24 : C      XMAT( 9) OMEGA 1 FACTOR
  25 : C      XMAT(10) MODULUS OF THE DESCEND. PART OF THE CURVE
  26 : C      XMAT(11) STRESS VALUE DEFINING THE COMPRESSION PLATEAU
  27 : C      XMAT(12) DESCENDING TRACTION CURVE FACTOR
  28 : C      XMAT(13) F2 (FOR LOADING/UNLOADING)  
  29 : C      XMAT(14) F1 (FOR LOADING/UNLOADING)
  30 : C      XMAT(15) RESIDUAL STRESS (TENSION)
  31 : C      XMAT(16) F1' (FOR LOADING/UNLOADING) 
  32 : C      XMAT(17) F2' (FOR LOADING/UNLOADING) 
  33 : C
  34 : C      SIG0 AND SIGF STRESS
  35 : C      DEPST,DEFP
  36 : C
  37 : C                      VAR0(1) VARF(1) FLAG
  38 : C                      VAR0(2) VARF(2) EPSO
  39 : C                      VAR0(3) VARF(3) MAX. STRAIN EVER REACHED
  40 : C                      VAR0(4) VARF(4) PLASTIC STRAIN
  41 : C                      VAR0(5) VARF(5) MOD. OF THE UNLOAD. CURVE
  42 : C                      VAR0(6) VARF(6) TANGENT MODULUS
  43 : C
  44 : C----------------------------------------------------------------------
  45 : C
  46 :        DEPS = DEPST
  47 :        FC___ = XMAT( 5)
  48 :        EZERO = XMAT( 6)
  49 :        FLT__ = XMAT( 7)
  50 :        ALFA1 = XMAT( 8)
  51 :        OMEG1 = XMAT( 9)
  52 :        ZETA_ = XMAT(10)
  53 :        STR85 = XMAT(11)
  54 :        TRAFA = XMAT(12)
  55 :        FTPL_ = XMAT(15)
  56 : C       FACT1 = 20.0
  57 : C       FACT2 = 20.0
  58 :        FACLO = XMAT(13)
  59 :        FACX2 = XMAT(14)
  60 :        FACT1 = XMAT(16)
  61 :        FACT2 = XMAT(17)
  62 : C
  63 :        IF (TRAFA .LT. 1.0) THEN
  64 :           WRITE (*,*) ' ERROR - TRAFA < 1.0 '
  65 :           STOP
  66 :        ENDIF
  67 : C
  68 :        ALFOM = ALFA1 * OMEG1
  69 : C
  70 :        IF (ZETA_ .LE. 0.) THEN
  71 : C
  72 : C ------- COMPUTING PROPERTIES OF CONFINED CONCRETE -------
  73 : C
  74 :           IF (ALFOM .LT. 1.0D-1) THEN
  75 :              BETAC = 1.000D0 + 2.50D0 * ALFOM
  76 :           ELSE
  77 :              BETAC = 1.125D0 + 1.25D0 * ALFOM
  78 :           ENDIF
  79 : C
  80 :           FC___ = BETAC * FC___
  81 :           EZERO = BETAC * BETAC * EZERO
  82 :           ZETA_ = (BETAC-0.85D0) / (BETAC*(0.1*ALFOM+0.0035-EZERO))
  83 :        ENDIF
  84 : C
  85 :        EPS85 = (1.0-STR85/FC___ ) / ZETA_ + EZERO
  86 :        EPS85 = -1.0 * EPS85 / EZERO
  87 : C
  88 :        ET___ = FLT__ / (2.0 * FC___ / EZERO)
  89 : C
  90 :        EPMXM = 0.0
  91 : C
  92 :        EPSON = (VAR0(1) + DEPS) / EZERO
  93 :        EPSO1 = VAR0(1) / EZERO
  94 :        EPMAX = VAR0(2) / EZERO
  95 :        EPMIN = VAR0(3) / EZERO
  96 :        STMIN = VAR0(4) / FC___
  97 :        TANGE = VAR0(5)
  98 :        TANGD = VAR0(6)
  99 :        STRE1 = SIG0 / FC___
 100 : C
 100 :        CALL NEWBE2
 101 :      .       ( FC___, EZERO, EPS85, ZETA_, FLT__, FTPL_,
 102 :      .         TRAFA, EPSON, EPSO1, EPMAX, EPMIN, STRES,
 103 :      .         STRE1, STMIN, TANGE, TANGD, FACT1, FACT2,
 104 :      .         FACX2, FACLO, EPLA1, ET___ )
 105 : 
 106 : C
 107 :        IF (EPSON .LT. EPLA1 .OR. STMIN .GT. 0.0) DEPS = 999E9
 108 : C
 109 :        VARF(1) = EPSON * EZERO
 110 :        VARF(2) = EPMAX * EZERO
 111 :        VARF(3) = EPMIN * EZERO
 112 :        VARF(4) = STMIN * FC___
 113 :        VARF(5) = TANGE
 114 :        VARF(6) = TANGD
 115 : C
 116 :        SIGF    = STRES * (FC___ * XCAR)
 117 :        DEFP    = 0.
 118 : C
 119 :        RETURN
 120 :        END
 121 :  
 122 :