Télécharger NSmchaml_nonreg.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : NSmchaml_nonreg.dgibi
  2. *** Jeu de données permettant de tester les spg des chaml
  3. *** suivant le type d'élément fini (vitesse/pression) choisi
  4. *** (option navi)
  5. *** Date de création : 16/03/2005
  6. *** Auteur : A. Bleyer
  7.  
  8. 'OPTI' 'DIME' 2 'ELEM' 'QUA8' ;
  9. 'OPTI' 'TRACE' 'X' ;
  10. *'OPTI' 'TRACE' 'PSC' ;
  11.  
  12. P0 = 0.0 0.0 ;
  13. P1 = 0.0 1.0 ;
  14.  
  15. *NX = 1 ;
  16. *NY = 1 ;
  17. NX = 3 ;
  18. NY = 3 ;
  19. EPSI = 1.E-6 ;
  20.  
  21. *D0 = P0 'DROI' (-1*NY) P1 'DINI' 0.1 'DFIN' 0.5 ;
  22. *D1 = P0 'DROI' (-1*NY) P1 'DINI' 0.5 'DFIN' 0.1 ;
  23. D0 = P0 'DROI' NY P1 ;
  24. D1 = P0 'DROI' NY P1 ;
  25. D1 = D1 'PLUS' (1.0 0.0) ;
  26. D2 = D0 'PLUS' (2.0 0.0) ;
  27.  
  28. S1 = D0 'REGL' NX D1 ;
  29. 'OPTI' 'ELEM' 'TRI6' ;
  30. S2 = D1 'REGL' NX D2 ;
  31. S0 = S1 'ET' S2 ;
  32. 'ELIM' S0 EPSI ;
  33.  
  34. M_S0 = 'CHAN' S0 'QUAF' ;
  35. M_S1 = 'CHAN' S1 'QUAF' ;
  36. M_S2 = 'CHAN' S2 'QUAF' ;
  37. 'ELIM' (M_S0 ET M_S1 ET M_S2) EPSI ;
  38.  
  39. LOPT = 'MOTS' 'LICE' 'LIMS' 'MCCE' 'MCP1' 'MCMS' 'QFP1' 'QFMS' ;
  40. ERROR = 0 ;
  41.  
  42. 'OPTI' 'ECHO' 0 ;
  43.  
  44. 'REPETER' BCL1 ('DIME' LOPT) ;
  45.  
  46. *** choix du couple Vitesse/Pression
  47. OPT1 = 'EXTR' LOPT &BCL1 ;
  48. 'OPTI' 'NAVI' OPT1 ;
  49.  
  50. 'MESS' 'Element fini en cours : ' OPT1 ;
  51.  
  52. *** création des objets modèles
  53. $S0 = 'MODE' M_S0 'NAVIER_STOKES' ;
  54. $S1 = 'MODE' M_S1 'NAVIER_STOKES' ;
  55. $S2 = 'MODE' M_S2 'NAVIER_STOKES' ;
  56.  
  57. ******************** traitement SOMMET ****************
  58. ********** champ scalaire à une seule composante ******
  59. MS0 = 'DOMA' $S0 'SOMMET' ;
  60. X1 Y1 = 'COORD' MS0 ;
  61.  
  62. CX1 = 'KCHT' $S0 'SCAL' 'SOMMET' X1 ;
  63. CY1 = 'KCHT' $S0 'SCAL' 'SOMMET' Y1 ;
  64. KX1 = 'CHAN' 'CHAM' CX1 $S0 ;
  65. KY1 = 'CHAN' 'CHAM' CY1 $S0 ;
  66.  
  67. *** verification bon passage SOMMET --> CENTREP1
  68. 'SI' (&BCL1 '>EG' 4) ;
  69. KKX2 = 'CHAN' KX1 $S0 'P1CENTRE' ;
  70. KKY2 = 'CHAN' KY1 $S0 'P1CENTRE' ;
  71.  
  72. SP1 = 'DOMA' $S0 'CENTREP1' ;
  73. XP1 YP1 = 'COORD' SP1 ;
  74.  
  75. CX2 = 'KCHT' $S0 'SCAL' 'CENTREP1' XP1 ;
  76. CY2 = 'KCHT' $S0 'SCAL' 'CENTREP1' YP1 ;
  77. KX2 = 'CHAN' 'CHAM' CX2 $S0 ;
  78. KY2 = 'CHAN' 'CHAM' CY2 $S0 ;
  79.  
  80. VALX_2 = 'MAXI' (KX2 - KKX2) 'ABS' ;
  81. VALY_2 = 'MAXI' (KY2 - KKY2) 'ABS' ;
  82.  
  83. 'MESS' 'SOMMET->CENTREP1 ' VALX_2 VALY_2 ;
  84. 'SI' ((VALX_2 '>' EPSI) 'OU' (VALY_2 '>' EPSI)) ;
  85. 'MESS' 'Probleme elt fini : ' OPT1 ' passage SOMMET->CENTREP1' ;
  86. ERROR = ERROR + 1 ;
  87. 'FINSI' ;
  88.  
  89. CKX2 = 'CHAN' KX2 $S0 'NOEUD' ;
  90. CKY2 = 'CHAN' KY2 $S0 'NOEUD' ;
  91.  
  92. VALX_2 = 'MAXI' (KX1 - CKX2) 'ABS' ;
  93. VALY_2 = 'MAXI' (KY1 - CKY2) 'ABS' ;
  94.  
  95. 'MESS' 'SOMMET<-CENTREP1 ' VALX_2 VALY_2 ;
  96. 'SI' ((VALX_2 '>' EPSI) 'OU' (VALY_2 '>' EPSI)) ;
  97. 'MESS' 'Probleme elt fini : ' OPT1 ' passage SOMMET<-CENTREP1' ;
  98. ERROR = ERROR + 1 ;
  99. 'FINSI' ;
  100.  
  101. 'FINSI' ;
  102.  
  103. *** verification bon passage SOMMET --> CENTRE
  104. 'SI' ((&BCL1 '<' 4) 'OU' (&BCL1 '>' 5)) ;
  105. KKX3 = 'CHAN' KX1 $S0 'GRAVITE' ;
  106. KKY3 = 'CHAN' KY1 $S0 'GRAVITE' ;
  107.  
  108. SC1 = 'DOMA' $S0 'CENTRE' ;
  109. XC1 YC1 = 'COORD' SC1 ;
  110.  
  111. CX3 = 'KCHT' $S0 'SCAL' 'CENTRE' XC1 ;
  112. CY3 = 'KCHT' $S0 'SCAL' 'CENTRE' YC1 ;
  113. KX3 = 'CHAN' 'CHAM' CX3 $S0 ;
  114. KY3 = 'CHAN' 'CHAM' CY3 $S0 ;
  115.  
  116. VALX_3 = 'MAXI' (KX3 - KKX3) 'ABS' ;
  117. VALY_3 = 'MAXI' (KY3 - KKY3) 'ABS' ;
  118.  
  119. 'MESS' 'SOMMET->CENTRE ' VALX_3 VALY_3 ;
  120. 'SI' ((VALX_3 '>' EPSI) 'OU' (VALY_3 '>' EPSI)) ;
  121. 'MESS' 'Probleme elt fini : ' OPT1 ' passage SOMMET->CENTRE' ;
  122. ERROR = ERROR + 1 ;
  123. 'FINSI' ;
  124.  
  125. 'FINSI' ;
  126.  
  127. *** verification bon passage SOMMET --> MSOMMET
  128. 'SI' (&BCL1 '>' 3) ;
  129. KKX4 = 'CHAN' KX1 $S0 'MSOMMET' ;
  130. KKY4 = 'CHAN' KY1 $S0 'MSOMMET' ;
  131.  
  132. SM1 = 'DOMA' $S0 'MSOMMET' ;
  133. XM1 YM1 = 'COORD' SM1 ;
  134.  
  135. CX4 = 'KCHT' $S0 'SCAL' 'MSOMMET' XM1 ;
  136. CY4 = 'KCHT' $S0 'SCAL' 'MSOMMET' YM1 ;
  137. KX4 = 'CHAN' 'CHAM' CX4 $S0 ;
  138. KY4 = 'CHAN' 'CHAM' CY4 $S0 ;
  139.  
  140. VALX_4 = 'MAXI' (KX4 - KKX4) 'ABS' ;
  141. VALY_4 = 'MAXI' (KY4 - KKY4) 'ABS' ;
  142.  
  143. 'MESS' 'SOMMET->MSOMMET ' VALX_4 VALY_4 ;
  144. 'SI' ((VALX_4 '>' EPSI) 'OU' (VALY_4 '>' EPSI)) ;
  145. 'MESS' 'Probleme elt fini : ' OPT1 ' passage SOMMET->MSOMMET' ;
  146. ERROR = ERROR + 1 ;
  147. 'FINSI' ;
  148.  
  149. CKX4 = 'CHAN' KX4 $S0 'NOEUD' ;
  150. CKY4 = 'CHAN' KY4 $S0 'NOEUD' ;
  151.  
  152. VALX_4 = 'MAXI' (KX1 - CKX4) 'ABS' ;
  153. VALY_4 = 'MAXI' (KY1 - CKY4) 'ABS' ;
  154.  
  155. 'MESS' 'SOMMET<-MSOMMET ' VALX_4 VALY_4 ;
  156. 'SI' ((VALX_4 '>' EPSI) 'OU' (VALY_4 '>' EPSI)) ;
  157. 'MESS' 'Probleme elt fini : ' OPT1 ' passage SOMMET<-MSOMMET' ;
  158. ERROR = ERROR + 1 ;
  159. 'FINSI' ;
  160.  
  161. 'FINSI' ;
  162.  
  163.  
  164.  
  165. 'FIN' BCL1 ;
  166.  
  167.  
  168. 'SI' (ERROR '>' 0) ;
  169. 'MESS' ERROR ' Probleme(s)' ;
  170. 'MESS' 'Necessite des verifications !' ;
  171. 'ERRE' 5 ;
  172. 'FINSI' ;
  173.  
  174. 'FIN' ;
  175.  
  176. 'OPTI' 'ECHO' 1 'DONN' 5 ;
  177.  
  178.  
  179.  

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