有限元编程的c++实现算例.docxVIP

下载本文档

14
0
约1.37万字
约 13页
2017-08-29 发布于河南
举报
版权申诉

有限元编程的c++实现算例.docx

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

有限元编程的c实现算例

有限元编程的c++实现算例/downloads76/doc/fileformat/290377/ganjian.cpp__.htm?? 1. #includestdio.h? ??? 2. #includemath.h? ??? 3.?? ??? 4.?? ??? 5. #define ne 3?????????????????????????????????????????????????? //单元数? ??? 6. #define nj 4?????????????????????????????????????????????????? //节点数? ??? 7. #define nz 6?????????????????????????????????????????????????? //支撑数? ??? 8. #define npj 0????????????????????????????????????????????????? //节点载荷数? ??? 9. #define npf 1????????????????????????????????????????????????? //非节点载荷数? ?? 10. #define nj3 12???????????????????????????????????????????????? //节点位移总数? ?? 11. #define dd 6?????????????????????????????????????????????????? //半带宽? ?? 12. #define e0 2.1E8????????????????????????????????????????????? //弹性模量? ?? 13. #define a0 0.008?????????????????????????????????????????????? //截面积? ?? 14. #define i0 1.22E-4???????????????????????????????????????????? //单元惯性距?? ?? 15. #define pi 3.141592654????????????????????????????????? ?? 16.?? ?? 17.?? ?? 18. int jm[ne+1][3]={{0,0,0},{0,1,2},{0,2,3},{0,4,3}};???????????? /*gghjghg*/?????????????? ?? 19. double gc[ne+1]={0.0,1.0,2.0,1.0};?????????????????????????????? ?? 20. double gj[ne+1]={0.0,90.0,0.0,90.0};? ?? 21. double mj[ne+1]={0.0,a0,a0,a0};? ?? 22. double gx[ne+1]={0.0,i0,i0,i0};? ?? 23. int zc[nz+1]={0,1,2,3,10,11,12};? ?? 24. double pj[npj+1][3]={{0.0,0.0,0.0}};? ?? 25. double pf[npf+1][5]={{0,0,0,0,0},{0,-20,1.0,2.0,2.0}};? ?? 26. double kz[nj3+1][dd+1],p[nj3+1];? ?? 27. double pe[7],f[7],f0[7],t[7][7];? ?? 28. double ke[7][7],kd[7][7];? ?? 29.?? ?? 30.?? ?? 31. //**kz[][]—整体刚度矩阵? ?? 32. //**ke[][]—整体坐标下的单元刚度矩阵? ?? 33. //**kd[][]—局部坐标下的单位刚度矩阵? ?? 34. //**t[][]—坐标变换矩阵? ?? 35.?? ?? 36. //**这是函数声明? ?? 37. void jdugd(int);? ?? 38. void zb(int);? ?? 39. void gdnl(int);? ?? 40. void dugd(int);? ?? 41.?? ?? 42.?? ?? 43. //**主程序开始? ?? 44. void main()? ?? 45. {? ?? 46.? int i,j,k,e,dh,h,ii,jj,hz,al,bl,m,l,dl,zl,z,j0;? ?? 47.? double cl,wy[7];? ?? 48.? int im,in,jn;? ?? 49.?? ?? 50. //*****