Télécharger barpar.eso

Retour à la liste

Numérotation des lignes :

  1. C BARPAR SOURCE CHAT 05/01/12 21:34:56 5004
  2. C BARPAR SOURCE DC 98/01/26
  3. SUBROUTINE BARPAR(XMAT,XCAR,DEPST,SIG0,VAR0,SIGF,VARF,DEFP)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION XMAT(6),VAR0(2),VARF(2)
  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 - Janvier 1998
  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) FY
  21. C XMAT( 6) HH
  22. C
  23. C SIGF STRESS
  24. C DEPST,DEFP
  25. C
  26. C VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  27. C
  28. C----------------------------------------------------------------------
  29. DEPS = DEPST
  30. C
  31. XSECT = XCAR
  32. C
  33. XSIG0 = SIG0/XSECT
  34. C
  35. GAMMA=XMAT(1)/2/(XMAT(2)+1)
  36. EYOUN = XMAT( 1)
  37. FYP = XMAT(5)
  38. FYN = (-1.)*XMAT(5)
  39. HHH = XMAT(6)
  40. FYPE = FYP + ( HHH * VAR0(1))
  41. FYNE = FYN + ( HHH * VAR0(1))
  42. C
  43. XEPSE = VAR0(1)
  44. C
  45. XSIGT = XSIG0 + EYOUN * DEPST
  46. C
  47. IF (XSIGT . GT . FYPE) THEN
  48. C
  49. XDEP = ((EYOUN/(EYOUN+HHH))*DEPST)
  50. & + ((XSIG0 - FYPE)/(EYOUN+HHH))
  51. XSIGF = XSIG0 + (EYOUN*(DEPST - XDEP))
  52. XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  53. C
  54. ELSEIF (XSIGT . LT . FYNE) THEN
  55. C
  56. XDEP = ((EYOUN/(EYOUN+HHH))*DEPST)
  57. & + ((XSIG0 - FYNE)/(EYOUN+HHH))
  58. XSIGF = XSIG0 + (EYOUN*(DEPST - XDEP))
  59. XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  60. C
  61. ELSE
  62. C
  63. XSIGF = XSIGT
  64. XDEP = 0.D0
  65. XTAN = EYOUN
  66. C
  67. ENDIF
  68. C
  69. C
  70. C
  71. SIGF = XSIGF * XSECT
  72. C
  73. VARF(1) = XEPSE + XDEP
  74. VARF(2) = XTAN
  75. C
  76. DEFP = XEPSE
  77. C
  78. RETURN
  79. END
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  

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