Télécharger dy_dev12.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dy_dev12.dgibi
  2. ************************************************************************
  3. *
  4. * Cas-Test de la liaison palier de l'operateur DYNE
  5. *
  6. *
  7. * Masse ponctuelle sur 1 squeeze-film (NL)
  8. * Chargement statique = poids propre = -m*g --> Wcharg
  9. * + Chargement tournant = balourd = m*u*w^2 --> Wbal
  10. * avec ou sans centering spring de raideur k
  11. *
  12. *
  13. * ref : JI Inayat-Hussain et al, On the bifurcations of a rigid rotor
  14. * response in squeeze-film dampers, Journal of Fluids and
  15. * Structures 17 (2003) 433-459
  16. *
  17. *
  18. * Auteur : BP, 2015-05-05
  19. *
  20. *
  21. * Mots-clés : Vibrations non-lineaires, DYNE, squeeze-film,
  22. * dynamique des machines tournantes
  23. *
  24. ************************************************************************
  25. *
  26. OPTI 'DIME' 3 'ELEM' SEG2 ;
  27. OPTI EPSI LINEAIRE;
  28. COMPLET = FAUX;
  29. * COMPLET = VRAI;
  30. * GRAPH = FAUX;
  31. GRAPH = VRAI;
  32.  
  33. * FLCOUR=VRAI pour utiliser le modele de palier_court (sinon rodeli)
  34. * FLCOUR = FAUX;
  35. FLCOUR = VRAI;
  36.  
  37. * devogelaere ou diff centrees ?
  38. FLDEVO = vrai;
  39. FLDEVO = faux;
  40.  
  41. * donnees du squeeze :
  42. Lpalier = 0.05;
  43. Rarbre = 0.02;
  44. * mu = 0.02;
  45. rhofl = 900.;
  46. xjeu = 0.50E-3; xjeum1 = 1./xjeu;
  47. * adimensionne
  48. B = 0.015; valB = mot'0015';
  49. * B = 0.09; valB = mot'009';
  50. * B = 0.40; valB = mot'040';
  51.  
  52. * donnees du rotor (rigide) :
  53. M1 = 1.;
  54.  
  55. *amplitude force statique
  56. g = 10.;
  57. * poids = M1*g;
  58.  
  59. * poids propre adimensionne
  60. W = 0.05; valW = mot '005';
  61. * W = 0.1; valW = mot '010';
  62. * W = 0.5; valW = mot '050';
  63.  
  64. * on deduit la vitesse de rotation
  65. si (EGA W 0.);
  66. poids = 0.;
  67. OMEGHZ = 100. ;
  68. OMEGA = OMEGHZ * (2*pi);
  69. sinon;
  70. poids = M1*g;
  71. OMEGA = (g / (xjeu * W))**0.5;
  72. OMEGHZ = OMEGA/(2*pi);
  73. mess '\W = ' OMEGA 'rad/s = ' OMEGHZ 'Hz' ;
  74.  
  75. * on deduit la viscosite
  76. mu = B * M1 * OMEGA * (xjeu**3) / (Rarbre * (Lpalier**3));
  77. mess '\m = ' mu ' unite SI ' ;
  78.  
  79. * amplitude balourd
  80. * Fbalad = 0.05; valU = mot'005';
  81. Fbalad = 0.10; valU = mot'010';
  82. Fbalad = 0.20; valU = mot'020';
  83. Fbalad = 0.30; valU = mot'030';
  84. Fbalad = 0.355; valU = mot'0355';
  85. * Fbalad = 0.363; valU = mot'0363';
  86. * Fbalad = 0.3635; valU = mot'03635';
  87. * * Fbalad = 0.364; valU = mot'0364';
  88. * * Fbalad = 0.3641; valU = mot'03641';
  89. * * Fbalad = 0.370; valU = mot'0370';
  90. * Fbalad = 0.385; valU = mot'0385';
  91. * * Fbalad = 0.40; valU = mot'040';
  92. * Fbal = Fbalad * M1 * xjeu * (OMEGA**2);
  93. Fbal = Fbalad * M1 * xjeu * (OMEGA**2);
  94. mess 'Fbal='Fbal 'N';
  95.  
  96. * centering spring
  97. * S = 0.10; valS = mot '010';
  98. S = 0.30; valS = mot '030';
  99. * S = 0.50; valS = mot '050';
  100. k1 = M1 * ((S * OMEGA)**2);
  101. mess 'k='k1 'N/m';
  102.  
  103. Y0 = W / (S**2);
  104. mess 'static eccentricity ratio Y0='Y0;
  105.  
  106.  
  107. * valeurs adimensionnées du squeeze :
  108. cha1 = chai 'Squeeze-film : B=' FORMAT '(F5.3)' B
  109. ' - W=' FORMAT '(F5.3)' W
  110. ' - S=' FORMAT '(F5.3)' S
  111. ' - U=' FORMAT '(F5.3)' Fbalad ;
  112. mess cha1;
  113. TITR cha1;
  114.  
  115. * prefixe
  116. si FLCOUR; mopal = mot 'court';
  117. sinon; mopal = mot 'rodeli';
  118. finsi;
  119. si FLDEVO; moalgo = mot 'DV';
  120. sinon; moalgo = mot 'DC';
  121. finsi;
  122.  
  123. prefix = chai 'dy_dev12-'mopal'-'moalgo
  124. '-B'valB'-W'valW'-S'valS'-U'valU;
  125. iCI = 0; prefix = chai prefix '-CI0';
  126. * iCI = 1; prefix = chai prefix '-CI1';
  127. * iCI = 2; prefix = chai prefix '-CI2';
  128. * iCI = 3; prefix = chai prefix '-CI3';
  129. mess prefix ;
  130.  
  131. ficxdr = chai prefix'.xdr';
  132. ficps = chai prefix'.ps';
  133. OPTI 'TRAC' 'PSC' 'POTR' 'HELVETICA_16' 'EPTR' 6 'FTRA' ficps;
  134.  
  135.  
  136. ************************************************
  137. * Nombre de pas de temps, pas de temps
  138. ************************************************
  139.  
  140. si (COMPLET);
  141. si (((Fbalad > 0.36 ) et (Fbalad <eg 0.37))
  142. ou (Fbalad ega 0.385)); Nperiod = 1000 ;
  143. sinon; Nperiod = 400 ;
  144. finsi;
  145. Nperiod = 400 ;
  146. sinon; Nperiod = 80 ;
  147. finsi;
  148. * nombre de points par periode, nombre de period
  149. nepp1 = 8;
  150. npp1 = 2**nepp1;
  151. NPas1 = npp1 * Nperiod;
  152. dt1 = 1. / (OMEGHZ * npp1);
  153. Tfinal = NPas1 * dt1;
  154. * Sortie tous les NSor1 pas de temps
  155. * NSor1 = 2**3 ;
  156. NSor1 = 1 ;
  157. NbMod1 = 2;
  158.  
  159.  
  160. ************************************************
  161. * Creation des MODES
  162. ************************************************
  163.  
  164. P1 = 0. 0. 0. ;
  165. Mail1 = MANU 'POI1' P1;
  166. *
  167. moprim = mots 'UY' 'UZ';
  168. MMass1= MANU 'RIGIDITE' 'TYPE' 'MASSE' Mail1 moprim (prog M1 0 0 M1);
  169. MRigi1= MANU 'RIGIDITE' 'TYPE' 'RIGIDITE' Mail1 moprim (prog k1 0 0 k1);
  170. *
  171. chpuy = manu 'CHPO' Mail1 2 'UY' 1. 'UZ' 0. 'NATURE' 'DIFFUS';
  172. chpuz = manu 'CHPO' Mail1 2 'UY' 0. 'UZ' 1. 'NATURE' 'DIFFUS';
  173. lchp = SUIT chpuy chpuz;
  174.  
  175. TBasR1 = VIBR 'PROC' (prog -0.1) (lect NbMod1)
  176. MRigi1 MMass1 'TBAS' lchp;
  177.  
  178. mess 'mode 1 : f=' TBasR1 . 'MODES' . 1 . 'FREQUENCE'
  179. 'M=' TBasR1 . 'MODES' . 1 . 'MASSE_GENERALISEE';
  180. mess 'mode 2 : f=' TBasR1 . 'MODES' . 2 . 'FREQUENCE'
  181. 'M*=' TBasR1 . 'MODES' . 2 . 'MASSE_GENERALISEE';
  182.  
  183. *
  184. ************************************************
  185. * Creation de la table des CONDITIONS INITIALES
  186. ************************************************
  187.  
  188. TInit1 = TABLE 'INITIAL' ;
  189. *
  190. si (ega iCI 0); UY0 = 0. ; UZ0 = 0. ; finsi;
  191. * * 1er jeu de CI :
  192. si (ega iCI 1); UY0 = 0.1*xjeu ; UZ0 = -0.1*xjeu ; finsi;
  193. * 2eme jeu de CI :
  194. si (ega iCI 2); UY0 = 0.0*xjeu ; UZ0 = 0.5*xjeu ; finsi;
  195. * 3eme jeu de CI :
  196. si (ega iCI 3); UY0 = 0.0*xjeu ; UZ0 = -0.8*xjeu ; finsi;
  197. * si (ega iCI 3); UY0 = 0.1*xjeu ; UZ0 = -0.95*xjeu ; finsi;
  198.  
  199. Dep0 = MANU 'CHPO' Mail1 2 'UY' UY0 'UZ' UZ0 NATURE DIFFUS ;
  200. *
  201. REPE bb1 NbMod1 ;
  202. i1 = &bb1 ;
  203. Def1 = TBasR1 . 'MODES' . i1 . 'DEFORMEE_MODALE' ;
  204. M_G1 = TBasR1 . 'MODES' . i1 . 'MASSE_GENERALISEE' ;
  205. P_R1 = TBasR1 . 'MODES' . i1 . 'POINT_REPERE' ;
  206. qi1 = (YTMX Def1 MMass1 Dep0) / M_G1 ;
  207. f_r1 = TBasR1 . 'MODES' . i1 . 'FREQUENCE';
  208. mess 'Mode' i1 'f=' f_r1 ' M*=' qi1;
  209. SI (EGA i1 1) ;
  210. TInit1.'DEPLACEMENT' = MANU 'CHPO' P_R1 1 'ALFA' qi1 NATURE DIFFUS ;
  211. TInit1.'VITESSE' = MANU 'CHPO' P_R1 1 'ALFA' 0. NATURE DIFFUS ;
  212. SINON ;
  213. TInit1.'DEPLACEMENT' = TInit1.'DEPLACEMENT' ET
  214. (MANU 'CHPO' P_R1 1 'ALFA' qi1 NATURE DIFFUS) ;
  215. TInit1.'VITESSE' = TInit1.'VITESSE' ET
  216. (MANU 'CHPO' P_R1 1 'ALFA' 0. NATURE DIFFUS) ;
  217. FINSI ;
  218. FIN bb1 ;
  219.  
  220. *
  221. **********************************************
  222. * Creation de la table des LIAISONS (paliers)
  223. **********************************************
  224. *
  225. TLiai1 = TABLE 'LIAISON' ;
  226. TLiai1 . LIAISON_B = TABLE 'LIAISON_B' ;
  227. PAL1 = TABLE 'LIAISON_ELEMENTAIRE' ;
  228. TLiai1 . LIAISON_B . 1 = PAL1 ;
  229. PAL1 . TYPE_LIAISON = 'PALIER_FLUIDE' ;
  230. PAL1 . POINT_SUPPORT = P1 ;
  231. PAL1 . LONGUEUR_PALIER = Lpalier ;
  232. PAL1 . RAYON_ARBRE = Rarbre ;
  233. PAL1 . VISCOSITE_FLUIDE = mu ;
  234. PAL1 . RHO_FLUIDE = rhofl ;
  235. PAL1 . PRESSION_ADMISSION = 0. ;
  236. PAL1 . VITESSE_ROTATION = 0.;
  237. si (FLCOUR);
  238. PAL1 . MODELE_PALIER = 'PALIER_COURT' ;
  239. PAL1 . JEU_USINAGE = xjeu ;
  240. sinon;
  241. PAL1 . MODELE_PALIER = 'RODELI' ;
  242. PAL1 . AFFICHAGE = 0 ;
  243. PAL1 . GEOMETRIE_PALIER = TABLE ;
  244. PAL1 . GEOMETRIE_PALIER . NOMBRE_LOBES = 1 ;
  245. PAL1 . GEOMETRIE_PALIER . CRITERE_ARRET = 1e-6 ;
  246. PAL1 . GEOMETRIE_PALIER . 1 = TABLE ;
  247. PAL1 . GEOMETRIE_PALIER . 1 . ASYMETRIE = 0. ;
  248. PAL1 . GEOMETRIE_PALIER . 1 . JEU_USINAGE = xjeu ;
  249. PAL1 . GEOMETRIE_PALIER . 1 . PRECHARGE = 0. ;
  250. PAL1 . GEOMETRIE_PALIER . 1 . ANGLE_DEBUT = 0. ;
  251. PAL1 . GEOMETRIE_PALIER . 1 . AMPL_ANGLE = 2 * pi ;
  252. PAL1 . GEOMETRIE_PALIER . 1 . COEF_SUR = 1.715 ;
  253. PAL1 . GEOMETRIE_PALIER . 1 . NB_MAILLES = 120 ;
  254. * PAL1 . GEOMETRIE_PALIER . 1 . NB_MAILLES = 360 ;
  255. finsi;
  256.  
  257.  
  258.  
  259. ***********************************************
  260. * Creation de la table des VARIABLES DE SORTIE
  261. ***********************************************
  262. TSort1 = TABLE 'SORTIE' ;
  263. TSort1 . 'VARIABLE' = TABLE 'VARIABLE' ;
  264. TSort1 . 'LIAISON_B' = TABLE 'LIAISON_B' ;
  265. TSort1 . 'LIAISON_B' . PAL1 = VRAI ;
  266.  
  267.  
  268. *************************************************************
  269. * Creation de la table du CHARGEMENT
  270. *************************************************************
  271.  
  272. LTps1 = PROG 0.0 'PAS' dt1 'NPAS' (NPas1+1) ;
  273.  
  274. *** evol temporelle ***
  275.  
  276. *gravite
  277. LFW1 = PROG (NPas1+2)*poids;
  278. evW1 = evol ROSE MANU LTps1 LFW1;
  279. Chp_FZ1 = MANU CHPO 2 Mail1 'FY' 0. 'FZ' -1. ;
  280.  
  281.  
  282. *balourd
  283. LFdY1 = PROG 'SINU' OMEGHZ 'AMPL' Fbal LTps1 ;
  284. LFdZ1 = PROG 'SINU' OMEGHZ 'PHAS' 90. 'AMPL' Fbal LTps1 ;
  285. Ev_FdY1 = EVOL 'BLEU' MANU 'LEGE' 'FY' 't(s)' LTps1 'FY' LFdY1 ;
  286. Ev_FdZ1 = EVOL 'OCEA' MANU 'LEGE' 'FZ' 't(s)' LTps1 'FZ' LFdZ1 ;
  287.  
  288. dess ((-1.*evW1) et Ev_FdY1 et (-1.*Ev_FdZ1))
  289. 'POSX' 'CENT' 'TITX' 't(s)' 'XBOR' 0. (1./OMEGHZ)
  290. 'POSY' 'CENT' 'TITY' 'F(N)' 'LEGE' 'NO';
  291. dess ((-1.*evW1) et Ev_FdY1 et (-1.*Ev_FdZ1))
  292. 'POSX' 'CENT' 'TITX' 't(s)' 'XBOR' 0. Tfinal
  293. 'POSY' 'CENT' 'TITY' 'F(N)' 'LEGE' 'NO';
  294.  
  295. *** chpoints ***
  296. Chp_FdY1 = MANU CHPO 2 Mail1 'FY' 1. 'FZ' 0. ;
  297. Chp_FdZ1 = MANU CHPO 2 Mail1 'FY' 0. 'FZ' -1. ;
  298.  
  299. *** chargement ***
  300. Chatot = (CHAR Chp_FdY1 Ev_FdY1)
  301. ET (CHAR Chp_FdZ1 Ev_FdZ1)
  302. ET (CHAR Chp_FZ1 evW1) ;
  303.  
  304. TChar1 = TABLE 'CHARGEMENT' ;
  305. TChar1 . 'BASE_A' = PJBA TBasR1 Chatot ;
  306.  
  307.  
  308. ****************************************************************
  309. * Cacul temporel 1 : Appel a DYNE
  310. ****************************************************************
  311. si FLDEVO; algodyne = mot 'DE_VOGELAERE';
  312. sinon; algodyne = mot 'DIFFERENCES_CENTREES';
  313. finsi;
  314.  
  315. TResu1 = DYNE algodyne TBasR1 TChar1 TInit1 TSort1
  316. TLiai1 dt1 NPas1 NSor1 ;
  317.  
  318. uy1 = TResu1 . PAL1 . 'UY_ARBRE';
  319. uz1 = TResu1 . PAL1 . 'UZ_ARBRE';
  320. vy1 = TResu1 . PAL1 . 'VY_ARBRE';
  321. vz1 = TResu1 . PAL1 . 'VZ_ARBRE';
  322. fy1 = TResu1 . PAL1 . 'FY_PALIER';
  323. fz1 = TResu1 . PAL1 . 'FZ_PALIER';
  324. tprog = TResu1 . 'TEMPS_DE_SORTIE' ;
  325. nt = dime tprog;
  326.  
  327.  
  328. ********************************
  329. * graphiques
  330. ********************************
  331. si (GRAPH);
  332.  
  333. evuy1 = EVOL BLEU 'MANU' 'LEGE' 'UY/e' 't' tprog 'UY/e' (uy1*xjeum1);
  334. evuz1 = EVOL TURQ 'MANU' 'LEGE' 'UZ/e' 't' tprog 'UZ/e' (uz1*xjeum1);
  335. evvy1 = EVOL BLEU 'MANU' 'LEGE' 'VY' 't' tprog 'VY' (vy1);
  336. evvz1 = EVOL TURQ 'MANU' 'LEGE' 'VZ' 't' tprog 'VZ' (vz1);
  337. evfy1 = EVOL BLEU 'MANU' 't' tprog 'FY' fy1;
  338. evfz1 = EVOL TURQ 'MANU' 't' tprog 'FZ' fz1;
  339. evyz1 = EVOL BLEU 'MANU' 'UY/e' (uy1*xjeum1) 'UZ/e' (uz1*xjeum1);
  340. t1 = prog -180. PAS 1. 180.;
  341. y1 = sin t1; z1 = -1.* (cos t1);
  342. evjeu = evol 'MANU' y1 z1;
  343. Tjeu = tabl; Tjeu . 2 = mot 'TIRR';
  344.  
  345. si (non complet);
  346. dess (evuy1 et evuz1)
  347. 'XBOR' 0. Tfinal TITX 'U/e' POSX CENT POSY CENT LEGE NO;
  348. dess (evvy1 et evvz1)
  349. 'XBOR' 0. Tfinal POSX CENT POSY CENT LEGE NO;
  350. dess (evfy1 et evfz1)
  351. 'XBOR' 0. Tfinal POSX CENT POSY CENT LEGE NO;
  352. finsi;
  353. dess (evyz1 et evjeu)
  354. POSX CENT POSY CENT 'CARR' Tjeu;
  355.  
  356. * derniere periode
  357. si complet; neplast = 7 ;
  358. sinon; neplast = 4 ;
  359. finsi;
  360. nplast = 2**neplast ;
  361. nlast = (Nperiod - nplast) * npp1 / NSor1;
  362. llast = lect nlast pas 1 (NPas1 / NSor1);
  363. tprogf = extr tprog llast;
  364. tprogf = tprogf - (extr tprogf 1);
  365. uy1f = extr uy1 llast;
  366. uz1f = extr uz1 llast;
  367. evuy1f = EVOL BLEU 'MANU' 'LEGE' 'UY/e'
  368. 't' tprogf 'UY/e' (uy1f*xjeum1);
  369. evuz1f = EVOL TURQ 'MANU' 'LEGE' 'UZ/e'
  370. 't' tprogf 'UZ/e' (uz1f*xjeum1);
  371. evyz1f = EVOL BLEU 'MANU' 'UY/e' (uy1f*xjeum1) 'UZ/e' (uz1f*xjeum1);
  372. dess (evyz1f et evjeu)
  373. POSX CENT POSY CENT 'CARR' Tjeu;
  374. dess (evyz1f)
  375. POSX CENT POSY CENT 'CARR' Tjeu;
  376.  
  377. * DSP
  378. ndsp = nepp1 + neplast ;
  379. dspy1 = DSPR ndsp evuy1f 'BLEU';
  380. dspz1 = DSPR ndsp evuz1f 'TURQ';
  381. * adimensionnement des frequences (+ pratique)
  382. dspy1 = EVOL 'BLEU' 'MANU' 'LEGE' 'UY/e'
  383. '\w/\W' ((EXTR dspy1 'ABSC')/OMEGHZ)
  384. 'UY/e' (EXTR dspy1 'ORDO');
  385. dspz1 = EVOL 'TURQ' 'MANU' 'LEGE' 'UZ/e'
  386. '\w/\W' ((EXTR dspz1 'ABSC')/OMEGHZ)
  387. 'UZ/e' (EXTR dspz1 'ORDO');
  388. si (non complet);
  389. dess (dspy1 et dspz1)
  390. POSX CENT POSY CENT LOGY YBOR 1.E-15 1. LEGE NE ;
  391. finsi;
  392. dess (dspy1 et dspz1)
  393. POSX CENT POSY CENT XBOR 0. 20. LOGY YBOR 1.E-15 1. LEGE NE ;
  394. dess (dspy1 et dspz1)
  395. POSX CENT POSY CENT XBOR 0. 4. XGRA 0.5 LOGY YBOR 1.E-15 1. LEGE NE;
  396. ymax = maxi (EXTR dspy1 'ORDO');
  397. * nymax = enti ((log ymax) / (log 10.)) + 1;
  398. nymax = enti 'SUPERIEUR' ((log ymax) / (log 10.));
  399. nymin = nymax - 6;
  400. ymin = 10.**nymin;
  401. ymax = 10.**nymax;
  402. dess (dspy1 et dspz1)
  403. POSX CENT POSY CENT XBOR 0. 4. XGRA 0.5 LOGY YBOR ymin ymax LEGE NE;
  404. dess (dspy1 et dspz1)
  405. POSX CENT POSY CENT XBOR 0. 2. XGRA 0.25 LOGY YBOR ymin ymax LEGE NE;
  406.  
  407. finsi;
  408.  
  409.  
  410. * *excentricite et angle de calage
  411. * prex1 = ((uy1**2) + (uz1**2))**0.5;
  412. * mas00 = (masq (abs uy1) 'EGINFE' 1.E-30 )
  413. * * (masq (abs uz1) 'EGINFE' 1.E-30 );
  414. * uy100 = (((prog nt*1.) - mas00) * uy1) + mas00;
  415. * prang1 = ATG uy100 (-1.*uz1);
  416. * mess 'excentricite ='; list prex1;
  417. * mess 'angle de calage ='; list prang1;
  418. * tfin = extr tprog nt;
  419. * an1fin = extr prang1 nt;
  420. * ex1fin = extr prex1 nt;
  421. * mess 'excentricite relative (t=' tfin ')=' (ex1fin*xjeum1);
  422. *
  423.  
  424.  
  425. ****************************************************************
  426. * SAUVEGARDE DES COURBES
  427. ****************************************************************
  428. si (FAUX);
  429.  
  430. TEVOL1 = INDEX '*EVOLUTIO';
  431. list TEVOL1;
  432. opti sauv ficxdr;
  433. sauv TEVOL1;
  434.  
  435. finsi ;
  436.  
  437.  
  438. ****************************************************************
  439. * Tests de Non Regression
  440. ****************************************************************
  441.  
  442. * on verifie que l'on obtient un regime 4T-Periodique
  443.  
  444. * recup des 4 dernieres periodes
  445. Tuy = tabl;
  446. Tuz = tabl;
  447. nfin = NPas1 / NSor1;
  448. repe b4 5;
  449. nlast = (Nperiod - &b4) * npp1 / NSor1;
  450. llast = lect nlast pas 1 nfin;
  451. Tuy . &b4 = extr uy1 llast;
  452. Tuz . &b4 = extr uz1 llast;
  453. mess 'nlast - nfin' nlast nfin (dime Tuy . &b4);
  454. nfin = nlast;
  455. fin b4;
  456.  
  457. * ecart uy(t)-uy(t-1*T)
  458. yecart = prog;
  459. zecart = prog;
  460. repe b4 4 ;
  461. yecart = yecart et (MAXI (Tuy . 1 - Tuy . (&b4 + 1)) 'ABS');
  462. zecart = zecart et (MAXI (Tuz . 1 - Tuz . (&b4 + 1)) 'ABS');
  463. fin b4;
  464. yecart = yecart*xjeum1;
  465. zecart = zecart*xjeum1;
  466. * resultats obtenus par BP en 2016-07-07 (avec COMPLET=vrai):
  467. * 0.84326 0.45235 0.71956 5.85832E-11
  468. * 1.3414 0.24801 1.2888 5.62050E-11
  469. * resultats obtenus par BP en 2016-07-07 (avec COMPLET=faux):
  470. * 0.84326 0.45235 0.71956 3.73157E-08
  471. * 1.3414 0.24801 1.2888 2.20570E-08
  472. list yecart;
  473. list zecart;
  474. opti echo 0;
  475. mess '-------------------------------------------------------------';
  476. repe b4 4 ;
  477. mess ' Max_t{ y(t) - y(t-'&b4'*T)} = ' (extr yecart &b4);
  478. fin b4;
  479. mess '-------------------------------------------------------------';
  480. opti echo 1;
  481.  
  482. SI ( ((extr yecart 1) < 0.1) ou ((extr yecart 2) < 0.1)
  483. ou ((extr yecart 3) < 0.1) ou ((extr yecart 4) > 1.E-6)
  484. ou ((extr zecart 1) < 0.1) ou ((extr zecart 2) < 0.1)
  485. ou ((extr zecart 3) < 0.1) ou ((extr zecart 4) > 1.E-6) );
  486. ERRE 5;
  487. SINON;
  488. ERRE 0;
  489. FINSI;
  490.  
  491.  
  492. fin ;
  493.  
  494.  
  495.  
  496.  
  497.  

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