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    CB215821  21/02/11    21:15:16     10886          
  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. *     IPVANO  (E)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  18. *                  VALEURS NODALES DE NOTRE COEFFICIENT
  19. *     IPGEOM  (E)  POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE
  20. *                  DU DOMAINE D'INTEGRATION
  21. *     IPVAEQ  (S)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS
  22. *                  NODALES EQUIVALENTES (APRES INTEGRATION NUMERIQUE)
  23. *
  24. * VARIABLES:
  25. * ----------
  26. *     XE(3,NBPTEL)  =  CORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  27. *
  28. C=======================================================================
  29.  
  30.       SUBROUTINE INTCN0 (IPVANO,IPGEOM,IPINTE,IPVAEQ)
  31.  
  32.       IMPLICIT INTEGER(I-N)
  33.       IMPLICIT REAL*8 (A-H,O-Z)
  34.  
  35.  
  36. -INC PPARAM
  37. -INC CCOPTIO
  38. -INC CCREEL
  39. -INC SMCHAML
  40. -INC SMELEME
  41. -INC SMINTE
  42. -INC SMCOORD
  43.  
  44.       PARAMETER (X2Pi=6.283185307179586476925286766559D0)
  45.       PARAMETER (X4Pi=12.566370614359172953850573533118D0)
  46.  
  47.       DIMENSION XE(3,1)
  48.  
  49. * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  50.       MELEME=IPGEOM
  51.       NBPTEL=NUM(/1)
  52.       NEL   =NUM(/2)
  53.  
  54. *  ON RECUPERE LES VALEURS NODALES DU COEFFICIENT
  55.       MELVA1=IPVANO
  56.       NBPTE1=MELVA1.VELCHE(/1)
  57.       NEL1  =MELVA1.VELCHE(/2)
  58.  
  59. *  INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  60.       N1PTEL=NBPTEL
  61.       N1EL  =NEL
  62.       N2PTEL=0
  63.       N2EL  =0
  64.       SEGINI,MELVAL
  65.  
  66. *     BOUCLE SUR LES ELEMENTS
  67.       DO iElt=1,NEL
  68.  
  69. *  COORDONNEES DES NOEUDS DE L'ELEMENT DANS LE REPERE GLOBAL
  70.         CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XE)
  71.         XValHT = MELVA1.VELCHE(1,MIN(NEL1,iElt))
  72.         IF (IFOUR.GE.12.AND.IFOUR.LE.14) THEN
  73. *         CAS DES ELEMENTS 1D AXISYMETRIQUES Q=2.Pi.R.H.T_ext
  74.           RR=XE(1,1)
  75.           XValQ=X2Pi*RR*XValHT
  76.         ELSE IF (IFOUR.EQ.15) THEN
  77. *         CAS DES ELEMENTS 1D SHPERIQUES  Q=4.Pi.R**2.H.T_ext
  78.           RR=XE(1,1)
  79.           XValQ=X4Pi*RR*RR*XValHT
  80.         ELSE
  81. *         CAS DES POI1 en 2D et 3D ==> Q=H.T_ext
  82.           XValQ=XValHT
  83.         ENDIF
  84.         VELCHE(1,iElt)=XValQ
  85.       ENDDO
  86.  
  87.       IPVAEQ=MELVAL
  88.  
  89.       END
  90.  
  91.  

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