Esope source files

   1 : C CAPA5     SOURCE    FANDEUR   11/07/19    21:15:21     7042           
   2 : 
   3 : C=======================================================================
   4 : C=                            C A P A 5                                =
   5 : C=                            ---------                                =
   6 : C=                                                                     =
   7 : C=  Fonction :                                                         =
   8 : C=  ----------                                                         =
   9 : C=  Calcul, en un point de Gauss, du terme Rho.Cp intervenant dans la  =
  10 : C=  contribution de ce point de Gauss a la matrice de CAPACITE d'un    =
  11 : C=  element fini. Ce terme prend egalement en compte la presence d'un  =
  12 : C=  eventuel changement de phase.                                      =
  13 : C=                                                                     =
  14 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  15 : C=  ------------                                                       =
  16 : C=   iGau    (E)  Numero du point de Gauss de l'element considere      =
  17 : C=   iElt    (E)  Numero de l'ELEMENT fini considere                   =
  18 : C=   IVAMAT  (E)  Segment MPTVAL                               (ACTIF) =
  19 : C=   NMATT   (E)  Nombre de composantes recherchees                    =
  20 : C=   IVAPHA  (E)  Segment MPTVAL pour changement de phase      (ACTIF) =
  21 : C=   NPHAT   (E)  Nombre de composantes recherchees                    =
  22 : C=   VACOMP  (S)  Valeur de Rho.Cp au point de Gauss...                =
  23 : C=                                                                     =
  24 : C=  Denis ROBERT, le 18 fevrier 1988.                                  =
  25 : C=======================================================================
  26 : 
  27 :       SUBROUTINE CAPA5 (iGau,iElt,IVAMAT,NMATT,IVAPHA,NPHAT, VACOMP)
  28 : 
  29 :       IMPLICIT INTEGER(I-N)
  30 :       IMPLICIT REAL*8 (A-H,O-Z)
  31 : 
  32 : -INC CCREEL
  33 : -INC SMCHAML
  34 : 
  35 :       SEGMENT MPTVAL
  36 :         INTEGER IPOS(NS),NSOF(NS),IVAL(NCOSOU)
  37 :         CHARACTER*16 TYVAL(NCOSOU)
  38 :       ENDSEGMENT
  39 : 
  40 :       DIMENSION VACOMP(*)
  41 : 
  42 : C ===
  43 : C  1 - Recuperation des caracteristiques materielles
  44 : C      Ces composantes sont obligatoires donc IVAL(i) n'est pas nul !
  45 : C ===
  46 :       VACOMP(1) = XZero
  47 :       VACOMP(2) = XZero
  48 :       VACOMP(3) = XZero
  49 :       MPTVAL = IVAMAT
  50 :       DO i = 1, NMATT
  51 :         MELVAL = IVAL(i)
  52 : c*      IF (IVAL(i).NE.0) THEN
  53 :           IGMN = MIN(iGau,VELCHE(/1))
  54 :           IEMN = MIN(iElt,VELCHE(/2))
  55 :           VACOMP(i) = VELCHE(IGMN,IEMN)
  56 : c*      ENDIF
  57 :       ENDDO
  58 :       VALRHO = VACOMP(1)
  59 :       VACOMP(1) = VALRHO * VACOMP(2)
  60 : C ===
  61 : C  2 - Ajout a la capacite calorifique du terme de chaleur latente lie
  62 : C      a un changement de phase
  63 : C ===
  64 :       IF (IVAPHA.NE.0) THEN
  65 :         MPTVAL = IVAPHA
  66 : c* NPHAT = 1 (normalement !)
  67 :         MELVAL = IVAL(1)
  68 :         IGMN = MIN(iGau,VELCHE(/1))
  69 :         IEMN = MIN(iElt,VELCHE(/2))
  70 :         VACOMP(1) = VACOMP(1) + VALRHO * VELCHE(iGau,iElt)
  71 :       ENDIF
  72 : 
  73 :       RETURN
  74 :       END
  75 : 
  76 :