Télécharger la plaquette Cast3M


fr en

Télécharger echang.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  echang.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************
  6. *                                                                      *
  7. * Echange par convection bilatéral                                     *
  8. *                                                                      *
  9. * Calcul 3D (ou 2D-plan) linéaire permanent                                         *
  10. *                                                                      *
  11. * on considère 3 plaques parallèles                                    *
  12. *                                                                      *
  13. *    -------- T3                                                       *
  14. *                      face à face                                     *
  15. *                                                                      *
  16. *    -------- T2                                                       *
  17. *                      face à face                                     *
  18. *                                                                      *
  19. *    -------- T1,e1                                                    *
  20. *                                                                      *
  21. * La température des plaques supérieure et inférieure est imposée et   *
  22. * uniforme (respectivement 1000K et 2000K) et on détermine la tempéra- *
  23. * ture d'équilibre de la plaque intermédiaire.                         *
  24. *                                                                      *
  25. ************************************************************************
  26.  
  27. OPTI ECHO 0;
  28.  
  29. xdim3 = VRAI;
  30.  
  31.  
  32. *-------------------------------
  33. * Données physiques du problème
  34. *-------------------------------
  35.  
  36. T1 = 1000. ; T3 = 2000. ;
  37. h12 = 200. ; h23 = 200. ;
  38. lamb = 2.e3 ;
  39.  
  40. *-------------------------------
  41. * Solution analytique
  42. *-------------------------------
  43.  
  44. T2_ana =  ((h12*T1) + (h23*T3))/(h12 + h23) ;
  45.  
  46. *-------------------------------
  47. *   Maillage 
  48. *------------------------------- 
  49.  
  50. * 3D
  51.  
  52. si xdim3 ;
  53. option  dime 3 elem cub8  ;
  54. p1=0 0 0 ; p2=1 0 0 ; p3=1 1 0 ; p4=0 1 0;
  55. n= 3 ;
  56. l1=p1 d n p2;
  57. l2=p2 d n p3;
  58. l3=p3 d n p4;
  59. l4=p4 d n p1;
  60. s1= (l1 et l2 et l3 et l4) surf plan ;
  61.  
  62. n=1 ;
  63. v = 0 0 0.1 ;  v1 = (-1)*v ;
  64. w = 0 0 0.001 ;  w1 = (-1)*w ;
  65.  
  66.  s2inf = s1 ;
  67.  s2sup = s2inf plus w ;
  68.  plaq2 = s2inf volu n s2sup ;
  69.  
  70.  s3inf = s2sup plus v ;
  71.  s3sup = s3inf plus w ;
  72.  plaq3 = s3inf volu n s3sup ;
  73.  
  74.  s1sup = s2inf plus v1;
  75.  s1inf = s1sup plus w1;
  76.  plaq1 = s1INF volu n s1sup ;
  77.  
  78.  rac12 = liai  0.1  s1sup s2inf ;
  79.  rac23 = liai  0.1  s2sup s3inf ;
  80.  
  81.  tout = plaq1 et plaq2 et plaq3 ;
  82.  
  83.  
  84. * 2D
  85.  
  86. sinon ;
  87.   option  dime 2 elem qua4  ;
  88.   p  = 0. 0. ; q  = 1. 0. ; s  = droi 1 p q ;
  89.   s1inf = s ;
  90.   plaq1 = s1inf tran 1 (0. 0.001) ;
  91.   s1sup = inve ( cote 3 plaq1) ;
  92.  
  93.   s2inf = s1inf plus (0. 0.1) ;
  94.   plaq2 = s2inf tran 1 (0. 0.001) ;
  95.   s2sup = inve ( cote 3 plaq2) ;
  96.  
  97.   s3inf = s2sup plus (0. 0.1) ;
  98.   plaq3 = s3inf tran 1 (0. 0.001) ;
  99.  
  100.   rac12 = racc 0.1 s1sup s2inf ;
  101.   rac23 = racc 0.1 s2sup s3inf ;
  102.  
  103.   tout = plaq1 et plaq2 et plaq3 ;
  104.  
  105. finsi ;
  106.  
  107. *------------------------------- 
  108. *  Conduction
  109. *------------------------------- 
  110.  
  111.   mcd = mode tout thermique ;
  112.   k = mate mcd 'K' lamb ;
  113.   cnd = cond mcd k ;
  114.  
  115. *------------------------------- 
  116. * Convection bilatérale
  117. *------------------------------- 
  118.  
  119.   mcv12 = mode rac12 thermique convection ;
  120.   mcv23 = mode rac23 thermique convection ;
  121.  
  122.   ma12 = MATE mcv12 'H' h12 ;
  123.   ma23 = MATE mcv23 'H' h23 ;
  124.  
  125.   cr12 = cond mcv12 ma12    ;
  126.   cr23 = cond mcv23 ma23    ;
  127.  
  128. *------------------------------- 
  129. *  Conditions aux limites
  130. *------------------------------- 
  131.  
  132.     c1 = bloq plaq1   'T' ;
  133.     tim1 = depi c1 T1 ;
  134.  
  135.     c3 = bloq plaq3   'T' ;
  136.     tim3 = depi c3 T3 ;
  137.  
  138. *------------------------------- 
  139. *  Résolution
  140. *------------------------------- 
  141.  
  142.  
  143.        crtot= cr12 et cr23  ;
  144.        cndtot = (crtot et cnd) et (c1 et c3) ;
  145.  
  146.        tt = resou cndtot (tim3 et tim1) ;
  147.  
  148.     Temp = exco 'T' tt 'T' ;
  149.     T2 = redu plaq2   Temp ;
  150.  
  151.  
  152. *------------------------------- 
  153. *  code fonctionnement
  154. *------------------------------- 
  155.    nno = nbno plaq2;
  156.    rr = resul t2 ;
  157.    rra = MAXI ( rr / nno );
  158.    ere = abs  (( rra - T2_ana ) / T2_ana );
  159.    MESS ' ERREUR RELATIVE '   ere;
  160.  
  161.  
  162.    SI (ere  <EG 1.d-8 ) ;
  163.      ERRE 0;
  164.    SINON;
  165.      ERRE 5;
  166.    FINSI;
  167. fin;
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  

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