Télécharger chaleu.eso

Retour à la liste

Numérotation des lignes :

  1. C CHALEU SOURCE CB215821 18/09/21 21:15:21 9930
  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. CALL ACTOBJ('MMODEL ',IPMODL,1)
  68. C =====
  69. C 1.3 - Lecture OBLIGATOIRE des sources volumiques donnees par :
  70. C 1) un CHPOINT (IPCHPO)
  71. C ou 2) d'un maillage (IPGEOM) et d'un flottant (S1)
  72. C ou 3) d'un MCHAML (IPCHEL)
  73. C =====
  74. IPCHPO = 0
  75. S1 = 0.D0
  76. IPGEOM = 0
  77. IPCHEL = 0
  78. CALL MESLIR(-138)
  79. CALL LIROBJ('CHPOINT',IPCHPO,0,iret)
  80. IF (iret.EQ.0) THEN
  81. CALL MESLIR(-139)
  82. CALL LIRREE(S1,0,iret)
  83. IF (iret.NE.0) THEN
  84. MOTERR(1:8)='MAILLAGE'
  85. CALL MESLIR(-137)
  86. CALL LIROBJ('MAILLAGE',IPGEOM,1,iret)
  87. IF (IERR.NE.0) RETURN
  88. ELSE
  89. CALL LIROBJ('MCHAML',IPCHEL,1,iret)
  90. IF (IERR.NE.0) RETURN
  91. ENDIF
  92. ENDIF
  93. C =====
  94. C 1.4 - Lecture FACULTATIVE du champ de caracteristiques (IPCARA)
  95. C =====
  96. IPCARA = 0
  97. CALL MESLIR(-145)
  98. CALL LIROBJ('MCHAML',IPCARA,0,iret)
  99. IF (IERR.NE.0) RETURN
  100. C =====
  101. C 1.5 - Lecture de la composante de la source (MOCOMP)
  102. C =====
  103. MOCOMP = ' '
  104. CALL LIRCHA(MOCOMP,0,iret)
  105. IF (IERR.NE.0) RETURN
  106.  
  107. C 2 - VERIFICATION DU MODELE A TRAITER
  108. C ======================================
  109. IPMOD1 = 0
  110. C IPMOD1 : Modele a traiter
  111. C IFORMU : indice correspondant a la formulation
  112. C ITYPEF : indicateur du type d'elements (massif, coque ou barre)
  113. CALL CHALVM(IPMODL,MOCOMP, IPMOD1,IFORMU,ITYPEF)
  114. IF (IERR.NE.0) GOTO 9900
  115.  
  116. C 3 - ANALYSE DU TERME SOURCE
  117. C =============================
  118. C Reduction du MCHAML de sources si fourni
  119. IF (IPCHEL.NE.0) THEN
  120. IPCHMZ = IPCHEL
  121. CALL ACTOBJ('MCHAML ',IPCHEL,1)
  122. CALL REDUAF(IPCHMZ,IPMODL,IPCHEL,0,iok,kerre)
  123. IF (iok.NE.1) CALL ERREUR(kerre)
  124. IF (IERR.NE.0) GOTO 9900
  125. ENDIF
  126. C Mise sous forme de MCHAML des donnees d'entrees et
  127. C Determination du support du MCHAML de sources
  128. CALL CHALVS(IPMOD1,IFORMU,ITYPEF, IPCHPO,S1,IPGEOM,IPCHEL,
  129. & IPCHSO,ISUPSO)
  130. IF (IERR.NE.0) GOTO 9900
  131.  
  132. C 4 - ANALYSE DU CHAMP DE CARACTERISTIQUES
  133. C ==========================================
  134. C Verification de l'existence des caracteristiques
  135. C dans le cas d'elements COQUEs et BARREs.
  136. IF (IPCARA.EQ.0) THEN
  137. IF (ITYPEF.EQ.2) THEN
  138. CALL ERREUR(514)
  139. ELSE IF (ITYPEF.EQ.3) THEN
  140. CALL ERREUR(518)
  141. ENDIF
  142. IF (IERR.NE.0) GOTO 9900
  143. ELSE
  144. C En elements MASSIFs, aucune caracteristique n'est necessaire.
  145. IF (ITYPEF.EQ.1) IPCARA = 0
  146. ENDIF
  147. IF (IPCARA.NE.0) THEN
  148. IPCHMZ = IPCARA
  149. mchelm = IPCHMZ
  150. SEGACT,mchelm
  151. IF (mchelm.titche.NE.'CARACTERISTIQUES') THEN
  152. MOTERR(1:16) = 'CARACTERISTIQUES'
  153. CALL ERREUR(291)
  154. SEGDES,mchelm
  155. GOTO 9900
  156. ENDIF
  157. CALL ACTOBJ('MCHAML ',IPCARA,1)
  158. CALL REDUAF(IPCHMZ,IPMODL,IPCARA,0,iok,kerre)
  159. IF (iok.NE.1) CALL ERREUR(kerre)
  160. IF (IERR.NE.0) GOTO 9900
  161. ENDIF
  162.  
  163. C 5 - CALCUL DES FLUX NODAUX EQUIVALENTS
  164. C ========================================
  165. CALL CHAL1(IPMOD1,IFORMU,ITYPEF, IPCHSO,ISUPSO, IPCARA, MOCOMP,
  166. & IPCHAL)
  167. IF (IERR.NE.0 .OR. IPCHAL.EQ.0) GOTO 9900
  168.  
  169. C 6 - ECRITURE DE L'OBJET RESULTAT (CHPOINT OU MCHAML)
  170. C ======================================================
  171. IF (ICLE.EQ.0) THEN
  172. CALL CHAMPO(IPCHAL,0,IPCHAP,iret)
  173. C= Le resultat est un CHPOINT de nature DISCRETE.
  174. CALL DTCHAM(IPCHAL)
  175. IF (iret.NE.1 .OR. IERR.NE.0) RETURN
  176. CALL ECROBJ('CHPOINT',IPCHAP)
  177. ELSE
  178. CALL ECROBJ('MCHAML',IPCHAL)
  179. ENDIF
  180.  
  181. C 7 - FIN
  182. C =========
  183. 9900 CONTINUE
  184. IF (IPMOD1.NE.0) THEN
  185. MMODE1 = IPMOD1
  186. SEGSUP,MMODE1
  187. ENDIF
  188.  
  189. RETURN
  190. END
  191.  
  192.  
  193.  
  194.  

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