Télécharger fibet2.eso

Retour à la liste

Numérotation des lignes :

  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
  102. CALL NEWBE2
  103. . ( FC___, EZERO, EPS85, ZETA_, FLT__, FTPL_,
  104. . TRAFA, EPSON, EPSO1, EPMAX, EPMIN, STRES,
  105. . STRE1, STMIN, TANGE, TANGD, FACT1, FACT2,
  106. . FACX2, FACLO, EPLA1, ET___ )
  107. C
  108. IF (EPSON .LT. EPLA1 .OR. STMIN .GT. 0.0D0) DEPS(1) = 999D9
  109. C
  110. VARF(1) = EPSON * EZERO
  111. VARF(2) = EPMAX * EZERO
  112. VARF(3) = EPMIN * EZERO
  113. VARF(4) = STMIN * FC___
  114. VARF(5) = TANGE
  115. VARF(6) = TANGD
  116. C
  117. SIGF(1) = STRES * FC___
  118. SIGF(2) = SIG0(2) + GAMMA * DEPS(2)
  119. SIGF(3) = SIG0(3) + GAMMA * DEPS(3)
  120. C
  121. RETURN
  122. END
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  

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