固支梁有限元课程设计71p.docVIP

下载本文档

1
0
约7.84万字
约 71页
2018-05-08 发布于河北
举报
版权申诉

固支梁有限元课程设计71p.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

固支梁有限元课程设计71p

题目如图所示固支梁，高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,d