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    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