Esope source files

   1 : C SAUFL2    SOURCE    AM        09/12/08    21:15:25     6582           
   2 :       SUBROUTINE SAUFL2(VARF,NVARI,CMAT,NMAT,SIGF,
   3 :      $           SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,SIGFI)
   4 : C
   5 : C_______________________________________________________________________
   6 : C
   7 : C SOUS PROGRAMME QUI CALCUL LA DEFORMATION 
   8 : C DE FLUAGE PROPRE A LA FIN DU PAS DE TEMPS
   9 : C_______________________________________________________________________
  10 : C
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)
  13 :       DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS)
  14 :       DIMENSION EPSFLU(6),SIGFI(NSTRS)
  15 :       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  16 : C
  17 : C RECUPERATION DES VARIABLES INTERNES INITIALES
  18 : C
  19 : C   FLUAGE PROPRE
  20 : C
  21 :       EPSFLU(1) = VARF(51)
  22 :       EPSFLU(2) = VARF(52)
  23 :       EPSFLU(3) = VARF(53)
  24 :       EPSFLU(4) = VARF(54)
  25 :       EPSFLU(5) = VARF(55)
  26 :       EPSFLU(6) = VARF(56)
  27 : C
  28 : C  TEST SI LE TEMPS EST NUL
  29 : C
  30 :       TDT = CMAT(NMAT)
  31 :       IF (TDT.EQ.0.D0) THEN
  32 :         GOTO 30
  33 :       ENDIF
  34 : C
  35 : C INITIALISATION DES MATRICES
  36 : C
  36 :       CALL ZERO(AN,NSTRS,1)
  36 :       CALL ZERO(BN,NSTRS,NSTRS)
  36 :       CALL ZERO(CN,NSTRS,NSTRS)
  36 :       CALL ZERO(SIGFI,NSTRS,1)
  37 : C
  38 :       DO 10 I=1,NSTRS
  39 :           SIGFI(I) = SIGF(I)
  40 :   10  CONTINUE
  41 : C
  42 : C_______________________________________________________________________
  43 : C
  44 : C  FLUAGE PROPRE AMORTISSEUR
  45 : C
  46 :       NDI=NSTRS
  46 :       CALL FLUDE1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  47 :      $                 AN,BN,CN,NDI)
  48 : C
  49 :       DO 20 I=1,NSTRS
  50 :         EPSFLU(I) = EPSFLU(I) + AN(I)
  51 :         DO 20 J=1,NSTRS
  52 :           EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  53 :   20  CONTINUE
  54 : C
  55 :   30  CONTINUE
  56 : C
  57 : C  SAUVEGARDE DES DEFORMATIONS
  58 : C
  59 : C    FLUAGE DESSICCATION
  60 : C
  61 :       VARF(45) = EPSFLU(1) 
  62 :       VARF(46) = EPSFLU(2)
  63 :       VARF(47) = EPSFLU(3)
  64 :       VARF(48) = EPSFLU(4)
  65 :       VARF(49) = EPSFLU(5)
  66 :       VARF(50) = EPSFLU(6)
  67 : C
  68 :       RETURN
  69 :       END
  70 :