Télécharger tran14.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : tran10.dgibi
  2. *
  3. * Copie du cas-test "tran8.dgibi" pour tester
  4. * la donnee a PASAPAS d'une :
  5. * - CAPACITE_CONSTANTE
  6. * - CONDUCTIVITE_CONSTRANTE
  7. * et resolution par la procedure DUPONT
  8.  
  9. * Calcul avec des EF TRI6 (au lieu de QUA4)
  10.  
  11. *******************************************
  12. *******************************************
  13.  
  14. *******************************************
  15. * Test tran10.dgibi: jeux de données *
  16. * --------------------------------- *
  17. *******************************************
  18.  
  19. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  20. * SINON SI GRAPH DIFFERENT DE N TOUS
  21. * LES GRAPHIQUES SONT AFFICHES
  22.  
  23. GRAPH = 'N' ;
  24.  
  25. SAUT PAGE;
  26. SI (NEG GRAPH 'N') ;
  27. OPTI ECHO 1 ;
  28. OPTI TRAC PSC ;
  29. SINO ;
  30. OPTI ECHO 0 ;
  31. FINSI ;
  32.  
  33. SAUT PAGE;
  34.  
  35. *******************************************
  36. * tran10 *
  37. * *
  38. * THERMIQUE TRANSITOIRE LINEAIRE EN 2D *
  39. * *
  40. * Test NAFEMS numero T3 *
  41. * *
  42. * *
  43. * description *
  44. * ----------- *
  45. * *
  46. * A C B *
  47. * <_0.02_> *
  48. * | | | *
  49. * |=================================| *
  50. * | | *
  51. * <______________0.1 m______________> *
  52. * ---> axe X *
  53. * *
  54. * *
  55. * conditions aux limites *
  56. * ---------------------- *
  57. * *
  58. * - temperature imposee au point A : *
  59. * TA = 0 *
  60. * *
  61. * - temperature imposee au point B : *
  62. * TB = 100 * sin (PI * temps / 40) *
  63. * *
  64. *******************************************
  65. OPTI ELEM TRI6 ;
  66.  
  67. *---------- geometrie : maillage ---------*
  68.  
  69. PA1 = 0. 0.; PB1 = 0.1 0. ;
  70. PA2 = 0. 0.01; PB2 = 0.1 0.01;
  71.  
  72. N = 10;
  73. D1 = PA1 DROI N PB1;
  74. D2 = PB1 DROI 1 PB2;
  75. D3 = PB2 DROI N PA2;
  76. D4 = PA2 DROI 1 PA1;
  77.  
  78. SUR1 = DALL D1 D2 D3 D4 PLAN;
  79.  
  80. PC = POIN SUR1 PROC (0.08 0.);
  81.  
  82.  
  83. SI (NEG GRAPH 'N');
  84. TITR 'TRAN10 : MAILLAGE';
  85. TRAC 'QUAL' SUR1;
  86. FINSI;
  87.  
  88. *--------- modeles - materiaux -----------*
  89. MODL1 = MODE SUR1 THERMIQUE ISOTROPE ;
  90. MATR1 = MATE MODL1 'RHO' 7200 'K' 20.0
  91. 'C' 240.5 ;
  92.  
  93. * test de l'option 'CAPACITE_CONSTANTE'
  94. * de PASAPAS :
  95. MATR2 = MATE MODL1 'RHO' 7200 'K' 15.0
  96. 'C' 200.0 ;
  97. CAP2 = CAPA MODL1 MATR2 ;
  98. CND2 = COND MODL1 MATR2 ;
  99.  
  100. * temperatures imposees : fonction du temps
  101. *
  102. * - Cote PA : temperature constante
  103. * de 0 degres celcius
  104. * - Cote PB : temperature variable
  105. * en fonction du temps
  106.  
  107. BLOCD4 = BLOQ D4 'T';
  108. BLOCD2 = BLOQ D2 'T';
  109.  
  110. TEMPD4 = DEPI BLOCD4 1.;
  111. TEMPD2 = DEPI BLOCD2 1.;
  112.  
  113. * Un pas toutes les secondes,
  114. * temps maximum 40. s.
  115. LTEMPS = PROG 0. PAS 1. 40.;
  116.  
  117. LD4 = PROG 41 * 0.;
  118. LD2 = PROG SINU (1. / 80.) AMPL 100 LTEMPS;
  119.  
  120. EVOLD4 = EVOL MANU TEMPS LTEMPS THETA LD4;
  121. EVOLD2 = EVOL MANU TEMPS LTEMPS THETA LD2;
  122.  
  123. CHAD4 = CHAR 'TIMP' TEMPD4 EVOLD4;
  124. CHAD2 = CHAR 'TIMP' TEMPD2 EVOLD2;
  125.  
  126. * Creation d'un flux nul (type chargement)
  127. FLU0 = MANU CHPO (D1 ET D3) 1 Q 0.;
  128.  
  129. EVOL3= EVOL MANU TEMP LTEMPS FLUX
  130. (PROG 41 * 1.);
  131. FLU1 = CHAR 'Q' FLU0 EVOL3;
  132.  
  133. *--- objets pour la procedure PASAPAS ----*
  134. BLOCT = (BLOCD2 ET BLOCD4);
  135. CHART = (FLU1 ET CHAD4 ET CHAD2);
  136.  
  137. TAB1 = TABL;
  138. TAB1.'TEMPERATURES' = TABL;
  139. TAB1.'TEMPERATURES' . 0 =
  140. MANU CHPO SUR1 1 T 0.;
  141. TAB1.'BLOCAGES_THERMIQUES' = BLOCT;
  142. TAB1.'CHARGEMENT' = CHART;
  143. TAB1.'MODELE' = MODL1;
  144. TAB1.'CARACTERISTIQUES' = MATR1;
  145. TAB1.'TEMPS_SAUVES' =
  146. PROG 0. PAS 1. 32.;
  147. TAB1.'TEMPS_CALCULES' =
  148. PROG 0. PAS 1. 32.;
  149. TAB1.'PROCEDURE_THERMIQUE' = DUPONT ;
  150. TAB1.'CAPACITE_CONSTANTE' = CAP2 ;
  151. TAB1.'CONDUCTIVITE_CONSTANTE' = CND2 ;
  152.  
  153. PASAPAS TAB1;
  154.  
  155. *------- extraction des resultats --------*
  156. * Temperature du point C a t = 32 s.
  157.  
  158. * Construction de l'evolution de T au cours
  159. * du temps
  160.  
  161. LTEMPE = VIDE 'LISTREEL';
  162. LTEMPS = VIDE 'LISTREEL';
  163. NBPAS = DIME (TAB1. 'TEMPS') ;
  164. REPE SURPAS NBPAS;
  165. INDICE = &SURPAS;
  166. LTEMPS = LTEMPS ET (PROG
  167. (TAB1.'TEMPS'. (INDICE - 1)));
  168. LTEMPE = LTEMPE ET (PROG
  169. (EXTR TAB1.'TEMPERATURES'. (INDICE - 1) T
  170. (SUR1 POIN PROC (0.08 0.))));
  171. FIN SURPAS;
  172.  
  173. EVTEMPE = EVOL BLEU MANU LTEMPS LTEMPE;
  174.  
  175.  
  176. THET1 = EXTR LTEMPE NBPAS;
  177. THET2 = 36.6;
  178.  
  179. ERG =100 * (ABS ((THET2 - THET1) / THET2));
  180.  
  181. * Trace facultatif de la repartition
  182. * de temperature a t = 32 s
  183. SI (NEG GRAPH 'N');
  184. TITR 'TRAN10 : Temperature a t = 32 s';
  185. CHPO3 = PECHE TAB1 'TEMPERATURES'
  186. TAB1. 'TEMPS' . (NBPAS - 1);
  187. TRAC SUR1 CHPO3;
  188.  
  189. TITR 'Evolution de la temperature en
  190. fonction du temps : ';
  191. DESS EVTEMPE ;
  192. FINSI;
  193.  
  194. *-------- Test Reactions dans solution -------*
  195. IREAC1 = 'EXIS' TAB1 'REACTIONS_THERMIQUES' ;
  196.  
  197. *-------- affichage des resultats --------*
  198.  
  199. MESS ' RESULTATS ';
  200. MESS ' --------- ';
  201.  
  202. MESS 'Temperature theorique :' THET2' C';
  203. MESS 'Temperature calculee :' THET1' C';
  204. MESS ' Soit un ecart de : ' ERG '%';
  205.  
  206. *------- code fonctionnement -------------*
  207. SI ((ERG &lt;EG 5) et IREAC1);
  208. ERRE 0;
  209. SINON;
  210. 'SI' ('NON' IREAC1) ;
  211. 'MESS' '***** ERREUR : il manque les reactions thermiques ! ' ;
  212. 'FINS'
  213. ERRE 5;
  214. FINSI;
  215.  
  216. TEMPS;
  217.  
  218. FIN;
  219.  
  220.  
  221.  
  222.  
  223.  

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