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

《振型分解反應(yīng)譜法matlab.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

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=[263770%各層框架側(cè)移剛度263770263770263770123582];m0=[1.904%各層質(zhì)量，重力荷載代表值/g2.6772.6772.6772.677]*1.0e3./9.8;n1=0.21712;%單榀框架地震力分配系數(shù)Tg=0.35;%特征

2、周期(按規(guī)范選取)s=0.05;%阻尼比(按規(guī)范選取)r=0.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)特征周期并作折

3、減，折減系數(shù)0.7fori=1:cn;[dl(i),j]=min(d);xgd(:,i)=x(:,j);d(j)=max(d)+1;endw=dl;%輸出結(jié)構(gòu)自振頻率x=xgd;forj=1:cn;%求解結(jié)構(gòu)振型參與系數(shù)和各質(zhì)點(diǎn)的水平相對(duì)位移xx(:,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ù)fori=1:cn;ifT(i)<0.1alpha(i)=5.5*amax*T(i)+0.45*amax;elseifT

4、(i)5*Tg;alpha(i)=[0.2.^r-y1*(T(i)-5*Tg)/y2]*amax;elsealpha(i)=(Tg/T(i)).^r*amax;endendendenda=alpha%輸出地震影響系數(shù)F=zeros(cn);%求解各振型下，各層水平地震力fori=1:cn;F(:,i)=x(:,i).*m0.*y(i).*a(i).*9.8;endf=F.*n1%求解分配后地震力，并排序z=zeros(cn,cn);b=cn+1;forn=1:cn;z(b-n,:)=f(n,:);endp=z%輸出地震力v=z

5、eros(cn,1);%求解水平地震力效應(yīng)值fori=1:cn;v(i)=sum(p(i,:).^2);endV=v.^0.5子程序functionkcju=matrixju(korc,cn)kcju=zeros(cn,cn);fori=1:cn-1;kcju(i,i)=korc(i)+korc(i+1);kcju(i,i+1)=-korc(i+1);kcju(i+1,i)=-korc(i+1);endkcju(cn,cn)=korc(cn);