本軟件使用Matlab函數(shù)和m文件編寫(xiě)而成。本程序使用有兩種方法

1、通過(guò)命令窗口輸入命令調(diào)用函數(shù)。
2、通過(guò)圖形用戶界面（模態(tài)圖，模態(tài)形狀，Campbell圖等）

使用過(guò)程：

1、運(yùn)行rotfeini.m文件，設(shè)定改程序的工作路徑（包括rotfe_ib文件夾及其所有的子文件夾）
2、在rotfe_ib文件夾下建立一個(gè)后綴名為.m的文件，建立轉(zhuǎn)子系統(tǒng)的模型（仿照tutorial.m和
   test.m文件）
    模型中包括

        a.節(jié)點(diǎn)
           NODES＝[x1 x2 ...],其中xi為節(jié)點(diǎn)的坐標(biāo)

         
        b.單元 
           ELEMENTS＝[1 2  30e-3 12e-3 1
                    2 3  32e-3 10e-3 1  格式為：[單元的第一個(gè)節(jié)點(diǎn) 單元的第二個(gè)節(jié)點(diǎn) 外直徑 內(nèi)直徑 單元材料表示號(hào)]
                    ...]
     

        c.材料
         MATERIALS＝[210e9 7800 0.3
                      70e9 3200 0.3      格式為：[楊氏模量 密度 泊松常數(shù)]
                      ... ...]  

        d.轉(zhuǎn)盤(pán)
         DISCS=[6  0.2 	12e-3 12e-3 1]   格式為：[圓盤(pán)所在位置的節(jié)點(diǎn)號(hào) 外直徑 內(nèi)直徑 厚度 材料編號(hào)]
  

        e.支承
         SPRINGS＝[b1 2e5 2.1e5 1e5 1e5
                   b2 2e5 2.1e5 1e5 1e5   格式為：[支承所在節(jié)點(diǎn)號(hào) Kxx Kyy Kxy Kyx]

        f.點(diǎn)質(zhì)量
        POINT_MASS=[5 0.01 0 0
                    ...  ...]       格式為：[節(jié)點(diǎn)編號(hào) 質(zhì)量 Jp Jd]
                   ...  ...]

        g.不平衡量
         UNBALANCE=[3 1e-3 0
                    ...  ...]   格式為：[不平衡所在的節(jié)點(diǎn) x方向的不平衡 y方向的不平衡]
                                其中x、y方向的不平衡分別為m*ex和m*ey 
                 
模型文件的格式如下：

  NODES=[0 0.3 0.5 0.6 0.8 1.0 ];                                                              
  ELEMENTS=[ 	...
      1 2  30e-3 12e-3 1                                                                                                                                            
      2 3  32e-3 10e-3 1
      3 4  22e-3 8e-3 1
      4 5  28e-3 6e-3 1
      5 6  30e-3 6e-3 1];                                                                          
  MATERIALS=[210e9 7800 0.3; 70e9 3200 0.3];     
  DISCS=[6 0.2 12e-3 12e-3 1];                                               
  SPRINGS=[1 2e3 2e3; 2 2e3 2e3];
  POINT_MASS=[5 0.01 0 0];

3、給出轉(zhuǎn)子的模型圖
命令格式:drawrot tutorial.m 其中tutorial.m為第二步中所建立的轉(zhuǎn)子系統(tǒng)模型
 另外：
     drawrot3D tutorial.m給出模型轉(zhuǎn)子的3維圖
     drawrot3D_sec tutorial.m 給出去掉軸1/4后的部分3維圖
     drawrot3D_spring tutorial.m 和drawrot3D_spring2 tutorial.m都給出帶有3維彈性支承的轉(zhuǎn)子模型

給出轉(zhuǎn)子的模型圖只是對(duì)建立轉(zhuǎn)子模型過(guò)程的校正，不是必要的一步。

4.建立有限元模型
 命令格式：
   R＝rotfe('tutorial');
   其中tutorial為轉(zhuǎn)子系統(tǒng)模型
   R為返回的轉(zhuǎn)子有限元結(jié)構(gòu)參數(shù)

5.畫(huà)出模態(tài)形狀并給出響應(yīng)的臨界轉(zhuǎn)速
  命令格式：
      modeplot('tutorial')        四者等價(jià)
      modeplot1('tutorial')
      modeplot2('tutorial')
      modeplot2f('tutorial')

6.給出campbell圖
  命令格式：
   campbell
  然后在彈出的對(duì)話框中將file改為tutotial即可


7.響應(yīng)的仿真分析
  利用matlab的simulink工具箱函數(shù)

   Rot=rotfe('tutorial')  
   [a,b,c,d]=rot2ss(ROt,w) 
   命令rot2ss將轉(zhuǎn)子的有限元模型變換成狀態(tài)空間模型
      空間狀態(tài)方程dx/dt=Ax+Bu
                  y=Cx+Du
      其中，前式表示轉(zhuǎn)子系統(tǒng)的狀態(tài)方程，u表示不平衡量，為該模塊的輸入量；
      y表示轉(zhuǎn)子的響應(yīng)，為該模塊的輸出量。


   rotsfuncw(t,x,u,flag,Rot)為解不平衡加速運(yùn)動(dòng)方程