原创力文档- 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语言
#include stdio.h
#include math.h
main()
{
FILE *fp;
int i,j,NN,NE,NF,n,m,z[20],a[20][2],g[20][4],c[20];
float E,A,Cx,Cy,l,T,Xi,Xj,Yi,Yj,K[20][20]={0},k[4][4]={0},F[4][4]={0},y,t,d,ui,uj,vi,vj,Ui,Uj,Vi,Vj,fxi,fxj,fyi,fyj;
clrscr();
fp=fopen(data3.txt,r);
fscanf(fp,%d,NN);
fscanf(fp,%d,NE);
fscanf(fp,%d,NF);
for(i=0;i=2*NF-1;i++)
fscanf(fp,%d,z[i]);
for(i=0;i=NN-1;i++)
for(j=0;j=1;j++)
fscanf(fp,%d,a[i][j]);
for(i=0;i=NN-1;i++)
for(j=0;j=3;j++)
fscanf(fp,%d,g[i][j]);
for(i=0;i=2*NN-1;i++)
fscanf(fp,%d,c[i]);
fclose(fp);
fp=fopen(answer3.txt,w);
fprintf(fp,NN\tNE\tNF\n);
fprintf(fp,%d\t%d\t%d\n,NN,NE,NF);
fprintf(fp,SUPPRESSED DISPLACIMENT POINT\n);
for(i=0;i=2*NF-1;i++)
{
fprintf(fp,%d\t,z[i]);
if(i%2!=0)
fprintf(fp,\n);
}
fprintf(fp,NO.\tX-COORD.\tY-COORD\n);
for(i=0;i=NN-1;i++)
{
fprintf(fp,%d\t,i+1);
for(j=0;j=1;j++)
fprintf(fp,%d\t ,a[i][j]);
fprintf(fp,\n);
}
fprintf(fp,I\tJ\tA\tE\n);
for(i=0;i=NN-1;i++)
{
for(j=0;j=3;j++)
fprintf(fp,%d\t,g[i][j]);
fprintf(fp,\n);
}
fprintf(fp,VECTOR OF NODAL LOAD\n);
fprintf(fp,NO.\tX-LOAD.\tY-LOAD.\n);
for(i=0;i=NN-1;i++)
{
fprintf(fp,(%d)\t,i+1);
fprintf(fp,%d\t%d,c[2*i],c[2*i+1]);
fprintf(fp,\n);
}
fprintf(fp,THE NODAL DISPLACEMENT ARE:\n);
for(i=1;i=NN;i++)
for(j=i+1;j=NN;j++)
for(t=0;t=3;t++)
{
if(i==g[t][0]j==g[t][1])
{
Xi=a[i-1][0];
Yi=a[i-1][1];
Xj=a[j-1][0];
Yj=a[j-1][1];
l=sqrt((Xi-Xj)*(Xi-Xj)+(Yi-Yj)*(Yi-Yj));
Cx=(Xj-Xi)/l;
Cy=(Yj-Yi)/l;
E=g[t][3];
A=g[t][2];
T=E*A/l;
k[0][0]=T*Cx*Cx;
k[0][1]=T*Cx*Cy;
k[0][2]=T*(-1)*Cx*Cx;
k[0][3]=T*(-1)*Cx*Cy;
k[1][0]=T*Cx*Cy;
k[1][1]=T*Cy*Cy;
k[1][2]=T*(-1)*Cx*Cy;
k[1][3]=T*(-1)*Cy*Cy;
k[2][0]=T*(-1)*Cx*Cx;
k[2][1]=T*(-1)*Cx*Cy;
k[2][2]=T*Cx*Cx;
k[2][3]=T*Cx*Cy;
k[3][0]=T*(-1)*Cx*Cy;
k[3][1]=T*(-1)*Cy*Cy;
k[3][2]=T*Cx*Cy;
k[3][3]=T*Cy*Cy;
K[2*i-1][2*i-1]=K[2*i-1][2*i-1]+k[0][0];
K[2*i-1][2*i]=K[2*i-1][2*i]+k[0][1];
K[2*i][2*i-1]=K[2*i][2*i-1]+k[1][0];
K[2*i][2*i]=K[2*i][2*i]+k[1][1];
K[2*j-1][2*j-1]=K[2*j-1][2*j-1]+k[2][2];
K[2*j-1][2*j]=K[2*j-1][2*j]+k[2][3];
K[2*j][2*j
文档评论（0）