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) (PROG 5. 1. 1. 1.));
  69.  
  70. * --- Matrice Raideur
  71. Ks = k*(MANU RIGI type RIGIDITE L1 (MOTS ALFA) DUAL (MOTS FALF) (PROG 2. -1. -1. 1.));
  72. * --- Matrice Raideur pour force suiveuse unitaire
  73. UK = 2*a*(MANU RIGI type RIGIDITE L1 (MOTS ALFA) DUAL (MOTS FALF) (PROG -1. 1. 0. 0.));
  74. * /////////////////////////////////////// *
  75. * Initialisation de l'algorithme *
  76. * /////////////////////////////////////// *
  77. *
  78. * Frequences propres
  79. Freqr = Table 'FREQ_REEL';
  80. Freqi = Table 'FREQ_IMAG';
  81. Repeter Freq1 4;
  82. Ifreq = &Freq1;
  83. Freqr.Ifreq = Prog;
  84. Freqi.Ifreq = Prog;
  85. Fin Freq1;
  86. *
  87. dlFs = dime lFs;
  88. *
  89. * /////////////////////////////////////// *
  90. * Itérations sur Fs *
  91. * /////////////////////////////////////// *
  92. *
  93. * --- Reperage de la force critique
  94. FSauv = 'N';
  95. *
  96. Repeter bloc1 dlFs;
  97. *
  98. i=&bloc1;
  99. Fs = extr lFs i;
  100. *
  101. * --- Matrice de Raideur totale
  102. Kt = Ks et (Fs*UK) ;
  103. *
  104. *******************************************
  105. * Résolution *
  106. *******************************************
  107. *
  108. BASC = VIBC Mt Kt ;
  109. *
  110. *******************************************
  111. * Traitement et stockage *
  112. *******************************************
  113. *
  114. *
  115. Si (AFFICH);
  116. Mess 'Iteration' I;
  117. Mess 'Force suiveuse : ' Fs ;
  118. Mess ' ';
  119. Mess '---------------------------------------------------------------- --------------';
  120. Mess ' Mode Frequence Amort. P. Reelle P. Imaginaire Stabilité ';
  121. Finsi;
  122. *
  123. BMOD = BASC.'MODES';
  124. Repeter Freq2 4;
  125. Ifreq=&Freq2;
  126. MOD = BMOD.Ifreq;
  127. Si (EGA MOD.SOUSTYPE 'MODE_ANNULE');
  128. f1 = 0.;
  129. if1 =1.;
  130. Sinon;
  131. f1 = MOD.'FREQUENCE_REELLE';
  132. if1 =MOD.'FREQUENCE_IMAGINAIRE';
  133. Finsi;
  134. Freqr.Ifreq = Freqr.Ifreq et (Prog f1);
  135. Freqi.Ifreq = Freqi.Ifreq et (Prog if1);
  136. Si (Affich);
  137. Msg = 'STAB.';
  138. Si ((< if1 -1.D-10) et (>EG f1 0.));
  139. Msg = 'INST.';
  140. Si (EGA FSauv 'N');
  141. Fcrit = Fs - (PFs/2.);
  142. FSauv = 'O';
  143. Finsi;
  144. Finsi;
  145. Si (> f1 0.);
  146. * Valeurs propres relatives a lambda = i omega = ix2PIxf
  147. Mess Ifreq f1 (if1/f1) (-2*PI*if1) (2*PI*f1) Msg;
  148. Sinon;
  149. Si (EGA f 0.);
  150. Mess Ifreq f1 ' ----- ' (-2*PI*if1) (2*PI*f1) Msg;
  151. Finsi;
  152. Finsi;
  153. Finsi;
  154. Fin Freq2;
  155. Mess ' ';
  156. Mess ' ';
  157. *
  158. *
  159. * /////////////////////////////////////// *
  160. * Fin du calcul *
  161. * /////////////////////////////////////// *
  162. *
  163. Fin bloc1;
  164. *
  165. Si ( EGA Graph 'O');
  166. *
  167. TabM = Table;
  168. TabM.1 = 'NOLI';
  169. Titr 'Frequence (partie reelle)';
  170. Evol1 = Evol BLEU manu lFs (Freqr.1);
  171. Titr 'Frequence (partie imaginaire)';
  172. Evol2 = Evol BLEU manu lFs (Freqi.1);
  173. Titr 'Frequence dans le plan complexe';
  174. Evol3 = Evol BLEU manu (Freqr.1) (Freqi.1);
  175. Repeter Trac1 3;
  176. i=&Trac1+1;
  177. TabM.i ='NOLI';
  178. Titr 'Frequence (partie reelle)';
  179. Evol1 = Evol1 et (Evol BLEU manu lFs (Freqr.i));
  180. Titr 'Frequence (partie imaginaire)';
  181. Evol2 = Evol2 et (Evol BLEU manu lFs (Freqi.i));
  182. Titr 'Frequence dans le plan complexe';
  183. Evol3 = Evol3 et (Evol BLEU manu (Freqr.i) (Freqi.i));
  184. Fin Trac1;
  185. Dess Evol1 Titx 'Force (N)' Tity 'Re(F) (Hz)' TabM;
  186. Dess Evol2 Titx 'Force (N)' Tity 'Im(F) (Hz)' TabM;
  187. Dess Evol3 Titx 'Re(F) (Hz)' Tity 'Im(F) (Hz)' TabM;
  188. Finsi;
  189. *
  190. Mess ' Theorie Calcul Erreur ';
  191. Mess Fs1 Fcrit (Abs(Fcrit-Fs1)/Fs1);
  192. *Opti donn 5;
  193. *
  194. Si ((Abs(Fcrit-Fs1)/Fs1) < 0.05);
  195. Erre 0;
  196. Sinon;
  197. Erre 5;
  198. Finsi;
  199. *
  200. FIN;
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  

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