Télécharger chama2.eso

Retour à la liste

Numérotation des lignes :

  1. C CHAMA2 SOURCE FANDEUR 19/10/08 21:15:10 10329
  2.  
  3. C=======================================================================
  4. C= C H A M A 2 =
  5. C= ----------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Calcul du MCHAML (champ par element aux noeuds) des flux nodaux =
  10. C= equivalents a des sources volumiques. Cas BIDIMENSIONNEL MASSIF =
  11. C= =
  12. C= Parametres : (E)=Entree (S)=Sortie =
  13. C= ------------ =
  14. C= IPSONO (E) Pointeur sur le segment MELVAL des SOURCES =
  15. C= IPGEOM (E) Pointeur sur un objet MAILLAGE elementaire =
  16. C= IPINTE (E) Pointeur sur un segment MINTE contenant les =
  17. C= caracteristiques d'integration =
  18. C= IPCHEQ (S) Pointeur sur un segment MELVAL contenant les =
  19. C= flux de chaleur nodaux equivalents =
  20. C= =
  21. C= Variables locales : =
  22. C= ------------------- =
  23. C= XEL(3,NBPTEL) Coordonnees GLOBALES des noeuds d'un element =
  24. C= Surf Surface elementaire au point de Gauss =
  25. C= AEL Vecteur de travail =
  26. C= =
  27. C= Denis ROBERT, le 16 fevrier 1988. =
  28. C=======================================================================
  29.  
  30. SUBROUTINE CHAMA2 (IPSONO,IPGEOM,IPINTE, IPCHEQ)
  31.  
  32. IMPLICIT INTEGER(I-N)
  33. IMPLICIT REAL*8 (A-H,O-Z)
  34.  
  35. -INC CCOPTIO
  36. -INC CCREEL
  37.  
  38. -INC SMCHAML
  39. -INC SMELEME
  40. -INC SMINTE
  41. -INC SMCOORD
  42.  
  43. SEGMENT MMAT1
  44. REAL*8 XEL(3,NBPTEL), AEL(NBPTEL)
  45. ENDSEGMENT
  46.  
  47. C= Quelques constantes (2.Pi)
  48. PARAMETER (X2Pi=6.283185307179586476925286766559D0)
  49.  
  50. C 1 - INITIALISATION
  51. C ====================
  52. C 1.1 - Recuperation d'informations sur le maillage elementaire
  53. C =====
  54. MELEME=IPGEOM
  55. SEGACT,MELEME
  56. NBPTEL=NUM(/1)
  57. NbElt=NUM(/2)
  58. C =====
  59. C 1.2 - Recuperation d'informations sur l'element fini du maillage
  60. C =====
  61. MINTE=IPINTE
  62. C* SEGACT,MINTE
  63. NBPGAU=POIGAU(/1)
  64. C* NBNO=SHPTOT(/2)
  65. C* On doit avoir -> NBNO = NBPTEL
  66. C =====
  67. C 1.3 - Activation du MCHAML contenant les valeurs des sources
  68. C =====
  69. MELVA1=IPSONO
  70. SEGACT,MELVA1
  71. NBPTE1=MELVA1.VELCHE(/1)
  72. NEL1=MELVA1.VELCHE(/2)
  73. * IPSONO s'appuie normalement sur IPINTE
  74. * On doit avoir : NBPGAU = NBPTE1 si NBPTE1 != 1 (champ constant par element)
  75. * NEL1 = NELT si NEL1 != 1 (champ uniforme)
  76. C =====
  77. C 1.4 - Initialisation du segment de travail MMAT1
  78. C =====
  79. SEGINI,MMAT1
  80. C =====
  81. C 1.5 - Initialisation du segment resultat MELVAL contenant les valeurs
  82. C des flux nodaux pour chaque element du maillage (IPGEOM)
  83. C =====
  84. N1PTEL=NBPTEL
  85. N1EL=NbElt
  86. N2PTEL=0
  87. N2EL=0
  88. SEGINI,MELVAL
  89. IPCHEQ=MELVAL
  90.  
  91. C 2 - BOUCLE SUR LES ELEMENTS DU MAILLAGE ELEMENTAIRE IMAMOD
  92. C ============================================================
  93. DO iElt=1,NbElt
  94.  
  95. IEMIN1 = MIN(NEL1,iElt)
  96. C =====
  97. C 2.1 - Mise a zero de la source nodale equivalente AEL
  98. C =====
  99. DO iNoe = 1, NBPTEL
  100. AEL(iNoe) = XZero
  101. ENDDO
  102. C =====
  103. C 2.2 - Recuperation des coordonnees GLOBALES des noeuds de l'element
  104. C =====
  105. CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XEL)
  106. C =====
  107. C 2.3 - Boucle sur les points de Gauss de l'element iElt
  108. C =====
  109. DO iGau=1,NBPGAU
  110. C =======
  111. C 2.3.1 - Calcul du volume associe au point de Gauss iGau
  112. C Pour les elements MASSIFs, traitement du cas axisymetrique
  113. C =======
  114. S1=XZERO
  115. S2=XZERO
  116. S3=XZERO
  117. S4=XZERO
  118. DO iNoe = 1, NBPTEL
  119. S1=S1+SHPTOT(2,iNoe,iGau)*XEL(1,iNoe)
  120. S2=S2+SHPTOT(3,iNoe,iGau)*XEL(1,iNoe)
  121. S3=S3+SHPTOT(3,iNoe,iGau)*XEL(2,iNoe)
  122. S4=S4+SHPTOT(2,iNoe,iGau)*XEL(2,iNoe)
  123. ENDDO
  124. VOLU=S1*S3-S2*S4
  125. IF (IFOMOD.EQ.0) THEN
  126. CALL DISTRR(XEL,SHPTOT(1,1,iGau),NBPTEL,RR)
  127. VOLU=X2Pi*RR*VOLU
  128. ENDIF
  129. VOLU=ABS(VOLU)*POIGAU(iGau)
  130. C =======
  131. C 2.3.2 - Calcul de la contribution du point de Gauss a AEL
  132. C pour cet element (MCHAML aux noeuds)
  133. C =======
  134. IGMIN1 = MIN(iGau,NBPTE1)
  135. R_z = VOLU * MELVA1.VELCHE(IGMIN1,IEMIN1)
  136. DO iNoe = 1, NBPTEL
  137. AEL(iNoe) = AEL(iNoe) + SHPTOT(1,iNoe,iGau) * R_z
  138. ENDDO
  139. ENDDO
  140. C =====
  141. C 2.4 - Calcul des flux equivalents sur cet element (MCHAML aux noeuds)
  142. C =====
  143. DO iNoe = 1,NBPTEL
  144. VELCHE(iNoe,iElt) = AEL(iNoe)
  145. ENDDO
  146. ENDDO
  147.  
  148. C 3 - MENAGE : DESACTIVATION/DESTRUCTION DE SEGMENTS
  149. C ====================================================
  150. SEGSUP,MMAT1
  151.  
  152. C* RETURN
  153. END
  154.  
  155.  
  156.  

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