Télécharger la plaquette Cast3M


fr en

Télécharger intcn0.eso

Retour à la liste

Numérotation des lignes : 

intcn0

  1. C INTCN0    SOURCE    OF166741  24/08/06    21:15:03     11982          
  2.  
  3. C=======================================================================
  4. C=                            I N T C N 0                              =
  5. *                             -----------                              =
  6. *
  7. * FONCTION:
  8. * ---------
  9. *     INTEGRATION NUMERIQUE DANS UN DOMAINE UNIDIMENSIONNEL DU PRODUIT:
  10. *                     COEF. * TRANSPOSEE( N )
  11. *     COEF. : GRANDEUR PHYSIQUE REPRESENTEE PAR UN CHAMELEM
  12. *     N     : FONCTIONS DE FORME DE L'ELEMENT MAILLANT LE DOMAINE
  13. *             CONSIDERE
  14. *
  15. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)
  16. * -----------
  17. *     IPVATE  (E)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  18. *                  VALEURS DE LA TEMPERATURE EXTERIEURE
  19. *     IPVAHC  (E)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  20. *                  VALEURS DU COEFFICIENT D'ECHANGE
  21. *     IPGEOM  (E)  POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE
  22. *                  DU DOMAINE D'INTEGRATION
  23. *     IPVAEQ  (S)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS
  24. *                  NODALES EQUIVALENTES (APRES INTEGRATION NUMERIQUE)
  25. *
  26. * VARIABLES:
  27. * ----------
  28. *     XE(3,NBPTEL)  =  CORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  29. *
  30. C=======================================================================
  31.  
  32.       SUBROUTINE INTCN0 (IPVATE,IPVAHC,IPGEOM,IPINTE,IPVAEQ)
  33.  
  34.       IMPLICIT INTEGER(I-N)
  35.       IMPLICIT REAL*8 (A-H,O-Z)
  36.  
  37. -INC PPARAM
  38. -INC CCOPTIO
  39. -INC CCREEL
  40.       PARAMETER (X2Pi= 6.283185307179586476925286766559D0)
  41.       PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  42.  
  43. -INC SMCHAML
  44. -INC SMELEME
  45. -INC SMINTE
  46. -INC SMCOORD
  47.  
  48.       DIMENSION XE(3,1)
  49.  
  50. * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  51.       MELEME = IPGEOM
  52.       NBPTEL = meleme.NUM(/1)
  53.       NEL    = meleme.NUM(/2)
  54.  
  55. *  ON RECUPERE LES VALEURS DE LA TEMPERATURE
  56.       MELVA1 = IPVATE
  57.       N1_1   = MELVA1.VELCHE(/1)
  58.       N1_2   = MELVA1.VELCHE(/2)
  59.  
  60. *  ON RECUPERE LES VALEURS DU COEFFICIENT
  61.       MELVA2 = IPVAHC
  62.       N2_1   = MELVA2.VELCHE(/1)
  63.       N2_2   = MELVA2.VELCHE(/2)
  64.  
  65. *  INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  66.       N1PTEL = NBPTEL
  67.       N1EL   = NEL
  68.       N2PTEL = 0
  69.       N2EL   = 0
  70.       SEGINI,MELVAL
  71.  
  72. *   BOUCLE SUR LES ELEMENTS
  73.       DO iElt = 1, NEL
  74.  
  75. *     COORDONNEES DES NOEUDS DE L'ELEMENT DANS LE REPERE GLOBAL
  76.         CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XE)
  77.  
  78. * Calcul de H().T_ext
  79.         IEMN1 = MIN(iElt,N1_2)
  80.         IEMN2 = MIN(iElt,N2_2)
  81.         XVal = MELVA1.VELCHE(1,IEMN1) * MELVA2.VELCHE(1,IEMN2)
  82.  
  83. *     CAS DES ELEMENTS 1D AXISYMETRIQUES Q = 2.Pi.R. H().T_ext
  84.         IF (IFOUR.GE.12.AND.IFOUR.LE.14) THEN
  85.           XVal = X2Pi * XE(1,1) * XVal
  86.  
  87. *     CAS DES ELEMENTS 1D SHPERIQUES  Q = 4.Pi.R**2. H().T_ext
  88.         ELSE IF (IFOUR.EQ.15) THEN
  89.           RR   = XE(1,1)
  90.           XVal = X4Pi * (RR*RR) * XVal
  91.  
  92. *     CAS DES POI1 en 2D et 3D ==> Q = H().T_ext
  93. c*        ELSE
  94. c*          XVal = XVal
  95.         ENDIF
  96.  
  97.         melval.VELCHE(1,iElt) = XVal
  98.  
  99.       ENDDO
  100.  
  101.       IPVAEQ = MELVAL
  102.  
  103. c      return
  104.       END
  105.  
  106.  
  107.  

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