Télécharger babael.eso

Retour à la liste

Numérotation des lignes :

  1. C BABAEL SOURCE CHAT 05/01/12 21:31:17 5004
  2. C BABAEL SOURCE DC 98/04/28
  3. SUBROUTINE BABAEL(XMAT,XCAR,DEPST,VAR0,SIGF,VARF,DEFP)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION XMAT(5),VAR0(3),VARF(3)
  7. C----------------------------------------------------------------------
  8. C
  9. C COMPORTEMENT DE FIBRE
  10. C MODELE BETON PLASTIQUE DU BAEL
  11. C ELEMENT DE BARRE
  12. C----------------------------------------------------------------------
  13. C Didier COMBESCURE - EMSI - Avril 1999
  14. C----------------------------------------------------------------------
  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) FC
  21. C
  22. C SIGF STRESS
  23. C DEPST,DEFP
  24. C
  25. C VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  26. C VAR0(2) VARF(2) EPSO DEFORMATION TOTALE
  27. C
  28. C----------------------------------------------------------------------
  29. C
  30. GAMMA=XMAT(1)/2/(XMAT(2)+1)
  31. EYOUN = XMAT( 1)
  32. FC = (-1.D0)*XMAT(5)
  33. XEPSLI = (2.D0)*(FC/EYOUN)
  34. C
  35. XEPSE = VAR0(1)
  36. XEPSO = VAR0(2)
  37. C
  38. XEPSOF = XEPSO + DEPST
  39. C
  40. IF (XEPSOF. GE . XEPSE) THEN
  41. C
  42. XSIGF = 0.D0
  43. XTAN = 0.D0
  44. C
  45. ELSE
  46. C
  47. XSIGF1 = EYOUN * (XEPSOF - XEPSE)
  48. C
  49. IF (XEPSOF . LT . XEPSLI) THEN
  50. XSIGF2 = FC
  51. XTAN = 0.D0
  52. ELSE
  53. XSIGF2 = FC*( ( (2.D0)*(XEPSOF/XEPSLI))
  54. & - ((XEPSOF/XEPSLI)**2) )
  55. XTAN = FC*( (2.D0/XEPSLI)
  56. & - (( 2.D0 * XEPSOF) /(XEPSLI**2) ) )
  57. ENDIF
  58. C
  59. IF (XSIGF2 . GE . XSIGF1) THEN
  60. XSIGF = XSIGF2
  61. XEPSE = XEPSOF - (XSIGF/EYOUN)
  62. ELSE
  63. XSIGF = XSIGF1
  64. XTAN = EYOUN
  65. ENDIF
  66. C
  67. ENDIF
  68. C
  69. C
  70. SIGF = XSIGF * XCAR
  71. C
  72. VARF(1) = XEPSE
  73. VARF(2) = XEPSOF
  74. VARF(3) = XTAN * XCAR
  75. C
  76. DEFP = XEPSE
  77. C
  78. RETURN
  79. END
  80.  
  81.  
  82.  
  83.  

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