Esope source files

   1 : C AREVA1    SOURCE    STRU      11/01/21    21:15:01     6855           
   2 :       SUBROUTINE AREVA1(SIG0,DEPST,VAR0,XMAT,XCAR,SIGF,VARF,DEFP)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    DRIVER OF THE BILINEAR PLASTIC BENDING MODEL
   6 : C
   7 : C  INPUT:
   8 : C
   9 : C    SIG0             INITIAL GENERALIZED STRESS
  10 : C    DEPST            INCREMENT OF GENERALIZED STRAIN
  11 : C    XMAT             MATERIAL CHARACTERISTICS
  12 : C    XCAR             GEOMETRIC CHARACTERISTICS
  13 : C    VAR0             INITIAL INTERNAL VARIABLE
  14 : C
  15 : C  OUTPUT:
  16 : C
  17 : C    SIGF             FINAL GENERALIZED STRESS
  18 : C    VARF             FINAL INTERNAL VARIABLE
  19 : C    DEFP             INCREMENT OF PLASTIC DEFORMATION
  20 : C
  21 : C  GENERALIZED STRESS (.=0 OR F):
  22 : C
  23 : C    MOMY = SIG.(5)
  24 : C
  25 : C  MATERIAL CHARACTERISTICS:
  26 : C
  27 : C    E    = XMAT(1)   ELASTIC YOUNG MODULUS
  28 : C    NU   = XMAT(2)   POISSON RATIO
  29 : C    ETAN = XMAT(5)   PLASTIC MODULUS
  30 : C    YIEF = XMAT(6)   YIELDING FORCE
  31 : C
  32 : C  GEOMETRIC CHARACTERISTICS:
  33 : C
  34 : C    INRX = XCAR(1)   X-INERTIAL TORSION
  35 : C    INRY = XCAR(2)   Y-INERTIAL BENDING        <----- 
  36 : C    INRZ = XCAR(3)   Z-INERTIAL BENDING
  37 : C    SECT = XCAR(4)   SECTION
  38 : C    SRDY = XCAR(5)   TRANSVERSE SHEAR'S REDUCED SECTION
  39 : C    SRDZ = XCAR(6)   TRANSVERSE SHEAR'S REDUCED SECTION
  40 : C
  41 : C  INTERNAL VARIABLE (.=0 OR F):
  42 : C
  43 : C    EPSE = VAR.(1)   CUMULATED PLASTIC AXIAL STRAIN
  44 : C    EITA = VAR.(2)   TANGENT X-AXIAL STIFFNESS
  45 : C    MOEY = VAR.(3)   ELASTIC AXIAL FORCE
  46 : C    EPSX = VAR.(4)   AXIAL STRAIN
  47 : C
  48 :       IMPLICIT INTEGER(I-N)
  49 :       IMPLICIT REAL*8(A-H,O-Z)
  50 : -INC CCOPTIO
  51 : C
  52 :       DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),
  53 :      .          SIGF(*),VARF(*) ,DEFP(*)
  54 : C
  55 : C-----------------------------------------------------------------------
  56 : C
  57 : C  ELASTIC STRESS INCREMENT FOR ALL STRESS COMPONENT EXCEPT EFFX
  58 : C
  59 :       CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2))
  60 :       SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2)
  61 :       SIGF(3)=SIG0(3) + CISA*XCAR(6)*DEPST(3)
  62 :       SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4)
  63 :       SIGF(5)=SIG0(5) + XMAT(1)*XCAR(2)*DEPST(5)
  64 :       SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6)
  65 : C
  66 : C  ELASTIC AND PLASTIC AXIAL STIFFNESS
  67 : C
  68 :       EA     = XMAT(1)*XCAR(4)
  69 :       ETA    = XMAT(5)*XCAR(4)
  70 : C
  71 : C
  72 : C  AXIAL STRAIN INCREMENT
  73 : C
  74 :       DEPSX = DEPST(1)
  75 : C
  75 : C  CALL TO THE SCALAR BILINEAR PLASTIC MODEL COMPUTING MOMY
  76 : C
  77 :       FEP0 = SIG0(1) - VAR0(3)
  78 :       FE0  = VAR0(3)
  78 :       CALL BILINE(DEPSX,XMAT(6),EA, ETA, VARF(2),
  79 :      >          FE0,FEP0,SIGF(1),DEPSE)
  80 : C
  81 : C  FINAL INTERNAL VARIABLE
  82 : C
  83 :       VARF(1)=VAR0(1)+DEPSE
  84 :       VARF(3)=FE0
  85 :       VARF(4)=VAR0(4)+DEPSX
  86 : C
  87 : C  PLASTIC DEFORMATION INCREMENT
  88 : C
  89 :       DEFP(1) = DEPSE
  90 :       DEFP(2) = 0.D0
  91 :       DEFP(3) = 0.D0
  92 :       DEFP(4) = 0.D0
  93 :       DEFP(5) = 0.D0
  94 :       DEFP(6) = 0.D0
  95 :       RETURN
  96 :       END
  97 :  
  98 :  
  99 :