Esope source files

   1 : C FIBPAA    SOURCE    CHAT      12/04/06    21:15:12     7348           
   2 : C FIBPAR    SOURCE    DC        98/08/31    21:16:01     3283           
   3 : C FIBFRA    SOURCE    DC        98/01/26            
   4 :       SUBROUTINE FIBPAA(XMAT,DEPS,SIG0,VAR0,SIGF,VARF)
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION XMAT(15),DEPS(3),SIG0(3),VAR0(8),SIGF(3),VARF(8)
   8 : C----------------------------------------------------------------------
   9 : C
  10 : C      DRIVER POUR LE MODELE A FIBRE 
  11 : C      MODELE ANCRAGE AVEC ACIER PARFAIT
  12 : C----------------------------------------------------------------------
  13 : C     Didier COMBESCURE - EMSI - MAI 2000 
  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      XMAT( 7) KS
  23 : C      XMAT( 8) S1T
  24 : C      XMAT( 9) S2T 
  25 : C      XMAT(10) S3T
  26 : C      XMAT(11) T1T 
  27 : C      XMAT(12) T3T
  28 : C      XMAT(13) ALFA 
  29 : C      XMAT(14) LANC
  30 : C      XMAT(15) SECT
  31 : C
  32 : C      DEPS(1) EPXX    SIG0(1) SIGF(1) SMXX
  33 : C      DEPS(2) GAXY    SIG0(2) SIGF(2) SMXY
  34 : C      DEPS(3) GAXZ    SIG0(3) SIGF(3) SMXZ
  35 : C                    
  36 : C                      VAR0(1) VARF(1) EPSE DEFORMATION PLASTIQUE
  37 : C                      VAR0( 3) VARF( 3) EPPP
  38 : C                      VAR0( 4) VARF( 4) EPPM
  39 : C                      VAR0( 5) VARF( 5) EPSOA
  40 : C                      VAR0( 6) VARF( 6) KSTAN
  41 : C                      VAR0( 7) VARF( 7) NITER
  42 : C                      VAR0( 8) VARF( 8) XLAMBDA
  43 : 
  44 : C
  45 : C----------------------------------------------------------------------
  46 : C
  47 : C
  48 :        XSIG0 = SIG0(1)
  49 :        DEPST = DEPS(1)
  50 :        XSECT = XMAT(15)
  51 :        GAMMA=XMAT(1)/2.D0/(XMAT(2)+1.D0)
  52 : C
  52 :        CALL PARANC(XMAT,XSECT,DEPST,XSIG0,VAR0,XSIGF,VARF,DEFP)   
  53 : C
  54 : C
  55 :        SIGF(1) = XSIGF
  56 : C
  57 :        SIGF(2) = SIG0(2) + GAMMA   * DEPS(2)
  58 :        SIGF(3) = SIG0(3) + GAMMA   * DEPS(3)
  59 : C
  60 :        RETURN
  61 :        END
  62 :  
  63 :  
  64 :  
  65 :  
  66 :