Télécharger dyna6.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna6.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna6.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test dyna6.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. * DYNA6
  30. *
  31. * CAS TEST OPERATEUR VIBRATION
  32. * MODES DE FLEXION D'UNE POUTRE
  33. * ENCASTREE-LIBRE
  34. *
  35. * Cet exemple permet de tester les trois options
  36. * de l'opérateur VIBR :
  37. *
  38. * - par intervalle (INTERVALLE) oº on sépare les
  39. * modes par dichotomie sur un intervalle donné
  40. * - par proximité (PROCHE) oº on recherche les
  41. * modes proches de valeurs données
  42. * - par la méthode de Lanczos (SIMULTANE) oº on
  43. * projette sur un sous-espace
  44. *
  45. *
  46. * Le problème consiste à trouver les fréquences pour
  47. * le premier et le deuxième mode, d'une poutre
  48. * encastrée à une extrémité et libre à l'autre.
  49. *
  50. * Les valeurs obtenues par les trois méthodes sont
  51. * comparées aux valeurs théoriques.
  52. *
  53. *******************************************************
  54. TITR 'FLEXION D UNE POUTRE ENCASTREE-LIBRE';
  55. OPTI DIME 3 ELEM SEG2 MODE TRID IMPI 0;
  56. OPTI EPSI LINEAIRE;
  57. TEMPS ;
  58. DENSITE 0.1;
  59.  
  60. *---------------------- MAILLAGE ----------------------
  61. *----------- Q2 EST UN NOEUD POUR LE MODE NUMERO 2 ----
  62. Q1 = 0. 0. 0.;Q2 = (2.*0.783) 0. 0. ;Q3 = 2. 0. 0. ;
  63. LIG= Q1 D Q2 D Q3 ;
  64. OEIL1 = 5 5 5 ;
  65.  
  66. SI(NEG GRAPH 'N');
  67. TRAC OEIL1 LIG ;
  68. FINSI;
  69.  
  70. STAN= MODE LIG MECANIQUE POUT ;
  71.  
  72. *------------ MATERIAU ET CARATERISTIQUES -------------
  73. MATPL1=MATE STAN YOUNG 2.E11 NU 0.3 RHO 7800;
  74. CARPL1=CARA STAN SECT 0.01 INRY 6.25E-6 INRZ 1.E-3
  75. TORS 1.E-2;
  76. MATPL1=MATPL1 ET CARPL1;
  77.  
  78. *------------------ RIGIDITE ET MASSE -----------------
  79. RIGPL1 =RIGI STAN MATPL1;
  80. MASPLA1=MASS STAN MATPL1 ;
  81.  
  82. *------------------ ON BLOQUE Q1 ----------------------
  83. BLOQ1 =BLOQ DEPL ROTA Q1 ;
  84. RIGPLA1=RIGPL1 ET BLOQ1 ;
  85.  
  86. *______________________________________________________
  87. * PREMIERE OPTION (INTERVALLE)
  88. * ON SEPARE LES MODES PAR DICHOTOMIE
  89. * PUIS ON LANCE LES ITERATIONS INVERSES
  90. *______________________________________________________
  91. MODPLA1=VIBR INTE 10. 120. RIGPLA1 MASPLA1 ;
  92.  
  93. *______________________________________________________
  94. * DEUXIEME OPTION (PROCHE)
  95. * ON LANCE LES ITERATIONS INVERSES
  96. * A PARTIR DES FREQUENCES CONTENUES
  97. * DANS LE LISTREEL
  98. *______________________________________________________
  99. PR=PROG 17. 110.;
  100. MODPLA2=VIBR PROC PR RIGPLA1 MASPLA1 ;
  101.  
  102. *______________________________________________________
  103. * TROISIEME OPTION (SIMULTANE)
  104. * ON UTILISE LA METHODE DE LANCZOS
  105. * (PROJECTION SUR UN SOUS-ESPACE)
  106. *______________________________________________________
  107. MODPLA3=VIBR SIMU 1. 2 RIGPLA1 MASPLA1 ;
  108. *
  109. *---------------------- R{sultats ---------------------
  110. *
  111. OPTI ECHO 0;
  112. SAUTER PAGE;
  113. F1=TIRE MODPLA1 FREQ RANG 1;
  114. G1=TIRE MODPLA1 FREQ RANG 2;
  115. F2=TIRE MODPLA2 FREQ RANG 1;
  116. G2=TIRE MODPLA2 FREQ RANG 2;
  117. F3=TIRE MODPLA3 FREQ RANG 1;
  118. G3=TIRE MODPLA3 FREQ RANG 2;
  119. *
  120. MESS '*';
  121. MESS '******** ON VERIFIE LES FREQUENCES';
  122. MESS '*';
  123. ' ANALYTIQUE * INTERVALLE * PROCHE
  124. * SIMULTANE';
  125. ' -------------*-----------------*-----------
  126. -----*-------------';
  127. MESS ' MODE 1 17.71 * ' F1 '*' F2 '* ' F3;
  128. MESS ' MODE 2 110.98 * ' G1 '*' G2 '* ' G3;
  129. MESS ' ';
  130. MESS ' ';
  131. MESS '*';
  132. '******** ON VERIFIE LA DEFORMEE DU DEUXIEME MODE';
  133. MESS '*';
  134. MESS ' ';
  135. X1=EXTR (TIRE MODPLA1 DEPL RANG 2) UZ Q2;
  136. X2=EXTR (TIRE MODPLA2 DEPL RANG 2) UZ Q2;
  137. X3=EXTR (TIRE MODPLA3 DEPL RANG 2) UZ Q2;
  138. ' ANALYTIQUE * INTERVALLE * PROCHE
  139. * SIMULTANE';
  140. ' -------------*-----------------*-----------
  141. -----*-------------';
  142. MESS 'MODE 2 0.00 * ' X1 '*' X2 '* ' X3;
  143. TEMPS ;
  144. * CODE FONCTIONNEMENT
  145. FREF=17.71;
  146. GREF=110.98;
  147.  
  148. RESI1=100 * (ABS((F1-FREF)/FREF));
  149. RESI2=100 * (ABS((G1-GREF)/GREF));
  150.  
  151. SI((RESI1 <EG 1.) ET (RESI2 <EG 1.));
  152. ERRE 0;
  153. SINO;
  154. ERRE 5;
  155. FINSI;
  156.  
  157. FIN;
  158.  
  159.  
  160.  
  161.  
  162.  

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