unbal2.m

来自「利用matalab编制的转子动力学计算软件。其中模型利用有限元建模」· M 代码 · 共 63 行

% ==> ELEMENTS 
%-------------------------------- 
% ==> dim 
%-------------------------------- 
% ==> M 
%-------------------------------- 
% ==> dof 
%-------------------------------- 
 dof =[3 4 5 6 9 10 11 12] ; 
%-------------------------------- 
% ==> W 
%-------------------------------- 
 W =10 ; 
%-------------------------------- 
% ==> NODES 
%-------------------------------- 
 NODES =[0:0.5:3] ; 
 NN=length(NODES);
  ELEMENTS=[];
 for q=1:NN-1,
     ELEMENTS=[  ELEMENTS ; q q+1 .02 0 1];
end
    
 %-------------------------------- 
% ==> MATERIALS 
%-------------------------------- 
 MATERIALS =[209000000000 7800 0.3;70e9 3000 .38 ; 40e9 5000 .29] ; 
%-------------------------------- 
% ==> DISCS 
%-------------------------------- 
 DISCS =[NN .4 .05 0.05 2] ; 
%-------------------------------- 
% ==> POINT_MASS 
%-------------------------------- 
 POINT_MASS =[] ; 
%-------------------------------- 
% ==> SPRINGS 
%-------------------------------- 
 SPRINGS =[] ; 
%-------------------------------- 
% ==> DASHPOTS 
%-------------------------------- 
 DASHPOTS =[] ; 
%-------------------------------- 
% ==> BCNodeDir 
%-------------------------------- 
 BCNodeDir =1 ; 
%-------------------------------- 
% ==> Bcord 
%-------------------------------- 
% ==> UNBALANCE 
%-------------------------------- 
 UNBALANCE =[NN 1.1 0] ; 
%-------------------------------- 
% ==> FNodeDir 
%-------------------------------- 
 FNodeDir =[NN+.1 NN+.3] ; 
%-------------------------------- 
% ==> RNodeDir 
%-------------------------------- 
 RNodeDir =[NN+.1 NN+.3] ; 
 %-------
 Reduct.Type='modal'; Reduct.Nmodes=4;