1 : C FIBETO SOURCE CHAT 12/04/06 21:15:10 7348 2 : SUBROUTINE FIBETO(XMAT,DEPS,SIG0,VAR0,SIGF,VARF) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : DIMENSION XMAT(18),DEPS(3),SIG0(3),VAR0(6),SIGF(3),VARF(6) 6 : C---------------------------------------------------------------------- 7 : C 8 : C COMPORTEMENT DE FIBRE DE BETON_LIBRE ET DE BETON_CONFINE 9 : C 10 : C---------------------------------------------------------------------- 11 : C Joao Guedes et Pierre Pegon (ISPRA) Juillet/Decembre 1993 12 : C---------------------------------------------------------------------- 13 : C 14 : C XMAT( 1) YOUNG 15 : C XMAT( 2) NU 16 : C XMAT( 3) RHO 17 : C XMAT( 4) ALPHA 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 27 : C Paramètres "inactifs" 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 DEPS(1) EPXX SIG0(1) SIGF(1) SMXX 35 : C DEPS(2) GAXY SIG0(2) SIGF(2) SMXY 36 : C DEPS(3) GAXZ SIG0(3) SIGF(3) SMXZ 37 : C 38 : C VAR0(1) VARF(1) FLAG 39 : C VAR0(2) VARF(2) EPSO 40 : C VAR0(3) VARF(3) MAX. STRAIN EVER REACHED 41 : C VAR0(4) VARF(4) PLASTIC STRAIN 42 : C VAR0(5) VARF(5) MOD. OF THE UNLOAD. CURVE 43 : C VAR0(6) VARF(6) TANGENT MODULUS 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 : C 57 : ALFOM = ALFA1 * OMEG1 58 : C 59 : IF (ZETA_ .LE. 0.D0) THEN 60 : C 61 : C ------- COMPUTING PROPERTIES OF CONFINED CONCRETE ------- 62 : C 63 : IF (ALFOM .LT. 1.0D-1) THEN 64 : BETAC = 1.000D0 + 2.50D0 * ALFOM 65 : ELSE 66 : BETAC = 1.125D0 + 1.25D0 * ALFOM 67 : ENDIF 68 : C 69 : FC___ = BETAC * FC___ 70 : EZERO = BETAC * BETAC * EZERO 71 : ZETA_ = (BETAC-0.85D0) / (BETAC*(0.1D0*ALFOM+0.0035D0-EZERO)) 72 : ENDIF 73 : C 74 : EPS85 = (1.0D0-STR85/FC___ ) / ZETA_ + EZERO 75 : EPS85 = -1.0D0 * EPS85 / EZERO 76 : C 77 : IFLAG = NINT(VAR0(1)) 78 : EPSON = (VAR0(2) + DEPS(1)) / EZERO 79 : EPSO1 = VAR0(2) / EZERO 80 : EPMAX = VAR0(3) / EZERO 81 : EPLAS = VAR0(4) / EZERO 82 : SLOPE = VAR0(5) 83 : TANGE = VAR0(6) 84 : STRE1 = SIG0(1) / FC___ 85 : C 85 : CALL NEWBET 86 : . ( FC___, EZERO, FLT__, ZETA_, IFLAG, EPSON, 87 : . EPSO1, EPMAX, EPLAS, EPS85, SLOPE, TANGE, 88 : . STRES, STRE1, TRAFA ) 89 : C 90 : VARF(1) = DBLE(IFLAG) 91 : VARF(2) = EPSON * EZERO 92 : VARF(3) = EPMAX * EZERO 93 : VARF(4) = EPLAS * EZERO 94 : VARF(5) = SLOPE 95 : VARF(6) = TANGE 96 : C 97 : SIGF(1) = STRES * FC___ 98 : SIGF(2) = SIG0(2) + GAMMA * DEPS(2) 99 : SIGF(3) = SIG0(3) + GAMMA * DEPS(3) 100 : C 101 : RETURN 102 : END 103 : 104 : 105 : 106 :
© Cast3M 2003 - All rights reserved.
Disclaimer