有限元课程设计.docVIP

PAGE 1 《有限元法》课程设计报告 班 级： 10级工程力学2班 姓 名： 刘海波 学 号： 1002060221 指导教师： 谢献忠 湖南科技大学土木工程学院 2013年 6月 12日  （一）问题描述 题2：图示薄板左边固定，右边受均布压力P=100Kn/m作用，板厚度为0.3cm； 三节点常应变单元；（2个和200个单元） 四节点矩形单元；（1个和50个单元） （二）有限元计算模型 因为该题荷载和边界都关于横向对称，故从中间取一半进行分析。网格划分方案如下图所示，将作用在单元上的均布载荷按虚功等效的原则移置到节点上，成为等效节点载荷。因取一半，故下方需在节点位置加上竖向链杆支座。  （三）基本理论 1.三角形单元的线性位移模式 单元应变矩阵 对于平面应力问题，弹性矩阵[D]为 总刚方程 [K]{d}={R} 2.矩形单元 位移模式 单元刚度矩阵 如果单元厚度t是常量，则 （四）计算结果及分析 两单元的计算结果及200单元的计算结果见Matlab中 每个单元的应变 2单元的应力δxδy Τxy sigm = 1.0e+009 * 0.0605 -0.3760 0.2419 1.1627 -0.4013 -0.2759 100单元的位移分析：（放大1000倍后） 结果数据过多，见matlab 矩形单元也进行了分析 放大1000倍的变形位移ans = -0.0000 0.0000 -0.0511 0.0112 -0.0000 0.0000 -0.0468 -0.0000 -0.0000 -0.0000 -0.0511 -0.0112 (五)多方案计算比较 取一半分析时，是三角形单元，2单元的没有100个单元的精度高,网格划分越多，结果也越精确，矩形双线性单元比三角形单元具有更高的精度。从定性分析来看，编程结果也基本符合实际物体受力。 （六）建议与体会 经过这些天的编程，我学到了很多关于matlab的技巧，对有限元也有了进一步的理解，当把所学的用程序表达出来时，虽然在所学过程中遇到了困难，花费了不少时间，当看到最后的编程成果时， 内心还是充满了成就感的，这些天的编程使我对三角形单元有了熟练的掌握，但对矩形双线性单元还有些困惑，深感自己在有限元课程上还有很长的道路要走。感谢老师的悉心教诲。回顾这个学期所学这门课程，有限元是门很有趣，很实用的课程。 （七）源程序 （1） %课程设计:第一问:2个单元的三角形单元，这里取一半来分析,然后对称 %生成原始数据 clc clear Gao=1.0;%板高1.0m Kuan=1.5;%板宽1.5m DG=1.0;% 沿高单元网格尺寸 DK=1.5;% 沿宽单元网格尺寸 NG=Gao/DG+1;% 沿高节点数(2) NK=Kuan/DK+1;% 沿寛节点数(2) t=0.003; %板厚3mm P=-100;%均布荷载100KN/m（100N/mm） E=1e9; v=0.25; for i=1:NG %NG=2 for j=1:NK %NK=2 x((i-1)*NK+j,:)=[(j-1)*DK Gao-(i-1)*DG];% 节点坐标 end end %节点整体编码与局部编码的对应关系 id(1,:)=[1 3 2]; id(2,:)=[2 3 4]; nelm=(NG-1)*(NK-1)*2;%单元数(2) nd1=NG*NK;%节点数(4) %计算单元的刚度矩阵并组集总刚 K=zeros(2*nd1,2*nd1);%整体刚度矩阵赋零初值 for i=1:nelm %单元常数(2) mA=[1 x(id(i,1),1) x(id(i,1),2) 1 x(id(i,2),1) x(id(i,2),2) 1 x(id(i,3),1) x(id(i,3),2)];%单元节点坐标矩阵 A(i)=0.5*det(mA);%单元面积 mb=[1 x(id(i,2),2) % bi 1 x(id(i,3),2)]; b(i,1)=-det(mb);