此程序用来计算单层隔振系统的固有频率

%此程序用来计算单层隔振系统的固有频率a=[2000,2000,2000,720,245,855];M=diag(a);%M表示柴油发电机组的质量矩阵b=[4040,2820,990,0,0,0];k=diag(b);%k表示单个弹性元件即橡胶支承的刚度矩阵%di表示第i个橡胶支承的安装点坐标矩阵d1=[0-0.33-0.7;0.3300.62;0.7-0.620];d2=[0-0.330;0.3300.62;0-0.620];d3=[0-0.330.7;0.3300.62;-0.7-0.620];d4=[00.330.7;-0.3300.62;-0.7-0.620];d5=[00.330;-0.3300.62;0-0.620];d6=[00.33-0.7;-0.3300.62;0.7-0.620];O=zeros(3,3);%产生一个3*3的零矩阵I=eye(3,3);%产生一个3*3的单位阵%ti表示第i个橡胶支承的局部坐标变换矩阵t1=[IO;d1I];t2=[IO;d2I];t3=[IO;d3I];t4=[IO;d4I];t5=[IO;d5I];t6=[IO;d6I];%K表示支承元件的刚度矩阵K=t1*k*t1'+t2*k*t2'+t3*k*t3'+t4*k*t4'+t5*k*t5'+t6*k*t6';G=inv(M);%对质量矩阵M求逆E=G*K;%E表示动力刚度矩阵[V,D]=eig(E);%计算特征值及特征向量DVforn=1:6f(n)=sqrt(D(n,n))/(2*pi);r(n)=50/f(n);endfr