【正文】 b41=。b33= 。b25= 。b17= 。b9= 。b1= 。 KE[7][2]=KE[2][7]。 KE[5][2]=KE[2][5]。 KE[3][2]=KE[2][3]。 KE[1][2]=C*(1+MU)/8。 B[3][5]=(a+x)/AE。j=8。 b=。 JDYL[r][2]=JDYL[r][2]+DYYL[f][n][2]。 fprintf(fp2,sx=% sy=% tou=%\n,SIG1,SIG2,SIG3)。j++) YL[i]=YL[i]+SZ[i][8*(n1)+j]*WY[j]。E=NE。j++) { h=j+i1。i=IM。 for(j=2。 if(NPJ0) { for(i=1。i=4。 double F,c,SIG1,SIG2,SIG3,PYL,RYL,MAYL,MIYL,CETA。 二、力學(xué)分析 LXM=0 q=10kN/ 三．程序框圖及程序 程序框圖 ： 開始 插入基本參數(shù) 插入其他參數(shù) 形成整體剛度矩陣 形成荷載列陣 引入支承條件 解方程輸出位移 求應(yīng)力輸出應(yīng)力 結(jié)束 子程序單元剛度矩陣 （ ASK=3）計算單元剛度矩陣 （ ASK=1）計算單元面積 （ ASK=2）計算 S 矩陣 四． 源程序 include include define NE 30 //單元數(shù) define NJ 44 //節(jié)點數(shù) define NZ 16 //支承數(shù) define NPJ 11 //節(jié)點載荷作用數(shù) define DD 26 //半帶寬 define NJ2 88 //節(jié)點位移數(shù) int LXM=0。 double AJZ[NJ+1][3]={{0,0,0}, {0,0,3},{0,3},{0,3,3},{0,3},{0,6,3},{0,3},{0,9,3},{0,3},{0,12,3},{0,3},{0,15,3}, {0,0,2},{0,2},{0,3,2},{0,2},{0,6,2},{0,2},{0,9,2},{0,2},{0,12,2},{0,2},{0,15,2}, {0,0,1},{0,1},{0,3,1},{0,1},{0,6,1},{0,1},{0,9,1},{0,1},{0,12,1},{0,1},{0,15,1}, {0,0,0},{0,0},{0,3,0},{0,0},{0,6,0},{0,0},{0,9,0},{0,0},{0,12,0},{0,0},{0,15,0}}。 if(LXM!=0) { E0=E0/(*MU)。 for(j=1。E=NE。 for(j=2。j=jn。 fprintf(fp1,JD U V\n)。j=2。 RYL=sqrt(pow((SIG1SIG2)/,2)+pow(SIG3,2))。 } } } fprintf(ab,各節(jié)點應(yīng)力： \n)。 JDYL[r][3]=(JDYL[r][0]+JDYL[r][1])*。 D[2][1]=D[1][2]。 B[1][3]=(by)/AE。i=3。 KE[1][7]=C*(b/(6*a)(1MU)*a/(6*b))。 KE[3][7]=KE[1][5]。 KE[5][7]=KE[1][3]。 KE[7][7]=KE[1][1]。 a4= 。 a12= 。 a20=。 a28= 。 a36= 。 a44=。 a39= 。 a31=。 a23= 。 a15= 。 a7= 。 KE[8][5]=KE[5][8]。 KE[6][5]=KE[5][6]。 KE[4][5]=KE[2][3]。 KE[2][5]=KE[1][6]。k++) S[i][j]=S[i][j]+D[i][k]*B[k][j]。 B[2][8]=(ax)/AE。 if(ASK1) { for(w=1。 fclose(fp2)。n++) if(JM[n+1][1]==r+1||JM[n+1][2]==r+1||JM[n+1][3]==r+1||JM[n+1][4]==r+1) { o=o+1。 printf(E=%d(%d)\n,E,n)。n++) { for(i=1。i=NJ。i=1。k=NJ1。 P[2*PE]=P[2*PE]+F。 dl=zldh+1。j++) KZ[i][j]=。 //節(jié)點應(yīng)力矩陣 double DYYL[4][NE][4]。 一． 題目 如圖所示固支梁，高 3 米，長 15m，承受均布載荷 q=10kN/， E=20GPa，u=，厚度 t=1m，忽略自重，按平面應(yīng)力問題分析。 double JDYL[NJ][6]。j=DD。 zl=2*(JM[E][j]1)+jj。 P[2*ME]=P[2*ME]+F。 for(k=1。 for(i=NJ1。 } for(i=1。n=5。 printf(\n)。nNE。 //關(guān)閉文件 fclose(fp1)。 D[3][3]=E0/(2*(1+MU))。 B[2][6]=(a+x)/AE。k=3。 KE[2][4]=C*(a/(6*b)(1MU)*b/(6*a))。 KE[4][4]=KE[2][2]。 KE[6][4]=KE[4][6]。 KE[8][4]=KE[4][8]。b6= 。b14= 。b22= 。b30= 。b38=。b44= 。b36=。b28= 。b20= 。b12= 。b4= 。 KE[7][8]=KE[1][6]。 KE[5][8]=KE[1][4]。 KE[3][8]=KE[1][2]。 KE[1][8]=C*((13*MU)/8)。i++) { for(j=1。 B[1][5]=(b+y)/AE。 D[2][2]=D[1][1]。 JDYL[r][4]=(float)sqrt(pow((JDYL[r][0]JDYL[r][1])*,2)+pow(JDYL[r][2],2))。 for(r=0。 MAYL=PYL+RYL。j++) { h=2*(i1)+j。 for(i=1。j++) { m=j+ik。j=JO。E++) { DUGD(E,1)。j=4。 MU=MU/()。 //節(jié)點坐標(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}}。 、應(yīng)力分布圖、單元劃分圖。 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。 for(i=1。i++) P[i]=。 KZ[z][1]=。 for(i=IN。j=JO。 } for(E=1。j=8。 fprintf(fp2,E=%d(%d)\n,E,n)。 JDYL[r][1]=JDYL[r][1]+DYYL[f][n][1]。 a=。i++) for(j=0。 B[3][4]=(by)/AE。 KE[1][1]=C*(b/(3*a)+(1MU)*a/(6*b))。 KE[3][1]=KE[1][3]。 KE[5][1]=KE[1][5]。 KE[7][1]=KE[1][7]。 } } 節(jié)點位移 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= 。 a9=。 a17= 。 a25= 。 a33=。 a41=。 a42=。 a34= 。 a26= 。 a18= 。 a10=。 a2= 。 KE[7][3]=KE[3][7]。 KE[5][3]=KE[3][5]。 KE[3][3]=KE[1][1]。 KE[1][3]=C*((b)/(3*a)+(1MU)*a/(12*b))。 B[3][6]=(b+y)/AE。j++) B[i][j]=。 AE=4*a*b。 } //printf(%f\n%f\n%f\n\n\n,JDYL[r][0],JDYL[r][1],JDYL[r][2])。 fprintf(fp2,s1=% s2=% theta=%\n,MAYL,MIYL,CETA)。 } SIG1=YL[1]。E++) { DUGD(E,2)。 P[i]=P[i]KZ[i][j]*P[h]。i++) { l=ik+1。j=DD。i=NPJ。i++