MATLAB的转子振动计算代码

转子振动计算的matlab代码，matlab网站上给出的。

%demo_rotfe
%
%Demonstrating  some of Rotfe abilities
%
% Rotfe :  a simple finite element program for rotordynamics
%
%
% Author:  Izhak Bucher
% Email:   bucher@technion.ac.il
% Where:   Mechanical Engineering
%          Technion Israel Institute of Technology
%          Haifa 32000
%          Israel
%


%----------> Must type the below command before running this
%
%     rotfeini   <<<<<<-------------------------------------
%
%      % include all relevant folders in path <----

close all bdclose all
echo demo_rotfe 
echo on



modeplot2 rz4
h=findobj（‘tag‘‘modeplot2critspd‘）;
h1=findobj（‘tag‘‘modeplot2solid‘）;
h2=findobj（‘tag‘‘modeplot2N‘）;
h3=findobj（‘tag‘‘modeplot2modeN‘）;
h4=findobj（‘tag‘‘modeplot2section‘）;

set（h‘value‘2）;  % do critical speed analysis （not modal analysis）

modeplot2 go       % compute critical speeds and corresponding motion（s）

pause（1）

for q=1:5
    modeplot2 modeplus  % move to next mode （like pressing the + pushbutton）
    pause（3）            % wait a bit
end


set（h1‘value‘1） % solid   % choose solid display （effective only in 3d mode）
set（h2‘value‘1） % 3D      % switch to 3D
set（h3‘value‘1） % mode 1  % choose model # 1
modeplot2 display           % and refresh display  -- like changing the displayed mode

for q=1:5
    modeplot2 modeplus      % advnace to visialize a few modes
    pause（.3）
end

%%%%%%%%%%%%%%%%%%%%%%%%555
% campbell
%  Compute the campbell diagram of a model
%
close all
figure      % new figure
 campbell simple1s      % choose campbell-diagram  of model=‘simple1s‘ -see under models 
 campbell go            % compute
 
set（gca‘ylim‘[0 350]） % choose part of the range （there are high frequencies that are not interesting）

pause （5）;


close all

figure
 S=dir（‘models\*.m‘）;  % use the drawrot function to show some exisitng models
 for q=1:length（S）
      try
          drawrot（S（q）.name）;
          pause（.3）
      catch
          disp（sprintf（‘ skipped %s  this one （not a valid model‘S（q）.name））
      end
  end
  
close all
  for q=1:length（S）
      try
          drawrot3d_sec（S（q）.name）;
          pause（.3）
      catch
          disp（sprintf（‘ skipped %s  this one （not a valid model‘S（q）.name））
      end
  end
  
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6148  2009-11-19 09:51  Rotfe25x\.DS_Store

     文件       2319  2003-07-19 12:21  Rotfe25x\demo_rotfe.m

     文件      27648  2010-08-26 17:17  Rotfe25x\doc\Alan_RRA_reply_reduct.doc

     文件      18432  1999-08-14 11:24  Rotfe25x\doc\Convergence of the FE model.doc

     文件      53248  1999-08-14 11:24  Rotfe25x\doc\coordinates.doc

     文件     146432  1999-08-14 11:24  Rotfe25x\doc\Model order reduction in ROTFE 2.doc

     文件      25600  2007-07-11 20:47  Rotfe25x\doc\model reduction.doc

     文件      20992  1999-08-14 11:23  Rotfe25x\doc\New features and modifications in RotFE 2.doc

     文件     147968  1999-08-14 11:24  Rotfe25x\doc\rotexample.doc

     文件     170496  2010-08-26 17:20  Rotfe25x\doc\ROTFE20a.doc

     文件      81920  1999-08-14 11:24  Rotfe25x\doc\Rotor Dynamics software in Matlab.doc

     目录          0  2011-09-07 13:37  Rotfe25x\doc

     文件        876  2002-09-24 10:32  Rotfe25x\draw_plane.m

     文件       1505  2002-09-24 10:30  Rotfe25x\make_asym_frf.m

     文件      13248  2000-05-25 20:45  Rotfe25x\models\asym51.m

     文件      11345  2000-11-09 21:53  Rotfe25x\models\beit_sun.m

     文件       2553  2002-02-16 17:29  Rotfe25x\models\cobra0.m

     文件       2600  2002-02-16 19:20  Rotfe25x\models\cobra0a.m

     文件       2632  2002-02-16 19:57  Rotfe25x\models\cobra0at.m

     文件       2631  2002-02-16 20:28  Rotfe25x\models\cobra0at2.m

     文件       2670  2002-02-16 21:30  Rotfe25x\models\cobra0at3.m

     文件       2580  2002-02-16 15:19  Rotfe25x\models\cobra1.m

     文件        411  2006-12-13 16:53  Rotfe25x\models\isqrt32.m

     文件       9273  2009-10-14 21:50  Rotfe25x\models\lacm.m

     文件       8579  2009-10-15 08:11  Rotfe25x\models\lacm2.m

     文件       8741  2009-10-15 08:26  Rotfe25x\models\lacm3.m

     文件       8586  2009-10-31 10:46  Rotfe25x\models\lacm4a.m

     文件       8741  2009-10-15 08:39  Rotfe25x\models\lacm4b.m

     文件       8741  2009-11-01 15:19  Rotfe25x\models\lacm4c.m

     文件       6575  2000-07-09 23:49  Rotfe25x\models\manydisks.m

............此处省略192个文件信息