Télécharger la plaquette Cast3M


fr en

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.  
  56. OPTI MODE AXIS;
  57. OPTI ELEM QUA8;
  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. *  je remplace le repeter par une condition sur les lignes
  111. *REPETER BLOC1 11;
  112. *  I = I + 1;
  113. *  J = EXTR LIST1 I;
  114. *  K = EXTR LIST2 I;
  115. *  CL3 = RELA 1 DEPL DIRECTION DIR4 (NOEUD J)
  116. *         -   1 DEPL DIRECTION DIR3 (NOEUD K);
  117. *  CLI1 = CLI1 ET CL3;
  118. *FIN BLOC1;
  119.   CL3 = RELA 1 DEPL DIRECTION DIR4 D23
  120.          -   1 DEPL DIRECTION DIR3 (D41 INVE);
  121.   CLI1 = CLI1 ET CL3;
  122.  
  123. *
  124. LIST1 = LECT 15 17 19  21 23 25;
  125. LIST2 = LECT 26 28 29 30 31 27;
  126. *
  127. I = 0;
  128. * la c'est trop compliquer. On le garde mais on est dependant de la numerotation
  129. REPETER BLOC2 6;
  130.   I = I + 1;
  131.   J = EXTR LIST1 I;
  132.   K = EXTR LIST2 I;
  133.   CL4 = RELA 1 DEPL DIRECTION DIR4 (NOEUD J)
  134.          -   1 UR (NOEUD K);
  135.   CLI1 = CLI1 ET CL4;
  136. FIN BLOC2;
  137. *
  138. RIG1 = RIG1 ET CLI1;
  139. *
  140. *-----------------------------------------------------*
  141. *            calcul des frequences propres            *
  142. *-----------------------------------------------------*
  143. *
  144. * Utilisation de l'operateur VIBR.
  145. *
  146. FRE1 = TABLE;
  147. FRE1.1 = 369.91;
  148. FRE1.2 = 838.03;
  149. FRE1.3 = 1451.2;
  150. FRE1.4 = 2117.0;
  151. FRE1.5 = 2795.8;
  152. *
  153. LIST1 = PROG FRE1.1 FRE1.2 FRE1.3 FRE1.4 FRE1.5;
  154. *
  155. RESUL = VIBR PROC LIST1 RIG1 MAS1;
  156. *
  157. *-----------------------------------------------------*
  158. *         extraction des resultats
  159. *-----------------------------------------------------*
  160. *
  161. MESS ' RESULTATS FV42.DGIBI';
  162. MESS ' --------------------';
  164. MESS ' SPHERE EN VIBRATION RADIALE';
  166. *
  167. FRE2 = TABLE; ERG = TABLE;
  168. DEP = TABLE;  DEF = TABLE; EVO = TABLE;
  169. *
  170. I = 0;
  171. REPETER BLOC3 5;
  172. *
  173.   I = I + 1;
  174.   FRE2.I =  RESUL . MODES . I . FREQUENCE;
  175.   ERG.I = 100 * (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  176.   MESS ' MODE ' I;
  177.   MESS ' ----------';
  178.   MESS ' Frequence propre theorique :' FRE1.I 'Hz';
  179.   MESS ' Frequence propre calculee  :' FRE2.I 'Hz';
  180.   MESS '     Soit un ecart de : ' ERG.I '%';
  182. *
  183. * Trace facultatif des modes propres
  184. *
  185.   SI (NEG GRAPH 'N');
  186.      DEP.I =  RESUL  . MODES . I . DEFORMEE_MODALE;
  187.      DEF0  = DEFO SUR1 DEP.I 0. NOIR;
  188.      DEF.I = DEFO SUR1 DEP.I    VERT;
  189.      TITR ' Mode propre : Deplacement radial';
  190.      TRAC (DEF0 ET DEF.I);
  191.      EVO.I = EVOL BLEU CHPO DEP.I UR D23;
  192.      DESSIN EVO.I;
  193.   FINSI;
  194. *
  195. FIN BLOC3;
  196. *
  197. *-----------------------------------------------------*
  198. *                 code fonctionnement                 *
  199. *-----------------------------------------------------*
  200. *
  201. ERGMAX = MAXI (PROG ERG.1 ERG.2 ERG.3 ERG.4 ERG.5);
  202. *
  204. SI (ERGMAX <EG 2.);
  205.    ERRE 0;
  206. SINON;
  207.    ERRE 5;
  208. FINSI;
  209. *
  211. TEMPS;
  213. *
  214. FIN; 
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  

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