Esope source files

   1 : C DJONL2    SOURCE    CHAT      05/01/12    22:52:31     5004           
   2 :       SUBROUTINE DJONL2(SIGI,DEPST,VARI,XMAT,SIGF,VARF,DEFP,KERRE)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C     ROUTINE DE CALCUL DE L'ECOULEMENT DANS LE JOINT
   6 : C
   7 : C     "DILATANT JOINT WITH COULOMBS FRICTION"
   8 : C     MODELE PROPOSED BY MF SNYMANN, WW BIRD AND JB MARTIN
   9 : C
  10 : C     WORK ONLY IN 2-D PLANE (STRAIN) STATE
  11 : C
  12 : C  INPUT
  13 : C     SIGI    INITIAL STRESS
  14 : C     DEPST   TRIAL INCREMENT OF ELASTIC STRAIN
  15 : C     XMAT    MATERIAL PROPERTIES
  16 : C     VARI    INITIAL INTERNAL VARIABLES
  17 : C
  18 : C  OUTPUT
  19 : C     SIGF    FINAL STRESS
  20 : C     VARF    FINAL INTERNAL VARIABLES
  21 : C     DEFP    INCREMENT OF PLASTIQUE STRAIN
  22 : C
  23 : C-----------------------------------------------------------------------
  24 : C  P.PEGON   OCTOBER 93
  25 : C-----------------------------------------------------------------------
  26 :       IMPLICIT INTEGER(I-N)
  27 :       IMPLICIT REAL*8(A-H,O-Z)
  28 : -INC CCOPTIO
  29 : -INC CCREEL
  30 : *
  31 :       DIMENSION SIGI(*),DEPST(*),VARI(*),XMAT(*),
  32 :      .          SIGF(*),VARF(*) ,DEFP(*)
  33 :       REAL*8 KS,KN,MU
  34 : C
  35 : C  QUICK QUIT IN CASE OF DIMENSION ERROR
  36 : C
  37 :       IF (IFOUR.NE.-3.AND.IFOUR.NE.-2.AND.IFOUR.NE.-1)THEN
  38 :         KERRE=99
  39 :         RETURN
  40 :       ELSE
  41 :         KERRE=0
  42 :       ENDIF
  43 : C
  44 : C-----------------------------------------------------------------------
  45 : C  MATERIAL PARAMETERS
  46 : C  ===================
  47 : C
  48 : C      KS    = XMAT(1)   Shear modulus
  49 : C      KN    = XMAT(2)   Traction modulus
  50 : C      PHI   = XMAT(5)   Friction angle
  51 : C      MU    = XMAT(6)   Dilatation angle
  52 : C      FRTC  = XMAT(7)   Traction maximum (tan(PHI)*FRTC = cohesion)
  53 : C-----------------------------------------------------------------------
  54 : C  STRESS
  55 : C  ======
  56 : C
  57 : C      SMSN.  = SIG.(1)
  58 : C      SMN.   = SIG.(2)
  59 : C-----------------------------------------------------------------------
  60 : C  STRAIN
  61 : C  ======
  62 : C
  63 : C      DRSN   = EPS(1)
  64 : C      DRN    = EPS(2)
  65 : C-----------------------------------------------------------------------
  66 : C  INTERNAL VARIABLES
  67 : C  ==================
  68 : C
  69 : C      DS0.   = VAR.(1) POSITION OF THE VERTEX      (RELATIVE SHEAR)
  70 : C      DN0.   = VAR.(2) POSITION OF THE VERTEX      (EXTENSION)
  71 : C      DRSN.  = VAR.(3) CURRENT DEFORMATION         (RELATIVE SHEAR)
  72 : C      DRN.   = VAR.(4) CURRENT DEFORMATION         (EXTENSION)
  73 : C      DPSN.  = VAR.(5) CURRENT PLASTIC DEFORMATION (RELATIVE SHEAR)
  74 : C      DPN.   = VAR.(6) CURRENT PLASTIC DEFORMATION (EXTENSION)
  75 : C
  76 : C      HC11   = VAR.(7)
  77 : C      HC12   = VAR.(8) CURRENT VALUES OF THE 
  78 : C      HC21   = VAR.(9) CONSISTENT HOOK MATRIX
  79 : C      HC22   = VAR.(10)
  80 : C-----------------------------------------------------------------------
  81 :        KS     = XMAT(1)
  82 :        KN     = XMAT(2)
  83 :        PSI    = ATAN(KN/KS*TAN(XMAT(5)*XPI/180))
  84 :        MU     = XMAT(6)*XPI/180
  85 :        FTRAC  = XMAT(7)
  86 :        DTRAC  = XMAT(7)/KN
  87 : C
  88 :        TANPSI = TAN(PSI)
  89 :        TANMU  = TAN(MU)
  90 : C
  91 : C      TOTAL "STRAIN" FROM THE STRESS INCREMENT
  92 : C
  93 :        DS0I   = VARI(1)
  94 :        DN0I   = VARI(2)
  95 :        DRSNF  = VARI(3) + DEPST(1)
  96 :        DRNF   = VARI(4) + DEPST(2)
  97 :        VARF(3)= DRSNF
  98 :        VARF(4)= DRNF
  99 : C
 100 : C      FIRST REGION
 101 : C
 102 :        IF ((DRNF-DTRAC).GE.(DN0I+ABS(DRSNF-DS0I)*TANMU))THEN
 103 :          VARF(5) = DRSNF
 104 :          VARF(6) = DRNF
 105 :          VARF(1) = DRSNF
 106 :          VARF(2) = DN0I+ABS(DRSNF-DS0I)*TANMU
 107 :          SIGF(1) = 0.D0
 108 :          SIGF(2) = FTRAC
 109 : C
 110 :          VARF(7) = 0.D0
 111 :          VARF(8) = 0.D0
 112 :          VARF(9) = 0.D0
 113 :          VARF(10)= 0.D0
 114 : C
 115 : C      SECOND REGION (ELASTIC CASE)
 116 : C
 117 :        ELSEIF((DRNF-DTRAC).LE.(DN0I-ABS(DRSNF-DS0I)/TANPSI))THEN
 118 :          VARF(5) = DS0I
 119 :          VARF(6) = DN0I
 120 :          VARF(1) = DS0I
 121 :          VARF(2) = DN0I
 122 :          SIGF(1) = KS * (DRSNF-DS0I)
 123 :          SIGF(2) = KN * (DRNF -DN0I)
 124 : C
 125 :          VARF(7) = KS
 126 :          VARF(8) = 0.D0
 127 :          VARF(9) = 0.D0
 128 :          VARF(10)= KN
 129 : C
 130 : C      THIRD REGION
 131 : C
 132 :        ELSE
 133 :          A=( (DRNF-DTRAC-DN0I)*SIN(PSI)+ABS(DRSNF-DS0I)*COS(PSI))
 134 :      >                                 /COS(PSI-MU)
 135 :          B=(-(DRNF-DTRAC-DN0I)*COS(MU )+ABS(DRSNF-DS0I)*SIN(MU ))
 136 :      >                                 /COS(PSI-MU)
 137 :          VARF(5) = DS0I + A * COS(MU) * SIGN(1.D0,DRSNF-DS0I)
 138 :          VARF(6) = DN0I + A * SIN(MU)
 139 :          VARF(1) = VARF(5)
 140 :          VARF(2) = VARF(6)
 141 :          SIGF(1) = KS * B * SIN(PSI) * SIGN(1.D0,DRSNF-DS0I)
 142 :          SIGF(2) = KN * B * COS(PSI) * (-1) + FTRAC
 143 : C
 144 :          COHOOK  = KN/(1+TANPSI*TANMU)
 145 :          TANPHI  = TANPSI*KS/KN
 146 :          VARF(7) = COHOOK*TANPHI*TANMU
 147 :          VARF(8) =-COHOOK*SIGN(1.D0,DRSNF-DS0I)*TANPHI
 148 :          VARF(9) =-COHOOK*SIGN(1.D0,DRSNF-DS0I)*TANMU
 149 :          VARF(10)= COHOOK
 150 : C
 151 :        ENDIF
 152 : C
 153 : C      INCREMENT OF PLASTIC STRAIN
 154 : C
 155 :        DEFP(1) = VARF(5) - VARI(5)
 156 :        DEFP(2) = VARF(6) - VARI(6)
 157 : C
 158 :        RETURN
 159 :        END
 160 :  
 161 :  
 162 :  
 163 :