Télécharger la plaquette Cast3M


fr en

Télécharger babael.eso

Retour à la liste

Numérotation des lignes : 

babael

  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