tst_reduce2.m

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

%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')