capa5
C CAPA5 SOURCE CB215821 17/01/16 21:15:07 9279 C======================================================================= C= C A P A 5 = C= --------- = C= = C= Fonction : = C= ---------- = C= Calcul, en un point de Gauss, du terme Rho.Cp intervenant dans la = C= contribution de ce point de Gauss a la matrice de CAPACITE d'un = C= element fini. Ce terme prend egalement en compte la presence d'un = C= eventuel changement de phase. = C= = C= Parametres : (E)=Entree (S)=Sortie = C= ------------ = C= iGau (E) Numero du point de Gauss de l'element considere = C= iElt (E) Numero de l'ELEMENT fini considere = C= IVAMAT (E) Segment MPTVAL (ACTIF) = C= NMATT (E) Nombre de composantes recherchees = C= 3 en THERMIQUE RHO CP EPAI = C= 2 en DIFFUSION CDIF EPAI = C= IVAPHA (E) Segment MPTVAL pour changement de phase (ACTIF) = C= NPHAT (E) Nombre de composantes recherchees = C= VACOMP (S) Valeur de Rho.Cp au point de Gauss... = C= = C= Denis ROBERT, le 18 fevrier 1988. = C======================================================================= SUBROUTINE CAPA5 (iGau,iElt,IVAMAT,NMATT,IVAPHA,NPHAT, VACOMP) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) -INC CCREEL -INC SMCHAML SEGMENT MPTVAL INTEGER IPOS(NS),NSOF(NS),IVAL(NCOSOU) CHARACTER*16 TYVAL(NCOSOU) ENDSEGMENT DIMENSION VACOMP(NMATT) C === C 1 - Recuperation des caracteristiques materielles C Ces composantes sont obligatoires donc IVAL(i) n'est pas nul ! C === MPTVAL = IVAMAT DO i = 1, NMATT VACOMP(i)=XZero MELVAL = IVAL(i) c* IF (IVAL(i).NE.0) THEN IGMN = MIN(iGau,VELCHE(/1)) IEMN = MIN(iElt,VELCHE(/2)) VACOMP(i) = VELCHE(IGMN,IEMN) c* ENDIF ENDDO VALRHO = VACOMP(1) C CAS THERMIQUE on fait RHO.CP VACOMP(1) = VALRHO * VACOMP(2) RETURN END C === C 2 - Ajout a la capacite calorifique du terme de chaleur latente lie C a un changement de phase C CB215821 - 2017 - Fait dans PROPHA et plus ici !!! C === C IF (IVAPHA.NE.0) THEN C MPTVAL = IVAPHA Cc* NPHAT = 1 (normalement !) C MELVAL = IVAL(1) C IGMN = MIN(iGau,VELCHE(/1)) C IEMN = MIN(iElt,VELCHE(/2)) C VACOMP(1) = VACOMP(1) + VALRHO * VELCHE(iGau,iElt) C ENDIF C C RETURN C END
© Cast3M 2003 - Tous droits réservés.
Mentions légales