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

有限元课程设计 姓名：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,