محاسبه فرکانس و شکل مودهای سازه

پس از محاسبه سختی و تشکیل ماتریس جرم می توانید با استفاده از کد زیر در متلب شکل مودها و فرکانس سازه را به صورت مرتب شده از کوچک به بزرگ بدست آورید

K=xlsread('stiffness Matrix ETABS.xlsx',1,'A1:L12');
M=xlsread('Mass Matrix dasti.xlsx',1,'A1:L12');
[V,D]=eig(K,M);
[W,k]=sort(diag(D));
V=V(:,k);
Factor=diag(V'*M*V);
Phi=V*inv(sqrt(diag(Factor)));
Omega=diag(sqrt(Phi'*K*Phi));

شکل مود های نرمال شده می باشد  Phi  همان فرکانس ها می باشد و  Omega