Esope source files

   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 :