Télécharger tran10.dgibi

Retour à la liste

Numérotation des lignes :

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

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