1 : C FIBSTE SOURCE CHAT 12/04/06 21:15:13 7348 2 : SUBROUTINE FIBSTE(XMAT,DEPS,SIG0,VAR0,SIGF,VARF) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : DIMENSION XMAT(*),DEPS(3),SIG0(3),VAR0(14),SIGF(3),VARF(14) 6 : DIMENSION EPSON(5),SIGMA(5) 7 : C---------------------------------------------------------------------- 8 : C 9 : C COMPORTEMENT DE FIBRE DE FERRAILLE 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 XMAT( 3) RHO 18 : C XMAT( 4) ALPHA 19 : C XMAT( 5) MAX. ELAST. STRESS 20 : C XMAT( 6) ULTIMATE STRESS 21 : C XMAT( 7) HARDENING STRAIN 22 : C XMAT( 8) ULTIMATE STRAIN 23 : C XMAT( 9) Ro FACTOR 24 : C XMAT(10) b FACTOR 25 : C XMAT(11) a1 FACTOR 26 : C XMAT(12) a2 FACTOR 27 : C XMAT(13) FACTOR L/D 28 : C XMAT(14) a6 FACTOR 29 : C XMAT(15) c FACTOR 30 : C XMAT(16) a FACTOR 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) FLAG 37 : C VAR0( 2) VARF( 2) EPSO 38 : C VAR0( 3) VARF( 3) R FACTOR 39 : C VAR0( 4) VARF( 4) TANGENT MODULUS 40 : C VAR0( 5) VARF( 5) LAST REV. STRAIN 41 : C VAR0( 6) VARF( 6) INTERSEPTION POINT 42 : C VAR0( 8) VARF( 7) MAX. REV. POST. STRAIN 43 : C VAR0( 9) VARF( 8) MIN. REV. POST. STRAIN 44 : C VAR0(10) VARF( 9) LAST REV. STRESS 45 : C VAR0(11) VARF(10) INTERSEPTION POINT 46 : C VAR0(13) VARF(11) MAX. REV. POST. STRESS 47 : C VAR0(14) VARF(12) MIN. REV. POST. STRESS 48 : C 49 : C---------------------------------------------------------------------- 50 : C 51 : GAMMA = XMAT( 1) / 2.0D0 / (XMAT(2)+1D0) 52 : EYOUN = XMAT( 1) 53 : FSY__ = XMAT( 5) 54 : FSU__ = XMAT( 6) 55 : EPSH_ = XMAT( 7) 56 : EPSU_ = XMAT( 8) 57 : ROFAC = XMAT( 9) 58 : BFACT = XMAT(10) 59 : A1FAC = XMAT(11) 60 : A2FAC = XMAT(12) 61 : FACLD = XMAT(13) 62 : A6FAC = XMAT(14) 63 : CFACT = XMAT(15) 64 : AFACT = XMAT(16) 65 : C 66 : IFLAG = NINT(VAR0(1)) 67 : STRAI = VAR0(2) + DEPS(1) 68 : RFACT = VAR0(3) 69 : TANGE = VAR0(4) 70 : C 71 : DO IFACT = 1, 2 72 : EPSON(IFACT) = VAR0(4+IFACT) 73 : SIGMA(IFACT) = VAR0(8+IFACT) 74 : END DO 75 : C 76 : EPSON(3) = VAR0(2) 77 : SIGMA(3) = SIG0(1) 78 : C 79 : DO IFACT = 4, 5 80 : EPSON(IFACT) = VAR0(3+IFACT) 81 : SIGMA(IFACT) = VAR0(7+IFACT) 82 : END DO 83 : C 83 : CALL STEEL1 84 : . ( EYOUN, FSY__, FSU__, EPSH_, EPSU_, ROFAC, BFACT, 85 : . A1FAC, A2FAC, FACLD, A6FAC, CFACT, AFACT, IFLAG, 86 : . STRAI, RFACT, EPSON, SIGMA, TANGE, STRES ) 87 : C 88 : VARF(1) = DBLE(IFLAG) 89 : VARF(2) = STRAI 90 : VARF(3) = RFACT 91 : VARF(4) = TANGE 92 : C 93 : DO IFACT = 1, 2 94 : VARF(4+IFACT) = EPSON(IFACT) 95 : VARF(8+IFACT) = SIGMA(IFACT) 96 : END DO 97 : C 98 : DO IFACT = 4, 5 99 : VARF(3+IFACT) = EPSON(IFACT) 100 : VARF(7+IFACT) = SIGMA(IFACT) 101 : END DO 102 : C 103 : SIGF(1) = STRES 104 : SIGF(2) = SIG0(2) + GAMMA * DEPS(2) 105 : SIGF(3) = SIG0(3) + GAMMA * DEPS(3) 106 : C 107 : RETURN 108 : END 109 : 110 : 111 : 112 :
© Cast3M 2003 - All rights reserved.
Disclaimer