固支梁有限元课程设计毕业设计.docVIP

题目 如图所示固支梁，高3米，长15m，承受均布载荷q=10kN/，E=20GPa，u=0.167，厚度t=1m，忽略自重，按平面应力问题分析。用有限元方法计算梁的变形及应力分布，要求用矩形单元。 要求： 1.单元数目不得少于20个。 2.采用矩形单元计算求解。 3.计算结果并给出变形图、应力分布图、单元划分图。 二、力学分析 1.题目可以看做是平面应力问题故有LXM=0 2.单元划分图 三．程序框图及程序 程序框图： 四．源程序 #includestdio.h #includemath.h #define NJ 44 //节点数 #define NZ 16 //支承数 #define NPJ 11 //节点载荷作用数 #define DD 26 //半带宽 #define NJ2 88 //节点位移数 int LXM=0; double E0=2e7; double MU=0.167; double LOU=0.0; double TE=1 ; double AJZ[NJ+1][3]={{0,0,0}, {0,0,3},{0,1.5,3},{0,3,3},{0,4.5,3},{0,6,3},{0,7.5,3},{0,9,3},{0,10.5,3},{0,12,3},{0,13.5,3},{0,15,