Télécharger la plaquette Cast3M


fr en

Télécharger fibael.eso

Retour à la liste

Numérotation des lignes : 

fibael

  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