Télécharger la plaquette Cast3M


fr en

Télécharger ther8.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ther8.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/06/13     PROVENANCE : TEST                                    
  5. SAUT PAGE;                                                                      
  6. *                                                                               
  7. ***************************************************************                 
  8. *                                                             *                 
  9. *   TEST THER8                                                *                 
  10. *                                                             *                 
  11. *   TRANSFERT DE CHALEUR AVEC CONVECTION EN 2D                *                 
  12. *                                                             *                 
  13. *   Test NAFEMS numero T4                                     *                 
  14. *                                                             *                 
  15. *                                                             *                 
  16. *  description                                                *                 
  17. *  -----------                                                *                 
  18. *                                                             *                 
  19. *            <- 0.6 m ->                                      *                 
  20. *           D           C  _____                              *                 
  21. *            .----------.       |                             *                 
  22. *           /|          |       |                             *                 
  23. *           /|          |       |                             *                 
  24. *           /|          |       |                             *                 
  25. *           /|          |       |                             *                 
  26. *           /|          |     1.0 m                           *                 
  27. *           /|          |       |                             *                 
  28. *           /|          | __    |                             *                 
  29. *           /|          |  |    |                             *                 
  30. *           /|          | 0.2 m |                             *                 
  31. *            .----------. _|____|                             *                 
  32. *           A           B                                     *                 
  33. *                                                             *                 
  34. *                                                             *                 
  35. *    conditions aux limites                                   *                 
  36. *    ----------------------                                   *                 
  37. *                                                             *                 
  38. *       - temperature imposee sur la droite AB :              *                 
  39. *                     T = 100 degres Celcius                  *                 
  40. *                                                             *                 
  41. *       - convection sur droites BC et CD :                   *                 
  42. *           Temperature exterieure de 0. degres Celcius       *                 
  43. *                                                             *                 
  44. *       - droite DA isolee : pas de convection                *                 
  45. *                                                             *                 
  46. *    On cherche la temperture au point E. La valeur trouvee   *                 
  47. *    est ensuite comparee a la valeur theorique.              *                 
  48. *                                                             *                 
  49. ***************************************************************                 
  50. *                                                                               
  51. OPTI ECHO 0;                                                                    
  52. *                                                                               
  54. OPTI ELEM QUA8;                                                                 
  55. *                                                                               
  56. GRAPH = 'N';                                                                    
  57. *                                                                               
  58. *-------------------- geometrie : maillage --------------------*                
  59. *                                                                               
  60. PA = 0.  0.;                                                                    
  61. PB = 0.6 0.;                                                                    
  62. PC = 0.6 1.0;                                                                   
  63. PD = 0.  1.0;                                                                   
  64. *                                                                               
  65. N1 = 5;                                                                         
  66. N2 = 3;                                                                         
  67. D1 = PA DROI N2 PB;                                                             
  68. D2 = PB DROI N1 PC;                                                             
  69. D3 = PC DROI N2 PD;                                                             
  70. D4 = PD DROI N1 PA;                                                             
  71. *                                                                               
  72. SUR1 = DALL D1 D2 D3 D4 PLAN;                                                   
  73. *                                                                               
  74. SI (NEG GRAPH 'N');                                                             
  75.   TITR 'THER8 : MAILLAGE';                                                      
  76.   TRAC QUAL SUR1;                                                               
  77. FINSI;                                                                          
  78. *                                                                               
  79. *------------------ modeles - materiaux -----------------------*                
  80. *                                                                               
  81. MODL1 = MODE SUR1 THERMIQUE ISOTROPE QUA8;                                      
  82. MATR1 = MATE MODL1 K 52.0;                                                      
  83. *                                                                               
  84. MODL2 = MODE (D2 ET D3) thermique CONVECTION;                                             
  85. MATR2 = MATE MODL2 H 750.0;                                                     
  86. *                                                                               
  87. *---------------- matrices de conductivite --------------------*                
  88. *                                                                               
  89. COND1 = CONDUCTIVITE MODL1 MATR1;                                               
  90. COND2 = CONDUCTIVITE MODL2 MATR2;                                               
  91. *                                                                               
  92. *------------------- temperatures imposees --------------------*                
  93. *                                                                               
  94. *    - Cote AB (droite D1) : temperature constante                              
  95. *                             de 100 degres celcius                             
  96. *                                                                               
  97. BLOCD1 = BLOQ D1 'T';                                                           
  98. *                                                                               
  99. TEMPD1 = DEPI BLOCD1 100.;                                                      
  100. *                                                                               
  101. *-------- flux equivalent a la condition de convection --------*                
  102. *                                                                               
  103. * Temperature ambiante de 0. degres Celcius                                     
  104. *                                                                               
  105. CNVBC = CONVECTION MODL2 MATR2 T 0.;                                            
  106. *                                                                               
  107. *------------------ assemblage des flux -----------------------*                
  108. *                                                                               
  109. FLU1 = TEMPD1 ET CNVBC;                                                         
  110. *                                                                               
  111. *---------------- assemblage des conductivites ----------------*                
  112. *                                                                               
  113. COND3 = COND1 ET COND2 ET BLOCD1;                                               
  114. *                                                                               
  115. *-------------------- resolution ------------------------------*                
  116. *                                                                               
  117. CHPO1 = RESO COND3 FLU1;                                                        
  118. *                                                                               
  119. *---------------- extraction des resultats --------------------*                
  120. *                                                                               
  121. * Trace facultatif de la repartition de temperature                             
  122. *                                                                               
  123. SI (NEG GRAPH 'N');                                                             
  124.    TITR 'THER8 : TEMPERATURE';                                                  
  125.    TRAC SUR1 CHPO1;                                                             
  126. FINSI;                                                                          
  127. *                                                                               
  128. PE = SUR1 POIN PROC (0.6 0.2);                                                  
  129. THETPE2 = EXTR CHPO1 T PE;                                                      
  130. THETPE1 = 18.3;                                                                 
  131. *                                                                               
  132. ERG = 100 *(ABS ((THETPE2 - THETPE1) / THETPE1));                               
  133. *                                                                               
  134. *---------------- affichage des resultats ----------------------*               
  135. *                                                                               
  136. MESS ' RESULTATS ';                                                             
  137. MESS ' --------- ';                                                             
  139. *                                                                               
  140. MESS '    Temperature theorique en E : ' THETPE1;                               
  141. MESS '    Temperature calculee en E  : ' THETPE2;                               
  142. MESS '       Soit un ecart de : 'ERG '%';                                       
  144. *                                                                               
  145. *---------------- code fonctionnement --------------------------*               
  146. *                                                                               
  147. SI (ERG &lt;EG 5);                                                                 
  148.   ERRE 0;                                                                       
  149. SINON;                                                                          
  150.   ERRE 5;                                                                       
  151. FINSI;                                                                          
  153. *                                                                               
  154. TEMPS;                                                                          
  155. *                                                                               
  156. FIN;                                                                            
  157. *                                                                               
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  

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