Télécharger la plaquette Cast3M


fr en

Télécharger barpar.eso

Retour à la liste

Numérotation des lignes : 

barpar

  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