Télécharger dy_devo4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dy_devo4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dy_devo4.dgibi: Jeux de donnees *
  5. * ------------------------------------ *
  6. * *
  7. * si GRAPH = N, les graphiques ne sont pas affiches
  8. * si GRAPH different de N, tous les graphiques
  9. * sont affiches
  10. *
  11. GRAPH = 'N' ; d_dyne = faux ;
  12. * test comparaison pasapas et dyn
  13. *
  14. SAUT PAGE ;
  15. *
  16. *-----------------------------------------------------*
  17. * DY_DEVO4 *
  18. *-----------------------------------------------------*
  19. * *
  20. * Tests de l'op{rateur DYNE option DE_VOGELAERE *
  21. * _____________________________________________ *
  22. * *
  23. * *
  24. * 1 poutre encastree, 10 elements finis, *
  25. * masse negligeable *
  26. * 1 masse concentree e son extremite *
  27. * Test PROFIL_PROFIL_EXTERIEUR avec des carres *
  28. * Comparaison avec la liaison POINT_PLAN *
  29. * NB_M: nombre de segments sur le cote du carre mobile*
  30. * NB_F: nombre de segments sur le cote du carre fixe *
  31. * *
  32. *-----------------------------------------------------*
  33.  
  34. OPTION ECHO 1 ;
  35. TEMPS ;
  36. *
  37. NB_M = 3 ;
  38. NB_F = 1 ;
  39. R_PPL = 1.E+05 ;
  40. R_PP = R_PPL ;
  41. EXPO = 0.5 ;
  42. *
  43. OPTION ELEM SEG2 DIME 3 ;
  44. OPTI EPSI LINEAIRE;
  45. *
  46. * Definition geometrie
  47. *
  48. P1 = 0. 0. 0. ;
  49. P2 = 0. 0. 0.9 ;
  50. LIG1 = P1 D 10 P2 ;
  51. *
  52. * Definition materiau
  53. *
  54. MODL1 = MODE LIG1 MECANIQUE ELASTIQUE POUT;
  55.  
  56. MA1 = MATE MODL1 YOUNG 2.1E11 NU 0.3 RHO 78. ;
  57. CA1 = 'CARA' MODL1 SECT 1.1E-6 INRY 8.33E-10
  58. INRZ 8.33E-10 TORS 5.0E-7 VECT (1. 0. 1. );
  59. *
  60. RIG1 = RIGI MODL1 (MA1 ET CA1) ;
  61. MAST = (MASSE MODL1 (MA1 ET CA1)) ET
  62. ( MASSE DEPL 2.0 P2 ) ;
  63. *
  64. * Conditions aux limites
  65. *
  66. ENC1 = BLOQUE DEPL ROTA P1 ;
  67. ENC2 = BLOQUE 'UY' LIG1 ;
  68. RIGT = RIG1 ET ENC1 ET ENC2 ;
  69. *
  70. * Calcul des modes
  71. *
  72. MOD1=VIBRATION PROCHE ( PROG 1.0 ) RIGT MAST ;
  73. *
  74. * Base modale
  75. *
  76. TB1 = TRADUIRE MOD1 ;
  77. TBAS = TABLE 'BASE_MODALE' ;
  78. TBAS.'MODES' = TB1 ;
  79.  
  80. ma1 = manu poi1 tb1 . 1 . point_repere ;
  81. mo1 = mode ma1 mecanique elastique modal ;
  82. ca1 = mate mo1 freq tb1 . 1 . frequence
  83. mass tb1 . 1 . masse_generalisee defo tb1 . 1 . deformee_modale
  84. amor 0. ;
  85.  
  86. *
  87. RIGI1 = RIGI TB1 ;
  88. *
  89. * Amortissement
  90. *
  91. AMOR1 = AMOR TBAS ( PROG 0.0 ) ;
  92. amor2 = amor mo1 ca1 ;
  93. *
  94. TA = TABLE 'AMORTISSEMENT' ;
  95. TA.'AMORTISSEMENT' = AMOR1 ;
  96. *
  97. * Deplacement initial
  98. *
  99. FE2 = FORCE FX 1. P2 ;
  100. XE2 = RESOU RIGT FE2 ;
  101. XMAX = MAXI XE2 AVEC (MOTS UX ) ;
  102. FE1 = FE2 * (0.139E-2 / XMAX );
  103. FGE1 = PJBA TBAS FE1 ;
  104. ALO = RESO RIGI1 FGE1 ;
  105. *
  106. TINIT = TABLE 'INITIAL' ;
  107. TINIT.'DEPLACEMENT' = ALO ;
  108. *
  109. *--- Definition des liaisons
  110. *
  111. TLIA = TABLE 'LIAISON' ;
  112. TLB = TABLE 'LIAISON_B' ;
  113. TLIA.'LIAISON_B' = TLB ;
  114. NORM1 = 0 0 1 ;
  115. *
  116. * Liaison PROFIL_PROFIL_EXTERIEUR
  117. *
  118. * profil mobile
  119. P1_1 = -0.5 0.0 0.9 ;
  120. P1_2 = 0.0 -0.5 0.9 ;
  121. P1_3 = 0.5 0.0 0.9 ;
  122. P1_4 = 0.0 0.5 0.9 ;
  123. MAIL_1 = P1_1 D NB_M P1_2 D NB_M P1_3 D NB_M P1_4 D
  124. NB_M P1_1 ;
  125. *
  126. * profil fixe
  127. P2_1 = 0.501 1.0 0.9 ;
  128. P2_2 = 0.501 -1.0 0.9 ;
  129. P2_3 = 0.601 -1.0 0.9 ;
  130. P2_4 = 0.601 1.0 0.9 ;
  131. MAIL_2 = P2_1 D NB_F P2_2 D NB_F P2_3 D NB_F P2_4 D
  132. NB_F P2_1 ;
  133. *
  134. SI ( NEG GRAPH 'N' ) ;
  135. TRAC ( 0 0 1000 ) ( MAIL_1 ET MAIL_2 ) ;
  136. FINSI ;
  137. *
  138. TL1 = TABLE 'LIAISON_ELEMENTAIRE' ;
  139. TL1.'TYPE_LIAISON' = MOT 'PROFIL_PROFIL_EXTERIEUR';
  140. TL1.'SUPPORT' = P2 ;
  141. TL1.'PROFIL_FIXE' = MAIL_2 ;
  142. TL1.'PROFIL_MOBILE' = MAIL_1 ;
  143. TL1.'NORMALE' = NORM1 ;
  144. TL1.'RAIDEUR' = R_PP ;
  145. TL1.'EXPOSANT_RAIDEUR' = EXPO ;
  146. *
  147. TLB . 1 = TL1 ;
  148. mappe1 = manu poi1 p2 ;
  149. moppe1 = modele mappe1 liaison PROFIL_PROFIL EXTERIEUR ;
  150. tuppe1 = table ; tuppe1 . moppe1 = vrai ;
  151. cappe1 = mate moppe1 'NORM' norm1 'RAID' r_pp 'ERAI' expo
  152. 'PFIX' mail_2 'PMOB' mail_1 'SORT' tuppe1 ;
  153.  
  154. *
  155. *--- Definition des resultats en sortie
  156. *
  157. TSORT = TABLE 'SORTIE' ;
  158. TSOR1 = TABLE 'SORTIE' ;
  159. TSORT . 'VARIABLE' = TSOR1 ;
  160. TSOR1 . 'DEPLACEMENT' = VRAI ;
  161. TSOR1 . 'VITESSE' = VRAI ;
  162. *
  163. TSOR2 = TABLE 'SORTIE' ;
  164. TSORT.'LIAISON_B' = TSOR2 ;
  165. TSOR2.TL1 = VRAI ;
  166. *
  167. * Temps
  168. *
  169. PDT = 0.1E-3 ;
  170. NP1 = 5000 ;
  171. NINS = 5 ;
  172. *
  173. ev1 = evol manu (prog 0. 1000.) (prog 1. 1.) ;
  174. cha1 = char force ev1 FGE1 ;
  175. *
  176. SI d_dyne ;
  177. tabdyn = table 'PASAPAS' ;
  178. tabdyn . modele = mo1 et moppe1 ;
  179. tabdyn . caracteristiques = ca1 et cappe1 ;
  180. tabdyn . sortie = table ;
  181. tabdyn . sortie . moppe1 = vrai ;
  182. tabdyn . sortie . vitesses = faux ;
  183. tabdyn . pas_de_temps = pdt ;
  184. tabdyn . nombre_pas = np1 ;
  185. tabdyn . pas_de_sortie = nins ;
  186. tabdyn . temps = table ;
  187. tabdyn . temps . 0 = 0. ;
  188. tabdyn . deplacements = table ;
  189. tabdyn . deplacements . 0 = ALO ;
  190. dyne de_vogelaere tabdyn ;
  191. SINON ;
  192. * pour pasapas
  193. * etape precontrainte
  194. t_prec = table ;
  195. t_prec . modele = mo1 ;
  196. t_prec . caracteristiques = ca1 ;
  197. t_prec . chargement = cha1 ;
  198. t_prec . temps_calcules = prog 0. 1000. ;
  199. t_prec . mes_sauvegardes = table ;
  200. t_prec . mes_sauvegardes . defto = vrai ;
  201. pasapas t_prec ;
  202.  
  203. *lacher
  204. tabdyn = table 'PASAPAS' ;
  205. tabdyn . modele = mo1 et moppe1 ;
  206. tabdyn . caracteristiques = ca1 et cappe1 ;
  207. *tabdyn . modele = mo1 ;
  208. *tabdyn . caracteristiques = ca1 ;
  209. tabdyn . temps = table ;
  210. tabdyn . temps . 0 = 0. ;
  211. tabdyn . deplacements = table ;
  212. tabdyn . deplacements . 0 = t_prec . deplacements . 1 ;
  213. tabdyn . contraintes = table ;
  214. tabdyn . contraintes . 0 = t_prec . contraintes . 1 ;
  215. tabdyn . deformations = table;
  216. tabdyn . deformations . 0 = t_prec . deformations . 1;
  217. tabdyn . dynamique = vrai ;
  218. tabdyn . temps_calcules = prog 0 pas (4*nins*pdt) (np1*pdt);
  219. pasapas tabdyn ;
  220.  
  221. FINSI ;
  222. *------ Op{rateur DYNE ------
  223. *
  224. TRESU1 = DYNE 'DE_VOGELAERE' TBAS TA TINIT TSORT TLIA
  225. NP1 PDT NINS ;
  226.  
  227. nsu1 = dime tresu1 . temps_de_sortie ;
  228. ovi1 = prog ;
  229. repeter bsu1 nsu1 ;
  230. ovi1 = ovi1 et
  231. (prog (extr tresu1 . &bsu1 . vitesse alfa tb1 . 1 . point_repere)) ;
  232. fin bsu1 ;
  233. evvi1 = evol manu 'temps' TRESU1 . 'TEMPS_DE_SORTIE'
  234. 'vitesse' ovi1 ;
  235. *
  236. * Trac{ des courbes
  237. *
  238. SI ( NEG GRAPH 'N' ) ;
  239. si d_dyne ;
  240. l_t2 = prog ;
  241. repeter bs (dime tabdyn . temps) ;
  242. si (ega &bs 1) ; iterer bs ; finsi ;
  243. l_t2 = l_t2 et (prog tabdyn . temps . (&bs - 1)) ;
  244. fin bs ;
  245. evfn1 = evol manu l_t2 tabdyn . force_de_choc ;
  246. evux1 = evol manu l_t2 tabdyn . ux ;
  247. TITRE '***** LIAISON PROFIL_PROFIL_EXTERIEUR *****';
  248. DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  249. 'FN' TRESU1 . TL1 . 'FORCE_DE_CHOC') et (evfn1 coul bleu))
  250. 'MIMA' ;
  251. *
  252. DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  253. 'UX' TRESU1 . TL1 . 'UX' ) et (evux1 coul bleu)) 'MIMA';
  254. *
  255. sinon ;
  256. *
  257. l_t20 = prog tabdyn . temps . 1 ;
  258. tt1 = 'EXTR' tabdyn . variables_internes . 1 'SLIA' 2 1 1 ;
  259. lutt1 = dime tt1 . force_de_choc ;
  260. l_fdc20 = extr tt1 . force_de_choc lutt1;
  261. l_ux20 = extr tt1 . ux lutt1 ;
  262. l_vi20 = prog
  263. (extr tabdyn . vitesses . 1 alfa tb1 . 1 . point_repere);
  264. indi1 = 1 ;
  265.  
  266. repeter bs ((dime tabdyn . temps) - 2) ;
  267. indi1 = indi1 + 1 ;
  268. ttu1 = 'EXTR' tabdyn . variables_internes . indi1 'SLIA' 2 1 1 ;
  269. l_ux20 = l_ux20 et ttu1 . ux ;
  270. l_fdc20 = l_fdc20 et ttu1 . force_de_choc ;
  271. l_vi20 = l_vi20 et
  272. (prog (extr tabdyn . vitesses . indi1 alfa tb1 . 1 . point_repere));
  273. l_t20 = l_t20 et (prog tabdyn . temps . indi1) ;
  274. fin bs ;
  275. evfdc20 = evol manu 'temps' l_t20 'force' l_fdc20 ;
  276. evux20 = evol manu 'temps' l_t20 'dep_pl' l_ux20 ;
  277. evit20 = evol manu 'temps' l_t20 'vitesse' l_vi20 ;
  278. evfdcr1 = evol manu 'temps' TRESU1 . 'TEMPS_DE_SORTIE'
  279. 'FN' TRESU1 . TL1 . 'FORCE_DE_CHOC' ;
  280. evuxr1 = EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  281. 'UX' TRESU1 . TL1 . 'UX' ;
  282. titre ' vitesse ' ;
  283. DESS (evvi1 et (evit20 coul bleu)) 'MIMA' ;
  284. titre 'force de choc' ;
  285. DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  286. 'FN' TRESU1 . TL1 . 'FORCE_DE_CHOC') et (evfdc20 coul bleu))
  287. 'MIMA' ;
  288. titre ' ux liaison ' ;
  289. DESS (( EVOL MANU 'TEMPS' TRESU1 . 'TEMPS_DE_SORTIE'
  290. 'UX' TRESU1 . TL1 . 'UX' ) et (evux20 coul bleu)) 'MIMA';
  291. *
  292. finsi ;
  293. FINSI ;
  294.  
  295.  
  296. * mess (extr tresu1 . temps_de_sortie 801) tabdyn . temps . 200 ;
  297. D_FC_1 = EXTRAIRE ( TRESU1 . TL1 . 'FORCE_DE_CHOC' )
  298. 801 ;
  299. D_DE_1 = EXTRAIRE ( TRESU1 . TL1 . 'UX' ) 801 ;
  300.  
  301. tt200 = 'EXTR' tabdyn . variables_internes . 200 'SLIA' 2 1 1 ;
  302. D_FC_2 =EXTRAIRE (tt200 . 'FORCE_DE_CHOC') 1;
  303. D_DE_2 =EXTRAIRE (tt200 . 'UX') 1;
  304.  
  305. *
  306. SAUT PAGE ;
  307. MESS '-----> liaison PROFIL_PROFIL_EXTERIEUR' ;
  308. MESS ' valeur de la force de choc :' D_FC_1 D_FC_2;
  309. MESS ' valeur du deplacement :' D_DE_1 D_DE_2;
  310. *
  311. R_FC = ABS ( ( D_FC_1 - D_FC_2 ) / D_FC_2 ) ;
  312. R_DE = ABS ( ( D_DE_1 - D_DE_2 ) / D_DE_2 ) ;
  313. R_DE0 = ABS((D_DE_1 - 1.29833E-03) / 1.29833E-03) ;
  314.  
  315. SI ( ( R_FC <EG (1.4497 * 1.2) ) ET ( R_DE <EG (0.15738 *1.2) )
  316. ET (R_DE0 < 1.e-2) ) ;
  317. ERRE 0 ;
  318. SINON ;
  319. ERRE 5 ;
  320. FINSI ;
  321.  
  322.  
  323. * Liaison POINT_PLAN
  324. *
  325. NORM1 = 1 0 0 ;
  326. TL1 = TABLE 'LIAISON_ELEMENTAIRE' ;
  327. TL1.'TYPE_LIAISON' = MOT 'POINT_PLAN' ;
  328. TL1.'SUPPORT' = P2 ;
  329. TL1.'NORMALE' = NORM1 ;
  330. TL1.'RAIDEUR' = R_PPL ;
  331. TL1.'JEU' = 0.001 ;
  332. *
  333. TLB . 1 = TL1 ;
  334. TSOR2.TL1 = VRAI ;
  335. mapp2 = manu poi1 p2 ;
  336. mopp2 = modele mapp2 liaison POINT_PLAN ;
  337. tuppe2 = table ; tuppe2 . mopp2 = vrai ;
  338. capp2 = mate mopp2 'NORM' norm1 'RAID' r_ppl 'JEU' 0.001 'SORT' tuppe2;
  339. *
  340. tabdyn2 = table 'PASAPAS' ;
  341. tabdyn2 . modele = mo1 et mopp2 ;
  342. tabdyn2 . caracteristiques = ca1 et capp2 ;
  343. tabdyn2 . temps = table ;
  344. tabdyn2 . temps . 0 = 0. ;
  345.  
  346. SI d_dyne ;
  347. tabdyn2 . deplacements = table ;
  348. tabdyn2 . deplacements . 0 = ALO ;
  349. tabdyn2 . sortie = table ;
  350. tabdyn2 . sortie . mopp2 = vrai ;
  351. tabdyn2 . sortie . vitesses = faux ;
  352. tabdyn2 . pas_de_temps = pdt ;
  353. tabdyn2 . nombre_pas = np1 ;
  354. tabdyn2 . pas_de_sortie = nins ;
  355. dyne de_vogelaere tabdyn2 ;
  356.  
  357. SINON ;
  358. tabdyn2 . deplacements = table ;
  359. tabdyn2 . deplacements . 0 = t_prec . deplacements . 1 ;
  360. tabdyn2 . contraintes = table ;
  361. tabdyn2 . contraintes . 0 = t_prec . contraintes . 1 ;
  362. tabdyn2 . deformations = table;
  363. tabdyn2 . deformations . 0 = t_prec . deformations . 1;
  364. tabdyn2 . dynamique = vrai ;
  365. tabdyn2 . temps_calcules = prog 0 pas (4*nins*pdt) (np1*pdt);
  366. *tabdyn2 . temps_calcules = prog 0 pas (pdt) .0005;
  367. *tabdyn2 . temps_sauves = prog 0 pas (nins*pdt) (np1*pdt) ;
  368. *tabdyn2 . temps_sauves = prog 0 pas (nins*pdt) .01 ;
  369. pasapas tabdyn2 ;
  370.  
  371. FINSI ;
  372.  
  373. *np1 = 5 ; nins = 1;
  374. TRESU2 = DYNE 'DE_VOGELAERE' TBAS TA TINIT TSORT TLIA
  375. NP1 PDT NINS ;
  376.  
  377. nsu2 = dime tresu2 . temps_de_sortie ;
  378. ovi2 = prog ;
  379. repeter bsu2 nsu2 ;
  380. ovi2 = ovi2 et
  381. (prog (extr tresu2 . &bsu2 . vitesse alfa tb1 . 1 . point_repere)) ;
  382. fin bsu2 ;
  383. evvi2 = evol manu 'temps' TRESU2 . 'TEMPS_DE_SORTIE'
  384. 'vitesse' ovi2 ;
  385. *
  386. * Trace des courbes
  387. *
  388. SI ( NEG GRAPH 'N' ) ;
  389. si d_dyne ;
  390. l_t22 = prog ;
  391. repeter bs (dime tabdyn2 . temps) ;
  392. si (&bs ega 1) ; iterer bs ; finsi ;
  393. l_t22 = l_t22 et (prog tabdyn2 . temps . (&bs - 1)) ;
  394. fin bs ;
  395. evfn2 = evol manu l_t22 tabdyn2 . force_de_choc;
  396. evux2 = evol manu l_t22 tabdyn2 . ux ;
  397. TITRE '***** LIAISON POINT_PLAN *****' ;
  398. DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  399. 'FN' TRESU2 . TL1 . 'FORCE_DE_CHOC') et (evfn2 coul
  400. vert)) 'MIMA' ;
  401. *
  402. DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  403. 'UX' TRESU2 . TL1 . 'UX' ) et (evux2 coul vert)) 'MIMA' ;
  404. sinon ;
  405. *
  406. l_t202 = prog tabdyn2 . temps . 1 ;
  407. tt2 = 'EXTR' tabdyn2 . variables_internes . 1 'SLIA' 2 1 1 ;
  408. lutt2 = dime tt2 . force_de_choc ;
  409. l_fdc202 = extr tt2 . force_de_choc lutt2;
  410. l_ux202 = extr tt2 . ux lutt2 ;
  411. l_vi202 = prog
  412. (extr tabdyn2 . vitesses . 1 alfa tb1 . 1 . point_repere);
  413. indi1 = 1 ;
  414.  
  415. repeter bs ((dime tabdyn2 . temps) - 2) ;
  416. indi1 = indi1 + 1 ;
  417. ttu2 = 'EXTR' tabdyn2 . variables_internes . indi1 'SLIA' 2 1 1 ;
  418. l_ux202 = l_ux202 et ttu2 . ux ;
  419. l_fdc202 = l_fdc202 et ttu2 . force_de_choc ;
  420. l_vi202 = l_vi202 et
  421. (prog (extr tabdyn2 . vitesses .indi1 alfa tb1 . 1 . point_repere)) ;
  422. l_t202 = l_t202 et (prog tabdyn2 . temps . indi1) ;
  423. fin bs ;
  424. evfdc202 = evol manu 'temps' l_t202 'force' l_fdc202 ;
  425. evux202 = evol manu 'temps' l_t202 'dep_pl' l_ux202 ;
  426. evit202 = evol manu 'temps' l_t202 'vitesse' l_vi202 ;
  427. evfdcr2 = evol manu 'temps' TRESU2 . 'TEMPS_DE_SORTIE'
  428. 'FN' TRESU2 . TL1 . 'FORCE_DE_CHOC' ;
  429. evuxr2 = EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  430. 'UX' TRESU2 . TL1 . 'UX' ;
  431. titre ' vitesse ' ;
  432. DESS (evvi2 et (evit202 coul bleu)) 'MIMA' ;
  433. titre 'force de choc' ;
  434. DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  435. 'FN' TRESU2 . TL1 . 'FORCE_DE_CHOC') et (evfdc202 coul bleu))
  436. 'MIMA' ;
  437. titre ' ux liaison ' ;
  438. DESS (( EVOL MANU 'TEMPS' TRESU2 . 'TEMPS_DE_SORTIE'
  439. 'UX' TRESU2 . TL1 . 'UX' ) et (evux202 coul bleu)) 'MIMA';
  440. *
  441. finsi ;
  442. FINSI ;
  443. *
  444. TEMPS ;
  445. *
  446. * Code de bon fonctionnement
  447. *
  448. D_FC_20 =EXTRAIRE ( TRESU2 . TL1 . 'FORCE_DE_CHOC' ) 801;
  449. D_DE_20 =EXTRAIRE ( TRESU2 . TL1 . 'UX' ) 801 ;
  450. tt2200 = 'EXTR' tabdyn2 . variables_internes . 200 'SLIA' 2 1 1 ;
  451. D_FC_21 =EXTRAIRE (tt2200 . 'FORCE_DE_CHOC') 1;
  452. D_DE_21 =EXTRAIRE (tt2200 . 'UX') 1;
  453.  
  454. SAUT 2 LIGNES ;
  455. MESS '-----> liaison POINT_PLAN ' ;
  456. MESS ' valeur de la force de choc :' D_FC_20 D_FC_21 ;
  457. MESS ' valeur du deplacement :' D_DE_20 D_DE_21 ;
  458. SAUT 3 LIGNES ;
  459. R_FC = ABS ( ( D_FC_20 - D_FC_21 ) / D_FC_20 ) ;
  460. R_DE = ABS ( ( D_DE_20 - D_DE_21 ) / D_DE_20 ) ;
  461. R_DE0 = ABS ((D_DE_20 - 1.29833E-03) / 1.29833E-03) ;
  462. SI ( ( R_FC &lt;EG (0.59178 *1.2)) ET ( R_DE &lt;EG (0.13598 * 1.2))
  463. ET (R_DE0 < 1.e-2) ) ;
  464. ERRE 0 ;
  465. SINON ;
  466. ERRE 5 ;
  467. FINSI ;
  468. *
  469.  
  470. FIN ;
  471.  
  472.  
  473.  
  474.  
  475.  
  476.  
  477.  
  478.  
  479.  
  480.  
  481.  
  482.  
  483.  
  484.  
  485.  
  486.  

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