Télécharger testlqm.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : testlqm.dgibi
  2. *---------------------------------------------------------
  3. * Ce jdd teste les changement LINE -> QUAF -> MACRO
  4. *--------------------------------------------------------
  5. COMPLET = FAUX ;
  6. *COMPLET = VRAI ;
  7. GRAPH = VRAI ;
  8. GRAPH = FAUX ;
  9. 'OPTI' 'TRAC' 'X' ;
  10. *'OPTI' 'TRAC' 'PSC' ;
  11.  
  12.  
  13. *--------------------------------------------------------
  14. * Definition du maillage de l'enceinte cylindrique
  15. *
  16. 'OPTI' 'DIME' 3 'ELEM' 'CUB8' ;
  17.  
  18. hj = 1. ; hc = 4.163 ; nc = 10 ; nz = 6 ;
  19. v1 = (0. 0. hj);
  20.  
  21. epsi = 1.000e-5 ; ;
  22.  
  23.  
  24. p0 = 0.000 0.000 0.000 ;
  25. pox = p0 plus (hj 0. 0.);
  26. pxy = pox plus (0. hj 0.);
  27. poy = p0 plus (0. hj 0.);
  28. poz = p0 plus (0. 0. hj);
  29. q0 = 0.000 0.000 hj ;
  30. qox = p0 plus (hj 0. hj);
  31. qxy = pox plus (0. hj hj);
  32. qoy = p0 plus (0. hj hj);
  33.  
  34.  
  35.  
  36. *2D
  37. q4=manu 'QUA4' p0 pox pxy poy;
  38. t3=manu 'TRI3' p0 pox pxy;
  39. Si Graph ;
  40. trace (q4 et (T3 plus (1.5 0. 0.))) 'NOEU' titr ' LINEAIRE: QUA4 TRI3';
  41. Finsi ;
  42. * On cree les QUAD
  43. qq4=chan quad q4;
  44. qt3=chan quad t3;
  45. Si Graph ;
  46. trace (qq4 et (qT3 plus (1.5 0. 0.))) 'NOEU' titr 'QUAD: QUA8 TRI6';
  47. Finsi ;
  48. * On cree les QUAF à partir de LINEAIRE
  49. m1q4=chan QUAF q4;
  50. m1t3=chan QUAF t3;
  51. Si Graph ;
  52. trace (m1q4 et (m1T3 plus (1.5 0. 0.))) 'NOEU' titr 'QUAF1: QUA9 TRI7';
  53. Finsi ;
  54. * On cree les QUAF à partir de QUAD
  55. m2q4=chan QUAF qq4;
  56. m2t3=chan QUAF qt3;
  57. Si Graph ;
  58. trace (m2q4 et (m2T3 plus (1.5 0. 0.))) 'NOEU' titr 'QUAF2: QUA9 TRI7';
  59. Finsi ;
  60.  
  61. $mod= mode (m2q4 et (m2T3 plus (1.5 0. 0.))) 'NAVIER_STOKES' LINE ;
  62. mf=doma $mod 'MAILLAGE';
  63. ms=(doma $mod 'MSOMMET') coul rouge;
  64. Si Graph ;
  65. trace (mf et ms) 'NOEU' titr '$mod line';
  66. Finsi ;
  67.  
  68. $mod= mode (m2q4 et (m2T3 plus (1.5 0. 0.))) 'NAVIER_STOKES' MACRO;
  69. mf=doma $mod 'MAILLAGE';
  70. ms=(doma $mod 'MSOMMET') coul rouge;
  71. Si Graph ;
  72. trace (mf et ms) 'NOEU' titr '$mod MACRO';
  73. Finsi ;
  74.  
  75. $mod= mode (m2q4 et (m2T3 plus (1.5 0. 0.))) 'NAVIER_STOKES' QUAF;
  76. mf=doma $mod 'MAILLAGE';
  77. ms=(doma $mod 'MSOMMET') coul rouge;
  78. Si Graph ;
  79. trace (mf et ms) 'NOEU' titr '$mod QUAF';
  80. Finsi ;
  81.  
  82.  
  83.  
  84. *3D
  85.  
  86. *Création des éléments linéaires
  87. t4=manu 'TET4' p0 pox pxy poz;
  88. c8=manu 'CUB8' p0 pox pxy poy q0 qox qxy qoy;
  89. c8=c8 plus (1.5 0. 0.);
  90. p6=manu 'PRI6' p0 pox pxy q0 qox qxy;
  91. p6=p6 plus (3. 0. 0.);
  92. p5=manu 'PYR5' p0 pox pxy poy q0;
  93. p5=p5 plus (4.5 0. 0.);
  94. Si Graph ;
  95. trace (t4 et c8 et p6 et p5) 'NOEU' titr 'LINE: TET4 CUB8 PRI6 PYR5';
  96. Finsi ;
  97.  
  98. *Création des éléments QUAD
  99. qad=chan quad (t4 et c8 et p6 et p5);
  100. Si Graph ;
  101. trace qad 'NOEU' titr 'QUAD: TE10 CU20 PR18 PY17';
  102. Finsi ;
  103.  
  104. *Création des éléments QUAF un à un
  105. mc8=chan quaf c8;
  106. mp6=chan quaf p6;
  107. mt4=chan quaf t4;
  108. mp5=chan quaf p5;
  109.  
  110. Si Graph ;
  111. trace (mc8 et mp6 et mt4 et mp5) 'NOEU'
  112. titr 'QUAF : CU27 PR18 TE15 PY19' ;
  113. Finsi ;
  114.  
  115. *Création des éléments MACRO un à un en passant par DOMA
  116.  
  117. $mod= mode mt4 'NAVIER_STOKES' 'MACRO' ;
  118. mf=doma $mod 'MAILLAGE';
  119. ms=(doma $mod 'MSOMMET') coul rouge;
  120. Si Graph ;
  121. trace (mf et ms) 'NOEU' titr '$mod MACRO TE15';
  122. Finsi ;
  123. dv1= (mesu t4) - (mesu mf);
  124. mess 'Volume t4=' (mesu t4) ' Volume mf=' (mesu mf) ' dv1=' dv1;
  125.  
  126. $mod= mode mc8 'NAVIER_STOKES' 'MACRO' ;
  127. mf=doma $mod 'MAILLAGE';
  128. ms=(doma $mod 'MSOMMET') coul rouge;
  129. Si Graph ;
  130. trace (mf et ms) 'NOEU' titr '$mod MACRO CU27';
  131. Finsi ;
  132. dv2= (mesu c8) - (mesu mf);
  133. mess 'Volume c8=' (mesu c8) ' Volume mf=' (mesu mf) ' dv2=' dv2;
  134.  
  135. $mod= mode mp6 'NAVIER_STOKES' 'MACRO' ;
  136. mf=doma $mod 'MAILLAGE';
  137. ms=(doma $mod 'MSOMMET') coul rouge;
  138. Si Graph ;
  139. trace (mf et ms) 'NOEU' titr '$mod MACRO PR21';
  140. Finsi ;
  141. dv3= (mesu p6) - (mesu mf);
  142. mess 'Volume p6=' (mesu p6) ' Volume mf=' (mesu mf) ' dv3=' dv3;
  143.  
  144. $mod= mode mp5 'NAVIER_STOKES' 'MACRO' ;
  145. mf=doma $mod 'MAILLAGE';
  146. ms=(doma $mod 'MSOMMET') coul rouge;
  147. Si Graph ;
  148. trace (mf et ms) 'NOEU' titr '$mod MACRO PY17';
  149. Finsi ;
  150. dv4= (mesu p5) - (mesu mf);
  151. mess 'Volume p5=' (mesu p6) ' Volume mf=' (mesu mf) ' dv4=' dv4;
  152.  
  153.  
  154. dvt=dv1 + dv2 + dv3 + dv4;
  155. mess ' dvt=' dvt;
  156.  
  157. Si (dvt > 1.e-15); erreur 5 ; Finsi ;
  158.  
  159. Si GRAPH;
  160. mf1=elem mf 'PYR5';
  161. mf2=elem mf 'TET4';
  162.  
  163.  
  164. a=(elem mf1 (lect 1)) coul rouge;
  165. b=(elem mf1 (lect 2)) coul verte;
  166. c=(elem mf1 (lect 3)) coul bleue;
  167. d=(elem mf1 (lect 4)) coul jaune;
  168.  
  169. trace (mf et a et b et c et d);
  170.  
  171. Finsi ;
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178. 'FIN' ;
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  

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