三维Laplace方程快速多极边界元软件用户手册.pdfVIP

三维 LAPLACE 方程快速多极边界元软件 FMBEM_LAPLACE3D 用户手册 王海涛 wanght@mail.tsinghua.edu.cn 1 关于本手册 本手册包括三维LAPLACE 方程快速多极边界元软件FMBEM_LAPLACE3D 的 基本说明和使用方法，请在使用软件前仔细阅读本手册。 本手册包括： FMBEM_LAPLACE3D 简介 使用方法 附录：前处理转换软件PAT_TRANSLATOR 使用说明 2 FMBEM_LAPLACE3D 简介 FMBEM_LAPLACE3D 是一个基于边界元法的三维大规模数值求解器，可用于 以Laplace 方程作为控制方程的三维物理问题求解，如稳态温度场、静电场等。 该软件的主要特点包括： 软件组成： 计算程序FMBEM_LAPLACE3D.EXE 前处理转换程序PAT_TRANSLATOR.EXE 运行环境：Windows 2000 ，Windows XP 单元类型：三结点三角形常值面单元 1 求解算法： 直接高斯消去 Direct Gauss Elimination 预处理GMRES 迭代 Preconditioned GMRES Iterator 树代码 Tree Code O(NlogN) 自适应快速多极 Adaptive Fast Multipole O(N) 新型自适应快速多极 New Version Adaptive Fast Multipole O(N) 计算功能：常规的直接高斯消去法和预处理GMRES 迭代法需要显式存储系 数矩阵，求解规模较小；树代码算法和快速多极算法隐式存储系数矩阵和迭 代向量信息，可以在单台普通PC 上完成大规模边界元计算，且精度可以控 制。 例：内存 1GB 的 PC，使用直接高斯消去法和预处理 GMRES 迭代法求解的上限约 为 1 万自由度；而使用树代码算法和快速多极算法的求解规模取决于计算精度的要 求和求解速度与存储量的平衡要求，通常能够达到几十万自由度。 3 使用方法 3.1 前处理 用户需要为本软件提供运行必需的参数卡片文件 BEM_DATACARD.DAT 和前 处理文件 BEM_INP.DAT，均是文本格式。BEM_DATACARD.DAT 包括求解器 选择、迭代参数和快速多极算法参数等信息；BEM_INP.DAT 包括三维结构的边 界单元节点坐标信息、单元节点编号信息和单元载荷信息。用户可以手工生成这 两个文件，也可以借助其他工具生成。本软件附带一个小的自动转换程序 PAT_TRANSLATOR.EXE ，该转换程序可以将通用有限元软件MSC.PATRAN 的 前处理文件*.bdf 自动转换成本软件可读的前处理文件BEM_INP.DAT ，具体使用 方法见附录。 3.1.1 BEM_DATACARD 参数卡片文件BEM_DATACARD.DAT 的具体格式如下： SOLVER_TYPE （第一行） 2 ITERATIONS （当SOLVER_TYPE = 2，3 时定义） ERRORS （当SOLVER_TYPE = 2，3 时定义） FMM_METHOD （当SOLVER_TYPE = 3 时定义） TRUNC_M （当SOLVER_TYPE = 3 时定义） TRUNC_E （当FMM_METHOD = 3 时定义） LEAF_MAX_ELE_COUNT （当SOLVER_TYPE