matlab振型分解法

基于matlab振型分解法求解结构地震响应摘要：建筑结构由地震引起的振动反应称为建筑结构的地震反应，它包括地震在结构中引起的内力，变形，位移，速度和加速度。结构地震反应可以通过结构动力学方法求解。本文利用Matlab程序求解结构地震作用引起的结构响应。关键词：动力特性振型分解法matlab地震响应Abstract:Architecturalstructurevibrationresponsescausedbyearthquakeknownastheseismicresponsesofstructures,includinginternal,deformation,displacement,velocityandaccelerationofstructurecausedbytheearthquake.Structuralseismicresponsecanbesolvedbystructuraldynamicsmethod.BasedonMatlabprogramsolvedynamicresponseofstructurecausedbyseismicresponsesinthisarticle.Keywords:dynamiccharacteristicsvibrationmodedecompositionmethodmatlabseismicresponse一．力学模型和运动方程在进行结构地震反应计算时，为了使问题简化，便于计算，需要把具体的结构体系，在满足工程计算精度要求的条件下，抽象为质点体系。根据达朗贝尔原理，可以推得单自由度体系的运动微分方程为）（tpkyycym可以写成为：）（ta22yyy其中，mc2/为单自由度结构体系的阻尼比，mk为单自由度结构体系的自振频率。当外部激励荷载为地震作用时，p（t）=-m)(tug，这时动力方程可以写成kyycym-m)(tug和yyy22)(tug根据平衡关系同样可以建立多自由度体系在地震荷载作用下的运动微分方程为)(tuMlKyyCyMgmm)(ty)(tym)(11tyk)(tyc其中M=Nmmm21为多自由度体系的质量矩阵，C=nnnnnnCCCCCCCCC212222111211为多自由度体系的阻尼矩阵，K=nnnnnnKKKKKKKKK212222111211为多自由度体系的刚度矩阵，l=111为n×1维单位列向量，y=nyyy21为多自由度体系的唯一向量，)(tug为输入到结构的地震速度。采用瑞雷阻尼假定，则阻尼矩阵为C=KM21，式中，1和2为与结构体系有关的常系数。二．动力特性的求解现考虑二自由度体系的情况运动方程为：该方程特解为：得到频率方程为：在matlab语言中，其求解过程是相当简单的。可以直接利用matlab数值工具箱中库函数命22222121ymykyk11212111ymykykmm1)(1tymm2)(2ty)sin()sin(2211tXytXy0222221122111mkkkmk令eig来实现。[x,d]=eig(ik,m);%ik和im分别为结构的刚度矩阵和质量矩阵d=diag(sqrt(d));fori=1:cn%cn为结构的层数，即质点数[d1(i),j]=min(d);xgd(:,i)=x(:,j);d(j)=max(d)+1;endw=d1；%所求自振频率x=xgd；%所求结构主振型三．振型分解法振型分解法就是利用各振型相互正交的特性，将原来藕联的微分方程组变成为若干相互独立的微分方程，从而使原来多自由度体系的动力计算变为若干单自由度体系的问题，在求得了各单自由度体系的解后，再将各个解进行组合，从而求得多自由度体系的地震反应。本文用matlab来计算一个三层的钢筋混凝土结构，结构的各层特性参数为：第一层到第三层质量m分别为2726kg，2760kg，2300kg，第一层到第三层刚度k分别为2.485e4N/M，1.921e4N/M，1.522e4N/M.地震波采用200galcentro波，采用周期为0.02s。经计算得该结构的自振频率为1=4.1041，2=10.4906，3=14.9514结构阵型矩阵为X=绘制该结构顶层位移，加速度时程曲线要调用plot函数即可。该过程实现如下：t=0:dt:ndzh*dt;subplot(2,2,1)plot(t,disp(3,:)*1000,'k-')Subplot(2,2,2)plot(t,accel(3,:),'k-')结构顶层的位移反应和加速度时程曲线如下图为了说明计算结果的可靠性，现用sap2000做了对比计算。计算得前三阶振型图如下计算所得的频率差别从下面图可以知道0246810121416123频率阶次频率值sap2000计算频率matlab计算频率从图上我们可以看到matlab计算的结果偏大。Sap2000计算得顶层位移时程曲线和顶层加速度时程曲线，可以看出时程曲线相差不大。四．小结对于对称性比较好的框架结构，把结构计算模型假设为层模型，利用matlab中的eig函数，能很方便的计算结构振型。利用matlab编制的振型分解法程序，可以计算结构的动力响应，通过与sap2000计算比较，动力响应基本吻合。参考文献：[1]徐赵东等matlab语言在建筑抗震工程中的应用，科学出版社。[2]北京金土木软件技术有限公司，中国建筑科学研究院，Sap2000中文版使用指南。人民交通出版社，2006.[3]陈以新，matlab在结构动力分析中的应用，《福建建筑》.1999.2.[4]徐荣桥，结构分析的有限元法与matlab程序设计，人民交通出版社，2001.12.