E=2.1e11; %elastic molulus
poisson =0.3; % poisson ratio
density=7.3e3; %density
t=0.05; %plate thickness
lx=2; %length in x direction
ly=2; %length in y direction
jdx=11; %number of nodes in x direction
jdy=11; %number of nodes in y direction
k(1:330,1:330)=0; %system stiffness matrix
m(1:330,1:330)=0; %system mass matrix
%prepare the arrays which are needed to describe this problem
en(1:100,1:4)=0; %element node
for ni=1:jdx-1
for nj=1:jdy-1
en(ni+(nj-1)*(jdx-1),1)=ni+(nj-1)*jdx;
en(ni+(nj-1)*(jdx-1),2)=ni+1+(nj-1)*jdx;
en(ni+(nj-1)*(jdx-1),4)=ni+nj*jdx;
en(ni+(nj-1)*(jdx-1),3)=ni+1+nj*jdx;
end
end
disp(1:jdx*jdy,1:3)=1; % node displacement
constraints=1:jdx:jdx*jdy; % constraints
disp(constraints,:)=0;
dof=0; %degree of freedom
for ni=1:jdx*jdy
for nj=1:3
if disp(ni,nj)~=0
dof=dof+1;
disp(ni,nj)=dof;
end
end
end
el=lx/(jdx-1); %element length
eh=ly/(jdy-1); %element height
[ek,dm]=km(el/2,eh/2,mu,poisson,E,density);
%km: function used to compute element stifness and mass matrix
%in this case, all elements have the same element stifness and mass matrix.
%built system stifness and mass matrix.
index(1:12)=0; % vector sontaining system dofs of nodes in each element.
for loopi=1:(jdx-1)*(jdy-1)
for zi=1:4
index((zi-1)*3+1)=disp(en(loopi,zi),1);
index((zi-1)*3+2)=disp(en(loopi,zi),2);
index((zi-1)*3+3)=disp(en(loopi,zi),3);
end
for jx=1:12
for jy=1:12
if(index(jx)*index(jy)~=0)
k(index(jx),index(jy))=k(index(jx),index(jy))+ek(jx,jy);
m(index(jx),index(jy))=m(index(jx),index(jy))+em(jx,jy);
end
end
end
end
%solve eigenvalue problem
[v,d] = eig(k,m);
tempd=diag(d);
[nd,sortindex]=sort(tempd);
v=v(:,sortindex);
mode_number=1:15;
frequency(mode_number)=sqrt(nd(mode_number))/(2*pi);
function [k,m]=km(a,b,poisson,t,E,density)
k=[E*t^3/(360-360*poisson^2)/a/b*(21-6*poisson+30*b^2/a^2+30*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(3*b+12*poisson*b+30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-3*a-12*poisson*a-30*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(-21+6*poisson-30*b^2/a^2+15*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*b-12*poisson*b+15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-3*a+3*poisson*a-30*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(21-6*poisson-15*b^2/a^2-15*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*b+3*poisson*b+15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(3*a-3*poisson*a-15*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(-21+6*poisson+15*b^2/a^2-30*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(3*b-3*poisson*b+30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(3*a+12*poisson*a-15*b^2/a);
E*t^3/(360-360*poisson^2)/a/b*(3*b+12*poisson*b+30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(8*b^2-8*poisson*b^2+40*a^2), -30*E*t^3/(360-360*poisson^2)*poisson, E*t^3/(360-360*poisson^2)/a/b*(-3*b-12*poisson*b+15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-8*b^2+8*poisson*b^2+20*a^2), 0,
E*t^3/(360-360*poisson^2)/a/b*(3*b-3*poisson*b-15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(2*b^2-2*poisson*b^2+10*a^2), 0,
E*t^3/(360-360*poisson^2)/a/b*(-3*b+3*poisson*b-30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-2*b^2+2*poisson*b^2+20*a^2), 0;
E*t^3/(360-360*poisson^2)/a/b*(-3*a-12*poisson*a-30*b^2/a), -30*E*t^3/(360-360*poisson^2)*poisson, E*t^3/(360-360*poisson^2)/a/b*(8*a^2-8*poisson*a^2+40*b^2), E*t^3/(360-360*poisson^2)/a/b*(3*a-3*poisson*a+30*b^2/a), 0,
E*t^3/(360-360*poisson^2)/a/b*(-2*a^2+2*poisson*a^2+20*b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*a+3*poisson*a+15*b^2/a), 0,
E*t^3/(360-360*poisson^2)/a/b*(2*a^2-2*poisson*a^2+10*b^2), E*t^3/(360-360*poisson^2)/a/b*(3*a+12*poisson*a-15*b^2/a), 0,
E*t^3/(360-360*poisson^2)/a/b*(-8*a^2+8*poisson*a^2+20*b^2);
E*t^3/(360-360*poisson^2)/a/b*(-21+6*poisson-30*b^2/a^2+15*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*b-12*poisson*b+15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(3*a-3*poisson*a+30*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(21-6*poisson+30*b^2/a^2+30*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(3*b+12*poisson*b+30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(3*a+12*poisson*a+30*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(-21+6*poisson+15*b^2/a^2-30*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(3*b-3*poisson*b+30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-3*a-12*poisson*a+15*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(21-6*poisson-15*b^2/a^2-15*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*b+3*poisson*b+15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-3*a+3*poisson*a+15*b^2/a);
E*t^3/(360-360*poisson^2)/a/b*(-3*b-12*poisson*b+15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-8*b^2+8*poisson*b^2+20*a^2), 0,
E*t^3/(360-360*poisson^2)/a/b*(3*b+12*poisson*b+30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(8*b^2-8*poisson*b^2+40*a^2), 30*E*t^3/(360-360*poisson^2)*poisson, E*t^3/(360-360*poisson^2)/a/b*(-3*b+3*poisson*b-30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-2*b^2+2*poisson*b^2+20*a^2), 0,
E*t^3/(360-360*poisson^2)/a/b*(3*b-3*poisson*b-15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(2*b^2-2*poisson*b^2+10*a^2), 0;
E*t^3/(360-360*poisson^2)/a/b*(-3*a+3*poisson*a-30*b^2/a), 0,
E*t^3/(360-360*poisson^2)/a/b*(-2*a^2+2*poisson*a^2+20*b^2), E*t^3/(360-360*poisson^2)/a/b*(3*a+12*poisson*a+30*b^2/a), 30*E*t^3/(360-360*poisson^2)*poisson, E*t^3/(360-360*poisson^2)/a/b*(8*a^2-8*poisson*a^2+40*b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*a-12*poisson*a+15*b^2/a), 0,
E*t^3/(360-360*poisson^2)/a/b*(-8*a^2+8*poisson*a^2+20*b^2), E*t^3/(360-360*poisson^2)/a/b*(3*a-3*poisson*a-15*b^2/a), 0,
E*t^3/(360-360*poisson^2)/a/b*(2*a^2-2*poisson*a^2+10*b^2);
E*t^3/(360-360*poisson^2)/a/b*(21-6*poisson-15*b^2/a^2-15*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(3*b-3*poisson*b-15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-3*a+3*poisson*a+15*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(-21+6*poisson+15*b^2/a^2-30*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*b+3*poisson*b-30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(-3*a-12*poisson*a+15*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(21-6*poisson+30*b^2/a^2+30*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(-3*b-12*poisson*b-30*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(3*a+12*poisson*a+30*b^2/a), E*t^3/(360-360*poisson^2)/a/b*(-21+6*poisson-30*b^2/a^2+15*a^2/b^2), E*t^3/(360-360*poisson^2)/a/b*(3*b+12*poisson*b-15*a^2/b), E*t^3/(360-360*poisson^2)/a/b*(3*a-3*poisson*a+30*b^2/a);
E*t^3/(360-360*poisson^2)/a/b*(-3*b+3*poisson*b+15*a^2/b), E*t^
pro6.zip_matlab 壳_shell element_壳 有限元_板壳_板壳有限元
版权申诉
5星 · 超过95%的资源 156 浏览量
2022-07-13
19:39:19
上传
评论
收藏 4KB ZIP 举报
JaniceLu
- 粉丝: 85
- 资源: 1万+
最新资源
- labelImg安装指导书.docx
- 2023AI自有光-她经济消费新图鉴(2023)-百度营销.pdf
- 2022中国新能源汽车内容生态趋势洞察(1).pdf
- Docker技术:Docker安装与配置教程+运维技术+超融合+虚拟技术+云计算
- ZeRO Memory Optimizations Toward Training LLM.pdf
- 高效SQL语句编写(how-to-write-efficient-sql)
- ZeroTermux-release_sign.apk
- 非弹性斜碰撞物理课件模拟-HTML网页制作
- 软件测试学习日志-测试基础-day02
- 制作一个简单的进销存(库存管理)页面.rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
评论2