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

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