Télécharger tran2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : tran2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/06/24 PROVENANCE : DELA
  5.  
  6. *******************************************
  7. * Test tran2.dgibi: jeux de données *
  8. * --------------------------------- *
  9. *******************************************
  10.  
  11. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  12. * SINON SI GRAPH DIFFERENT DE N TOUS
  13. * LES GRAPHIQUES SONT AFFICHES
  14.  
  15. GRAPH = 'N' ;
  16.  
  17. SAUT PAGE;
  18. SI (NEG GRAPH 'N') ;
  19. OPTI ECHO 1 ;
  20. OPTI TRAC X ;
  21. SINO ;
  22. OPTI ECHO 1 ;
  23. FINSI ;
  24.  
  25. SAUT PAGE;
  26.  
  27. *-----------------------------------------*
  28. * *
  29. * TEST TRAN2 *
  30. * *
  31. *--- Problème : *
  32. * *
  33. * Problème de Valentin : *
  34. * thermique transitoire linéaire *
  35. * *
  36. * *
  37. *--- Description de la géométrie : *
  38. * *
  39. * A ----- B *
  40. * | | *
  41. * | | AB = 1. *
  42. * | | BC = 10. *
  43. * | | EC = 4. *
  44. * | | *
  45. * | | *
  46. * F |-----| E *
  47. * | | *
  48. * D | | C *
  49. * X------ -----> R *
  50. * R=0. *
  51. * *
  52. * Conductivité isotrope : K = 1.*
  53. * Capacité du matériau : C = 1.*
  54. * RHO= 1.*
  55. * Sur AB , température imposée: T = 1.*
  56. * Sur BC , convection forcée : TE= 1.*
  57. * H= 10.*
  58. * Dans CDFE , source volumique: Q = 1.*
  59. * Ailleurs, flux imposé : F = 0.*
  60. * *
  61. *--- Calcul : *
  62. * Calcul axisymétrique *
  63. * Recherche de la température maximale*
  64. * dans le cylindre épais par une *
  65. * méthode pas à pas. L'évolution de la*
  66. * température est comparée avec la *
  67. * solution analytique *
  68. * *
  69. *--- Résultats : *
  70. * *
  71. * Comparer avec le cas-test THERMX *
  72. * 'PLI45' *
  73. * La température maximale doit tendre *
  74. * vers 1.3 *
  75. * *
  76. *-----------------------------------------*
  77.  
  78. TITR 'PROBLEME DE VALENTIN - THERMIQUE
  79. TRANSITOIRE LINEAIRE' ;
  80. OPTI DIME 2 ELEM TRI6 MODE AXIS ;
  81. TEMPS ;
  82.  
  83. *------- Création de la géométrie --------*
  84. PD = 0. 0. ;
  85. PC = 1. 0. ;
  86. PB = 1. 10. ;
  87. PA = 0. 10. ;
  88. PE = 1. 4. ;
  89. PF = 0. 4. ;
  90. D1 = DROITE 4 PD PC ;
  91. D2 = DROITE 5 PC PE ;
  92. D3 = DROITE 4 PE PF ;
  93. D4 = DROITE 5 PF PD ;
  94. LIG1 = D1 ET D2 ET D3 ET D4 ;
  95. D5 = DROITE 8 PF PA ;
  96. D6 = DROITE 4 PA PB ;
  97. D7 = DROITE 8 PB PE ;
  98. LIG2 = D3 ET D5 ET D6 ET D7 ;
  99. SURF1 = SURF LIG1 PLANE ;
  100. SURF2 = SURF LIG2 PLANE ;
  101. SURF3 = SURF1 ET SURF2 ;
  102.  
  103. SI(NEG GRAPH 'N');
  104. TRAC 'QUAL' SURF3 ;
  105. FINSI;
  106.  
  107. *--- DONNEES DU PROBLEME DE THERMIQUE ----*
  108.  
  109. *------ 1 - Modélisation du maillage SURF3:
  110. * formulation + matériau
  111. MOD1 = MODE SURF3 THERMIQUE ISOTROPE ;
  112.  
  113. *------ 2 - Création d'un CHAMELEM de
  114. * sous-type CARACTERISTIQUES
  115. MAT1 = MATE MOD1 K 1. 'C' 1. RHO 1. ;
  116.  
  117. *------ 3 - Condition de source volumique :
  118. * calcul du CHPOINT des flux
  119. * nodaux équivalents
  120. S1 = SOURCE MOD1 1. SURF1 ;
  121.  
  122. *------ 4 - Modèle pour la formulation
  123. * convection et type d'élément
  124. MOD2 = MODE (D2 ET D7) THERMIQUE
  125. CONVECTION;
  126. MOD4 = MOD1 ET MOD2;
  127.  
  128. *------ 5 - Création de CHAMELEM de
  129. * caractéristiques pour la
  130. * convection
  131. MAT2 = MATE MOD2 H 10. ;
  132. MAT4 = MAT1 ET MAT2 ;
  133.  
  134. *------ 6 - Blocage des températures
  135. B1 = BLOQ T D6 ;
  136.  
  137. *------ 7 - Valeur de T imposée
  138. E1 = DEPI B1 1. ;
  139.  
  140. *------ 8 - Condition de flux nul : calcul
  141. * du CHPOINT de flux équivalent
  142. FF1 = FLUX MOD1 0. D1 ;
  143.  
  144. *------ 9 - Assemblage des CHPOINTs de flux
  145. * équivalents
  146. QQQ1 = S1 ET FF1 ;
  147.  
  148. *------ 10 - chargement
  149. D69 = D7 'ET' D2;
  150. CH_TE = 'MANU' 'CHPO' 1 D69 'T' 1.;
  151.  
  152. PP1 = PROG 0. PAS 0.05 NPAS 50 ;
  153. PP2 = PROG 51. * 1. ;
  154. MEVOL1 = EVOL MANU 'TEMPS' PP1 'COEFF' PP2;
  155. MEVOL2 = EVOL MANU 'TEMPS' (PROG 0. 2.) 'T' (PROG 1. 1.);
  156. CHA1 = CHAR 'Q ' QQQ1 MEVOL1 ;
  157. CHA2 = CHAR 'TIMP' E1 MEVOL2;
  158. CHA3 = CHAR 'TECO' CH_TE MEVOL2;
  159. CHA4 = CHA1 ET CHA2 ET CHA3;
  160.  
  161. * ----- 11 - remplissage pour PASAPAS
  162. CHTER0 = MANU CHPO SURF3 1 T 1. ;
  163.  
  164. ETAB = TABL ;
  165. ETAB.'TEMPERATURES' = TABL ;
  166. ETAB.'TEMPERATURES' . 0 = CHTER0 ;
  167. ETAB.'MODELE' = MOD4 ;
  168. ETAB.'CARACTERISTIQUES' = MAT4 ;
  169. ETAB.'BLOCAGES_THERMIQUES' = B1 ;
  170. ETAB.'CHARGEMENT' = CHA4 ;
  171. ETAB.'PROCEDURE_THERMIQUE' = LINEAIRE ;
  172. ETAB.'TEMPS_CALCULES' =
  173. PROG 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
  174. 1. 1.1 1.2 1.3;
  175.  
  176.  
  177. PASAPAS ETAB ;
  178.  
  179. LISTT = PROG ;
  180. TEMP1 = PROG ;
  181. REPETER BOUCSOR (DIME ETAB.'TEMPERATURES');
  182. ITERA = &BOUCSOR;
  183. LISTT = INSE LISTT ITERA
  184. (ETAB.'TEMPS'. (ITERA-1)) ;
  185. TEMP1 = INSE TEMP1 ITERA
  186. ( MAXI ( ETAB.'TEMPERATURES'. (ITERA-1) ));
  187. FIN BOUCSOR ;
  188.  
  189. *----- 12 - SORTIE DES RESULTATS
  190. EV1 = EVOL BLEU MANU 'TEMPS' LISTT
  191. 'TEMPERATURE' TEMP1 ;
  192. TIME = PROG 0.1 0.200 0.300 0.400 0.500
  193. 0.600 0.700 0.800 0.900 1.000
  194. 1.100 1.200 1.300 ;
  195. LTEM = PROG 1.076 1.135 1.179 1.212 1.236
  196. 1.253 1.266 1.276 1.282 1.287
  197. 1.290 1.293 1.295 ;
  198. EV2 = EVOL ROUGE MANU TIME LTEM ;
  199.  
  200. SI(NEG GRAPH 'N');
  201. 'EVOLUTION DE LA TEMPERATURE MAXIMALE';
  202. DESS (EV1 ET EV2) ;
  203. FINSI;
  204.  
  205. TEMPS ;
  206.  
  207. * CODE DE FONCTIONNEMENT BASé SUR UN ECART
  208. * MILIEU DE COURBE
  209. * TEST RUDIMENTAIRE POUVANT ETRE OPTIMISé
  210. TREF = 1.253;
  211. ITER1A= 6;
  212. TMAX = MAXI (ETAB.'TEMPERATURES'.ITER1A);
  213. RESI = 100. * (ABS((TMAX-TREF)/TREF));
  214. SI(RESI <EG 5.);
  215. ERRE 0;
  216. SINO;
  217. ERRE 5;
  218. FINSI;
  219.  
  220.  
  221. FIN;
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  

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