Code covered by the BSD License  

Highlights from
Rotor Dynamics toolbox (RotFE)

image thumbnail

Rotor Dynamics toolbox (RotFE)

by

Izhak Bucher (view profile)

 

Toolbox models rotating, elastic shafts with disks

tst_reduce2.m
%tst_sim1.m
%

Rot=rotfe('unbal2');
Rot.UNBALANCE=[5 0.3 0];
Rot=rmfield(Rot,'Reduct'); % no reduction
zeta=[ones(1,4)*.03 ];%ones(1,4)*.12 ones(1,4)*.12];
Rot.PROP_DAMP=zeta;		 % add 2% damping
%Rot.D=addpdamp(Rot.M,Rot.K,zeta);
Rot=rotfe(Rot);   % build new model withour reduction
Rot.W=00;
[a1 b1 c1 d1 bu1 du1]=rot2ss(Rot);

%sim_tutorial


Rot.Reduct.Type='modal';Rot.Reduct.Nmodes=4;
Rot=rotfe(Rot);   % build new model this time reduce order
Rot.W=00;
[a2 b2 c2 d2 bu2 du2]=rot2ss(Rot);

Rot.Reduct.Type='mixed';Rot.Reduct.Nmodes=2;
Rot=rotfe(Rot);   % build new model this time reduce order
Rot.W=00;
[a3 b3 c3 d3 bu3 du3]=rot2ss(Rot);
 
[y x t]=impulse(a1,b1,c1,d1,1);
[y2 x2 t2]=impulse(a2,b2,c2,d2,1,t);
[y3 x3 t3]=impulse(a3,b3,c3,d3,1,t);

plot(t,[y(:,1) y2(:,1) y3(:,1)],t,[y(:,2) y2(:,2) y3(:,2)]+1e-2)

legend('y1','y1-modal','y1-mixed','y2','y2-modal','y2-mixed')

Contact us