Télécharger la plaquette Cast3M


fr en

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