Télécharger vibr2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : vibr2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test vibr2.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test vibr2.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 PSC ;
  23. SINO ;
  24. OPTI ECHO 0 ;
  25. FINSI ;
  26.  
  27. SAUT PAGE;
  28. *********************************************************************
  29. * VIBR2
  30. *********************************************************************
  31. * *
  32. * Mots-clés : Vibrations, calcul modal, axisymetrique *
  33. * *
  34. * TEST VIBR2 *
  35. * *
  36. * CAS TEST DU 91/06/13 PROVENANCE : TEST *
  37. * *
  38. * SPHERE EPAISSE : VIBRATION RADIALE UNIFORME *
  39. * *
  40. * Test NAFEMS , test numero fv42 *
  41. * *
  42. * On se propose de calculer les frequences et *
  43. * modes propres d'une sphere epaisse. On ne s'inte- *
  44. * resse qu'aux modes de vibration radiale uniforme. *
  45. * *
  46. * On ne maille qu'un secteur de la sphere et l'on *
  47. * travaille en mode axisymetrique. *
  48. * *
  49. * Pour obtenir les conditions aux limites corres- *
  50. * pondantes, on utilise l'operateur RELA (les noeuds *
  51. * ayant le meme rayon ont le meme deplacement radial) *
  52. * *
  53. *********************************************************************
  54.  
  55. OPTI MODE AXIS;
  56. OPTI ELEM QUA8;
  57. OPTI EPSI LINEAIRE;
  58.  
  59. *-----------------------------------------------------*
  60. * geometrie - maillage *
  61. *-----------------------------------------------------*
  62. *
  63. ALPHA = 5.; COSALPHA = COS ALPHA; SINALPHA = SIN ALPHA;
  64. *
  65. P0 = 0. 0.;
  66. P1 = 1.8 ( 1.8 * SINALPHA);
  67. P2 = 1.8 (-1.8 * SINALPHA);
  68. P3 = 6.0 (-6.0 * SINALPHA);
  69. P4 = 6.0 ( 6.0 * SINALPHA);
  70. *
  71. D23 = P2 DROI 5 P3;
  72. *
  73. SUR1 = D23 ROTA 1 (2 * ALPHA) P0;
  74. D41 = COTE 3 SUR1;
  75. *
  76. ELIM (SUR1 ET D23 ET D41) 0.001;
  77. *
  78. SI (NEG GRAPH 'N');
  79. TITR 'VIBR2 : MAILLAGE';
  80. TRAC 'QUAL' SUR1;
  81. FINSI;
  82. *
  83. *-----------------------------------------------------*
  84. * MODE - materiau - rigidite - masse *
  85. *-----------------------------------------------------*
  86. *
  87. MOD1 = MODE SUR1 MECANIQUE QUA8;
  88.  
  89. MAT1 = MATE MOD1 YOUN 2.E11 NU 0.3 RHO 8E3;
  90. RIG1 = RIGI MOD1 MAT1;
  91. MAS1 = MASS MOD1 MAT1;
  92. *
  93. *-----------------------------------------------------*
  94. * conditions aux limites *
  95. *-----------------------------------------------------*
  96. DIR1 = (-1 * SINALPHA) (-1 * COSALPHA);
  97. DIR2 = (-1 * SINALPHA) COSALPHA;
  98. *
  99. CL1 = BLOQ DEPL DIRECTION DIR1 D23;
  100. CL2 = BLOQ DEPL DIRECTION DIR2 D41;
  101. CLI1 = CL1 ET CL2;
  102. *
  103. DIR3 = COSALPHA (-1 * SINALPHA);
  104. DIR4 = COSALPHA SINALPHA;
  105. *
  106. LIST1 = LECT 27 22 17 23 18 24 19 25 20 26 21;
  107. LIST2 = LECT 1 6 2 7 3 8 4 9 5 10 11;
  108. *
  109. I = 0;
  110. REPETER BLOC1 11;
  111. I = I + 1;
  112. J = EXTR LIST1 I;
  113. K = EXTR LIST2 I;
  114. CL3 = RELA 1 DEPL DIRECTION DIR4 (NOEUD J)
  115. - 1 DEPL DIRECTION DIR3 (NOEUD K);
  116. CLI1 = CLI1 ET CL3;
  117. FIN BLOC1;
  118. *
  119. LIST1 = LECT 27 17 18 19 20 21;
  120. LIST2 = LECT 28 12 13 14 15 16;
  121. *
  122. I = 0;
  123. REPETER BLOC2 6;
  124. I = I + 1;
  125. J = EXTR LIST1 I;
  126. K = EXTR LIST2 I;
  127. CL4 = RELA 1 DEPL DIRECTION DIR4 (NOEUD J)
  128. - 1 UR (NOEUD K);
  129. CLI1 = CLI1 ET CL4;
  130. FIN BLOC2;
  131. *
  132. RIG1 = RIG1 ET CLI1;
  133. *
  134. *-----------------------------------------------------*
  135. * calcul des frequences propres *
  136. *-----------------------------------------------------*
  137. *
  138. * Utilisation de l'operateur VIBR.
  139. *
  140. FRE1 = TABLE;
  141. FRE1.1 = 369.91;
  142. FRE1.2 = 838.03;
  143. FRE1.3 = 1451.2;
  144. FRE1.4 = 2117.0;
  145. FRE1.5 = 2795.8;
  146. *
  147. LIST1 = PROG FRE1.1 FRE1.2 FRE1.3 FRE1.4 FRE1.5;
  148. *
  149. RESUL = VIBR PROC LIST1 RIG1 MAS1;
  150. *
  151. *-----------------------------------------------------*
  152. * extraction des resultats
  153. *-----------------------------------------------------*
  154. *
  155. MESS ' RESULTATS FV42.DGIBI';
  156. MESS ' --------------------';
  157. MESS ' SPHERE EN VIBRATION RADIALE';
  158. *
  159. FRE2 = TABLE; ERG = TABLE;
  160. DEP = TABLE; DEF = TABLE; EVO = TABLE;
  161. *
  162. I = 0;
  163. REPETER BLOC3 5;
  164. *
  165. I = I + 1;
  166. FRE2.I = RESUL . MODES . I . FREQUENCE;
  167. ERG.I = 100 * (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  168. MESS ' MODE ' I;
  169. MESS ' ----------';
  170. MESS ' Frequence propre theorique :' FRE1.I 'Hz';
  171. MESS ' Frequence propre calculee :' FRE2.I 'Hz';
  172. MESS ' Soit un ecart de : ' ERG.I '%';
  173. *
  174. * Trace facultatif des modes propres
  175. *
  176. SI (NEG GRAPH 'N');
  177. DEP.I = RESUL . MODES . I . DEFORMEE_MODALE;
  178. DEF0 = DEFO SUR1 DEP.I 0. NOIR;
  179. DEF.I = DEFO SUR1 DEP.I VERT;
  180. TITR ' Mode propre : Deplacement radial';
  181. TRAC (DEF0 ET DEF.I);
  182. EVO.I = EVOL BLEU CHPO DEP.I UR D23;
  183. DESSIN EVO.I;
  184. FINSI;
  185. *
  186. FIN BLOC3;
  187. *
  188. *-----------------------------------------------------*
  189. * code fonctionnement *
  190. *-----------------------------------------------------*
  191. *
  192. ERGMAX = MAXI (PROG ERG.1 ERG.2 ERG.3 ERG.4 ERG.5);
  193. *
  194. SI (ERGMAX <EG 2.);
  195. ERRE 0;
  196. SINON;
  197. ERRE 5;
  198. FINSI;
  199. *
  200. TEMPS;
  201. *
  202. FIN;
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  

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