Esope source files

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