Télécharger chaleu.eso

Retour à la liste

Numérotation des lignes :

  1. C CHALEU SOURCE FANDEUR 16/11/30 21:15:06 9222
  2.  
  3. C=======================================================================
  4. C= C H A L E U =
  5. C= ----------- =
  6. C= =
  7. C= OPERATEUR CAST3M "SOURCE" : =
  8. C= --------------------------- =
  9. C= FF1 = 'SOURCE' MODL1 | S1 MAIL1 | ( CARA1 ) ( 'ELEM' ) ; =
  10. C= | CHP1 | =
  11. C= | CHEL1 | =
  12. C= =
  13. C= Cet operateur sert a calculer les flux nodaux equivalents a une =
  14. C= source volumique de chaleur (CHPOINT au second membre) =
  15. C= =
  16. C= ARGUMENTS : =
  17. C= ----------- =
  18. C= MODL1 (MMODEL) Modele (global) associe a la structure =
  19. C= S1 (FLOTTANT) Valeur algebrique de la source (constante) =
  20. C= MAI1 (MAILLAGE) Partie de la structure ou on impose la source =
  21. C= de chaleur de valeur S1. =
  22. C= CHP1 (CHPOINT) Valeurs algebriques des sources aux NOEUDS =
  23. C= CHEL1 (MCHAML) Valeurs algebriques des sources par ELEMENT =
  24. C= CARA1 (MCHAML) Caracteristiques geometriques (utilisees dans =
  25. C= le cas des elements de type COQUE et BARRe) =
  26. C= Sous-type 'CARACTERISTIQUES' =
  27. C= =
  28. C= RESULTAT : =
  29. C= ---------- =
  30. C= FF1 (CHPOINT) Flux nodaux equivalents =
  31. C= =
  32. C= Remarques : =
  33. C= ----------- =
  34. C= Le CHPOINT resultat FF1 est de nature DISCRETE. =
  35. C= Si le MOT 'ELEM' (facultatif) est present, le champ resultat FF1 =
  36. C= est alors un MCHAML aux NOEUDS. =
  37. C= =
  38. C= Denis ROBERT, le 26 fevrier 1988. =
  39. C=======================================================================
  40.  
  41. SUBROUTINE CHALEU
  42.  
  43. IMPLICIT INTEGER(I-N)
  44. IMPLICIT REAL*8 (A-H,O-Z)
  45.  
  46. -INC CCOPTIO
  47.  
  48. -INC SMMODEL
  49. -INC SMCHAML
  50.  
  51. CHARACTER*4 MOTELE(1), MOCOMP
  52. DATA MOTELE / 'ELEM' /
  53.  
  54. C 1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  55. C ==========================================
  56. C 1.1 - Lecture FACULTATIVE du mot 'ELEM'
  57. C =====
  58. ICLE = 0
  59. CALL LIRMOT(MOTELE,1,ICLE,0)
  60. C =====
  61. C 1.2 - Lecture OBLIGATOIRE du modele (IPMODL)
  62. C =====
  63. MOTERR(1:8)='MODELE'
  64. CALL MESLIR(-137)
  65. CALL LIROBJ('MMODEL',IPMODL,1,IRet)
  66. IF (IERR.NE.0) RETURN
  67. C =====
  68. C 1.3 - Lecture OBLIGATOIRE des sources volumiques donnees par :
  69. C 1) un CHPOINT (IPCHPO)
  70. C ou 2) d'un maillage (IPGEOM) et d'un flottant (S1)
  71. C ou 3) d'un MCHAML (IPCHEL)
  72. C =====
  73. IPCHPO = 0
  74. S1 = 0.D0
  75. IPGEOM = 0
  76. IPCHEL = 0
  77. CALL MESLIR(-138)
  78. CALL LIROBJ('CHPOINT',IPCHPO,0,iret)
  79. IF (iret.EQ.0) THEN
  80. CALL MESLIR(-139)
  81. CALL LIRREE(S1,0,iret)
  82. IF (iret.NE.0) THEN
  83. MOTERR(1:8)='MAILLAGE'
  84. CALL MESLIR(-137)
  85. CALL LIROBJ('MAILLAGE',IPGEOM,1,iret)
  86. IF (IERR.NE.0) RETURN
  87. ELSE
  88. CALL LIROBJ('MCHAML',IPCHEL,1,iret)
  89. IF (IERR.NE.0) RETURN
  90. ENDIF
  91. ENDIF
  92. C =====
  93. C 1.4 - Lecture FACULTATIVE du champ de caracteristiques (IPCARA)
  94. C =====
  95. IPCARA = 0
  96. CALL MESLIR(-145)
  97. CALL LIROBJ('MCHAML',IPCARA,0,iret)
  98. IF (IERR.NE.0) RETURN
  99. C =====
  100. C 1.5 - Lecture de la composante de la source (MOCOMP)
  101. C =====
  102. MOCOMP = ' '
  103. CALL LIRCHA(MOCOMP,0,iret)
  104. IF (IERR.NE.0) RETURN
  105.  
  106. C 2 - VERIFICATION DU MODELE A TRAITER
  107. C ======================================
  108. IPMOD1 = 0
  109. C IPMOD1 : Modele a traiter
  110. C IFORMU : indice correspondant a la formulation
  111. C ITYPEF : indicateur du type d'elements (massif, coque ou barre)
  112. CALL CHALVM(IPMODL,MOCOMP, IPMOD1,IFORMU,ITYPEF)
  113. IF (IERR.NE.0) GOTO 9900
  114.  
  115. C 3 - ANALYSE DU TERME SOURCE
  116. C =============================
  117. C Reduction du MCHAML de sources si fourni
  118. IF (IPCHEL.NE.0) THEN
  119. IPCHMZ = IPCHEL
  120. CALL REDUAF(IPCHMZ,IPMODL,IPCHEL,0,iok,kerre)
  121. IF (iok.NE.1) CALL ERREUR(kerre)
  122. IF (IERR.NE.0) GOTO 9900
  123. ENDIF
  124. C Mise sous forme de MCHAML des donnees d'entrees et
  125. C Determination du support du MCHAML de sources
  126. CALL CHALVS(IPMOD1,IFORMU,ITYPEF, IPCHPO,S1,IPGEOM,IPCHEL,
  127. & IPCHSO,ISUPSO)
  128. IF (IERR.NE.0) GOTO 9900
  129.  
  130. C 4 - ANALYSE DU CHAMP DE CARACTERISTIQUES
  131. C ==========================================
  132. C Verification de l'existence des caracteristiques
  133. C dans le cas d'elements COQUEs et BARREs.
  134. IF (IPCARA.EQ.0) THEN
  135. IF (ITYPEF.EQ.2) THEN
  136. CALL ERREUR(514)
  137. ELSE IF (ITYPEF.EQ.3) THEN
  138. CALL ERREUR(518)
  139. ENDIF
  140. IF (IERR.NE.0) GOTO 9900
  141. ELSE
  142. C En elements MASSIFs, aucune caracteristique n'est necessaire.
  143. IF (ITYPEF.EQ.1) IPCARA = 0
  144. ENDIF
  145. IF (IPCARA.NE.0) THEN
  146. IPCHMZ = IPCARA
  147. mchelm = IPCHMZ
  148. SEGACT,mchelm
  149. IF (mchelm.titche.NE.'CARACTERISTIQUES') THEN
  150. MOTERR(1:16) = 'CARACTERISTIQUES'
  151. CALL ERREUR(291)
  152. SEGDES,mchelm
  153. GOTO 9900
  154. ENDIF
  155. CALL REDUAF(IPCHMZ,IPMODL,IPCARA,0,iok,kerre)
  156. IF (iok.NE.1) CALL ERREUR(kerre)
  157. IF (IERR.NE.0) GOTO 9900
  158. ENDIF
  159.  
  160. C 5 - CALCUL DES FLUX NODAUX EQUIVALENTS
  161. C ========================================
  162. CALL CHAL1(IPMOD1,IFORMU,ITYPEF, IPCHSO,ISUPSO, IPCARA, MOCOMP,
  163. & IPCHAL)
  164. IF (IERR.NE.0 .OR. IPCHAL.EQ.0) GOTO 9900
  165.  
  166. C 6 - ECRITURE DE L'OBJET RESULTAT (CHPOINT OU MCHAML)
  167. C ======================================================
  168. IF (ICLE.EQ.0) THEN
  169. CALL CHAMPO(IPCHAL,0,IPCHAP,iret)
  170. C= Le resultat est un CHPOINT de nature DISCRETE.
  171. CALL DTCHAM(IPCHAL)
  172. IF (iret.NE.1 .OR. IERR.NE.0) RETURN
  173. CALL ECROBJ('CHPOINT',IPCHAP)
  174. ELSE
  175. CALL ECROBJ('MCHAML',IPCHAL)
  176. ENDIF
  177.  
  178. C 7 - FIN
  179. C =========
  180. 9900 CONTINUE
  181. IF (IPMOD1.NE.0) THEN
  182. MMODE1 = IPMOD1
  183. SEGSUP,MMODE1
  184. ENDIF
  185.  
  186. RETURN
  187. END
  188.  
  189.  
  190.  

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