1 : C FIBET2 SOURCE CHAT 12/04/06 21:15:10 7348 2 : C FIBET2 SOURCE AM 94/02/28 21:24:59 1061 3 : SUBROUTINE FIBET2(XMAT,DEPS,SIG0,VAR0,SIGF,VARF) 4 : IMPLICIT INTEGER(I-N) 5 : IMPLICIT REAL*8(A-H,O-Z) 6 : DIMENSION XMAT(18),DEPS(3),SIG0(3),VAR0(6),SIGF(3),VARF(6) 7 : C---------------------------------------------------------------------- 8 : C 9 : C COMPORTEMENT DE FIBRE DE BETON_LIBRE ET DE BETON_CONFINE 10 : C 11 : C---------------------------------------------------------------------- 12 : C Joao Guedes et Pierre Pegon (ISPRA) Juillet/Decembre 1993 13 : C---------------------------------------------------------------------- 14 : C 15 : C XMAT( 1) YOUNG 16 : C XMAT( 2) NU 17 : C 18 : C XMAT( 5) MAXIMUM COMPRESSION STRESS 19 : C XMAT( 6) STRAIN OF THE MAXIMUM COMPRESSION STRESS 20 : C XMAT( 7) MAXIMUM TENSION STRESS 21 : C XMAT( 8) ALPHA 1 FACTOR 22 : C XMAT( 9) OMEGA 1 FACTOR 23 : C XMAT(10) MODULUS OF THE DESCEND. PART OF THE CURVE 24 : C XMAT(11) STRESS VALUE DEFINING THE COMPRESSION PLATEAU 25 : C XMAT(12) DESCENDING TRACTION CURVE FACTOR 26 : C XMAT(13) F2 (FOR LOADING/UNLOADING) 27 : C XMAT(14) F1 (FOR LOADING/UNLOADING) 28 : C XMAT(15) RESIDUAL STRESS (TENSION) 29 : C XMAT(16) F1' (FOR LOADING/UNLOADING) 30 : C XMAT(17) F2' (FOR LOADING/UNLOADING) 31 : C 32 : C DEPS(1) EPXX SIG0(1) SIGF(1) SMXX 33 : C DEPS(2) GAXY SIG0(2) SIGF(2) SMXY 34 : C DEPS(3) GAXZ SIG0(3) SIGF(3) SMXZ 35 : C 36 : C VAR0(1) VARF(1) EPSO 37 : C VAR0(2) VARF(2) MAX. STRAIN EVER REACHED 38 : C VAR0(3) VARF(3) MIN. STRAIN EVER REACHED 39 : C VAR0(4) VARF(4) MIN. STRESS EVER REACHED 40 : C VAR0(5) VARF(5) TANGENT MODULUS 41 : C VAR0(6) VARF(6) DISCHARGE TANGENT MODULUS 42 : C 43 : C---------------------------------------------------------------------- 44 : C 45 : C 46 : C 47 : GAMMA = XMAT( 1) / 2.0D0 / (XMAT(2)+1.D0) 48 : FC___ = XMAT( 5) 49 : EZERO = XMAT( 6) 50 : FLT__ = XMAT( 7) 51 : ALFA1 = XMAT( 8) 52 : OMEG1 = XMAT( 9) 53 : ZETA_ = XMAT(10) 54 : STR85 = XMAT(11) 55 : TRAFA = XMAT(12) 56 : FTPL_ = XMAT(15) 57 : C FACT1 = 20.0 58 : C FACT2 = 20.0 59 : FACLO = XMAT(13) 60 : FACX2 = XMAT(14) 61 : FACT1 = XMAT(16) 62 : FACT2 = XMAT(17) 63 : C 64 : IF (TRAFA .LT. 1.0d0) THEN 65 : WRITE (*,*) ' ERROR - TRAFA < 1.0 ' 66 : STOP 67 : ENDIF 68 : C 69 : ALFOM = ALFA1 * OMEG1 70 : C 71 : IF (ZETA_ .LE. 0.D0) THEN 72 : C 73 : C ------- COMPUTING PROPERTIES OF CONFINED CONCRETE ------- 74 : C 75 : IF (ALFOM .LT. 1.0D-1) THEN 76 : BETAC = 1.000D0 + 2.50D0 * ALFOM 77 : ELSE 78 : BETAC = 1.125D0 + 1.25D0 * ALFOM 79 : ENDIF 80 : C 81 : FC___ = BETAC * FC___ 82 : EZERO = BETAC * BETAC * EZERO 83 : ZETA_ = (BETAC-0.85D0) / (BETAC*(0.1D0*ALFOM+0.0035D0-EZERO)) 84 : ENDIF 85 : C 86 : EPS85 = (1.0D0-STR85/FC___ ) / ZETA_ + EZERO 87 : EPS85 = -1.0D0 * EPS85 / EZERO 88 : C 89 : ET___ = FLT__ / (2.0D0 * FC___ / EZERO) 90 : C 91 : EPMXM = 0.0D0 92 : C 93 : EPSON = (VAR0(1) + DEPS(1)) / EZERO 94 : EPSO1 = VAR0(1) / EZERO 95 : EPMAX = VAR0(2) / EZERO 96 : EPMIN = VAR0(3) / EZERO 97 : STMIN = VAR0(4) / FC___ 98 : TANGE = VAR0(5) 99 : TANGD = VAR0(6) 100 : STRE1 = SIG0(1) / FC___ 101 : C 101 : CALL NEWBE2 102 : . ( FC___, EZERO, EPS85, ZETA_, FLT__, FTPL_, 103 : . TRAFA, EPSON, EPSO1, EPMAX, EPMIN, STRES, 104 : . STRE1, STMIN, TANGE, TANGD, FACT1, FACT2, 105 : . FACX2, FACLO, EPLA1, ET___ ) 106 : C 107 : IF (EPSON .LT. EPLA1 .OR. STMIN .GT. 0.0D0) DEPS(1) = 999D9 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(1) = STRES * FC___ 117 : SIGF(2) = SIG0(2) + GAMMA * DEPS(2) 118 : SIGF(3) = SIG0(3) + GAMMA * DEPS(3) 119 : C 120 : RETURN 121 : END 122 : 123 : 124 : 125 : 126 :
© Cast3M 2003 - All rights reserved.
Disclaimer