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