Esope source files

   1 : C BILIFO    SOURCE    CHAT      05/01/12    21:40:04     5004           
   2 :       SUBROUTINE BILIFO(SIG0,DEPST,VAR0,XMAT,XCAR,SIGF,VARF,DEFP)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    DRIVER OF THE BILINEAR PLASTIC TRANSVERSE SHEAR 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    EFFZ = SIG.(3)
  24 : C
  25 : C  MATERIAL CHARACTERISTICS:
  26 : C
  27 : C    E    = XMAT(1)   ELASTIC YOUNG MODULUS
  28 : C    NU   = XMAT(2)   POISSON'S RATIO
  29 : C    EAYI = XMAT(5)   PLASTIC MODULUS
  30 : C    YEFF = XMAT(6)   YIELDING TRANSVERSE SHEAR 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    SECZ = XCAR(6)   Z-TRANSVERSE SHEAR SECTION  <------
  40 : C
  41 : C  INTERNAL VARIABLE (.=0 OR F):
  42 : C
  43 : C    EPSE = VAR.(1)   CUMULATED PLASTIC TRANSVERSE SHEAR STRAIN
  44 : C    GATA = VAR.(2)   TANGENT Z-SHEAR STIFFNESS
  45 : C    EFEZ = VAR.(3)   ELASTIC SHEAR FORCE
  46 : C    CISZ = VAR.(4)   TRANSVERSE SHEAR STRAIN
  47 : C
  48 : C-----------------------------------------------------------------------
  49 : C  P. PEGON      26/11/92 renewed new ECOU 24/11/95
  50 : C-----------------------------------------------------------------------
  51 :       IMPLICIT INTEGER(I-N)
  52 :       IMPLICIT REAL*8(A-H,O-Z)
  53 : -INC CCOPTIO
  54 : C
  55 :       DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),
  56 :      .          SIGF(*),VARF(*) ,DEFP(*)
  57 : C
  58 : C-----------------------------------------------------------------------
  59 : C
  60 : C  ELASTIC STRESS INCREMENT FOR ALL STRESS COMPONENT EXCEPT MOMY
  61 : C
  62 :       CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2))
  63 : C     SIGF(1)=SIG0(1) + DSIGT(1)
  64 :       SIGF(1)=SIG0(1) + XMAT(1)*XCAR(4)*DEPST(1)
  65 : C     SIGF(2)=SIG0(2) + DSIGT(2)
  66 :       SIGF(2)=SIG0(2) + CISA*XCAR(5)*DEPST(2)
  67 : C     SIGF(4)=SIG0(4) + DSIGT(4)
  68 :       SIGF(4)=SIG0(4) + CISA*XCAR(1)*DEPST(4)
  69 : C     SIGF(5)=SIG0(5) + DSIGT(5)
  70 :       SIGF(5)=SIG0(5) + XMAT(1)*XCAR(2)*DEPST(5)
  71 : C     SIGF(6)=SIG0(6) + DSIGT(6)
  72 :       SIGF(6)=SIG0(6) + XMAT(1)*XCAR(3)*DEPST(6)
  73 : C
  74 : C  ELASTIC AND PLASTIC SHEAR STIFNESS
  75 : C
  76 :       XMODG=0.5D0 / (1.D0 + XMAT(2))
  77 :       GAZ    = XMAT(1)*XCAR(6)*XMODG
  78 :       GYAZ   = XMAT(5)*XCAR(6)*XMODG
  79 : C
  80 : C  Z-TRANSVERSE SHEAR INCREMENT
  81 : C
  82 : C     DCISAZ = DSIGT(3) / GAZ
  83 :       DCISAZ = DEPST(3)
  84 : C
  84 : C  CALL TO THE SCALAR BILINEAR PLASTIC MODEL COMPUTING MOMY
  85 : C
  86 :       FEP0 = SIG0(3) - VAR0(3)
  87 :       FE0  = VAR0(3)
  87 :       CALL BILINE(DCISAZ,XMAT(6),GAZ, GYAZ, VARF(2),
  88 :      >          FE0,FEP0,SIGF(3),DEPSE)
  89 : C
  90 : C  FINAL INTERNAL VARIABLE
  91 : C
  92 :       VARF(1)=VAR0(1)+DEPSE
  93 :       VARF(3)=FE0
  94 :       VARF(4)=VAR0(4)+DCISAZ
  95 : C
  96 : C  PLASTIC DEFORMATION INCREMENT
  97 : C
  98 :       DEFP(1) = 0.D0
  99 :       DEFP(2) = 0.D0
 100 :       DEFP(3) = DEPSE
 101 :       DEFP(4) = 0.D0
 102 :       DEFP(5) = 0.D0
 103 :       DEFP(6) = 0.D0
 104 :       RETURN
 105 :       END
 106 :  
 107 :