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

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