【正文】 = 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= 。 a44=。 a42=。 a40= 。 a38= 。 a36= 。 a34= 。 a32=。 a30=。 a28= 。 a26= 。 a24= 。 a22=。 a20=。 a18= 。 a16= 。 a14= 。 a12= 。 a10=。 a8=。 a6= 。 a4= 。 a2= 。 KE[8][7]=KE[7][8]。 KE[8][3]=KE[3][8]。 KE[7][7]=KE[1][1]。 KE[7][3]=KE[3][7]。 KE[6][7]=KE[1][8]。 KE[6][3]=KE[3][6]。 KE[5][7]=KE[1][3]。 KE[5][3]=KE[3][5]。 KE[4][7]=KE[1][2]。 KE[4][3]=KE[3][4]。 KE[3][7]=KE[1][5]。 KE[3][3]=KE[1][1]。 KE[2][7]=KE[1][4]。 KE[2][3]=KE[1][8]。 KE[1][7]=C*(b/(6*a)(1MU)*a/(6*b))。 KE[1][3]=C*((b)/(3*a)+(1MU)*a/(12*b))。 SZ[i][n]=S[i][j]。 for(k=1。i=3。 B[3][6]=(b+y)/AE。 B[3][2]=(yb)/AE。 B[2][4]=(a+x)/AE。 B[1][3]=(by)/AE。j++) B[i][j]=。w++) { for(i=0。 D[3][2]=0。 D[2][1]=D[1][2]。 AE=4*a*b。 int i, j, k,n,w。 //輸出節(jié)點(diǎn)應(yīng)力至文件 } fclose(ab)。 JDYL[r][3]=(JDYL[r][0]+JDYL[r][1])*。 } //printf(%f\n%f\n%f\n\n\n,JDYL[r][0],JDYL[r][1],JDYL[r][2])。f5。 for(n=0。 } } } fprintf(ab,各節(jié)點(diǎn)應(yīng)力： \n)。 fprintf(fp2,s1=% s2=% theta=%\n,MAYL,MIYL,CETA)。 printf(s1=% s2=% theta=%\n,MAYL,MIYL,CETA)。 else CETA=*atan2(SIG3,(SIG2MIYL))。 RYL=sqrt(pow((SIG1SIG2)/,2)+pow(SIG3,2))。 } SIG1=YL[1]。i++) { YL[i]=0。 } } for(n=1。j=2。E++) { DUGD(E,2)。b%d=%。 fprintf(fp1,%d % %\n,i,P[2*i1],P[2*i])。 fprintf(fp1,JD U V\n)。 P[i]=P[i]KZ[i][j]*P[h]。 else JO=DD。 } } P[NJ2]=P[NJ2]/KZ[NJ2][1]。j=jn。i++) { l=ik+1。 else IM=NJ2。 } NJ1=NJ21。 for(j=2。j=DD。i=NZ。 P[2*JE]=P[2*JE]+F。E=NE。i=NPJ。 } } } } } for(i=1。jj++) { l=2*(j1)+jj。 for(j=1。i++) { for(ii=1。E=NE。i++) { for(j=0。 if(LXM!=0) { E0=E0/(*MU)。 double WY[8+1],YL[3+1]。 void DUGD(int,int)。 //節(jié)點(diǎn)載荷數(shù)組 //節(jié)點(diǎn)載荷數(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 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}}。 double E0=2e7。 要求： 20 個(gè)。用有限元方法計(jì)算梁的變形及應(yīng)力分布，要求用矩形單元。 二、力學(xué)分析 LXM=0 q=10kN/ 三．程序框圖及程序 程序框圖 ： 開始 插入基本參數(shù) 插入其他參數(shù) 形成整體剛度矩陣 形成荷載列陣 引入支承條件 解方程輸出位移 求應(yīng)力輸出應(yīng)力 結(jié)束 子程序單元?jiǎng)偠染仃? （ ASK=3）計(jì)算單元?jiǎng)偠染仃? （ ASK=1）計(jì)算單元面積 （ ASK=2）計(jì)算 S 矩陣 四． 源程序 include include define NE 30 //單元數(shù) define NJ 44 //節(jié)點(diǎn)數(shù) define NZ 16 //支承數(shù) define NPJ 11 //節(jié)點(diǎn)載荷作用數(shù) define DD 26 //半帶寬 define NJ2 88 //節(jié)點(diǎn)位移數(shù) int LXM=0。 double TE=1 。 //支撐數(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}}。 //單元應(yīng)力矩陣 int IE,JE,ME,PE。 double F,c,SIG1,SIG2,SIG3,PYL,RYL,MAYL,MIYL,CETA。 fp2=fopen(單元應(yīng)力 .txt,w)。i=NJ2。 } for(E=1。i=4。 dh=2*(JM[E][i]1)+ii。jj=2。 if(dl0) KZ[dh][dl]=KZ[dh][dl]+KE[h][l]。 if(NPJ0) { for(i=1。 } } if(LOU0) { for(E=1。 P[2*IE]=P[2*IE]+F。 } } for(i=1。 for(j=2。 else JO=z。 } P[z]=。k++) { if(NJ2k+DD1) IM=k+DD1。i=IM。 for(j=1。 } P[i]=P[i]c*P[k]。i) { if(DD=NJ2i+1) JO=NJ2i+1。j++) { h=j+i1。 printf(JD U V\n)。i++) { printf(%d % %\n,i,P[2*i1],P[2*i])。i++) { fprintf(fp1,a%d=%。E=NE。i++) { for(j=1。 WY[h]=P[dh]。i=3。j++) YL[i]=YL[i]+SZ[i][8*(n1)+j]*WY[j]。 PYL=(SIG1+SIG2)/2。 if(SIG2==MIYL) CETA=0。 printf(sx=% sy=% tou=%\n,SIG1,SIG2,SIG3)。 fprintf(fp2,sx=% sy=% tou=%\n,SIG1,SIG2,SIG3)。 DYYL[n1][E1][2]=SIG3。r++) //計(jì)算節(jié)點(diǎn)應(yīng)力（始） { o=0。 for(f=0。 JDYL[r][2]=JDYL[r][2]+DYYL[f][n][2]。 JDYL[r][2]=JDYL[r][2]*。 //計(jì)算節(jié)點(diǎn)應(yīng)力（終） fprintf(ab,%2d\t%f\n,r+1,JDYL[r][5]/1000000)。 } /**************DUGD 子程序 ******************/ void DUGD(int E,int ASK) { double B[3+1][8+1],D[3+1][3+1],x,y,a,b,C。 b=。 D[1][3]=0。 D[3][1]=0。w=5。j=8。 B[1][1]=(yb)/AE。 B[2][2]=(xa)/AE。 B[3][1]=(xa)/AE。 B[3][5]=(a+x)/AE。 for(i=1。j++) { S[i][j]=。 n=8*(w1)+j。 KE[1][2]=C*(1+MU)/8。 KE[1][6]=(C)*(1+MU)/8。 KE[2][2]=C*(a/(3*b)+(1MU)*b/(6*a))。 KE[2][6]=C*((a)/(6*b)(1MU)*b/(12*a))。 KE[3][2]=KE[2][3]。 KE[3][6]=KE[1][4]。 KE[4][2]=KE[2][4]。 KE[4][6]=KE[2][8]。 KE[5][2]=KE[2][5]。 KE[5][6]=KE[1][2]。 KE[6][2]=KE[2][6]。 KE[6][6]=KE[2][2]。 KE[7][2]=KE[2][7]。 KE[7][6]=KE[6][7]。 KE[8][2]=KE[2][8]。 KE[8][6]=KE[6][8]。b1= 。b3= 。b5= 。b7= 。b9= 。b11= 。b13= 。b15= 。b17= 。b19= 。b21= 。b23= 。b25= 。b27= 。b29= 。b31= 。b33= 。b35=。b37=。b39=。b41=。b