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 LCAR LCARA G GRA;
  94. *
  95. *--------------- RIGIDITE ET MASSE --------------------
  96. *
  97. RAI=RIGI MODLIQ MAT;
  98. MAS=MASS MODLIQ MAT;
  99. *
  100. *______________________________________________________
  101. * ON CALCULE LE PREMIER MODE AVEC L'OPTION INTERVALLE
  102. * DE VIBR. ON RENORMALISE LE CHAMP DE DEPLACEMENT
  103. * (ET DE PRESSION) DE MANIERE A CE QUE LA PRESSION SOIT
  104. * EGALE A 1 A L'ALTITUDE 0 ET ON EXTRAIT LES PRESSIONS
  105. * AUX POINTS A2, A5 ET A3
  106. *______________________________________________________
  107. *
  108. MOD1=VIBR INTE 0. 1. RAI MAS;
  109. F1 =MOD1 . MODES . 1 . FREQUENCE;
  110. PRE=MOD1 . MODES . 1 . DEFORMEE_MODALE;
  111. PBAS=EXTR PRE P A2;
  112. PP=PRE/PBAS;
  113. X1=EXTR PP P A2;Y1=EXTR PP P A5;Z1=EXTR PP P A3;
  114. *
  115. *______________________________________________________
  116. * ON CALCULE LE PREMIER MODE AVEC L'OPTION PROCHE DE VIBR
  117. *______________________________________________________
  118. *
  119. MOD2=VIBR PROC (PROG 0.5) RAI MAS;
  120. F2= MOD2 . MODES . 1 . FREQUENCE;
  121. PRE= MOD2 . MODES . 1 . DEFORMEE_MODALE;
  122. PBAS=EXTR PRE P A2;
  123. PP=PRE/PBAS;
  124. X2=EXTR PP P A2;Y2=EXTR PP P A5;Z2=EXTR PP P A3;
  125. *
  126. *______________________________________________________
  127. * ON CALCULE LE PREMIER MODE AVEC L'OPTION SIMULTANE DE VIBR
  128. *______________________________________________________
  129. *
  130. MOD3=VIBR SIMU 0.5 1 RAI MAS;
  131. F3= MOD3 . MODES . 1 . FREQUENCE;
  132. PRE= MOD3 . MODES . 1 . DEFORMEE_MODALE;
  133. PBAS=EXTR PRE P A2;
  134. PP=PRE/PBAS;
  135. X3=EXTR PP P A2;Y3=EXTR PP P A5;Z3=EXTR PP P A3;
  136. *______________________________________________________
  137. * ON CALCULE LE PREMIER MODE AVEC L'OPTION IRAM DE VIBR
  138. *______________________________________________________
  139. *
  140. MOD4=VIBR SIMU 0.5 1 RAI MAS;
  141. F4= MOD4 . MODES . 1 . FREQUENCE;
  142. PRE= MOD4 . MODES . 1 . DEFORMEE_MODALE;
  143. PBAS=EXTR PRE P A2;
  144. PP=PRE/PBAS;
  145. X4=EXTR PP P A2;Y4=EXTR PP P A5;Z4=EXTR PP P A3;
  146. *
  147. *______________________________________________________
  148. * ON RECOMMENCE LES OPERATIONS CETTE FOIS
  149. * AVEC L'HARMONIQUE 3
  150. *______________________________________________________
  151. *
  152. OPTI MODE FOUR 3;
  153. MODLIQ1=MODE FLUID LIQUIDE LQU4;
  154. MODLIQ2=MODE SURFA LIQUIDE LSU2;
  155. *
  156. MODLIQ=MODLIQ1 Et MODLIQ2;
  157. *
  158. MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF LCAR LCARA G GRA;
  159. *
  160. RAI=RIGI MODLIQ MAT;
  161. MAS=MASS MODLIQ MAT;
  162. *
  163. MOD1=VIBR INTE 0. 1. RAI MAS;
  164. *------------------------------------------------------
  165. G1= MOD1 . MODES . 1 . FREQUENCE;
  166. PRE= MOD1 . MODES . 1 . DEFORMEE_MODALE;
  167. PBAS=EXTR PRE P A2;
  168. PP=PRE/PBAS;
  169. T1=EXTR PP P A2;U1=EXTR PP P A5;V1=EXTR PP P A3;
  170. *------------------------------------------------------
  171. MOD2=VIBR PROC (PROG 0.9) RAI MAS;
  172. *
  173. G2=MOD2 . MODES . 1 . FREQUENCE;
  174. PRE= MOD2 . MODES . 1 . DEFORMEE_MODALE;
  175. PBAS=EXTR PRE P A2;
  176. PP=PRE/PBAS;
  177. T2=EXTR PP P A2;U2=EXTR PP P A5;V2=EXTR PP P A3;
  178. *------------------------------------------------------
  179. MOD3=VIBR SIMU 0.9 1 RAI MAS;
  180. G3=MOD3 . MODES . 1 . FREQUENCE;
  181. PRE= MOD3 . MODES . 1 . DEFORMEE_MODALE;
  182. PBAS=EXTR PRE P A2;
  183. PP=PRE/PBAS;
  184. T3=EXTR PP P A2;U3=EXTR PP P A5;V3=EXTR PP P A3;
  185. *------------------------------------------------------
  186. MOD4=VIBR SIMU 0.9 1 RAI MAS;
  187. G4=MOD4 . MODES . 1 . FREQUENCE;
  188. PRE= MOD4 . MODES . 1 . DEFORMEE_MODALE;
  189. PBAS=EXTR PRE P A2;
  190. PP=PRE/PBAS;
  191. T4=EXTR PP P A2;U4=EXTR PP P A5;V4=EXTR PP P A3;
  192. *------------------------------------------------------
  193. SAUT PAGE; SAUT 2 LIGNE ;
  194. MESS '*';
  195. MESS '******** ON VERIFIE LES FREQUENCES';
  196. MESS '*';
  197. MESS ' ';
  198. MESS ' ANALYTIQUE * INTERVALLE * PROCHE * SIMULTANE * IRAM';
  199. MESS ' -------------*-----------------*----------- -----*-------------';
  200. MESS ' HARM 1 0.489 * ' F1 ' * ' F2 ' * ' F3 ' * ' F4;
  201. MESS ' HARM 3 0.863 * ' G1 ' * ' G2 ' * ' G3 ' * ' G4;
  202. MESS ' ';
  203. MESS '*';
  204. MESS '******** ON VERIFIE LE CHAMP DE PRESSION AUX POINTS A1,A5 ET A3';
  205. MESS '*';
  206. MESS ' ';
  207. MESS ' ANALYTIQUE * INTERVALLE * PROCHE * SIMULTANE * IRAM';
  208. MESS ' -------------*-----------------*----------- -----*-------------';
  209. MESS 'HARM 1 1.00 * ' X1 ' * ' X2 ' * ' X3' * ' X4 ;
  210. MESS 'HARM 1 1.81 * ' Y1 ' * ' Y2 ' * ' Y3' * ' Y4 ;
  211. MESS 'HARM 1 3.76 * ' Z1 ' * ' Z2 ' * ' Z3' * ' Z4 ;
  212. MESS 'HARM 3 1.00 * ' T1 ' * ' T2 ' * ' T3' * ' T4 ;
  213. MESS 'HARM 3 18.31 * ' U1 ' * ' U2 ' * ' U3' * ' U4 ;
  214. MESS 'HARM 3 201.72 * ' V1 ' * ' V2 ' * ' V3' * ' V4 ;
  215. TEMPS ;
  216.  
  217. *------------- code de bon fonctionnement ------------*
  218. SAUT 1 PAGE;
  219. ERR = TABL;
  220. ERR.1 = 100 * (ABS(0.489 -F1)/0.489);
  221. ERR.2 = 100 * (ABS(0.489 -F2)/0.489);
  222. ERR.3 = 100 * (ABS(0.489 -F3)/0.489);
  223. ERR.4 = 100 * (ABS(0.863 -G1)/0.863);
  224. ERR.5 = 100 * (ABS(0.863 -G2)/0.863);
  225. ERR.6 = 100 * (ABS(0.863 -G3)/0.863);
  226. ERR.7 = 100 * (ABS(1-X1)/1);
  227. ERR.8 = 100 * (ABS(1-X2)/1);
  228. ERR.9 = 100 * (ABS(1-X3)/1);
  229. ERR.10= 100 * (ABS(1.81 -Y1)/1.81);
  230. ERR.11= 100 * (ABS(1.81 -Y2)/1.81);
  231. ERR.12= 100 * (ABS(1.81 -Y3)/1.81);
  232. ERR.13= 100 * (ABS(3.76 -Z1)/3.76);
  233. ERR.14= 100 * (ABS(3.76 -Z2)/3.76);
  234. ERR.15= 100 * (ABS(3.76 -Z3)/3.76);
  235. ERR.16= 100 * (ABS(1-T1));
  236. ERR.17= 100 * (ABS(1-T2));
  237. ERR.18= 100 * (ABS(1-T3));
  238. ERR.19= 100 * (ABS(18.31 -U1)/18.31);
  239. ERR.20= 100 * (ABS(18.31 -U2)/18.31);
  240. ERR.21= 100 * (ABS(18.31 -U3)/18.31);
  241. ERR.22= 100 * (ABS(201.72-V1)/201.72);
  242. ERR.23= 100 * (ABS(201.72-V2)/201.72);
  243. ERR.24= 100 * (ABS(201.72-V3)/201.72);
  244. * AJOUT IRAM
  245. ERR.25 = 100 * (ABS(0.489 -F4)/0.489);
  246. ERR.26 = 100 * (ABS(0.863 -G4)/0.863);
  247. ERR.27 = 100 * (ABS(1 -X4)/1);
  248. ERR.28 = 100 * (ABS(1.81 -Y4)/1.81);
  249. ERR.29 = 100 * (ABS(3.76 -Z4)/3.76);
  250. ERR.30 = 100 * (ABS(1 -T4));
  251. ERR.31 = 100 * (ABS(18.31 -U4)/18.31);
  252. ERR.32 = 100 * (ABS(201.72-V4)/201.72);
  253.  
  254. I=0;
  255. REPETER BOUC1 32;
  256. I=I+1;
  257. SI (ERR.I <EG 1);
  258. ERRE 0;
  259. SINON;
  260. MESS 'erreur valeur ' I ' du tableau :' ERR.I;
  261. ERRE 5;
  262. FINSI;
  263. FIN BOUC1;
  264.  
  265. FIN;
  266.  
  267.  
  268.  
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  

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