发布时间 : 星期日 文章matlaB程序的有限元法解泊松方程更新完毕开始阅读176f49d649649b6648d747f4
NE(1,n1)=NN(i,i); NE(2,n1)=NN(i-1,i); NE(3,n1)=NN(i-1,i-1); n1=n1+1;
NE(1,n1)=NN(i,i); NE(2,n1)=NN(i-1,i+1); NE(3,n1)=NN(i-1,i); n1=n1+1;
NE(1,n1)=NN(i,i); NE(2,n1)=NN(i,i+1); NE(3,n1)=NN(i-1,i+1);
end %下斜边界节点要对应三个单元，上左，左上，左下 n1=n1+1;
NE(1,n1)=NN(Imax+1,Imax+1); NE(2,n1)=NN(Imax,Imax+1);
NE(3,n1)=NN(Imax,Imax);%右顶点的左下 n1=n1+1;
NE(1,n1)=NN(Imax+1,Imax+1); NE(2,n1)=NN(Imax,Imax+2);
NE(3,n1)=NN(Imax,Imax+1);%右顶点的左上 K=0;
for i=Imax:-1:2 K=K+2; n1=n1+1;
NE(1,n1)=NN(i,i+K); NE(2,n1)=NN(i-1,i+K+1); NE(3,n1)=NN(i-1,i+K); end %右上边界的左上
nel=n1;%此时n1值为总的单元个数
六、输出文本文件
*********有限元法求解三角形区域上Possion方程的结果**********
节点编号 坐标分量x 坐标分量y u（x,y）的值
1 0.01333 -0.98667 0.00000 2 0.01333 -0.98000 0.00001 3 0.02667 -0.98000 0.00001
4 0.01333 -0.97333 0.00002 5 0.02667 -0.97333 0.00003 6 0.04000 -0.97333 0.00002 7 0.01333 -0.96667 0.00002 8 0.02667 -0.96667 0.00004 9 0.04000 -0.96667 0.00005 10 0.05333 -0.96667 0.00003 11 0.01333 -0.96000 0.00003 12 0.02667 -0.96000 0.00006 13 0.04000 14 0.05333 15 0.06667 16 0.01333 17 0.02667 18 0.04000 19 0.05333 20 0.06667 21 0.08000 22 0.01333 23 0.02667 24 0.04000 25 0.05333 26 0.06667 27 0.08000 28 0.09333 29 0.01333 30 0.02667 31 0.04000 32 0.05333 33 0.06667 34 0.08000 35 0.09333 36 0.10667 37 0.01333 38 0.02667 39 0.04000 40 0.05333 41 0.06667 42 0.08000 43 0.09333 44 0.10667 45 0.12000 46 0.01333 47 0.02667 48 0.04000 -0.96000 0.00007 -0.96000 0.00007 -0.96000 0.00005 -0.95333 0.00004 -0.95333 0.00008 -0.95333 0.00010 -0.95333 0.00011 -0.95333 0.00010 -0.95333 0.00007 -0.94667 0.00005 -0.94667 0.00010 -0.94667 0.00014 -0.94667 0.00016 -0.94667 0.00016 -0.94667 0.00013 -0.94667 0.00008 -0.94000 0.00006 -0.94000 0.00012 -0.94000 0.00017 -0.94000 0.00020 -0.94000 0.00022 -0.94000 0.00021 -0.94000 0.00017 -0.94000 0.00010 -0.93333 0.00007 -0.93333 0.00015 -0.93333 0.00021 -0.93333 0.00025 -0.93333 0.00028 -0.93333 0.00028 -0.93333 0.00026 -0.93333 0.00021 -0.93333 0.00012 -0.92667 0.00009 -0.92667 0.00017 -0.92667 0.00024
49 0.05333 -0.92667 0.00030 50 0.06667 -0.92667 0.00034 51 0.08000 -0.92667 0.00036 52 0.09333 -0.92667 0.00036 53 0.10667 -0.92667 0.00032 54 0.12000 -0.92667 0.00025 55 0.13333 -0.92667 0.00015 56 0.01333 -0.92000 0.00010 57 0.02667 -0.92000 0.00020 58 0.04000 59 0.05333 60 0.06667 61 0.08000 62 0.09333 63 0.10667 64 0.12000 65 0.13333 66 0.14667 67 0.01333 68 0.02667 69 0.04000 70 0.05333 71 0.06667 72 0.08000 73 0.09333 74 0.10667 75 0.12000 76 0.13333 77 0.14667 78 0.16000 79 0.01333 80 0.02667 81 0.04000 82 0.05333 83 0.06667 84 0.08000 85 0.09333 86 0.10667 87 0.12000 88 0.13333 89 0.14667 90 0.16000 91 0.17333 92 0.01333 93 0.02667 -0.92000 0.00028 -0.92000 0.00036 -0.92000 0.00041 -0.92000 0.00045 -0.92000 0.00046 -0.92000 0.00044 -0.92000 0.00038 -0.92000 0.00030 -0.92000 0.00017 -0.91333 0.00011 -0.91333 0.00022 -0.91333 0.00032 -0.91333 0.00041 -0.91333 0.00048 -0.91333 0.00053 -0.91333 0.00056 -0.91333 0.00056 -0.91333 0.00052 -0.91333 0.00045 -0.91333 0.00034 -0.91333 0.00019 -0.90667 0.00013 -0.90667 0.00025 -0.90667 0.00037 -0.90667 0.00047 -0.90667 0.00055 -0.90667 0.00062 -0.90667 0.00066 -0.90667 0.00068 -0.90667 0.00066 -0.90667 0.00061 -0.90667 0.00052 -0.90667 0.00039 -0.90667 0.00022 -0.90000 0.00014 -0.90000 0.00028
94 0.04000 -0.90000 0.00041 95 0.05333 -0.90000 0.00053 96 0.06667 -0.90000 0.00063 97 0.08000 -0.90000 0.00071 98 0.09333 -0.90000 0.00077 99 0.10667 -0.90000 0.00080 100 0.12000 -0.90000 0.00080 101 0.13333 -0.90000 0.00077 102 0.14667 -0.90000 0.00070 103 0.16000 104 0.17333 105 0.18667 106 0.01333 107 0.02667 108 0.04000 109 0.05333 110 0.06667 111 0.08000 112 0.09333 113 0.10667 114 0.12000 115 0.13333 116 0.14667 117 0.16000 118 0.17333 119 0.18667 120 0.20000 121 0.01333 122 0.02667 123 0.04000 124 0.05333 125 0.06667 126 0.08000 127 0.09333 128 0.10667 129 0.12000 130 0.13333 131 0.14667 132 0.16000 133 0.17333 134 0.18667 135 0.20000 136 0.21333 137 0.01333 138 0.02667 -0.90000 0.00059 -0.90000 0.00044 -0.90000 0.00024 -0.89333 0.00016 -0.89333 0.00031 -0.89333 0.00045 -0.89333 0.00059 -0.89333 0.00071 -0.89333 0.00080 -0.89333 0.00088 -0.89333 0.00093 -0.89333 0.00095 -0.89333 0.00093 -0.89333 0.00089 -0.89333 0.00080 -0.89333 0.00067 -0.89333 0.00049 -0.89333 0.00027 -0.88667 0.00017 -0.88667 0.00034 -0.88667 0.00050 -0.88667 0.00065 -0.88667 0.00078 -0.88667 0.00090 -0.88667 0.00099 -0.88667 0.00106 -0.88667 0.00110 -0.88667 0.00110 -0.88667 0.00107 -0.88667 0.00101 -0.88667 0.00090 -0.88667 0.00074 -0.88667 0.00055 -0.88667 0.00030 -0.88000 0.00019 -0.88000 0.00037