CN106842951A - 面向电性能与控制的空间网状天线状态空间建模方法 - Google Patents
面向电性能与控制的空间网状天线状态空间建模方法 Download PDFInfo
- Publication number
- CN106842951A CN106842951A CN201710135148.XA CN201710135148A CN106842951A CN 106842951 A CN106842951 A CN 106842951A CN 201710135148 A CN201710135148 A CN 201710135148A CN 106842951 A CN106842951 A CN 106842951A
- Authority
- CN
- China
- Prior art keywords
- matrix
- modal
- unit
- subscript
- electrical property
- 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
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Aerials With Secondary Devices (AREA)
Abstract
本发明公开了一种面向电性能与控制的空间网状天线状态空间建模方法,包括:输入天线几何参数、材料参数、电参数与控制矩阵、控制向量;确定系统状态变量与系统输出量;计算理想天线远区电场;建立天线结构有限元模型;进行天线结构模态分析;输出固有频率矩阵、模态质量阵、模态振型阵;计算模态刚度阵、模态阻尼阵、模态输入矩阵;生成系统状态矩阵、系统控制矩阵;提取节点、单元与形函数信息;计算电性能对节点位移一阶、二阶系数矩阵;计算电性能对模态坐标一阶、二阶系数矩阵;生成系统一阶、二阶输出矩阵;建立系统状态空间模型。本发明克服了传统建模方法的不足,实现了同时面向电性能与控制的空间网状天线状态空间建模。
Description
技术领域
本发明属于雷达天线技术领域,具体涉及雷达天线领域中的一种面向电性能与控制的空间网状天线状态空间建模方法。
背景技术
网状天线由于其质量轻、收拢体积小等优点被逐渐应用于空间天线设计中。网状天线在轨运行过程中,周期性地受到太空辐射热、冲击等载荷的影响。在载荷影响下,天线形面发生柔性变形,导致电性能恶化。为了降低外部载荷对天线电性能的影响,一方面需要针对天线结构所受载荷进行分析并提出相应的设计策略,另一方面就是天线柔性结构进行精密主被动控制。对天线结构进行控制需要建立较为精确的天线状态空间模型;更进一步,为了实现电性能指标要求,其状态空间模型需要同时面向电性能与控制,也就是建立面向电性能与控制的天线结构状态空间模型,以此进行面向电性能的天线精密控制。
Gawaronski在文献“Advanced structural dynamics and active control ofstructures”(NewYork,Springer,2004年出版)中,公开了一种面向结构性能的天线结构状态空间模型,该方法虽然可以面向结构响应,但其对电性能的影响没有考虑在内,若要计算天线电性能,需要更进一步进行电性能分析。张洁、黄进、宋瑞雪、邱丽丽在中国专利“一种面向控制的大型天线建模方法”(授权号:ZL201310496650.5)中提出了一种面向指向精度与控制的大型天线动力学建模方法,该方法同样没有直接以天线电性能为系统输出,也无法在控制中考虑天线电性能。因此,针对空间网状天线的控制问题,需要一种面向电性能与控制的状态空间建模方法,能够克服传统建模方法无法直接建立与电性能联系的不足,实现同时面向电性能与控制的空间网状天线状态空间建模。
发明内容
本发明的目的是克服上述现有技术的不足,提供一种面向电性能与控制的空间网状天线状态空间建模方法,该方法从机械-电磁与控制集成的角度出发,以天线电性能为系统输出变量,进行空间网状天线面向电性能与控制的状态空间建模。
本发明的技术方案是:面向电性能与控制的空间网状天线状态空间建模方法,包括如下步骤:
(1)输入天线几何参数、材料参数、电参数与控制矩阵、控制向量
输入用户提供的空间网状天线的几何参数、材料参数、电参数与与控制矩阵、控制向量;其中几何参数包括口径、焦距、偏置距离、前后网面最小距离;材料参数包括索结构、桁架结构以及金属丝网结构的材料密度、横截面积、杨氏弹性模量、泊松比和瑞利阻尼系数;电参数包括工作波长、馈源参数、馈源初级方向图以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求;
(2)确定系统状态变量与系统输出量
2a)按照下式确定以模态坐标系下的模态坐标与模态坐标的导数为系统状态变量:
其中,x为系统状态变量,q表示模态坐标系下的模态坐标列向量,表示模态坐标导数列向量,上标·表示求导数运算;
2b)按照下式确定以天线远区电场为系统输出量:
y=E
其中,y表示系统输出量,E表示天线远区电场;
(3)计算理想天线远区电场
根据天线几何参数中的口径、焦距、偏置高度,电参数中的工作波长、馈源参数、馈源初级方向图,采用物理光学法计算理想天线远区电场;
(4)建立天线结构有限元模型
根据用户提供的天线几何参数、材料参数建立天线结构有限元模型,其中索结构采用只受拉的杆单元进行建模,桁架结构采用梁单元进行建模,金属丝网结构采用壳单元进行建模;
(5)进行天线结构模态分析
根据天线结构有限元模型,进行天线结构模态分析,确定天线结构各阶模态质量、各阶振动频率与对应的振型模态向量;
(6)输出固有频率矩阵、模态质量阵、模态振型阵
分别以天线结构模态分析得到的各阶模态质量、各阶振动频率为对角线元素,输出固有频率矩阵与模态质量阵,将天线结构各阶振动频率对应的振型模态向量按列放置在矩阵中,输出模态振型阵;
(7)计算模态刚度阵、模态阻尼阵、模态输入矩阵;
(8)生成系统状态矩阵、系统控制矩阵
8a)按照下式生成系统状态矩阵:
其中,Z表示天线结构的模态阻尼比矩阵,Mm为步骤(6)输出的模态质量阵,下标m表示模态,上标-1表示矩阵求逆运算,Dm为天线结构模态阻尼阵,Ω为步骤(6)输出的固有频率阵,A表示系统状态矩阵,I表示单位矩阵;
8b)按照下式生成系统控制矩阵:
其中,B为系统控制矩阵,Bm天线结构的模态输入矩阵,下标m表示模态;
(9)提取节点、单元与形函数信息
以建立的结构有限元模型为基础,提取有限元模型中处于电磁波照射下的反射面部分的节点、单元和形函数信息;
(10)计算电性能对节点位移一阶、二阶系数矩阵;
(11)计算电性能对模态坐标一阶、二阶系数矩阵;
(12)生成系统一阶、二阶输出矩阵
12a)按照下式生成系统一阶输出矩阵:
其中,C1表示系统一阶输出矩阵,下标1表示一阶,Gq为步骤(11)得到的电性能对模态坐标一阶系数矩阵,下标q表示模态坐标,上标T表示转置运算;
12b)按照下式生成系统二阶输出矩阵:
其中,C2表示系统二阶输出矩阵,下标2表示二阶,Hq为步骤(11)得到的电性能对模态坐标二阶系数矩阵,下标q表示模态坐标;
(13)建立系统状态空间模型
按照下式建立系统状态空间模型:
其中,x表示以模态坐标与模态坐标导数构成的系统状态变量,上标·表示求导数运算,A为步骤(8)得到的系统状态矩阵,B表示步骤(8)得到的系统控制矩阵,u表示步骤(1)输入的控制向量,y表示系统输出量,C1表示步骤(12)得到的系统一阶输出矩阵,下标1表示一阶,C2表示步骤(12)得到的系统二阶输出矩阵,下标2表示二阶,E0表示步骤(3)得到的理想天线远区电场。
上述步骤(3)中所述的物理光学法是一种基于面电流分布的高频近似方法,计算公式如下:
其中,E0表示理想天线远区电场,表示远场观察点位置矢量,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Σ表示反射曲面,表示反射面上位置矢量处的面电流密度,表示反射面位置矢量,表示远场观察点的单位矢量,σ表示投影口面,表示单位法向矢量,表示反射面位置矢量处的入射磁场。
上述步骤(7)中按照下式计算天线结构模态刚度阵:
Km=MmΩ2
其中,Km为天线结构模态刚度矩阵,下标m表示模态,Mm为步骤(6)输出的模态质量阵,Ω为步骤(6)输出的固有频率阵;
按照下式计算天线结构模态阻尼阵:
Dm=αKm+βMm
其中,Dm为天线结构模态阻尼阵,下标m表示模态,Km为天线结构模态刚度阵,Mm为步骤(6)输出的模态质量阵,α、β分别步骤(1)中输入的瑞利阻尼系数;
按照下式计算天线结构的模态输入矩阵:
其中,Bm为天线结构的模态输入矩阵,下标m表示模态,Mm为步骤(6)输出的模态质量阵,上标-1表示矩阵求逆运算,Φ为步骤(6)输出的模态振型阵,上标T表示矩阵转置运算,B0为步骤(1)输入的控制矩阵,下标0表示初始状态。
上述步骤(10)中通过下式计算电性能对节点位移一阶系数矩阵:
其中,表示单元e的电性能一阶系数矩阵,上标e表示从步骤(9)中提取的结构有限元模型中某一单元,下标i表示位于单元e上的节点编号,表示单元e的电性能一阶系数矩阵的第i个分量,符号∈表示从属关系,NUM表示单元e上的节点总数,表示单元e的法向矢量,表示反射面位置矢量处的入射磁场,表示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,k表示自由空间波数,表示远场观察点的单位矢量,Qi表示步骤(9)中提取的相对于第i个节点的形函数,θs表示位置矢量在馈源坐标系下的俯仰角,下标s表示馈源坐标系,θ表示远场观察点俯仰角,σe表示单元e在口径面内的投影面积;G表示电性能对节点位移一阶系数矩阵,m表示单元总数,Σ表示有限元组集运算;
通过下式计算电性能对节点位移二阶系数矩阵:
其中,表示单元e的电性能二阶系数矩阵,上标e表示从步骤(9)中提取的结构有限元模型中某一单元,u和v分别表示位于单元e上的节点编号,表示由节点u和v构成的单元e的电性能二阶系数矩阵分量,符号∈表示从属关系,NUM表示单元e上的节点总数,k表示自由空间波数,表示单元e的法向矢量,表示反射面位置矢量处的入射磁场,表示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,表示远场观察点的单位矢量,Qu表示步骤(9)中提取的相对于第u个节点的形函数,Qv表示步骤(9)中提取的相对于第v个节点的形函数,θs表示位置矢量在馈源坐标系下的俯仰角,下标s表示馈源坐标系,θ表示远场观察点俯仰角,σe表示单元e在口径面内的投影面积;H表示电性能对节点位移二阶系数矩阵,m表示单元总数,Σ表示有限元组集运算。
上述步骤(11)中通过下式计算电性能对模态坐标一阶系数矩阵:
其中,Gq为电性能对模态坐标一阶系数矩阵,下标q表示模态坐标,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Φ为步骤(6)输出的模态振型阵,上标T表示转置运算,G为步骤(10)得到的电性能对节点位移一阶系数矩阵;
通过下式计算电性能对模态坐标二阶系数矩阵:
其中,Hq为电性能对模态坐标二阶系数矩阵,下标q表示模态坐标,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Φ为步骤(6)输出的模态振型阵,上标T表示转置运算,H为步骤(10)得到的电性能对节点位移二阶系数矩阵。
本发明的有益效果:本发明首先输入天线几何参数、材料参数、电参数与控制矩阵、控制向量,确定天线系统状态变量与系统输出量,采用物理光学法计算理想天线远区电场;与此同时,根据几何参数、材料参数信息建立天线结构有限元模型;其次,在结构有限元模型的基础上,进行天线结构模态分析,输出固有频率矩阵、模态质量阵和模态振型阵,并依次计算模态刚度阵、模态阻尼阵与模态输入矩阵,生成系统状态矩阵与系统控制矩阵;再次,从结构有限元模型中提取节点、单元与形函数信息,计算电性能对节点位移的一阶、二阶系数矩阵,结合模态振型阵,获得电性能对模态坐标的一阶、二阶系数矩阵,生成系统一阶、二阶输出矩阵;最后,以系统状态矩阵、系统控制矩阵与系统一阶、二阶输出矩阵为基础,建立面向电性能与控制的系统状态空间模型。
与现有技术相比,本发明具有以下优点:
1.本发明从天线结构有限元模型出发,进行天线结构模态分析,将模态分析得到的模态振型阵与电性能对节点位移的一阶、二阶系数矩阵相结合,获得了电性能对模态坐标的一阶、二阶系数矩阵,从而在系统状态空间模型中实现了对天线电性能的考虑。
2.本发明以天线电性能为系统输出量,克服了传统状态空间模型建模方法无法直接建立与电性能联系的不足,实现了同时面向电性能与控制的空间网状天线状态空间建模。
以下将结合附图对本发明做进一步详细说明。
附图说明
图1为本发明的流程图;
图2为网状天线结构示意图;
图3为网状天线投影示意图;
图4为冲击激励作用下的电性能响应时间历程图。
具体实施方式
下面结合附图1,对本发明具体实施方式作进一步的详细描述:
本发明提供了一种面向电性能与控制的空间网状天线状态空间建模方法,包括如下步骤:
步骤1,输入用户提供的空间网状天线的几何参数、材料参数、电参数与与控制矩阵、控制向量;其中几何参数包括口径、焦距、偏置距离、前后网面最小距离;材料参数包括索结构、桁架结构和金属丝网结构的材料密度、横截面积、杨氏弹性模量、泊松比以及瑞利阻尼系数;电参数包括工作波长、馈源参数、馈源初级方向图以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求;
步骤2,确定系统状态变量与系统输出量
2a)按照下式确定以模态坐标系下的模态坐标与模态坐标的导数为系统状态变量:
其中,x为系统状态变量,q表示模态坐标系下的模态坐标列向量,表示模态坐标导数列向量,上标·表示求导数运算;
2b)按照下式确定以天线远区电场为系统输出量:
y=E
其中,y表示系统输出量,E表示天线远区电场;
步骤3,根据天线几何参数中的口径、焦距、偏置高度,电参数中的工作波长、馈源参数、馈源初级方向图,采用下式物理光学法计算理想天线远区电场:
其中,E0表示远区电场,表示远场观察点位置矢量,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Σ表示反射曲面,表示反射面上位置矢量处的面电流密度,表示反射面位置矢量,表示远场观察点的单位矢量,σ表示投影口面,表示单位法向矢量,表示反射面位置矢量处的入射磁场;
步骤4,根据用户提供的天线几何参数、材料参数建立天线结构有限元模型,其中索结构采用只受拉的杆单元进行建模,桁架结构采用梁单元进行建模,金属丝网结构采用壳单元进行建模;
步骤5,根据天线结构有限元模型,进行天线结构模态分析,确定天线结构各阶模态质量、各阶振动频率与对应的振型模态向量;
步骤6,分别以天线结构模态分析得到的各阶模态质量、各阶振动频率为对角线元素,输出固有频率矩阵与模态质量阵,将天线结构各阶振动频率对应的振型模态向量按列放置在矩阵中,输出模态振型阵;
步骤7,计算模态刚度阵、模态阻尼阵、模态输入矩阵
7a)按照下式计算天线结构模态刚度阵:
Km=MmΩ2
其中,Km为天线结构模态刚度矩阵,下标m表示模态,Mm为步骤6输出的模态质量阵,Ω为步骤6输出的固有频率阵;
7b)按照下式计算天线结构模态阻尼阵:
Dm=αKm+βMm
其中,Dm为天线结构模态阻尼阵,下标m表示模态,Km为天线结构模态刚度阵,Mm为步骤6输出的模态质量阵,α、β分别步骤1中输入的瑞利阻尼系数;
7c)按照下式计算天线结构的模态输入矩阵:
其中,Bm为天线结构的模态输入矩阵,下标m表示模态,Mm为步骤6输出的模态质量阵,上标-1表示矩阵求逆运算,Φ为步骤6输出的模态振型阵,上标T表示矩阵转置运算,B0为步骤1输入的控制矩阵,下标0表示初始状态;
步骤8,生成系统状态矩阵、系统控制矩阵
8a)按照下式生成系统状态矩阵:
其中,Z表示天线结构的模态阻尼比矩阵,Mm为步骤6输出的模态质量阵,下标m表示模态,上标-1表示矩阵求逆运算,Dm为天线结构模态阻尼阵,Ω为步骤6输出的固有频率阵,A表示系统状态矩阵,I表示单位矩阵;
8b)按照下式生成系统控制矩阵:
其中,B为系统控制矩阵,Bm天线结构的模态输入矩阵,下标m表示模态;
步骤9,以建立的结构有限元模型为基础,提取有限元模型中处于电磁波照射下的反射面部分的节点、单元和形函数信息;
步骤10,计算电性能对节点位移一阶、二阶系数矩阵
10a)通过下式计算电性能对节点位移一阶系数矩阵:
其中,表示单元e的电性能一阶系数矩阵,上标e表示从步骤9中提取的结构有限元模型中某一单元,下标i表示位于单元e上的节点编号,表示单元e的电性能一阶系数矩阵的第i个分量,符号∈表示从属关系,NUM表示单元e上的节点总数,表示单元e的法向矢量,表示反射面位置矢量处的入射磁场,表示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,k表示自由空间波数,表示远场观察点的单位矢量,Qi表示步骤9中提取的相对于第i个节点的形函数,θs表示位置矢量在馈源坐标系下的俯仰角,下标s表示馈源坐标系,θ表示远场观察点俯仰角,σe表示单元e在口径面内的投影面积;G表示电性能对节点位移一阶系数矩阵,m表示单元总数,Σ表示有限元组集运算;
10b)通过下式计算电性能对节点位移二阶系数矩阵:
其中,表示单元e的电性能二阶系数矩阵,上标e表示从步骤9中提取的结构有限元模型中某一单元,u和v分别表示位于单元e上的节点编号,表示由节点u和v构成的单元e的电性能二阶系数矩阵分量,符号∈表示从属关系,NUM表示单元e上的节点总数,k表示自由空间波数,表示单元e的法向矢量,表示反射面位置矢量处的入射磁场,表示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,表示远场观察点的单位矢量,Qu表示步骤9中提取的相对于第u个节点的形函数,Qv表示步骤9中提取的相对于第v个节点的形函数,θs表示位置矢量在馈源坐标系下的俯仰角,下标s表示馈源坐标系,θ表示远场观察点俯仰角,σe表示单元e在口径面内的投影面积;H表示电性能对节点位移二阶系数矩阵,m表示单元总数,Σ表示有限元组集运算;
步骤11,计算电性能对模态坐标一阶、二阶系数矩阵
11a)通过下式计算电性能对模态坐标一阶系数矩阵:
其中,Gq为电性能对模态坐标一阶系数矩阵,下标q表示模态坐标,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Φ为步骤6输出的模态振型阵,上标T表示转置运算,G为步骤10得到的电性能对节点位移一阶系数矩阵;
11b)通过下式计算电性能对模态坐标二阶系数矩阵:
其中,Hq为电性能对模态坐标二阶系数矩阵,下标q表示模态坐标,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Φ为步骤6输出的模态振型阵,上标T表示转置运算,H为步骤10得到的电性能对节点位移二阶系数矩阵;
步骤12,生成系统一阶、二阶输出矩阵
12a)按照下式生成系统一阶输出矩阵:
其中,C1表示系统一阶输出矩阵,下标1表示一阶,Gq为步骤11得到的电性能对模态坐标一阶系数矩阵,下标q表示模态坐标,上标T表示转置运算;
12b)按照下式生成系统二阶输出矩阵:
其中,C2表示系统二阶输出矩阵,下标2表示二阶,Hq为步骤11得到的电性能对模态坐标二阶系数矩阵,下标q表示模态坐标;
步骤13,按照下式建立系统状态空间模型:
其中,x表示以模态坐标与模态坐标导数构成的系统状态变量,上标·表示求导数运算,A为步骤8得到的系统状态矩阵,B表示步骤8得到的系统控制矩阵,u表示步骤1输入的控制向量,y表示系统输出量,C1表示步骤12得到的系统一阶输出矩阵,下标1表示一阶,C2表示步骤12得到的系统二阶输出矩阵,下标2表示二阶,E0表示步骤3得到的理想天线远区电场。
本发明的优点可通过以下仿真实验进一步说明:
1.仿真条件:
网状天线最大投影口径9.23m、焦距6m,偏置高度5m,前后网面最小间距0.2m,如图2所示。工作频率2GHz,馈源采用Cosine-Q类型馈源,极化方式为右旋圆极化,馈源参数为Qx=Qy=8.338,馈源倾斜角41.64度。网状反射面前网面在口径面内沿半径方向等分6段,如图3所示,其中虚线代表周边桁架,实线代表索网。
2.仿真结果:
根据网状天线参数,对天线结构上网面最高点施加大小为25N、方向与主轴z方向相反的冲击载荷。针对天线结构受到冲击载荷的电性能响应时间历程,分别采用传统由结构到电磁的间接建模方法与本发明面向电性能与控制的直接建模方法进行对比计算。图4分别为采用传统间接建模方法与本发明直接建模方法得到的电性能动力响应时间历程图。表1为最小主轴方向系数。结合图4与表1可以看出,本发明方法得到的电性能动力响应时间历程曲线与传统间接建模方法吻合,验证了方法的有效性。
表1本方法与准确模型比较表
综上所述,本发明首先输入天线几何参数、材料参数、电参数与控制矩阵、控制向量,确定天线系统状态变量与系统输出量,采用物理光学法计算理想天线远区电场;与此同时,根据几何参数、材料参数信息建立天线结构有限元模型;其次,在结构有限元模型的基础上,进行天线结构模态分析,输出固有频率矩阵、模态质量阵和模态振型阵,并依次计算模态刚度阵、模态阻尼阵与模态输入矩阵,生成系统状态矩阵与系统控制矩阵;再次,从结构有限元模型中提取节点、单元与形函数信息,计算电性能对节点位移的一阶、二阶系数矩阵,结合模态振型阵,获得电性能对模态坐标的一阶、二阶系数矩阵,生成系统一阶、二阶输出矩阵;最后,以系统状态矩阵、系统控制矩阵与系统一阶、二阶输出矩阵为基础,建立面向电性能与控制的系统状态空间模型。
与现有技术相比,本发明具有以下优点:
1.本发明从天线结构有限元模型出发,进行天线结构模态分析,将模态分析得到的模态振型阵与电性能对节点位移的一阶、二阶系数矩阵相结合,获得了电性能对模态坐标的一阶、二阶系数矩阵,从而在系统状态空间模型中实现了对天线电性能的考虑。
2.本发明以天线电性能为系统输出量,克服了传统状态空间模型建模方法无法直接建立与电性能联系的不足,实现了同时面向电性能与控制的空间网状天线状态空间建模。
本实施方式中没有详细叙述的部分属本行业的公知的常用手段,这里不一一叙述。以上例举仅仅是对本发明的举例说明,并不构成对本发明的保护范围的限制,凡是与本发明相同或相似的设计均属于本发明的保护范围之内。
Claims (5)
1.面向电性能与控制的空间网状天线状态空间建模方法,其特征在于,包括如下步骤:
(1)输入天线几何参数、材料参数、电参数与控制矩阵、控制向量
输入用户提供的空间网状天线的几何参数、材料参数、电参数与控制矩阵、控制向量;其中几何参数包括口径、焦距、偏置距离、前后网面最小距离;材料参数包括索结构、桁架结构以及金属丝网结构的材料密度、横截面积、杨氏弹性模量、泊松比和瑞利阻尼系数;电参数包括工作波长、馈源参数、馈源初级方向图以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求;
(2)确定系统状态变量与系统输出量
2a)按照下式确定以模态坐标系下的模态坐标与模态坐标的导数为系统状态变量:
其中,x为系统状态变量,q表示模态坐标系下的模态坐标列向量,表示模态坐标导数列向量,上标·表示求导数运算;
2b)按照下式确定以天线远区电场为系统输出量:
y=E
其中,y表示系统输出量,E表示天线远区电场;
(3)计算理想天线远区电场
根据天线几何参数中的口径、焦距、偏置高度,电参数中的工作波长、馈源参数、馈源初级方向图,采用物理光学法计算理想天线远区电场;
(4)建立天线结构有限元模型
根据用户提供的天线几何参数、材料参数建立天线结构有限元模型,其中索结构采用只受拉的杆单元进行建模,桁架结构采用梁单元进行建模,金属丝网结构采用壳单元进行建模;
(5)进行天线结构模态分析
根据天线结构有限元模型,进行天线结构模态分析,确定天线结构各阶模态质量、各阶振动频率与对应的振型模态向量;
(6)输出固有频率矩阵、模态质量阵、模态振型阵
分别以天线结构模态分析得到的各阶模态质量、各阶振动频率为对角线元素,输出固有频率矩阵与模态质量阵,将天线结构各阶振动频率对应的振型模态向量按列放置在矩阵中,输出模态振型阵;
(7)计算模态刚度阵、模态阻尼阵、模态输入矩阵;
(8)生成系统状态矩阵、系统控制矩阵
8a)按照下式生成系统状态矩阵:
其中,Z表示天线结构的模态阻尼比矩阵,Mm为步骤(6)输出的模态质量阵,下标m表示模态,上标-1表示矩阵求逆运算,Dm为天线结构模态阻尼阵,Ω为步骤(6)输出的固有频率阵,A表示系统状态矩阵,I表示单位矩阵;
8b)按照下式生成系统控制矩阵:
其中,B为系统控制矩阵,Bm天线结构的模态输入矩阵,下标m表示模态;
(9)提取节点、单元与形函数信息
以建立的结构有限元模型为基础,提取有限元模型中处于电磁波照射下的反射面部分的节点、单元和形函数信息;
(10)计算电性能对节点位移一阶、二阶系数矩阵;
(11)计算电性能对模态坐标一阶、二阶系数矩阵;
(12)生成系统一阶、二阶输出矩阵
12a)按照下式生成系统一阶输出矩阵:
其中,C1表示系统一阶输出矩阵,下标1表示一阶,Gq为步骤(11)得到的电性能对模态坐标一阶系数矩阵,下标q表示模态坐标,上标T表示转置运算;
12b)按照下式生成系统二阶输出矩阵:
其中,C2表示系统二阶输出矩阵,下标2表示二阶,Hq为步骤(11)得到的电性能对模态坐标二阶系数矩阵,下标q表示模态坐标;
(13)建立系统状态空间模型
按照下式建立系统状态空间模型:
其中,x表示以模态坐标与模态坐标导数构成的系统状态变量,上标·表示求导数运算,A为步骤(8)得到的系统状态矩阵,B表示步骤(8)得到的系统控制矩阵,u表示步骤(1)输入的控制向量,y表示系统输出量,C1表示步骤(12)得到的系统一阶输出矩阵,下标1表示一阶,C2表示步骤(12)得到的系统二阶输出矩阵,下标2表示二阶,E0表示步骤(3)得到的理想天线远区电场。
2.根据权利要求1所述的面向电性能与控制的空间网状天线状态空间建模方法,其特征在于:步骤(3)中所述的物理光学法是一种基于面电流分布的高频近似方法,计算公式如下:
其中,E0表示理想天线远区电场,表示远场观察点位置矢量,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Σ表示反射曲面,表示反射面上位置矢量处的面电流密度,表示反射面位置矢量,表示远场观察点的单位矢量,σ表示投影口面,表示单位法向矢量,表示反射面位置矢量处的入射磁场。
3.根据权利要求1所述的面向电性能与控制的空间网状天线状态空间建模方法,其特征在于:步骤(7)中按照下式计算天线结构模态刚度阵:
Km=MmΩ2
其中,Km为天线结构模态刚度矩阵,下标m表示模态,Mm为步骤(6)输出的模态质量阵,Ω为步骤(6)输出的固有频率阵;
按照下式计算天线结构模态阻尼阵:
Dm=αKm+βMm
其中,Dm为天线结构模态阻尼阵,下标m表示模态,Km为天线结构模态刚度阵,Mm为步骤(6)输出的模态质量阵,α、β分别步骤(1)中输入的瑞利阻尼系数;
按照下式计算天线结构的模态输入矩阵:
其中,Bm为天线结构的模态输入矩阵,下标m表示模态,Mm为步骤(6)输出的模态质量阵,上标-1表示矩阵求逆运算,Φ为步骤(6)输出的模态振型阵,上标T表示矩阵转置运算,B0为步骤(1)输入的控制矩阵,下标0表示初始状态。
4.根据权利要求1所述的面向电性能与控制的空间网状天线状态空间建模方法,其特征在于:步骤(10)中通过下式计算电性能对节点位移一阶系数矩阵:
其中,表示单元e的电性能一阶系数矩阵,上标e表示从步骤(9)中提取的结构有限元模型中某一单元,下标i表示位于单元e上的节点编号,表示单元e的电性能一阶系数矩阵的第i个分量,符号∈表示从属关系,NUM表示单元e上的节点总数,表示单元e的法向矢量,表示反射面位置矢量处的入射磁场,表示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,k表示自由空间波数,表示远场观察点的单位矢量,Qi表示步骤(9)中提取的相对于第i个节点的形函数,θs表示位置矢量在馈源坐标系下的俯仰角,下标s表示馈源坐标系,θ表示远场观察点俯仰角,σe表示单元e在口径面内的投影面积;G表示电性能对节点位移一阶系数矩阵,m表示单元总数,Σ表示有限元组集运算;
通过下式计算电性能对节点位移二阶系数矩阵:
其中,表示单元e的电性能二阶系数矩阵,上标e表示从步骤(9)中提取的结构有限元模型中某一单元,u和v分别表示位于单元e上的节点编号,表示由节点u和v构成的单元e的电性能二阶系数矩阵分量,符号∈表示从属关系,NUM表示单元e上的节点总数,k表示自由空间波数,表示单元e的法向矢量,表示反射面位置矢量处的入射磁场,表示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,表示远场观察点的单位矢量,Qu表示步骤(9)中提取的相对于第u个节点的形函数,Qv表示步骤(9)中提取的相对于第v个节点的形函数,θs表示位置矢量在馈源坐标系下的俯仰角,下标s表示馈源坐标系,θ表示远场观察点俯仰角,σe表示单元e在口径面内的投影面积;H表示电性能对节点位移二阶系数矩阵,m表示单元总数,Σ表示有限元组集运算。
5.根据权利要求1所述的面向电性能与控制的空间网状天线状态空间建模方法,其特征在于:步骤(11)中通过下式计算电性能对模态坐标一阶系数矩阵:
其中,Gq为电性能对模态坐标一阶系数矩阵,下标q表示模态坐标,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Φ为步骤(6)输出的模态振型阵,上标T表示转置运算,G为步骤(10)得到的电性能对节点位移一阶系数矩阵;
通过下式计算电性能对模态坐标二阶系数矩阵:
其中,Hq为电性能对模态坐标二阶系数矩阵,下标q表示模态坐标,j表示虚数单位,k表示自由空间波数,↑表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,π表示圆周率,表示单位并矢,表示单位矢量的并矢,Φ为步骤(6)输出的模态振型阵,上标T表示转置运算,H为步骤(10)得到的电性能对节点位移二阶系数矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710135148.XA CN106842951B (zh) | 2017-03-08 | 2017-03-08 | 面向电性能与控制的空间网状天线状态空间建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710135148.XA CN106842951B (zh) | 2017-03-08 | 2017-03-08 | 面向电性能与控制的空间网状天线状态空间建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842951A true CN106842951A (zh) | 2017-06-13 |
CN106842951B CN106842951B (zh) | 2020-02-07 |
Family
ID=59143419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710135148.XA Active CN106842951B (zh) | 2017-03-08 | 2017-03-08 | 面向电性能与控制的空间网状天线状态空间建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842951B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108628345A (zh) * | 2018-06-28 | 2018-10-09 | 北京航空航天大学 | 一种电磁航天器编队悬停协同控制方法及系统 |
CN109472066A (zh) * | 2018-10-23 | 2019-03-15 | 西安电子科技大学 | 基于单元中心点位移的反射面天线随机误差分析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104112051A (zh) * | 2014-07-25 | 2014-10-22 | 西安电子科技大学 | 基于单元形函数的反射面天线机电集成设计方法 |
CN105844018A (zh) * | 2016-03-23 | 2016-08-10 | 西安电子科技大学 | 一种大型反射面天线反射体俯仰关键模态的选取方法 |
CN106021743A (zh) * | 2016-05-24 | 2016-10-12 | 西安电子科技大学 | 基于结构-电磁集成灵敏度信息的星载网状天线波束赋形设计方法 |
CN106354910A (zh) * | 2016-08-22 | 2017-01-25 | 西安电子科技大学 | 一种面向有源相控阵天线的结构主模态确定方法 |
CN106354906A (zh) * | 2016-08-19 | 2017-01-25 | 西安电子科技大学 | 一种基于三角形面片误差的网状可展开天线找形方法 |
CN106446469A (zh) * | 2016-11-14 | 2017-02-22 | 西安电子科技大学 | 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 |
-
2017
- 2017-03-08 CN CN201710135148.XA patent/CN106842951B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104112051A (zh) * | 2014-07-25 | 2014-10-22 | 西安电子科技大学 | 基于单元形函数的反射面天线机电集成设计方法 |
CN105844018A (zh) * | 2016-03-23 | 2016-08-10 | 西安电子科技大学 | 一种大型反射面天线反射体俯仰关键模态的选取方法 |
CN106021743A (zh) * | 2016-05-24 | 2016-10-12 | 西安电子科技大学 | 基于结构-电磁集成灵敏度信息的星载网状天线波束赋形设计方法 |
CN106354906A (zh) * | 2016-08-19 | 2017-01-25 | 西安电子科技大学 | 一种基于三角形面片误差的网状可展开天线找形方法 |
CN106354910A (zh) * | 2016-08-22 | 2017-01-25 | 西安电子科技大学 | 一种面向有源相控阵天线的结构主模态确定方法 |
CN106446469A (zh) * | 2016-11-14 | 2017-02-22 | 西安电子科技大学 | 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108628345A (zh) * | 2018-06-28 | 2018-10-09 | 北京航空航天大学 | 一种电磁航天器编队悬停协同控制方法及系统 |
CN108628345B (zh) * | 2018-06-28 | 2020-03-27 | 北京航空航天大学 | 一种电磁航天器编队悬停协同控制方法及系统 |
CN109472066A (zh) * | 2018-10-23 | 2019-03-15 | 西安电子科技大学 | 基于单元中心点位移的反射面天线随机误差分析方法 |
CN109472066B (zh) * | 2018-10-23 | 2023-04-18 | 西安电子科技大学 | 基于单元中心点位移的反射面天线随机误差分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106842951B (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106446469B (zh) | 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 | |
Denli et al. | Structural-acoustic optimization of sandwich structures with cellular cores for minimum sound radiation | |
CN102203782B (zh) | 求解拉格朗日形式的欧拉方程的数值方法 | |
CN104392034B (zh) | 一种基于l1/2范数的稀疏线性阵列优化方法 | |
CN107038299A (zh) | 一种考虑互耦效应的变形阵列天线远场方向图补偿方法 | |
Sarkar et al. | Wave farm modelling of oscillating wave surge converters | |
CN106940739A (zh) | 一种振动对机翼共形相控阵天线电性能影响的快速预测方法 | |
Noad et al. | Optimisation of arrays of flap-type oscillating wave surge converters | |
CN104615836B (zh) | 星载微带阵列天线热变形对电性能影响的快速预测方法 | |
CN106599504B (zh) | 基于机电耦合模型的空间网状天线动力载荷分析方法 | |
CN105160115B (zh) | 基于近似与灵敏度分析的反射面天线机电集成优化设计方法 | |
CN104112051B (zh) | 基于单元形函数的反射面天线机电集成设计方法 | |
CN107657070A (zh) | 基于遗传算法频率分集共形阵列波束的解耦合方法及天线 | |
CN107103124B (zh) | 基于机电耦合理论的变形阵列天线远场方向图分析方法 | |
CN105302962A (zh) | 基于结构-电磁混合单元的反射面天线机电集成优化设计方法 | |
CN106021743A (zh) | 基于结构-电磁集成灵敏度信息的星载网状天线波束赋形设计方法 | |
CN106842951A (zh) | 面向电性能与控制的空间网状天线状态空间建模方法 | |
Neshat et al. | New insights into position optimisation of wave energy converters using hybrid local search | |
CN104166764B (zh) | 面向多层次性能的导弹总体多学科设计方法 | |
CN104701637B (zh) | 基于机电耦合与最小二乘法的变形阵列天线电性能补偿方法 | |
La Magna | Bending-active plates: strategies for the induction of curvature through the means of elastic bending of plate-based structures | |
CN106202657A (zh) | 冲击波载荷作用下平面阵列天线的电性能预测方法 | |
Larsen et al. | Wake modeling and simulation | |
CN107515956A (zh) | 一种基于hfss单元法的大型有限平面阵列分析方法 | |
CN104794289A (zh) | 一种扩展直角坐标系下完全匹配吸收边界的实现方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |