Télécharger chama3.eso

Retour à la liste

Numérotation des lignes :

  1. C CHAMA3 SOURCE FANDEUR 19/10/08 21:15:11 10329
  2.  
  3. C=======================================================================
  4. C= C H A M A 3 =
  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 TRIDIMENSIONNEL 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= VOLU Volume elementaire au point de Gauss =
  25. C= AEL Vecteur de travail =
  26. C=======================================================================
  27.  
  28. SUBROUTINE CHAMA3 (IPSONO,IPGEOM,IPINTE, IPCHEQ)
  29.  
  30. IMPLICIT INTEGER(I-N)
  31. IMPLICIT REAL*8 (A-H,O-Z)
  32.  
  33. -INC CCOPTIO
  34. -INC CCREEL
  35.  
  36. -INC SMCHAML
  37. -INC SMELEME
  38. -INC SMINTE
  39. -INC SMCOORD
  40.  
  41. SEGMENT MMAT1
  42. REAL*8 XEL(3,NBPTEL), AEL(NBPTEL)
  43. ENDSEGMENT
  44.  
  45. C 1 - INITIALISATION
  46. C ====================
  47. C 1.1 - Recuperation d'informations sur le maillage elementaire
  48. C =====
  49. MELEME=IPGEOM
  50. SEGACT,MELEME
  51. NBPTEL=NUM(/1)
  52. NbElt=NUM(/2)
  53. C =====
  54. C 1.2 - Recuperation d'informations sur l'element fini du maillage
  55. C =====
  56. MINTE=IPINTE
  57. c* SEGACT,MINTE
  58. NBPGAU=POIGAU(/1)
  59. C* NBNO=SHPTOT(/2)
  60. C* On doit avoir -> NBNO = NBPTEL
  61. C =====
  62. C 1.3 - Activation du MCHAML contenant les valeurs des sources
  63. C =====
  64. MELVA1=IPSONO
  65. SEGACT,MELVA1
  66. NBPTE1=MELVA1.VELCHE(/1)
  67. NEL1=MELVA1.VELCHE(/2)
  68. * IPSONO s'appuie normalement sur IPINTE
  69. * On doit avoir : NBPGAU = NBPTE1 si NBPTE1 != 1 (champ constant par element)
  70. * NEL1 = NELT si NEL1 != 1 (champ uniforme)
  71. C =====
  72. C 1.4 - Initialisation du segment de travail MMAT1
  73. C =====
  74. SEGINI,MMAT1
  75. C =====
  76. C 1.5 - Initialisation du segment resultat MELVAL contenant les valeurs
  77. C des flux nodaux pour chaque element du maillage (IPGEOM)
  78. C =====
  79. N1PTEL=NBPTEL
  80. N1EL=NbElt
  81. N2PTEL=0
  82. N2EL=0
  83. SEGINI,MELVAL
  84. IPCHEQ=MELVAL
  85.  
  86. C 2 - BOUCLE SUR LES ELEMENTS DU MAILLAGE ELEMENTAIRE IMAMOD
  87. C ============================================================
  88. DO iElt=1,NbElt
  89.  
  90. IEMIN1 = MIN(NEL1,iElt)
  91. C =====
  92. C 2.1 - Mise a zero de la source nodale equivalente AEL
  93. C =====
  94. DO iNoe = 1, NBPTEL
  95. AEL(iNoe) = XZero
  96. ENDDO
  97. C =====
  98. C 2.2 - Recuperation des coordonnees GLOBALES des noeuds de l'element
  99. C =====
  100. CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,iElt,XEL)
  101. C =====
  102. C 2.3 - Boucle sur les points de Gauss de l'element iElt
  103. C =====
  104. DO iGau=1,NBPGAU
  105. C =======
  106. C 2.3.1 - Calcul du volume associe au point de Gauss iGau
  107. C =======
  108. V1=XZERO
  109. V2=XZERO
  110. V3=XZERO
  111. V4=XZERO
  112. V5=XZERO
  113. V6=XZERO
  114. V7=XZERO
  115. V8=XZERO
  116. V9=XZERO
  117. DO iNoe = 1, NBPTEL
  118. V1=V1+SHPTOT(2,iNoe,iGau)*XEL(2,iNoe)
  119. V2=V2+SHPTOT(3,iNoe,iGau)*XEL(3,iNoe)
  120. V3=V3+SHPTOT(3,iNoe,iGau)*XEL(2,iNoe)
  121. V4=V4+SHPTOT(2,iNoe,iGau)*XEL(3,iNoe)
  122. V5=V5+SHPTOT(3,iNoe,iGau)*XEL(1,iNoe)
  123. V6=V6+SHPTOT(2,iNoe,iGau)*XEL(1,iNoe)
  124. V7=V7+SHPTOT(4,iNoe,iGau)*XEL(1,iNoe)
  125. V8=V8+SHPTOT(4,iNoe,iGau)*XEL(2,iNoe)
  126. V9=V9+SHPTOT(4,iNoe,iGau)*XEL(3,iNoe)
  127. ENDDO
  128. VOLU=(V3*V9-V2*V8)*V6+(V2*V7-V5*V9)*V1+(V5*V8-V3*V7)*V4
  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