猴子YWY 猴子YWY
关注数: 134 粉丝数: 146 发帖数: 1,739 关注贴吧数: 26
球各种电脑高手解决一下这程序 跪求 。。。。顶 #include "stdio.h" #include <math.h> #include <stdlib.h> double areas (double,double, double,double,double,double); //在此定义函数原型,计算三角形单元面积 //double fabs(double x); main() { int i,j,k,r,s; const double T=1000; const double E=30000000; const double nou=0.3; const double t=1; const int dy_n=2; const int jd_n=4; const int jd_n2=2*jd_n; const int dy_jdn=3; double jd_x[100],jd_y[100];//jd_z[100]; double dyjd_x[dy_n+1][dy_jdn],dyjd_y[dy_n+1][dy_jdn];//dyjd_z[dy_n+1][dy_jdn]; double A[dy_n+1]; double K_e[jd_n2+1][jd_n2+1][dy_n+1], K_ea[jd_n2+1][jd_n2+1][dy_n+1], K_ea1[jd_n2+1][jd_n2+1][dy_n+1]; double K_all[jd_n2+1][jd_n2+1]; double F_all[jd_n2+1]; double d_all[jd_n2+1]; //double d_all_0[jd_n2+1]; double d_dy[2*dy_jdn+1][dy_n+1]; double d_bj; double nn_a[dy_n+1][dy_jdn],nn_b[dy_n+1][dy_jdn],nn_c[dy_n+1][dy_jdn]; double DY_B[4][7][dy_n+1]; double DY_DE[4][4]; double DY_YL[4][dy_n+1]; double n_max=10000000000; int n_num=6; double n_min=0.001; double n_min1=0.00000001; char buf[40]=""; int ndigit=8; int ii,jj,kk; double ff1,ff2; int mm1; double ss[4],ss1[jd_n2+1][jd_n2+1]; jd_x[1]=0;jd_x[2]=0;jd_x[3]=20;jd_x[4]=20; jd_y[1]=0;jd_y[2]=10;jd_y[3]=10;jd_y[4]=0; F_all[5]=100*t*T/20;F_all[6]=0; F_all[7]=100*t*T/20;F_all[8]=0; for(j=1;j<=dy_n;j++) d_all[j]=10000; d_all[1]=d_all[2]=d_all[3]=d_all[4]=0; for (i=1;i<=dy_n;i++) {for (j=0;j<=jd_n2;j++) for(k=0;k<=jd_n2;k++) K_e[j][k][i]= K_ea[j][k][i]= K_ea1[j][k][i]=0; switch(i) { case 1: ii=1;jj=3;kk=2; K_e[1][0][i]=K_e[2][0][i]=ii;//用0行0列标记每行每列所属节点号 K_e[3][0][i]=K_e[4][0][i]=jj; K_e[5][0][i]=K_e[6][0][i]=kk; for(j=1;j<=2*dy_jdn;j++) K_e[0][j][i]=K_e[j][0][i]; break; case 2: ii=1;jj=4;kk=3; K_e[1][0][i]=K_e[2][0][i]=ii; K_e[3][0][i]=K_e[4][0][i]=jj; K_e[5][0][i]=K_e[6][0][i]=kk; for(j=1;j<=2*dy_jdn;j++) K_e[0][j][i]=K_e[j][0][i]; break; } for(j=1;j<=2*dy_jdn;j++) K_ea[j][0][i]=K_e[j][0][i]; for(k=1;k<=2*dy_jdn;k++) K_ea[0][k][i]=K_e[0][k][i]; dyjd_x[i][0]=jd_x[ii]; dyjd_y[i][0]=jd_y[ii]; dyjd_x[i][1]=jd_x[jj]; dyjd_y[i][1]=jd_y[jj]; dyjd_x[i][2]=jd_x[kk]; dyjd_y[i][2]=jd_y[kk]; A[i]=areas(dyjd_x[i][0],dyjd_y[i][0],dyjd_x[i][1],dyjd_y[i][1],dyjd_x[i][2],dyjd_y[i][0]); nn_a[i][0]=dyjd_x[i][1]*dyjd_y[i][2]-dyjd_x[i][2]*dyjd_y[i][1]; nn_a[i][1]=dyjd_x[i][2]*dyjd_y[i][0]-dyjd_x[i][0]*dyjd_y[i][2]; nn_a[i][2]=dyjd_x[i][0]*dyjd_y[i][1]-dyjd_x[i][1]*dyjd_y[i][0]; nn_b[i][0]=dyjd_y[i][1]-dyjd_y[i][1]; nn_b[i][1]=dyjd_y[i][2]-dyjd_y[i][1]; nn_b[i][2]=dyjd_y[i][0]-dyjd_y[i][1]; nn_c[i][0]=dyjd_x[i][2]-dyjd_x[i][1]; nn_c[i][1]=dyjd_x[i][0]-dyjd_x[i][2]; nn_c[i][2]=dyjd_x[i][1]-dyjd_x[i][0]; for (j=0;j<=3;j++) for(k=0;k<=3;k++) DY_DE[j][k]=0; for (j=0;j<=3;j++) for(k=0;k<=6;k++) DY_B[j][k][i]=0; ff1=E/(1-nou*nou); DY_DE[1][1]=ff1;DY_DE[1][2]=ff1*nou; DY_DE[2][1]=ff1*nou;DY_DE[2][2]=ff1; DY_DE[3][3]=ff1*0.5*(1-nou); ff1=0.5/A[i]; DY_B[1][1][i]=ff1*nn_b[i][0]; DY_B[1][3][i]=ff1*nn_b[i][1]; DY_B[1][5][i]=ff1*nn_b[i][2]; DY_B[2][2][i]=ff1*nn_c[i][0]; DY_B[2][4][i]=ff1*nn_c[i][1]; DY_B[2][6][i]=ff1*nn_c[i][2]; DY_B[3][1][i]=ff1*nn_c[i][0]; DY_B[3][3][i]=ff1*nn_c[i][1]; DY_B[3][5][i]=ff1*nn_c[i][2]; DY_B[3][2][i]=ff1*nn_b[i][0]; DY_B[3][4][i]=ff1*nn_b[i][1]; DY_B[3][6][i]=ff1*nn_b[i][2]; ff1=t*E/(4*A[i]*(1-nou*nou)); ff2=0.5*(1-nou); for (r=0;r<=2;r++) for(s=0;s<=2;s++) {
1 下一页