Télécharger ccar3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ccar3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **
  5. ** --- 10 OCTOBRE 1997 ---
  6. **
  7. ** TEST CAVITE CARRE
  8. **
  9. ** teste LAPL KONV(SUPG) KMAB KMBT en IMPL et DFDT (centrep0 et p1)
  10. ** Formulation MACRO P0 et P1
  11. ** Formulation QUADR P0 et P1
  12.  
  13. *opti trace ps ;
  14. OPTION ISOV 'SULI' ;
  15. GRAPH = 'N' ;
  16. COMPLET = FAUX ;
  17.  
  18. SI ( COMPLET ) ;
  19. ds1=0.02 ; ds2=0.2 ;
  20. *ds1=0.02 ; ds2=0.1 ;
  21. ds1=0.01 ; ds2=0.05 ;
  22. SINON ;
  23. err1=2.e-2 ;
  24. ds1=0.04 ; ds2=0.3 ;
  25. FINSI ;
  26.  
  27.  
  28.  
  29. DEBPROC TEST KPRESS*MOT TYPK*MOT MACRO*MOT GRAPH*MOT ;
  30. option dime 2 elem TYPK ;
  31.  
  32. p1= 0 0 ; p12=0.5 0. ; p2= 1 0 ;
  33.  
  34. ab=p1 d dini ds1 dfin ds2 p12 d dini ds2 dfin ds1 p2 ;
  35. mt= ab trans dini ds1 dfin ds2 (0 0.5) trans dini ds2
  36. dfin ds1 (0 0.5) ;
  37. bc=cote 2 mt ; cd=cote 3 mt ; da=cote 4 mt ;
  38. AA=bc moin (0.5 0.) ;
  39. BB=ab plus (0. 0.5) ;
  40. elim (AA et BB et mt) 1.e-3 ;
  41.  
  42. mt= chan mt QUAF ;
  43. $mt=mode mt 'NAVIER_STOKES' MACRO ;
  44.  
  45. EPSS=1.e-8;
  46. MU =1. ;
  47. RO= 400. ;
  48. DT=1. ;
  49.  
  50. * La cavité est fermée il faut imposer la pression en 1 point !
  51. prep1=doma $mt kpress;
  52. bcp=elem prep1 POI1 (lect 1) ;
  53. *
  54. CD=CHAN CD POI1;
  55. CD=ELEM CD (LECT 2 PAS 1 ((NBEL CD) - 1)) ;
  56. *
  57.  
  58. RV= eqex 'OMEGA' 0.7 'NITER' 5
  59. 'OPTI' 'EF' 'IMPL' 'SUPG' KPRESS
  60. ZONE $mt OPER LAPN MU INCO 'UN'
  61. ZONE $mt OPER KONV RO 'UN' MU DT INCO 'UN'
  62. ZONE $mt OPER KMAB (-1.) 1.E-2 INCO 'UN' 'PRES'
  63. ZONE $mt OPER KMBT (-1.) INCO 'PRES' 'UN'
  64. 'OPTI' 'EFMC' INCOD KPRESS
  65. ZONE $mt OPER DFDT EPSS 'PRES' DT 'UN' MU INCO 'PRES'
  66. ;
  67. RV= eqex RV CLIM
  68. PRES TIMP bcp 0.
  69. UN UIMP CD 1. UN VIMP CD 0. UN UIMP DA 0. UN VIMP DA 0.
  70. UN UIMP AB 0. UN VIMP AB 0. UN UIMP BC 0. UN VIMP BC 0. ;
  71.  
  72. rv.inco= table inco ;
  73. rv.inco.un = kcht $mt vect sommet (1.e-5 1.e-5) ;
  74. rv.inco.pres = kcht $mt scal KPRESS 1.e-5 ;
  75.  
  76. exec rv ;
  77.  
  78. AA = chan AA quaf ;
  79. BB = chan BB quaf ;
  80. elim (mt et aa et bb ) 1.e-3 ;
  81. $AA=mode AA 'NAVIER_STOKES' MACRO ;
  82. $BB=mode BB 'NAVIER_STOKES' MACRO ;
  83. srti=doma $AA 'MAILLAGE' ;
  84. srth=doma $BB 'MAILLAGE' ;
  85. evolV = EVOL 'CHPO' (rv.'INCO'.'UN') UX (srti ) ;
  86. evolH = EVOL 'CHPO' (rv.'INCO'.'UN') UY (srth ) ;
  87. evx=extr evolV 'ORDO' ;
  88. list evx ;
  89. evy=extr evolV 'ABSC' ;
  90. evolV= evol 'MANU' 'Vitesse' evx 'Hauteur' evy ;
  91. rv.'EVOLV'=evolV ;
  92. rv.'EVOLH'=evolH ;
  93.  
  94. si ('EGA' graph 'O' );
  95. TAB1=TABLE;
  96. TAB1.'TITRE'=TABLE ;
  97. TAB1 . 1 ='MARQ REGU ' ;
  98. TAB1.'TITRE' . 1 = mot 'Composante_UX ' ;
  99. DESS evolV 'TITX' 'R (m)' 'TITY' 'V (m/s)' LEGE TAB1 ;
  100.  
  101. TAB1.'TITRE' . 1 = mot 'Composante_UY ' ;
  102. DESS evolH 'TITX' 'R (m)' 'TITY' 'V (m/s)' LEGE TAB1 ;
  103. c1=vect (rv.inco.'UN') 0.3 ux uy jaune ;
  104. trace c1 mt ;
  105. pn= ELNO $mt (rv.inco.'PRES') KPRESS ;
  106. trace pn mt ;
  107. mess ' MAX P ' (maxi pn) 'MIN P ' (mini pn) ;
  108. finsi ;
  109. FINPROC RV ;
  110.  
  111.  
  112. RV= test 'CENTRE' 'QUA8' 'MACRO' GRAPH ;
  113. evv= (rv.'EVOLV');
  114. evh= (rv.'EVOLH');
  115. SI ( NON COMPLET ) ;
  116. lrr=prog
  117. 2.43000E-08 +2.99486E-02 +5.66717E-02 +9.69482E-02 +.13621
  118. +.19572 +.23763 +.21778 +.11062 +2.22662E-03
  119. 8.41545E-02 .14460 .21077 .29282 .49443
  120. .71545 .99757 ;
  121. lr=extr (rv.'EVOLV') 'ABSC' ;
  122. lr = abs lr ;
  123. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  124. mess ' Ecart sur CENTRE QUA8 MACRO ' er ;
  125. si ( er > err1 ) ; erreur 5 ; finsi ;
  126. FINSI ;
  127.  
  128. RV= test 'CENTREP1' 'QUA8' 'QUAF' GRAPH ;
  129. evv= evv et (rv.'EVOLV');
  130. evh= evh et (rv.'EVOLH');
  131. SI ( NON COMPLET ) ;
  132. lrr=prog
  133. 2.43000E-08 +2.16450E-02 +4.15949E-02 +7.32300E-02 +.10680
  134. +.16017 +.21000 +.19926 +9.28972E-02 +2.03073E-02
  135. 6.24749E-02 .11395 .18806 .27293 .47598
  136. .70673 .99757 ;
  137. lr=extr (rv.'EVOLV') 'ABSC' ;
  138. lr = abs lr ;
  139. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  140. mess ' Ecart sur CENTREP1 QUA8 QUADR ' er ;
  141. si ( er > err1 ) ; erreur 5 ; finsi ;
  142. FINSI ;
  143.  
  144. RV= test 'CENTREP1' 'QUA8' 'MACRO' GRAPH ;
  145. evv=evv et (rv.'EVOLV');
  146. evh= evh et (rv.'EVOLH');
  147. SI ( NON COMPLET ) ;
  148. lrr=prog
  149. 2.43000E-08 +2.75896E-02 +5.30539E-02 +9.28867E-02 +.13248
  150. +.19214 +.24900 +.23503 +.10435 +1.59022E-02
  151. 9.16737E-02 .15820 .22408 .30042 .50088
  152. .72370 .99757;
  153. lr=extr (rv.'EVOLV') 'ABSC' ;
  154. lr = abs lr ;
  155. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  156. mess ' Ecart sur CENTREP1 QUA8 MACRO ' er ;
  157. si ( er > err1 ) ; erreur 5 ; finsi ;
  158. FINSI ;
  159.  
  160. RV= test 'CENTREP0' 'QUA8' 'QUAF' GRAPH ;
  161. evv= evv et (rv.'EVOLV');
  162. evh= evh et (rv.'EVOLH');
  163. SI ( NON COMPLET ) ;
  164. lrr=prog
  165. 2.43000E-08 +3.29351E-02 +5.56572E-02 +8.29802E-02 +.10512
  166. +.12415 +.12389 +6.57575E-02 +4.16763E-02 +2.59546E-02
  167. +1.09653E-02 +5.26515E-03 5.43831E-02 .15114 .38804
  168. .64474 .99757;
  169. lr=extr (rv.'EVOLV') 'ABSC' ;
  170. lr = abs lr ;
  171. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  172. mess ' Ecart sur CENTREP0 QUA8 QUADR ' er ;
  173. si ( er > err1 ) ; erreur 5 ; finsi ;
  174. FINSI ;
  175.  
  176. RV= test 'CENTREP0' 'QUA8' 'MACRO' GRAPH ;
  177. evv=evv et (rv.'EVOLV');
  178. evh= evh et (rv.'EVOLH');
  179. SI ( NON COMPLET ) ;
  180. lrr=prog
  181. 2.43000E-08 3.88683E-02 6.63672E-02 9.75319E-02 +.12062
  182. .13600 .13108 6.96810E-02 4.52764E-02 2.02565E-02
  183. 4.81434E-03 3.33747E-03 6.00279E-02 .15440 .39504
  184. .65279 .99757;
  185. lr=extr (rv.'EVOLV') 'ABSC' ;
  186. lr = abs lr ;
  187. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  188. mess ' Ecart sur CENTREP0 QUA8 MACRO ' er ;
  189. si ( er > err1 ) ; erreur 5 ; finsi ;
  190. FINSI ;
  191.  
  192. si ('EGA' graph 'O' );
  193. titre 'Coupe oy à x=1/2 ' ;
  194. TAB1=TABLE;
  195. TAB1.'TITRE'= table ;
  196. TAB1 . 1 ='TIRR MARQ CROI REGU ';
  197. TAB1.'TITRE' . 1 ='MACRO STAB QUA4 ';
  198. TAB1 . 2 ='TIRM MARQ PLUS REGU ';
  199. TAB1.'TITRE' . 2 =' P1-Q2 ';
  200. TAB1 . 3 ='TIRR MARQ ETOI REGU ';
  201. TAB1.'TITRE' . 3 =' P1 iso Q2';
  202. TAB1 . 4 ='TIRM MARQ CARR REGU ';
  203. TAB1.'TITRE' . 4 =' P0-Q2 ';
  204. TAB1 . 5 ='TIRR MARQ LOSA REGU ';
  205. TAB1.'TITRE' . 5 =' P0 iso Q2 ';
  206. DESS EVV 'TITX' 'Oy' LEGE TAB1 ;
  207. titre 'Coupe ox à y=1/2 ' ;
  208. DESS EVH 'TITX' 'Ox' LEGE TAB1 ;
  209. FINSI ;
  210.  
  211. RV= test 'CENTRE' 'TRI6' 'MACRO' GRAPH ;
  212. evv= (rv.'EVOLV');
  213. evh= (rv.'EVOLH');
  214. SI ( NON COMPLET ) ;
  215. lrr=prog
  216. 2.43000E-08 3.78634E-02 7.36738E-02 .10991 .15048
  217. .17902 .21175 .19267 .12321 1.24406E-02
  218. 8.49699E-02 .14204 .19784 .27675 .49157
  219. .70911 .99757;
  220. lr=extr (rv.'EVOLV') 'ABSC' ;
  221. lr = abs lr ;
  222. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  223. mess ' Ecart sur CENTRE TRI6 MACR ' er ;
  224. si ( er > err1 ) ; erreur 5 ; finsi ;
  225. FINSI ;
  226.  
  227. RV= test 'CENTREP1' 'TRI6' 'QUAF' GRAPH ;
  228. evv=evv et (rv.'EVOLV');
  229. evh= evh et (rv.'EVOLH');
  230. SI ( NON COMPLET ) ;
  231. lrr=prog
  232. 2.43000E-08 2.63551E-02 5.27287E-02 8.52686E-02 .12436
  233. .17501 .22987 .20393 .11337 1.39453E-02
  234. 8.24727E-02 .13427 .21546 .29014 .49608
  235. .71022 .99757;
  236. lr=extr (rv.'EVOLV') 'ABSC' ;
  237. lr = abs lr ;
  238. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  239. mess ' Ecart sur CENTREP1 TRI6 QUADR ' er ;
  240. si ( er > err1 ) ; erreur 5 ; finsi ;
  241. FINSI ;
  242.  
  243. RV= test 'CENTREP1' 'TRI6' 'MACRO' GRAPH ;
  244. evv=evv et (rv.'EVOLV');
  245. evh= evh et (rv.'EVOLH');
  246. SI ( NON COMPLET ) ;
  247. lrr=prog
  248. 2.43000E-08 2.84128E-02 4.97722E-02 8.16473E-02 .10407
  249. .12931 .14750 9.07169E-02 5.61807E-02 1.24650E-02
  250. 2.90796E-03 2.01045E-02 6.12103E-02 .15432 .40073
  251. .67195 .99757 ;
  252. lr=extr (rv.'EVOLV') 'ABSC' ;
  253. lr = abs lr ;
  254. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  255. mess ' Ecart sur CENTREP1 TRI6 MACRO ' er ;
  256. si ( er > err1 ) ; erreur 5 ; finsi ;
  257. FINSI ;
  258.  
  259. RV= test 'CENTREP0' 'TRI6' 'MACRO' GRAPH ;
  260. evv=evv et (rv.'EVOLV');
  261. evh= evh et (rv.'EVOLH');
  262. SI ( NON COMPLET ) ;
  263. lrr=prog
  264. 2.43000E-08 2.67197E-02 4.69013E-02 7.77459E-02 9.90798E-02
  265. .12505 .14073 9.63472E-02 6.05324E-02 1.80841E-02
  266. 1.57786E-03 2.34364E-02 6.58036E-02 .16034 .40541
  267. .67444 .99757;
  268. lr=extr (rv.'EVOLV') 'ABSC' ;
  269. lr = abs lr ;
  270. ER=(SOMM( abs (lr - lrr) ))*0.0625 ;
  271. mess ' Ecart sur CENTREP0 TRI6 MACRO ' er ;
  272. si ( er > err1 ) ; erreur 5 ; finsi ;
  273. FINSI ;
  274.  
  275.  
  276. si ('EGA' graph 'O' );
  277. titre 'Coupe oy à x=1/2 ' ;
  278. TAB1=TABLE;
  279. TAB1.'TITRE'= table ;
  280. TAB1 . 1 ='TIRR MARQ CROI REGU ';
  281. TAB1.'TITRE' . 1 ='P0 P1 Stab ';
  282. TAB1 . 2 ='TIRM MARQ PLUS REGU ';
  283. TAB1.'TITRE' . 2 =' P1 P2+Bulle ';
  284. TAB1 . 3 ='TIRR MARQ ETOI REGU ';
  285. TAB1.'TITRE' . 3 =' P1 iso P2';
  286. TAB1 . 4 ='TIRM MARQ CARR REGU ';
  287. TAB1.'TITRE' . 4 =' P0 iso P2';
  288. DESS EVV 'TITX' 'Oy' LEGE TAB1 ;
  289. titre 'Coupe ox à y=1/2 ' ;
  290. DESS EVH 'TITX' 'Ox' LEGE TAB1 ;
  291. FINSI ;
  292.  
  293. FIN ;
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  
  313.  
  314.  
  315.  
  316.  

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