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) ou on sépare les
  39. * modes par dichotomie sur un intervalle donné
  40. * - par proximité (PROCHE) ou on recherche les
  41. * modes proches de valeurs données
  42. * - par la méthode de Lanczos (SIMULTANE) ou on
  43. * projette sur un sous-espace
  44. * - par la méthode d Arnoldi (IRAM) ou on
  45. * projette sur un sous-espace
  46. *
  47. *
  48. * Le problème consiste à trouver les fréquences pour
  49. * le premier et le deuxième mode, d'une poutre
  50. * encastrée à une extrémité et libre à l'autre.
  51. *
  52. * Les valeurs obtenues par les trois méthodes sont
  53. * comparées aux valeurs théoriques.
  54. *
  55. *******************************************************
  56. TITR 'FLEXION D UNE POUTRE ENCASTREE-LIBRE';
  57. OPTI DIME 3 ELEM SEG2 MODE TRID IMPI 0;
  58. TEMPS ;
  59. DENSITE 0.1;
  60.  
  61. *---------------------- MAILLAGE ----------------------
  62. *----------- Q2 EST UN NOEUD POUR LE MODE NUMERO 2 ----
  63. Q1 = 0. 0. 0.;Q2 = (2.*0.783) 0. 0. ;Q3 = 2. 0. 0. ;
  64. LIG= Q1 D Q2 D Q3 ;
  65. OEIL1 = 5 5 5 ;
  66.  
  67. SI(NEG GRAPH 'N');
  68. TRAC OEIL1 LIG ;
  69. FINSI;
  70.  
  71. STAN= MODE LIG MECANIQUE POUT ;
  72.  
  73. *------------ MATERIAU ET CARATERISTIQUES -------------
  74. MATPL1=MATE STAN YOUNG 2.E11 NU 0.3 RHO 7800;
  75. CARPL1=CARA STAN SECT 0.01 INRY 6.25E-6 INRZ 1.E-3 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. * 4EME OPTION (IRAM)
  110. *______________________________________________________
  111. MODPLA4=VIBR IRAM 1. 2 RIGPLA1 MASPLA1 ;
  112. *
  113. *---------------------- Resultats ---------------------
  114. *
  115. OPTI ECHO 0;
  116. SAUTER PAGE;
  117. F1= MODPLA1 . 'MODES' . 1 . 'FREQUENCE';
  118. G1= MODPLA1 . 'MODES' . 2 . 'FREQUENCE' ;
  119. F2= MODPLA2 . 'MODES' . 1 . 'FREQUENCE' ;
  120. G2= MODPLA2 . 'MODES' . 2 . 'FREQUENCE' ;
  121. F3= MODPLA3 . 'MODES' . 1 . 'FREQUENCE' ;
  122. G3= MODPLA3 . 'MODES' . 2 . 'FREQUENCE' ;
  123. F4= MODPLA4 . 'MODES' . 1 . 'FREQUENCE' ;
  124. G4= MODPLA4 . 'MODES' . 2 . 'FREQUENCE' ;
  125. *
  126. MESS '*';
  127. MESS '******** ON VERIFIE LES FREQUENCES';
  128. MESS '*';
  129. MESS ' ANALYTIQUE * INTERVALLE * PROCHE * SIMULTANE * IRAM ';
  130. MESS ' -------------*-----------------*----------- -----*-------------';
  131. MESS ' MODE 1 17.71 * ' F1 ' * ' F2 ' * ' F3 ' * ' F4;
  132. MESS ' MODE 2 110.98 * ' G1 ' * ' G2 ' * ' G3 ' * ' G4;
  133. MESS ' ';
  134. MESS ' ';
  135. MESS '*';
  136. MESS '******** ON VERIFIE LA DEFORMEE DU DEUXIEME MODE';
  137. MESS '*';
  138. MESS ' ';
  139. X1=EXTR (MODPLA1 . 'MODES' . 2 . 'DEFORMEE_MODALE') UZ Q2;
  140. X2=EXTR (MODPLA2 . 'MODES' . 2 . 'DEFORMEE_MODALE') UZ Q2;
  141. X3=EXTR (MODPLA3 . 'MODES' . 2 . 'DEFORMEE_MODALE') UZ Q2;
  142. X4=EXTR (MODPLA4 . 'MODES' . 2 . 'DEFORMEE_MODALE') UZ Q2;
  143. MESS ' ANALYTIQUE * INTERVALLE * PROCHE * SIMULTANE * IRAM';
  144. MESS ' -------------*-----------------*----------- -----*-------------';
  145. MESS 'MODE 2 0.00 * ' X1 ' * ' X2 ' * ' X3 ' *' X4;
  146. TEMPS ;
  147. * CODE FONCTIONNEMENT
  148. FREF=17.71;
  149. GREF=110.98;
  150. * le point Q2 est proche du "noeud" du 2eme mode
  151. XREF= 0.;
  152. F = prog F1 F2 F3 F4;
  153. G = prog G1 G2 G3 G4;
  154. X = prog X1 X2 X3 X4;
  155. RES_F = 100 * (MAXI 'ABS' ((F-FREF)/FREF));
  156. RES_G = 100 * (MAXI 'ABS' ((G-GREF)/GREF));
  157. RES_X = 100 * (MAXI 'ABS' X);
  158.  
  159. SI((RES_F <EG 1.) ET (RES_G <EG 1.) ET (RES_X <EG 1.));
  160. ERRE 0;
  161. SINO;
  162. ERRE 5;
  163. FINSI;
  164.  
  165. FIN;
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  

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