Esope source files

   1 : C FIBETA    SOURCE    CHAT      11/03/16    21:22:27     6902           
   2 : C FIBETA      SOURCE    AM        96/12/24    21:19:32     2448           
   3 : C======================================================================
   4 : C       FIBETA  - D. COMBESCURE et P. PEGON - ELSA- 1996
   5 : C======================================================================
   6 : C
   7 : C                    MODELE GLOBAL D'ETAGE 
   8 : C           (Sur des elements de poutre TIMO - Fibre
   9 : C                 (Contrainte de cisaillement/Cisail.)
  10 : C
  11 :       SUBROUTINE FIBETA(XMAT,XCAR,SIG0,VAR0,SIGF,VARF,DEPS,
  12 :      . WRK2,NCURVP,NCURVN,KERRE)
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8(A-H,O-Z)
  15 : C
  16 : C=======================================================================
  17 : C             CETTE ROUTINE EST APPELE DANS ECO100
  18 : C
  19 : C
  20 : C WRK0      = Segment materiaux
  21 : C WRK1      = Segment contraintes
  22 : C WRK2      = Segment courbe
  23 : C NCURVP    = Longueur courbe sens positif
  24 : C NCURVN    = Longueur courbe sens negatif
  25 : C
  26 : C======================================================================
  27 : C XCAR      = Caracteristique de la section 
  28 : C DEPS     = Increment de deformation axiale
  29 : C SIG0      = Contrainte initiale
  30 : C VAR0      = Variables internes initiales
  31 : C SIGF      = Contrainte finale
  32 : C VARF      = Variables internes finales
  33 : C DEFP      = Deformation plastique
  34 : C
  35 : C=======================================================================
  36 : C
  37 :       PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
  38 :        DIMENSION XMAT(*),XCAR(*),SIG0(*),VAR0(*),SIGF(*),VARF(*),DEPS(*)
  39 : C
  40 : C=======================================================================
  41 : C      VARIABLES ET SEGMENTS NECESSAIRES
  42 : C=========================================================================
  43 : -INC CCOPTIO
  44 : C  Segment de la courbe
  45 :        SEGMENT WRK2
  46 :          REAL*8 TRAC(LTRAC)
  47 :        ENDSEGMENT
  48 : C
  49 : C
  50 : C  Lecture materiau
  51 : C
  52 :        YOUN = XMAT(1)
  53 :        XNU = XMAT(2)
  54 :        XDELAP = XMAT(5)
  55 :        XDELAN = (-1.D0)*XMAT(6)
  56 :        XDMAXP = XMAT(7)
  57 :        XDMAXN = XMAT(8)
  58 :        XBETA  = XMAT(9)
  59 :        XALPH  = XMAT(10)
  60 :        XTETA  = XMAT(11)
  61 : C
  62 : C
  63 :        IF (IDIM.EQ.3) THEN
  64 :         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(2)
  65 : C
  66 :         XDDEP = DEPS(3)
  67 :         XFOR0 = SIG0(3)
  68 : C
  69 :        ELSE
  70 : C
  71 :         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(1)
  72 : C
  73 :         XDDEP = DEPS(2)
  74 :         XFOR0 = SIG0(2)
  75 :        ENDIF
  76 : C
  77 : C Appel de la loi 
  78 : C
  78 :         CALL ETAGEL(XDDEP,XFOR0,XFORF,VAR0,VARF,SECZ,
  79 :      & XDELAP,XDELAN,XDMAXP,XDMAXN,XBETA,XALPH,XTETA,
  80 :      & WRK2,NCURVP,NCURVN,KERRE)
  81 : C
  82 : C==================================================================
  83 : C
  84 : C  On remplit les tableaux avant de sortir
  85 : C
  86 : C==================================================================
  87 :        IF (IDIM.EQ.3) THEN
  88 :         SIGF(1) = XZER 
  89 :         SIGF(2) = XZER
  90 :         SIGF(3) = XFORF
  91 :         SIGF(4) = XZER
  92 :         SIGF(5) = XZER
  93 :         SIGF(6) = XZER
  94 :        ELSE
  95 :         SIGF(1) = XZER 
  96 :         SIGF(2) = XFORF
  97 :         SIGF(3) = XZER
  98 :        ENDIF
  99 : C
 100 : C=====================================================
 101 : C            FIN DE LA ROUTINE DU MODELE
 102 : C======================================================
 103 :       RETURN
 104 : 
 105 :       END
 106 :  
 107 :  
 108 :  
 109 :