基于移动相似中心的比例边界有限元方法

在传统的比例边界有限元中,相似中心是固定的,难以用其求解关于偏心域的场问题。本文引入移动相似中心的概念,建立新的比例边界坐标变换,并利用加权余量法将控制方程半弱化为关于径向坐标的二阶常微分方程,引入对偶变量,将其降为系数矩阵为Hamilton矩阵的一阶常微分方程。对Hamilton矩阵进行Schur分解,得到微分方程的通解,代入边值条件可得关于积分常数的代数方程。此方法将比例边界有限元扩展到偏心域的边值问题,同时在径向是半解析的,解的精度高;仅需要离散求解域的一个边界,数据量小;在计算中仅需要对Hamilton矩阵进行Schur分解以及求解关于积分常数的代数方程,运算量少。将偏心环形域静电场边值问题的算例与解析解或其他数值方法计算结果的比较,表明此方法具有精度高、数据量小及运算量小的优点。