zhxcyxsh.m

matlab编写结构计算振型分解法

function [qx,a]=zhxcyxsh(B,b,ww,dn)
%求解振型参与系数qx，地震影响系数a---输入变量：振型矩阵B，
%                                              质量b
%                                              结构层数dn
%                                    返回变量：振型参与系数qx，地震影响系数a

Tg=0.4;amax=0.16;%求解地震影响系数所需参数,查规范

r=0.9;n1=0.02;n2=1.0; %一般结构的阻尼参数

l=diag(ones(dn));b=diag(b);

ww=diag(ww);T=2*pi*inv(ww);T=diag(T);%求解自振周期T

for j=1:dn
    qx(j)=(B(:,j))'*diag(b)/((B(:,j))'*diag(B(:,j))*diag(b));
    if T(j)<=0.1
        a(j)=(n2*amax-0.45*amax)*10*T(j)+0.45*amax;
    elseif T(j)>0.1&T(j)<=Tg
        a(j)=n2*amax;
    elseif T(j)>Tg&T(j)<=(5*Tg)
        a(j)=(Tg/T(j))^r*n2*amax;
    else
        a(j)=(n2*0.2^r-n1*(T(j)-5*Tg))*amax;
    end
end