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

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