Esope source files

   1 : C RISPL3    SOURCE    CHAT      05/01/13    03:00:58     5004           
   2 : C RISPL3.eso     SOURCE     bald     MAR 28/03/95     00:00:00
   3 :         SUBROUTINE RISPL3(RI0,SIG0,DCON,DU1,DU2,NSTRS,SGMT,DU1p,DU2p)
   4 : C
   5 : C-------------------------------------------------------------
   6 : C
   7 : C   ENTREES
   8 : C
   9 : C      RI0(NSTRS,NSTRS)   = MATRICE DE RIGIDITE ELASTIQUE DU JOINT
  10 : C      SIG0(NSTRS)        = CONTRAINTES INITIALES (AU PAS PRECEDENT)
  11 : C      DCON(NSTRS)        = VECTEUR DES INCREMENTS DE CONTRAINTE
  12 : C      DU1                = INCREMENT DE LA U1
  13 : C      DU2                = INCREMENT DE LA U2
  14 : C      NSRTS              = NOMBRE DE COMPOSANTES DE CONTRAINTES
  15 : C      SGMT               = CONTRAINTE LIMITE NORMALE  
  16 : C
  17 : C   SORTIES  
  18 : C
  19 : C      DU1p               = INCREMENT PLASTIQUE DE LA DU1
  20 : C      DU2p               = INCREMENT PLASTIQUE DE LA DU2
  21 : C
  22 : C-------------------------------------------------------------
  23 : C
  24 : c
  25 :       IMPLICIT INTEGER(I-N)
  26 :        IMPLICIT REAL*8(A-H,O-Z)
  27 :        DIMENSION RI0(NSTRS,NSTRS), DCON(NSTRS), SIG0(NSTRS)
  28 : c
  29 :        TOL1=ABS(1.D-3*SGMT)
  30 : c
  31 :        DETERM=0.D0 
  32 :        DETERM=DETERM+RI0(1,1)*RI0(2,2)*RI0(3,3)
  33 :        DETERM=DETERM+RI0(1,2)*RI0(2,3)*RI0(3,1)
  34 :        DETERM=DETERM+RI0(1,3)*RI0(2,1)*RI0(3,2)
  35 :        DETERM=DETERM-RI0(1,3)*RI0(2,2)*RI0(3,1)
  36 :        DETERM=DETERM-RI0(1,2)*RI0(2,1)*RI0(3,3)
  37 :        DETERM=DETERM-RI0(1,1)*RI0(3,2)*RI0(2,3)
  38 : c
  39 :        IF(ABS(SIG0(3)).LT.TOL1)THEN
  40 : c
  41 : c    le detérminant est nul
  42 : c
  43 :        DU1e=0.D0
  44 :        DU2e=0.D0
  45 : c
  46 :        ELSE
  47 : c
  48 :        DETER1=0.D0
  49 :        DETER1=DETER1+DCON(1)*RI0(2,2)*RI0(3,3)
  50 :        DETER1=DETER1+RI0(1,2)*RI0(2,3)*DCON(3)
  51 :        DETER1=DETER1+RI0(1,3)*DCON(2)*RI0(3,2)
  52 :        DETER1=DETER1-RI0(1,3)*RI0(2,2)*DCON(3)
  53 :        DETER1=DETER1-RI0(1,2)*DCON(2)*RI0(3,3)
  54 :        DETER1=DETER1-DCON(1)*RI0(3,2)*RI0(2,3)
  55 : c
  56 :        DETER2=0.D0
  57 :        DETER2=DETER2+RI0(1,1)*DCON(2)*RI0(3,3)
  58 :        DETER2=DETER2+DCON(1)*RI0(2,3)*RI0(3,1)
  59 :        DETER2=DETER2+RI0(1,3)*RI0(2,1)*DCON(3)
  60 :        DETER2=DETER2-RI0(1,3)*DCON(2)*RI0(3,1)
  61 :        DETER2=DETER2-DCON(1)*RI0(2,1)*RI0(3,3)
  62 :        DETER2=DETER2-RI0(1,1)*DCON(3)*RI0(2,3)
  63 : c
  64 :        DU1e=DETER1/DETERM
  65 :        DU2e=DETER2/DETERM
  66 : c
  67 :        END IF
  68 : c
  69 :        DU1p=DU1-DU1e
  70 :        DU2p=DU2-DU2e
  71 : c
  72 :        RETURN
  73 : c
  74 :        END
  75 : 
  76 : 
  77 :