一种计算蜂窝复合材料目标电磁散射的方法

技术领域

本发明涉及雷达技术领域，特别涉及一种计算蜂窝复合材料目标电磁散射的方法。

背景技术

如今，蜂窝复合材料因其高强度-重量比，耐高温性和高微波吸收性而广泛用于航空航天和飞行器设计领域。因此，蜂窝结构覆盖的目标电磁散射成为重要的课题。克服此问题有三个挑战。第一，蜂窝复合材料的成分复杂度高，所以很难描述其电磁参数。第二，蜂窝的精细结构需要大量的网格来描述其几何轮廓，这些网格将导致电磁仿真中巨大的计算量。第三，每个六边形蜂窝是电小尺寸，而目标通常是电大尺寸，解决具有不同网格尺寸的多尺度问题是一个巨大的挑战。所有这些问题使人们非常关注如何从具有复合材料的目标中有效地建立散射模型。

在电磁散射领域中，所有的模拟方法都可以分为两类：数值方法(例如矩量法(MoM)，时域有限差分法(FDTD)，有限元方法(FEM)等)和高频方法(例如物理光学(PO)，弹跳射线法(SBR)，迭代物理光学(IPO)等)。目标和复杂的蜂窝都将在网格划分中产生大量未知，处理这些未知数远远超出了数值方法的能力。因此，高频方法更为实用，在高频法中，SBR被广泛用于解决复杂目标的散射。然而，在这个问题中，蜂窝单元电尺寸很小，每个单元的长度仅为约2-4毫米，高频方法在解决电小问题时不够精确。而且，蜂窝内壁中通常涂有损耗材料。每当射线撞击界面时，它将分成两束射线(反射射线和透射射线)。这将大大增加计算成本。

为了克服上述问题，我们提出了一种从蜂窝复合材料覆盖的目标上获得散射特性的新方法。根据等效原理，周期性蜂窝结构可以用强干涉理论等效为具有一定介电常数和磁导率的介质层。然后通过广义传播矩阵法推导了各向异性材料在每个入射角的反射系数。然后将得到的反射系数代入SBR方法中，得到覆盖蜂窝复合材料的目标雷达散射截面。这样，可以简化复杂蜂窝状复合材料的散射计算，在良好的精度下大大降低计算成本。

发明内容

本发明提供一种计算蜂窝复合材料目标电磁散射的方法，能够简化复杂蜂窝状复合材料的散射场的计算，在良好的精度下大大降低计算成本。

本发明提供了一种计算蜂窝复合材料目标电磁散射的方法，包括以下步骤：

S1、根据强干扰理论确定蜂窝复合材料的等效介电常数和等效磁导率；

S2、根据等效介电常数和等效磁导率，采用广义传播矩阵法确定带有PEC(PerfectlyElectricalConductor)基板的蜂窝复合材料在不同入射角下的反射系数；

S3、根据反射系数计算入射电磁波在蜂窝复合材料目标表面产生的感应电流，采用物理光学积分法确定每次反射时产生感应电流的散射场，叠加得到总散射场。

上述步骤S1计算等效介电常数和等效磁导率的方法为：

将蜂窝看作具有等效介电常数ε＝diag(ε

ε

其中，ε

ε

其中，ε

相应地，垂直于蜂窝轴线方向的等效磁导率μ

μ

其中，μ

μ

其中，μ

上述步骤S2计算蜂窝复合材料反射系数的方法为：

蜂窝和PEC基板的反射系数R

R

其中，exp(-j2kh)为空间相位延迟，h为蜂窝单元格的高度，θ为入射角，

上述步骤S3计算入射电磁波在蜂窝复合材料目标表面产生感应电流的散射场的方法为：

将入射电磁波认为是一系列平行传播的射线管，射线从射线管中发出，多次反射追踪N阶反射射线的方向，射线管初始的传播方向，即入射波传播方向为

第N次反射的反射方向为：

其中

第N次反射的反射电场强度为：

其中

第N次反射磁场强度为：

其中

射线在目标表面产生感应电磁流分别为：

其中R是蜂窝复合材料的反射系数，进而通过对所有感应电流的辐射求和来获得散射场；

其中

根据物理光学积分，得到每次反射感应电流散射场为：

其中，r为射线与目标相交点到接收雷达的距离，

目标总散射场为

目标雷达散射截面(RCS)为：

其中，

与现有技术相比，本发明的有益效果在于：

本发明利用强干扰理论和广义传播矩阵法得到蜂窝复合材料的等效介电常数，等效磁导率和反射系数，进而利用物理光学积分法计算得每次反射产生的散射场再叠加得到总散射场。实验表明，本发明的方法具有较高的仿真精度，而且对计算机的资源占用也比较小，能够简化复杂蜂窝状复合材料的散射场的计算，在良好的精度下大大降低计算成本。

附图说明

图1为蜂窝单元的几何形状。

图1(a)是单位单元的俯视图；图1(b)展示出了单位单元的侧视图。

图2为蜂窝模型及其等效模型。

图2(a)显示了PEC基板上蜂窝的实际模型；图2(b)显示了等效模型。

图3为发射面与入射方向示意图。

图4为射线追踪示意图。

图5为MoM和等效模型在不同入射角下获得的反射系数的比较。

图6为MoM和等效模型在不同频率下获得的反射系数的比较。

图7为实际的蜂窝模型及其等效模型。

图8为实际模型和等效模型的双站RCS的比较。

图9为实测模型和环境。

图9(a)显示了实测模型和环境；图9(b)是模型的照片；图9(c)示出了蜂窝布置的细节。

图10为等效模型与实测数据的RCS比较。

图11为带蜂窝和不带蜂窝的飞机的单站散射。

具体实施方式

下面结合附图1-11，对本发明的一个具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。

如图1所示，本发明实施例提供的一种计算蜂窝复合材料目标电磁散射的方法，包括以下步骤：

步骤1，根据强干扰理论确定蜂窝复合材料的等效介电常数和等效磁导率；

蜂窝的几何形状如图1所示。根据强干扰理论，当入射波的波长远大于单元格的大小时，可以忽略蜂窝的每个单元格之间的耦合。图1(a)是蜂窝单元的俯视图。图1(b)展示出了蜂窝单元的侧视图。材料1是蜂窝壁，蜂窝壁的厚度为t，两个相邻的蜂窝单元格之间的距离为2a。材料2是涂覆在每个蜂窝内壁表面上的有损耗材料，d是材料2的厚度。材料3是其他填充材料(通常是空气)，每个单元格的高度是h。

蜂窝可以看作是具有等效介电常数ε＝diag(ε

ε

其中，ε

ε

其中，ε

将以上公式中ε全部替换为μ即可得到平行于蜂窝轴线方向和垂直于蜂窝轴线方向的等效磁导率μ

步骤2，通过建立蜂窝模型的等效模型以及步骤1得到的等效介电常数和等效磁导率，采用广义传播矩阵法确定带有PEC基底的蜂窝复合材料在不同入射角下的反射系数；

基于广义传播矩阵法，我们推导了带有PEC基底的蜂窝复合材料的反射系数。图2(a)显示了PEC基板上蜂窝的实际模型。图2(b)显示了等效模型。区域0是充满空气的半空间。区域1是蜂窝复合材料，可以等效于具有特定介电常数和磁导率的材料。

根据麦克斯韦方程，我们可以获得以下矩阵方程：

其中S＝[E

S(z)＝a·exp(jλz)·B

其中，a＝[a

入射角是θ。在区域0中，特征值λ＝±k

a

a

a

其中，

在区域0中，

在区域1中，

其中Γ是具有PEC基底的蜂窝的广义反射系数，其定义为：

下标o和e分别表示寻常波和非寻常波。下标+和-表示波沿z轴或与z轴相反传播。Γ

根据区域0和区域1的界面上的边界条件，我们有：

然后可以得到：

其中，Γ

代入上式，可以获得广义反射系数和透射系数。

根据PEC表面的边界条件，电场的切向分量为零，E

代入上式，可获得蜂窝和PEC基板的广义反射系数。

在各向同性介质中，寻常波和非寻常波会退化为水平(h)极化和垂直(v)极化。广义反射系数与反射系数之间的关系如下：

总之，R

R

其中，

步骤3，将入射电磁波认为是一系列平行传播的射线管，射线从射线管中发出，多次反射追踪N阶反射射线的方向，通过步骤2得到的反射系数计算射线在目标表面产生的感应电流，进而采用物理光学积分法确定每次反射产生的感应电流散射场，叠加得到总散射场。

将发射至目标的入射电磁波认为是一系列平行传播的射线管，射线从射线管中发出，确定射线在初始传播方向上的磁场矢量。如图3所示，射线管从发射面发射出来，发射面是与入射方向垂直的虚拟平面，该平面与目标中心的距离需要足够远，使得任意延扩该平面都不会与目标相交。每一根射线管都从发射面发出，为保证精度，射线管的边长设置为入射波长的十分之一。

射线管初始的传播方向，即入射波传播方向为

如图4所示，第N次反射的反射方向为：

其中

第N次反射的反射电场强度为：

其中

第N次反射磁场强度为：

其中

射线在目标表面产生感应电磁流分别为：

其中R是蜂窝复合材料的反射系数。进而可以通过对所有感应电流的辐射求和来获得散射场。

其中

根据物理光学积分，得到每次反射感应电流散射场为：

其中，r为射线与目标相交点到接收雷达的距离，

目标总场为

目标雷达散射截面(RCS)为：

其中，一般情况下

仿真结果：

首先，我们用Feko中的MoM验证了我们理论得到的反射系数。在Feko中，仅建立一个单位单元的物理模型，然后使用周期性边界条件来模拟无限周期蜂窝结构的反射。在此模拟中，蜂窝的高度h＝5mm，蜂窝框架的厚度t＝0.1mm，两个相邻蜂窝单元之间的距离2a＝4.763mm，涂层材料的厚度d＝0.04mm，蜂窝框架的介电常数和磁导率为ε

然后，我们用MoM计算带有PEC基板的实际蜂窝中的散射，以及用SBR-PO方法计算等效平板(如图7所示)的散射，入射频率为10GHz，平板的大小为0.4m×0.4m，蜂窝的高度为h＝5mm，入射角θ

接着，通过与实测数据进行比较来证实所提出的模型。我们测量了蜂窝和PEC基底的单站散射。图9(a)显示了实测模型和环境，图9(b)是模型的照片，图9(c)示出了蜂窝布置的细节。被测模型的大小为0.6m×0.6m，蜂窝的高度h＝5mm，蜂窝框架的厚度t＝0.04mm，两个相邻蜂窝单元之间的距离2a＝5.1962mm。图10显示了我们的模型和实测数据的RCS比较，入射角从0°改变为90°。两条曲线显示出良好的一致性，证明了所提出模型的准确性。

然后，将所提出的模型用于通过SBR-PO方法计算覆盖有蜂窝的飞机的散射。将模拟结果与相同尺寸的PEC机翼进行比较。图11显示了带有和不带有蜂窝涂层的飞机的RCS与不同入射角的关系。蜂窝参数与图5相同。仿真结果表明，蜂窝涂层可以大大降低目标的RCS。

本发明公开了一种计算蜂窝复合材料目标电磁散射的有效方法，涉及雷达技术领域，通过将等效反射系数模型与高频方法相结合，提出了一种新的蜂窝复合材料目标电磁散射模型。基于强干涉理论，将蜂窝复合材料视为单轴介质。然后，采用广义传播矩阵法推导了带有PEC基底的蜂窝复合材料的反射系数。最后，将所得到的反射系数引入弹跳射线法和物理光学(SBR-PO)方法中，以获得蜂窝材料覆盖的复杂目标的散射特性。通过分别与数值方法和实测数据进行比较，证明了我们方法的准确性。本发明模拟并讨论了蜂窝材料覆盖目标的散射。

以上公开的仅为本发明的几个具体实施例，但是，本发明实施例并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。