摘要

用有限元法数值求解时,定义在流形曲面上的偏微分方程的数值解精度会因为传统多边形单元的几何逼近误差而严重降低,为此提出基于有理Bernstein多项式的几何精确有限元法.首先插入重复节点从NURBS曲面直接生成有理Bézier单元,这一过程保持原有几何不变;然后通过Galerkin法建立参数曲面上包含Laplace-Beltrami微分算子的二阶椭圆偏微分方程的等效弱形式;针对Bernstein基函数的非插值性,通过配点法施加Dirichlet类型的边界约束,得到最优收敛的离散格式.数值算例结果表明,该方法能有效地减少网格离散误差,提高分析结果精度.