固支梁有限元课程研究与设计开发.docVIP

q=10kN/ q=10kN/ 要求： 1.单元数目不得少于20个。 2.采用矩形单元计算求解。 3.计算结果并给出变形图、应力分布图、单元划分图。 二、力学分析 1.题目可以看做是平面应力问题故有LXM=0 2.单元划分图 三．程序框图及程序 程序框图： 开始插入基本参数插入其他参数形成整体刚度矩阵 开始 插入基本参数 插入其他参数 形成整体刚度矩阵 形成荷载列阵 引入支承条件 解方程输出位移 求应力输出应力 结束 子程序单元刚度矩阵 （ASK=3）计算单元刚度矩阵 （ASK=1）计算单元面积 （ASK=2）计算S矩阵 四．源程序 #includestdio.h #includemath.h #define NE 30 //单元数 #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,3}, {0,0,2},{0,1.5,2},{0,3,2},{0,4.5,2},{0,6,2},{0,7.5,2},{0,9,2},{0,10.5,2},{0,12,2},{0,13.5,2},{0,15,2}, {0,0,1},{0,1.5,1},{0,3,1},{0,4.5,1},{0,6,1},{0,7.5,1},{0,9,1},{0,10.5,1},{0,12,1},{0,13.5,1},{0,15,1}, {0,0,0},{0,1.5,0},{0,3,0},{0,4.5,0},{0,6,0},{0,7.5,0},{0,9,0},{0,10.5,0},{0,12,0},{0,13.5,0},{0,15,0}}; //节点坐标 int JM[NE+1][5]={{0,0,0,0,0},{0,10,11,22,21},{0,9,10,21,20},{0,8,9,20,19},{0,7,8,19,18},{0,6,7,18,17},{0,5,6,17,16},{0,4,5,16,15},{0,3,4,15,14},{0,2,3,14,13},{0,1,2,13,12},{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,32,33,44,43},{0,31,32,43,42},{0,30,31,42,41},{0,29,30,41,40},{0,28,29,40,39},{0,27,28,39,38},{0,26,27,38,37},{0,25,26,37,36},{0,24,25,36,35},{0,23,24,35,34}}; int NZC[NZ+1]={0,1,2,21,22,23,24,43,44,45,46,65,66,67,68,87,88 }; //支撑数组 double PJ[NPJ+1][2+1]={{0,0,0},{0,-7500,68},{0,-15000,70},{0,-15000,72},{0,-15000,74},{0,-15000,76},{0,-15000,78},{0,-15000,80},{0,-15000,82},{0,-15000,84},{0,-15000,86},{0,-7500,88}}; //节点载荷数组 //节点载荷数组 double AE,KZ[NJ2+1][DD+1],P[NJ2+1],S[3+1][8+1],KE[8+1][8+1],SZ[3+1][5*8+1]; double JDYL[NJ][6]; //节点应力矩阵 double DYYL[4][NE][4]; //单元应力矩阵 int IE,JE,ME,PE; void DUGD(int,int); FILE