Télécharger fibeto.eso

Retour à la liste

Numérotation des lignes :

  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
  86. CALL NEWBET
  87. . ( FC___, EZERO, FLT__, ZETA_, IFLAG, EPSON,
  88. . EPSO1, EPMAX, EPLAS, EPS85, SLOPE, TANGE,
  89. . STRES, STRE1, TRAFA )
  90. C
  91. VARF(1) = DBLE(IFLAG)
  92. VARF(2) = EPSON * EZERO
  93. VARF(3) = EPMAX * EZERO
  94. VARF(4) = EPLAS * EZERO
  95. VARF(5) = SLOPE
  96. VARF(6) = TANGE
  97. C
  98. SIGF(1) = STRES * FC___
  99. SIGF(2) = SIG0(2) + GAMMA * DEPS(2)
  100. SIGF(3) = SIG0(3) + GAMMA * DEPS(3)
  101. C
  102. RETURN
  103. END
  104.  
  105.  
  106.  
  107.  
  108.  

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