Télécharger vibr9.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : vibr9.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. **********************************************************************
  6. * *
  7. * Mots-clés : flambage, modes complexes, *
  8. * forces suiveuses, flottement *
  9. * *
  10. * TEST : VIBR9 *
  11. * *
  12. * Calcul des modes propres complexes d'une *
  13. * structure soumise a une force suiveuse *
  14. * (la structure est composée de deux barres *
  15. * reliees par deux ressorts spiraux) *
  16. * *
  17. **********************************************************************
  18. *
  19. Opti Echo 0;
  20. Opti dime 2 elem seg2 impi 0;
  21. *rem bp : "OPTI EPSI LINEAIRE;" inutile car matrices entrées a la main
  22.  
  23. *--- Affichage
  24. Graph = 'N';
  25. Affich = VRAI;
  26. *******************************************
  27. * Points reperes *
  28. *******************************************
  29. P1= 0 0;
  30. P2= 1 0;
  31. L1 = D 1 P1 P2;
  32. L2 = D 1 P2 P1;
  33. *
  34. *******************************************
  35. * Description du système *
  36. *******************************************
  37. *
  38. * Demi-longueur d'une barre
  39. a = 1;
  40. * Masse d'une barre
  41. m = 10;
  42. * Raideur des ressorts spiraux
  43. k = 20;
  44. * Variation de la force suiveuse
  45. Fmin = 0.;
  46. Fmax = 100.;
  47. PFs = 1.;
  48. *
  49. lFs= Prog Fmin PAS PFs Fmax;
  50. *
  51. *******************************************
  52. * Frontiere des domaines *
  53. *******************************************
  54. *
  55. Fs1 = 5*k/(4*a);
  56. Fs2 = 13*k/(4*a);
  57. *
  58. Si (Affich);
  59. Mess 'Fs 1 = ' Fs1;
  60. Mess 'Fs 2 = ' Fs2;
  61. Finsi;
  62. *
  63. *******************************************
  64. * Définition des grandeurs physiques *
  65. *******************************************
  66. *
  67. * --- Matrice Masse totale
  68. Mt = m*a*a*(MANU RIGI type MASSE L1 (MOTS ALFA) DUAL (MOTS FALF)
  69. (PROG 5. 1. 1. 1.));
  70.  
  71. * --- Matrice Raideur
  72. Ks = k*(MANU RIGI type RIGIDITE L1 (MOTS ALFA) DUAL (MOTS FALF)
  73. (PROG 2. -1. -1. 1.));
  74. * --- Matrice Raideur pour force suiveuse unitaire
  75. UK = 2*a*(MANU RIGI type RIGIDITE L1 (MOTS ALFA) DUAL (MOTS FALF)
  76. (PROG -1. 1. 0. 0.));
  77. * /////////////////////////////////////// *
  78. * Initialisation de l'algorithme *
  79. * /////////////////////////////////////// *
  80. *
  81. * Frequences propres
  82. Freqr = Table 'FREQ_REEL';
  83. Freqi = Table 'FREQ_IMAG';
  84. Repeter Freq1 4;
  85. Ifreq = &Freq1;
  86. Freqr.Ifreq = Prog;
  87. Freqi.Ifreq = Prog;
  88. Fin Freq1;
  89. *
  90. dlFs = dime lFs;
  91. *
  92. * /////////////////////////////////////// *
  93. * Itérations sur Fs *
  94. * /////////////////////////////////////// *
  95. *
  96. * --- Reperage de la force critique
  97. FSauv = 'N';
  98. *
  99. Repeter bloc1 dlFs;
  100. *
  101. i=&bloc1;
  102. Fs = extr lFs i;
  103. *
  104. * --- Matrice de Raideur totale
  105. Kt = Ks et (Fs*UK) ;
  106. *
  107. *******************************************
  108. * Résolution *
  109. *******************************************
  110. *
  111. BASC = VIBC Mt Kt ;
  112. *
  113. *******************************************
  114. * Traitement et stockage *
  115. *******************************************
  116. *
  117. *
  118. Si (AFFICH);
  119. Mess 'Iteration' I;
  120. Mess 'Force suiveuse : ' Fs ;
  121. Mess ' ';
  122. Mess '----------------------------------------------------------------
  123. --------------';
  124. Mess ' Mode Frequence Amort. P. Reelle P. Imaginaire
  125. Stabilité ';
  126. Finsi;
  127. *
  128. BMOD = BASC.'MODES';
  129. Repeter Freq2 4;
  130. Ifreq=&Freq2;
  131. MOD = BMOD.Ifreq;
  132. Si (EGA MOD.SOUSTYPE 'MODE_ANNULE');
  133. f1 = 0.;
  134. if1 =1.;
  135. Sinon;
  136. f1 = MOD.'FREQUENCE_REELLE';
  137. if1 =MOD.'FREQUENCE_IMAGINAIRE';
  138. Finsi;
  139. Freqr.Ifreq = Freqr.Ifreq et (Prog f1);
  140. Freqi.Ifreq = Freqi.Ifreq et (Prog if1);
  141. Si (Affich);
  142. Msg = 'STAB.';
  143. Si ((< if1 -1.D-10) et (>EG f1 0.));
  144. Msg = 'INST.';
  145. Si (EGA FSauv 'N');
  146. Fcrit = Fs - (PFs/2.);
  147. FSauv = 'O';
  148. Finsi;
  149. Finsi;
  150. Si (> f1 0.);
  151. * Valeurs propres relatives a lambda = i omega = ix2PIxf
  152. Mess Ifreq f1 (if1/f1) (-2*PI*if1) (2*PI*f1) Msg;
  153. Sinon;
  154. Si (EGA f 0.);
  155. Mess Ifreq f1 ' ----- ' (-2*PI*if1) (2*PI*f1) Msg;
  156. Finsi;
  157. Finsi;
  158. Finsi;
  159. Fin Freq2;
  160. Mess ' ';
  161. Mess ' ';
  162. *
  163. *
  164. * /////////////////////////////////////// *
  165. * Fin du calcul *
  166. * /////////////////////////////////////// *
  167. *
  168. Fin bloc1;
  169. *
  170. Si ( EGA Graph 'O');
  171. *
  172. TabM = Table;
  173. TabM.1 = 'NOLI';
  174. Titr 'Frequence (partie reelle)';
  175. Evol1 = Evol BLEU manu lFs (Freqr.1);
  176. Titr 'Frequence (partie imaginaire)';
  177. Evol2 = Evol BLEU manu lFs (Freqi.1);
  178. Titr 'Frequence dans le plan complexe';
  179. Evol3 = Evol BLEU manu (Freqr.1) (Freqi.1);
  180. Repeter Trac1 3;
  181. i=&Trac1+1;
  182. TabM.i ='NOLI';
  183. Titr 'Frequence (partie reelle)';
  184. Evol1 = Evol1 et (Evol BLEU manu lFs (Freqr.i));
  185. Titr 'Frequence (partie imaginaire)';
  186. Evol2 = Evol2 et (Evol BLEU manu lFs (Freqi.i));
  187. Titr 'Frequence dans le plan complexe';
  188. Evol3 = Evol3 et (Evol BLEU manu (Freqr.i) (Freqi.i));
  189. Fin Trac1;
  190. Dess Evol1 Titx 'Force (N)'
  191. Tity 'Re(F) (Hz)' TabM;
  192. Dess Evol2 Titx 'Force (N)'
  193. Tity 'Im(F) (Hz)' TabM;
  194. Dess Evol3 Titx 'Re(F) (Hz)'
  195. Tity 'Im(F) (Hz)' TabM;
  196. Finsi;
  197. *
  198. Mess ' Theorie Calcul Erreur ';
  199. Mess Fs1 Fcrit (Abs(Fcrit-Fs1)/Fs1);
  200. *Opti donn 5;
  201. *
  202. Si ((Abs(Fcrit-Fs1)/Fs1) < 0.05);
  203. Erre 0;
  204. Sinon;
  205. Erre 5;
  206. Finsi;
  207. *
  208. FIN;
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  

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