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 VIMT CD >s`an style="color: #cc66cc;">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