Télécharger la plaquette Cast3M


fr en

Télécharger intcn2.eso

Retour à la liste

Numérotation des lignes : 

intcn2

  1. C INTCN2    SOURCE    OF166741  24/08/06    21:15:03     11982          
  2.  
  3. ************************************************************************
  4. *
  5. *                             I N T C N 2
  6. *                             -----------
  7. *
  8. * FONCTION:
  9. * ---------
  10. *     INTEGRATION NUMERIQUE DANS UN DOMAINE BIDIMENSIONNEL DU PRODUIT:
  11. *                     COEF. * TRANSPOSEE( N )
  12. *     COEF. : GRANDEUR PHYSIQUE REPRESENTEE PAR UN CHAMELEM
  13. *     N     : FONCTIONS DE FORME DE L'ELEMENT MAILLANT LE DOMAINE
  14. *             CONSIDERE
  15. *
  16. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  17. * -----------
  18. *
  19. *     IPVATE  (E)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  20. *                  VALEURS DE LA TEMPERATURE EXTERIEURE
  21. *     IPVAHC  (E)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES
  22. *                  VALEURS DU COEFFICIENT D'ECHANGE
  23. *     IPGEOM  (E)  POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE
  24. *                  DU DOMAINE D'INTEGRATION
  25. *     IPINTE  (E)  POINTEUR SUR UN SEGMENT MINTE CONTENANT LES
  26. *                  CARACTERISTIQUES D'INTEGRATION
  27. *     IPVAEQ  (S)  POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS
  28. *                  NODALES EQUIVALENTES
  29. * VARIABLES:
  30. * ----------
  31. *
  32. *     XE(3,NBPTEL)  =  COORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL
  33. *     A ET S        =  TABLEAUX DE TRAVAIL
  34. *
  35. * REMARQUES:
  36. * ----------
  37. *     L'UTILISATION DE CE S-P PRESUPPOSE UN PRE ET POST-TRAITEMENT
  38. *     DES SEGMENTS MELVAL PASSES EN TANT QUE PARAMETRES
  39. *
  40. * AUTEUR, DATE DE CREATION:
  41. * -------------------------
  42. *     DENIS ROBERT,LE 15 AVRIL 1988.
  43. *
  44. * LANGAGE:
  45. * --------
  46. *     ESOPE + FORTRAN77
  47. *
  48. ************************************************************************
  49.  
  50.       SUBROUTINE INTCN2 (IPVATE,IPVAHC,IPGEOM,IPINTE,IPVAEQ)
  51.  
  52.       IMPLICIT INTEGER(I-N)
  53.       IMPLICIT REAL*8 (A-H,O-Z)
  54.  
  55. -INC PPARAM
  56. -INC CCOPTIO
  57. -INC CCREEL
  58.       PARAMETER (X2Pi=6.283185307179586476925286766559D0)
  59.  
  60. -INC SMCOORD
  61. -INC SMCHAML
  62. -INC SMELEME
  63. -INC SMINTE
  64.  
  65.       SEGMENT,MMAT1
  66.         REAL*8 S(2,3),XEL(3,NBPTEL),AEL(NBPTEL)
  67.       ENDSEGMENT
  68.  
  69. C ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE
  70.       MELEME = IPGEOM
  71.       NBPTEL = meleme.NUM(/1)
  72.       NEL    = meleme.NUM(/2)
  73.  
  74. C  ON RECUPERE LES CARACTERISTIQUES D'INTEGRATION DES FACES
  75.       MINTE  = IPINTE
  76.       NBPGAU = minte.POIGAU(/1)
  77.  
  78. C  ON RECUPERE LES VALEURS DE LA TEMPERATURE
  79.       MELVA1 = IPVATE
  80.       N1_1   = MELVA1.VELCHE(/1)
  81.       N2_1   = MELVA1.VELCHE(/2)
  82.  
  83. C  ON RECUPERE LES VALEURS DU COEFFICIENT
  84.       MELVA2 = IPVAHC
  85.       N1_2   = MELVA2.VELCHE(/1)
  86.       N2_2   = MELVA2.VELCHE(/2)
  87.  
  88. C  INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES
  89.       N1PTEL = NBPTEL
  90.       N1EL   = NEL
  91.       N2PTEL = 0
  92.       N2EL   = 0
  93.       SEGINI,MELVAL
  94.  
  95.       SEGINI,MMAT1
  96.  
  97. * ========
  98. *  BOUCLE (1) SUR LES ELEMENTS
  99. * ========
  100.       DO IEL = 1, NEL
  101.  
  102.         IEMN1 = MIN(IEL,N2_1)
  103.         IEMN2 = MIN(IEL,N2_2)
  104.  
  105.         DO INOE = 1, NBPTEL
  106.           AEL(INOE) = XZERO
  107.         ENDDO
  108.  
  109. *        ON CHERCHE LES COORDONNEES DES NOEUDS DANS LE REPERE GLOBAL
  110. *
  111.         CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,IEL,XEL)
  112.  
  113. * ========
  114. *  BOUCLE (10) SUR LES POINTS D'INTEGRATION
  115. * ========
  116.         DO IGAU = 1, NBPGAU
  117.  
  118. *           CALCUL DU JACOBIEN AU POINT DE GAUSS CONSIDERE
  119.           DO IP = 1, 2
  120.             IP1 = IP + 1
  121.             DO IQ = 1, 3
  122.               r_z = XZERO
  123.               DO INOE = 1, NBPTEL
  124.                 r_z = r_z + SHPTOT(IP1,INOE,IGAU)*XEL(IQ,INOE)
  125.               END DO
  126.               S(IP,IQ) = r_z
  127.             END DO
  128.           END DO
  129.           S1 = S(1,2)*S(2,3)-S(1,3)*S(2,2)
  130.           S2 = S(1,3)*S(2,1)-S(1,1)*S(2,3)
  131.           S3 = S(1,1)*S(2,2)-S(1,2)*S(2,1)
  132.           DJAC = POIGAU(IGAU) * SQRT ( S1*S1 + S2*S2 + S3*S3 )
  133.  
  134. *      CAS DES ELEMENTS AXISYMETRIQUES
  135.           IF (IFOMOD.EQ.0) THEN
  136.             RR = XZERO
  137.             DO INOE = 1, NBPTEL
  138.               RR = RR + SHPTOT(1,INOE,IGAU)*XEL(1,INOEL)
  139.             END DO
  140.             DJAC = X2Pi * RR * DJAC
  141.           ENDIF
  142.  
  143.           INMN1 = MIN(IGAU,N1_1)
  144.           INMN2 = MIN(IGAU,N1_2)
  145.           r_z = MELVA1.VELCHE(INMN1,IEMN1) * MELVA2.VELCHE(INMN2,IEMN2)
  146.  
  147.           r_z = r_z * DJAC
  148.           DO INOE = 1, NBPTEL
  149.             AEL(INOE) = AEL(INOE) + SHPTOT(1,INOE,IGAU) * r_z
  150.           END DO
  151.  
  152.         END DO
  153. * ========
  154. *  BOUCLE (10) FIN
  155. * ========
  156.  
  157.         DO INOE = 1, NBPTEL
  158.           VELCHE(INOE,IEL) = AEL(INOE)
  159.         END DO
  160.  
  161.       END DO
  162. * ========
  163. *  BOUCLE (1) - FIN
  164. * ========
  165.  
  166.       SEGSUP,MMAT1
  167.  
  168.       IPVAEQ = MELVAL
  169.  
  170. c      return
  171.       END
  172.  
  173.  
  174.  

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