有限元和边界元方法.pptVIP

计算物理有限元和边界元方法物理问题的变分原理(1/3)物理问题的变分原理(2/3)物理问题的变分原理(3/3)泊松方程的有限元方法(1/11)泊松方程的有限元方法(2/11)泊松方程的有限元方法(3/11)泊松方程的有限元方法(4/11)泊松方程的有限元方法(5/11)泊松方程的有限元方法(6/11)泊松方程的有限元方法(7/11)泊松方程的有限元方法(8/11)泊松方程的有限元方法(9/11)泊松方程的有限元方法(10/11)泊松方程的有限元方法(11/11)扩散方程的有限元方法(1/1)波动方程的有限元方法(1/1)边界积分方程(1/2)边界积分方程(2/2)边界元近似(1/4)边界元近似(2/4)边界元近似(3/4)边界元近似(4/4)单一边界下的边界元法(1/6)单一边界下的边界元法(2/6)单一边界下的边界元法(3/6)单一边界下的边界元法(4/6)单一边界下的边界元法(5/6)单一边界下的边界元法(6/6)两种介质的边界元方法(1/1)作业**3/lesson/ComputationalPhysics有限元和边界元方法物理问题的变分原理泊松方程的有限元方法扩散方程的有限元方法波动方程的有限元方法边界积分方程边界元近似单一边界下的边界元法两种介质的边界元方法√有限元方法基于变分原理的离散化方法——部分逼近地离散化划分整体区域为有限个基本块(单元)在单元上插值逼近，得到结构简单的函数集(有限元空间，是泛函J(y)的定义域的子集)将边值问题转化为泛函的极值问题在有限元空间中寻找泛函J(y)的极小值，作为近似解物理中的变分例：力学体系的最小作用量原理体系的特性可以用拉格朗日函数L(q,q,t)描写在时刻t1和t2之间体系按照以下积分取最小值的方式运动(即，运动轨道由泛函的的极小值决定)√例：光学的费马原理光从点A到点B的传播路径是使光程L取极值由dL=0得到几何光学的折射定律和反射定律例：电磁学的麦克斯韦方程组电磁场的拉格朗日函数L是空间积分电磁学的作用量是时间积分运动方程由泛函的的极小值决定(即dS=0)√例：静电场的泊松方程第一类边界条件等价的变分问题为求解泛函的极值问题泛函的求解必须在边界条件下：条件变分问题第二类和第三类边界条件等价的变分问题为求解泛函的极值问题边界条件包含泛函中：自然边界条件√静电场中二维泊松方程的有限元方法G1G2ABDG1阳极G2G1阴极G2例：阴极射线管(如右图)，在两极上(边界G1)的电势u是已知的，在左右两侧(边界G2)的q是已知的。如果管中的自由电荷密度分布r(x,y)已知，则以上的泊松方程等价为求解以下泛函J(u)的极值问题√有限元方法的具体步骤划分整体区域为有限个单元和编号划分要点三角形的顶点相连避免钝角(因引入较大误差)每个三角形不跨越不同的介质每个三角形最多只有一条边在G2上(方便计算)三角形覆盖尽量多的区域编号约定三角形单元的编号：e=①,②,③,…顶点的编号：按逆时针为1,2,3顶点的坐标：(x1,y1),(x2,y2),(x3,y3)单元的泛函：Je(u)③①②e123e整体的泛函：√3(x3,y3)e1(x1,y1)2(x2,y2)构造线性插值函数假设每个单元内u(x,y)是x和y的线性函数每个的三个基函数u(x,y)的插值表达式中，a,b,c,d可由三角形的顶点坐标确定，只剩余u1,u2,u3未知√建立单元的矩阵单元的泛函123eG2第一项积分与单元刚度矩阵(zij)第二项积分与单元矩阵(rfj)第三项积分与单元矩阵(rqj)√建立顶点和结点的(V-n)对应关系单元编号：有一条边在G2上且q?0的单元编号为1,2,…,e1，其余的单元编号为e1+1,e1+2,…,e0G2G1③①②e1e1+1e0顶点编号：用V(e,i)表示，逆时针方向，2和3在G2上结点编号：内部和G2上的结点编号为1,2,…,n1，G1上的结点编号为n1+1,n1+2,…,n0建立顶点和结点的对应关系：V(e,i)=n集成泛函和建立方程泛函的离散化K为总体刚度矩阵，由单元刚度矩阵(zij)合成Rf由单元矩阵(rfj)合成，Rq