原创力文档- 1、本文档共71页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
固支梁有限元课程设计毕业设计
题目
如图所示固支梁，高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,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 *fp1,*fp2,*ab;
void main()
{
int NJ1,k,IN,IM,jn,m,i,j,z,JO,ii,jj,h,dh,E,l,zl,
文档评论（0）