Code covered by the BSD License  

Highlights from
Rotor Dynamics toolbox (RotFE)

image thumbnail

Rotor Dynamics toolbox (RotFE)

by

 

Toolbox models rotating, elastic shafts with disks

ofer1.m
% model.m - creat the database for the FE model                                     
%====================================================================================           
%   (1)                                                                                         
% define nodal locations from left side to right                                                                       
% NODES=[z1 z2 z3 ... ]                                                                         
% units=[Length] ([Meters])                                                                     
%                                                                                               
%====================================================================================           
  
 NODES=[0 12 29.775 47.55 (57.1:9.55:143.05) 1.60825e+002 178.6 190.6] ;   
 NODES=NODES*1e-3; %in [m]
%====================================================================================           
%   (2)                                                                                         
% define ELEMENTS                                                                               
% elements=[node1 node2 d_out d_in material_no ;  % element #1                                  
%   ...    % element #2                                                                         
%    ... ]     % etc.                                                                           
%                                                                                               
% node1, node2 : integer indices of entries in NODES                                            
% d_out, d_in  : [length][M] outer and inner diameters of a hollow shaft                        
% material _no : integer index. defines approp. row in the MATERIAL matrix (table)              
%====================================================================================           
s=size(NODES,2)
dout=1.9050e-002;  % shaft diameter [m] ;  
din=0;
mat=1;
for i=1:s-1
    ELEMENTS(i,:)=[i (i+1) dout din mat ];                                                                          
end
     
                                      
%====================================================================================           
%   (3)                                                                                         
% define material sets                                                                          
% MATERIALS=[E1 rho1 nu1; E2 rho2 nu2; ..... ]                                                          
% units=   E - [Pa] [psi] (young modulus)                                                       
%  rho [Kg/m^3] [lb/in^3] (density)  
%  nu [.] poison ration (0.3 for most metals)                                                           
%                                                                                               
%====================================================================================           
                                                                                                
 MATERIALS=[
                        210e9 7800 0.3;         %shaft
                        210e9 7800 0.3            %disc
                    ];     
%====================================================================================           
% DISCS=[node1 d_out d_in width material_no cylinder_mass;....]
%====================================================================================           
 
 DISCS=[ 2 0.062  1.9050e-002 0.012 2 0.1128;...
        16 0.062  1.9050e-002 0.012 2 0.1128];                                                                         

 %====================================================================================           
% SPRINGS=[node1 Kxx1 Kzz1 Kxz1 Kzx1;....]                     
%====================================================================================           
% spr1=size([p1 p2 p3],2)+1;
% spr2=size([p1 p2 p3 p4 p5 p6 p7 p8],2)+1;

  SPRINGS=[	4 596 	596   0 0 ;
           	14 596 	596  0 0];   

%====================================================================================           
% DASHPOTS=[node1 xi;.....]
%====================================================================================           
 DASHPOTS=[2  0.03;      
           15  0.03];   


%====================================================================================           
% UNBALANCE=[node1 ux1 uz1 ; node2  ux2 uz2; ... nodeQ uxQ uzQ]   node1=1,2, ..
  UNBALANCE=[2 1 0;
             16 1 0];

  %==================================================================================== 
           
% AMB=[node1 node2]   if node2==0 point force otherwise node2-node1==1 and force is spread between
  AMB=[11 12];
            

  %==================================================================================== 
%   (11)                                                                                         
% Point mass specification
%
% POINT_MASS=[node1 m Jp Jd]
%  node1=1,2, ..  m-mass [kg] , Jp,Jd - moment of inertia [Kg-m^2]
% nd1=POINT_MASS(:,1); 	% find node
% m=POINT_MASS(:,2); 	% masses
% Jp=POINT_MASS(:,3);  Jd=POINT_MASS(:,4); 

%  POINT_MASS=[2 1e-3 1e3 1e3];
  POINT_MASS=[];

Contact us