Télécharger capa5.eso

Retour à la liste

Numérotation des lignes :

capa5
  1. C CAPA5 SOURCE CB215821 17/01/16 21:15:07 9279
  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= 3 en THERMIQUE RHO CP EPAI =
  21. C= 2 en DIFFUSION CDIF EPAI =
  22. C= IVAPHA (E) Segment MPTVAL pour changement de phase (ACTIF) =
  23. C= NPHAT (E) Nombre de composantes recherchees =
  24. C= VACOMP (S) Valeur de Rho.Cp au point de Gauss... =
  25. C= =
  26. C= Denis ROBERT, le 18 fevrier 1988. =
  27. C=======================================================================
  28.  
  29. SUBROUTINE CAPA5 (iGau,iElt,IVAMAT,NMATT,IVAPHA,NPHAT, VACOMP)
  30.  
  31. IMPLICIT INTEGER(I-N)
  32. IMPLICIT REAL*8 (A-H,O-Z)
  33.  
  34. -INC CCREEL
  35. -INC SMCHAML
  36.  
  37. SEGMENT MPTVAL
  38. INTEGER IPOS(NS),NSOF(NS),IVAL(NCOSOU)
  39. CHARACTER*16 TYVAL(NCOSOU)
  40. ENDSEGMENT
  41.  
  42. DIMENSION VACOMP(NMATT)
  43.  
  44. C ===
  45. C 1 - Recuperation des caracteristiques materielles
  46. C Ces composantes sont obligatoires donc IVAL(i) n'est pas nul !
  47. C ===
  48. MPTVAL = IVAMAT
  49. DO i = 1, NMATT
  50. VACOMP(i)=XZero
  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.  
  60. C CAS THERMIQUE on fait RHO.CP
  61. VACOMP(1) = VALRHO * VACOMP(2)
  62. RETURN
  63. END
  64.  
  65.  
  66. C ===
  67. C 2 - Ajout a la capacite calorifique du terme de chaleur latente lie
  68. C a un changement de phase
  69. C CB215821 - 2017 - Fait dans PROPHA et plus ici !!!
  70. C ===
  71. C IF (IVAPHA.NE.0) THEN
  72. C MPTVAL = IVAPHA
  73. Cc* NPHAT = 1 (normalement !)
  74. C MELVAL = IVAL(1)
  75. C IGMN = MIN(iGau,VELCHE(/1))
  76. C IEMN = MIN(iElt,VELCHE(/2))
  77. C VACOMP(1) = VACOMP(1) + VALRHO * VELCHE(iGau,iElt)
  78. C ENDIF
  79. C
  80. C RETURN
  81. C END
  82.  
  83.  
  84.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales