Esope source files

   1 : C FIBFRA    SOURCE    CHAT      12/04/06    21:15:11     7348           
   2 : C FIBFRA    SOURCE    DC        98/01/26            
   3 :       SUBROUTINE FIBFRA(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(4),SIGF(3),VARF(4)
   7 : C----------------------------------------------------------------------
   8 : C
   9 : C     COMPORTEMENT DE FIBRE FRAGILE EN COMPRESSION ET TRACTION
  10 : C     MODELE PLASTIQUE POUR LE RENFORCEMENT TFC
  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      XMAT( 3) RHO
  18 : C      XMAT( 4) ALPHA
  19 : C      XMAT( 5) FC
  20 : C      XMAT( 6) FC_R
  21 : C      XMAT( 7) STC
  22 : C      XMAT( 8) FT
  23 : C      XMAT( 9) FT_R
  24 : C      XMAT(10) STT
  25 : C
  26 : C      DEPS(1) EPXX    SIG0(1) SIGF(1) SMXX
  27 : C      DEPS(2) GAXY    SIG0(2) SIGF(2) SMXY
  28 : C      DEPS(3) GAXZ    SIG0(3) SIGF(3) SMXZ
  29 : C                    
  30 : C                      VAR0(1) VARF(1) EPSO ALLONGEMENT TOTAL
  31 : C                      VAR0(2) VARF(2) DAMC ENDOMMAGEMENT EN COMP.
  32 : C                      VAR0(3) VARF(3) DAMT ENDOMMAGEMENT EN TRACT.
  33 : C
  34 : C----------------------------------------------------------------------
  35 :        GAMMA=XMAT(1)/2.D0/(XMAT(2)+1.D0)
  36 :        EYOUN = XMAT( 1)
  37 :        FC = (-1.D0)*XMAT(5)
  38 :        FC_R = (-1.D0)*XMAT(6)
  39 :        STC = (-1.D0)*XMAT(7)
  40 :        XEPSOC = FC/EYOUN
  41 :        FT = XMAT(8)
  42 :        FT_R = XMAT(9)
  43 :        STT = XMAT(10)
  44 :        XEPSOT = FT/EYOUN
  45 : C
  46 :        XDAMC =  VAR0(2) 
  47 :        XDAMT =  VAR0(3) 
  48 : C
  49 :        XEPSO = VAR0(1) + DEPS(1)     
  50 : C
  51 : C        CALCUL DE LA VARIABLE D'ENDOMMAGEMENT
  52 : C
  53 :        IF (XEPSO.GT. 0.D0) THEN
  54 : C
  55 : C         CAS DE LA TRACTION
  56 : C
  57 :           XDAMT0 = 1.0D0  - ( (XEPSOT/XEPSO)* ( (FT_R / FT) +
  58 :      & ((1.D0 - (FT_R / FT))*EXP((XEPSOT - XEPSO)/ STT) )) )
  59 : C  
  60 :          IF (XDAMT0. GT . XDAMT)  THEN
  61 :            XDAMT = XDAMT0
  62 :          ENDIF
  63 : C
  64 :          XDAM = XDAMT
  65 : C
  66 :        ELSEIF (XEPSO.LT. 0.D0) THEN
  67 : C
  68 : C        CAS DE LA COMPRESSION
  69 : C
  70 :           XDAMC0 = 1.0D0  - ( (XEPSOC/XEPSO)* ( (FC_R /FC) +
  71 :      & ((1.D0 - (FC_R /FC))*EXP((XEPSOC - XEPSO)/STC) )) )
  72 : C  
  73 :          IF (XDAMC0. GT . XDAMC)  THEN
  74 :            XDAMC = XDAMC0
  75 :          ENDIF        
  76 : C
  77 :          XDAM = XDAMC
  78 : C
  79 :        ELSE
  80 :          XDAM= 0.D0
  81 :        ENDIF
  82 : C
  83 : C
  84 :        SIGF(1)= EYOUN* (1.D0 - XDAM) * XEPSO
  85 : C
  86 :        SIGF(2) = SIG0(2) + GAMMA   * DEPS(2)
  87 :        SIGF(3) = SIG0(3) + GAMMA   * DEPS(3)
  88 : C
  89 :        VARF(1) = XEPSO
  90 :        VARF(2) = XDAMC
  91 :        VARF(3) = XDAMT      
  92 :        VARF(4) = (1.D0 - XDAM) * EYOUN 
  93 : C
  94 :        RETURN
  95 :        END
  96 :  
  97 :  
  98 :  
  99 :