Télécharger shearkei.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : shearkei.dgibi
  2. 'OPTI' dime 2 elem QUA4;
  3. 'OPTI' 'EPTR' 6 ;
  4.  
  5. COMPLET = FAUX ;
  6. *COMPLET = VRAI ;
  7. GRAPH = FAUX;
  8. *GRAPH = VRAI;
  9.  
  10.  
  11. DISCR = 'LINE';
  12. DISCR = 'QUAF';
  13. DISCR = 'MACRO';
  14. KPRES = 'CENTREP1';
  15. KSUPG = 'SUPG' ;
  16.  
  17. Ntauo=1. ;
  18. Tadim = 100.;
  19. Tadim = 20. ;
  20.  
  21. IKAS ='L1n';
  22. IKAS ='M02n';
  23.  
  24. IKAS ='M002n';
  25. IKAS ='M002n+';
  26. IKAS ='M002m+';
  27. IKAS ='M002n-';
  28. IKAS ='Q01f';
  29. IKAS ='M01f';
  30. IKAS ='M1n';
  31. IKAS ='Q1n';
  32. IKAS ='Q01f';
  33. IKAS ='M01n';
  34. IKAS ='M05n';
  35.  
  36. *'OPTI' 'TRAC' 'PSC';
  37. *'OPTI' 'FTRA' (chai '/test2/magn/SHEARKEI1/' ikas '.ps');
  38.  
  39. Si(EGA IKAS 'M002n');
  40. Ntauo=0.02;Tadim = 80.;DISCR = 'MACRO';
  41. ALGOKEPS= mots 'Nut' ;
  42. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  43. Finsi ;
  44.  
  45. Si(EGA IKAS 'M002n+');
  46. Ntauo=0.02;Tadim = 100.;DISCR = 'MACRO';
  47. ALGOKEPS= mots 'Nut' ;
  48. KEA= 0.03; KEB=KEA/2. ; EE = 1.e0 ;
  49. Finsi ;
  50.  
  51. Si(EGA IKAS 'M002m+');
  52. Ntauo=0.02;Tadim = 100.;DISCR = 'MACRO';
  53. ALGOKEPS= mots 'Nut' ;
  54. KEA= 0.015; KEB=KEA/2. ; EE = 1.e0 ;
  55. Finsi ;
  56.  
  57. Si(EGA IKAS 'M002n-');
  58. Ntauo=0.02;Tadim = 80.;DISCR = 'MACRO';
  59. ALGOKEPS= mots 'Nut' ;
  60. KEA= 0.003; KEB=KEA/2. ; EE = 1.e0 ;
  61. Finsi ;
  62.  
  63. Si(EGA IKAS 'M02n');
  64. Ntauo=0.2;Tadim = 100.;DISCR = 'MACRO';
  65. ALGOKEPS= mots 'Nut' ;
  66. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  67. Finsi ;
  68.  
  69. Si(EGA IKAS 'M05n');
  70. Ntauo=0.5;Tadim = 100.;DISCR = 'MACRO';
  71. ALGOKEPS= mots 'Nut' ;
  72. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  73. Finsi ;
  74.  
  75. Si(EGA IKAS 'M1n');
  76. Ntauo=1.;Tadim = 100.;DISCR = 'MACRO';
  77. ALGOKEPS= mots 'Nut' ;
  78. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  79. Finsi ;
  80.  
  81. Si(EGA IKAS 'Q1n');
  82. Ntauo=1.;Tadim = 80.;DISCR = 'QUAF';
  83. ALGOKEPS= mots 'Nut' ;
  84. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  85. Finsi ;
  86.  
  87. Si(EGA IKAS 'Q01n');
  88. Ntauo=0.1;Tadim =20.;DISCR = 'QUAF';
  89. ALGOKEPS= mots 'Nut' ;
  90. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  91. *KEA= 0.03; KEB=KEA/2. ; EE = 1.e0 ;
  92. Finsi ;
  93.  
  94. Si(EGA IKAS 'Q01f');
  95. Ntauo=0.1;Tadim =20.;DISCR = 'QUAF';
  96. ALGOKEPS= mots 'Nut' 'Filtre';
  97. Echl=0.01;
  98. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  99. *KEA= 0.03; KEB=KEA/2. ; EE = 1.e0 ;
  100. Finsi ;
  101.  
  102. Si(EGA IKAS 'M01n');
  103. Ntauo=0.1;Tadim = 100.;DISCR = 'MACRO';
  104. ALGOKEPS= mots 'Nut' ;
  105. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  106. *KEA= 0.03; KEB=KEA/2. ; EE = 1.e0 ;
  107. Finsi ;
  108.  
  109. Si(EGA IKAS 'M01f');
  110. Ntauo=0.1;Tadim = 100.;DISCR = 'MACRO';
  111. ALGOKEPS= mots 'Nut' 'Filtre';
  112. Echl=0.001;
  113. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  114. *KEA= 0.03; KEB=KEA/2. ; EE = 1.e0 ;
  115. Finsi ;
  116.  
  117. Si(EGA IKAS 'L1n');
  118. KPRES = 'CENTRE';
  119. Ntauo=0.5;Tadim = 60.;DISCR = 'LINE';
  120. ALGOKEPS= mots 'Nut' ;
  121. KE = 0.01; EE = 1.e0 ;
  122. *KE = 0.03;
  123. Finsi ;
  124.  
  125.  
  126. Si(NON COMPLET);
  127. DISCR = 'LINE'; KPRES = 'CENTRE';
  128. Ntauo=0.1;Tadim = 20.;
  129. ALGOKEPS= mots 'Nut' ;
  130. KEA= 0.01; KEB=KEA/2. ; EE = 1.e0 ;
  131. Finsi ;
  132.  
  133.  
  134. NITMA =Tadim/Ntauo ENTI ;
  135.  
  136. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  137. *
  138. * CALCUL D'UNE COUCHE DE MELANGE
  139. * ------------------------------
  140. *
  141.  
  142. Ua = 42.8 ;
  143. Ub = 25.2 ;
  144. DeltaU= Ua - Ub ;
  145. DeltaU2= DeltaU * DeltaU ;
  146. Tauo = 0.03 / DeltaU ;
  147. DT = Ntauo * Tauo;
  148. mess '--------------------------------------- ' ;
  149. mess 'Ua=' Ua ' Ub=' Ub 'DELTAU=' DELTAU
  150. 'DELTAU2=' DELTAU2 ' Tauo=' Tauo ;
  151. mess ' Tauo = ' Tauo ' DT = ' DT ;
  152. mess ' NITMA= ' NITMA ;
  153. mess '--------------------------------------- ' ;
  154. *-------------------------------------------------------------------
  155.  
  156.  
  157. LX = 0.75;
  158. LY0= 0.18;
  159. LY0= 0.09;
  160. LY =LY0*2.;
  161.  
  162. p0=0. 0. ;
  163. p1=0. LY0;
  164. p2=0. ((-1.)*LY0);
  165. q0=LX 0.;
  166. q1=LX LY0;
  167. q2=LX ((-1.)*LY0);
  168.  
  169. dn1=0.03*LY0; dn2=0.3*LY0;
  170.  
  171. ent1=p0 d dini dn1 dfin dn2 p1 ;
  172. ent2 = ent1 syme 'POINT' p0;
  173. m1= ent1 trans dini dn1 dfin dn2 (q1 moins p1);
  174. m2= m1 syme 'DROITE' p0 q0 ;
  175. mf = m1 et m2;
  176. elim mf 1.e-5;
  177. *trace mf ;
  178. *
  179. *
  180. haut = cote 2 m1;
  181. bas = cote 2 m2 ;
  182. axe = cote 4 m2 ;
  183. elim (mf et ent2) 1.e-5 ;
  184. *trace (ent1 et ent2 et haut et bas et axe);
  185. Mmf = chan QUAF mf ;
  186. Ment1 = chan QUAF ent1;
  187. Ment2 = chan QUAF ent2;
  188. Mhaut = chan QUAF haut;
  189. Mbas = chan QUAF bas ;
  190. Maxe = chan QUAF axe;
  191. elim (Mmf et Ment1 et Ment2 et Mhaut et Mbas et MAxe)1.e-5;
  192.  
  193. $mf = mode Mmf 'NAVIER_STOKES' DISCR;
  194. $ent1=mode Ment1 'NAVIER_STOKES' DISCR;
  195. ent1 = doma $ent1 maillage ;
  196. $ent2=mode Ment2 'NAVIER_STOKES' DISCR;
  197. ent2 = doma $ent2 maillage ;
  198.  
  199. $bas=mode Mbas 'NAVIER_STOKES' DISCR;
  200. bas = doma $bas maillage ;
  201. $haut=mode Mhaut 'NAVIER_STOKES' DISCR;
  202. haut= doma $haut maillage ;
  203. $axe=mode Maxe 'NAVIER_STOKES' DISCR;
  204. axe = doma $axe maillage inve ;
  205.  
  206. doma $mf 'IMPR';
  207. mf = doma $mf maillage;
  208. cmf = cont mf ;
  209.  
  210. entt=ent1 et ent2 ;
  211. enta = diff (chan ent2 poi1) (chan entt poi1);
  212. entb = diff (chan ent1 poi1) (chan entt poi1);
  213.  
  214. axe1 = chan POI1 axe ;
  215. cflstab = (coor 1 (poin (elem axe1 (lect 2))1))/Ua;
  216. CFL = DT/cflstab;
  217. CFLI = ENTI CFL ;
  218. mess '--------------------------------------- ' ;
  219. mess '--------------------------------------- ' ;
  220. mess ' CFLstab = ' cflstab ' CFL = ' CFL ;
  221. mess '--------------------------------------- ' ;
  222. mess '--------------------------------------- ' ;
  223.  
  224. nba = nbno axe1;
  225. tbzone= table ;
  226. nz=0 ;
  227. tbaret= table ;
  228. tbaret.1 = entt;
  229. repeter BCL (nba-1);
  230. pi = poin axe1 (&bcl+1);
  231. entti=entt plus (pi moins p0);
  232. elim (mf et entti) 1.e-5 ;
  233. tbaret.(&bcl+1) = entti;
  234. *trace (axe et entti);
  235. Si ((coor 1 pi) > 0.35);
  236. nz = nz+1;
  237. tbzone.nz= entti;
  238. Finsi ;
  239. FIN BCL ;
  240. xaxe = evol chpo (coor 1 axe) axe;
  241. xaxe = extr xaxe ABSC;
  242.  
  243. NU =1.E-5;
  244. cnu = 0.09 ;
  245.  
  246. MUF = KCHT $mf SCAL SOMMET 1.e-2;
  247.  
  248. RV= EQEX $mf 'ITMA' NITMA 'OMEGA' 1. 'NITER' 1 'FIDT' 20
  249. 'OPTI' 'IMPL' EF KSUPG KPRES
  250. 'ZONE' $mf 'OPER' 'KEPSILON' 1. 'UN' NU DT 'INCO' 'KN' 'EN'
  251. 'ZONE' $mf 'OPER' 'NS' 1. 'UN' 'MUF' 'INCO' 'UN'
  252. 'OPTI' 'EFM1'
  253. 'ZONE' $mf 'OPER' 'DFDT' 1. 'UN' DT 'UN' 'MUF' 'INCO' 'UN'
  254. ;
  255. RV= EQEX RV
  256. CLIM 'UN' VIMP (HAUT et BAS et ent1 et ent2 ) 0.
  257. CLIM 'UN' UIMP ENT1 Ua 'UN' UIMP ENTB Ub
  258. * CLIM 'UN' VIMP (ENTA et ENT2 et HAUT et BAS) 0.
  259. * CLIM 'KN' TIMP (ENTA et ENT2 et HAUT et BAS) KE
  260. * CLIM 'EN' TIMP (ENTA et ENT2 et HAUT et BAS) EE
  261. * CLIM 'UN' UIMP ENT1 Ua 'UN' UIMP ENT2 Ub
  262. CLIM 'UN' VIMP (ENTA et ENT2 ) 0.
  263. CLIM 'KN' TIMP ENT1 KEA 'KN' TIMP ENTB KEB
  264. CLIM 'EN' TIMP (ENTA et ENT2 ) EE
  265. ;
  266.  
  267. RVP = EQEX
  268. 'OPTI' 'EF' KPRES
  269. ZONE $mf OPER KBBT (-1.) INCO 'UN' 'PRES'
  270. * CLIM 'PRES' TIMP Sortiep 0.
  271. ;
  272. rvp.'METHINV'.TYPINV=2 ;
  273. rvp.'METHINV'.IMPINV=0 ;
  274. rvp.'METHINV'.NITMAX=250;
  275. rvp.'METHINV'.PRECOND=3 ;
  276. rvp.'METHINV'.RESID =1.e-10 ;
  277. rvp.'METHINV' . 'FCPRECT'=300 ;
  278. rvp.'METHINV' . 'FCPRECI'=300 ;
  279.  
  280. RV.'PROJ'= RVP ;
  281.  
  282. RV.INCO=TABLE 'INCO' ;
  283. RV.INCO.'PRES'=kcht $mf SCAL KPRES 0.;
  284. RV.INCO.'UN'=kcht $mf VECT SOMMET (1. 0.);
  285. RV.INCO.'KN'=kcht $mf SCAL SOMMET KEB ;
  286. RV.INCO.'EN'=kcht $mf SCAL SOMMET EE ;
  287. RV.INCO.'MUF'=kcht $mf SCAL SOMMET NU ;
  288. * RV.INCO.'Ro' =kcht $mf SCAL CENTRE 1. ;
  289. RV.INCO.'Ro' = 1. ;
  290. RV.INCO.'Mu' =kcht $mf SCAL SOMMET NU ;
  291. RV.INCO.'Echl'=Echl;
  292.  
  293. RV.'NOMVI' = 'UN' ;
  294. RV.'ALGO_KEPSILON'=ALGOKEPS;
  295.  
  296. na = 10;
  297. nna=nba/na ;
  298. lh= elem axe1 (lect 1 pas na (nna*na+1)) ;
  299. his = khis 'UN' 1 lh
  300. 'UN' 2 lh
  301. 'KN' lh
  302. 'EN' lh ;
  303. his.'KFIH'=1;
  304. rv.'HIST'=his ;
  305.  
  306. * temp place ;
  307. *option donn 5 ;
  308. temp place ;
  309. exec rv ;
  310. temp place ;
  311.  
  312. Si GRAPH;
  313. ************** Controle sur les historiques *****************************
  314. his=rv.'HIST' ;
  315. dessin his.'TABD' his.'1UN' ;
  316. dessin his.'TABD' his.'2UN' ;
  317. dessin his.'TABD' his.'KN' ;
  318. dessin his.'TABD' his.'EN' ;
  319. ************** Controle sur les historiques *****************************
  320. Finsi;
  321.  
  322.  
  323.  
  324. *************************************************************************
  325. un = rv.inco.'UN';
  326.  
  327. ************** Evolution de l'épaisseur de la couche de mélange **********************
  328. uux = EXCO un 'UX';
  329. dudy = exco (kops 'GRADS' uux $mf) 'UY';
  330.  
  331. deltaw = prog ;
  332. repeter BCL nba;
  333. ai = tbaret.&bcl;
  334. u1 = redu uux ai;
  335. du1dy = redu dudy ai;
  336. delti = ((maxi u1) - (mini u1))/ (maxi du1dy);
  337. deltaw = deltaw et (prog delti);
  338. FIN BCL ;
  339.  
  340. evdelta= evol manu xaxe deltaw ;
  341.  
  342. TAB1 = TABLE ;
  343. tab1.1 = 'MARQ LOSA NOLI' ;
  344. tab1.2 = 'MARQ TRIA REGU' ;
  345. tab1.3 = 'MARQ CROI REGU' ;
  346. tab1.4 = 'MARQ TRIB REGU' ;
  347. tab1.5 = 'MARQ PLUS REGU' ;
  348. tab1.'TITRE' = TABLE ;
  349. tab1.'TITRE'. 1 = 'MOT' 'Expérience CEAT' ;
  350. tab1.'TITRE'. 2 = 'MOT' (chai 'CAST3M K-eps: CFL=' CFLI);
  351.  
  352. dltexpx= prog 0.55 0.7 1. 1.4 1.8 2.35 2.8 3.5 4.15 4.85 5.5 6.25
  353. 6.9 7.6 8.3 9. 10.1 10.4;
  354. dltexpx= dltexpx* (0.8/11.1);
  355. dltexp = prog 1.4 1.5 1.65 1.85 2.05 2.4 2.35 2.7 3.2 3.6 4. 4.35
  356. 4.7 5.3 5.75 6.2 6.95 7.05;
  357. dltexp = dltexp * (0.05/8.65);
  358.  
  359. dlexp = prog 1.4 1.5 1.65 1.85 2.05 2.4 2.35 2.7 3.2 3.6 4. 4.35 4.7
  360. 5.3 5.75 6.2 6.95 7.05;
  361. dlexp = dlexp * (0.05/8.65);
  362. xexp = prog 0.55 0.7 1. 1.4 1.8 2.35 2.8 3.5 4.15 4.85 5.5 6.25 6.9
  363. 7.6 8.3 9. 10.1 10.4 ;
  364. xexp = xexp * (0.8/11.1);
  365. edltexp = evol manu xexp dlexp;
  366. evd = edltexp et evdelta;
  367. ******** FIN * Evolution de l'épaisseur de la couche de mélange **********************
  368.  
  369. * Graphiques
  370. Si GRAPH;
  371.  
  372. dess evd TITR (chai ' Delta_w (m) Tadim=' Tadim ' CFL=' CFLI)
  373. TITX ' X (m) ' TITY 'Delta_w (m)' LEGE tab1;
  374.  
  375. kn = rv.inco.'KN';
  376. ung = vect un 0.00055 ux uy jaune ;
  377. trace ung mf ;
  378.  
  379. *trace rv.inco.'Lw' mf cmf (prog 1.e-2 pas 0.5 10.) TITR 'Lw' ;
  380.  
  381. trace (exco un 'UX') mf TITR ' Composante UX de la vitesse';
  382.  
  383. trace (exco un 'UY') mf cmf TITR ' Composante UY de la vitesse';
  384.  
  385.  
  386. ************** Autosimilarité du profil transversal de vitesse ***********************
  387. ysec = prog ;
  388. delt = prog ;
  389. nbz = dime tbzone;
  390. nbt = dime tbaret;
  391.  
  392. teta= rv.inco.'TKTE';
  393. Ksi = rv.inco.'Ksi';
  394. unxa=(1./DeltaU)*((exco un 'UX') - Ub);
  395. knxa=(1./DeltaU/DeltaU)*kn;
  396. tetaxa=teta;
  397. ksia=ksi;
  398.  
  399.  
  400. repeter BCLZ (dime tbzone);
  401. ar = tbzone.&bclz inve ;
  402. deltawsi=extr (nbt - nbz + &bclz) deltaw ;
  403. taui =(extr (nbt - nbz + &bclz) deltaw)/DeltaU ;
  404.  
  405. ev1 = evol chpo unxa ar;
  406. evk = evol chpo knxa ar;
  407. evt = evol chpo tetaxa ar;
  408. evksi = evol chpo ksia ar;
  409. ysec = (extr ev1 'ABSC');
  410. ny= dime ysec ;
  411. ysec0 = prog ny*LY0;
  412. yseci = (ysec - ysec0)/deltawsi;
  413. uxa=extr ev1 'ORDO';
  414. kxa=extr evk 'ORDO';
  415. txa=(extr evt 'ORDO')/(prog ny*taui);
  416. ksixa=extr evksi 'ORDO';
  417. eva = evol manu yseci uxa;
  418. evak= evol manu yseci kxa;
  419. evat= evol manu yseci txa;
  420. evaksi= evol manu yseci ksixa;
  421.  
  422. Si (EGA &bclz 1);
  423. evux = eva;
  424. evkx = evak;
  425. evtx = evat;
  426. evksix = evaksi;
  427. Sinon ;
  428. evux = evux et eva;
  429. evkx = evkx et evak;
  430. evtx = evtx et evat;
  431. evksix = evksix et evaksi;
  432. Finsi ;
  433. FIN BCLZ;
  434.  
  435. kex = prog 5. 13.5 52. 116. 166. 166. 115. 50. 2.;
  436. xkex = prog -37. -29. -21. -12. -3.5 3.5 12. 21. 45.;
  437. kex = kex * (0.5e-2/183.*7.);
  438. xkex = xkex* (3./114.);
  439. evkex = evol manu xkex kex ;
  440. xuex1= prog 65. 56. 50. 44. 38. 35. 31.5 28.5 25. 22. 19. 16. 13. 10. 6.5
  441. 3. 0. ;
  442. xuex1=xuex1*(-3./115.);
  443. uex1=prog 181. 181. 181.5 182.5 182. 182. 181.5 181. 175. 172. 166. 156.5
  444. 145. 133.5 119. 108.5 95.;
  445. uex1 = ((prog 17*183.) - uex1)*(1./183.);
  446. xuex2= prog 3. 5.5 9.5 12.5 15. 18. 22. 25. 28. 31. 34.5 37. 41. 44.
  447. 47. 50. 53. 56. ;
  448. xuex2=xuex2*(3./115.);
  449. uex2= prog 80. 70. 54. 41. 32.5 24. 18.5 12.5 10. 6.5 6.5 3. 3. 3.
  450. 1. 1. 1. 1.;
  451. uex2 = ((prog 18*183.) - uex2)*(1./183.);
  452. uex = uex1 et uex2;
  453. xuex = xuex1 et xuex2;
  454. evuex = evol manu xuex uex;
  455.  
  456. TAB1='TABLE' ;
  457. TAB1 . 1 = 'MOT' 'MARQ LOSA NOLI';
  458. *TAB1 . 2 = 'MOT' 'TIRR ';
  459. TAB1 . 'TITRE' = 'TABLE' ;
  460. TAB1 . 'TITRE' . 1 = 'MOT' 'Expérience CEAT/LEAT';
  461.  
  462. dess (evuex et evux)
  463. TITR (chai '(U-Ub)/(Ua-Ub) : Tadim=' Tadim ' CFL=' CFLI)
  464. TITX '(Y-Y0)/delta_w' TITY '(U-Ub)/(Ua-Ub)' 'LEGE' TAB1 ;
  465.  
  466. dess (evkex et evkx)
  467. TITR (chai 'k/(Ua-Ub)**2. : Tadim=' Tadim ' CFL=' CFLI)
  468. TITX '(Y-Y0)/delta_w' TITY 'k/(Ua-Ub)**2.' 'LEGE' TAB1 ;
  469.  
  470. dess evTx TITR (chai 'TETA : Tadim=' Tadim ' CFL=' CFLI)
  471. TITX '(Y-Y0)/delta_w' TITY 'TETA' ;
  472.  
  473. dess evksix TITR (chai 'Ksi : Tadim=' Tadim ' CFL=' CFLI)
  474. TITX '(Y-Y0)/delta_w' TITY 'Ksi' ;
  475.  
  476. ********** FIN Autosimilarité du profil transversal de vitesse ***********************
  477.  
  478. pe = rv.inco.'PRESSION' ;
  479. pn = elno pe $mf KPRES;
  480. trace pn mf cmf TITR ' Pression' ;
  481.  
  482.  
  483. trace kn mf cmf TITR 'KN';
  484. en = rv.inco.'EN';
  485. trace en mf cmf TITR 'EN';
  486. nut = cnu * kn * kn * (inve en) ;
  487. nutsnu= nut* (1./NU);
  488. trace nutsnu mf cmf TITR ' Nut / Nu';
  489.  
  490. ksi = rv.inco.'Ksi';
  491. tkti= rv.inco.'TKTI' ;
  492. To = rv.inco.'TKTE' ;
  493. trace tkti mf cmf TITR ' TKTI';
  494. trace To mf cmf TITR ' To ';
  495. trace ksi mf cmf TITR ' ksi ';
  496. Finsi;
  497. * Fin Graphiques
  498.  
  499. Si (NON COMPLET);
  500. xd= extr edltexp 'ABSC';
  501. yd= extr edltexp 'ORDO';
  502. xc= extr evdelta 'ABSC';
  503. yc= extr evdelta 'ORDO';
  504. list yc;
  505.  
  506. ycr=prog
  507. 5.41126E-03 5.41313E-03 5.41535E-03 5.41784E-03 5.42038E-03
  508. 5.42327E-03 5.42664E-03 5.43045E-03 5.43462E-03 5.43922E-03
  509. 5.44424E-03 5.44978E-03 5.45584E-03 5.46250E-03 5.46982E-03
  510. 5.47788E-03 5.48679E-03 5.49666E-03 5.50763E-03 5.51989E-03
  511. 5.53364E-03 5.54914E-03 5.56668E-03 5.58681E-03 5.60970E-03
  512. 5.63609E-03 5.66645E-03 5.70168E-03 5.74247E-03 5.79062E-03;
  513. ycr=ycr et (prog
  514. 5.84654E-03 5.91233E-03 5.98995E-03 6.08212E-03 6.18944E-03
  515. 6.31623E-03 6.46651E-03 6.64144E-03 6.84154E-03 7.06852E-03
  516. 7.32907E-03 7.61250E-03 7.90744E-03 8.19769E-03 8.48321E-03
  517. 8.68052E-03 8.83170E-03 8.95027E-03 9.34207E-03 9.83521E-03
  518. 1.03896E-02 1.10250E-02 1.18099E-02 1.26858E-02 1.36029E-02);
  519. ycr=ycr et (prog
  520. 1.44548E-02 1.53243E-02 1.61430E-02 1.68317E-02 1.77894E-02
  521. 1.89391E-02 1.98979E-02 2.07947E-02 2.15067E-02 2.25427E-02
  522. 2.38273E-02 2.52953E-02 2.62061E-02 2.72421E-02 3.12650E-02
  523. 3.40889E-02 3.30795E-02 );
  524.  
  525. xd= (prog (mini xc)) et xd et (prog (maxi xc));
  526. yd= (prog (mini yd)) et yd et (prog (maxi yd));
  527. yd=yd - (prog 20*0.0106);
  528.  
  529. *edltexpr= evol manu xd yd;
  530. edltexpr= evol manu xc ycr;
  531.  
  532.  
  533. *evd = edltexpr et evdelta ; dess evd;
  534. m=(INTG edltexpr 'ABSO') SOMM ;
  535. ddeltaw= SOMM (INTG (edltexpr - evdelta) 'ABSO')/m;
  536. mess ' Deltaw : m=' m ' ddeltaw=' ddeltaw;
  537. Si (ddeltaw > 1.e-3 );ERREUR 5;Finsi;
  538. FINSI;
  539.  
  540.  
  541. FIN;
  542.  
  543.  
  544.  
  545.  
  546.  
  547.  
  548.  
  549.  

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