Télécharger la plaquette Cast3M


fr en

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') ;
  535. ddeltaw= (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.  
  550.  
  551.  

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