Télécharger barbet.eso

Retour à la liste

Numérotation des lignes :

  1. C BARBET SOURCE CHAT 05/01/12 21:33:18 5004
  2. SUBROUTINE BARBET(XMAT,XCAR,DEPST,SIG0,VAR0,SIGF,VARF,DEFP)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. DIMENSION XMAT(17),VAR0(6),VARF(6)
  6. C----------------------------------------------------------------------
  7. C
  8. C COMPORTEMENT UNIAXIAL DE BETON_LIBRE ET DE BETON_CONFINE
  9. C
  10. C----------------------------------------------------------------------
  11. C Pierre Pegon (ISPRA) Juillet 1994 +- clone de FIBETO
  12. C----------------------------------------------------------------------
  13. C
  14. C XCAR SECT
  15. C
  16. C XMAT( 1) YOUNG
  17. C XMAT( 2) NU
  18. C XMAT( 3) RHO
  19. C XMAT( 4) ALPHA
  20. C XMAT( 5) MAXIMUM COMPRESSION STRESS
  21. C XMAT( 6) STRAIN OF THE MAXIMUM COMPRESSION STRESS
  22. C XMAT( 7) MAXIMUM TENSION STRESS
  23. C XMAT( 8) ALPHA 1 FACTOR
  24. C XMAT( 9) OMEGA 1 FACTOR
  25. C XMAT(10) MODULUS OF THE DESCEND. PART OF THE CURVE
  26. C XMAT(11) STRESS VALUE DEFINING THE COMPRESSION PLATEAU
  27. C XMAT(12) DESCENDING TRACTION CURVE FACTOR
  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 SIG0 AND SIGF STRESS
  35. C DEPST,DEFP
  36. C
  37. C VAR0(1) VARF(1) FLAG
  38. C VAR0(2) VARF(2) EPSO
  39. C VAR0(3) VARF(3) MAX. STRAIN EVER REACHED
  40. C VAR0(4) VARF(4) PLASTIC STRAIN
  41. C VAR0(5) VARF(5) MOD. OF THE UNLOAD. CURVE
  42. C VAR0(6) VARF(6) TANGENT MODULUS
  43. C
  44. C----------------------------------------------------------------------
  45. C
  46. DEPS = DEPST
  47. FC___ = XMAT( 5)
  48. EZERO = XMAT( 6)
  49. FLT__ = XMAT( 7)
  50. ALFA1 = XMAT( 8)
  51. OMEG1 = XMAT( 9)
  52. ZETA_ = XMAT(10)
  53. STR85 = XMAT(11)
  54. TRAFA = XMAT(12)
  55. FTPL_ = XMAT(15)
  56. C FACT1 = 20.0
  57. C FACT2 = 20.0
  58. FACLO = XMAT(13)
  59. FACX2 = XMAT(14)
  60. FACT1 = XMAT(16)
  61. FACT2 = XMAT(17)
  62. C
  63. IF (TRAFA .LT. 1.0) THEN
  64. WRITE (*,*) ' ERROR - TRAFA < 1.0 '
  65. STOP
  66. ENDIF
  67. C
  68. ALFOM = ALFA1 * OMEG1
  69. C
  70. IF (ZETA_ .LE. 0.) THEN
  71. C
  72. C ------- COMPUTING PROPERTIES OF CONFINED CONCRETE -------
  73. C
  74. IF (ALFOM .LT. 1.0D-1) THEN
  75. BETAC = 1.000D0 + 2.50D0 * ALFOM
  76. ELSE
  77. BETAC = 1.125D0 + 1.25D0 * ALFOM
  78. ENDIF
  79. C
  80. FC___ = BETAC * FC___
  81. EZERO = BETAC * BETAC * EZERO
  82. ZETA_ = (BETAC-0.85D0) / (BETAC*(0.1*ALFOM+0.0035-EZERO))
  83. ENDIF
  84. C
  85. EPS85 = (1.0-STR85/FC___ ) / ZETA_ + EZERO
  86. EPS85 = -1.0 * EPS85 / EZERO
  87. C
  88. ET___ = FLT__ / (2.0 * FC___ / EZERO)
  89. C
  90. EPMXM = 0.0
  91. C
  92. EPSON = (VAR0(1) + DEPS) / EZERO
  93. EPSO1 = VAR0(1) / EZERO
  94. EPMAX = VAR0(2) / EZERO
  95. EPMIN = VAR0(3) / EZERO
  96. STMIN = VAR0(4) / FC___
  97. TANGE = VAR0(5)
  98. TANGD = VAR0(6)
  99. STRE1 = SIG0 / FC___
  100. 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.  
  107. C
  108. IF (EPSON .LT. EPLA1 .OR. STMIN .GT. 0.0) DEPS = 999E9
  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 = STRES * (FC___ * XCAR)
  118. DEFP = 0.
  119. C
  120. RETURN
  121. END
  122.  
  123.  
  124.  

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