[发明专利]一种蜂窝吸波结构电磁快速数值建模方法

权利要求书

1.一种蜂窝吸波结构电磁数值建模方法，其特征在于，包括如下步骤：

步骤1、将整个计算区域分为有限元计算区域与边界元计算区域两部分；其中，有限元计算区域为蜂窝吸波结构的三维模型实体部分，边界元计算区域为蜂窝吸波结构中包围有限元计算区域的闭合曲面；

然后将限元计算区域分成全波精确电磁建模区和均质等效建模区；全波精确电磁建模区为蜂窝吸波结构外侧设定宽度的环形区域；对精确电磁建模区中每个蜂窝单元都要进行实体建模，且每个蜂窝单元的边界为原蜂窝璧的外表面，均质等效建模区为将蜂窝吸波结构除全波精确电磁建模区的中间部分，将其替换成均匀实心介质，最终得到混合等效蜂窝模型；

步骤二：对均质等效建模区使用四面体网格剖分，同时均质等效建模区与全波精确电磁建模区最内侧的实心六棱柱接触的面定义为有限元交界面；全波精确电磁建模区中的实心六棱柱结构同样使用四面体网格剖分；然后对边界元计算区域进行三角形网格剖分；

步骤三、对步骤二剖分后的模型进行计算，获得磁场分布。

2.如权利要求1所述的一种蜂窝吸波结构电磁数值建模方法，其特征在于，所述步骤三的具体方法为：

在全波精确建模区和均质等效建模区的有限元交界面采用完全二阶Robin传输条件；

在有限元计算区域与边界元计算区域的交界面采用一阶Robin传输条件；

对均质等效建模区的电磁场通过有限元泛函变分模拟；

在全波精确建模区域，各蜂窝单元外表面两端的电磁场不连续性条件阻抗边界条件等效；

将全波精确电磁建模区的场通过泛函变分模拟；

将边界元计算区域的边界条件采用联合积分方程描述；

求解得到模型的磁场分布。

3.如权利要求2所述的一种蜂窝吸波结构电磁数值建模方法，其特征在于，所述在全波精确建模区和均质等效建模区的有限元交界面采用完全二阶Robin传输条件具体为：

其中，

表示有限元交界面Γ_mn上最小的网格边尺寸，为虚数单位，为自由空间的波数，其中，μ₀和ε₀分别表示空气中磁导率和介电常数；ω＝2πf表示为角频率，f表示计算的工作频率；和分布表示交界面上的电流和电场，i＝m，n，m和n分别代表全波精确电磁建模区和均质等效建模区。

4.如权利要求3所述的一种蜂窝吸波结构电磁数值建模方法，其特征在于，有限元计算区域外表面与边界元计算区域的交界面采用一阶Robin传输条件，具体为：

和表示有限元计算区域外表面上的电流和电场；和表示边界元计算区域表面上的电流和电场。

5.如权利要求4所述的一种蜂窝吸波结构电磁数值建模方法，其特征在于，均质等效建模区的电磁场通过有限元泛函变分模拟为：

其中，Z₀为自由空间阻抗，均质等效建模区内部使用Ω_m表示，dV表示单位体积，dS对应积分面的单位面积；E_m表示Ω_m中的电场，和分别表示Γ_m和面上的磁场，和分别表示Γ_m和面上的法向；ε_m,r表示Ω_m中的相对介电常数。

6.如权利要求5所述的一种蜂窝吸波结构电磁数值建模方法，其特征在于，所述阻抗边界条件等效为：

其中，是单位法向量；z_s＝jω(ε_r-ε₀)d，ε_r为蜂窝壁的介电常数，d为等效前的蜂窝壁厚度；E⁺表示等效前的蜂窝单元内壁上的电场，E^-表示等效前的蜂窝单元外壁上的电场，E表示等效前的蜂窝壁内的电场；H⁺表示等效前的蜂窝单元内壁上的磁场，H^-表示等效前的蜂窝单元外壁上的磁场。