Télécharger tran8.dgibi

Retour à la liste

Numérotation des lignes :

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

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