Télécharger la plaquette Cast3M


fr en

Télécharger capaj1.eso

Retour à la liste

Numérotation des lignes : 

capaj1

  1. C CAPAJ1    SOURCE    OF166741  25/02/21    21:15:24     12166          
  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. -INC TMPTVAL
  25.  
  26.       MELEME = IPMAIL
  27.       NBNN   = NUM(/1)
  28.       NBELEM = NUM(/2)
  29.  
  30.       XMATRI = IPMATR
  31.       MPTVAL = IVAMAT
  32.  
  33. C     BOUCLE SUR LES ELEMENTS DU MAILLAGE ELEMENTAIRE IPMAIL
  34. C     ======================================================
  35.       DO iElt = 1, NBELEM
  36.  
  37. C       2 - BOUCLE SUR LES NOEUDS DES ELEMENTS SUPPORT (SEG2)
  38. C       =====================================================
  39.         DO iNOEU = 1, NBNN
  40.           IF    (INFOR .EQ. 1)THEN
  41. C           Cas THERMIQUE : Calcul du terme M * Cp
  42.             MELVAL = MPTVAL.IVAL(1)
  43.             IGMN   = MIN(iNOEU,VELCHE(/1))
  44.             IEMN   = MIN(iElt,VELCHE(/2))
  45.             XM     = VELCHE(IGMN,IEMN)
  46.  
  47.             MELVAL = MPTVAL.IVAL(2)
  48.             IGMN   = MIN(iNOEU,VELCHE(/1))
  49.             IEMN   = MIN(iElt,VELCHE(/2))
  50.             XC     = VELCHE(IGMN,IEMN)
  51.  
  52.             XRIG = XM * XC
  53.  
  54.           ELSEIF(INFOR .EQ. 2)THEN
  55. C           Cas DIFFUSION : Directement le CDG (CD généralise)
  56.             MELVAL = MPTVAL.IVAL(1)
  57.             IGMN   = MIN(iNOEU,VELCHE(/1))
  58.             IEMN   = MIN(iElt,VELCHE(/2))
  59.             XRIG   = VELCHE(IGMN,IEMN)
  60.           ENDIF
  61.  
  62. C         Remplissage direct de la matrice
  63.           RE(iNOEU,iNOEU,iElt) = XRIG
  64.         ENDDO
  65.       ENDDO
  66.  
  67. C      RETURN
  68.       END 
  69.  
  70.  

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