【正文】 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= 。 a41=。 a37= 。 a33=。 a29= 。 a25= 。 a21=。 a17= 。 a13= 。 a9=。 a5= 。 } } 節(jié)點(diǎn)位移 JD U V 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 a1= 。 KE[8][1]=KE[1][8]。 KE[7][1]=KE[1][7]。 KE[6][1]=KE[1][6]。 KE[5][1]=KE[1][5]。 KE[4][1]=KE[1][4]。 KE[3][1]=KE[1][3]。 KE[2][1]=KE[1][2]。 KE[1][1]=C*(b/(3*a)+(1MU)*a/(6*b))。j=8。 B[3][4]=(by)/AE。 B[1][7]=(b+y)/AE。i++) for(j=0。 D[2][3]=0。 a=。 JDYL[r][5]=JDYL[r][3]+JDYL[r][4]。 JDYL[r][1]=JDYL[r][1]+DYYL[f][n][1]。rNJ。 fprintf(fp2,E=%d(%d)\n,E,n)。 MIYL=PYLRYL。j=8。 dh=2*(JM[E][i]1)+j。 } for(E=1。i=NJ。j=JO。 KZ[i][j]=KZ[i][j]c*KZ[k][m]。 for(i=IN。j++) KZ[zj+1][j]=。 KZ[z][1]=。 F=LOU*(AE)*TE/4。i++) P[i]=。j++) { for(jj=1。 for(i=1。 } for(i=0。 void main() { int NJ1,k,IN,IM,jn,m,i,j,z,JO,ii,jj,h,dh,E,l,zl,dl,r,n,o,f。 int NZC[NZ+1]={0,1,2,21,22,23,24,43,44,45,46,65,66,67,68,87,88 }。 、應(yīng)力分布圖、單元?jiǎng)澐謭D。 。 //節(jié)點(diǎn)坐標(biāo) int JM[NE+1][5]={{0,0,0,0,0},{0,10,11,22,21},{0,9,10,21,20},{0,8,9,20,19},{0,7,8,19,18},{0,6,7,18,17},{0,5,6,17,16},{0,4,5,16,15},{0,3,4,15,14},{0,2,3,14,13},{0,1,2,13,12},{0,21,22,33,32},{0,20,21,32,31},{0,19,20,31,30},{0,18,19,30,29},{0,17,18,29,28},{0,16,17,28,27},{0,15,16,27,26},{0,14,15,26,25},{0,13,14,25,24},{0,12,13,24,23},{0,32,33,44,43},{0,31,32,43,42},{0,30,31,42,41},{0,29,30,41,40},{0,28,29,40,39},{0,27,28,39,38},{0,26,27,38,37},{0,25,26,37,36},{0,24,25,36,35},{0,23,24,35,34}}。 FILE *fp1,*fp2,*ab。 MU=MU/()。E++) { DUGD(E,3)。j=4。i=NJ2。E++) { DUGD(E,1)。i++) { z=NZC[i]。j=JO。 IN=k+1。j++) { m=j+ik。 for(j=2。 for(i=1。\n,i,AJZ[i][1]+P[2*i1],i,AJZ[i][2]+P[2*i])。j++) { h=2*(i1)+j。 for(j=1。 MAYL=PYL+RYL。 fprintf(fp2,\n)。 for(r=0。f++) if(JM[n+1][f+1]==r+1)