【文章內容簡介】 enf = sig 。returnfunction DisplayResults% 顯示計算結果，并與解析解結果比較 global gNode gElement gMaterial gDelta gElementStress outfile gBC1 global gDF n m r p b Stress outfile=fopen(39。39。,39。w39。)。 Stress = zeros((m+1)*(n+1),3)。 %集成節(jié)點應力矩陣 for k=1:1:3 for i=1:1:m Stress ((n+1)*(i1)+1,k)=gElementStress(n*(i1)+1,4,k)。 end for j=1:1:n Stress ((n+1)*m+j,k)=gElementStress(n*(m1)+j,1,k)。 end Stress ((n+1)*(m+1),k)=gElementStress(n*m,2,k)。 for i=1:1:m for j=1:1:n Stress ((i1)*(n+1)+j+1,k)=gElementStress((i1)*n+j,3,k)。 end end end xx_max=max(max(Stress(:,1)))。 N_xx_max=find(Stress==xx_max)。 xx_min=min(min(Stress(:,1)))。 N_xx_min=find(Stress==xx_min)。 yy_max=max(max(Stress(:,2)))。 N_yy_max=find(Stress==yy_max)(m+1)*(n+1)。 yy_min=min(min(Stress(:,2)))。 N_yy_min=find(Stress==yy_min)(m+1)*(n+1)。 xy_max=max(max(Stress(:,3)))。 N_xy_max=find(Stress==xy_max)2*(m+1)*(n+1)。 xy_min=min(min(Stress(:,3)))。 N_xy_min=find(Stress==xy_min)2*(m+1)*(n+1)。 fprintf(outfile, 39。 表0 節(jié)點坐標（mm) \n39。 ) 。 fprintf(outfile, 39。============================\n39。 ) 。 fprintf(outfile, 39。結點號 X坐標 Y坐標 \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 for i=1:size(gNode(:,1)) fprintf(outfile, 39。%4d , % , % \n39。,i,gNode(i,1),gNode(i,2)) 。 end fprintf( outfile,39。============================\n\n\n39。 ) 。 fprintf(outfile, 39。 表1 節(jié)點位移、應力 \n39。 ) 。 fprintf(outfile, 39。============================================================================\n39。 ) 。 fprintf(outfile, 39。結點號 x位移 y位移 Stressxx Stressyy Stressxy \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 for i=1:size(gNode(:,1)) fprintf(outfile, 39。%4d % % % % %\n39。,... i, full(gDelta(2*i1)),full(gDelta(2*i)),Stress(i,1),Stress(i,2),Stress(i,3)) 。 end fprintf( outfile,39。============================================================================\n\n\n39。 ) 。 fprintf(outfile, 39。 表2 單元節(jié)點序列 \n39。 ) 。 fprintf(outfile, 39。==========================================\n39。 ) 。 fprintf(outfile, 39。單元號 材料組 節(jié)點號（逆時針排列） \n39。 ) 。 fprintf(outfile, 39。 1 2 3 4 \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 for i=1:size(gElement(:,1)) fprintf(outfile, 39。%4d , %4d , %4d , %4d , %4d , %4d\n39。,... i, gElement(i,5), gElement(i,1:4)) 。 end fprintf( outfile,39。==========================================\n\n\n39。 ) 。 fprintf(outfile, 39。 表3 材料性質 \n39。 ) 。 fprintf(outfile, 39。==========================================\n39。 ) 。 fprintf(outfile, 39。材料號 彈性模量 泊松比 比重 \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 for i=1:size(gMaterial(:,1)) fprintf(outfile, 39。%4d % % %\n39。,... i,gMaterial(i,1),gMaterial(i,2),gMaterial(i,3)) 。 end fprintf( outfile,39。==========================================\n\n\n39。 ) 。 fprintf(outfile, 39。 表4 邊界條件 \n39。 ) 。 fprintf(outfile, 39。==========================================\n39。 ) 。 fprintf(outfile, 39。 序號 結點號 約束自由度 給定值 \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 for i=1:size(gBC1(:,1)) fprintf(outfile, 39。%4d %4d %4d %4d\n39。,... i,gBC1(i,1),gBC1(i,2),gBC1(i,3)) 。 end fprintf(outfile, 39。==========================================\n\n\n39。 ) 。 fprintf(outfile, 39。 表5 分布力 \n39。 ) 。 fprintf(outfile, 39。==========================================\n39。 ) 。 fprintf(outfile, 39。邊界 結點1 結點2 均布力1 均布力2 \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 for i=1:size(gDF(:,1)) fprintf(outfile, 39。%4d %4d %4d % %\n39。, ... i,gDF(i,1),gDF(i,2),gDF(i,3),gDF(i,4)) 。 end fprintf( outfile,39。==========================================\n\n\n39。 ) 。 fprintf( outfile,39。 表6 xx向節(jié)點應力極值 \n39。 ) 。 fprintf(outfile, 39。=======================================================\n39。 ) 。 fprintf(outfile, 39。極大點 Stress_xx_max 極小點 Stress_xx_min \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 fprintf(outfile,39。%4d % %4d %\n39。, ... N_xx_max,xx_max,N_xx_min,xx_min) 。 fprintf(outfile, 39。=======================================================\n\n\n39。) 。 fprintf( outfile,39。 表7 yy向節(jié)點應力極值 \n39。 ) 。 fprintf(outfile, 39。=======================================================\n39。 ) 。 fprintf(outfile, 39。極大點 Stress_yy_max 極小點 Stress_yy_min \n39。 ) 。 fprintf(outfile, 39。\n39。 ) 。 fprintf(outfile,39。%4d % %4d %\n39。, ... N_yy_max,yy_max,N_yy_min,yy_min) 。 fprintf(outfile, 39。=======================================================\n\n\n39。) 。 fprintf( outfile,39。 表8 xy向節(jié)點應力極值 \n39。 ) 。 fprintf(outfile, 39。==========