微纳系统仿真大作业.docx

微纳系统仿真大作业班级：机硕1302姓名：刘智勇学号：M201370382指导老师：甘志银1、用有限差分法和有线元方法把以下问题变成数值方程，并说明两种方法的异同：?2?(x,y)??0边界条件：?(x,0)???(x,1)??0?(0,y)???(1,y)??0解：（1）有限差分法有限差分法的基本思想是将问题求解域划分为均匀的差分网格，用有限个网格节点代替连续的求解域。基于Taylor级数展开等方法，把描述问题的微分方程中的微分用网格节点上的函数值的差分来代替，从而将微分方程转化为以网格节点上函数值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法?2?(x,y)??0该方程为拉普拉斯方程?2u??uuyy将运算符离散化，则f?(x)的近似值公式为f?(x)??f(x??h)??2f(x)??f(x??h)??O(h2)h2（1-1）由于定义域为0??x??1,0??y??1,将定义域划分为均匀网域，x，y方向的步长都记为h，则?2u??u(x??h,y)??u(x??h,y)??u(x,y??h)??u(x,y??h)??4u(x,y)??O(h2)h2将O(h2)这个高阶无穷小舍去，则（1-2）?2u??u(x??h,y)??u(x??h,y)??u(x,y??h)??u(x,y??h)??4u(x,y)??0h2（1-3）其中在所有的内部网格点(x,y)??(xi,yj),i??2,?,n?1;j??2,?,m?1,精度为O(h2)，用ui,j近似表示u(xi,yj)，则式（1-3）可以表示为u??u??uu??4u?2ui,j????i?1, ji?1,ji,j?1i,j?1i,j ??0h2（1-4）这就是拉普拉斯方程的五点差分公式，将h2消去可得ui?1,j ??ui?1,j ??ui,j?1??ui,j?1??ui,j?1??4ui,j ??0（1-5）根据边界条件?(x,0)???(x,1)???(0,y)???(1,y)??0,将定义域划分为5×5网格，如下图所示：最终得到由9个方程组成的AP=B线性方程组，表示为：??4p1??p2p1??4p2??p3p2??4p3p4p5??0??0??p6??0p1??4p4??p5p2??p4??4p5??p6p7??0??p8??0p3??p5??4p6p4??4p7??p8??p9??0??0p5??p7??4p8??p9??0p6??p8??4p9??0通过高斯消去法可以得到解向量P??[p1, p2,?,p9]（2）有限元方法使用有限元发的计算流程为：A、求解区域离散化；B、构造插值函数形成分段光滑的坐标函数系；C、用Ritz方法求解微分方程。对??x,y??构造函数1???????2????????????x,y????????????????x??????????y??dxdy??首先将整个区域离散为三角形的子区域如图1，三角形微小子区域中的值由三角形节点值的插值结果表示，即图1求解域离散为n三角形??x,y????Ni?i ??Nj?j ??Nk?k其中Ni,Nj,Nk为三角形的节点插值函数，?i,?j,?k为函数??x,y??在节点i,j,k处的函数值。可知每一个三角形子区域中泛函可由节点插值函数和节点函数值表示，那么对于整个求解区域的泛函表达式为,????x,y???????1,?2,???,?n?由变分原理可知，????????????????????0??1??2??n计算这n个方程即可得到整个求解区域的值。（3）有限差分法与有限元法的异同相同的地方：两种方法都使用了离散化的思想，对求解域进行离散。并且都是将求解域划分成有限个网格进行近似求解。不同的地方：有限差分法是利用级数的概念将连续函数离散化，正如高等数学上所学的连续函数用泰勒级数表达一样，网格上的结点就是级数中的一个取值点，这样以级数和的形式求得最终的解，这个解是近似解，其余项就是误差。有限元法是利用插值原理对求域进行近似求解，将求解域划分网格，每个网格看作一个单元进行求解，这样可以得到若干有限个单元的解，这些解的集和构成整体函数的解。就是说每个单元一个解，这些解分布在整个求解域上，构成不同区域解的变化，如力的变化，温度的变化，这样就可以宏观上看到在不同点上不同的值了。其中具体的表现在以下几个方面：1、有限元单元的形状没有限制，可以做到单元形状与边界重合，处理无规边界的问题；2、有限元统一对待区域内的节点与边界上的节点，因此节点的计算精度总体上协调；而有限差分必须分别处理微分方程与定解条件；3、有限元法是各种单元可以混合使用，所以写不出方程也可以求解；差分法划分的网格是规则的，对方程进行离散化，就是用很