【正文】 { JDYL[r][0]=JDYL[r][0]+DYYL[f][n][0]。 JDYL[r][4]=(float)sqrt(pow((JDYL[r][0]JDYL[r][1])*,2)+pow(JDYL[r][2],2))。 double xy[6][3]={{0,0,0},{0,1,1},{0,1,1},{0,1,1},{0,1,1},{0,0,0}}。 D[2][2]=D[1][1]。i=3。 B[1][5]=(b+y)/AE。 B[3][3]=(a+x)/AE。i++) { for(j=1。 } } } } if(ASK2) { C=E0*TE/(1MU*MU)。 KE[1][8]=C*((13*MU)/8)。 KE[2][8]=C*((a)/(3*b)+(1MU)*b/(12*a))。 KE[3][8]=KE[1][2]。 KE[4][8]=KE[2][6]。 KE[5][8]=KE[1][4]。 KE[6][8]=KE[2][4]。 KE[7][8]=KE[1][6]。 KE[8][8]=KE[2][2]。b4= 。b8= 。b12= 。b16= 。b20= 。b24= 。b28= 。b32= 。b36=。b40=。b44= 。b42=。b38=。b34= 。b30= 。b26= 。b22= 。b18= 。b14= 。b10= 。b6= 。b2= 。 KE[8][4]=KE[4][8]。 KE[7][4]=KE[4][7]。 KE[6][4]=KE[4][6]。 KE[5][4]=KE[4][5]。 KE[4][4]=KE[2][2]。 KE[3][4]=KE[1][6]。 KE[2][4]=C*(a/(6*b)(1MU)*b/(6*a))。 KE[1][4]=C*((3*MU1)/8)。k=3。 B[3][7]=(ax)/AE。 B[2][6]=(a+x)/AE。 x=xy[w][1]*a。 D[3][3]=E0/(2*(1+MU))。 D[1][1]=E0/(1MU*MU)。 //關閉文件 fclose(fp1)。 //檢查單元應力各項是否正確求和 } JDYL[r][0]=JDYL[r][0]*。nNE。 if(n5) { DYYL[n1][E1][0]=SIG1。 printf(\n)。 SIG2=YL[2]。n=5。 for(i=1。 } for(i=1。 } P[i]=P[i]/KZ[i][1]。 for(i=NJ1。 c=KZ[k][l]/KZ[k][1]。 for(k=1。j++) KZ[z][j]=。 P[2*ME]=P[2*ME]+F。i++) { j=(int)PJ[i][2]。 zl=2*(JM[E][j]1)+jj。ii=2。j=DD。 ab=fopen(節(jié)點應力 .txt,w)。 double JDYL[NJ][6]。 double MU=。 一． 題目 如圖所示固支梁，高 3 米，長 15m，承受均布載荷 q=10kN/， E=20GPa，u=，厚度 t=1m，忽略自重，按平面應力問題分析。 double LOU=。 //節(jié)點應力矩陣 double DYYL[4][NE][4]。 fp1=fopen(節(jié)點位移 .txt,w)。j++) KZ[i][j]=。ii++) { h=2*(i1)+ii。 dl=zldh+1。 P[j]=PJ[i][1]。 P[2*PE]=P[2*PE]+F。 if(z!=1) { if(zDD) JO=DD。k=NJ1。 jn=DDl+1。i=1。 } printf(\n)。i=NJ。i=4。n++) { for(i=1。 SIG3=YL[3]。 printf(E=%d(%d)\n,E,n)。 //將各單元應力記錄到單元應力矩陣中 DYYL[n1][E1][1]=SIG2。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。 JDYL[r][1]=JDYL[r][1]*。 fclose(fp2)。 D[1][2]=MU*E0/(1MU*MU)。 if(ASK1) { for(w=1。 y=xy[w][2]*b。 B[2][8]=(ax)/AE。 B[3][8]=(b+y)/AE。k++) S[i][j]=S[i][j]+D[i][k]*B[k][j]。 KE[1][5]=C*((b)/(6*a)(1MU)*a/(12*b))。 KE[2][5]=KE[1][6]。 KE[3][5]=KE[1][7]。 KE[4][5]=KE[2][3]。 KE[5][5]=KE[1][1]。 KE[6][5]=KE[5][6]。 KE[7][5]=KE[5][7]。 KE[8][5]=KE[5][8]。 a3= 。 a7= 。 a11=。 a15= 。 a19=。 a23= 。 a27= 。 a31=。 a35= 。 a39= 。 a43=。 a44=。 a40= 。 a36= 。 a32=。 a28= 。 a24= 。 a20=。 a16= 。 a12= 。 a8=。 a4= 。 KE[8][7]=KE[7][8]。 KE[7][7]=KE[1][1]。 KE[6][7]=KE[1][8]。 KE[5][7]=KE[1][3]。 KE[4][7]=KE[1][2]。 KE[3][7]=KE[1][5]。 KE[2][7]=KE[1][4]。 KE[1][7]=C*(b/(6*a)(1MU)*a/(6*b))。 SZ[i][n]=S[i][j]。i=3。 B[3][2]=(yb)/AE。 B[1][3]=(by)/AE。w++) { for(i=0。 D[2][1]=D[1][2]。 int i, j, k,n,w。 JDYL[r][3]=(JDYL[r][0]+JDYL[r][1])*。f5。 } } } fprintf(ab,各節(jié)點應力： \n)。 printf(s1=% s2=% theta=%\n,MAYL,MIYL,CETA)。 RYL=sqrt(pow((SIG1SIG2)/,2)+pow(SIG3,2))。i++) { YL[i]=0。j=2。b%d=%。 fprintf(fp1,JD U V\n)。 else JO=DD。j=jn。 else IM=NJ2。 for(j=2。i=NZ。E=NE。 } } } } } for(i=1。 for(j=1。E=NE。 if(LXM!=0) { E0=E0/(*MU)。 void DUGD(int,int)。 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}}。 要求： 20 個。 二、力學分析 LXM=0 q=10kN/ 三．程序框圖及程序 程序框圖 ： 開始 插入基本參數(shù) 插入其他參數(shù) 形成整體剛度矩陣 形成荷載列陣 引入支承條件 解方程輸出位移 求應力輸出應力 結(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。 //支撐數(shù)組 double PJ[NPJ+1][2+1]={{0,0,0},{0,7500,68},{0,15000,70},{0,15000,72},{0,15000,74},{0,15000,76},{0,15000,78},{0,15000,80},{0,15000,82},{0,15000,84},{0,15000,86},{0,7500,88}}。 double F,c,SIG1,SIG2,SIG3,PYL,RYL,MAYL,MIYL,CETA。i=NJ2。i=4。jj=2。 if(NPJ0) { for(i=1。 P[2*IE]=P[2*IE]+F。 for(j=2。 } P[z]=。i=IM。 } P[i]=P[i]c*P[k]。j++) { h=j+i1。i++) { printf(%d % %\n,i,P[2*i1],P[2*i])。E=NE。 WY[h]=P[dh]。j++) YL[i]=YL[i]+SZ[i][8*(n1)+j]*WY[j]。 if(SIG2==MIYL) CETA=0。 fprintf(fp2,sx=% sy=% tou=%\n,SIG1,SIG2,SIG3)。r++) //計算節(jié)點應力（始） { o=0。 JDYL[r][2]=JDYL[r][2]+DYYL[f][n][2]。 //計算節(jié)點應力（終） fprintf(ab,%2d\t%f\n,r+1,JDYL[r][5]/1000000)。 b=。 D[3][1]=0。j=8。 B[2][2]=(xa)/AE。 B[3][5]=(a+x)/AE。j++) { S[i][j]=。 KE[1][2]=C*(1+MU)/8。 KE[2][2]=C*(a/(3*b)+(1MU)*b/(6*a))。 KE[3][2]=KE[2][3]。 KE[4][2]=KE[2][4]。 KE[5][2]=KE[2][5]。 KE[6][2]=KE[2][6]。 KE[7][2]=KE[2][7]。 KE[8][2]=KE[2][8]。b1= 。b5= 。b9= 。b13= 。b17= 。b21= 。b25= 。b29= 。b33= 。b37=。b