CN106021764A - 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 - Google Patents
面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 Download PDFInfo
- Publication number
- CN106021764A CN106021764A CN201610369219.8A CN201610369219A CN106021764A CN 106021764 A CN106021764 A CN 106021764A CN 201610369219 A CN201610369219 A CN 201610369219A CN 106021764 A CN106021764 A CN 106021764A
- Authority
- CN
- China
- Prior art keywords
- array antenna
- active phase
- antenna
- confidence level
- phase array
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,包括确定有源相控阵天线性能的影响因素;建立天线有限元模型;进行温度场分析,计算热参数的置信度;基于阵面温度均方根误差计算网格大小的置信度;计算天线结构热变形,提取天线单元几何中心节点的位置偏移量;计算结构位移提取的置信度;使用机电耦合模型计算天线的电性能;将计算结果与电磁仿真软件计算结果进行对比,计算机电耦合模型的置信度;基于层次分析法,确定加权系数;建立置信度计算公式,计算并最终确定有源相控阵天线性能仿真的置信度。本发明建立了天线性能仿真置信度计算方法,给出了有源相控阵天线结构热变形对电性能影响仿真计算结果的评判标准。
Description
技术领域
本发明属于天线技术领域,具体涉及面向机电耦合的有源相控阵天线性能仿真置信度的计算方法。
背景技术
随着电子科学技术的飞速发展,有源相控阵天线因其具有快速改变天线波束指向和波束形状、可用于多部发射机在空间进行功率合成、易于形成多个发射与接收波束、可使有源相控阵天线与雷达平台共形等特点,在观测高速运动目标、实现多种雷达功能和多目标跟踪等方面具有无可比拟的优势,其在雷达系统中得到广泛的应用。目前,有源相控阵天线已应用到预警、引导、制导、火控和超视距探测等几乎所有雷达应用领域。可以说有源相控阵天线已经成为当今实用雷达天线发展的主流。
有源相控阵天线阵面是雷达的核心结构部分,天线电性能在很大程度上依赖于其机械结构,天线阵面作为电磁信号传输的载体和边界条件,其位移场直接影响着电磁场在空间中的幅度和相位分布。有源相控阵天线在服役阶段,受到高功率器件T/R组件发热导致有源相控阵天线阵面发生结构热变形,从而导致天线的电性能恶化。
为了预测实际工况下有源相控阵天线结构热变形对电性能的影响,有必要采用仿真计算软件模拟结构热变形对电性能的影响。可首先采用有限元分析软件ANSYS计算天线的温度场和结构位移场的分布,得到天线结构的热变形,再利用有源相控阵天线机电耦合公式计算其电性能的变化,快速预测电性能。然而,在整个仿真计算过程中,模型热参数的选取、有限元网格的尺寸、结构位移提取、机电耦合模型计算电性能的精确度等均会直接影响计算结果,为了准确的衡量仿真计算结果,有必要建立有源相控阵天线性能仿真置信度的计算方法,给出有源相控阵天线性能仿真计算结果的评判标准。
发明内容
针对上述问题,本发明建立了有源相控阵天线性能仿真置信度的计算方法,确定有源相控阵天线性能仿真置信度的影响因素,给出了仿真计算结果的评判标准,从而确保仿真结果的可靠性。
实现本发明目的的技术解决方案是,面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,该方法包括下述步骤:
(1)根据有源相控阵天线的工作要求,确定有源相控阵天线性能的影响因素,以及天线结构参数、材料属性和电磁参数;
(2)根据有源相控阵天线的结构参数和材料属性,建立有源相控阵天线有限元模型;
(3)根据有源相控阵天线的工作要求和有限元模型,施加边界条件,计算天线阵面温度场分布;
(4)根据有源相控阵天线的材料热属性,基于阵面温度场分布,计算有源相控阵天线热参数的置信度P1;
(5)根据天线阵面温度场分布,基于阵面温度均方根误差RMS,计算网格划分大小的置信度P2;
(6)将阵面温度场分布作为载荷,施加到天线有限元模型,计算有源相控阵天线阵面的结构热变形;
(7)提取天线单元几何中心节点的位置偏移量,确定天线单元位置偏移量的均值,计算结构位移提取的置信度P3;
(8)根据有源相控阵天线的电磁工作参数以及单元几何中心的位置偏移量,使用有源相控阵天线机电耦合模型,计算天线的电性能;
(9)将机电耦合模型的计算结果与电磁仿真软件的计算结果进行对比,计算机电耦合模型的置信度P4;
(10)基于层次分析法,确定热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1、α2、α3、α4;
(11)建立置信度计算公式P=α1·P1+α2·P2+α3·P3+α4·P4,计算并最终确定有源相控阵天线性能仿真的置信度。
所述步骤(1)中,有源相控阵天线性能影响因素包括热参数、网格大小、结构位移提取和机电耦合模型;所述有源相控阵天线的结构参数,包括天线单元的行数M、列数N及单元间距dx和dy、印制板、铝蜂窝上面板、铝蜂窝、铝蜂窝下面板及T/R组件的几何参数;所述有源相控阵天线的材料属性包括弹性模量、泊松比、剪切模量、密度、导热系数和热膨胀系数;所述有源相控阵天线的电磁参数包括天线的单元形式和中心工作频率f。
所述步骤(2)中,确定有源相控阵天线单元的几何中心,并在几何中心处设置硬点;在ANSYS中建立有源相控阵天线有限元模型;
所述步骤(4)中,根据有源相控阵天线的材料热属性,基于阵面温度场分布,计算有源相控阵天线在建模过程中热参数的置信度,公式如下:
式中,P1为热参数的置信度,Tideal为建模初始温度时其热膨胀系数,Treal为最终稳定状态下的热膨胀系数。
所述步骤(5)计算网格划分的置信度,按照以下步骤进行:
(5a)选取一组网格大小Mi(1≤i≤n),确定边界条件,并施加天线环境热载荷,得到天线阵面的温度场分布,计算天线阵面温度均方根误差RMSi,当RMSi+1与RMSi的绝对值误差时,则将第i+1组网格大小Mi+1和均方根误差RMSi+1作为标准,分别记为Mc和RMSc。
(5b)当选取网格大小为Md时,计算天线阵面温度均方根误差RMSd,得到网格划分大小的置信度,公式如下:
所述步骤(6)中,将T/R组件热功耗作用下的天线温度场分布作为载荷,施加于天线有限元模型,确定结构约束,计算有源相控阵天线阵面的结构热变形。
所述步骤(7)中,根据有源相控阵天线的结构热变形,提取第(m,n)个天线单元几何中心的位置偏移量δmn,用天线单元几何中心的位移代替天线单元整体位移,计算结构位移提取的置信度P3,公式如下:
式中,为天线单元位置偏移量的均值。
所述步骤(8)利用机电耦合模型,计算天线电性能按照以下步骤进行:
(8a)设有源相控阵天线中共有M×N个天线单元,第(m,n)个天线单元的位置矢量为观察点P相对于坐标系O-xyz所在的方向(θ,φ)上的方位矢量为其中:
式中,xmn、ymn、zmn分别为位置矢量在坐标轴x、y、z方向的投影; 分别为坐标轴x、y、z方向的单位矢量;αx、αy、αz分别为位置矢量与坐标轴x、y、z方向的夹角;
(8b)当天线阵面产生变形时,第(m,n)个天线单元产生的位置偏移量为则变形后第(m,n)个天线单元的位置矢量为:
式中,Δxmn、Δymn、Δzmn分别为位置偏移量在坐标轴x、y、z方向的投影;
(8c)根据有源相控阵天线电磁叠加原理,理想情况下有源相控阵天线的方向图函数为:
结合(8b)得到变形有源相控阵天线方向图函数为:
式中,为有源相控阵天线的波常数,λ为有源相控阵天线工作波长,Ee为天线单元方向图函数,Amn是天线单元激励电流幅度,j为虚数单位,是天线单元激励电流相位;
(8d)根据有源相控阵天线方向图函数绘制的天线方向图,计算得到变形后有源相控阵天线的增益Gc。
所述步骤(9)计算机电耦合模型的置信度按照以下步骤进行:
(9a)利用电磁分析软件HFSS计算有源相控阵天线的增益Gd;
(9b)将步骤(8d)中利用机电耦合模型计算得到的增益Gc和利用HFSS软件仿真得到的增益Gd进行对比,计算机电耦合模型的置信度P4,公式如下:
所述步骤(10)中,确定热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数,通过下述方法进行:
(10a)建立有源相控阵天线性能仿真置信度的层次结构模型,以有源相控阵天线性能仿真的置信度P为标准,分别比较热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4相对于有源相控阵天线性能仿真置信度P的重要性,从而构造出判断矩阵R;
(10b)设热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1,α2,α3,α4,记为W=[α1 α2 α3 α4];
(10c)由判断矩阵R计算判断矩阵R最大的特征值λ,再由公式RW=λW,求得W=[α1 α2 α3 α4]。
本发明与现有技术相比,具有以下特点:
1.本发明可同时研究多因素对有源相控阵天线性能仿真置信度的影响,克服了传统方法中只考虑单一因素下有源相控阵天线性能仿真置信度的不足。
2.本发明综合分析了热参数、网格大小、结构位移提取和机电耦合模型对有源相控阵天线性能仿真置信度的影响,基于层次分析法,确定加权系数,建立有源相控阵天线性能仿真置信度的计算公式,给出了有源相控阵天线结构热变形对电性能影响仿真计算结果的评判标准,具有明显的工程应用价值。本发明中的计算方法对结构热变形对电性能影响这一类仿真置信度的分析具有普遍适用性。
附图说明
图1是本发明面向机电耦合的有源相控阵天线性能仿真置信度的计算方法的流程图;
图2是有源相控阵天线的结构示意图;
图3是有源相控阵天线的有限元模型;
图4是有源相控阵天线温度场分布云图;
图5是有源相控阵天线的热变形分布云图;
图6是φ=0°时,有源相控阵天线变形后HFSS仿真得到的和机电耦合模型计算得到的增益方向图;
图7是φ=90°时,有源相控阵天线变形后HFSS仿真得到的和机电耦合模型计算得到的增益方向图。
图8是有源相控阵天线性能仿真置信度的层次结构模型。
具体实施方式
下面结合附图及实施例对本发明做进一步说明
参照图1,本发明为面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,具体步骤如下:
步骤1,确定有源相控阵天线性能影响因素,以及天线结构参数、材料属性、电磁参数。
有源相控阵天线性能影响因素包括热参数、网格大小、结构位移提取、机电耦合模型;如图2所示,有源相控阵天线的结构参数包括天线中单元的行数M、列数N及单元间距dx和dy,印制板2、铝蜂窝上面板3、铝蜂窝4、铝蜂窝下面板5及T/R组件6的几何参数;微带贴片天线1在印制板2上方,T/R组件6在铝蜂窝下面板5下方;有源相控阵天线的材料属性包括弹性模量、泊松比、剪切模量、密度、导热系数和热膨胀系数;有源相控阵天线的电磁参数,包括天线的单元形式和中心工作频率f。
步骤2,建立有源相控阵天线有限元模型
确定有源相控阵天线单元的几何中心,并在几何中心处设置硬点;在ANSYS中建立有源相控阵天线有限元模型。
步骤3,计算有源相控阵天线温度场分布
根据有源相控阵天线的工作要求和有限元模型,施加边界条件,计算天线阵面温度场分布。
步骤4,计算有源相控阵天线热参数的置信度
根据有源相控阵天线的材料热属性,基于阵面温度场分布,计算有源相控阵天线在建模过程中热参数的置信度,公式如下:
式中,P1为热参数的置信度,Tideal为建模初始温度时其热膨胀系数,Treal为最终稳定状态下的热膨胀系数。
步骤5,计算网格大小的置信度
(5a)选取一组网格大小Mi(1≤i≤n),确定边界条件,并施加天线环境热载荷,得到天线阵面的温度场分布,计算天线阵面温度均方根误差RMSi,当RMSi+1与RMSi的绝对值误差时,则将第i+1组网格大小Mi+1和均方根误差RMSi+1作为标准,分别记为Mc和RMSc;
(5b)当选取网格大小为Md时,计算天线阵面温度均方根误差RMSd,得到网格划分的置信度,公式如下:
步骤6,计算有源相控阵天线结构热变形
将T/R组件热功耗作用下的天线温度场分布作为载荷,施加于天线有限元模型,确定结构约束,计算有源相控阵天线阵面的结构热变形。
步骤7,计算结构位移提取的置信度
根据有源相控阵天线的结构热变形,提取第(m,n)个天线单元几何中心的位置偏移量δmn,用天线单元几何中心的位移代替天线单元整体位移,计算结构位移提取的置信度P3,公式如下:
式中,为天线单元位置偏移量的均值。
步骤8,使用机电耦合模型,计算天线的电性能
(8a)设有源相控阵天线中共有M×N个天线单元,第(m,n)个天线单元的位置矢量为观察点P相对于坐标系O-xyz所在的方向(θ,φ)上的方位矢量为其中:
式中,xmn、ymn、zmn分别为位置矢量在坐标轴x、y、z方向的投影; 分别为坐标轴x、y、z方向的单位矢量;αx、αy、αz分别为位置矢量与坐标轴x、y、z方向的夹角;
(8b)当天线阵面产生变形时,第(m,n)个天线单元产生的位置偏移量为则变形后第(m,n)个天线单元的位置矢量为:
式中,Δxmn、Δymn、Δzmn分别为位置偏移量在坐标轴x、y、z方向的投影;
(8c)根据有源相控阵天线电磁叠加原理,理想情况下有源相控阵天线的方向图函数为:
结合(8b)得到变形有源相控阵天线方向图函数为:
式中,为有源相控阵天线的波常数,λ为有源相控阵天线工作波长,Ee为天线单元方向图函数,Amn是天线单元激励电流幅度,j为虚数单位,是天线单元激励电流相位;
(8d)根据有源相控阵天线方向图函数绘制的天线方向图,计算得到变形后有源相控阵天线的增益Gc。
步骤9,计算机电耦合模型的置信度
(9a)利用电磁分析软件HFSS计算有源相控阵天线的增益Gd;
(9b)将步骤(8d)中利用机电耦合模型计算得到的增益Gc和利用HFSS软件仿真得到的增益Gd进行对比,计算机电耦合模型的置信度P4,公式如下:
步骤10,基于层次分析法,确定加权系数
基于层次分析法,确定热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1、α2、α3、α4。
(10a)建立有源相控阵天线性能仿真置信度的层次结构模型,以有源相控阵天线性能仿真的置信度P为标准,分别比较热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4相对于有源相控阵天线性能仿真置信度P的重要性,从而构造出判断矩阵R;
(10b)设热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1,α2,α3,α4,记为W=[α1 α2 α3 α4];
(10c)由判断矩阵R计算判断矩阵R最大的特征值λ,再由公式RW=λW,求得W=[α1 α2 α3 α4]。
步骤11,确定有源相控阵天线性能仿真置信度
建立置信度计算公式P=α1·P1+α2·P2+α3·P3+α4·P4,计算并最终确定有源相控阵天线性能仿真的置信度。
本发明的优点可通过以下仿真实验进一步说明:
一、确定有源相控阵天线结构参数、热参数和电磁参数
本实例以工作频率为2.25GHZ的微带贴片天线,x方向阵元数M=5,y方向阵元数N=5,排列间距为dx×dy=60mm×60mm组成的有源相控阵天线为对象。其结构参数、材料属性如表1和表2所示,T/R组件的热功耗Q=5W。
表1有源相控阵天线的结构参数
表2有源相控阵天线的材料属性
二、计算有源相控阵天线性能仿真的置信度
1.计算天线温度场分布
研究对象为微带贴片天线,在天线单元的几何中心处建立硬点,根据有源相控阵天线的结构参数、材料属性,在ANSYS中采用SOLID278和SHELL131建立有源相控阵天线的有限元模型,如图3所示,施加边界条件,计算天线温度场分布,如图4所示。
2.计算有源相控阵天线热参数的置信度
根据材料热属性,基于阵面温度场分布,计算有源相控阵天线在建模过程中热参数的置信度,公式如下:
式中,P1为热参数的置信度,Tideal为建模初始温度时其热膨胀系数,Treal为最终稳定状态下的热膨胀系数,则Tideal=2.3e-5m/K,Tideal=2.41e-5m/K,由得出P1=0.954。
3.计算网格大小的置信度
3.1选取一组网格大小Mi(1≤i≤n),确定边界条件,并施加天线环境热载荷,得到天线阵面的温度分布,计算天线阵面温度均方根误差RMSi,当RMSi+1与RMSi的绝对值误差时,则将第i+1组网格大小Mi+1和均方根误差RMSi+1作为标准,记为Mc和RMSc,计算结果如表3。
表3网格大小及对应的均方根误差
由表3可知,则Mc=2,RMSc=5.04。
3.2当选取网格大小为Md时,计算天线阵面温度均方根误差RMSd,得到网格划分的置信度,公式如下:
划分网格时,当选取网格大小Md=4,则RMSd=5.32,由公式计算得P2=0.944。
4.计算天线结构热变形
将T/R组件热功耗作用下的天线温度场分布作为载荷,施加于天线有限元模型,确定结构约束,计算有源相控阵天线阵面的结构热变形,位移云图如图5所示。
5.计算结构位移提取的置信度
根据有源相控阵天线的结构热变形,提取第(m,n)(1≤m≤5,1≤n≤5)个天线单元几何中心的位置偏移量δmn,用天线单元几何中心的位移代替天线单元整体位移,计算结构位移提取的置信度P3,公式如下:
式中,为天线单元位置偏移量的均值,计算得到结构位移提取的置信度P3=0.948。
6.使用机电耦合模型,计算天线的电性能
6.1设有源相控阵天线中共有M×N个天线单元,第(m,n)个天线单元的位置矢量为观察点P相对于坐标系O-xyz所在的方向(θ,φ)上的方位矢量为其中:
式中,xmn、ymn、zmn分别为位置矢量在坐标轴x、y、z方向的投影; 分别为坐标轴x、y、z方向的单位矢量;αx、αy、αz分别为位置矢量与坐标轴x、y、z方向的夹角;
6.2当天线阵面产生变形时,第(m,n)个天线单元产生的位置偏移量为 则变形后第(m,n)个天线单元的位置矢量为:
式中,Δxmn、Δymn、Δzmn分别为位置偏移量在坐标轴x、y、z方向的投影;
6.3根据有源相控阵天线电磁叠加原理,理想情况下有源相控阵天线的方向图函数为:
结合6.2得到变形有源相控阵天线方向图函数为:
式中,为有源相控阵天线的波常数,λ为有源相控阵天线工作波长,Ee为天线单元方向图函数,Amn是天线单元激励电流幅度,是天线单元激励电流相位;
6.4根据有源相控阵天线方向图函数绘制的天线方向图,如图6和图7所示,计算得到变形后有源相控阵天线的增益Gc,增益Gc=14.087dB。
7.计算机电耦合模型的置信度
7.1利用电磁分析软件HFSS计算变形后有源相控阵天线的增益Gd,得Gd=14.168dB;
7.2将步骤6.4中利用机电耦合模型计算得到的增益Gc和利用HFSS软件仿真得到的增益Gd进行对比,计算机电耦合模型的置信度P4,公式如下:
代入数值计算得到机电耦合模型的置信度P4=0.992。
8.基于层次分析法,确定加权系数
8.1建立有源相控阵天线性能仿真置信度的层次结构模型,如图8所示,其中P代表有源相控阵天线性能仿真的置信度,P1表示热参数的置信度,P2表示网格划分大小的置信度,P3表示结构位移提取的置信度,P4表示机电耦合模型的置信度;
8.2根据层次结构模型两两比较P1~P4相对于P的重要性,从而构造出判断矩阵
8.3设热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1,α2,α3,α4,记为W=[α1 α2 α3 α4],由8.2中判断矩阵R计算判断矩阵R最大的特征值λ,再由公式RW=λW,求得则
9.确定有源相控阵天线性能仿真置信度
建立置信度计算公式P=α1·P1+α2·P2+α3·P3+α4·P4,由以上步骤可知,P1=0.954,P2=0.944,P3=0.948,P4=0.992,最终计算得到有源相控阵天线性能仿真的置信度P=0.967=96.7%。
Claims (10)
1.面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,包括下述步骤:
(1)根据有源相控阵天线的工作要求,确定有源相控阵天线性能的影响因素,以及天线结构参数、材料属性和电磁参数;
(2)根据有源相控阵天线的结构参数和材料属性,建立有源相控阵天线有限元模型;
(3)根据有源相控阵天线的工作要求和有限元模型,施加边界条件,计算天线阵面温度场分布;
(4)根据有源相控阵天线的材料热属性,基于阵面温度场分布,计算有源相控阵天线热参数的置信度P1;
(5)根据天线阵面温度场分布,基于阵面温度均方根误差RMS,计算网格划分大小的置信度P2;
(6)将阵面温度场分布作为载荷,施加到天线有限元模型,计算有源相控阵天线阵面的结构热变形;
(7)提取天线单元几何中心节点的位置偏移量,确定天线单元位置偏移量的均值,计算结构位移提取的置信度P3;
(8)根据有源相控阵天线的电磁工作参数以及单元几何中心的位置偏移量,使用有源相控阵天线机电耦合模型,计算天线的电性能;
(9)将机电耦合模型的计算结果与电磁仿真软件的计算结果进行对比,计算机电耦合模型的置信度P4;
(10)基于层次分析法,确定热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1、α2、α3、α4;
(11)建立置信度计算公式P=α1·P1+α2·P2+α3·P3+α4·P4,计算并最终确定有源相控阵天线性能仿真的置信度。
2.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(1)中,有源相控阵天线性能影响因素包括热参数、网格大小、结构位移提取和机电耦合模型;所述有源相控阵天线的结构参数,包括天线单元的行数M、列数N及单元间距dx和dy、印制板、铝蜂窝上面板、铝蜂窝、铝蜂窝下面板及T/R组件的几何参数;所述有源相控阵天线的材料属性包括弹性模量、泊松比、剪切模量、密度、导热系数和热膨胀系数;所述有源相控阵天线的电磁参数包括天线的单元形式和中心工作频率f。
3.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(2)中,确定有源相控阵天线单元的几何中心,并在几何中心处设置硬点;在ANSYS中建立有源相控阵天线有限元模型。
4.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(4)中,根据有源相控阵天线的材料热属性,基于阵面温度场分布,计算有源相控阵天线在建模过程中热参数的置信度,公式如下:
式中,P1为热参数的置信度,Tideal为建模初始温度时其热膨胀系数,Treal为最终稳定状态下的热膨胀系数。
5.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(5)按如下过程进行:
(5a)选取一组网格大小Mi(1≤i≤n),确定边界条件,并施加天线环境热载荷,得到天线阵面的温度场分布,计算天线阵面温度均方根误差RMSi,当RMSi+1与RMSi的绝对值误差时,则将第i+1组网格大小Mi+1和均方根误差RMSi+1作为标准,分别记为Mc和RMSc;
(5b)当选取网格大小为Md时,计算天线阵面温度均方根误差RMSd,得到网格划分大小的置信度,公式如下:
6.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(6)中,将T/R组件热功耗作用下的天线温度场分布作为载荷,施加于天线有限元模型,确定结构约束,计算有源相控阵天线阵面的结构热变形。
7.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(7)中,根据有源相控阵天线的结构热变形,提取第(m,n)个天线单元几何中心的位置偏移量δmn,用天线单元几何中心的位移代替天线单元整体位移,计算结构位移提取的置信度P3,公式如下:
式中,为天线单元位置偏移量的均值。
8.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(8)按如下过程进行:
(8a)设有源相控阵天线中共有M×N个天线单元,第(m,n)个天线单元的位置矢量为 观察点P相对于坐标系O-xyz所在的方向(θ,φ)上的方位矢量为 其中:
式中,xmn、ymn、zmn分别为位置矢量在坐标轴x、y、z方向的投影; 分别为坐标轴x、y、z方向的单位矢量;αx、αy、αz分别为位置矢量与坐标轴x、y、z方向的夹角;
(8b)当天线阵面产生变形时,第(m,n)个天线单元产生的位置偏移量为 则变形后第(m,n)个天线单元的位置矢量为:
式中,Δxmn、Δymn、Δzmn分别为位置偏移量在坐标轴x、y、z方向的投影;
(8c)根据有源相控阵天线电磁叠加原理,理想情况下有源相控阵天线的方向图函数为:
结合(8b)得到变形有源相控阵天线方向图函数为:
式中,为有源相控阵天线的波常数,λ为有源相控阵天线工作波长,Ee为天线单元方向图函数,Amn是天线单元激励电流幅度,j为虚数单位,是天线单元激励电流相位;
(8d)根据有源相控阵天线方向图函数绘制的天线方向图,计算得到变形后有源相控阵天线的增益Gc。
9.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(9)按如下过程进行:
(9a)利用电磁分析软件HFSS计算有源相控阵天线的增益Gd;
(9b)将步骤(8d)中利用机电耦合模型计算得到的增益Gc和利用HFSS软件仿真得到的增益Gd进行对比,计算机电耦合模型的置信度P4,公式如下:
10.根据权利要求1中所述的面向机电耦合的有源相控阵天线性能仿真置信度的计算方法,其特征在于,所述步骤(10)中,确定热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数,通过下述方法进行:
(10a)建立有源相控阵天线性能仿真置信度的层次结构模型,以有源相控阵天线性能仿真的置信度P为标准,分别比较热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4相对于有源相控阵天线性能仿真置信度P的重要性,从而构造出判断矩阵
(10b)设热参数的置信度P1、网格划分大小的置信度P2、结构位移提取的置信度P3和机电耦合模型的置信度P4的加权系数分别为α1,α2,α3,α4,记为W=[α1 α2 α3 α4];
(10c)由判断矩阵R计算判断矩阵R最大的特征值λ,再由公式RW=λW,求得W=[α1 α2 α3 α4]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610369219.8A CN106021764B (zh) | 2016-05-30 | 2016-05-30 | 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610369219.8A CN106021764B (zh) | 2016-05-30 | 2016-05-30 | 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021764A true CN106021764A (zh) | 2016-10-12 |
CN106021764B CN106021764B (zh) | 2019-04-02 |
Family
ID=57091430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610369219.8A Active CN106021764B (zh) | 2016-05-30 | 2016-05-30 | 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021764B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107918696A (zh) * | 2017-10-20 | 2018-04-17 | 西安电子科技大学 | 相控阵天线的多场耦合分析方法及计算机程序 |
CN108614916A (zh) * | 2018-03-29 | 2018-10-02 | 西安电子科技大学 | 一种快速补偿大型轻薄有源相控阵天线阵面变形的方法 |
CN109543334A (zh) * | 2018-12-02 | 2019-03-29 | 湖南大学 | 一种随机风载荷作用下相控阵雷达动态电性能分析方法 |
CN111123229A (zh) * | 2020-01-07 | 2020-05-08 | 浙江大学 | 基于电性能幅值加权的阵元安装位置的测量方法 |
CN111273240A (zh) * | 2020-02-21 | 2020-06-12 | 中国西安卫星测控中心 | 一种面向北斗三号立体四面相控阵天线测量基准确定方法 |
CN111428369A (zh) * | 2020-03-26 | 2020-07-17 | 中国人民解放军32035部队 | 一种空间目标碰撞预警结果置信度计算方法 |
CN112597703A (zh) * | 2020-12-22 | 2021-04-02 | 浙江大学 | 基于双分支深度神经网络的变尺寸阵列结构性能预测方法 |
CN117724049A (zh) * | 2024-02-07 | 2024-03-19 | 四川省华盾防务科技股份有限公司 | 相控阵系统阵元运行控制方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077485A (zh) * | 2014-06-30 | 2014-10-01 | 西安电子科技大学 | 基于吻合度的模型正确性的评价方法 |
CN104615836A (zh) * | 2015-02-12 | 2015-05-13 | 西安电子科技大学 | 星载微带阵列天线热变形对电性能影响的快速预测方法 |
-
2016
- 2016-05-30 CN CN201610369219.8A patent/CN106021764B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104077485A (zh) * | 2014-06-30 | 2014-10-01 | 西安电子科技大学 | 基于吻合度的模型正确性的评价方法 |
CN104615836A (zh) * | 2015-02-12 | 2015-05-13 | 西安电子科技大学 | 星载微带阵列天线热变形对电性能影响的快速预测方法 |
Non-Patent Citations (7)
Title |
---|
FANG KE 等: "AHP ultra weight for simulation credibility evaluation", 《JOURNAL OF BEIJING UNIVERSITY OF AERONAUTICS AND ASTRONAUTICS》 * |
刘聪锋 等: "自适应调零天线对抗效能层次分析评估方法", 《西安电子科技大学学报(自然科学版)》 * |
刘荷明 等: "层次分析法在仿真系统可信度评估中的应用", 《航天电子对抗》 * |
刘超: "面向机电热耦合的有源相控阵天线结构振动与热分析软件", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
王从思 等: "基于热—结构耦合的有源相控阵天线综合分析方法", 《2009年全国天线年会论文集》 * |
王国平 等: "仿真系统可信度分析", 《现代防御技术》 * |
郁正德 等: "层次分析法在雷达侦察机仿真试验置信度评估中的应用", 《情报指挥控制系统与仿真技术》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107918696A (zh) * | 2017-10-20 | 2018-04-17 | 西安电子科技大学 | 相控阵天线的多场耦合分析方法及计算机程序 |
CN107918696B (zh) * | 2017-10-20 | 2021-05-11 | 西安电子科技大学 | 相控阵天线的多场耦合分析方法及计算机程序 |
CN108614916A (zh) * | 2018-03-29 | 2018-10-02 | 西安电子科技大学 | 一种快速补偿大型轻薄有源相控阵天线阵面变形的方法 |
CN109543334A (zh) * | 2018-12-02 | 2019-03-29 | 湖南大学 | 一种随机风载荷作用下相控阵雷达动态电性能分析方法 |
CN111123229A (zh) * | 2020-01-07 | 2020-05-08 | 浙江大学 | 基于电性能幅值加权的阵元安装位置的测量方法 |
CN111123229B (zh) * | 2020-01-07 | 2021-12-07 | 浙江大学 | 基于电性能幅值加权的阵元安装位置的测量方法 |
CN111273240B (zh) * | 2020-02-21 | 2023-05-16 | 中国西安卫星测控中心 | 一种面向北斗三号立体四面相控阵天线测量基准确定方法 |
CN111273240A (zh) * | 2020-02-21 | 2020-06-12 | 中国西安卫星测控中心 | 一种面向北斗三号立体四面相控阵天线测量基准确定方法 |
CN111428369A (zh) * | 2020-03-26 | 2020-07-17 | 中国人民解放军32035部队 | 一种空间目标碰撞预警结果置信度计算方法 |
CN112597703A (zh) * | 2020-12-22 | 2021-04-02 | 浙江大学 | 基于双分支深度神经网络的变尺寸阵列结构性能预测方法 |
CN112597703B (zh) * | 2020-12-22 | 2022-08-23 | 浙江大学 | 基于双分支深度神经网络的变尺寸阵列结构性能预测方法 |
CN117724049A (zh) * | 2024-02-07 | 2024-03-19 | 四川省华盾防务科技股份有限公司 | 相控阵系统阵元运行控制方法、装置、设备及存储介质 |
CN117724049B (zh) * | 2024-02-07 | 2024-04-09 | 四川省华盾防务科技股份有限公司 | 相控阵系统阵元运行控制方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106021764B (zh) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021764A (zh) | 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法 | |
CN103353904B (zh) | 有源夹层微带天线与电磁综合的数据驱动设计方法及天线 | |
CN102890741B (zh) | 一种振动变形对阵列天线电性能影响的预测方法 | |
Dubuc et al. | Solution of the unsteady Euler equations using an implicit dual-time method | |
Le et al. | Effect of hydrofoil flexibility on the power extraction of a flapping tidal generator via two-and three-dimensional flow simulations | |
CN106940739B (zh) | 振动对机翼共形相控阵天线电性能影响的快速预测方法 | |
CN104036093B (zh) | 一种基于阵元互耦的大型变形阵列天线副瓣性能预测方法 | |
CN104615836B (zh) | 星载微带阵列天线热变形对电性能影响的快速预测方法 | |
CN107103124B (zh) | 基于机电耦合理论的变形阵列天线远场方向图分析方法 | |
CN102998540B (zh) | 共形承载微带天线阵面形貌对电性能影响的预测方法 | |
Yoon et al. | Effects of camber angle on aerodynamic performance of flapping-wing micro air vehicle | |
CN101344564A (zh) | 基于机电热三场耦合的有源相控阵天线电性能预测方法 | |
CN106446469B (zh) | 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 | |
Liu et al. | Discrete vortex method with flow separation corrections for flapping-foil power generators | |
CN105740502A (zh) | 一种载荷变形下的圆柱共形阵列天线电性能预测方法 | |
CN104036114A (zh) | 一种基于机电耦合的六边形有源相控阵天线结构公差的快速确定方法 | |
CN102788920B (zh) | 基于机电耦合模型的偏置反射面天线电性能预测方法 | |
Maheri et al. | WTAB, a computer program for predicting the performance of horizontal axis wind turbines with adaptive blades | |
Zhang et al. | Linearized Euler solver for rapid frequency-domain aeroelastic analysis | |
Li et al. | An efficient implementation of aeroelastic tailoring based on efficient computational fluid dynamics-based reduced order model | |
Shahverdi et al. | An efficient reduced-order modelling approach based on fluid eigenmodes and boundary element method | |
CN104063533B (zh) | 一种机翼蒙皮天线机电耦合分析方法 | |
CN106528954A (zh) | 空间网状天线温度载荷作用下机电集成优化设计方法 | |
CN109543334A (zh) | 一种随机风载荷作用下相控阵雷达动态电性能分析方法 | |
CN107590309A (zh) | 基于近似计算公式的网状天线电性能时域特性分析方法 |
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 |