有限元单元刚度矩阵汇.docVIP

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

实现功能:输入欲求的单元号,得到该单元的面积和该单元刚度矩阵;输入完所有的单元号,得到所有的单元刚度矩阵后,直接得到半带宽存储的数组。(所有的单元刚度矩阵和半带宽数组分别存在所有单元的刚度矩阵.txt和SK矩阵.txt) 第1单元面积为: 0.5000 弹性模量、泊松比和厚度分别为:100.0000 0.3000 0.1000 单元1单元的应力矩阵 -109.8901 -32.9670 109.8901 0.0000 0.0000 32.9670 -32.9670 -109.8901 32.9670 0.0000 0.0000 109.8901 -38.4615 -38.4615 0.0000 38.4615 38.4615 0.0000 第1单元的单元刚度阵为 7.4176 3.5714 -5.4945 -1.9231 -1.9231 -1.6484 3.5714 7.4176 -1.6484 -1.9231 -1.9231 -5.4945 -5.4945 -1.6484 5.4945 0.0000 0.0000 1.6484 -1.9231 -1.9231 0.0000 1.9231 1.9231 0.0000 -1.9231 -1.9231 0.0000 1.9231 1.9231 0.0000 -1.6484 -5.4945 1.6484 0.0000 0.0000 5.4945 半带宽存储的数组 7.4176 3.5714 -1.9231 -1.6484 -5.4945 -1.9231 7.4176 -1.9231 -5.4945 -1.6484 -1.9231 0.0000 9.8901 -2.3810 -4.9451 4.7619 -3.0220 -0.4579 9.8901 4.7619 -4.9451 -0.7326 0.5495 0.0000 20.8791 -2.3810 0.5495 -2.3810 -10.989 1.6484 13.7363 -2.3810 -3.0220 1.9231 -3.8462 0.0000 3.4341 1.1905 -0.9615 1.9231 0.0000 0.0000 5.2198 1.6484 -2.7473 0.0000 0.0000 0.0000 11.9506 -3.5714 0.0000 0.0000 0.0000 0.0000 6.5934 0.0000 0.0000 0.0000 0.0000 0.0000 #include stdio.h #define jds 5 #define d 2 main() {int nj,ne,i,j,k,a,n,*lnd; int JDS[3];xj,yj,xo,yo; float B[3][6],b1[3][6];D[3][3],E,u,t;S[3][6],s1[6][3];K[6][6];Kk[2*jds][2*jds]={{0}}; float SK[2*jds][(d+1)*2]={{0}}; float x1,x2,x3,y1,y2,y3,ae,*xy; FILE *fp1,*fp2,*fp3,*fp4,*fp5,*fp7; fp1=fopen(节点坐标.txt,r); fp2=fopen(单元节点编号.txt,r);fp3=fopen(材料参数.txt,r); fp4=fopen(输出结果.txt,a); fp5=fopen(节点总数和单元总数.txt,r);fp7=fopen(SK矩阵.txt,w); fscanf(fp5,%d,nj); fscanf(fp5,%d,ne); xy=(float*)malloc(nj*2*sizeof(float)); lnd=(int*)malloc(ne*3*sizeof(int)); for(i=0;inj;i++) for(j=0;j2;j++) fscanf(fp1,%f,xy[i*2+j]); for (i=0;ine;i++) for (j=0;j3;j++) fscanf(fp2,%d,lnd[i*3+j]); for(n=0;nne;n++){ printf(请输入所求的单元号:); scanf(%d,a); i=lnd[(a-1)*3+0],j=lnd[(a-1)*3+1],k=lnd[(a-1)*3+2]; JDS[0]=i;JDS[1]=j;JDS[2]=k

文档评论(0)

hhuiws1482 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

版权声明书
用户编号:5024214302000003

1亿VIP精品文档

相关文档