CN106294894B - 快速分析非均匀目标电磁散射特性的有限元边界积分方法 - Google Patents
快速分析非均匀目标电磁散射特性的有限元边界积分方法 Download PDFInfo
- Publication number
- CN106294894B CN106294894B CN201510251514.9A CN201510251514A CN106294894B CN 106294894 B CN106294894 B CN 106294894B CN 201510251514 A CN201510251514 A CN 201510251514A CN 106294894 B CN106294894 B CN 106294894B
- Authority
- CN
- China
- Prior art keywords
- group
- complex point
- field
- point source
- finite element
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种快速分析非均匀目标电磁散射特性的有限元边界积分方法。该方法结合了多层复点源快速方法与有限元边界积分方法分析非均匀目标电磁散射特性。采用八叉树分组技术对所分析问题进行分组,将组内基函数对远场的作用利用等效原理等效为所在组等效面上数量远小于组内基函数个数的复点源对远场的作用。本发明在计算密网格剖分的非均匀目标的电磁散射特性分析中能够有效地节省计算时间和内存需求。
Description
技术领域
本发明属于目标电磁散射特性数值计算技术领域,特别是结合了多层复点源快速方法和有限元边界积分方法的电磁散射特性快速仿真方法。
背景技术
对于非均匀目标的电磁散射特性分析一直是比较热门的研究领域。许多算法可以分析此类问题。对于微分类方法,有限元方法(FEM)和时域有限差分方法(FDTD)可以用来分析非均匀目标。然而,对于微分类算法,为了得到唯一解,需要使用吸收边界条件(ABC)或完美匹配层(PML)进行边界截断,这样会增大计算规模。对于积分类方法,边界积分方法(BIE)和体面积分方法(VSIE)可以用来分析介质目标,特别是体面积分方法可以以用来分析非均匀问题。这些积分方法自然满足辐射边界条件。然而积分方法的阻抗矩阵是稠密矩阵,相比于微分类方法的稀疏矩阵,这样限制了矩阵方程的求解效率。有限元边界积分方法(X.Q.Sheng,J.M.Jin,J.M.Song,C.C.Lu,and W.C.Chew,“On the formulation of hybridfinite element and boundary integral methods for 3-D scattering,”IEEETrans.Antennas Propagat,vol.46,no.3,pp.303-311,Mar.1998)将边界积分方法引入有限元中替换吸收边界条件或者是完美匹配层,结合了微分类方法和积分类方法的优点。然而,对于电大尺寸的非均匀目标问题,为了保证计算精度需要密网格剖分,这样会导致边界面上较大的未知量,这样会大大限制求解效率。而多层复点源技术能够有效的加速大未知量密网格问题的求解,从而能够结合有限元边界积分方法有效的分析电大尺寸非均匀目标电磁散射特性问题。
发明内容
本发明的目的在于提供一种快速分析非均匀目标电磁散射特性的有限元边界积分方法。
实现本发明目的的技术解决方案为:一种快速分析非均匀目标电磁散射特性的有限元边界积分方法,步骤如下:
步骤1、建立非均匀目标的物理模型,使用剖分软件对介质区域用四面体剖分,对物体表面采用三角形剖分,得到物理模型的结构信息,即四面体的编号及各节点坐标和三角形的编号及各节点坐标;
步骤2、将所分析问题利用八叉树分组技术进行分组,得到一种多层结构,在每一层中,以组的中心为球心,建立等效球面将组全部包围,按照各层的复点源个数将复点源均匀分布在等效球面上;
步骤3、将组内基函数对远场的作用利用等效原理用所在组等效球面上复点源进行展开,利用复点源展开系数表示远场组的相互作用矩阵;
步骤4、将复点源表示的远场矩阵矢量乘引入有限元边界积分方程的迭代求解器中对矩阵方程进行迭代求解,得到有限元计算区域的电场系数和边界积分区域的电流系数,即可得到分析物体的远场散射场信息,从而可以计算的到目标的雷达散射截面。
本发明与现有技术相比,其显著优点为:(1)将组内基函数的远场相互作用转化为等效面上数量小于组内基函数个数的复点源的远场相互作用,可以有效的加快矩阵矢量乘及降低远场内存需求,从而降低算法复杂度。(2)与格林函数无关,可以灵活有效地分析非均匀目标,程序实现起来比较容易。
附图(表)说明
图1是等效面和测试面位置示意图。
图2是本发明某散射体示意图。
图3是本发明实施例中某散射体双站RCS曲线图。
表1是本发明与其它方法时间及内存上的比较。
具体实施方式
下面结合附图对本发明作进一步详细描述。
步骤1、建立非均匀目标的物理模型,使用剖分软件对介质区域用四面体剖分,对物体表面采用三角形剖分,得到物理模型的结构信息,即四面体的编号及各节点坐标和三角形的编号及各节点坐标;
步骤2、将所分析问题利用八叉树分组技术进行分组,得到一种多层结构,在每一层中,以组的中心为球心,建立等效球面将组全部包围,按照各层的复点源个数将复点源均匀分布在等效球面上;
步骤3、将组内基函数对远场的作用利用等效原理用所在组等效球面上复点源进行展开,利用复点源展开系数表示远场组的相互作用矩阵;具体步骤如下:
有限元边界积分中的边界积分方程表达式如下所示:
J和M分别表示感应电流和感应磁流,Einc和Hinc表示入射电场和入射磁场。η表示自由空间波阻抗,算子L和K的表达式为
电流J和磁流M用RWG基函数展开,用迦辽金方法对边界积分方程进行测试。公式(5)中远场相互作用的运算符可以表示为
在公式(10)和公式(11)中,积分Sn可以用复点源表示为
这里α表示运算符L或K。Q是复点源的个数,r′nq为第q个复点源的位置,wnq为对于fn(r′)第q个复点源的展开系数。所以公式(10)和公式(11)可以表示为
这里C表示与积分有关的常数,然后公式(15)可以变为
积分Sm可以用复点源表示为
公式(12)和公式(13)可以变为
这里D表示与积分无关的常数。将公式(14)带入公式(19)并且利用并失格林函数的对称性,公式(19)可以表示为
积分Sm可以用复点源表示为
因此公式(12)和(13)表示的阻抗矩阵可以表示为
将公式(18)和(22)按照公式(5)的形式组合即可得到边界积分方程远场相互作用矩阵复点源的表示形式。
步骤4、将复点源表示的远场矩阵矢量乘引入有限元边界积分方程的迭代求解器中对矩阵方程进行迭代求解,得到有限元计算区域的电场系数和边界积分区域的电流系数,即可得到分析物体的远场散射场信息,从而可以计算的到目标的雷达散射截面。具体表达式如下:
ES(r)为远场r点处散射场,j是虚数符号,k表示的波数,η表示自由空间波阻抗,r′表示电流系数所在基函数的位置,表示r′的单位方向向量,Jθ和Jφ分别表示电流系数的及分量。目标雷达散射截面的表达式为:
σ表示的是雷达散射截面,Einc表示的是入射场。
实施例
本实施例进行了电磁散射的典型仿真,仿真在主频2.0GHz、内存512GB的服务器上实现,以x,y,z三个方向尺寸分别为3.49m,0.698m,0.698m的尖顶拱模型为例,如图2所示,其表面有一层厚度为0.03m,相对介电常数为(4.0,0.5)的涂覆。以0.03λ剖分尺寸对模型进行剖分,得到41570个四面体和13138个三角形单元。在最细层采用0.2λ波长分组得到三层结构。入射波频率为300MHz,入射波的方向θ=0°,观察角度为分别采用本发明方法和传统的有限元边界积分方法对目标进行电磁散射分析。图3为两种电磁散射特性仿真的RCS曲线图,从图中的曲线可以看出,两种方法得到的RCS曲线吻合的很好。在表一中给出了两种方法在分析过程中内存和时间的使用情况,从表中可以看出本发明提出的方法相比于传统的有限元边界积分方法在内存和时间上都有很大程度节省。
表1。
Claims (1)
1.一种快速分析非均匀目标电磁散射特性的有限元边界积分方法,其特征在于步骤如下:
步骤1,建立非均匀目标的物理模型,对介质区域用四面体剖分,对物体表面采用三角形剖分,得到物理模型的结构信息,即四面体的编号及各节点坐标、三角形的编号及各节点坐标;
步骤2,将所分析问题利用八叉树分组方法进行分组,得到一种多层结构,在每一层中,以组的中心为球心,建立等效球面将组全部包围,按照各层的复点源个数将复点源均匀分布在等效球面上;
步骤3,将组内基函数对远场的作用利用等效原理用所在组等效球面上复点源进行展开,利用复点源展开系数表示远场组的相互作用矩阵;
步骤4,将复点源表示的远场矩阵矢量乘引入有限元边界积分方程的迭代求解器中,对矩阵方程进行迭代求解,得到有限元计算区域的电场系数和边界积分区域的电流系数,即可得到分析目标的远场散射场信息,从而可以计算的到目标的雷达散射截面RCS;
步骤2的具体实现如下:
用一个立方体将目标体包围住,该立方体定义为第零层的第一个且是最后一个组,把该立方体等分为八个子立方体形成第一层组,这八个子立方体称为第零层立方体组的子层组,第零层立方体组则称为八个子层组的父层组,然后再对每个子立方体组进行与上一步相同的细分得到第二层组,重复上述对每个子立方体组进行的操作得到一个多层分组的结构,并以第一步的网格信息来判断最底层立方体的尺寸;
步骤3所述将组内基函数对远场的作用利用等效原理用所在组等效球面上复点源进行展开具体为:
步骤3所述利用复点源展开系数表示远场组的相互作用矩阵,具体为:
C为与积分无关的常数,P为m基函数所在组复点源个数,ump为对于m基函数第p个复点源展开的系数;
步骤4所述得到有限元计算区域的电场系数和边界积分区域的电流系数,即可得到分析目标的远场散射场信息,从而计算得到目标的雷达散射截面RCS,具体表示形式如下所示:
ES(r)为远场r点处散射场,j是虚数符号,k表示的波数,η表示自由空间波阻抗,r′表示电流系数所在基函数的位置,表示r′的单位方向向量,Jθ和Jφ分别表示电流系数的及分量;目标雷达散射截面的表达式为:
σ表示的是雷达散射截面,Einc表示的是入射场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510251514.9A CN106294894B (zh) | 2015-05-15 | 2015-05-15 | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510251514.9A CN106294894B (zh) | 2015-05-15 | 2015-05-15 | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106294894A CN106294894A (zh) | 2017-01-04 |
CN106294894B true CN106294894B (zh) | 2020-06-19 |
Family
ID=57632278
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510251514.9A Active CN106294894B (zh) | 2015-05-15 | 2015-05-15 | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106294894B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107609274B (zh) * | 2017-09-14 | 2020-03-13 | 哈尔滨工业大学 | 基于传输线与级别调度法的二维静磁场并行有限元方法 |
CN108230414A (zh) * | 2017-12-13 | 2018-06-29 | 南京理工大学 | 一种低高频混合多尺度目标低秩压缩分析方法 |
CN108519524A (zh) * | 2018-03-19 | 2018-09-11 | 江苏伏波海洋探测科技有限公司 | 舰船静电场的基于点电源的深度换算方法 |
CN108575055B (zh) * | 2018-04-24 | 2020-11-13 | 广东工业大学 | 一种差分过孔分析方法 |
CN112016232B (zh) * | 2020-08-31 | 2024-06-14 | 中国原子能科学研究院 | 一种撕裂有限元过程处理方法及系统 |
CN112989677B (zh) * | 2021-04-20 | 2021-08-20 | 北京智芯仿真科技有限公司 | 大规模集成电路层间耦合部分累加计算的迭代方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331995A (zh) * | 2011-07-08 | 2012-01-25 | 华东师范大学 | 一种基于点源的地球信息三维模型的信息采集方法和系统 |
CN104200074A (zh) * | 2014-08-18 | 2014-12-10 | 南京理工大学 | 快速获取目标电磁散射特性的多层复波束方法 |
-
2015
- 2015-05-15 CN CN201510251514.9A patent/CN106294894B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102331995A (zh) * | 2011-07-08 | 2012-01-25 | 华东师范大学 | 一种基于点源的地球信息三维模型的信息采集方法和系统 |
CN104200074A (zh) * | 2014-08-18 | 2014-12-10 | 南京理工大学 | 快速获取目标电磁散射特性的多层复波束方法 |
Non-Patent Citations (2)
Title |
---|
Nested Equivalent Source Approximation for the Modeling of Multiscale Structures;Mengmeng Li等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20140731;第3664-3678页 * |
自适应快速多级子边界面算法分析大规模声学问题;王现辉;《中国博士学位论文全文数据库-信息科技辑》;20140115;第I138-70页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106294894A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106294894B (zh) | 快速分析非均匀目标电磁散射特性的有限元边界积分方法 | |
Sumithra et al. | Review on computational electromagnetics | |
CN110737873B (zh) | 一种大规模阵列天线散射的快速分析方法 | |
CN111046603A (zh) | 基于gpu并行加速特征基函数算法的电磁散射特性分析方法 | |
CN109765538B (zh) | 非均匀介质目标体的散射场确定方法 | |
Zhang et al. | Performance of a massively parallel higher-order method of moments code using thousands of CPUs and its applications | |
CN102081690A (zh) | 复杂电路的矩阵分解结合新奇异值分解方法 | |
CN102708229A (zh) | 复杂分层媒质结构的矩阵分解结合新奇异值分解方法 | |
CN112329303A (zh) | 基于有限元区域分解的阵列天线电磁特性求解方法 | |
CN108052738A (zh) | 色散媒质的高阶局部无条件稳定时域间断伽辽金分析方法 | |
Shi et al. | Comparison of interpolating functions and interpolating points in full-wave multilevel Green's function interpolation method | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN107305536B (zh) | 混合阶时域不连续伽略金方法 | |
Xu et al. | Improved synthetic basis functions method for nonperiodic scaling structures with arbitrary spatial attitudes | |
CN111832183B (zh) | 基于八叉树结构mlfma区域分解的电磁仿真方法 | |
CN115408972A (zh) | 多层印刷电路板电路参数求解方法 | |
CN109063386B (zh) | 一种密集旋转对称体群电磁散射仿真方法 | |
CN104714929B (zh) | 一种实现ah‑fdtd算法按阶并行求解的方法 | |
CN114491992A (zh) | 一种基于等效偶极矩和物理光学法的高效电磁散射方法 | |
CN114239239A (zh) | 用于弹目交会目标电磁特性快速模拟的直接稀疏求解方法 | |
CN106294898B (zh) | 一种加速分析介质目标电磁散射特性的复点源求解方法 | |
CN111931353A (zh) | 一种应用于仿真fss结构的散射场求解方法 | |
CN106156431A (zh) | 基于nurbs曲面建模的导体目标电磁散射仿真方法 | |
Ong et al. | Full-Wave analysis of large-scale interconnects using the multilevel UV method with the sparse matrix iterative approach (SMIA) | |
Chen et al. | Analysis of electromagnetic scattering problems by means of a VSIE-ODDM-MLFMA method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |