Télécharger chaleu.eso

Retour à la liste

Numérotation des lignes :

chaleu
  1. C CHALEU SOURCE PV 20/03/31 14:33:12 10567
  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.  
  47. -INC PPARAM
  48. -INC CCOPTIO
  49. -INC SMCOORD
  50.  
  51. -INC SMMODEL
  52. -INC SMCHAML
  53.  
  54. CHARACTER*4 MOTELE(1), MOCOMP
  55. DATA MOTELE / 'ELEM' /
  56.  
  57. segact mcoord
  58. C 1 - LECTURE DES ARGUMENTS DE L'OPERATEUR
  59. C ==========================================
  60. C 1.1 - Lecture FACULTATIVE du mot 'ELEM'
  61. C =====
  62. ICLE = 0
  63. CALL LIRMOT(MOTELE,1,ICLE,0)
  64. C =====
  65. C 1.2 - Lecture OBLIGATOIRE du modele (IPMODL)
  66. C =====
  67. MOTERR(1:8)='MODELE'
  68. CALL MESLIR(-137)
  69. CALL LIROBJ('MMODEL ',IPMODL,1,IRet)
  70. CALL ACTOBJ('MMODEL ',IPMODL,1)
  71. IF (IERR.NE.0) RETURN
  72. C =====
  73. C 1.3 - Lecture OBLIGATOIRE des sources volumiques donnees par :
  74. C 1) un CHPOINT (IPCHPO)
  75. C ou 2) d'un maillage (IPGEOM) et d'un flottant (S1)
  76. C ou 3) d'un MCHAML (IPCHEL)
  77. C =====
  78. IPCHPO = 0
  79. S1 = 0.D0
  80. IPGEOM = 0
  81. IPCHEL = 0
  82. CALL MESLIR(-138)
  83. CALL LIROBJ('CHPOINT ',IPCHPO,0,iret)
  84. IF (iret.EQ.0) THEN
  85. CALL MESLIR(-139)
  86. CALL LIRREE(S1,0,iret)
  87. IF (iret.NE.0) THEN
  88. MOTERR(1:8)='MAILLAGE'
  89. CALL MESLIR(-137)
  90. CALL LIROBJ('MAILLAGE',IPGEOM,1,iret)
  91. CALL ACTOBJ('MAILLAGE',IPGEOM,1)
  92. ELSE
  93. CALL LIROBJ('MCHAML ',IPCHEL,1,iret)
  94. CALL ACTOBJ('MCHAML ',IPCHEL,1)
  95. ENDIF
  96.  
  97. ELSE
  98. CALL ACTOBJ('CHPOINT ',IPCHPO,1)
  99. ENDIF
  100. IF (IERR.NE.0) RETURN
  101. C =====
  102. C 1.4 - Lecture FACULTATIVE du champ de caracteristiques (IPCARA)
  103. C =====
  104. IPCARA = 0
  105. IPCAR1 = 0
  106. CALL MESLIR(-145)
  107. CALL LIROBJ('MCHAML ',IPCARA,0,iret)
  108. IF (iret .EQ. 1) CALL ACTOBJ('MCHAML ',IPCARA,1)
  109. IF (IERR.NE.0) RETURN
  110. C =====
  111. C 1.5 - Lecture de la composante de la source (MOCOMP)
  112. C =====
  113. MOCOMP = ' '
  114. CALL LIRCHA(MOCOMP,0,iret)
  115. IF (IERR.NE.0) RETURN
  116.  
  117. C 2 - VERIFICATION DU MODELE A TRAITER
  118. C ======================================
  119. IPMOD1 = 0
  120. C IPMOD1 : Modele a traiter
  121. C IFORMU : indice correspondant a la formulation
  122. C ITYPEF : indicateur du type d'elements (massif, coque ou barre)
  123. CALL CHALVM(IPMODL,MOCOMP, IPMOD1,IFORMU,ITYPEF)
  124. IF (IERR.NE.0) GOTO 9900
  125.  
  126. C 3 - ANALYSE DU TERME SOURCE
  127. C =============================
  128. C Reduction du MCHAML de sources si fourni
  129. IF (IPCHEL.NE.0) THEN
  130. IPCHMZ = IPCHEL
  131. CALL REDUAF(IPCHMZ,IPMOD1,IPCHEL,0,iok,kerre)
  132. IF (iok.NE.1) CALL ERREUR(kerre)
  133. IF (IERR.NE.0) GOTO 9900
  134. C Si IPCHEL, verif. si caracteristiques sont dans MCHAML source
  135. IF (IPCARA.EQ.0.AND.(ITYPEF.EQ.2.OR.ITYPEF.EQ.3)) THEN
  136. CALL ECRCHA('GEOM')
  137. CALL ECROBJ('MMODEL ',IPMOD1)
  138. CALL EXTRAI
  139. IF (IERR.NE.0) GOTO 9900
  140. CALL LIROBJ('LISTMOTS',MLMOT1,1,IRET)
  141. IF (IERR.NE.0) GOTO 9900
  142. CALL EXCOC2(IPCHEL,MLMOT1,IPCAR1,MLMOT1,1)
  143. IF (IERR.NE.0) GOTO 9900
  144. IPCARA = IPCAR1
  145. ENDIF
  146. ENDIF
  147. C Mise sous forme de MCHAML des donnees d'entrees et
  148. C Determination du support du MCHAML de sources
  149. CALL CHALVS(IPMOD1,IFORMU,ITYPEF, IPCHPO,S1,IPGEOM,IPCHEL,
  150. & IPCHSO,ISUPSO)
  151. IF (IERR.NE.0) GOTO 9900
  152.  
  153. C 4 - ANALYSE DU CHAMP DE CARACTERISTIQUES
  154. C ==========================================
  155. C Verification de l'existence des caracteristiques
  156. C dans le cas d'elements COQUEs et BARREs.
  157. C Si IPCHEL, verif. si caracteristiques sont dans MCHAML source
  158. IF (IPCARA.EQ.0) THEN
  159. IF (ITYPEF.EQ.2) THEN
  160. CALL ERREUR(514)
  161. ELSE IF (ITYPEF.EQ.3) THEN
  162. CALL ERREUR(518)
  163. ENDIF
  164. IF (IERR.NE.0) GOTO 9900
  165. ELSE
  166. C En elements MASSIFs, aucune caracteristique n'est necessaire.
  167. IF (ITYPEF.EQ.1) IPCARA = 0
  168. ENDIF
  169. IF (IPCARA.NE.0) THEN
  170. IPCHMZ = IPCARA
  171. CALL REDUAF(IPCHMZ,IPMOD1,IPCARA,0,iok,kerre)
  172. IF (iok.NE.1) CALL ERREUR(kerre)
  173. IF (IERR.NE.0) GOTO 9900
  174. C
  175. mchelm = IPCARA
  176. SEGACT,mchelm
  177. IF (mchelm.titche.NE.'CARACTERISTIQUES') THEN
  178. MOTERR(1:16) = 'CARACTERISTIQUES'
  179. CALL ERREUR(291)
  180. C SEGDES,mchelm
  181. GOTO 9900
  182. ENDIF
  183. ENDIF
  184.  
  185. C 5 - CALCUL DES FLUX NODAUX EQUIVALENTS
  186. C ========================================
  187. CALL CHAL1(IPMOD1,IFORMU,ITYPEF, IPCHSO,ISUPSO, IPCARA, MOCOMP,
  188. & IPCHAL)
  189. IF (IERR.NE.0 .OR. IPCHAL.EQ.0) GOTO 9900
  190.  
  191. C 6 - ECRITURE DE L'OBJET RESULTAT (CHPOINT OU MCHAML)
  192. C ======================================================
  193. IF (ICLE.EQ.0) THEN
  194. CALL CHAMPO(IPCHAL,0,IPCHAP,iret)
  195. C= Le resultat est un CHPOINT de nature DISCRETE.
  196. CALL DTCHAM(IPCHAL)
  197. IF (iret.NE.1 .OR. IERR.NE.0) RETURN
  198. CALL ACTOBJ('CHPOINT ',IPCHAP,1)
  199. CALL ECROBJ('CHPOINT ',IPCHAP)
  200. ELSE
  201. CALL ACTOBJ('MCHAML ',IPCHAL,1)
  202. CALL ECROBJ('MCHAML ',IPCHAL)
  203. ENDIF
  204.  
  205. C 7 - FIN
  206. C =========
  207. 9900 CONTINUE
  208. IF (IPMOD1.NE.0) THEN
  209. MMODE1 = IPMOD1
  210. SEGSUP,MMODE1
  211. ENDIF
  212. IF (IPCAR1.NE.0) THEN
  213. MCHEL1 = IPCAR1
  214. SEGSUP,MCHEL1
  215. ENDIF
  216.  
  217. END
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  

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