Esope source files

   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 :