(最新整理)结构动力学振型分析

1、(完整)结构动力学振型分析(完整)结构动力学振型分析 编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)结构动力学振型分析）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快 业绩进步，以下为(完整)结构动力学振型分析的全部内容。matalab作业某三层钢筋混凝土结构，结构的各层特性参数为：第一层到第三层质量m分别为2400kg，1200

2、kg，1200kg，第一层到第三层刚度k分别为3.3104n/m,1.1104n/m,0。664n/m。地震采用acc_elcentro_0.34g,采样周期为0.02。m3=1200kgk3=0.66104n/m。m2=1200kgk2=1.1104n/mm1=2400kgk1=3。3*104n/m用振型分解法求解结构地震反应的matlab层序如下，编制该程序的程序框图以下所示开始输入地震参数和结构参数结束求解结构的地震反应计算振型参与系数输出结果计算单自由度体系的地震反应计算结构振型与自振型频率%振型分解法求解结构地震反应；主程序clearclc%地震波数据xs=20。287； dzhbo

3、=load(acc_elcentro_0。34g_0.02s.txt); ag=dzhbo0.01xs; dt=0.02； ndzh=400； cn=3; %cn为结构的层数，即质点数m0=2。4 1.2 1.2*1e+3； %结构各层质量 k0=3。3 1。1 0.66*1e+5； 结构各层刚度l=diag(ones（cn); m=diag（m0）； %计算质量矩阵ik=matrixju（k0,cn）; 计算刚度矩阵x,d=eig（ik，m）; 结构动力特性求解d=diag（sqrt（d）; 求解结构圆频率for i=1：cn; d1（i)，j=min(d)； xgd（:,i)=x（：,j）

4、； d(j)=max(d)+1； end %以此循环对所求频率和振型进行排序 w=d1; 所求自振频率x=xgd; 所求结构主振型a1=2*w(1）w（2）*（0.05*w（2）0。07w（1）)/（w（2）2w（1)2)； a2=2*（0.07w（2）-0.05w(1）/(w（2）2w（1)2); for j=1：cn x（：，j）=x(：，j)/x(cn,j）; znb0（j)=（a1+a2*w（j)2)/2/w(j）； zhcan（j）=(x(：,j）)m*l/（x（:,j)）m*x（：,j）； 求解振型参数 dlt（j,：)，dltacceler（j,:）=zxzj（znb0（j),w

5、（j），ag); end%求解结构各层的地震反应for i=1：cn； disp1=0; accel1=0; for j=1：cn disp0=zhcan(j)*dlt(j,:）x（i，j); accel0=zhcan（j）dltacceler（j,：）*x（i，j）； disp1=disp1+disp0; accel1=accel1+accel0; end disp（i,：）=disp1; accel（i，:）=accel1； end t=0:dt:ndzhdt；subplot（2,2，1） plot(t，disp（3，:)1e+3，k-) %subplot（2,2,2） plot（t，ac

6、cel(3，:）,k） 子程序用于求解单自由度结构体系的地震反应function bx，acceler=zxzj(znb，w，dag)dt=0.02;n=400;x(1）=0；dx（1)=0;ddx(1)=0;s=1+znb*dtw+w2*dt2/6; 中间参数sfor i=1:n a（i）=x（i)+dx(i）dt+ddx(i)dt2/3； b(i)=dx（i)+ddx(i)*dt/2; ddx(i+1)=-1*（dag（i+1)+1*znbwb(i)+w2*a（i)）/s;%加速度 dx(i+1)=b（i）+ddx（i+1）dt/2;速度 x(i+1）=a(i)+ddx(i+1)*dt2/

7、6;%位移endbx=x;acceler=ddx;%子程序刚度和阻尼矩阵的聚合function kcju=matrixju（korc,cn）kcju=zeros（cn)；for i=1：cn1 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);经程序求解，该结构的自振频率为w= 4.8683 11。0666 15。4541结构的振型矩阵为 x=0.1634 0.9238 2。76040.5691 -1.2267 -3。34231。0000 1.0000 1.