基于电磁数值计算的几何模型细小结构消除转化方法

摘要

本发明涉及基于电磁数值计算的几何模型细小结构消除转化方法，步骤如下：步骤101：将模型中信息读入内存；步骤102：设变量i＝0；步骤103：从内存中读取模型中的第一个面片；步骤104：计算出该面片的顶点数n；步骤106：将面片的四条边读出并按长度排序；步骤107：判断a是否小于m；步骤108：令j＝j+1；步骤109：对短边的处理；步骤110：j是否小于n；步骤111：获取次短边的信息；步骤112：i是否小于num1；步骤113：从内存中读取模型中的下一个面片；步骤114：处理折面以及天线；步骤115：从内存中读取第一个面片；步骤116：判断该面的顶点数n；步骤117：从模型中删除该面片；步骤118：令i＝i+1；步骤119：判断i是否小于num1；步骤120：从内存中读取模型中的下一个面片。

权利要求

1.基于电磁数值计算的几何模型细小结构消除转化方法，其特征是：包括如下步骤：步骤101：将模型中信息读入内存，并计算出模型中的面片数num1和天线数num2；

步骤102：设变量i=0；

步骤103：从内存中读取模型中的第一个面片；

步骤104：计算出该面片的顶点数n，并设变量j=0；

步骤105：i=i+1；

步骤106：将面片的四条边读出并按长度排序，获取最短边及其边长a；

步骤107：判断a是否小于m，m=kλ，λ为波长，如果a小于m则转到步骤108；如果a大于或等于m，则转到步骤112；

步骤108：令j=j+1；

步骤109：对短边的处理；包括如下步骤：步骤201：读取该短边的两个端点A与B的信息；

步骤202：找到模型中点A的相关面片；

步骤203：将相关面片中与A点重合的顶点坐标值修改为B点的坐标值，处理短边结束；

步骤110：j是否小于n，如果是，则转到步骤111；如果不是，则转到步骤112；

步骤111：获取次短边的信息及其边长a并转到步骤107；

步骤112：i是否小于num1，如果是，则转到步骤113；如果不是，则转到步骤114；

步骤113：从内存中读取模型中的下一个面片，转到步骤104；

步骤114：处理折面以及天线；包括如下步骤：步骤301：从内存中读取模型中的第一个面片并设i=0；

步骤302：判断该面片的顶点数n是否为4，如果是，则转到步骤303；如果不是，则转到步骤305；

步骤303：判断面片的四个顶点是否共面，如果共面，则转到步骤305；如果不共面，则转到步骤304；

步骤304：用较短的对角线将该面片划分为两个三角形；

步骤305：判断i是否小于num1，如果是，则转到步骤306；如果不是，则转到步骤307；

步骤306：从内存中读取下一个面片，转到步骤302；

步骤307：从内存中读取模型中的第一根天线并设i=0；

步骤308：令i=i+1；

步骤309：判断该天线是否落点在模型的表面上，如果是，则转到步骤313；如果不是，则转到步骤310；

步骤310：找到与该天线所在直线相交且与其落点距离最近的面；

步骤311：将该天线所在直线与该面求交，求得的交点即为天线的新落点；

步骤312：新的落点与原来的落点相比坐标的增量为Δ，则天线原来的顶点坐标加上Δ即为该天线新的顶点坐标；

步骤313：判断i是否小于num2，如果是，则转到步骤314；如果不是，则对折面以及天线的处理结束；

步骤314：从内存中读取下一根天线的信息，转到步骤308；

步骤115：从内存中读取模型中的第一个面片并令i=0；

步骤116：判断该面片的顶点数n是否小于3，如果是，则转到步骤117；如果不是，则转到步骤119；

步骤117：从模型中删除该面片；

步骤118：令i=i+1；

步骤119：判断i是否小于num1，如果是，则转到步骤120；如果不是，算法结束；

步骤120：从内存中读取模型中的下一个面片，返回步骤116。

2.根据权利要求1所述的基于电磁数值计算的几何模型细小结构消除转化方法，其特征是：所述的找到模型中点A的相关面片是对所有以点A为顶点的面片。

3.根据权利要求1所述的基于电磁数值计算的几何模型细小结构消除转化方法，其特征是：所述的步骤303四点共面的判断方法为：设面片的四个顶点分别为A、B、C、D，可以组成三个向量 和 如果满足则可确定A、B、C、D四点共面。

说明书

基于电磁数值计算的几何模型细小结构消除转化方法

技术领域

[0001] 本发明属于计算机载体模型建模领域，基于电磁数值计算的几何模型细小结构消除转化方法。

背景技术

[0002] 车载、机载、舰载、船载和星载通信系统均可归结为载体通信系统。载体以运载工具不同可分为陆基的车体、空基的机体、海基的舰(船)体、天基的星体，载体为电磁学上的金属导体，其上安装的天线有短波天线、超短波天线、超高频天线、全球通天线、卫星天线、微波天线、集群车载天线、无线接入天线等。载体天线通信系统在抗震救灾、抗洪抢险、突发事件应急以及军事指挥和作战方面有着广泛的应用，已成为重要的民用和军用装备。

[0003] 载体天线通信系统的载体上安装有少则几部天线，多则十几部天线，这些天线分布在较狭小的空间，天线之间往往存在严重的邻道干扰、谐波干扰、互调干扰等，载体结构和天线位置也严重影响天线间的隔离度及天线的方向图，这些电磁兼容方面的问题轻者会影响系统的性能，重者会导致通信的中断。采用电磁兼容仿真设计和预测分析软件是解决这些问题的有效途径，并能缩短产品研发周期、降低产品成本、优化产品设计。电磁兼容仿真设计和预测分析软件的核心是电磁数值计算以及为数值计算提供离散数据的网格划分。现阶段主要的电磁数值计算方法为有限元法、时域有限差分法、矩量法等。有限元法和时域有限差分法是体求解方法，而对于载体天线结构，由于高频电流的趋肤效应，感应电流分布在载体金属的表层，采用基于表面电流技术的矩量法进行求解，其适应性要好于有限元法和时域有限差分法，而且矩量法的计算结果精度比较高、可求解复杂形状目标的电磁计算，并可与快速多极子方法、快速傅里叶变换、波形渐进预估技术、时域方法、有限元方法等相结合，以解决高频计算、电大尺寸目标计算、宽频带计算、非均匀介质计算等，因而它获得了越来越广泛的应用，成为计算电磁学中最为重要的方法。

[0004] 对于载体天线结构中的鞭天线而言，由于其天线较细，可认为电流只沿轴线分布；载体可用封闭的金属导体面近似。在求解这类问题时，需要在导线上建立电场积分方程，在载体面上建立磁场积分方程，且在天线与载体面相连接的区域，其积分域包括直线段和导体面，也要通过建立电场和磁场混合积分方程来求解。求解这些积分方程可得到天线上的线电流和载体面的面电流，进而可计算出天线的隔离度以及电场分布。目前求解这些积分方程的数值方法主要是矩量法。矩量法求解的前提是载体面的网格划分和天线的网格划分，天线的网格划分简单易行，而载体面的网格划分由于载体结构的复杂往往不能得到理想的结果。

[0005] 对于单元形状为三角形的网格划分方法而言，在符合数值计算要求的前提下，划分结果中狭长三角形要尽量少，因为狭长三角形会使数值计算的结果发生畸变，使计算结果与实际相差很远。目前，已经有不少成熟的网格划分方法可以自动调整网格单元的质量，使划分结果尽量理想。但是在某些情况下，如当载体中出现一些细小结构时，划分结果中就会不可避免得出现狭长三角形，从而严重地影响计算结果的准确性。

[0006] 为了得到更好的网格划分结果，人们往往追求更完善的网格划分方法，而很少会从模型简化的角度去解决问题。在电磁仿真中，载体模型的建立是为了利用计算机得到与实际相符合的仿真结果，以解决实际问题，这样既可以缩短工作时间，又可以节约成本。而为了得到更加符合实际的计算结果而一味地追求模型的精细，不但不能达到本来的目的，还会使结果发生更大的偏差。因此，我们可以在适当的范围内对载体模型进行简化，从而解决网格划分中出现狭长三角形的问题。

[0007] 目前对于三维网格模型的划分方法有很多种。论文“基于多边形顶点法矢量网络模型简化算法”(中国图像图形学报，第七卷第六期，2002年6月，601-605，作者：周石琳等)中提出了一种基于多边形顶点法矢量网络模型简化算法；论文“基于特征的三维模型简化算法研究”(系统仿真学报，第19卷第11期，2007年6月，2434-2414，作者：李建军等)提出了一种半边折叠算法；论文“一种基于区域分割的几何模型简化方法”(计算机学报，第29卷第10期，2006年10月，1834-1842，作者：全红艳等)提出了先对模型进行区域划分，然后各个区域按照三角形数目的比例进行简化的方法；论文“基于三角形网格模型简化的研究”(计算机工程与科学，第32卷第12期，2010年，69-72，作者：闫涛等)设计了一种基于三角网格删除的简化算法；论文“基于相似性的点模型简化算法”(浙江大学学报(工学版)，第43卷第3期，2009年3月，448-454，作者：王仁芳等)提出了一种基于相似性的曲率自适应点模型简化算法。这些方法都是对于三维网格模型的简化方法但是其初衷都是在保持原模型形状的前提下尽量减少模型的数据量，即减少网格数。而对于用于电磁计算的网格模型的简化方法还有待研究。

[0008] 目前，还没有专门针对用于电磁计算的网格模型提出的成熟的模型简化方法。在实际仿真中，研究人员多是根据经验，针对某个模型提出简化方案，这要求仿真软件的使用者既要非常了解软件中数值计算的方法原理，又要熟悉软件的建模方法，而且在待仿真的模型较多时，模型简化的工作量就会变得非常大。

发明内容

[0009] 本发明的目的在于克服上述现有的模型中存在细小结构所引起的问题，提供一种基于电磁数值计算的几何模型细小结构消除转化方法，通过计算机对模型进行简化处理。根据天线工作频率的大小，自动识别并简化细小结构，增加仿真结果的准确性，从而减少工作量，节约时间。

[0010] 本发明的目的是这样实现的：基于电磁数值计算的几何模型细小结构消除转化方法，其特征是：包括如下步骤：

[0011] 步骤101：将模型中信息读入内存，并计算出模型中的面片数num1和天线数num2；

[0012] 步骤102：设变量i＝0；

[0013] 步骤103：从内存中读取模型中的第一个面片；

[0014] 步骤104：计算出该面片的顶点数n，并设变量j＝0；

[0015] 步骤105：i＝i+1；

[0016] 步骤106：将面片的四条边读出并按长度排序，获取最短边及其边长a；

[0017] 步骤107：判断a是否小于m(m＝kλ)，如果a小于m则转到步骤108；如果a大于或等于m则转到步骤112；

[0018] 步骤108：令j＝j+1；

[0019] 步骤109：对短边的处理；

[0020] 步骤110：j是否小于n，如果是，则转到步骤111；如果不是，则转到步骤112；

[0021] 步骤111：获取次短边的信息及其边长a并转到步骤107；

[0022] 步骤112：i是否小于num1，如果是，则转到步骤113，如果不是，则转到步骤114；

[0023] 步骤113：从内存中读取模型中的下一个面片，转到步骤104；

[0024] 步骤114：处理折面以及天线；

[0025] 步骤115：从内存中读取模型文件中的第一个面片并令i＝0；

[0026] 步骤116：判断该面的顶点数n是否小于3，如果是，则转到步骤117；如果不是，则转到步骤119；

[0027] 步骤117：从模型中删除该面片；

[0028] 步骤118：令i＝i+1；

[0029] 步骤119：判断i是否小于num1，如果是，则转到步骤120；如果不是，算法结束；

[0030] 步骤120：从内存中读取模型中的下一个面片，返回步骤116。

[0031] 所述的步骤109对短边的处理；包括如下步骤：

[0032] 步骤201：读取该短边的两个端点A与B的信息；

[0033] 步骤202：找到模型中点A的相关面片；

[0034] 步骤203：将相关面中与A点重合的顶点坐标值修改为B点的坐标值，处理短边结束。

[0035] 所述的找到模型中点A的相关面片是对所有以点A为顶点的面片。

[0036] 所述的步骤114处理折面以及天线包括如下步骤：

[0037] 步骤301：从内存中读取模型中的第一个面片并设i＝0；

[0038] 步骤302：判断该面的顶点数n是否为4，如果是，则转到步骤303；如果不是，则转到步骤305；

[0039] 步骤303：判断面片的四个顶点是否共面，如果共面，则转到步骤305；如果不共面，则转到步骤304。假设面片的四个顶点分别为A、B、C、D，可以组成三个向量 和如果 则可确定A、B、C、D四点共面，而在实际计算中，由于计算机中储存的数据所保留的有效数字位数是有限的，计算后的结果会有误差，因此本发明实际用的四点共面的判断公式为：

[0040]

[0041] 步骤304：用较短的对角线将该面片划分为两个三角形；

[0042] 步骤305：判断i是否小于num1，如果是，则转到步骤306；如果不是，则转到步骤307；

[0043] 步骤306：从内存中读取下一个面片，转到步骤302；

[0044] 步骤307：从内存中读取模型中的第一根天线并设i＝0；

[0045] 步骤308：令i＝i+1；

[0046] 步骤309：判断该天线是否落点在模型的表面上，如果是，则转到步骤313；如果不是，则转到步骤310；

[0047] 步骤310：找到与该天线所在直线相交且与其落点距离最近的面；

[0048] 步骤311：将该天线所在直线与该面求交，求得的交点即为天线的新落点；

[0049] 步骤312：新的落点与原来的落点相比坐标的增量为Δ，则天线原来的顶点坐标加上Δ即为该天线新的顶点坐标；

[0050] 步骤313：判断i是否小于num2，如果是，转到步骤314；如果不是，则对折面以及天线的处理结束；

[0051] 步骤314：从内存中读取下一根天线的信息，转到步骤308。

[0052] 本发明对表面模型的细小结构进行简化处理，从而得到便于进行电磁计算的模型，有如下优点：

[0053] (1)本发明对载体模型中的细小结构进行了简化处理，使对其进行网格划分的结果更加理想，使电磁计算的结果更加准确；

[0054] (2)本发明对载体模型中细小结构的识别以及简化处理都是由程序自动完成的，提高了工作的效率以及准确性。

附图说明

[0055] 图1是本发明的总流程图；

[0056] 图2是本发明处理短边的流程图；

[0057] 图3是本发明处理折面及天线的流程图；

[0058] 图4是本发明模拟的载体模型示例图；

[0059] 图5是经过本发明简化后的载体模型示意图。

具体实施方式

[0060] 本发明重点研究载体模型中细小结构的简化处理。本发明对载体模型的要求是：载体模型必须为表面封闭的几何体，其表面由顶点对顶点的三角形或四边形平面片围成，且三角形和四边形平面片的顶点顺序为逆时针方向。

[0061] 以下结合附图对本发明的模型简化方法作进一步详细描述：

[0062] 参照图1，本发明的模型简化包括如下步骤：

[0063] 步骤101：将模型中信息读入内存，并计算出模型中的面片数num1和天线数num2；

[0064] 步骤102：设变量i＝0；

[0065] 步骤103：从内存中读取模型中的第一个面片；

[0066] 步骤104：计算出该面片的顶点数n，并设变量j＝0；

[0067] 步骤105：i＝i+1；

[0068] 步骤106：将面片的四条边读出并按长度排序，获取最短边及其边长a；

[0069] 步骤107：判断a是否小于m(m＝kλ)，如果a小于m则转到步骤108；如果a大于或等于m则转到步骤112；

[0070] 步骤108：令j＝j+1；

[0071] 步骤109：对短边的处理，参照图2，对短边的处理包括如下步骤：

[0072] 步骤201：读取该短边的两个端点A与B的信息；

[0073] 步骤202：找到模型中点A的相关面片(即所有以点A为顶点的面片)；

[0074] 步骤203：将相关面中与A点重合的顶点坐标值修改为B点的坐标值，处理短边结束。

[0075] 步骤110：j是否小于n，如果是，则转到步骤111；如果不是，则转到步骤112；

[0076] 步骤111：获取次短边的信息及其边长a并转到步骤107；

[0077] 步骤112：i是否小于num1，如果是，则转到步骤113，如果不是，则转到步骤114；

[0078] 步骤113：从内存中读取模型中的下一个面片，转到步骤104；

[0079] 步骤114：处理折面以及天线，参照图3，处理折面以及天线包括如下步骤：

[0080] 步骤301：从内存中读取模型中的第一个面片并设i＝0；

[0081] 步骤302：判断该面的顶点数n是否为4，如果是，则转到步骤303；如果不是，则转到步骤305；

[0082] 步骤303：判断面片的四个顶点是否共面，如果共面，则转到步骤305；如果不共面，则转到步骤304。假设面片的四个顶点分别为A、B、C、D，可以组成三个向量 和如果 则可确定A、B、C、D四点共面，而在实际计算中，由于计算机中储存的数据所保留的有效数字位数是有限的，计算后的结果会有误差，因此本发明实际用的四点共面的判断公式为

[0083]

[0084] 步骤304：用较短的对角线将该面片划分为两个三角形；

[0085] 步骤305：判断i是否小于num1，如果是，则转到步骤306；如果不是，则转到步骤307；

[0086] 步骤306：从内存中读取下一个面片，转到步骤302；

[0087] 步骤307：从内存中读取模型中的第一根天线并设i＝0；

[0088] 步骤308：令i＝i+1；

[0089] 步骤309：判断该天线是否落点在模型的表面上，如果是，则转到步骤313；如果不是，则转到步骤310；

[0090] 步骤310：找到与该天线所在直线相交且与其落点距离最近的面；

[0091] 步骤311：将该天线所在直线与该面求交，求得的交点即为天线的新落点；

[0092] 步骤312：新的落点与原来的落点相比坐标的增量为Δ，则天线原来的顶点坐标加上Δ即为该天线新的顶点坐标；

[0093] 步骤313：判断i是否小于num2，如果是，转到步骤314；如果不是，则对折面以及天线的处理结束；

[0094] 步骤314：从内存中读取下一根天线的信息，转到步骤308。

[0095] 步骤115：从内存中读取模型文件中的第一个面片并令i＝0；

[0096] 步骤116：判断该面的顶点数n是否小于3，如果是，则转到步骤117，如果不是，则转到步骤119；

[0097] 步骤117：从模型中删除该面片；

[0098] 步骤118：令i＝i+1；

[0099] 步骤119：判断i是否小于num1，如果是，则转到步骤120；如果不是，算法结束；

[0100] 步骤120：从内存中读取模型中的下一个面片，返回步骤116。

[0101] 经过以上步骤对于拥有细小结构的载体模型的简化就完成了。

[0102] 仿真实例

[0103] 利用本发明对图4所示的载体模型进行简化处理，该模型包含两根天线，结构a与结构b为两处细小结构。

[0104] 在天线发射频率为100MHz即波长λ＝3米的条件下对该模型进行简化处理，此时kλ为0.017379米。当读取到面片1时(如图4中结构a所示)，判断最短边为边AB，长度为0.010米，小于kλ，因此要对其进行处理；首先找到A点的相关面片，即面片1、面片2、面片3和面片4，将这四个面片上与A点坐标相同的顶点的坐标值改为与B点相同，该边的处理结束；找到面片1的次短边即边CD，其长度为0.010米，也小于kλ，应做同样的处理，而剩下的两个边长度大于kλ，该面片处理结束；对于结构b处的细小结构做同样处理；
然后判断模型中是否有折面以及落点不在模型表面上的天线，发现天线I落点不在模型表面，按照算法将其下移到模型表面；最后判断模型中是否有n＜3的面片，并将其删除，简化结束。

[0105] 图5为简化处理后的模型，可以看出两处细小结构都经过了处理，符合要求。