Télécharger formation_debutant_2_thermique.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : formation_debutant_2_thermique.dgibi
  2. *
  3. ************************************************************************
  4. * FORMATION DEBUTER AVEC CAST3M - CALCULS THERMO-MECANIQUES *
  5. * *
  6. * Modelisation du comportement thermo-mecanique d'une plaque perforee *
  7. * *
  8. * Ce fichier est la partie 2 sur 3 de la formation debutant a Cast3M *
  9. * *
  10. * Partie 2 : calculs thermiques *
  11. * restitution des objets *
  12. * conduction, convection, rayonnement *
  13. * regimes stationnaires et transitoires *
  14. * creation d'une procedure Gibiane *
  15. * instructions conditionelles, gestion des erreurs *
  16. * sauvegarde des objets *
  17. * *
  18. ************************************************************************
  19.  
  20.  
  21.  
  22.  
  23. * RESTITUTION DES OBJETS CREES
  24. * ----------------------------
  25. OPTI 'REST' 'formation_debutant_1_maillage.sauv' ;
  26.  
  27.  
  28. * OBJET LOGIQUE POUR INDIQUER SI L'ON FAIT LES TRACES OU NON
  29. * ----------------------------------------------------------
  30. ITRAC = FAUX ;
  31.  
  32.  
  33.  
  34.  
  35.  
  36.  
  37. ************************************************************************
  38. * CHAPITRE 2 : THERMIQUE LINEAIRE STATIONNAIRE *
  39. ************************************************************************
  40.  
  41.  
  42. * MODELE THERMIQUE (CONDUCTION) A MATERIAU UNIFORME ET CONSTANT
  43. * -------------------------------------------------------------
  44. MOT = MODE SU 'THERMIQUE' 'ISOTROPE' 'CONDUCTION' ;
  45. MAT = MATE MOT 'K' CONDUMAT 'C' CAPAMAT 'RHO' RHOMAT ;
  46.  
  47.  
  48. * CALCUL DE LA MATRICE DE CONDUCTIVITE
  49. * ------------------------------------
  50. CON = COND MOT MAT ;
  51.  
  52.  
  53. * MATRICE DE BLOCAGE : TEMPERATURE IMPOSEE SUR LE TROU
  54. * ----------------------------------------------------
  55. BLT = BLOQ CE 'T' ;
  56.  
  57.  
  58. * FLUX NODAUX ASSOCIES AU BLOCAGE
  59. * -------------------------------
  60. TMAX = 250. ;
  61. FLT1 = DEPI BLT TMAX ;
  62.  
  63.  
  64. * FLUX IMPOSE SUR LA LIGNE GAUCHE
  65. * -------------------------------
  66. FLT2 = FLUX MOT LIGA -3.5E4 ;
  67.  
  68.  
  69. * CALCUL DU CHAMP DE TEMPERATURES PAR APPEL AU SOLVEUR
  70. * ----------------------------------------------------
  71. TCON1 = RESO (CON ET BLT) (FLT1 ET FLT2) ;
  72.  
  73.  
  74. * AFFICHAGE DU CHAMP DE TEMPERATURE
  75. * ---------------------------------
  76. LISO1 = PROG 0. 'PAS' (TMAX / 15.) TMAX ;
  77. SI ITRAC ;
  78. TRAC TCON1 SU ;
  79. TRAC TCON1 SU CSU ;
  80. TRAC TCON1 SU CSU LISO1 ;
  81. TRAC TCON1 SU CSU LISO1 'TITR' '[2] Temperature' ;
  82. FINSI ;
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93. ************************************************************************
  94. * CHAPITRE 3 : THERMIQUE LINEAIRE TRANSITOIRE *
  95. * UTILISATION DE PASAPAS *
  96. ************************************************************************
  97.  
  98.  
  99. * INSTANT FINAL DU CALCUL THERMIQUE
  100. * ---------------------------------
  101. TPSFIN = 5.E4 ;
  102.  
  103.  
  104. * CHAMP DE TEMPERATURES INITIALES (UNIFORME)
  105. * ------------------------------------------
  106. T0 = 25. ;
  107. T_INI = MANU 'CHPO' 1 SU 'T' T0 ;
  108.  
  109.  
  110. * CHARGEMENTS THERMIQUES POUR LES TEMPERATURES ET FLUX IMPOSES
  111. * MONTEE RAPIDE PUIS MAINTENUE AU COURS DU TEMPS
  112. * ------------------------------------------------------------
  113. LIST1 = PROG 0. (0.1 * TPSFIN) TPSFIN ;
  114. LIST2 = PROG (T0 / TMAX) 1. 1. ;
  115. EVT1 = EVOL 'MANU' LIST1 LIST2 ;
  116. CHT1 = CHAR 'TIMP' FLT1 EVT1 ;
  117.  
  118. LIST3 = PROG 0. (0.2 * TPSFIN) (0.5 * TPSFIN) TPSFIN ;
  119. LIST4 = PROG 0. 0. 1. 1. ;
  120. EVT2 = EVOL 'MANU' LIST3 LIST4 ;
  121. CHT2 = CHAR 'Q' FLT2 EVT2 ;
  122. CHT = CHT1 ET CHT2 ;
  123.  
  124.  
  125. * DEFINITION DE LA TABLE D'ARGUMENTS A FOURNIR EN DONNEE D'ENTREE A LA
  126. * PROCEDURE PASAPAS
  127. * --------------------------------------------------------------------
  128. TAB1 = TABL ;
  129. TAB1 . 'MODELE' = MOT ;
  130. TAB1 . 'CARACTERISTIQUES' = MAT ;
  131. TAB1 . 'BLOCAGES_THERMIQUES' = BLT ;
  132. TAB1 . 'CHARGEMENT' = CHT ;
  133. TAB1 . 'TEMPS_CALCULES' = PROG 0. 'PAS' (0.02 * TPSFIN) TPSFIN ;
  134. TAB1 . 'TEMPERATURES' = TABL ;
  135. TAB1 . 'TEMPERATURES' . 0 = T_INI ;
  136.  
  137.  
  138. * APPEL A PASAPAS
  139. * ---------------
  140. PASAPAS TAB1 ;
  141.  
  142.  
  143. * POST TRAITEMENT
  144. * ---------------
  145. * EVOLUTION TEMPORELLE DE LE TEMPERATURE EN DEUX POINTS POUR VERIFIER
  146. * QUE L'ON A ATTEINT L'ETAT STATIONNAIRE
  147. PMIL = SU POIN 'PROC' ((0.5 * LONG) (0.5 * HAUT)) ;
  148. EV1 = EVOL 'ROUG' 'TEMP' TAB1 'TEMPERATURES' 'T' PMIL ;
  149. EV2 = EVOL 'BLEU' 'TEMP' TAB1 'TEMPERATURES' 'T' PG ;
  150. SI ITRAC ;
  151. DESS (EV1 ET EV2) 'TITR' '[3] Temperature vs Temps' ;
  152. FINSI ;
  153.  
  154. * BOUCLE SUR TOUS LES PAS DE TEMPS ET TRACE DU CHAMP DE TEMPERATURES
  155. N1 = DIME (TAB1 . 'TEMPERATURES') ;
  156. SI ITRAC ;
  157. REPE B1 N1 ;
  158. T_I = TAB1 . 'TEMPERATURES' . (&B1 - 1) ;
  159. TPS_I = TAB1 . 'TEMPS' . (&B1 - 1) ;
  160. MOT_I = CHAI '[3] Temperature au temps ' TPS_I ;
  161. TRAC T_I SU CSU 'TITR' MOT_I LISO1 ;
  162. FIN B1 ;
  163. FINSI ;
  164.  
  165. * VECTEUR FLUX DE CHALEUR ==> CREATION D'UNE PROCEDURE
  166. DEBP VECFLU CHP1*'CHPOINT' MOD1*'MMODEL' MAT1*'MCHAML' MOT1/'MOT' ;
  167. * CALCUL DU GRADIENT DE T ET CHANGEMENT DU TYPE
  168. G1 = GRAD CHP1 MOD1 ;
  169. G2 = CHAN 'TYPE' G1 'CARACTERISTIQUES' ;
  170. * MULTIPLICATION DES CHAMPS ENTRE EUX
  171. Q = MAT1 * G2 (MOTS 'K' 'K') (MOTS 'T,X' 'T,Y') (MOTS 'QX' 'QY') ;
  172. Q = -1. * Q ;
  173. * CREATION D'UN OBJET VECTEUR
  174. SI (EGA (TYPE MOT1) 'MOT') ;
  175. SI (EGA MOT1 'COUL') ;
  176. QP = CHAN 'CHPO' Q MOD1 ;
  177. VEC1 = @VECOUL QP 2.E-6 ;
  178. SINON ;
  179. ERRE '--> ERREUR dans la procedure VECFLU, le mot clef attendu est COUL <--' ;
  180. FINSI ;
  181. SINON ;
  182. VEC1 = VECT Q MOD1 (MOTS 'QX' 'QY') 2.E-6 ;
  183. FINSI ;
  184. FINP VEC1 ;
  185.  
  186. * BOUCLE SUR TOUS LES PAS DE TEMPS ET TRACE VECTEUR FLUX DE CHALEUR
  187. SI ITRAC ;
  188. REPE B1 N1 ;
  189. T_I = TAB1 . 'TEMPERATURES' . (&B1 - 1) ;
  190. VF_I = VECFLU T_I MOT MAT 'COUL' ;
  191. TPS_I = TAB1 . 'TEMPS' . (&B1 - 1) ;
  192. MOT_I = CHAI '[3] Vecteur flux de chaleur au temps ' TPS_I ;
  193. TRAC VF_I CSU 'TITR' MOT_I ;
  194. FIN B1 ;
  195. FINSI ;
  196.  
  197. * VECTEUR FLUX ET CHAMP DE TEMPRATURE SOUS FORME DE LIGNES D'ISOVALEURS
  198. OPTI 'ISOV' 'LIGN' ;
  199. SI ITRAC ;
  200. TRAC VF_I T_I SU CSU 15 'TITR' MOT_I ;
  201. FINSI ;
  202. OPTI 'ISOV' 'SURF' ;
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212. ************************************************************************
  213. * CHAPITRE 4 : THERMIQUE NON LINEAIRE TRANSITOIRE *
  214. * AJOUT DE CONVECTION ET DE RAYONNEMENT *
  215. ************************************************************************
  216.  
  217.  
  218. * RECUPERATION DE LA LIGNE A DROITE
  219. * ---------------------------------
  220. X1 = COOR SU 1 ;
  221. PLD = X1 POIN 'EGAL' LONG ;
  222. LD = CSU ELEM 'APPU' 'STRI' PLD ;
  223.  
  224.  
  225. * MODELE DE CONVECTION SUR LA LIGNE DE DROITE
  226. * ET MATERIAU REPRESENTANT LE COEFFICIENT D'ECHANGE
  227. * -------------------------------------------------
  228. MOC = MODE LD 'THERMIQUE' 'CONVECTION' ;
  229. MAC = MATE MOC 'H' 100. ;
  230.  
  231.  
  232. * CHARGEMENT DE CONVECTION REPRESENTANT LA TEMPERATURE
  233. * EXTERIEURE (MAINTENUE)
  234. * ----------------------------------------------------
  235. CHTC = MANU 'CHPO' LD 1 'T' T0 ;
  236. CHACONV = CHAR 'TECO' CHTC EVT2 ;
  237.  
  238.  
  239. * MODELE DE RAYONNEMENT SUR LA LIGNE DU HAUT
  240. * ET MATERIAU REPRESENTANT L'EMISSIVITE
  241. * ------------------------------------------
  242. MOR = MODE LHAUT 'THERMIQUE' 'RAYONNEMENT' 'INFINI' ;
  243. MAR = MATE MOR 'EMIS' 0.8 ;
  244.  
  245.  
  246. * CHARGEMENT DE RAYONNEMENT REPRESENTANT LA TEMPERATURE
  247. * EXTERIEURE (MAINTENUE)
  248. * -----------------------------------------------------
  249. CHTR = MANU 'CHPO' LHAUT 1 'T' T0 ;
  250. CHARAYE = CHAR 'TERA' CHTR EVT2 ;
  251.  
  252.  
  253. * REDEFINITION DE LA TABLE TAB1 POUR PASAPAS
  254. * ------------------------------------------
  255. TAB1 = TABL ;
  256. TAB1 . 'MODELE' = MOT ET MOC ET MOR ;
  257. TAB1 . 'CARACTERISTIQUES' = MAT ET MAC ET MAR ;
  258. TAB1 . 'BLOCAGES_THERMIQUES' = BLT ;
  259. TAB1 . 'CHARGEMENT' = CHT ET CHACONV ET CHARAYE ;
  260. TAB1 . 'TEMPS_CALCULES' = PROG 0. 'PAS' (0.02 * TPSFIN) TPSFIN ;
  261. TAB1 . 'TEMPERATURES' = TABL ;
  262. TAB1 . 'TEMPERATURES' . 0 = T_INI ;
  263. TAB1 . 'CELSIUS' = VRAI ;
  264.  
  265.  
  266. * APPEL A PASAPAS
  267. * ---------------
  268. PASAPAS TAB1 ;
  269.  
  270.  
  271. * POST TRAITEMENT
  272. * ---------------
  273. * EVOLUTION TEMPORELLE DE LA TEMPERATURE
  274. EV11 = EVOL 'ROUG' 'TEMP' TAB1 'TEMPERATURES' 'T' PMIL ;
  275. EV22 = EVOL 'BLEU' 'TEMP' TAB1 'TEMPERATURES' 'T' PG ;
  276. TL = TABL ;
  277. TL . 1 = 'TIRR' ;
  278. TL . 2 = 'TIRR' ;
  279. TL . 'TITRE' = TABL ;
  280. TL . 'TITRE' . 1 = 'PMIL cond seule' ;
  281. TL . 'TITRE' . 2 = 'PG cond seule' ;
  282. TL . 'TITRE' . 3 = 'PMIL cond+conv+ray' ;
  283. TL . 'TITRE' . 4 = 'PG cond+conv+ray' ;
  284. SI ITRAC ;
  285. DESS (EV1 ET EV2 ET EV11 ET EV22) 'LEGE' TL 'TITR' '[4] Temperature vs Temps' ;
  286. FINSI ;
  287.  
  288. * BOUCLE SUR TOUS LES PAS DE TEMPS ET TRACE DU CHAMP DE TEMPERATURES
  289. N1 = DIME (TAB1 . 'TEMPERATURES') ;
  290. SI ITRAC ;
  291. REPE B1 N1 ;
  292. T_I = TAB1 . 'TEMPERATURES' . (&B1 - 1) ;
  293. TPS_I = TAB1 . 'TEMPS' . (&B1 - 1) ;
  294. MOT_I = CHAI '[4] Temperature au temps ' TPS_I ' (conv. + ray.)' ;
  295. TRAC T_I SU CSU 'TITR' MOT_I LISO1 ;
  296. FIN B1 ;
  297. FINSI ;
  298.  
  299. * BOUCLE SUR TOUS LES PAS DE TEMPS ET TRACE VECTEUR FLUX DE CHALEUR
  300. SI ITRAC ;
  301. REPE B1 N1 ;
  302. T_I = TAB1 . 'TEMPERATURES' . (&B1 - 1) ;
  303. VF_I = VECFLU T_I MOT MAT 'COUL' ;
  304. TPS_I = TAB1 . 'TEMPS' . (&B1 - 1) ;
  305. MOT_I = CHAI '[4] Vecteur flux de chaleur au temps ' TPS_I ;
  306. TRAC VF_I CSU 'TITR' MOT_I ;
  307. FIN B1 ;
  308. FINSI ;
  309.  
  310.  
  311. * SAUVEGARDE DES OBJETS CREES
  312. * ---------------------------
  313. OPTI 'SAUV' 'formation_debutant_2_thermique.sauv' ;
  314.  
  315.  
  316. FIN ;
  317.  
  318.  
  319.  
  320.  
  321.  
  322.  
  323.  

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