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

cobra0at3.m
%cobra0.m                                                                                     
%                                                                                               
% by: I. Bucher 16-Feb-2002
                                     
                                                                                         
 NODES=[linspace(0,1.46,10) ] ;                                                              
                                                                                                
 N_NODES=length(NODES);      
 
  NODES=[linspace(0,3*1.46,30) ] ;                                                              
 N_NODES1=length(NODES);      

  TR35=3.;    
 for q=1:N_NODES1-1, 
	ELEMENTS(q,:)=[q q+1 TR35*25.4e-3 (TR35-2*0.052)*25.4e-3 1]; 
 end                                       
                                                                                                
                                                                                                
                                                                        
% MATERIALS=[E1 rho1 nu1; E2 rho2 nu2; ..... ]                                                          
                                                                                                
 MATERIALS=[72400000000 2770 0.3;
		210e9 7800 0.3; 70e9 3200 0.3];     
% DISCS=[node1 d_out d_in width material_no; % define disc #1                                   
                                                                  
 DISCS=[	1       TR35*[2.54e-2 2.54e-2-4e-3*2]   2.22*25.4e-3 1                                                                     
 		N_NODES TR35*[2.54e-2 2.54e-2-4e-3*2]	2.22*25.4e-3 1    ];                                                                         

   DISCS=[];                                                                                          
   
% SPRINGS=[node1 Kxx1 Kyy1 Kxy1 Kyx1 Ktt1 Kpp1 Ktp1 Kpt1   % bearing #1                         
                                                                                                
 SPRINGS=[];                                                                                      
 DASHPOTS=[];                                                                                   

 
 BCNodeDir=[1+[1:4]/10 N_NODES*[1:3]+.1 N_NODES*[1:3]+.3   ]';

  UNBALANCE=[2 1e-3 0];

  POINT_MASS = [ 1  1.0 2e-2 1.6e-2
                N_NODES 1.0 2e-2  1.6e-2
                2*N_NODES 1.0 2e-2  1.6e-2
                3*N_NODES 1.0 2e-2  1.6e-2
            ];
   
TorsionBC=1;                


Contact us