Esope source files

   1 : C ANCREL    SOURCE    CHAT      05/01/12    21:22:18     5004           
   2 : C      
   3 :       SUBROUTINE ANCREL(SIGI,DEPST,VARI,XMAT,
   4 :      .                       SIGF,VARF,DEFP,KERRE)
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION XMAT(*),VARI(*),VARF(*),SIGI(*)
   8 :      . ,DEPST(*),DEFP(*),SIGF(*)
   9 : C----------------------------------------------------------------------
  10 : C
  11 : C     LOI D'ANCRAGE ELIGEHAUSEN SUR UN ELEMENT DE JOINT 2D
  12 : C----------------------------------------------------------------------
  13 : C     Didier COMBESCURE - EMSI - MAI 2000
  14 : C----------------------------------------------------------------------
  15 : C
  16 : C      XMAT( 1) XKS
  17 : C      XMAT( 2) XKN
  18 : C      XMAT( 3) RHO
  19 : C      XMAT( 4) ALPHA
  20 : C      XMAT( 5) S1T
  21 : C      XMAT( 6) S2T 
  22 : C      XMAT( 7) S3T
  23 : C      XMAT( 8) T1T 
  24 : C      XMAT( 9) T3T
  25 : C      XMAT(10) ALFA 
  26 : C      XMAT(11) PERI 
  27 : C
  28 : C  STRESS
  29 : C  ======
  30 : C
  31 : C      SMSN.  = SIG.(1)
  32 : C      SMN.   = SIG.(2)
  33 : C-----------------------------------------------------------------------
  34 : C  STRAIN
  35 : C  ======
  36 : C
  37 : C      DRSN   = EPS(1)
  38 : C      DRN    = EPS(2)
  39 : C-----------------------------------------------------------------------
  40 : C  INTERNAL VARIABLES
  41 : C  ==================
  42 : C                    
  43 : C      EPPP = VAR0(1) VARF(1) EPPP GLISSEMENT PLASTIQUE POSITIF
  44 : C      EPPM = VAR0(2) VARF(2) EPPM GLISSEMENT PLASTIQUE NEGATIF
  45 : C      EPSO = VAR0(3) VARF(3) EPSO GLISSEMENT TOTAL 
  46 : C
  47 : C  MATRICE TANGENTE
  48 : C
  49 : C      KSTAN   = VARF(4)
  50 : C----------------------------------------------------------------------
  51 : C
  52 :        XKS  = XMAT(1)
  53 :        XKN  = XMAT(2)
  54 :        S1T = XMAT(5)
  55 :        S2T = XMAT(6)
  56 :        S3T = XMAT(7)
  57 :        T1T = XMAT(8)
  58 :        T3T = XMAT(9)
  59 :        ALFA = XMAT(10)
  60 :        PERI = XMAT(11) 
  61 : C
  62 :        G12 = XKS/PERI 
  63 : C
  64 : C  EFFORT NORMAL (ELASTIQUE)
  65 : C
  66 :         SIGF(2) = SIGI(2)+(XKN*DEPST(2))     
  67 : C
  68 : C  CISAILLEMENT - GLISSEMENT
  69 : C
  70 :        EPPP   = VARI(1)
  71 :        EPPM   = VARI(2)
  72 : C
  73 :        EPSO = VARI(3) + DEPST(1)
  74 : C
  74 :        CALL ELIGEH(EPSO,G12,S1T,S2T,S3T,T1T,T3T,ALFA,
  75 :      .            EPPP,EPPM,EPPPF,EPPMF,XKTAN,TAUF)
  76 : C
  77 :        SIGF(1) = PERI * TAUF
  78 : C
  79 :        VARF(1) = EPPPF
  80 :        VARF(2) = EPPMF
  81 :        VARF(3) = EPSO      
  82 : C
  83 :        VARF(4) = PERI * XKTAN
  84 :        VARF(5) = XKN
  85 : C
  86 :        RETURN
  87 :        END
  88 :  
  89 :  
  90 :  
  91 :  
  92 :