Esope source files

   1 : C BILINE    SOURCE    CHAT      05/01/12    21:40:08     5004           
   2 :       SUBROUTINE BILINE(DELTAD,FY,STIFE,STIFP,RIG,FE0,FEP0,F,DELDPL)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C SUBROUTINE FOR ELASTIC PLASTIC (WITH HARDENING) RELATIONSHIP.
   6 : C
   7 : C
   8 : C A PERFECT E.P. SCALAR MODEL WITH INITIAL STIFNESS STIFEQ AND
   9 : C YIELDING FYEQ IS ADDED IN PARALEL TO AN ELASTIC SCALAR MODEL
  10 : C OF STIFNESS STIFP.
  11 : C
  12 : C INPUT
  13 : C
  14 : C    DELTAD  DISPLACEMENT INCREMENT
  15 : C    FY      YIELDING FORCE (FY=STIFE*FYEQ/STIFEQ)
  16 : C    STIFE   ELASTIC STIFNESS (STIFE=STIFEQ+STIFP)
  17 : C    STIFP   PLASTIC STIFNESS
  18 : C    FE0     INITIAL FORCE IN THE ELASTIC MODEL
  19 : C    FEP0    INITIAL FORCE IN THE E.P. MODEL
  20 : C
  21 : C OUTPUT
  22 : C
  23 : C    RIG     FINAL STIFFNESS
  24 : C    F       FINAL FORCE
  25 : C    FE0     FINAL FORCE IN THE ELASTIC MODEL
  26 : C    FEP0    FINAL FORCE IN THE ELASTIC MODEL
  27 : C    DELDPL  PLASTIC DEFORMATION INCREMENT
  28 : C
  29 : C WARNING FE0 AND FEP0 SHOULD BE INITIALIZED !!
  30 : C
  31 : C-----------------------------------------------------------------------
  32 : C VERY OLD (CRAZY ?) STUFF SLIGHTLY MODIFIED FOR DELDPL
  33 : C AND TAKEN AS AN EXEMPLE OF "EXISTING" MODEL
  34 : C-----------------------------------------------------------------------
  35 :       IMPLICIT INTEGER(I-N)
  36 :       IMPLICIT REAL*8(A-H,O-Z)
  37 : C
  38 : C  COMPUTATION OF THE STIFNESS AND YIELDING OF THE E.P. MODEL
  39 : C
  40 :       STIFEQ=STIFE-STIFP
  41 :       FYEQ=FY*STIFEQ/STIFE
  42 : C
  43 : C  ELASTIC PREDICTION OF THE FORCE IN THE E.P. MODEL ...
  44 : C
  45 :       FEP=FEP0+STIFEQ*DELTAD
  46 : C
  47 : C  ... AND POSSIBLE RETURN
  48 : C
  49 :       IF (ABS(FEP).GT.FYEQ)THEN
  50 :         FEP=SIGN(FYEQ,FEP)
  51 :         RIG=STIFP
  52 :         DELDPL=DELTAD-(FEP-FEP0)/STIFEQ
  53 :       ELSE
  54 :         RIG=STIFE
  55 :         DELDPL=0.D0
  56 :       ENDIF
  57 :       FEP0=FEP
  58 : C
  59 : C  FORCE IN THE ELASTIC MODEL
  60 : C
  61 :       FE=FE0+STIFP*DELTAD
  62 :       FE0=FE
  63 : C
  64 : C  TOTAL FORCE
  65 : C
  66 :       F=FE+FEP
  67 :       RETURN
  68 :       END
  69 :  
  70 :