Esope source files

   1 : C RISPL2    SOURCE    CHAT      05/01/13    03:00:54     5004           
   2 : C RISPL2.eso     SOURCE     bald     MAR 28/03/95     00:00:00
   3 :         SUBROUTINE RISPL2(RI0,SIG0,DCON,DU,NSTRS,SGMT,DUp)
   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      DU                 = INCREMENT DE LA U
  13 : C      NSRTS              = NOMBRE DE COMPOSANTES DE CONTRAINTES
  14 : C      SGMT               = CONTRAINTE LIMITE NORMALE 
  15 : C
  16 : C   SORTIES  
  17 : C
  18 : C      DUp                = INCREMENT PLASTIQUE DE LA U
  19 : C
  20 : C-------------------------------------------------------------
  21 : C
  22 : c
  23 :       IMPLICIT INTEGER(I-N)
  24 :        IMPLICIT REAL*8(A-H,O-Z)
  25 :        DIMENSION RI0(NSTRS,NSTRS), DCON(NSTRS), SIG0(NSTRS)
  26 : c
  27 :        TOL1=ABS(1.D-3*SGMT)
  28 : c
  29 :        DETERM=0.D0 
  30 :        DETERM=DETERM+RI0(1,1)*RI0(2,2)
  31 :        DETERM=DETERM-RI0(1,2)*RI0(2,1)
  32 : c
  33 :        IF(ABS(SIG0(2)).LT.TOL1)THEN
  34 : c
  35 : c    le detérminant est nul
  36 : c
  37 :        DUe=0.D0
  38 : c
  39 :        ELSE
  40 : c
  41 :        DETER=0.D0
  42 :        DETER=DETER+DCON(1)*RI0(2,2)
  43 :        DETER=DETER-DCON(2)*RI0(1,2)
  44 : c
  45 :        DUe=DETER/DETERM
  46 : c
  47 :        END IF
  48 : c
  49 :        DUp=DU-DUe
  50 : c
  51 :        RETURN
  52 : c
  53 :        END
  54 : 
  55 : 
  56 :