Esope source files

   1 : C BARPAR    SOURCE    CHAT      05/01/12    21:34:56     5004           
   2 : C BARPAR    SOURCE    DC        98/01/26            
   3 :       SUBROUTINE BARPAR(XMAT,XCAR,DEPST,SIG0,VAR0,SIGF,VARF,DEFP)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 :       DIMENSION XMAT(6),VAR0(2),VARF(2)
   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 - Janvier 1998 
  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) FY
  21 : C      XMAT( 6) HH
  22 : C
  23 : C      SIGF STRESS
  24 : C      DEPST,DEFP
  25 : C                    
  26 : C                      VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  27 : C
  28 : C----------------------------------------------------------------------
  29 :        DEPS = DEPST
  30 : C
  31 :        XSECT = XCAR
  32 : C
  33 :        XSIG0 = SIG0/XSECT
  34 : C
  35 :        GAMMA=XMAT(1)/2/(XMAT(2)+1)
  36 :        EYOUN = XMAT( 1)
  37 :        FYP = XMAT(5)
  38 :        FYN = (-1.)*XMAT(5)
  39 :        HHH = XMAT(6)
  40 :        FYPE = FYP + ( HHH * VAR0(1))
  41 :        FYNE = FYN + ( HHH * VAR0(1))
  42 : C
  43 :        XEPSE =  VAR0(1)
  44 : C
  45 :        XSIGT = XSIG0 + EYOUN * DEPST
  46 : C 
  47 :        IF (XSIGT . GT . FYPE) THEN
  48 : C
  49 :          XDEP = ((EYOUN/(EYOUN+HHH))*DEPST) 
  50 :      &         + ((XSIG0 - FYPE)/(EYOUN+HHH))
  51 :          XSIGF =  XSIG0 + (EYOUN*(DEPST - XDEP))
  52 :          XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  53 : C
  54 :        ELSEIF (XSIGT . LT . FYNE) THEN
  55 : C
  56 :          XDEP = ((EYOUN/(EYOUN+HHH))*DEPST) 
  57 :      &         + ((XSIG0 - FYNE)/(EYOUN+HHH))
  58 :          XSIGF =  XSIG0 + (EYOUN*(DEPST - XDEP))
  59 :          XTAN = (EYOUN*HHH)/(EYOUN+HHH)
  60 : C
  61 :        ELSE
  62 : C
  63 :          XSIGF = XSIGT
  64 :          XDEP = 0.D0        
  65 :          XTAN = EYOUN
  66 : C
  67 :        ENDIF
  68 : C
  69 : C
  70 : C
  71 :        SIGF = XSIGF * XSECT
  72 : C
  73 :        VARF(1) = XEPSE + XDEP
  74 :        VARF(2) = XTAN
  75 : C
  76 :        DEFP    = XEPSE
  77 : C
  78 :        RETURN
  79 :        END
  80 :  
  81 :  
  82 :  
  83 :  
  84 :