有限元分析简单实例之平面矩形薄板（matlab）

最新推荐文章于 2021-10-30 17:22:45 发布

monarch-_-

最新推荐文章于 2021-10-30 17:22:45 发布

阅读量9.9k

点赞数 16

本文链接：https://blog.csdn.net/w1029126938/article/details/114676929

版权

本文介绍了使用matlab进行平面矩形薄板的有限元分析，包括问题描述、所需函数、解答步骤。通过计算单元刚度矩阵、组装刚度矩阵、处理边界条件来求解节点位移和支座反力。

摘要由CSDN通过智能技术生成

有限元分析简单实例之平面矩形薄板（matlab）

问题描述

在这里插入图片描述
对于如图所示的一个平面矩形薄板结构，施加如右图所示的几个方向力，对其进行有限元分析，计算各个节点的位移及支座反力。（其中F是合力，E是弹性模量，μ是泊松比，t是厚度）

要用到的函数

（1）计算单元的刚度矩阵

function k = Triangle2D3Node_Stiffness(E,NU,t,xi,yi,xj,yj,xm,ym,ID)
% 	计算单元的刚度矩阵
%   输入弹性模量E，泊松比NU，厚度t
%   输入三个节点的坐标xi,yi,xj,yj,xm,ym
%  	输入平面问题性质指标参数ID（1为平面应力，2为平面应变）
%  	输出单元刚度矩阵k(6*6)
A = (xi*(yj-ym)+xj*(ym-yi)+xm*(yi-yj))/2;
betai = yj-ym;
betaj = ym-yi;
betam = yi-yj;
gammai = xm-xj;
gammaj = xi-xm;
gammam = xj-xi;
B = [betai 0 betaj 0 betam 0;
    0 gammai 0 gammaj 0 gammam;
    gammai betai gammaj betaj gammam betam]/(2*A);
if ID==1
    D=(E/(1-NU*NU))*[1 NU 0;NU 1 0;0 0 (1-NU)/2];
elseif ID==2
    D =(E/(1+NU)/(1-2*NU))*[1-NU NU 0;NU 1-NU 0;0 0 (1-2*NU)/2];
end
k = t*A*B'*D*B;
end