振型分解反應(yīng)譜法matlab

1、%本程序采用振型分解反應(yīng)譜法計(jì)算框架結(jié)構(gòu)水平地震力%采用KN.M單位%運(yùn)行本程序之前請(qǐng)運(yùn)行CYGD1.M和CYGD2.M求解框架水平側(cè)移剛度%本程序未考慮扭轉(zhuǎn)耦聯(lián)振動(dòng)，只能用于平面框架計(jì)算。求解所有振型。%結(jié)構(gòu)地震影響系數(shù)按高規(guī)3.3.8選取%地震作用和作用效應(yīng)按高規(guī)3.3.10計(jì)算clear %清理WORKSPACEk0= %各層框架側(cè)移剛度;m0=1.904 %各層質(zhì)量，重力荷載代表值/g2.6772.6772.6772.677*1.0e3./9.8;n1= 0.21712; %單榀框架地震力分配系數(shù) Tg=0.35; %特征周期(按規(guī)范選取)s=0.05; %阻尼比(按規(guī)范選取)r=0.

2、9; %衰減系數(shù)(按規(guī)范選取)y1=0.02; %阻尼比調(diào)整系數(shù)1(按規(guī)范選取)y2=1; %阻尼比調(diào)整系數(shù)2(按規(guī)范選取)amax=0.08; %水平地震最大影響系數(shù)(按規(guī)范選取)zjxu=0.7 %周期折減系數(shù)(按規(guī)范選取) cn=length(m0) %計(jì)算樓層數(shù)l=diag(ones(cn); m=diag(m0); %計(jì)算質(zhì)量矩陣ik=matrixju(k0,cn); %計(jì)算剛度矩陣x,d=eig(ik,m) %求解特征值和特征向量d=diag(sqrt(d) %求解結(jié)構(gòu)圓頻率T=zjxu*2*pi./d %求解結(jié)構(gòu)特征周期并作折減，折減系數(shù)0.7for i=1:cn; dl(i),

3、j=min(d); xgd(:,i)=x(:,j); d(j)=max(d)+1;endw=dl; %輸出結(jié)構(gòu)自振頻率x=xgd; for j=1:cn; %求解結(jié)構(gòu)振型參與系數(shù)和各質(zhì)點(diǎn)的水平相對(duì)位移x x(:,j)=x(:,j)/x(cn,j); zhcan(j)=(x(:,j)*m*l/(x(:,j)*m*x(:,j);endy=zhcan %輸出振型參與系數(shù)alpha=zeros(cn,1); %求解地震影響系數(shù)for i=1:cn;if T(i)0.1 alpha(i)=5.5*amax*T(i)+0.45*amax;else if T(i)5*Tg; alpha(i)=0.2.r-y

4、1*(T(i)-5*Tg)/y2*amax; else alpha(i)=(Tg/T(i).r*amax; end end endenda=alpha %輸出地震影響系數(shù)F=zeros(cn); %求解各振型下，各層水平地震力for i=1:cn; F(:,i)=x(:,i).*m0.*y(i).*a(i).*9.8;endf=F.*n1 %求解分配后地震力，并排序z=zeros(cn,cn);b=cn+1;for n=1:cn; z(b-n,:)=f(n,:);endp=z %輸出地震力v=zeros(cn,1); %求解水平地震力效應(yīng)值 for i=1:cn; v(i)=sum(p(i,:).2); end V=v.0.5 子程序 function kcju=matrixju(korc,cn)kcju=zeros(cn,cn);for i=1:cn-1; kcju(i,i)=korc(i)+korc(i+1); kc