Esope source files

   1 : C FLUPR1    SOURCE    AM        09/12/08    21:15:18     6582           
   2 :       SUBROUTINE FLUPR1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
   3 :      $                  AN,BN,CN,KERRE)
   4 : C
   5 : C ROUTINE APPELE DANS CAMECA
   6 : C_______________________________________________________________________
   7 : C
   8 : C SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATION 
   9 : C  DE FLUAGE PROPRE 
  10 : C
  11 : C_______________________________________________________________________
  12 : C
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8(A-H,O-Z)
  15 :       DIMENSION VARF(NVARI),CMAT(NMAT)
  16 :       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  17 :       DIMENSION A1(10),B1(10,10),C1(10,10)
  18 :       DIMENSION A2(10),B2(10,10),C2(10,10)
  19 :       DIMENSION A3(10),B3(10,10),C3(10,10)
  20 :       DIMENSION A4(10),B4(10,10),C4(10,10)
  21 :       DIMENSION A5(10),B5(10,10),C5(10,10)
  22 :       DIMENSION A6(10),B6(10,10),C6(10,10)
  23 :       DIMENSION AA(10),BA(10,10),CA(10,10)
  24 :       DIMENSION AD(10),BD(10,10),CD(10,10)
  25 : 
  26 : C  TEST SUR LA DIMENSION
  27 : 
  28 :       NDI = 10
  29 :       IF(NSTRS.GT.NDI) THEN
  30 :          PRINT *, ' FLUD - DIMENSION INSUFFISANTE'
  31 :          KERRE=99
  32 :          RETURN
  33 :       ENDIF
  34 : C
  35 : C CALCUL DES DEFORMATIONS DE FLUAGE
  36 : C   POUR LES 6 CHAINES DE KELVIN VOIGT
  37 : C
  38 : C   Chaine 1
  39 : C
  40 :       NN=1
  40 :       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  41 :      $                       A1,B1,C1,NN,NDI)
  42 : C
  43 : C   Chaine 2
  44 : C
  45 :       NN=2
  45 :       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  46 :      $                       A2,B2,C2,NN,NDI)
  47 : C
  48 : C   Chaine 3
  49 : C
  50 :       NN=3
  50 :       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  51 :      $                       A3,B3,C3,NN,NDI)
  52 : C
  53 : C   Chaine 4
  54 : C
  55 :       NN=4
  55 :       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  56 :      $                       A4,B4,C4,NN,NDI)
  57 : C
  58 : C   Chaine 5
  59 : C
  60 :       NN=5
  60 :       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  61 :      $                       A5,B5,C5,NN,NDI)
  62 : C
  63 : C   Chaine 6
  64 : C
  65 :       NN=6
  65 :       CALL FLUD(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  66 :      $                       A6,B6,C6,NN,NDI)
  67 : C
  68 : C   Amortisseur
  69 : C
  69 :       CALL FLUA(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  70 :      $                       AA,BA,CA,NDI)
  71 : C
  72 : C   fluage dessiccation
  73 : C
  73 :       CALL FLUDE1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  74 :      $                       AD,BD,CD,NDI)
  75 : C
  76 : C
  77 : C   ASSEMBLAGE DES MATRICES
  78 : C
  79 :       DO 10 I=1,NSTRS
  80 :         AN(I) = A1(I) + A2(I) + A3(I) + A4(I) + A5(I)
  81 :      $          + A6(I) + AA(I) + AD(I)
  82 :         DO 10 J=1,NSTRS
  83 :           BN(I,J) = B1(I,J) + B2(I,J) + B3(I,J)
  84 :      $              + B4(I,J) + B5(I,J) + B6(I,J)
  85 :      $              + BA(I,J) + BD(I,J)
  86 :           CN(I,J) = C1(I,J) + C2(I,J) + C3(I,J)
  87 :      $              + C4(I,J) + C5(I,J) + C6(I,J)
  88 :      $              + CA(I,J) + CD(I,J)
  89 :    10   CONTINUE 
  90 : C
  91 :       END
  92 :  
  93 :