有限元和边界元方法.ppt

计算物理 有限元和边界元方法 物理问题的变分原理(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 取极值 由 d L = 0 得到几何光学的折射定律和反射定律 例：电磁学的麦克斯韦方程组 电磁场的拉格朗日函数 L 是空间积分 电磁学的作用量是时间积分 运动方程由泛函的的极小值决定(即 d S = 0 ) √ 例：静电场的泊松方程 第一类边界条件 等价的变分问题为求解泛函的极值问题 泛函的求解必须在边界条件下：条件变分问题 第二类和第三类边界条件 等价的变分问题为求解泛函的极值问题 边界条件包含泛函中：自然边界条件 √ 静电场中二维泊松方程的有限元方法 G1 G2 A B D G1 阳极 G2 G1 阴极 G2 例：阴极射线管(如右图) ，在两极上(边界 G1)的电势 u 是已知的，在左右两侧(边界 G2)的 q 是已知的。如果管中的自由电荷密度分布 r (x, y) 已知，则 以上的泊松方程等价为求解以下泛函 J(u) 的极值问题 √ 有限元方法的具体步骤 划分整体区域为有限个单元和编号 划分要点 三角形的顶点相连 避免钝角(因引入较大误差) 每个三角形不跨越不同的介质 每个三角形最多只有一条边在 G2 上(方便计算) 三角形覆盖尽量多的区域 编号约定 三角形单元的编号：e = ①, ②, ③,… 顶点的编号：按逆时针为 1, 2, 3 顶点的坐标：(x1, y1), (x2, y2), (x3, y3) 单元的泛函：Je(u) ③ ① ② e 1 2 3 e 整体的泛函： √ 3 (x3, y3) e 1 (x1, y1) 2 (x2, y2) 构造线性插值函数 假设每个单元内 u(x, y) 是 x 和 y 的线性函数 每个的三个基函数 u(x, y) 的插值表达式中，a, b, c, d 可由三角形的顶点坐标确定，只剩余 u1, u2, u3 未知 √ 建立单元的矩阵 单元的泛函 1 2 3 e G2 第一项积分与单元刚度矩阵 (z i j) 第二项积分与单元矩阵 (r f j) 第三项积分与单元矩阵 (r q j) √ 建立顶点和结点的(V-n)对应关系 单元编号：有一条边在 G2 上且 q?0 的单元编号为 1, 2, …, e1，其余的单元编号为 e1+1, e1+2, …, e0 G2 G1 ③ ① ② e1 e1+1 e0 顶点编号：用 V(e, i) 表示，逆时针方向，2和 3在 G2 上 结点编号：内部和 G2 上的结点编号为 1, 2, …, n1，G1 上的结点编号为 n1+1, n1+2, …, n0 建立顶点和结点的对应关系：V(e, i) = n 集成泛函和建立方程 泛函的离散化 K 为总体刚度矩阵，由单元刚度矩阵 (z i j) 合成 Rf 由单元矩阵 (r f j) 合成，Rq