Télécharger capaj1.eso

Retour à la liste

Numérotation des lignes :

  1. C CAPAJ1 SOURCE CB215821 21/02/11 21:15:05 10886
  2. SUBROUTINE CAPAJ1(IPMAIL,IVAMAT,NMATT,IPMATR,INFOR)
  3.  
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6.  
  7. C=======================================================================
  8. C= C A P A J 1 =
  9. C= ----------- =
  10. C= =
  11. C= Fonction : =
  12. C= ---------- =
  13. C= Calcul de la matrice de CAPACITE CALORIFIQUE d'un element JOI1 =
  14. C= =
  15. C= Clement BERTHINIER, le 03 fevrier 2021. =
  16. C=======================================================================
  17.  
  18. -INC PPARAM
  19. -INC CCOPTIO
  20. -INC SMELEME
  21. -INC SMRIGID
  22. -INC SMCHAML
  23.  
  24. SEGMENT MPTVAL
  25. INTEGER IPOS(NS),NSOF(NS),IVAL(NCOSOU)
  26. CHARACTER*16 TYVAL(NCOSOU)
  27. ENDSEGMENT
  28.  
  29. MELEME = IPMAIL
  30. NBNN = NUM(/1)
  31. NBELEM = NUM(/2)
  32.  
  33. XMATRI = IPMATR
  34. MPTVAL = IVAMAT
  35.  
  36.  
  37. C BOUCLE SUR LES ELEMENTS DU MAILLAGE ELEMENTAIRE IPMAIL
  38. C ======================================================
  39. DO iElt = 1, NBELEM
  40.  
  41. C 2 - BOUCLE SUR LES NOEUDS DES ELEMENTS SUPPORT (SEG2)
  42. C =====================================================
  43. DO iNOEU = 1, NBNN
  44. IF (INFOR .EQ. 1)THEN
  45. C Cas THERMIQUE : Calcul du terme M * Cp
  46. MELVAL = MPTVAL.IVAL(1)
  47. IGMN = MIN(iNOEU,VELCHE(/1))
  48. IEMN = MIN(iElt,VELCHE(/2))
  49. XM = VELCHE(IGMN,IEMN)
  50.  
  51. MELVAL = MPTVAL.IVAL(2)
  52. IGMN = MIN(iNOEU,VELCHE(/1))
  53. IEMN = MIN(iElt,VELCHE(/2))
  54. XC = VELCHE(IGMN,IEMN)
  55.  
  56. XRIG = XM * XC
  57.  
  58. ELSEIF(INFOR .EQ. 2)THEN
  59. C Cas DIFFUSION : Directement le CDG (CD généralise)
  60. MELVAL = MPTVAL.IVAL(1)
  61. IGMN = MIN(iNOEU,VELCHE(/1))
  62. IEMN = MIN(iElt,VELCHE(/2))
  63. XRIG = VELCHE(IGMN,IEMN)
  64. ENDIF
  65.  
  66. C Remplissage direct de la matrice
  67. RE(iNOEU,iNOEU,iElt) = XRIG
  68. ENDDO
  69. ENDDO
  70.  
  71. END
  72.  

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