Télécharger dyna7.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna7.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna7.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test dyna7.dgibi: jeux de données *
  10. * --------------------------------- *
  11. *******************************************************
  12.  
  13. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  14. * SINON SI GRAPH DIFFERENT DE N TOUS
  15. * LES GRAPHIQUES SONT AFFICHES
  16.  
  17. GRAPH = 'N' ;
  18.  
  19. SAUT PAGE;
  20. SI (NEG GRAPH 'N') ;
  21. OPTI ECHO 1 ;
  22. OPTI TRAC X ;
  23. SINO ;
  24. OPTI ECHO 0 ;
  25. FINSI ;
  26.  
  27. SAUT PAGE;
  28. *******************************************************
  29. * DYNA7
  30. * CAS TEST OPERATEUR VIBRATION
  31. *
  32. *
  33. * Cet exemple permet de tester les trois options de
  34. * l'opérateur VIBR
  35. *
  36. * On calcule les modes de ballotement d'une lame mince
  37. * de fluide dans un espace annulaire en mode de
  38. * Fourier 1 et 3.
  39. *
  40. * On suppose la cavité indéformable. L'effet de la
  41. * gravité est pris en compte.
  42. *
  43. *
  44. *******************************************************
  45. TITR 'MODES DE BALLOTEMENT D UNE LAME MINCE DE FLUIDE';
  46. OPTI DIME 2 ELEM QUA4;
  47. OPTI EPSI LINEAIRE;
  48. TEMPS ;
  49.  
  50. *--------------------------- MAILLAGE -----------------
  51. *----------------- ON NE MAILLE QUE LE FLUIDE ---------
  52. *
  53. A1=0.995 0.;A2=1.005 0.;A3=1.005 2.;A4=0.995 2.;
  54. A5=1.005 1.2;A6=0.995 1.2;
  55. L1=D 1 A1 A2;
  56. L2=A2 D 24 A5 D 16 A3;
  57. L3=D 1 A3 A4;
  58. L4=A4 D 16 A6 D 24 A1;
  59. FLUID=DALL L1 L2 L3 L4;
  60. *
  61. SI (NEG GRAPH 'N');
  62. TRAC 'QUAL' FLUID ;
  63. FINSI;
  64. *
  65. OPTI ELEM SEG2;
  66. SURFA=D 1 A3 A4;
  67. *
  68. SI (NEG GRAPH 'N');
  69. TRAC 'QUAL' SURFA ;
  70. FINSI;
  71. *
  72. MAIL=FLUID ET SURFA;
  73. *
  74. SI (NEG GRAPH 'N');
  75. TRAC 'QUAL' MAIL;
  76. FINSI;
  77. *
  78. *---------- DEFINITION DE QUELQUES CONSTANTES ---------
  79. *
  80. ROREF=1.;CEREF=1.;LCARA=1.;
  81. RO0=1000.;CS=1200.;GRA=9.81;
  82. *
  83. *-- MODELE ET OBJET AFFECTE FLUIDE ET SURFACE LIBRE ---
  84. *-------- POUR L'HARMONIQUE DE FOURIER 1 --------------
  85. *
  86. OPTI MODE FOUR 1;
  87. MODLIQ1=MODE FLUID LIQUIDE LQU4;
  88. MODLIQ2=MODE SURFA LIQUIDE LSU2;
  89. MODLIQ=MODLIQ1 ET MODLIQ2;
  90. *
  91. *--------------- MATERIAU LIQUIDE ---------------------
  92. *
  93. MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF
  94. LCAR LCARA G GRA;
  95. *
  96. *--------------- RIGIDITE ET MASSE --------------------
  97. *
  98. RAI=RIGI MODLIQ MAT;
  99. MAS=MASS MODLIQ MAT;
  100. *
  101. *______________________________________________________
  102. * ON CALCULE LE PREMIER MODE AVEC L'OPTION INTERVALLE
  103. * DE VIBR. ON RENORMALISE LE CHAMP DE DEPLACEMENT
  104. * (ET DE PRESSION) DE MANIERE A CE QUE LA PRESSION SOIT
  105. * EGALE A 1 A L'ALTITUDE 0 ET ON EXTRAIT LES PRESSIONS
  106. * AUX POINTS A2, A5 ET A3
  107. *______________________________________________________
  108. *
  109. MOD1=VIBR INTE 0. 1. RAI MAS;
  110. F1 =MOD1 . MODES . 1 . FREQUENCE;
  111. PRE=MOD1 . MODES . 1 . DEFORMEE_MODALE;
  112. PBAS=EXTR PRE P A2;
  113. PP=PRE/PBAS;
  114. X1=EXTR PP P A2;Y1=EXTR PP P A5;Z1=EXTR PP P A3;
  115. *
  116. *______________________________________________________
  117. * ON CALCULE LE PREMIER MODE AVEC L'OPTION PROCHE DE VIBR
  118. *______________________________________________________
  119. *
  120. MOD2=VIBR PROC (PROG 0.5) RAI MAS;
  121. F2= MOD2 . MODES . 1 . FREQUENCE;
  122. PRE= MOD2 . MODES . 1 . DEFORMEE_MODALE;
  123. PBAS=EXTR PRE P A2;
  124. PP=PRE/PBAS;
  125. X2=EXTR PP P A2;Y2=EXTR PP P A5;Z2=EXTR PP P A3;
  126. *
  127. *______________________________________________________
  128. * ON CALCULE LE PREMIER MODE AVEC L'OPTION SIMULTANE DE VIBR
  129. *______________________________________________________
  130. *
  131. MOD3=VIBR SIMU 0.5 1 RAI MAS;
  132. F3= MOD3 . MODES . 1 . FREQUENCE;
  133. PRE= MOD3 . MODES . 1 . DEFORMEE_MODALE;
  134. PBAS=EXTR PRE P A2;
  135. PP=PRE/PBAS;
  136. X3=EXTR PP P A2;Y3=EXTR PP P A5;Z3=EXTR PP P A3;
  137. *______________________________________________________
  138. * ON CALCULE LE PREMIER MODE AVEC L'OPTION IRAM DE VIBR
  139. *______________________________________________________
  140. *
  141. MOD4=VIBR SIMU 0.5 1 RAI MAS;
  142. F4= MOD4 . MODES . 1 . FREQUENCE;
  143. PRE= MOD4 . MODES . 1 . DEFORMEE_MODALE;
  144. PBAS=EXTR PRE P A2;
  145. PP=PRE/PBAS;
  146. X4=EXTR PP P A2;Y4=EXTR PP P A5;Z4=EXTR PP P A3;
  147. *
  148. *______________________________________________________
  149. * ON RECOMMENCE LES OPERATIONS CETTE FOIS
  150. * AVEC L'HARMONIQUE 3
  151. *______________________________________________________
  152. *
  153. OPTI MODE FOUR 3;
  154. MODLIQ1=MODE FLUID LIQUIDE LQU4;
  155. MODLIQ2=MODE SURFA LIQUIDE LSU2;
  156. *
  157. MODLIQ=MODLIQ1 Et MODLIQ2;
  158. *
  159. MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF
  160. LCAR LCARA G GRA;
  161. *
  162. RAI=RIGI MODLIQ MAT;
  163. MAS=MASS MODLIQ MAT;
  164. *
  165. MOD1=VIBR INTE 0. 1. RAI MAS;
  166. *------------------------------------------------------
  167. G1= MOD1 . MODES . 1 . FREQUENCE;
  168. PRE= MOD1 . MODES . 1 . DEFORMEE_MODALE;
  169. PBAS=EXTR PRE P A2;
  170. PP=PRE/PBAS;
  171. T1=EXTR PP P A2;U1=EXTR PP P A5;V1=EXTR PP P A3;
  172. *------------------------------------------------------
  173. MOD2=VIBR PROC (PROG 0.9) RAI MAS;
  174. *
  175. G2=MOD2 . MODES . 1 . FREQUENCE;
  176. PRE= MOD2 . MODES . 1 . DEFORMEE_MODALE;
  177. PBAS=EXTR PRE P A2;
  178. PP=PRE/PBAS;
  179. T2=EXTR PP P A2;U2=EXTR PP P A5;V2=EXTR PP P A3;
  180. *------------------------------------------------------
  181. MOD3=VIBR SIMU 0.9 1 RAI MAS;
  182. G3=MOD3 . MODES . 1 . FREQUENCE;
  183. PRE= MOD3 . MODES . 1 . DEFORMEE_MODALE;
  184. PBAS=EXTR PRE P A2;
  185. PP=PRE/PBAS;
  186. T3=EXTR PP P A2;U3=EXTR PP P A5;V3=EXTR PP P A3;
  187. *------------------------------------------------------
  188. MOD4=VIBR SIMU 0.9 1 RAI MAS;
  189. G4=MOD4 . MODES . 1 . FREQUENCE;
  190. PRE= MOD4 . MODES . 1 . DEFORMEE_MODALE;
  191. PBAS=EXTR PRE P A2;
  192. PP=PRE/PBAS;
  193. T4=EXTR PP P A2;U4=EXTR PP P A5;V4=EXTR PP P A3;
  194. *------------------------------------------------------
  195. SAUT PAGE; SAUT 2 LIGNE ;
  196. MESS '*';
  197. MESS '******** ON VERIFIE LES FREQUENCES';
  198. MESS '*';
  199. MESS ' ';
  200. ' ANALYTIQUE * INTERVALLE * PROCHE
  201. * SIMULTANE * IRAM';
  202. ' -------------*-----------------*-----------
  203. -----*-------------';
  204. MESS ' HARM 1 0.489 * ' F1 ' * ' F2 ' * ' F3 ' * ' F4;
  205. MESS ' HARM 3 0.863 * ' G1 ' * ' G2 ' * ' G3 ' * ' G4;
  206. MESS ' ';
  207. MESS '*';
  208. '******** ON VERIFIE LE CHAMP DE PRESSION AUX POINTS
  209. A1,A5 ET A3';
  210. MESS '*';
  211. MESS ' ';
  212. ' ANALYTIQUE * INTERVALLE * PROCHE
  213. * SIMULTANE * IRAM';
  214. ' -------------*-----------------*-----------
  215. -----*-------------';
  216. MESS 'HARM 1 1.00 * ' X1 ' * ' X2 ' * ' X3' * ' X4 ;
  217. MESS 'HARM 1 1.81 * ' Y1 ' * ' Y2 ' * ' Y3' * ' Y4 ;
  218. MESS 'HARM 1 3.76 * ' Z1 ' * ' Z2 ' * ' Z3' * ' Z4 ;
  219. MESS 'HARM 3 1.00 * ' T1 ' * ' T2 ' * ' T3' * ' T4 ;
  220. MESS 'HARM 3 18.31 * ' U1 ' * ' U2 ' * ' U3' * ' U4 ;
  221. MESS 'HARM 3 201.72 * ' V1 ' * ' V2 ' * ' V3' * ' V4 ;
  222. TEMPS ;
  223.  
  224. *------------- code de bon fonctionnement ------------*
  225. SAUT 1 PAGE;
  226. ERR = TABL;
  227. ERR.1 = 100 * (ABS(0.489 -F1)/0.489);
  228. ERR.2 = 100 * (ABS(0.489 -F2)/0.489);
  229. ERR.3 = 100 * (ABS(0.489 -F3)/0.489);
  230. ERR.4 = 100 * (ABS(0.863 -G1)/0.863);
  231. ERR.5 = 100 * (ABS(0.863 -G2)/0.863);
  232. ERR.6 = 100 * (ABS(0.863 -G3)/0.863);
  233. ERR.7 = 100 * (ABS(1-X1)/1);
  234. ERR.8 = 100 * (ABS(1-X2)/1);
  235. ERR.9 = 100 * (ABS(1-X3)/1);
  236. ERR.10= 100 * (ABS(1.81 -Y1)/1.81);
  237. ERR.11= 100 * (ABS(1.81 -Y2)/1.81);
  238. ERR.12= 100 * (ABS(1.81 -Y3)/1.81);
  239. ERR.13= 100 * (ABS(3.76 -Z1)/3.76);
  240. ERR.14= 100 * (ABS(3.76 -Z2)/3.76);
  241. ERR.15= 100 * (ABS(3.76 -Z3)/3.76);
  242. ERR.16= 100 * (ABS(1-T1));
  243. ERR.17= 100 * (ABS(1-T2));
  244. ERR.18= 100 * (ABS(1-T3));
  245. ERR.19= 100 * (ABS(18.31 -U1)/18.31);
  246. ERR.20= 100 * (ABS(18.31 -U2)/18.31);
  247. ERR.21= 100 * (ABS(18.31 -U3)/18.31);
  248. ERR.22= 100 * (ABS(201.72-V1)/201.72);
  249. ERR.23= 100 * (ABS(201.72-V2)/201.72);
  250. ERR.24= 100 * (ABS(201.72-V3)/201.72);
  251. * AJOUT IRAM
  252. ERR.25 = 100 * (ABS(0.489 -F4)/0.489);
  253. ERR.26 = 100 * (ABS(0.863 -G4)/0.863);
  254. ERR.27 = 100 * (ABS(1 -X4)/1);
  255. ERR.28 = 100 * (ABS(1.81 -Y4)/1.81);
  256. ERR.29 = 100 * (ABS(3.76 -Z4)/3.76);
  257. ERR.30 = 100 * (ABS(1 -T4));
  258. ERR.31 = 100 * (ABS(18.31 -U4)/18.31);
  259. ERR.32 = 100 * (ABS(201.72-V4)/201.72);
  260.  
  261. I=0;
  262. REPETER BOUC1 32;
  263. I=I+1;
  264. SI (ERR.I <EG 1);
  265. ERRE 0;
  266. SINON;
  267. MESS 'erreur valeur ' I ' du tableau :' ERR.I;
  268. ERRE 5;
  269. FINSI;
  270. FIN BOUC1;
  271.  
  272. FIN;
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  
  279.  

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