Esope source files

   1 : C BARSTE    SOURCE    CHAT      05/01/12    21:35:29     5004           
   2 :       SUBROUTINE BARSTE(XMAT,XCAR,DEPST,SIG0,VAR0,SIGF,VARF,DEFP)
   3 : *
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 :       DIMENSION XMAT(16),VAR0(14),VARF(14)
   7 :       DIMENSION EPSON(5),SIGMA(5),DEPST(*)
   8 :       CHARACTER*8 CMATE
   9 : C----------------------------------------------------------------------
  10 : C
  11 : C     COMPORTEMENT DE FIBRE DE FERRAILLE
  12 : C
  13 : C----------------------------------------------------------------------
  14 : C     Pierre Pegon (ISPRA) Juillet 1994 +- clone de FIBSTE
  15 : C----------------------------------------------------------------------
  16 : C
  17 : C      XMAT( 1) YOUNG
  18 : C      XMAT( 2) NU
  19 : C      XMAT( 3) RHO
  20 : C      XMAT( 4) ALPHA
  21 : C      XMAT( 5) MAX. ELAST. STRESS
  22 : C      XMAT( 6) ULTIMATE STRESS
  23 : C      XMAT( 7) HARDENING STRAIN
  24 : C      XMAT( 8) ULTIMATE STRAIN
  25 : C      XMAT( 9) Ro FACTOR
  26 : C      XMAT(10) b  FACTOR
  27 : C      XMAT(11) a1 FACTOR
  28 : C      XMAT(12) a2 FACTOR
  29 : C      XMAT(13) FACTOR L/D
  30 : C      XMAT(14) a6 FACTOR
  31 : C      XMAT(15) c FACTOR
  32 : C      XMAT(16) a FACTOR
  33 : C
  34 : C      SIG0 AND SIGF STRESS
  35 : C      DEPST,DEFP
  36 : C
  37 : C                      VAR0( 1) VARF( 1) FLAG
  38 : C                      VAR0( 2) VARF( 2) EPSO
  39 : C                      VAR0( 3) VARF( 3) R FACTOR
  40 : C                      VAR0( 4) VARF( 4) TANGENT MODULUS
  41 : C                      VAR0( 5) VARF( 5) LAST REV. STRAIN
  42 : C                      VAR0( 6) VARF( 6) INTERSEPTION POINT
  43 : C                      VAR0( 8) VARF( 7) MAX. REV. POST. STRAIN
  44 : C                      VAR0( 9) VARF( 8) MIN. REV. POST. STRAIN
  45 : C                      VAR0(10) VARF( 9) LAST REV. STRESS
  46 : C                      VAR0(11) VARF(10) INTERSEPTION POINT
  47 : C                      VAR0(13) VARF(11) MAX. REV. POST. STRESS
  48 : C                      VAR0(14) VARF(12) MIN. REV. POST. STRESS
  49 : C
  50 : C----------------------------------------------------------------------
  51 : *
  52 :        DEPS = DEPST(1)
  53 :        EYOUN = XMAT( 1)
  54 :        FSY__ = XMAT( 5)
  55 :        FSU__ = XMAT( 6)
  56 :        EPSH_ = XMAT( 7)
  57 :        EPSU_ = XMAT( 8)
  58 :        ROFAC = XMAT( 9)
  59 :        BFACT = XMAT(10)
  60 :        A1FAC = XMAT(11)
  61 :        A2FAC = XMAT(12)
  62 :        FACLD = XMAT(13)
  63 :        A6FAC = XMAT(14)
  64 :        CFACT = XMAT(15)
  65 :        AFACT = XMAT(16)
  66 : C
  67 :        IFLAG = NINT(VAR0(1))
  68 :        STRAI = VAR0(2) + DEPS
  69 :        RFACT = VAR0(3)
  70 :        TANGE = VAR0(4)
  71 : C
  72 :        DO IFACT = 1, 2
  73 :           EPSON(IFACT) = VAR0(4+IFACT)
  74 :           SIGMA(IFACT) = VAR0(8+IFACT)
  75 :        END DO
  76 : C
  77 :        EPSON(3) = VAR0(2)
  78 :        SIGMA(3) = SIG0 / XCAR
  79 : C
  80 :        DO IFACT = 4, 5
  81 :           EPSON(IFACT) = VAR0(3+IFACT)
  82 :           SIGMA(IFACT) = VAR0(7+IFACT)
  83 :        END DO
  84 : C
  84 :        CALL STEEL1
  85 :      .       ( EYOUN, FSY__, FSU__, EPSH_, EPSU_, ROFAC, BFACT,
  86 :      .         A1FAC, A2FAC, FACLD, A6FAC, CFACT, AFACT, IFLAG,
  87 :      .         STRAI, RFACT, EPSON, SIGMA, TANGE, STRES )
  88 : C
  89 :        VARF(1) = DBLE(IFLAG)
  90 :        VARF(2) = STRAI
  91 :        VARF(3) = RFACT
  92 :        VARF(4) = TANGE
  93 : C
  94 :        DO IFACT = 1, 2
  95 :           VARF(4+IFACT) = EPSON(IFACT)
  96 :           VARF(8+IFACT) = SIGMA(IFACT)
  97 :        END DO
  98 : C
  99 :        DO IFACT = 4, 5
 100 :           VARF(3+IFACT) = EPSON(IFACT)
 101 :           VARF(7+IFACT) = SIGMA(IFACT)
 102 :        END DO
 103 : C
 104 :        SIGF    = STRES * XCAR
 105 :        DEFP    = DEPS - ((SIGF - SIG0) / (EYOUN * XCAR))
 106 : C
 107 :        RETURN
 108 :        END
 109 :  
 110 :  
 111 :