matlaB程序的有限元法解泊松方程(共138页).doc

基于matlaB编程的有限元法一、待求问题：泛定方程：边界条件：以（0,-1），（0,1），（1,0）为顶点的三角形区域边界上2、 编程思路及方法1、 给节点和三角形单元编号，并设定节点坐标画出以（0,-1），（0,1），（1,0）为顶点的三角形区域figure1由于积分区域规则，故采用特殊剖分单元，将区域沿水平竖直方向分等份，此时所有单元都是等腰直角三角形，剖分单元个数由自己输入，但竖直方向份数（用Jmax表示）必须是水平方向份数（Imax）的两倍，所以用户只需输入水平方向的份数Imax。采用上述剖分方法，节点位置也比较规则。然后利用循环从区域内部（非边界）的节点开始编号，格式为NN(i,j)=n1，i，j分别表示节点所在列数与行数，并将节点坐标存入相应矩阵X(n1)，Y(n1)。由于区域上下两部分形状不同因此，分两个循环分别编号赋