原创力文档- 1、本文档共7页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
【MATLAB 算例】基于 3节点三角形单元的 矩形薄板分析
将此结构按三角形单元 划分成 432个三角形 （X 方向分成 18段， Y 方向分成 12
段），总共分成 19X13=247个结点 的有限元模型，具体步骤详细程序如下：
tic;
Initial_info=[0.09 0.06 18 12];
disp(该程序计算的是 ,num2str(Initial_info(3)+1),X,num2str(Initial_info(4)+1),=,...
num2str((Initial_info(3)+1)*(Initial_info(4)+1)), 个结点的有限元模型 ]);
LX=Initial_info(1);
LY=Initial_info(2);
nx=Initial_info(3);
ny=Initial_info(4);
ne=2*nx*ny;
np=(nx+1)*(ny+1);
for i=1:nx+1;
j=1:ny+1;
Np(i,j)=j+(i-1)*(ny+1);
end
生成节点编号矩阵 Np
for i=1:nx+1;
j=1:ny+1;
XX(i,j)=(i-1)*LX/nx;
YY(i,j)=(j-1)*LY/ny;
end
XY=[reshape(XX,np,1),reshape(YY,np,1)];
nx2=nx/2;
Np1=Np(1:nx2+1,:);
Np2=Np(nx2+1:end,:);
for i=1:nx2*ny;
if rem(i,nx2)==0
xp=nx2;
yp=i/nx2;
else
xp=rem(i,nx2);
yp=fix(i/nx2)+1;
end
Dof1(i,:)=[Np1(xp,yp),Np1(xp+1,yp),Np1(xp,yp+1)];
Dof1(i+nx2*ny,:)=[Np1(xp+1,yp),Np1(xp+1,yp+1),Np1(xp,yp+1)];
Dof2(i,:)=[Np2(xp,yp),Np2(xp+1,yp),Np2(xp+1,yp+1)];
Dof2(i+nx2*ny,:)=[Np2(xp,yp),Np2(xp+1,yp+1),Np2(xp,yp+1)];
end
Dof=[Dof1;Dof2];
for i=1:ne
unit(i,:)=[XY(Dof(i,1),1),XY(Dof(i,2),1),XY(Dof(i,3),1),...
XY(Dof(i,1),2),XY(Dof(i,2),2),XY(Dof(i,3),2)];
end
disp(前处理完成 );
前处理完成
单元刚度矩阵
E=2*10^11;u=0.3; 平面应力问题
D=E/(1-u^2)*[1 u 0;u 1 0;0 0 (1-u)/2];
for i=1:ne
xi=unit(i,1);
yi=unit(i,4);
xj=unit(i,2);
yj=unit(i,5);
xm=unit(i,3);
ym=unit(i,6);
ai=xj*ym-xm*yj;
aj=xm*yi-xi*ym;
am=xi*yj-xj*yi;
bi=yj-ym;
bj=ym-yi;
bm=yi-yj;
ci=-(xj-xm);
cj=-(xm-xi);
cm=-(xi-xj);
area=abs((ai+aj+am)/2);
B = [bi 0 bj 0 bm 0
0 ci 0 cj 0 cm
ci bi cj bj cm bm];
Be{i,1} = B/2/area ;
ke{i,1}=[Be{i,1}]*D*Be{i,1}*area;
end
总刚度矩阵叠加
KK=sparse(2*n
文档评论（0）