有限元分析课程设计.docx

文档介绍：有限元分析
课程设计

学院:土木建筑工程学院
专业: 工程力学
班级: 力学131
学号: 3121631023
姓名: 崔晨露


单元划分与节点编号

C语言程序
#include<>
#include<>
#define NE 32 //单元数
#define NJ 45 //节点数
#define NZ 10 //支承数
#define NPJ 7 //节点荷载数
#define NJ2 90 //节点位移数
#define DD 14 //半带宽
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,,},{0,,},{0,,3},{0,,0},{0,,},{0,,},{0,,},{0,,3},{0,,0},{0,,},{0,,},{0,,},{0,,3},{0,15,0},{0,15,},{0,15,},{0,15,},{0,15,3}};//共36个节点
int JM [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,12,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个单元
int NZC [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,-,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];
int IE,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,CETA;
double WY[8+