梁的有限元C++程序设计.docVIP

计算机编写程序实验课及实践课程 指导教师：孙先念（教授） 目标掌握利用c++语言编写有限元程序。 例如：计算某货轮第三货舱肋骨刚架（如下图所示：） 程序代码：（肖有才创作） #include iostream #include iomanip #include cmath #include fstream using namespace std; int main() {int i,j,ng,nj,WZSL[100][2],n=1,jad[7]; double tk[20]; float I[100],XY[50][2],KE[7][7],KK[50][50],AREA[100],zancun[100],E,PE[7][4],PP[300],pe[7],CJ[20][4],x[15],PJ[5][4]; ifstream infile(jiegoulixuein.txt,ios::in|ios::binary); ofstream outfile(datagangjiaout.txt,ios::in|ios::binary|ios::out); for(i=1;i=100;i++) infilezancun[i]; ng=zancun[1]; nj=zancun[2]; E=zancun[3]; coutng=ng nj=nj E= Eendl; for(i=1;i=nj*2;i++) for(j=1;j=2;j++) { XY[i][j]=zancun[n+3]; n++;} for(j=1;j=nj;j++) for(i=1;i=2;i++) { if(i==1) coutendl; coutsetw(4)setiosflags(ios::left)setiosflags(ios::fixed)setprecision(1)XY[j][i]; } n=1; for(i=1;i=nj*2;i++) for(j=1;j=2;j++) { WZSL[i][j]=zancun[n+3+2*nj]; n++;} for(j=1;j=ng;j++) for(i=1;i=2;i++) { if(i==1) coutendl; coutsetw(4)setiosflags(ios::left)setiosflags(ios::fixed)setprecision(0)WZSL[j][i]; } n=1; for(i=1;i=ng;i++) {I[i]=zancun[n+3+2*nj+2*ng]; n++;} for(j=1;j=ng;j++) { coutsetw(10)setiosflags(ios::left)setiosflags(ios::fixed)setprecision(8)setiosflags(ios::scientific)I[j]; } coutendl; n=1; for(i=1;i=ng;i++) {AREA[i]=zancun[n+3+2*nj+2*ng+ng]; n++;} for(j=1;j=ng;j++) { coutendl; coutsetw(10)setiosflags(ios::left)setiosflags(ios::fixed)setprecision(8)setiosflags(ios::scientific)AREA[j]; } coutendl; n=1; for(i=1;i=3;i++) for(j=1;j=4;j++) {PE[i][j]=zancun[n+3+2*nj+2*ng+ng+ng]; n++;} for(i=1;i=3;i++) for(j=1;j=4;j++) { if(j==1) coutendl; coutsetw(8)setiosflags(ios::left)setiosflags(ios::fixed)setprecision(2)PE[i][j]; } coutendl; n=1; for(i=1;i=3;i++) for(j=1;j=4;j++) {PJ[i][j]=zancun[n+3+2*nj+2*ng+ng+