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=TIRE MOD1 FREQ RANG 1;
  111. PRE=TIRE MOD1 DEPL RANG 1;
  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
  118. * PROCHE DE VIBR
  119. *______________________________________________________
  120. *
  121. MOD2=VIBR PROC (PROG 0.5) RAI MAS;
  122. F2=TIRE MOD2 FREQ RANG 1;
  123. PRE=TIRE MOD2 DEPL RANG 1;
  124. PBAS=EXTR PRE P A2;
  125. PP=PRE/PBAS;
  126. X2=EXTR PP P A2;Y2=EXTR PP P A5;Z2=EXTR PP P A3;
  127. *
  128. *______________________________________________________
  129. * ON CALCULE LE PREMIER MODE AVEC L'OPTION SIMULTANE
  130. * DE VIBR
  131. *______________________________________________________
  132. *
  133. MOD3=VIBR SIMU 0.5 1 RAI MAS;
  134. F3=TIRE MOD3 FREQ RANG 1;
  135. PRE=TIRE MOD2 DEPL RANG 1;
  136. PBAS=EXTR PRE P A2;
  137. PP=PRE/PBAS;
  138. X3=EXTR PP P A2;Y3=EXTR PP P A5;Z3=EXTR PP P A3;
  139. *
  140. *______________________________________________________
  141. * ON RECOMMENCE LES OPERATIONS CETTE FOIS
  142. * AVEC L'HARMONIQUE 3
  143. *______________________________________________________
  144. *
  145. OPTI MODE FOUR 3;
  146. MODLIQ1=MODE FLUID LIQUIDE LQU4;
  147. MODLIQ2=MODE SURFA LIQUIDE LSU2;
  148. *
  149. MODLIQ=MODLIQ1 Et MODLIQ2;
  150. *
  151. MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF
  152. LCAR LCARA G GRA;
  153. *
  154. RAI=RIGI MODLIQ MAT;
  155. MAS=MASS MODLIQ MAT;
  156. *
  157. MOD1=VIBR INTE 0. 1. RAI MAS;
  158. *------------------------------------------------------
  159. G1=TIRE MOD1 FREQ RANG 1;
  160. PRE=TIRE MOD1 DEPL RANG 1;
  161. PBAS=EXTR PRE P A2;
  162. PP=PRE/PBAS;
  163. T1=EXTR PP P A2;U1=EXTR PP P A5;V1=EXTR PP P A3;
  164. *------------------------------------------------------
  165. MOD2=VIBR PROC (PROG 0.9) RAI MAS;
  166. *
  167. G2=TIRE MOD2 FREQ RANG 1;
  168. PRE=TIRE MOD2 DEPL RANG 1;
  169. PBAS=EXTR PRE P A2;
  170. PP=PRE/PBAS;
  171. T2=EXTR PP P A2;U2=EXTR PP P A5;V2=EXTR PP P A3;
  172. *------------------------------------------------------
  173. MOD3=VIBR SIMU 0.9 1 RAI MAS;
  174. G3=TIRE MOD3 FREQ RANG 1;
  175. PRE=TIRE MOD3 DEPL RANG 1;
  176. PBAS=EXTR PRE P A2;
  177. PP=PRE/PBAS;
  178. T3=EXTR PP P A2;U3=EXTR PP P A5;V3=EXTR PP P A3;
  179. *------------------------------------------------------
  180. SAUT PAGE; SAUT 2 LIGNE ;
  181. MESS '*';
  182. MESS '******** ON VERIFIE LES FREQUENCES';
  183. MESS '*';
  184. MESS ' ';
  185. ' ANALYTIQUE * INTERVALLE * PROCHE
  186. * SIMULTANE';
  187. ' -------------*-----------------*-----------
  188. -----*-------------';
  189. MESS ' HARM 1 0.489 * ' F1 '*' F2 '* ' F3;
  190. MESS ' HARM 3 0.863 * ' G1 '*' G2 '* ' G3;
  191. MESS ' ';
  192. MESS '*';
  193. '******** ON VERIFIE LE CHAMP DE PRESSION AUX POINTS
  194. A1,A5 ET A3';
  195. MESS '*';
  196. MESS ' ';
  197. ' ANALYTIQUE * INTERVALLE * PROCHE
  198. * SIMULTANE';
  199. ' -------------*-----------------*-----------
  200. -----*-------------';
  201. MESS 'HARM 1 1.00 * ' X1 '*' X2 '* ' X3 ;
  202. MESS 'HARM 1 1.81 * ' Y1 '*' Y2 '* ' Y3 ;
  203. MESS 'HARM 1 3.76 * ' Z1 '*' Z2 '* ' Z3 ;
  204. MESS 'HARM 3 1.00 * ' T1 '*' T2 '* ' T3 ;
  205. MESS 'HARM 3 18.31 * ' U1 '*' U2 '* ' U3 ;
  206. MESS 'HARM 3 201.72 * ' V1 '*' V2 '* ' V3 ;
  207. TEMPS ;
  208.  
  209. *------------- code de bon fonctionnement ------------*
  210. SAUT 1 PAGE;
  211. ERR = TABL;
  212. ERR.1 = 100 * (ABS(0.489 -F1)/0.489);
  213. ERR.2 = 100 * (ABS(0.489-F2)/0.489);
  214. ERR.3 = 100 * (ABS(0.489 -F3)/0.489);
  215. ERR.4 = 100 * (ABS(0.863 -G1)/0.863);
  216. ERR.5 = 100 * (ABS(0.863 -G2)/0.863);
  217. ERR.6 = 100 * (ABS(0.863 -G3)/0.863);
  218. ERR.7 = 100 * (ABS(1-X1)/1);
  219. ERR.8 = 100 * (ABS(1-X2)/1);
  220. ERR.9 = 100 * (ABS(1-X3)/1);
  221. ERR.10= 100 * (ABS(1.81 -Y1)/1.81);
  222. ERR.11= 100 * (ABS(1.81 -Y2)/1.81);
  223. ERR.12= 100 * (ABS(1.81 -Y3)/1.81);
  224. ERR.13= 100 * (ABS(3.76 -Z1)/3.76);
  225. ERR.14= 100 * (ABS(3.76 -Z2)/3.76);
  226. ERR.15= 100 * (ABS(3.76 -Z3)/3.76);
  227. ERR.16= 100 * (ABS(1-T1));
  228. ERR.17= 100 * (ABS(1-T2));
  229. ERR.18= 100 * (ABS(1-T3));
  230. ERR.19= 100 * (ABS(18.31 -U1)/18.31);
  231. ERR.20= 100 * (ABS(18.31 -U2)/18.31);
  232. ERR.21= 100 * (ABS(18.31 -U3)/18.31);
  233. ERR.22= 100 * (ABS(201.72-V1)/201.72);
  234. ERR.23= 100 * (ABS(201.72-V2)/201.72);
  235. ERR.24= 100 * (ABS(201.72-V3)/201.72);
  236.  
  237. I=0;
  238. REPETER BOUC1 24;
  239. I=I+1;
  240. SI (ERR.I <EG 1);
  241. ERRE 0;
  242. SINON;
  243. MESS 'erreur valeur ' I ' du tableau :' ERR.I;
  244. ERRE 5;
  245. FINSI;
  246. FIN BOUC1;
  247.  
  248. FIN;
  249.  
  250.  
  251.  
  252.  
  253.  

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