Esope source files

   1 : C CALC2     SOURCE    CHAT      05/01/12    21:45:37     5004           
   2 : C CALC2.eso     SOURCE     bald     MAR 28/03/95     00:00:00
   3 :         SUBROUTINE CALC2(KNN,KNT,KTN,KTT,SIG0,DU,DV,NSTRS,
   4 :      $                   DELTA,RI0,DCON)
   5 : C
   6 : C-------------------------------------------------------------
   7 : C
   8 : C   ENTREES
   9 : C
  10 : C     KNN,KNT,KTN,KTT    = RAIDEURS DU JOINT 2D
  11 : C     SIG0(NSTRS)        = CONTRAINTES INITIALES (AU PAS PRECEDENT)
  12 : C     DU                 = INCREMENT DE LA U
  13 : C     DV                 = INCREMENT DE LA V
  14 : C     NSTRS              = NOMBRE DE COMPOSANTES DE CONTRAINTES
  15 : C
  16 : C   SORTIES  
  17 : C
  18 : C     DELTA              = INCREMENT DE LA U EQUIVALENTE
  19 : C     RI0(NSTRS,NSTRS)   = MATRICE DE RIGIDITE DU JOINT
  20 : C     DCON(NSTRS)        = VECTEUR DES INCREMENTS DE CONTRAINTE
  21 : C
  22 : C-------------------------------------------------------------
  23 : C
  24 : c
  25 :       IMPLICIT INTEGER(I-N)
  26 :        IMPLICIT REAL*8(A-H,O-Z)
  27 :        DIMENSION SIG0(NSTRS), RI0(NSTRS,NSTRS), DCON(NSTRS)
  28 :        REAL*8 KNN, KNT, KTN, KTT
  29 : c
  30 :        TOL1=1.D-20
  31 :        TOL2=1.D-8
  32 : c
  33 : c    calcul de la matrice Dt reale (2D)
  34 : c
  35 :        IF(ABS(SIG0(1)).LE.TOL1)THEN
  36 : c
  37 : c    la tau est nulle
  38 : c       
  39 :           IF(ABS(DU).LE.TOL2)THEN
  40 : c
  41 : c       l' incrément de déformation est nul
  42 : c
  43 :              DELTA=0.D0
  44 : c
  45 :              RI0(1,1)=0.D0
  46 :              RI0(1,2)=0.D0
  47 : c
  48 :              RI0(2,1)=0.D0
  49 :              RI0(2,2)=KNN       
  50 : c
  51 :           ELSE
  52 : c
  53 :              DELTA=ABS(DU)
  54 : c
  55 :              RI0(1,1)=KTT
  56 :              RI0(1,2)=KTN*(DU/ABS(DU))
  57 : c
  58 :              RI0(2,1)=KNT*(DU/ABS(DU))
  59 :              RI0(2,2)=KNN 
  60 : c
  61 :           END IF
  62 : c
  63 :        ELSE
  64 : c
  65 :        DELTA=DU*(SIG0(1)/ABS(SIG0(1)))
  66 : c 
  67 :        RI0(1,1)=KTT
  68 :        RI0(1,2)=KTN*(SIG0(1)/ABS(SIG0(1)))
  69 : c
  70 :        RI0(2,1)=KNT*(SIG0(1)/ABS(SIG0(1)))
  71 :        RI0(2,2)=KNN
  72 : c
  73 :        END IF
  74 : c
  75 :        DCON(1)=RI0(1,1)*DU+RI0(1,2)*DV
  76 :        DCON(2)=RI0(2,1)*DU+RI0(2,2)*DV
  77 : c
  78 :        RETURN
  79 : c
  80 :        END
  81 : 
  82 :