【正文】 有限元分析 課程設計 學 院： 土木建筑工程學院 專 業(yè)： 工程力學 班 級： 力學 131 學 號： 3121631023 姓 名 ： 崔晨露 一、 單元 劃分 與節(jié)點編號 二、 C語言程序 include include define NE 32 //單元數(shù) define NJ 45 //節(jié)點數(shù) define NZ 10 //支承數(shù) define NPJ 7 //節(jié)點荷載數(shù) define NJ2 90 //節(jié)點位移數(shù) define DD 14 //半帶寬 int LXM=0。 //類型判別碼 double EO=170e9。 //楊氏模量 double MU=。 //泊松比 double LOU=0。 //容重 double TE=。 //厚度 double AJZ [NJ+1][3]= {{0,0,0},{0,0,0},{0,0,},{0,0,},{0,0,},{0,0,3},{0,0},{0,},{0,},{0,},{0,3},{0,0},{0,},{0,},{0,},{0,3},{0,0},{0,5,},{0,},{0,},{0,3},{0,0},{0,5},{0,},{0,},{0,3},{0,0},{0,},{0,9.375,},{0,},{0,3},{0,0},{0,},{0,11.25,},{0,},{0,3},{0,0},{0,},{0,1,},{0,},{0,3},{0,15,0},{0,15,},{0,15,1.5},{0,15,},{0,15,3}}。//共 36 個節(jié)點 int JM [NE+1][5]= {{0,0,0,0},{0,1,6,7,2},{0,2,7,8,3},{0,3,8,9,4},{0,4,9,10,5},{0,6,11,12,7},{0,7,12,13,8},{0,8,13,14,9},{0,9,14,15,10},{0,11,16,17,12},{0,12,17,18,13},{0,13,18,19,14},{0,14,19,20,15},{0,16,21,22,17},{0,17,22,23,18},{0,18,23,24,19},{0,19,24,25,20},{0,21,26,27,22},{0,22,27,28,23},{0,23,28,29,24},{0,24,29,30,25},{0,26,31,32,27},{0,27,32,33,28},{0,28,33,34,29},{0,29,34,35,30},{0,31,36,37,32},{0,32,37,38,33},{0,33,38,39,34},{0,34,39,40,35},{0,36,41,42,37},{0,37,42,43,38},{0,38,43,44,39},{0,39,44,45,40}}。//共 32個單元 int NZC [NZ+1]={0,1,2,3,4,5,6,7,8,9,10}。//15號節(jié)點的 x,y 被約束 double PJ [NPJ+1][2+1]={{0,0,0},{0,50},{0,,60},{0,70},{0,80},{0,8e4,83},{0,8e4,87},{0,4,90}}。 double AE,JS[NJ*4][4],js[NJ+1][4],KZ[NJ2+1][DD+1],P[NJ2+1],S[3+1][8+1],KE[8+1][8+1],SZ[3+1][32+1]。 int IE,JE,ME,LE。 void DUGD(int,int)。//生成 S矩陣， KE 矩陣 void main() { int NJ1,k,IN,IM,jn,m,i,j,z,JO,ii,jj,h,dh,E,l,zl,dl,n。 double PE,c,SIG1,SIG2,SIG3,PYL,RYL,MAYL,MIYL,CETA。 double WY[8+1],YL[3+1]。 if(LXM!=0)//平面應力問題與平面應變問題的判別 { EO=EO/(*MU)。 MU=MU/()。 } for(i=0。i=NJ2。i++) { for(j=0。j=DD。j++) KZ[i][j]=。 } for(E=1。E=NE。E++) { DUGD(E,3)。 for(i=1。i=4。i++) { for(ii=1。ii=2。ii++) { h=2*(i1)+ii。 dh=2*(JM[E][i]1)+ii。 for(j=1。j=4。j++) { for(jj=1。jj=2。jj++) { l=2*(j1)+jj。 zl=2*(JM[E][j]1)+jj。 dl=zldh+1。 if(dl0) KZ[dh][dl]=KZ[dh][dl]+KE[h][l]。 } } } } } //********形成 P 矩陣 ******** for(i=1。i=NJ2。i++) P[i]=。 if(NPJ0) { for(i=1。i=NPJ。i++) { j=(int)PJ[i][2]。 P[j]=PJ[i][1]。 } } if(LOU0) { for(E=1。E=NE。E++) { DUGD(E,1)。 PE=LOU*(AE)*TE*。 P[2*IE]=P[2*IE]+PE。 P[2*JE]=P[2*JE]+PE。 P[2*ME]=P[2*ME]+PE。 P[2*LE]=P[2*LE]+PE。 } } //********邊界條件 ******** for(i=1。i=NZ。i++) { z=NZC[i]。 KZ[z][1]=。 for(j=2。j=DD。j++) KZ[z][j]=。 if(z!=1) { if(zDD) JO=DD。 else JO=z。 for(j=2。j=JO。j++) KZ[zj+1][j]=。 } P[z]=。 } //********求解方程 ******** NJ1=NJ21。 for(k=1。k=NJ1。k++) { if(NJ2k+DD1) IM=k+DD1。 else IM=NJ2。 IN=k+1。 for(i=IN。i=IM。i++) { l=ik+1。 c=KZ[k][l]/KZ[k][1]。 jn=DDl+1。 for(j=1。j=jn。j++) { m=j+ik。 KZ[i][j]=KZ[i][j]c*KZ[k][m]。 } P[i]=P[i]c*P[k]。 } } //*******矩陣回代 ******** P[NJ2]=P[NJ2]/KZ[NJ2][1]。 for(i=NJ1。i=1。i) { if(DDNJ2i+1) JO=NJ2i+1。 else JO=DD。 for(j=2。j=JO。j++) { h=j+i1。 P[i]=P[i]KZ[i][j]*P[h]。 } P[i]=P[i]/KZ[i][1]。 } printf(\n)。 printf(JD U V\n)。 for(i=1