悬臂梁受载荷的 矩形单位 无限元课程设计[整理版].docVIP

有限元课程设计姓名：XXX专业：工程力学081学号：XXXXXXXXX指导老师：XXX2011年6月29日一、题目已知：悬臂梁的材料参数为E=190GPa =0.25 t=10 cm =0 尺寸参数如图所示。要求：1.单元数目不得少于20个。2.采用矩形单元计算求解。3.计算结果并给出变形图、应力分布图、单元划分图。二、力学分析1.题目可以看做是平面应力问题故有LXM=02.单元划分图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 26int 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,