Esope source files

   1 : C SJONC2    SOURCE    CHAT      05/01/13    03:20:33     5004           
   2 :       SUBROUTINE SJONC2(SIGI,DEPST,VARI,XMAT,TRAS,NTRAS,
   3 :      >                                     SIGF,VARF,DEFP,KERRE)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C     ROUTINE DE CALCUL DE L'ECOULEMENT DANS LE JOINT
   7 : C
   8 : C     WORK ONLY IN 2-D PLANE STATE
   9 : C
  10 : C  INPUT
  11 : C     SIGI    INITIAL STRESS
  12 : C     DEPST   INCREMENT OF STRAIN
  13 : C     XMAT    MATERIAL PROPERTIES
  14 : C     TRAS    SHEAR CURVE
  15 : C     VARI    INITIAL INTERNAL VARIABLES
  16 : C
  17 : C  OUTPUT
  18 : C     SIGF    FINAL STRESS
  19 : C     VARF    FINAL INTERNAL VARIABLES
  20 : C     DEFP    INCREMENT OF PLASTIQUE STRAIN
  21 : C
  22 : C-----------------------------------------------------------------------
  23 : C  P.PEGON   SEPTEMBER 2001
  24 : C-----------------------------------------------------------------------
  25 :       IMPLICIT INTEGER(I-N)
  26 :       IMPLICIT REAL*8(A-H,O-Z)
  27 : *
  28 :       DIMENSION SIGI(*),DEPST(*),VARI(*),
  29 :      .          SIGF(*),VARF(*) ,DEFP(*),
  30 :      .          XMAT(*),TRAS(2,NTRAS)
  31 : C
  32 : C-----------------------------------------------------------------------
  33 : C  MATERIAL PARAMETERS
  34 : C  ===================
  35 : C
  36 : C      KS    = XMAT(1)   Shear modulus
  37 : C      TRAS  = X/Y       PLASTIC Shear curve
  38 : C      BETA  = XMAT(5)   Unloading Coefficient
  39 : C-----------------------------------------------------------------------
  40 : C  STRESS
  41 : C  ======
  42 : C
  43 : C      SMS1.  = SIG.(1)
  44 : C      SMS2   = SIG.(2)
  45 : C-----------------------------------------------------------------------
  46 : C  STRAIN
  47 : C  ======
  48 : C
  49 : C      DRS1   = EPS(1)
  50 : C      DRS2   = EPS(2)
  51 : C-----------------------------------------------------------------------
  52 : C  INTERNAL VARIABLES
  53 : C  ==================
  54 : C
  55 : C
  56 : C      EPSS   = VAR.(1) PLASTIC SHEAR COEFFICIENT
  57 : C
  58 : C
  59 : C-----------------------------------------------------------------------
  60 : C
  61 : C        CURRENT POSITION ON THE LOADING CURVE
  62 : C
  63 :       EPSS=VARI(1)
  63 :       CALL YOFXCU(EPSS,TRAS,NTRAS,SMSNP,DSMSNP,KERRE)
  64 :       IF (KERRE.NE.0) RETURN
  65 : C
  66 : C        VALUE OF THE LOADING/UNLOADING MODULUS
  67 : C
  68 :       DMODU=SMSNP/(EPSS+(TRAS(2,1)/XMAT(1)))
  69 :       DMODU=(1-XMAT(5))*XMAT(1)+DMODU*XMAT(5)
  70 : C
  71 : C        ELASTIC STRESS PREDICTION
  72 : C
  73 :       SIGF(1)=SIGI(1)+DMODU*DEPST(1)
  74 :       SIGF(2)=SIGI(2)+DMODU*DEPST(2)
  75 :       SMSN=SQRT(SIGF(1)**2+SIGF(2)**2)
  76 : C
  77 : C        TEST ON THE FLOW
  78 : C
  79 :       IF (SMSN.GT.SMSNP) THEN
  80 : C
  81 : C        -> YIELDING IN SHEAR ...
  82 : C
  83 :         EPSS= EPSS+(SMSN-SMSNP)/DMODU
  83 :         CALL YOFXCU(EPSS,TRAS,NTRAS,SMSNP,DSMSNP,KERRE)
  84 :         IF (KERRE.NE.0) RETURN
  85 : C
  86 :         SIGF(1)=SIGF(1)*SMSNP/SMSN
  87 :         SIGF(2)=SIGF(2)*SMSNP/SMSN
  88 : C
  89 :       ENDIF 
  90 : C
  91 :       VARF(1)= EPSS
  92 :       DEFP(1) = VARF(1) - VARI(1)
  93 : C
  94 :       RETURN
  95 :       END
  96 :  
  97 :