Esope source files

   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 :