有限元分析课程设计.docx

文档介绍：有限元分析课程设计学院: 土木建筑工程学院专业: 工程力学班级: 力学 131 学号: 3121631023 姓名: 崔晨露一、单元划分与节点编号二、 C语言程序#include<> #include<> #define NE32//单元数#define NJ45//节点数#define NZ10//支承数#define NPJ 7//节点荷载数#define NJ2 90//节点位移数#define DD14//半带宽 int LXM=0; //类型判别码 double EO=170e9; //杨氏模量 double MU=; //泊松比 double LOU=0; //容重 double TE=; //厚度 double AJZ [NJ+1][3]= {{0,0,0},{0,0,0},{0,0,},{0,0,},{0,0,},{0,0,3},{0,,0}, {0,,},{0,,},{0,,},{0,,3},{0,,0},{0, ,},{0,,},{0,,},{0,,3},{0,,0},{0,, },{0,,},{0,,},{0,,3},{0,,0},{0,,}, {0,,},{0,,},{0,,3},{0,,0},{0,,},{0, 5,},{0,,},{0,,3},{0,,0},{0,,},{0,, },{0,,},{0,,3},{0,,0},{0,,},{0, 5,},{0,,},{0,,3},{0,15,0},{0,15,},{0,15,},{ 0,15,},{0,15,3}};// 共36个节点 intJM[NE+1][5]= {{0,0,0,0},{0,1,6,7,2},{0,2,7,8,3},{0,3,8,9,4},{0,4,9,10,5},{0,6,11,1 2,7},{0,7,12,13,8},{0,8,13,14,9},{0,9,14,15,10},{0,11,16,17,12},{0,12, 17,18,13},{0,13,18,19,14},{0,14,19,20,15},{0,16,21,22,17},{0,17,22,23, 18},{0,18,23,24,19},{0,19,24,25,20},{0,21,26,27,22},{0,22,27,28,23},{ 0,23,28,29,24},{0,24,29,30,25},{0,26,31,32,27},{0,27,32,33,28},{0,28, 33,34,29},{0,29,34,35,30},{0,31,36,37,32},{0,32,37,38,33},{0,33,38,39, 34},{0,34,39,40,35},{0,36,41,42,37},{0,37,42,43,38},{0,38,43,44,39},{ 0,39,44,45,40}};// 共32个单元 intNZC [NZ+1]={0,1,2,3,4,5,6,7,8,9,10};//1- 5 号节点的 x,y 被约束 double PJ[NPJ+1][2+1]={{0,0,0},{0,-,50},{0,- ,60},{0,-,70},{0,-,80},{0,8e4,83},{0,8e4,87},{0,- 4,90}}; double AE,JS[NJ*4][4],js[NJ+1][4],KZ[NJ2+1][DD+1] ,P[NJ2+1],S[3+1][8+1],KE[8+1][8+1],SZ[3+1][32+1]; intIE,JE,ME,LE; void DUGD(int,int);// 生成 S矩阵, KE矩阵 void main() {int NJ1,k,IN,IM,jn,m,i,j,z,JO,ii,jj,h,dh,E,l, zl,dl,n; double PE,c,SIG1,SIG2,SIG3,PYL,RYL,MAYL,MIYL,CET A;double WY[8+1],YL[3+1]; if(LXM!=0)// 平面应力问题与平面应变问题的判别{EO=EO/(