Esope source files

   1 : C SAUFLA    SOURCE    AM        09/12/08    21:15:26     6582           
   2 :       SUBROUTINE SAUFLA(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 VALEURS DES PARAMETRES MATERIAU
  18 : C
  19 :       TDT  = CMAT(NMAT)
  20 : C
  21 : C RECUPERATION DES VARIABLES INTERNES INITIALES
  22 : C
  23 : C   FLUAGE PROPRE
  24 : C
  25 :       EPSFLU(1) = VARF(45)
  26 :       EPSFLU(2) = VARF(46)
  27 :       EPSFLU(3) = VARF(47)
  28 :       EPSFLU(4) = VARF(48)
  29 :       EPSFLU(5) = VARF(49)
  30 :       EPSFLU(6) = VARF(50)
  31 : C
  32 : C  TEST SI LE TEMPS EST NUL
  33 : C
  34 :       IF (TDT.EQ.0.D0) THEN
  35 :         GOTO 30
  36 :       ENDIF
  37 : C
  38 : C INITIALISATION DES MATRICES
  39 : C
  39 :       CALL ZERO(AN,NSTRS,1)
  39 :       CALL ZERO(BN,NSTRS,NSTRS)
  39 :       CALL ZERO(CN,NSTRS,NSTRS)
  39 :       CALL ZERO(SIGFI,NSTRS,1)
  40 : C
  41 :       DO 10 I=1,NSTRS
  42 :           SIGFI(I) = SIGF(I)
  43 :   10  CONTINUE
  44 : C
  45 : C_______________________________________________________________________
  46 : C
  47 : C  FLUAGE PROPRE AMORTISSEUR
  48 : C
  49 :       NDI=NSTRS
  49 :       CALL FLUA(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  50 :      $                 AN,BN,CN,NDI)
  51 : C
  52 :       DO 20 I=1,NSTRS
  53 :         EPSFLU(I) = EPSFLU(I) + AN(I)
  54 :         DO 20 J=1,NSTRS
  55 :           EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  56 :   20  CONTINUE
  57 : C
  58 :   30  CONTINUE
  59 : C
  60 : C  SAUVEGARDE DES DEFORMATIONS
  61 : C
  62 : C    FLUAGE PROPRE
  63 : C
  64 :       VARF(45) = EPSFLU(1) 
  65 :       VARF(46) = EPSFLU(2)
  66 :       VARF(47) = EPSFLU(3)
  67 :       VARF(48) = EPSFLU(4)
  68 :       VARF(49) = EPSFLU(5)
  69 :       VARF(50) = EPSFLU(6)
  70 : C
  71 :       END
  72 :