Télécharger fibael.eso

Retour à la liste

Numérotation des lignes :

  1. C FIBAEL SOURCE CHAT 12/04/06 21:15:09 7348
  2. C FIBFRA SOURCE DC 98/01/26
  3. SUBROUTINE FIBAEL(XMAT,DEPS,SIG0,VAR0,SIGF,VARF)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION XMAT(*),DEPS(3),SIG0(3),VAR0(3),SIGF(3),VARF(3)
  7. C----------------------------------------------------------------------
  8. C
  9. C COMPORTEMENT DE FIBRE
  10. C MODELE BETON PLASTIQUE DU BAEL
  11. C----------------------------------------------------------------------
  12. C Didier COMBESCURE - EMSI - Janvier 1998
  13. C----------------------------------------------------------------------
  14. C
  15. C XMAT( 1) YOUNG
  16. C XMAT( 2) NU
  17. C
  18. C XMAT( 5) FC
  19. C
  20. C DEPS(1) EPXX SIG0(1) SIGF(1) SMXX
  21. C DEPS(2) GAXY SIG0(2) SIGF(2) SMXY
  22. C DEPS(3) GAXZ SIG0(3) SIGF(3) SMXZ
  23. C
  24. C VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  25. C VAR0(2) VARF(2) EPSO DEFORMATION TOTALE
  26. C
  27. C----------------------------------------------------------------------
  28. GAMMA=XMAT(1)/2.D0/(XMAT(2)+1.D0)
  29. EYOUN = XMAT( 1)
  30. FC = (-1.D0)*XMAT(5)
  31. XEPSLI = (2.D0)*(FC/EYOUN)
  32. C
  33. XEPSE = VAR0(1)
  34. XEPSO = VAR0(2)
  35. C
  36. XEPSOF = XEPSO + DEPS(1)
  37. C
  38. IF (XEPSOF. GE . XEPSE) THEN
  39. C
  40. XSIGF = 0.D0
  41. XTAN = 0.D0
  42. C
  43. ELSE
  44. C
  45. XSIGF1 = EYOUN * (XEPSOF - XEPSE)
  46. C
  47. IF (XEPSOF . LT . XEPSLI) THEN
  48. XSIGF2 = FC
  49. XTAN = 0.D0
  50. ELSE
  51. XSIGF2 = FC*( ( (2.D0)*(XEPSOF/XEPSLI))
  52. & - ((XEPSOF/XEPSLI)**2) )
  53. XTAN = FC*( (2.D0/XEPSLI)
  54. & - (( 2.D0 * XEPSOF) /(XEPSLI**2) ) )
  55. ENDIF
  56. C
  57. IF (XSIGF2 . GE . XSIGF1) THEN
  58. XSIGF = XSIGF2
  59. XEPSE = XEPSOF - (XSIGF/EYOUN)
  60. ELSE
  61. XSIGF = XSIGF1
  62. XTAN = EYOUN
  63. ENDIF
  64. C
  65. ENDIF
  66. C
  67. C
  68. SIGF(1) = XSIGF
  69. C
  70. SIGF(2) = SIG0(2) + GAMMA * DEPS(2)
  71. SIGF(3) = SIG0(3) + GAMMA * DEPS(3)
  72. C
  73. VARF(1) = XEPSE
  74. VARF(2) = XEPSOF
  75. VARF(3) = XTAN
  76. C
  77. RETURN
  78. END
  79.  
  80.  
  81.  
  82.  
  83.  

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