4.边界元方法.pptVIP

计算物理 边界元方法 边界积分方程(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 边界元方法 边界积分方程 边界元近似 单一边界下的边界元法 两种介质的边界元方法 √ 边界元方法的特点 基于边界积分方程的近似方法 结点仅分布在区域的边界 以边界积分方程为控制方程，将边界离散插值，转化为代数方程组 未知量的个数少，求解的计算量少 面积分?边界积分 定理：如果 u(x, y) 和 v(x, y) 是定义在平面域 D 上的两个任意函数，并且它们在边界外法线上的导数为 x y O ds ds √ 证明：利用格林公式(略) d 函数 定义： 性质： 平面域： 边界积分方程 v 和 p 是已知的，上式给出 D + G 内任意一点(等式左边)与 G 上某点(等式右边)的 u 和 q 之间的关系 上式称为边界积分方程，是边界元法的基础 √ 边界元方程 当定点 i 在边界 G 上 常数边界离散化 边界积分方程的右边：将边界 G 分成 N 段，以每段中点的 u 和 q 近似该段的函数值 例：泊松方程 当定点 i 不在边界 G 上：利用前面的结果，代入最上公式 √ 对角元 Hii 和 Gii 的计算(定点 i 在 Gi 上) 非对角元 Hij 和 Gij 的计算(定点 i 不在 Gj 上) √ Bi 的计算：将区域 D 划分为有限个三角形单元 例：三角区域的电势和电量 x O y 1 1 泊松方程的混合边界问题 常数边界离散化 对角元 Hii 和 Gii A ① G2 G1 G1 ② ③ 非对角元 Hij 和 Gij ① ② ③ ① ② ③ A ① G2 G1 G1 ② ③ A ① G2 G1 G1 ② ③ A ① G2 G1 G1 ② ③ A ① G2 G1 G1 ② ③ √ Bi ：因为 f = 0，所以 Bi = 0 构造方程 AX = R 内部(三角形中心(1/3, 0))：u = 0.623 √ 例：三角区域的电势和电量 主程序 x O y 1 1 √ 边界离散化和初始化 √ 构造矩阵 H、 G、A 和 R √ 非对角元 Hij 和 Gij √ 求解线性方程 √ 计算内点 作图 √ 二维区域、无自由电荷的泊松方程 G1 G2 GI D1 D2 区域拆分和边界元方程 区域拆分约定 G1 G2 G2I D1 D2 G1I 边界元方程 求解 根据边界条件，未知量移左边，已知量移右：AX = R 编号规则：按逆时针方向，两个区域的编号连接 √ *