有限元程序设计作业.docVIP

有限元程序设计作业 2.1 平面刚架结构的程序设计 1、程序源代码 注明程序每一部分的功能 2、计算对象，单元及节点编码在图上标注，已知条件，需要求解的内容 计算题目可以自行选择 3、所选桁架算例如下： 确定节点1的水平、竖直位移，以及每一个杆的应力。 弹性模量E=200Gpa,截面面积为A=6×10-4m2。 4、总刚度矩阵以及程序 #include stdio.h #include math.h #define N 2 #define M 2 #define X 4 #define P 3 #define Q 2 /*N:结点数 M:自由度 X：总刚节点数 P:杆数量 Q:杆：2，梁：3*/ int MATNO1[P]={1,1,1};//定义 float PROPS[P][2]={{210,400},{210,400},{210,400}}; float IA1[P][N]={{1,2},{1,3},{1,4}}; float XY[X][2] ={{0,0},{3,0},{44,{0,3}}; float ekzong[M*X][M*X]={0}; float ek1[M*2][M*2]={0}; float ek2[M*2][M*2]={0}; float ek3[M*2][M*2]={0}; void main() {