CN103246827B - 复杂外形金属目标电磁散射的无网格仿真方法 - Google Patents
复杂外形金属目标电磁散射的无网格仿真方法 Download PDFInfo
- Publication number
- CN103246827B CN103246827B CN201310216541.3A CN201310216541A CN103246827B CN 103246827 B CN103246827 B CN 103246827B CN 201310216541 A CN201310216541 A CN 201310216541A CN 103246827 B CN103246827 B CN 103246827B
- Authority
- CN
- China
- Prior art keywords
- node
- edge
- point
- linked list
- triangle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002184 metal Substances 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000005684 electric field Effects 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000004088 simulation Methods 0.000 claims description 18
- 238000001914 filtration Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 2
- 238000005056 compaction Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 9
- 230000014509 gene expression Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000004613 tight binding model Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种复杂外形金属目标电磁散射的无网格仿真方法,步骤如下:节点局部积分域的建立:根据复杂外形金属目标表面的离散节点的分布信息,通过区域增长算法,构造出每个节点周围各自的局部积分域;建立金属目标电场积分方程;展开离散节点处的面电流,构造出每个局部积分域中心位置处的面电流;对电场积分方程采用点匹配进行测试得到阻抗矩阵方程;求解矩阵方程,得到电流系数,根据电流系数确定雷达散射截面积。本发明不依赖于目标表面的模型化网格剖分,仅需事先知道目标表面离散节点的分布信息,便可对其进行快速的电磁散射仿真,其实现过程灵活自由,具有很强的实际工程应用价值。
Description
一技术领域
本发明属于电磁散射特性仿真技术领域,特别是一种复杂外形金属目标电磁散射的无网格仿真方法。
二背景技术
随着军工领域的发展,越来越多具有复杂外形的金属目标迫切需要对其进行快速的电磁散射仿真。传统的方法分析金属目标电磁散射问题的时候,都需要事先采用商业软件对金属目标进行结构化网格建模,以获得金属目标表面的网格剖分信息。
虽然此传统方法已在学术领域应用普遍,但是在实际工程中处理具有复杂外形结构的金属目标时并不适用。因为对于复杂外形的金属目标,对其进行高质量的网格生成往往是十分困难的,尤其是当金属目标出现剧烈形变的时候,很难获得满足要求的网格信息,并且传统意义上的网格面片还必须成对出现,所以当复杂外形的金属目标中出现三面共线的情况时,便无法对其进行所需的网格生成,从而不能进行金属目标的电磁散射仿真分析。
三发明内容
本发明的目的在于提供一种稳定、高效的复杂外形金属目标电磁散射的无网格仿真方法,该方法不依赖于金属目标表面的模型化网格剖分,实现过程灵活自由,具有很强的实际工程应用价值。
实现本发明目的的技术解决方案为:一种复杂外形金属目标电磁散射的无网格仿真方法,其特征在于,步骤如下:
第1步,节点局部积分域的建立:根据复杂外形金属目标表面的离散节点的分布信息构造出每个节点周围各自的局部积分域;
第2步,建立金属目标电场积分方程;
第3步,展开离散节点处的面电流,构造出每个局部积分域中心位置处的面电流;
第4步,对第2步中的电场积分方程采用点匹配进行测试得到阻抗矩阵方程;
第5步,求解矩阵方程,得到电流系数,根据电流系数确定雷达散射截面积。
本发明与现有技术相比,其显著优点:(1)工程中适用范围更广泛:因为本方法可以处理各种实际中出现的复杂外形金属目标,对目标模型没有严格的网格化要求,所以本方法理论上可以分析几乎所有的金属目标;(2)方法实现过程灵活简单:由于本方法只需要实际金属目标表面的离散节点的分布信息,不需要事先知道目标表面的网格分布情况,前处理过程可以做到自适应性构造局部积分域,所以本方法更加自由灵活简单;(3)形成矩阵方程性态较好:本方法中的方程属于第二类积分方程,而第二类积分方程具有良好的迭代求解性态,当采用迭代求解时,可以快速收敛至所需要的计算精度。
下面结合附图对本发明作进一步详细描述。
四附图说明
图1是本发明金属目标表面离散节点分布示意图。
图2是本发明区域增长的基本概念示意图。
图3是本发明最佳匹配点为外点时区域增长示意图。
图4是本发明最佳匹配点为当前边的起点邻接点时区域增长示意图。
图5是本发明最佳匹配点为当前边的终点邻接点时区域增长示意图。
图6是本发明最佳匹配点为其它边界点时区域增长示意图。
图7是本发明最佳匹配点为左、右相连点时区域增长示意图。
图8是本发明节点周围局部积分域的重构过程示意图。
图9是本发明金属目标表面电流展开成的两个正交单位矢量示意图。
图10是本发明实施例中金属目标双站RCS曲线图。
五具体实施方式
下面结合附图对本发明作进一步详细描述。
首先对如下的基本概念作出定义:
内边:有两个邻接三角形的边。
外边:只被一个三角形拥有的边。
活边:每一条新生成的、没有经过处理的边,即此边没有经历找点生成新边的过程。
死边:经历过找点生成新边过程,但没有成功找到新匹配点的边。
内点:若一个节点的所有邻接边都为内边。
外点:没有被选择成为匹配点的节点。
边界边:如果一条边只属于一个三角形,则这条边称为边界边。
边界环:由边界边首尾相连组成的空间多边形。
边界点:边界边上的点。
结合图1~图8,本发明复杂外形金属目标电磁散射的无网格仿真方法,步骤如下:
第1步,节点局部积分域的建立:根据图1复杂外形金属目标表面的离散节点的分布信息构造出每个节点周围各自的局部积分域,具体步骤如下:
(1.1)建立点云数据集合,并使用小立方盒对集合中所有点云数据进行分割,得到小立方盒一维链表。
(1.1.1)首先在空间建立一个笛卡尔直角坐标系xyz,空间中任意一点的坐标值为(x,y,z),设点云数据集合为整个金属目标体被三维立方体区域分割,每个立方体区域称之为一个小立方盒,则有:
nK=[RK/Δu]+1,K=x,y,z
其中Δu表示小立方盒的边长,K表示x轴、y轴或z轴,nK表示K轴方向上的小立方盒的个数,RK表示K轴方向上的距离长度,Ki表示K轴方向上第i个节点;i代表点云数据的节点编号且i=1,2,...,N,N代表点云数据的节点总个数,[]为取整函数,沿K轴方向分别作平面:
其中x0、y0、z0表示的是金属目标体分别在x、y、z三个轴向的最小值, 分别表示金属目标体分别在x、y、z三个轴向的最大值,则六个平面x=x0,y=y0,z=z0,形成的包围盒将点云数据全部包含,而分别平行于三个坐标平面的三组平行平面xj、yj、zj将包围盒分割成了边长相等的小立方盒,此时K轴上的区间:
被分成了nK(K=x,y,z)等份;这些小立方盒共有nx×ny×nz个,所有的小立方盒排成了一个一维链表,这里的一维链表指的是一个一维数组,数组中的每个位置代表一个小立方盒的序号,通过对数组的索引查找出对应的小立方盒的空间位置信息,任意一个小立方盒(ix,iy,iz)在此一维链表中的位置为:
i=iz×nx×ny+iy×nx+ix (4)
在小立方盒的各个坐标分量上分别加上或者减去1,并利用式(4)搜索将找到与该小立方盒相邻的小立方盒。
(1.1.2)将点云数据中的节点分别放入到对应小立方盒中,对任意点pi=(xi,yi,zi),使得:
iK表示第i个节点所在的小立方盒在K轴方向上的编号,这样每一个节点都唯一的对应一个有序数组(ix,iy,iz),即此节点在相应的小立方盒中,从而完成了对整个点云数据的分割。
(1.2)对一维链表中的点云数据进行精简过滤,形成新的一维链表;
(1.2.1)将点云数据文件中的节点读到步骤1.1中的一维链表中;
(1.2.2)随机采样一维链表中5~8个节点,对每个节点计算点云数据中与其最近节点的距离,并求出这些距离的平均值l作为点云数据精简的参考值;
(1.2.3)将PNMOV赋值为一维链表的头节点,在一维链表中删掉PNMOV后面到PNMOV的距离小于DIST=αl的所有节点,所述距离指节点之间的欧式距离,其中α代表节点过滤因子,α的范围为0~1,特别是α=0.5,其值越大,剩余的节点越少,这一步利用步骤(1.1)的点云数据分割,使得PNMOV点的附近点的查找方便快捷,避免了对整个采样一维链表的遍历;
(1.2.4)将PNMOV赋值为一维链表中PNMOV的下一个节点;
(1.2.5)若PNMOV不是一维链表的尾节点,对新的点云数据按照步骤(1.1)中的方法进行分割,然后返回步骤(1.2.3),如果PNMOV是一维链表的尾节点,完成对点云数据的精简过滤,形成新的一维链表。
(1.3)根据步骤(1.2)中新形成的一维链表,构造种子三角形形成剖分区域,建立剖分区域的三角形链表和边界边链表。
(1.3.1)取出步骤(1.2)中新的一维链表中的首个节点,记作为firstpoint;
(1.3.2)取点云数据一维链表中firsitpoint之后的距离firstpoint最近的节点,作为第二个节点,记为secondpoint,这两节点构成一条边;
(1.3.3)在点云数据一维链表中选取firstpoint之后不等于secondpoint,且与firstpoint和secondpoint之间的距离之和最小的节点,作为第三个节点,记为thirdpoint,这三个节点构成一个三角形;
(1.3.4)如果这个三角形的最小内角小于π/6,选择点云一维链表中firstpoint节点的下一个节点作为firstpoint并返回步骤(1.3.2),如果这个三角形的最小内角不小于π/6,将该三角形作为种子三角形,形成剖分区域;
(1.3.5)建立三角形链表,将种子三角形编号设为1并置于三角形链表的首位;建立边界边链表,将剖分区域的边界边置于边界边链表中。
(1.4)寻找边界边的最佳匹配点,将剖分区域进行边界更新得到新的剖分区域,并根据最佳匹配点的情况更新三角形链表和边界边链表,重构出所有节点周围的局部积分域;图2为区域增长的基本概念示意图。
(1.4.1)选边界边链表中一条边界边为当前边,以当前边的两个端点所在的小立方盒为中心,分别从一维链表中找出这两个小立方盒周围邻接的小立方盒;
(1.4.2)对邻接的小立方盒内的节点进行检测,检测各节点是否满足最小内角原则、二面角原则、边长限制原则,将满足以上原则的节点作为候选点;
检测候选点是否满足最小内角原则:即候选点与当前边两端点所构成的三角形的最小内角不小于设定阈值。阈值的确定方式如下:(a)当候选点为当前边的邻接边界点时,它与当前边所构成三角形的最小内角的阈值设为30度;(b)当候选点为其它边界点时,它与当前边所构成三角形的最小内角的阈值设为10度;
检测候选点是否满足二面角原则:分为如下两种情况,(a)当候选点为当前边的邻接边界点时,候选点与当前边所构成的三角形所在面与其有公共边的三角形所在平面形成的二面角要求不小于5π/7;(b)当候选点为其它边界点时,候选点与当前边所构成的三角形所在面与其有公共边的三角形所在平面形成的二面角要求不小于π/2;
检测候选点是否满足边长限制原则:即候选点与当前边所构成三角形的边的边长要小于βDIST,β的范围为1~5,特别是β=3。这避免了自交现象的产生,特别是边界附近不正常边界边的产生;
(1.4.3)如果经过以上检测所得候选点为0个时,将此当前边标注为死边,进入步骤(1.4.5);如果经过以上检测所得到的候选点不止1时,选择候选点中与当前边两端点所构成两边的夹角最大者作为当前边的最佳匹配点;所得候选点只有1个时,将该节点作为当前边的最佳匹配点;所得最佳匹配点与当前边构成新的三角形;
(1.4.4)根据新的三角形,进行边界更新得到新的剖分区域,并根据最佳匹配点的情况更新三角形链表和边界边链表:
情况一:如图3所示,当最佳匹配点为外点时,产生两条新边,新边一和新边二。新边一由当前边的起点与最佳匹配点相连构成,新边二由最佳匹配点与当前边的终点相连构成,在边界边链表中将新边一插入到当前边的后面,新边二插入到新边一的后面,然后删除当前边,最后将新生成的三角形加入到三角形链表中;
情况二:如图4、5所示,当最佳匹配点为当前边的起点或者终点的邻接点时,产生一条新边,此边是由与当前边的起点相连接的边界边的起点和当前边的终点相连而构成,同时产生两条内边和一个新内点:内边一和内边二,内边一是与当前边的起点相连接的边界边,内边二是当前边,新内点是当前边的起点,将新边插入到边界边链表中当前边之后,并将当前边和与当前边的起点相连接的边界边从边界边链表中删除,将新产生的三角形加入到三角形链表中;
情况三:如图6所示,当最佳匹配点为其它边界点时,此时最佳匹配点是与当前边的起、终点不相邻接的边界点,这时产生两条新边,新边一和新边二。新边一由当前边的起点与最佳匹配点相连构成,新边二由最佳匹配点与当前边的终点相连构成,且边界环(在最佳匹配点处分裂成两个,内边界环和外边界环,为了保证剖分过程中边界环始终只有一个,需要检查新边一和新边二是否与边界环中边界边反向重合,若有一条新边与某条边界边重合,则将此新边和边界边标记为死边,并在边界边链表中将新边一插入到当前边的后面,新边二插入到新边一的后面,然后删除当前边,否则在边界边链表中将新边一插入到当前边的后面,新边二插入到新边一的后面,然后删除当前边。但是候选点在几何上既是左相连点同时又是右相连点,如图7所示,这时候应该删掉当前边,以及相关连的另外两条边。
(1.4.5)如果边界边链表中不全为死边时,按照步骤(1.4.1)~(1.4.4)的方法对边界边链表中未被标记为死边的各条边界边不断寻找最佳匹配点,直到边界边链表中没有活边时结束,通过不断更新剖分区域的边界,使得剖分区域不断增长,最终重构出所有节点周围的局部积分域,如图8所示。
第2步,建立金属目标电场积分方程,具体步骤如下:
令均匀平面波照射到复杂外形金属目标上,金属目标的表面产生感应电流J(r′),根据理想导体的电场边界条件,即金属表面的总场切向分量为0,得到金属目标的电场积分方程(EFIE)和磁场积分方程(MFIE),如下:
[Einc(r)+Esca(r)]tan=0 (6)
[Hinc(r)+Hsca(r)]tan=0 (7)
其中,tan表示切向分量,Einc表示入射电场,Hinc表示入射磁场,Esca表示散射电场,Hsca表示散射磁场,具体表达式为:
其中,表示单位并矢,将上述两式代入到式(8)、式(9)中,并且联立方程(6)和(7),得金属目标电场积分方程和磁场积分方程:
其中表示金属目标表面的单位切向分量,P.V.表示主值积分项,表示金属目标表面的单位外法向量,J(r)表示目标表面上场点r处的感应电流。
第3步,离散节点处面电流的展开,构造出每个局部积分域中心位置处的面电流。即将每个局部积分域中心处的面电流用局部积分域周围的节点处的面电流进行近似展开,从而构造出每个局部积分域中心处的面电流。
根据金属目标表面的离散节点的分布信息,将金属目标表面的电流展开,待求的等效表面电流表示如下:
其中,N表示金属目标表面上所有节点的总个数,i表示点云数据中节点编号,Ji(r)为场点r处的未知电流函数在第i个节点处的值,(r)是第i个节点电流函数在场点r处的值,和是金属目标表面场点r处的两个正交的单位切向矢量,和分别表示第i个节点处方向和方向上待求解的相关节点参数,如图9所示。
第4步,对第2步中的电场积分方程采用点匹配进行测试得到阻抗矩阵方程。即将第3步中用作面电流展开时候的系数作为矩阵方程的未知量,并采用点匹配进行测试。
对方程(12)和(13)进行点匹配测试,测试函数选择为场点rm处的和则可得矩阵方程:
其中,rm为节点编号为m的场点, 和是阻抗矩阵,下标mi表示第m个节点与第i个节点之间的作用,上标uu表示方向的测试函数与方向电流函数之间的作用、uv表示方向的测试函数与方向电流函数之间的作用、vu表示方向的测试函数与方向电流函数之间的作用、vv表示方向的测试函数与方向电流函数之间的作用,表示第i个节点处方向上的待求节点参数,表示第i个节点处方向上的待求节点参数,表示第m个节点处方向测试函数的平面波激励,表示第m个节点处方向测试函数的平面波激励。
以MFIE中矩阵元素为例,具体表达式如下:
其中是第i个节点电流函数在场点rm处的值,表示的是场点rm处的方向分量,表示场点rm处的单位法向分量,r′表示源点,表示的是源点r′处的方向分量,是第i个节点电流函数在源点r′处的值,R表示场点rm和源点r′之间的距离。
第5步,联合式(12)和式(13)求解矩阵方程(15),得到电流系数和为了突出地反映金属目标的电磁散射特征,用均匀散射场能量来归一化实际散射场能量,根据电流系数确定雷达散射截面积σ,雷达散射截面RCS常表示成对数形式,单位为分贝(dB):
RCS=10×lg(σ)(18)其中,r表示观察场点的位置,Einc和Esca表示入射电场和散射电场。
实施例1
为了验证本文方法的正确性与有效性,进行了具有金属目标电磁散射的典型仿真,仿真在主频2.83GHz、内存3.5GB的个人计算机上实现,以直径为2m的金属球为例,入射波频率为300MHz,如图9所示,入射波的方向θ=0°,为了验证本发明方法的正确性,以解析方法Mie级数仿真结果作为参照。图10为两种电磁散射特性仿真的RCS曲线图,从图中的曲线可以看出,本文方法不管是在入射波的水平极化还是入射波的垂直极化照射的情况下,在所有的观察角度上都与正确的解析解的数值结果吻合,说明本文方法能够快速仿真分析复杂外形金属目标的电磁散射特性。
综上所述,本发明仅需知道复杂外形金属目标表面离散节点的分布信息,便可对其进行快速的电磁散射仿真,而不依赖于金属目标表面的模型化网格剖分,其实现过程灵活自由,具有很强的实际工程应用价值。
Claims (5)
1.一种复杂外形金属目标电磁散射的无网格仿真方法,其特征在于,步骤如下:
第1步,节点局部积分域的建立:根据复杂外形金属目标表面的离散节点的分布信息构造出每个节点周围各自的局部积分域;所述的节点局部积分域的建立,分为以下步骤:
(1.1)建立点云数据集合,并使用小立方盒对集合中所有点云数据进行分割,得到小立方盒一维链表;
(1.2)对一维链表中的点云数据进行精简过滤,形成新的一维链表;
(1.3)根据新形成的一维链表,构造种子三角形形成剖分区域,建立剖分区域的三角形链表和边界边链表;
(1.4)从一维链表中寻找边界边的最佳匹配点,将剖分区域进行边界更新得到新的剖分区域,并根据最佳匹配点的情况更新三角形链表和边界边链表,重构出所有节点周围的局部积分域;
第2步,建立金属目标电场积分方程;
第3步,展开离散节点处的面电流,构造出每个局部积分域中心位置处的面电流;
第4步,对第2步中的电场积分方程采用点匹配进行测试得到阻抗矩阵方程;
第5步,求解矩阵方程,得到电流系数,根据电流系数确定雷达散射截面积。
2.根据权利要求1所述的复杂外形金属目标电磁散射的无网格仿真方法,其特征在于,步骤(1.1)中所述小立方盒一维链表的具体构建包括以下步骤:
(1.1.1)首先在空间建立一个笛卡尔直角坐标系xyz,空间中任意一点的坐标值为(x,y,z),设点云数据集合为整个金属目标体被三维立方体区域分割,每个立方体区域称之为一个小立方盒,则有:
其中Δu表示小立方盒的边长,K表示x轴、y轴或z轴,nK表示K轴方向上的小立方盒的个数,RK表示K轴方向上的距离长度,Ki表示K轴方向上第i个节点;i代表点云数据的节点编号且i=1,2,...,N,N代表点云数据的节点总个数,[]为取整函数,沿K轴方向分别作平面:
其中x0、y0、z0表示的是金属目标体分别在x、y、z三个轴向的最小值,分别表示金属目标体分别在x、y、z三个轴向的最大值,则六个平面x=x0,y=y0,z=z0,形成的包围盒将点云数据全部包含,而分别平行于三个坐标平面的三组平行平面xj、yj、zj将包围盒分割成了边长相等的小立方盒,此时K轴上的区间:
被分成了nK(K=x,y,z)等份;这些小立方盒共有nx×ny×nz个,所有的小立方盒排成了一个一维链表,这里的一维链表指的是一个一维数组,数组中的每个位置代表一个小立方盒的序号,通过对数组的索引查找出对应的小立方盒的空间位置信息,任意一个小立方盒(ix,iy,iz)在此一维链表中的位置为:
i=iz×nx×ny+iy×nx+ix (4)
在小立方盒的各个坐标分量上分别加上或者减去1,并利用式(4)搜索将找到与该小立方盒相邻的小立方盒;
(1.1.2)将点云数据中的节点分别放入到对应小立方盒中,对任意点pi=(xi,yi,zi),得:
iK表示第i个节点所在的小立方盒在K轴方向上的编号,每一个节点都唯一的对应一个有序数组(ix,iy,iz),即此节点在相应的小立方盒中,从而完成了对整个点云数据的分割。
3.根据权利要求1所述的复杂外形金属目标电磁散射的无网格仿真方法,其特征在于,步骤(1.2)中所述对一维链表中的点云数据进行精简过滤,具体过程如下:
(1.2.1)将点云数据文件中的节点读到步骤(1.1)中的一维链表中;
(1.2.2)随机采样一维链表中5~8个节点,对每个节点计算点云数据中与其最近节点的距离,并求出这些距离的平均值l作为点云数据精简的参考值;
(1.2.3)将PNMOV赋值为一维链表的头节点,在一维链表中删掉PNMOV后面到PNMOV的欧式距离小于DIST=α×l的所有节点,其中α代表节点过滤因子;
(1.2.4)将PNMOV赋值为一维链表中PNMOV的下一个节点;
(1.2.5)若PNMOV不是一维链表的尾节点,对新的点云数据按照步骤(1.1)中的方法进行分割,然后返回步骤(1.2.3),如果PNMOV是一维链表的尾节点,完成对点云数据的精简过滤,更新一维链表。
4.根据权利要求1所述的复杂外形金属目标电磁散射的无网格仿真方法,其特征在于,步骤(1.3)中所述根据新形成的一维链表,构造种子三角形形成剖分区域,建立剖分区域的三角形链表和边界边链表,具体步骤如下:
(1.3.1)取出新形成的一维链表中的首个节点,记作为firstpoint;
(1.3.2)取firsitpoint之后的距离firstpoint最近的节点,作为第二个节点,记为secondpoint,这两节点构成一条边;
(1.3.3)选取firstpoint之后不等于secondpoint,且与firstpoint和secondpoint之间的距离之和最小的节点,作为第三个节点,记为thirdpoint,这三个节点构成一个三角形;
(1.3.4)如果这个三角形的最小内角小于π/6,选择firstpoint节点的下一个节点作为firstpoint并返回(1.3.2),如果这个三角形的最小内角不小于π/6,将该三角形作为种子三角形,形成剖分区域;
(1.3.5)建立三角形链表,将种子三角形置于三角形链表的首位;建立边界边链表,将剖分区域的边界边置于边界边链表中。
5.根据权利要求1所述的复杂外形金属目标电磁散射的无网格仿真方法,其特征在于,步骤(1.4)中所述重构出所有节点周围的局部积分域,具体过程如下:
(1.4.1)选边界边链表中一条边界边为当前边,以当前边的两个端点所在的小立方盒为中心,分别从一维链表中找出这两个小立方盒周围邻接的小立方盒;
(1.4.2)对邻接的小立方盒内的节点进行检测,检测各节点是否满足最小内角原则、二面角原则、边长限制原则,将同时满足以上原则的节点作为候选点;
检测候选点是否满足最小内角原则:即候选点与当前边两端点所构成的三角形的最小内角不小于设定阈值;阈值的确定方式如下:(a)当候选点为当前边的邻接边界点时,它与当前边所构成三角形的最小内角的阈值设为30度;(b)当候选点为其它边界点时,它与当前边所构成三角形的最小内角的阈值设为10度;
检测候选点是否满足二面角原则:分为如下两种情况,(a)当候选点为当前边的邻接边界点时,候选点与当前边所构成的三角形所在面与其有公共边的三角形所在平面形成的二面角要求不小于5π/7;(b)当候选点为其它边界点时,候选点与当前边所构成的三角形所在面与其有公共边的三角形所在平面形成的二面角要求不小于π/2;
检测候选点是否满足边长限制原则:即候选点与当前边所构成三角形的边的边长要小于βDIST,β的范围为1~5;
(1.4.3)如果经过以上检测所得候选点为0个时,将此当前边标注为死边,进入步骤(1.4.5);如果经过以上检测所得到的候选点不止1时,选择候选点中与当前边两端点所构成两边的夹角最大者作为当前边的最佳匹配点;所得候选点只有1个时,将该节点作为当前边的最佳匹配点;所得最佳匹配点与当前边构成新的三角形;
(1.4.4)根据新的三角形,进行边界更新得到新的剖分区域,并更新三角形链表和边界边链表;
(1.4.5)如果边界边链表中不全为死边时,按照步骤(1.4.1)~(1.4.4)的方法对边界边链表中未被标记为死边的各条边界边不断寻找最佳匹配点,直到边界边链表中没有活边时结束,通过不断更新剖分区域的边界,使得剖分区域不断增长,最终重构出所有节点周围的局部积分域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310216541.3A CN103246827B (zh) | 2013-05-31 | 2013-05-31 | 复杂外形金属目标电磁散射的无网格仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310216541.3A CN103246827B (zh) | 2013-05-31 | 2013-05-31 | 复杂外形金属目标电磁散射的无网格仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103246827A CN103246827A (zh) | 2013-08-14 |
CN103246827B true CN103246827B (zh) | 2016-09-21 |
Family
ID=48926343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310216541.3A Active CN103246827B (zh) | 2013-05-31 | 2013-05-31 | 复杂外形金属目标电磁散射的无网格仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103246827B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732050B (zh) * | 2013-12-20 | 2018-04-03 | 南京理工大学 | 雷电脉冲下碳纤维材料飞行目标内的电磁分布预估方法 |
CN104778151B (zh) * | 2014-01-09 | 2017-11-14 | 南京理工大学 | 基于矩量法和抛物线方程的含腔目标电磁散射分析方法 |
CN114299240A (zh) * | 2021-12-20 | 2022-04-08 | 重庆市勘测院 | 一种基于距离阈值的并行点云抽稀方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129523A (zh) * | 2011-03-23 | 2011-07-20 | 南京理工大学 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
-
2013
- 2013-05-31 CN CN201310216541.3A patent/CN103246827B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129523A (zh) * | 2011-03-23 | 2011-07-20 | 南京理工大学 | 基于mda和mlssm的分析复杂目标电磁散射的方法 |
Non-Patent Citations (3)
Title |
---|
An efficient surface integral equation solution to em scattering by chiral objects above a lossy half space;RS Chen等;《Antennas and Propagation》;20090527;第57卷(第11期);3586-3593页 * |
由点云生成网格曲面算法的研究;高福顺;《中国博士学位论文全文数据库基础科学辑》;20100915(第9期);A002-55 * |
电磁散射的无网格方法研究;李灵侠;《中国优秀博硕士学位论文全文数据库 (硕士)工程科技Ⅱ辑 》;20070615(第6期);C042-10 * |
Also Published As
Publication number | Publication date |
---|---|
CN103246827A (zh) | 2013-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106295120A (zh) | 一种目标电磁散射特性快速计算方法 | |
CN104778151B (zh) | 基于矩量法和抛物线方程的含腔目标电磁散射分析方法 | |
EP2887092A1 (en) | Computing radar cross section | |
CN103218487B (zh) | 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法 | |
CN105930567A (zh) | 一种基于子区域自适应积分的电磁散射特性获取方法 | |
Botha et al. | On the variational formulation of hybrid finite element-boundary integral techniques for electromagnetic analysis | |
Liu et al. | Topological properties and fractal analysis of a recurrence network constructed from fractional Brownian motions | |
CN104200074A (zh) | 快速获取目标电磁散射特性的多层复波束方法 | |
CN103246827B (zh) | 复杂外形金属目标电磁散射的无网格仿真方法 | |
Zhang et al. | Broadband monostatic RCS and ISAR computation of large and deep open cavities | |
CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
CN106777536A (zh) | 基于精细电磁仿真的电磁远场二、三维可视化处理方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN104346488B (zh) | 电大复杂外形金属目标混合建模及电磁散射快速仿真方法 | |
CN105277927B (zh) | 飞行器编队瞬态电磁特性时域阶数步进分析方法 | |
CN104731996A (zh) | 一种快速提取电大尺寸金属腔体目标瞬态散射信号的仿真方法 | |
CN106156475A (zh) | 电大尺寸目标的瞬态电磁特性快速提取方法 | |
Walton et al. | Advances in co-volume mesh generation and mesh optimisation techniques | |
CN104915324A (zh) | 腔体含介质目标电磁散射混合分析方法 | |
CN115510690A (zh) | 基于awe技术的电大不确定外形金属目标的新型电磁特性计算方法 | |
CN105205299A (zh) | 电大目标电磁散射特性快速降维分析方法 | |
CN110196961B (zh) | 非合作不确定外形的飞机雷达回波预测方法 | |
CN105184033A (zh) | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 | |
Ikehata | On finding an obstacle with the Leontovich boundary condition via the time domain enclosure method | |
Warecka et al. | Hybrid method analysis of unshielded guiding structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |