Télécharger la plaquette Cast3M


fr en

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. TEMPS ;
  48.  
  49. *--------------------------- MAILLAGE -----------------
  50. *----------------- ON NE MAILLE QUE LE FLUIDE ---------
  51. *
  52. A1=0.995 0.;A2=1.005 0.;A3=1.005 2.;A4=0.995 2.;
  53. A5=1.005 1.2;A6=0.995 1.2;
  54. L1=D 1 A1 A2;
  55. L2=A2 D 24 A5 D 16 A3;
  56. L3=D 1 A3 A4;
  57. L4=A4 D 16 A6 D 24 A1;
  58. FLUID=DALL L1 L2 L3 L4;
  59. *
  60. SI (NEG GRAPH 'N');
  61.   TRAC 'QUAL' FLUID ;
  62. FINSI;
  63. *
  64. OPTI ELEM SEG2;
  65. SURFA=D 1 A3 A4;
  66. *
  67. SI (NEG GRAPH 'N');
  68.   TRAC 'QUAL' SURFA ;
  69. FINSI;
  70. *
  71. MAIL=FLUID ET SURFA;
  72. *
  73. SI (NEG GRAPH 'N');
  74.   TRAC 'QUAL' MAIL;
  75. FINSI;
  76. *
  77. *---------- DEFINITION DE QUELQUES CONSTANTES ---------
  78. *
  79. ROREF=1.;CEREF=1.;LCARA=1.;
  80. RO0=1000.;CS=1200.;GRA=9.81;
  81. *
  82. *-- MODELE ET OBJET AFFECTE FLUIDE ET SURFACE LIBRE ---
  83. *-------- POUR L'HARMONIQUE DE FOURIER 1 --------------
  84. *
  85. OPTI MODE FOUR 1;
  86. MODLIQ1=MODE FLUID LIQUIDE LQU4;
  87. MODLIQ2=MODE SURFA LIQUIDE LSU2;
  88. MODLIQ=MODLIQ1 ET MODLIQ2;
  89. *
  90. *--------------- MATERIAU LIQUIDE ---------------------
  91. *
  92. MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF LCAR LCARA G GRA;
  93. *
  94. *--------------- RIGIDITE ET MASSE --------------------
  95. *
  96. RAI=RIGI MODLIQ MAT;
  97. MAS=MASS MODLIQ MAT;
  98. *
  99. *______________________________________________________
  100. * ON CALCULE LE PREMIER MODE AVEC L'OPTION INTERVALLE
  101. * DE VIBR. ON RENORMALISE LE CHAMP DE DEPLACEMENT
  102. * (ET DE PRESSION) DE MANIERE A CE QUE LA PRESSION SOIT
  103. * EGALE A 1 A L'ALTITUDE 0 ET ON EXTRAIT LES PRESSIONS
  104. * AUX POINTS A2, A5 ET A3
  105. *______________________________________________________
  106. *
  107. MOD1=VIBR INTE 0. 1. RAI MAS;
  108. F1 =MOD1 . MODES . 1 . FREQUENCE;
  109. PRE=MOD1 . MODES . 1 . DEFORMEE_MODALE;
  110. PBAS=EXTR PRE P A2;
  111. PP=PRE/PBAS;
  112. X1=EXTR PP P A2;Y1=EXTR PP P A5;Z1=EXTR PP P A3;
  113. *
  114. *______________________________________________________
  115. * ON CALCULE LE PREMIER MODE AVEC L'OPTION PROCHE DE VIBR
  116. *______________________________________________________
  117. *
  118. MOD2=VIBR PROC (PROG 0.5) RAI MAS;
  119. F2= MOD2 . MODES . 1 . FREQUENCE;
  120. PRE= MOD2 . MODES . 1 . DEFORMEE_MODALE;
  121. PBAS=EXTR PRE P A2;
  122. PP=PRE/PBAS;
  123. X2=EXTR PP P A2;Y2=EXTR PP P A5;Z2=EXTR PP P A3;
  124. *
  125. *______________________________________________________
  126. * ON CALCULE LE PREMIER MODE AVEC L'OPTION SIMULTANE DE VIBR
  127. *______________________________________________________
  128. *
  129. MOD3=VIBR SIMU 0.5 1 RAI MAS;
  130. F3= MOD3 . MODES . 1 . FREQUENCE;
  131. PRE= MOD3 . MODES . 1 . DEFORMEE_MODALE;
  132. PBAS=EXTR PRE P A2;
  133. PP=PRE/PBAS;
  134. X3=EXTR PP P A2;Y3=EXTR PP P A5;Z3=EXTR PP P A3;
  135. *______________________________________________________
  136. * ON CALCULE LE PREMIER MODE AVEC L'OPTION IRAM DE VIBR
  137. *______________________________________________________
  138. *
  139. MOD4=VIBR IRAM 0.5 1 RAI MAS;
  140. F4= MOD4 . MODES . 1 . FREQUENCE;
  141. PRE= MOD4 . MODES . 1 . DEFORMEE_MODALE;
  142. PBAS=EXTR PRE P A2;
  143. PP=PRE/PBAS;
  144. X4=EXTR PP P A2;Y4=EXTR PP P A5;Z4=EXTR PP P A3;
  145. *
  146. *______________________________________________________
  147. * ON RECOMMENCE LES OPERATIONS CETTE FOIS
  148. * AVEC L'HARMONIQUE 3
  149. *______________________________________________________
  150. *
  151. OPTI MODE FOUR 3;
  152. MODLIQ1=MODE FLUID LIQUIDE LQU4;
  153. MODLIQ2=MODE SURFA LIQUIDE LSU2;
  154. *
  155. MODLIQ=MODLIQ1 Et MODLIQ2;
  156. *
  157. MAT=MATE MODLIQ RHO RO0 RORF ROREF CSON CS CREF CEREF LCAR LCARA G GRA;
  158. *
  159. RAI=RIGI MODLIQ MAT;
  160. MAS=MASS MODLIQ MAT;
  161. *
  162. MOD1=VIBR INTE 0. 1. RAI MAS;
  163. *------------------------------------------------------
  164. G1= MOD1 . MODES . 1 . FREQUENCE;
  165. PRE= MOD1  . MODES . 1 . DEFORMEE_MODALE;
  166. PBAS=EXTR PRE P A2;
  167. PP=PRE/PBAS;
  168. T1=EXTR PP P A2;U1=EXTR PP P A5;V1=EXTR PP P A3;
  169. *------------------------------------------------------
  170. MOD2=VIBR PROC (PROG 0.9) RAI MAS;
  171. *
  172. G2=MOD2 . MODES . 1 . FREQUENCE;
  173. PRE= MOD2  . MODES . 1 . DEFORMEE_MODALE;
  174. PBAS=EXTR PRE P A2;
  175. PP=PRE/PBAS;
  176. T2=EXTR PP P A2;U2=EXTR PP P A5;V2=EXTR PP P A3;
  177. *------------------------------------------------------
  178. MOD3=VIBR SIMU 0.9 1 RAI MAS;
  179. G3=MOD3 . MODES . 1 . FREQUENCE;
  180. PRE= MOD3  . MODES . 1 . DEFORMEE_MODALE;
  181. PBAS=EXTR PRE P A2;
  182. PP=PRE/PBAS;
  183. T3=EXTR PP P A2;U3=EXTR PP P A5;V3=EXTR PP P A3;
  184. *------------------------------------------------------
  185. MOD4=VIBR IRAM 0.9 1 RAI MAS;
  186. G4=MOD4 . MODES . 1 . FREQUENCE;
  187. PRE= MOD4  . MODES . 1 . DEFORMEE_MODALE;
  188. PBAS=EXTR PRE P A2;
  189. PP=PRE/PBAS;
  190. T4=EXTR PP P A2;U4=EXTR PP P A5;V4=EXTR PP P A3;
  191. *------------------------------------------------------
  192. SAUT PAGE; SAUT 2 LIGNE ;
  193. MESS '*';
  194. MESS '********    ON VERIFIE LES FREQUENCES';
  195. MESS '*';
  196. MESS ' ';
  197. MESS '         ANALYTIQUE   *   INTERVALLE    *     PROCHE *    SIMULTANE * IRAM';
  198. MESS '         -------------*-----------------*----------- -----*-------------';
  199. MESS ' HARM 1    0.489      *   ' F1 ' * ' F2 ' * ' F3 ' * ' F4;
  200. MESS ' HARM 3    0.863      *   ' G1 ' * ' G2 ' * ' G3 ' * ' G4;
  201. MESS ' ';
  202. MESS '*';
  203. MESS '********   ON VERIFIE LE CHAMP DE PRESSION AUX POINTS A1,A5 ET A3';
  204. MESS '*';
  205. MESS ' ';
  206. MESS '        ANALYTIQUE   *   INTERVALLE    *     PROCHE *    SIMULTANE * IRAM';
  207. MESS '        -------------*-----------------*----------- -----*-------------';
  208. MESS 'HARM 1     1.00      *  ' X1 ' * ' X2 ' * ' X3' * ' X4 ;
  209. MESS 'HARM 1     1.81      *  ' Y1 ' * ' Y2 ' * ' Y3' * ' Y4 ;
  210. MESS 'HARM 1     3.76      *  ' Z1 ' * ' Z2 ' * ' Z3' * ' Z4 ;
  211. MESS 'HARM 3     1.00      *  ' T1 ' * ' T2 ' * ' T3' * ' T4 ;
  212. MESS 'HARM 3    18.31      *  ' U1 ' * ' U2 ' * ' U3' * ' U4 ;
  213. MESS 'HARM 3   201.72      *  ' V1 ' * ' V2 ' * ' V3' * ' V4 ;
  214. TEMPS ;
  215.  
  216. *------------- code de bon fonctionnement ------------*
  217. SAUT 1 PAGE;
  218. ERR   = TABL;
  219. ERR.1 = 100 * (ABS(0.489 -F1)/0.489);
  220. ERR.2 = 100 * (ABS(0.489 -F2)/0.489);
  221. ERR.3 = 100 * (ABS(0.489 -F3)/0.489);
  222. ERR.4 = 100 * (ABS(0.863 -G1)/0.863);
  223. ERR.5 = 100 * (ABS(0.863 -G2)/0.863);
  224. ERR.6 = 100 * (ABS(0.863 -G3)/0.863);
  225. ERR.7 = 100 * (ABS(1-X1)/1);
  226. ERR.8 = 100 * (ABS(1-X2)/1);
  227. ERR.9 = 100 * (ABS(1-X3)/1);
  228. ERR.10= 100 * (ABS(1.81 -Y1)/1.81);
  229. ERR.11= 100 * (ABS(1.81 -Y2)/1.81);
  230. ERR.12= 100 * (ABS(1.81 -Y3)/1.81);
  231. ERR.13= 100 * (ABS(3.76 -Z1)/3.76);
  232. ERR.14= 100 * (ABS(3.76 -Z2)/3.76);
  233. ERR.15= 100 * (ABS(3.76 -Z3)/3.76);
  234. ERR.16= 100 * (ABS(1-T1));
  235. ERR.17= 100 * (ABS(1-T2));
  236. ERR.18= 100 * (ABS(1-T3));
  237. ERR.19= 100 * (ABS(18.31 -U1)/18.31);
  238. ERR.20= 100 * (ABS(18.31 -U2)/18.31);
  239. ERR.21= 100 * (ABS(18.31 -U3)/18.31);
  240. ERR.22= 100 * (ABS(201.72-V1)/201.72);
  241. ERR.23= 100 * (ABS(201.72-V2)/201.72);
  242. ERR.24= 100 * (ABS(201.72-V3)/201.72);
  243. * AJOUT IRAM
  244. ERR.25 = 100 * (ABS(0.489 -F4)/0.489);
  245. ERR.26 = 100 * (ABS(0.863 -G4)/0.863);
  246. ERR.27 = 100 * (ABS(1     -X4)/1);
  247. ERR.28 = 100 * (ABS(1.81  -Y4)/1.81);
  248. ERR.29 = 100 * (ABS(3.76  -Z4)/3.76);
  249. ERR.30 = 100 * (ABS(1     -T4));
  250. ERR.31 = 100 * (ABS(18.31 -U4)/18.31);
  251. ERR.32 = 100 * (ABS(201.72-V4)/201.72);
  252.  
  253. I=0;
  254. REPETER BOUC1 32;
  255. I=I+1;
  256.   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.  
  275.  
  276.  
  277.  
  278.  

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