悬臂梁受载荷的矩形单元有限元课程设计.doc

有限元课程设计 姓名：XXX 专业：工程力学081 学号：XXXXXXXXX 指导老师：XXX 2011年6月29日 一、题目 已知：悬臂梁的材料参数为E=190GPa =0.25 t=10 cm =0 尺寸参数如图所示。 要求： 1.单元数目不得少于20个。 2.采用矩形单元计算求解。 3.计算结果并给出变形图、应力分布图、单元划分图。 二、力学分析 1.题目可以看做是平面应力问题故有LXM=0 2.单元划分图 3.简化成力学模型 根据力学模型可以写出支承数组NZC、节点载荷数组PJ 三、程序代码 #include stdio.h #include math.h #define NE 40 #define NJ 55 #define NZ 10 #define NPJ 14 #define NJ2 110 #define DD 26 int LXM=0; double E0=190000000000; double MU=0.25; double LOU=0.0; double TE=0.1; double AJZ[NJ+1][3]={{0,0,0}, {0,0,0},{0,0.5,0},{0,1,0},{0,1.5,0},{0,2,0},{0,2.5,0},{0,3,0},{0,3.5,0},{0,4,0},{0,4.5,0},{0,5,0}, {0,0,0.25},{0,0.5,0.25},{0,1,0.25},{0,1.5,0.25},{0,2,0.25},{0,2.5,0.25},{0,3,0.25},{0,3.5,0.25},{0,4,0.25},{0,4.5,0.25},{0,5,0.25}, {0,0,0.5},{0,0.5,0.5},{0,1,0.5},{0,1.5,0.5},{0,2,0.5},{0,2.5,0.5},{0,3,0.5},{0,3.5,0.5},{0,4,0.5},{0,4.5,0.5},{0,5,0.5}, {0,0,0.75},{0,0.5,0.75},{0,1,0.75},{0,1.5,0.75},{0,2,0.75},{0,2.5,0.75},{0,3,0.75},{0,3.5,0.75},{0,4,0.75},{0,4.5,0.75},{0,5,0.75}, {0,0,1},{0,0.5,1},{0,1,1},{0,1.5,1},{0,2,1},{0,2.5,1},{0,3,1},{0,3.5,1},{0,4,1},{0,4.5,1},{0,5,1}}; int JM[NE+1][5]={{0,0,0,0,0}, {0,1,2,13,12},{0,2,3,14,13},{0,3,4,15,14},{0,4,5,16,15},{0,5,6,17,16}, {0,6,7,18,17},{0,7,8,19,18},{0,8,9,20,19},{0,9,10,21,20},{0,10,11,22,21}, {0,21,22,33,32},{0,20,21,32,31},{0,19,20,31,30},{0,18,19,30,29},{0,17,18,29,28}, {0,16,17,28,27},{0,15,16,27,26},{0,14,15,26,25},{0,13,14,25,24},{0,12,13,24,23}, {0,23,24,35,34},{0,24,25,36,35},{0,25,26,37,36},{0,26,27,38,37},{0,27,28,39,38}, {0,28,29,40,39},{0,29,30,41,40},{0,30,31,42,41},{0,31,32,43,42},{0,32,33,44,43}, {0,43,44,55,54},{0,42,43,54,53},{0,41,42,53,52},{0,40,41,52,51},{0,39,40,51,50}, {0,38,39,50,49},{0,37,38,49,48},{0,36,37,48,47},{0,35,36,47,46},{0,34,35,46,45}}; int NZC[NZ+1]={0,1,2,23,24,45,46,67,68,89,90}; double PJ[NPJ+1][2+1]={{0,0,0}, {0,-104166.7,21},{0,-125000,43},{0,125000,87},{0,250000,89},{0,500000,91}, {0,500000,93},{0,500000,95},{0,500000,97},{0,500000,99},{0,50