Télécharger fibste.eso

Retour à la liste

Numérotation des lignes :

  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
  84. CALL STEEL1
  85. . ( EYOUN, FSY__, FSU__, EPSH_, EPSU_, ROFAC, BFACT,
  86. . A1FAC, A2FAC, FACLD, A6FAC, CFACT, AFACT, IFLAG,
  87. . STRAI, RFACT, EPSON, SIGMA, TANGE, STRES )
  88. C
  89. VARF(1) = DBLE(IFLAG)
  90. VARF(2) = STRAI
  91. VARF(3) = RFACT
  92. VARF(4) = TANGE
  93. C
  94. DO IFACT = 1, 2
  95. VARF(4+IFACT) = EPSON(IFACT)
  96. VARF(8+IFACT) = SIGMA(IFACT)
  97. END DO
  98. C
  99. DO IFACT = 4, 5
  100. VARF(3+IFACT) = EPSON(IFACT)
  101. VARF(7+IFACT) = SIGMA(IFACT)
  102. END DO
  103. C
  104. SIGF(1) = STRES
  105. SIGF(2) = SIG0(2) + GAMMA * DEPS(2)
  106. SIGF(3) = SIG0(3) + GAMMA * DEPS(3)
  107. C
  108. RETURN
  109. END
  110.  
  111.  
  112.  
  113.  
  114.  

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