Télécharger la plaquette Cast3M


fr en

Télécharger fibpar.eso

Retour à la liste

Numérotation des lignes : 

fibpar

  1. C FIBPAR    SOURCE    CHAT      12/04/06    21:15:12     7348
  2. C FIBFRA    SOURCE    DC        98/01/26
  3.       SUBROUTINE FIBPAR(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(2),SIGF(3),VARF(2)
  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) FY
  19. C      XMAT( 6) HH
  20. C
  21. C      DEPS(1) EPXX    SIG0(1) SIGF(1) SMXX
  22. C      DEPS(2) GAXY    SIG0(2) SIGF(2) SMXY
  23. C      DEPS(3) GAXZ    SIG0(3) SIGF(3) SMXZ
  24. C
  25. C                      VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  26. C
  27. C----------------------------------------------------------------------
  28.        GAMMA=XMAT(1)/2.D0/(XMAT(2)+1.D0)
  29.        EYOUN = XMAT( 1)
  30.        FYP = XMAT(5)
  31.        FYN = (-1.D0)*XMAT(5)
  32.        HHH = XMAT(6)
  33.        FYPE = FYP + ( HHH * VAR0(1))
  34.        FYNE = FYN + ( HHH * VAR0(1))
  35. C
  36.        XEPSE =  VAR0(1)
  37. C
  38.        XSIGT = SIG0(1) + EYOUN * DEPS(1)
  39. C
  40.        IF (XSIGT . GT . FYPE) THEN
  41. C
  42.          XDEP = ((EYOUN/(EYOUN+HHH))*DEPS(1))
  43.      &         + ((SIG0(1) - FYPE)/(EYOUN+HHH))
  44.          XSIGF =  SIG0(1) + (EYOUN*(DEPS(1) - XDEP))
  45.          XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  46. C
  47.        ELSEIF (XSIGT . LT . FYNE) THEN
  48. C
  49.          XDEP = ((EYOUN/(EYOUN+HHH))*DEPS(1))
  50.      &         + ((SIG0(1) - FYNE)/(EYOUN+HHH))
  51.          XSIGF =  SIG0(1) + (EYOUN*(DEPS(1) - XDEP))
  52.          XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  53. C
  54.        ELSE
  55. C
  56.          XSIGF = XSIGT
  57.          XDEP = 0.D0
  58.          XTAN = EYOUN
  59. C
  60.        ENDIF
  61. C
  62. C
  63.        SIGF(1) = XSIGF
  64. C
  65.        SIGF(2) = SIG0(2) + GAMMA   * DEPS(2)
  66.        SIGF(3) = SIG0(3) + GAMMA   * DEPS(3)
  67. C
  68.        VARF(1) = XEPSE + XDEP
  69.        VARF(2) = XTAN
  70. C
  71.        RETURN
  72.        END
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  

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