Télécharger vibr3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : vibr3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test vibr3.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test vibr3.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. * VIBR3
  30. ******************************************************************
  31. * *
  32. * Mots-clés : Vibrations, calcul modal, 2D *
  33. * *
  34. * TEST VIBR3 *
  35. * *
  36. * CAS TEST DU 91/06/13 PROVENANCE : TEST *
  37. * *
  38. * POUTRE DE SECTION VARIABLE : VIBRATION DANS LE PLAN *
  39. * *
  40. * Test NAFEMS numero fv32 *
  41. * *
  42. * La structure est constituee d'une poutre de section *
  43. * variable. Sa longueur est de 10.0 m, sa plus grande largeur *
  44. * est de 5.0 m et sa plus petite largeur est de 1.0 m. *
  45. * *
  46. * On s'interesse aux vibrations dans le plan xOy, c'est-a- *
  47. * dire que le deplacemnt uz de tous les noeuds est nul. On *
  48. * travaille donc en mode plan (dimension 2). *
  49. * *
  50. * La poutre est bloquee le long de l'axe des ordonnees pour *
  51. * les deplacements ux et uy. *
  52. * *
  53. ******************************************************************
  54. OPTI ELEM QUA8;
  55. OPTI EPSI LINEAIRE;
  56.  
  57. TEMPS;
  58. *
  59. *================================================================*
  60. * geometries et maillages *
  61. *================================================================*
  62. *
  63. * Dimensions en metres
  64. *
  65. * Points
  66. *
  67. P1 = 0. -2.5;
  68. P2 = 0. 2.5;
  69. P3 = 10. 0.5;
  70. P4 = 10. -0.5;
  71. *
  72. * Droites
  73. *
  74. K = 1 ;
  75. N1 = 4 * K; N2 = 8 * K;
  76. *
  77. D12 = P1 DROI N1 P2;
  78. D23 = P2 DROI N2 P3;
  79. D34 = P3 DROI N1 P4;
  80. D41 = P4 DROI N2 P1;
  81. *
  82. SUR1 = DALL D12 D23 D34 D41 PLAN;
  83. *
  84. * Trace facultatif du maillage
  85. *
  86. SI (NEG GRAPH 'N');
  87. TITR ' VIBR3 : MAILLAGE';
  88. TRAC 'QUAL' SUR1;
  89. FINSI;
  90. *
  91. *================================================================*
  92. * modl- materiau - rigidite - masse *
  93. *================================================================*
  94. *
  95. MOD1 = MODE SUR1 MECANIQUE QUA8;
  96. MAT1 = MATE MOD1 YOUN 2.E11 NU 0.3 RHO 8.E3;
  97. RIG1 = RIGI MOD1 MAT1;
  98. MAS1 = MASS MOD1 MAT1;
  99. *
  100. *==============================================================*
  101. * conditions aux limites *
  102. *==============================================================*
  103. *
  104. CLI1 = BLOQ DEPL D12;
  105. RIG1 = RIG1 ET CLI1;
  106. *
  107. *==============================================================*
  108. * calcul des frequences propres *
  109. * et *
  110. * extraction des resultats *
  111. *==============================================================*
  112. *
  113. * Utilisation de l'operateur VIBR.
  114. *
  115. FRE1 = TABLE;
  116. FRE1.1 = 44.623;
  117. FRE1.2 = 130.03;
  118. FRE1.3 = 162.70;
  119. FRE1.4 = 246.05;
  120. FRE1.5 = 379.90;
  121. FRE1.6 = 391.44;
  122. *
  123. RESUL = VIBR INTERVALLE 40 450 BASSE 6 RIG1 MAS1;
  124. *
  125. *==============================================================*
  126. * affichache des resultats *
  127. *==============================================================*
  128. *
  129. MESS ' RESULTATS ';
  130. MESS ' --------- ';
  131. FRE2 = TABLE; MOD = TABLE; DEF = TABLE;ERG = TABLE;
  132. *
  133. I = 0;
  134. REPETER BLOC1 6;
  135. I = I + 1;
  136. FRE2.I = RESUL . MODES . I . FREQUENCE;
  137. ERG.I = 100 * (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  138. MESS ' MODE ' I;
  139. MESS ' ----------';
  140. MESS ' Frequence theorique :' FRE1.I 'Hz';
  141. MESS ' Frequence calculee :' FRE2.I 'Hz';
  142. MESS ' Soit un ecart de :' ERG.I '%';
  143. *
  144. * Extraction et trace des modes propres et calcul des deformees.
  145. * (facultatif)
  146. SI (NEG GRAPH 'N');
  147. MOD.I = RESUL . MODES . I . DEFORMEE_MODALE;
  148. DEF0 = DEFO SUR1 0. MOD.I;
  149. DEF.I = DEFO SUR1 MOD.I BLEU;
  150. TITR 'VIBR3 : MODES PROPRES';
  151. TRAC (DEF0 ET DEF.I);
  152. FINSI;
  153. FIN BLOC1;
  154. *
  155. *==============================================================*
  156. * code fonctionnement *
  157. *==============================================================*
  158. *
  159. ERGMAX = MAXI (PROG ERG.1 ERG.2 ERG.3 ERG.4 ERG.5 ERG.6);
  160.  
  161. SI (ERGMAX <EG 5.5);
  162. ERRE 0;
  163. SINON;
  164. ERRE 5;
  165. FINSI;
  166.  
  167. TEMPS;
  168.  
  169. FIN;
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  

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