NURBS在电场计算中研究.docVIP

NURBS在电场计算中研究 　　摘要：在运用边界元方法计算开域场中的电场时，为提高传统一阶边界元法的精度问题，提出NURBS曲面四边形边界元法。此方法是对计算模型进行二阶四边形剖分，采用剖分信息构造出NURBS曲面参数方程，然后用面积比值法定义曲面上顶点的形状函数，曲面四边形的四个顶点作为求解点，计算NURBS曲面四边形单元上的积分值。 　　关键词：NURBS；电场强度；提高精度；电场计算；边界元法 文献标识码：A 　　中图分类号：TM151 文章编号：1009-2374（2016）29-0018-04 DOI：10.13535/j.cnki.11-4406/n.2016.29.008 　　边界元法在求解电磁数值计算中，是利用离散技术将求解的边界面离散为边界单元，将边界问题等价转化为边界积分方程问题，同时把边界积分方程离散为代数方程组，最后利用数值方法计算方程组，得到边界积分方程的解。 　　本文提出NURBS曲面四边形边界元方法。利用二阶剖分单元的节点信息，构造出NURBS曲面四边形单元，以曲面四边形单元的顶点为求解点，计算NURBS曲面四边形单元上的积分，最后形成NURBS曲面四边形边界元方法求解三维电场边值问题。 　　1 NURBS曲面边界元的基本原理 　　1.1 有理B样条基函数 　　NURBS（Non-Uniform Rational B-Spline）方法，是在B样条方法基础上发展起来的，为了能进一步精确地描述有理数，推广得到有理多项式基函数。将B样条基函数推广到有理的基函数，就得到了NURBS方法的基函数即有理B样条基函数，而且NURBS曲面在中心投影下得到B样条曲面，所以研究NURBS方法就必须了解B样条方法基础。 　　B样条函数的概念最初是由Schoenberg于1946年提出的。由de Boor-Cox公式得出提出B样条的递推公式为： 　　式中：称为相应于序列的次，即阶B样条基函数。其中，下标表示序号，下标表示次数。可以看出想要得到第个次B样条基函数，就得共个节点序列，而且的支撑区间为。由该定义出发，可以导出任意阶次的B样条基函数。文章中是采用高斯积分方法求解数值积分，由于高斯积分法对应的积分区间为。而在区间上B样条基函数为均匀的基函数，根据B样条基函数的性质，次B样条基函数只有个基函数在上非零，而且在上都为多项式函数。取，则由式（1）可得均匀2次B样条基函数表达式为： 　　有理B样条基函数作为对B样条基函数的推广，一个基函数变为有理分式首先需要相应的一个多项式作为分母，然而考虑到相同的分母将会给有理分式之间的计算带来简便。再注意到B样条基函数具有的“权”性，即在区间的非负性和单位分解性，因此可以采取下面的方法：给每一个B样条基函数乘以一个“权因子”，再用它们的“和”作为统一分母进行“平均”，这样就得到了次有理B样条基函数为： 　　1.2 NURBS曲面 　　NURBS曲面是非有理张量积型性B样条曲面的有理推广，类似于NURBS曲线定义，一张次NURBS曲面有理式可以表示为： 　　其中是双变量的有理基函数： 　　式中：为控制点；为与控制点相应的权因子；和分别为参数向次和参数向次的定义在和规范有理B样条基函数。 　　根据式（3），当时，得到广泛应用的双二次NURBS曲面的数学模型为： 　　上式中是由式（2）推导出，为的转置矩阵。 　　根据式（7），要得出NURBS曲面的参数方程，必须就算出控制点P11，P12，…，P33。根据有理式表示的NURBS曲面参数方程还可以表示为，其中为曲面数据点列，维数为；为有理基函数，维数为，为方阵，则逆矩阵存在；是控制点，维数为。在已知曲面上的数据点列，假定权因子时，经过矩阵求逆变换可得到双二次NURBS曲面控制点矩阵求逆变换可得。 　　1.3 四边形单元上的形状函数 　　本文是以四边形单元的四个顶点为求解点，运用面积比值法构造曲面四边形顶点上的形状函数。曲面单元如图1所示，曲面单元上任意一点处的函数值由4个顶点节点i、j、k、l上的函数值和点相对于4个顶点节点的形状函数插值得到。在图1中，点处的两条参数等值线将曲面单元分成、、、共4部分，其中每一部分的面积都可以通过曲面参数方程得到。设曲面单元总的面积为。则点处关于4个顶点的形状函数定义为： 　　1.4 边界元法 　　三维静电场光滑边界的积分方程为： 　　式中：为边界上的场点；为边界上的源点；为场点电位；为源点电位；为源点到场点的距离，；为边界外法线方向；S为求解区域的边界。 　　2 NURBS曲面边界元算例分析 　　2.1 算例一 　　如图2是高为2、半径为1的四分之一圆柱面。运用面积计算可知其表面积为。将其剖分一个单元，得出节点信息后，导入一阶平面四边形边界元程序计算其表