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

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