薄板单元例题matlab133209140.pdfVIP

薄板问题有限元分析例题 例题：正方形薄板（如图 1 所示）边长 2l=2m ，厚度t=0.01m ，材料弹性模量E=200GPa ， 泊松比 μ=0.3 。薄板四边固支，其中心受垂直集中载荷 P=400N 作用，试求薄板中心挠 2 P (2l ) Et 3 度 w 。【理论解 D ）】 w =−0.0056 =−0m （ D 12 1− μ 2 ( ) 图1 薄板结构图 解：求解中采用国际单位制。 完整Matlab 程序 1. 单元划分： % Step 1 划分网格 Enode=[1 1 2 5 4; 2 2 3 6 5; 3 5 6 9 8; 4 4 5 8 7]; % 每一行表示一个单元的信息，第1 个数字表示单元编号，第2-5 个数字 表示该单元的结点编号。 ndof=3; % 薄板单元每个结点的自由度是3 。 Edof=caldof(Enode,ndof); % 计算各单元结点所对应的自由度 编号。 单元1 （结点1、2 、5、4 ） 单元2 （结点2 、3、6、5 ） ex=[0 1 1 0; 单元3 （结点5、6、9、8 ） 1 2 2 1; 单元4 （结点4 、5、8、7 ） 第 1 页 共 16 页 1 2 2 1; 0 1 1 0]; % 每一行表示一个单元4 个结点的x 坐标值。 ey=[0 0 1 1; 0 0 1 1; 1 1 2 2; 1 1 2 2]; % 每一行表示一个单元4 个结点的y 坐标值。 % Step 2 求单元刚度矩阵 2. 单元刚度矩阵分析： E=2.0e11; v=0.3; D=hooke(E,v); % E 是弹性模量，v 是泊松比，D 是 弹性矩阵。 t=0.01; ep=[t]; % ep=[t]表示薄板的厚度。 nie=size(Enode,1); ’ ’ ’ % nie 表示单元总数。 薄板单元在局部坐标系 X O Y 下的刚度矩阵按照以下公式计 算 for i=1:nie 3 e t T Ke(:,:,i)=platre(ex(i,:),ey(i,:),ep,D); K B D B dxdy [ ] 12 ∫∫ [ ] [ ][ ] End A % Ke(:,:,i)表示第 i 个单元的单元刚 式中 B 和 D 的计算公式在附录1 中。