资源描述：
《哈工大结构动力学大作业2012春》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、结构动力学大作业对于如下结构，是研究质量块的质量变化和在简支梁上位置的变化对整个系统模态的影响。以上为一个简支梁结构。集中质量块放于梁上，质量块距简支梁的左端点距离为L.将该简支梁简化为欧拉伯努利梁，并离散为N个单元。每个单元有两个节点，四个自由度。单元的节点位移可表示为：则单元内一点的挠度可计作：带入边界条件：建立了单元位移模式后，其动能势能均可用节点位移表示。单元的动能为：其中m为单元质量阵，并有：带入公式后积分可得：单元势能可表示为其中K为单元刚度矩阵，并有以上为单元类型矩阵，通过定义全局位移矩阵，可以得到系统刚度矩阵和系统质量矩
2、阵。当集中质量块加到简支梁上时，可以认为系统的刚度矩阵不变，而质量矩阵在有基础上有所增加。当集中质量块位于第i个单元内，距单元左节点的位移为，则其动能可表示为：则质量块附加到总体刚度矩阵的子矩阵尾为：例如当质量为的质量块位于单元的右节点上时，其质量阵可表示为根据以上理论，编写matlab程序。（相关参数：密度7860KG/M3，长度L=1m，截面尺寸0.02m*0.02m）以下为matlab源代码：主程序clear;clc;Beam_InputData541;%输入相关参数，划分了40个单%元k=zeros(No_nel*No_dof,
3、No_nel*No_dof);%初始化单元刚度阵和质量阵m=zeros(No_nel*No_dof,No_nel*No_dof);kk=zeros(Sys_dof,Sys_dof);mm=zeros(Sys_dof,Sys_dof);%初始化总体质量阵和刚度阵index=zeros(No_nel*No_dof,1);foriel=1:No_elnd(1)=iel;nd(2)=iel+1;leng=0.025;[k,m]=BeamElement11(prop,leng);%调用子函数形成单元刚度矩阵%质量矩阵index=femEldof(
4、nd,No_nel,No_dof);kk=femAssemble1(kk,k,index);mm=femAssemble1(mm,m,index);%调用子函数形成总体刚度矩阵end%和质量矩阵Mzhiliangkuai=1/5*prop(3)*prop(6)*1;%将集中质量块加进总体质量矩Ms=zeros(82,82);%之中Ms(41,41)=Mzhiliangkuai;%此为质量块为于梁中央时的处mm=Ms+mm;理方法kk(1,:)=[];%引入边界条件，划去相应的kk(80,:)=[];%行和列kk(:,1)=[];kk(:
5、,80)=[];mm(1,:)=[];mm(80,:)=[];mm(:,1)=[];mm(:,80)=[];[V,D]=eig(kk,mm);%求解特征值和特征向量%或使用雅克比迭代求解特征值与特征向量，以下为代码%function[aa,v]=jac(a)%此程序用jacobi方法求实对称矩阵的全部特征值和特征向量%输入x：nxn矩阵%----------------------------------------------------%%n=length(a);%aa=zeros(1,n);%u=zeros(n,n);%l=0;%
6、v是n阶单位矩阵%v=eye(n,n);%while(1)%l=l+1;%fm=0.0;%fori=1:n%forj=1:n%ifi==j%continue;%else%d=abs(a(i,j));%ifd>fm%fm=d;%p=i;%q=j;%end%end%end%end%iffm<0.0001%break;%end%x=2*a(p,q)*(sign(a(p,p)-a(q,q)));%y=abs(a(p,p)-a(q,q));%ify~=0%cn=sqrt((1+y/sqrt(x^2+y^2))/2.0);%else%cn=1/sq
7、rt(2);%end%ify~=0%sn=x/sqrt(x^2+y^2)/(2*cn);%elseify==0&a(p,q)>0%sn=1/sqrt(2);%elseify==0&a(p,q)<0%sn=-1/sqrt(2);%end%%更新a%a_t=a;%fori=1:n%ifi~=p&i~=q%a(i,p)=a_t(i,p)*cn+a_t(i,q)*sn;%a(i,q)=-a_t(i,p)*sn+a_t(i,q)*cn;%a(p,i)=a(i,p);%a(q,i)=a(i,q);%end%end%a(p,p)=a_t(p,p)*c
8、n*cn+2*a_t(p,q)*cn*sn+a_t(q,q)*sn*sn;%a(q,q)=a_t(p,p)*sn*sn-2*a_t(p,q)*sn*cn+a_t(q,q)*cn*cn;%a(p,q)=0.0
