【正文】 ) { for(ii=1。 double WY[8+1],YL[3+1]。 double E0=2e7。 double TE=1 。 fp2=fopen(單元應力 .txt,w)。 dh=2*(JM[E][i]1)+ii。 } } if(LOU0) { for(E=1。 else JO=z。 for(j=1。 printf(JD U V\n)。i++) { for(j=1。 PYL=(SIG1+SIG2)/2。 DYYL[n1][E1][2]=SIG3。 JDYL[r][2]=JDYL[r][2]*。 D[1][3]=0。 B[1][1]=(yb)/AE。 for(i=1。 KE[1][6]=(C)*(1+MU)/8。 KE[3][6]=KE[1][4]。 KE[5][6]=KE[1][2]。 KE[7][6]=KE[6][7]。b3= 。b11= 。b19= 。b27= 。b35=。b43=。b39=。b31= 。b23= 。b15= 。b7= 。 KE[8][6]=KE[6][8]。 KE[6][6]=KE[2][2]。 KE[4][6]=KE[2][8]。 KE[2][6]=C*((a)/(6*b)(1MU)*b/(12*a))。 n=8*(w1)+j。 B[3][1]=(xa)/AE。w=5。 } /**************DUGD 子程序 ******************/ void DUGD(int E,int ASK) { double B[3+1][8+1],D[3+1][3+1],x,y,a,b,C。 for(f=0。 printf(sx=% sy=% tou=%\n,SIG1,SIG2,SIG3)。i=3。i++) { fprintf(fp1,a%d=%。i) { if(DD=NJ2i+1) JO=NJ2i+1。k++) { if(NJ2k+DD1) IM=k+DD1。 } } for(i=1。 if(dl0) KZ[dh][dl]=KZ[dh][dl]+KE[h][l]。 } for(E=1。 //單元應力矩陣 int IE,JE,ME,PE。用有限元方法計算梁的變形及應力分布，要求用矩形單元。 //節(jié)點載荷數(shù)組 //節(jié)點載荷數(shù)組 double AE,KZ[NJ2+1][DD+1],P[NJ2+1],S[3+1][8+1],KE[8+1][8+1],SZ[3+1][5*8+1]。i++) { for(j=0。jj++) { l=2*(j1)+jj。 P[2*JE]=P[2*JE]+F。 } NJ1=NJ21。 } } P[NJ2]=P[NJ2]/KZ[NJ2][1]。 fprintf(fp1,%d % %\n,i,P[2*i1],P[2*i])。 } } for(n=1。 else CETA=*atan2(SIG3,(SIG2MIYL))。 for(n=0。 //輸出節(jié)點應力至文件 } fclose(ab)。 D[3][2]=0。 B[2][4]=(a+x)/AE。 for(k=1。 KE[2][3]=KE[1][8]。 KE[4][3]=KE[3][4]。 KE[6][3]=KE[3][6]。 KE[8][3]=KE[3][8]。 a6= 。 a14= 。 a22=。 a30=。 a38= 。 單元應力 E=1(1) sx= sy= tou= s1= s2= theta= E=1(2) sx= sy= tou= s1= s2= theta= E=1(3) sx= sy= tou= s1= s2= theta= E=1(4) sx= sy= tou= s1= s2= theta= E=1(5) sx= sy= tou= s1= s2= theta= E=2(1) sx= sy= tou= s1= s2= theta= E=2(2) sx= sy= tou= s1= s2= theta= E=2(3) sx= sy= tou= s1= s2= theta= E=2(4) sx= sy= tou= s1= s2= theta= E=2(5) sx= sy= tou= s1= s2= theta= E=3(1) sx= sy= tou= s1= s2= theta= E=3(2) sx= sy= tou= s1= s2= theta= E=3(3) sx= sy= tou= s1= s2= theta= E=3(4) sx= sy= tou= s1= s2= theta= E=3(5) sx= sy= tou= s1= s2= theta= E=4(1) sx= sy= tou= s1= s2= theta= E=4(2) sx= sy= tou= s1= s2= theta= E=4(3) sx= sy= tou= s1= s2= theta= E=4(4) sx= sy= tou= s1= s2= theta= E=4(5) sx= sy= tou= s1= s2= theta= E=5(1) sx= sy= tou= s1= s2= theta= E=5(2) sx= sy= tou= s1= s2= theta= E=5(3) sx= sy= tou= s1= s2= theta= E=5(4) sx= sy= tou= s1= s2= theta= E=5(5) sx= sy= tou= s1= s2= theta= E=6(1) sx= sy= tou= s1= s2= theta= E=6(2) sx= sy= tou= s1= s2= theta= E=6(3) sx= sy= tou= s1= s2= theta= E=6(4) sx= sy= tou= s1= s2= theta= E=6(5) sx= sy= tou= s1= s2= theta= E=7(1) sx= sy= tou= s1= s2= theta= E=7(2) sx= sy= tou= s1= s2= theta= E=7(3) sx= sy= tou= s1= s2= theta= E=7(4) sx= sy= tou= s1= s2= theta= E=7(5) sx= sy= tou= s1= s2= theta= E=8(1) sx= sy= tou= s1= s2= theta= E=8(2) sx= sy= tou= s1= s2= theta= E=8(3) sx= sy= tou= s1= s2= theta= E=8(4) sx= sy= tou= s1= s2= theta= E=8(5) sx= sy= tou= s1= s2= theta= E=9(1) sx= sy= tou= s1= s2= theta= E=9(2) sx= sy= tou= s1= s2= theta= E=9(3) sx= sy= tou= s1= s2= theta= E=9(4) sx= sy= tou= s1= s2= theta= E=9(5) sx= sy= tou= s1= s2= theta= E=10(1) sx= sy= tou= s1= s2= theta= E=10(2) sx= sy= tou= s1= s2= theta= E=10(3) sx= sy= tou= s1= s2= theta= E=10(4) sx= sy= tou= s1= s2= theta= E=10(5) sx= sy= tou= s1= s2= theta= E=11(1) sx= sy= tou= s1= s2= theta= E=11(2) sx= sy= tou= s1= s2= theta= E=11(3) sx= sy= tou= s1= s2= theta= E=11(4) sx= sy= tou= s1= s2= theta= E=11(5) sx= sy= tou= s1= s2= theta= E=12(1) sx= sy= tou= s1= s2= theta= E=12(2) sx= sy= tou= s1= s2= theta= E=12(3) sx= sy= tou= s1= s2= theta= E=12(4) sx= sy= tou= s1= s2= theta= E=12(5) sx= sy= tou= s1= s2= theta= E=13(1) sx= sy= tou= s1= s2= theta= E=13(2) sx= sy= tou= s1= s2= theta= E=13(3) sx= sy= tou= s1= s2= theta= E=13(4) sx= sy= tou= s1= s2= theta= E=13(5) sx= sy= tou= s1= s2= theta= E=14(1) sx= sy= tou= s1= s2= theta= E=14(2) sx= sy= tou= s1= s2= theta= E=14(3) sx= sy= tou= s1= s2= theta= E=14(4) sx= sy= tou= s1= s2= theta= E=14(5) sx= sy= tou= s1= s2= theta= E=15(1) sx= sy= tou= s1= s2= theta= E=15(2) sx= sy= tou= s1= s2= theta= E=15(3) sx= sy= tou= s1= s2= theta= E=15(4) sx= sy= tou= s1= s2= theta= E=15(5) sx= sy= tou= s1= s2= theta= E=16(1) sx= sy= tou= s1= s2= theta= E=16(2) sx= sy= tou= s1= s2= theta= E=16(3) sx= sy= tou= s1= s2= theta= E=16(4) sx= sy= tou= s1= s2= theta= E=16(5) sx= sy= tou= s1= s2= theta= E=17(1) sx= sy= tou= s1= s2= theta= E=17(2) sx= sy= tou= s1= s2= theta= E=17(3) sx= sy= tou= s1= s2= theta=