CN106446469A - 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 - Google Patents

基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 Download PDF

Info

Publication number
CN106446469A
CN106446469A CN201610997422.XA CN201610997422A CN106446469A CN 106446469 A CN106446469 A CN 106446469A CN 201610997422 A CN201610997422 A CN 201610997422A CN 106446469 A CN106446469 A CN 106446469A
Authority
CN
China
Prior art keywords
represent
antenna
unit
electrical property
order
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
Application number
CN201610997422.XA
Other languages
English (en)
Other versions
CN106446469B (zh
Inventor
张树新
杜敬利
张岳震
宋立伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian Univ
Original Assignee
Xidian Univ
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian Univ filed Critical Xidian Univ
Priority to CN201610997422.XA priority Critical patent/CN106446469B/zh
Publication of CN106446469A publication Critical patent/CN106446469A/zh
Application granted granted Critical
Publication of CN106446469B publication Critical patent/CN106446469B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Abstract

本发明公开了一种基于二阶近似计算公式的空间网状天线电性能动力响应分析方法,包括:输入天线几何参数、材料参数与电参数;计算理想天线远区电场;建立天线结构有限元模型;提取节点、单元与形函数信息;计算单元电性能一阶、二阶系数矩阵;组集总体电性能一阶、二阶系数矩阵;进行天线结构模态分析;输出模态矩阵;计算电性能动力响应一阶导数列向量、二阶Hessian阵;施加动力载荷;计算模态坐标;近似计算天线远区电场变化量;计算远区电场动力响应;判断电性能是否满足要求;若满足输出天线结构设计方案;否则更新天线参数。本发明在保证计算精度的前提下可有效减少电性能动力响应计算时间,实现天线结构机电集成优化设计。

Description

基于二阶近似计算公式的空间网状天线电性能动力响应分析 方法
技术领域
[0001] 本发明属于雷达天线技术领域,具体涉及雷达天线领域中的一种基于二阶近似计 算公式的空间网状天线电性能动力响应分析方法。
背景技术
[0002] 网状天线由于其质量轻、收拢体积小等优点被逐渐应用于空间天线设计中。网状 天线在轨运行过程中,周期性地受到太空辐射热、冲击等动力载荷的影响,动力载荷引起的 天线表面变形对天线的电性能产生很大影响。动力载荷带来的天线表面变形将引起天线结 构发生振动,电性能增益下降、副瓣电平上升,波束倾斜,严重影响天线的工作性能。与静力 载荷不同,动力载荷随时间变化,其对天线电性能的影响程度也是随时间变化。因此,有必 要针对空间网状天线在轨运行受到的动力载荷进行分析,提出一种分析动力载荷作用下的 空间网状天线电性能动力响应计算公式,并进行天线结构机电集成优化设计。
[0003] 王从思等人在中国专利"一种振动变形对阵列天线电性能影响的预测方法"中提 出了一种分析振动变形对阵列天线电性能影响的预测方法。该方法以阵列天线为对象,分 析了随机振动对阵列天线电性能的影响;虽然该分析方法可以为空间网状天线提供借鉴, 但由于天线电性能计算的非线性特点,对于动力载荷作用下的天线电性能响应,该分析方 法出现计算量大的问题。周金柱、宋立伟等在文献"动载荷对结构功能一体化天线力电性能 的影响,机械工程学报,第52卷第9期,2016年5月,105-115"分析了动力载荷对结构功能一 体化天线的力电性能的影响,该分析方法同样难以解决计算量大,费时的问题。因此,针对 空间网状天线受动力载荷影响下的电性能动力响应问题,提出了一种基于二阶近似计算公 式的电性能动力响应分析方法,在保证计算精度的前提下,可以减少动力响应计算时间,并 以此进行天线结构机电集成优化设计。
发明内容
[0004] 本发明的目的是克服上述现有技术的不足,提供一种基于二阶近似计算公式的空 间网状天线电性能动力响应分析方法。该方法基于二阶近似计算公式,分析动力载荷作用 下的空间网状天线电性能动力响应特性,实现快速电性能动力响应计算,并进行天线结构 机电集成优化设计。
[0005] 本发明的技术方案是:基于二阶近似计算公式的空间网状天线电性能动力响应分 析方法,包括如下步骤:
[0006] (1)输入天线几何参数、材料参数与电参数
[0007] 输入用户提供的空间网状天线的几何参数、材料参数以及电参数;其中几何参数 包括口径、焦距、偏置距离、前后网面最小距离;材料参数包括索结构、桁架结构和金属丝网 结构的材料密度、横截面积、杨氏弹性模量、泊松比;电参数包括工作波长、馈源参数、馈源 初级方向图以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求;
[0008] (2)计算理想天线远区电场
[0009] 根据天线几何参数中的口径、焦距、偏置高度,电参数中的工作波长、馈源参数、馈 源初级方向图,采用物理光学法计算理想天线远区电场;
[0010] (3)建立天线结构有限元模型
[0011] 根据用户提供的天线几何参数、材料参数建立天线结构有限元模型,其中索结构 采用只受拉的杆单元进行建模,桁架结构采用梁单元进行建模,金属丝网结构采用壳单元 进行建模;
[0012] ⑷提取节点、单元与形函数信息
[0013] 以建立的结构有限元模型为基础,提取有限元模型中处于电磁波照射下的反射面 部分的节点、单元和形函数信息;
[0014] (5)计算单元电性能一阶、二阶系数矩阵
[0015] 5a)通过下式计算单元电性能一阶系数矩阵:
[0016] ^ = {η:157;:2,ieNUM
[0017] = |2/V x H(r)Qxp( jkr · R)jkQi{cos6s +cos^)c/a,
[0018] 其中,Af表示单元e的电性能一阶系数矩阵,上标e表示从步骤⑷中提取的结构有 限元模型中某一单元,下标i表示位于单元e上的节点编号表示单元e的电性能一阶系 数矩阵的第i个分量,符号e表示从属关系,NUM表示单元e上的节点总数,#表示单元e的法 向矢量,及表示反射面位置矢量F处的入射磁场,_?表示反射面位置矢量,exp表示自然对 数的指数运算,j表示虚数单位,k表示自由空间波数,左表示远场观察点的单位矢量,Qi表示 步骤⑷中提取的相对于第i个节点的形函数,9 5表示位置矢量?在馈源坐标系下的俯仰角, 下标s表示馈源坐标系,Θ表示远场观察点俯仰角,表示单元e在口径面内的投影面积;
[0019] 5b)通过下式计算单元电性能二阶系数矩阵:
[0020] h\ = ,…X:_uv},託 NUM,v e ΝϋΜ:
[0021] f^nr --k1^2Nx H(r)QKip(jkr · R)QnQv(cos0s + cos&fdae
[0022] 其中,/?表示单元e的电性能二阶系数矩阵,上标e表示从步骤(4)中提取的结构有 限元模型中某一单元,u和v分别表示位于单元e上的节点编号,表示由节点u和v构成的 单元e的电性能二阶系数矩阵分量,符号e表示从属关系,NUM表示单元e上的节点总数,k表 示自由空间波数,#表示单元e的法向矢量,#(?)表示反射面位置矢量F处的入射磁场,尹表 示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,及表示远场观察点的单 位矢量,Q u表示步骤(4)中提取的相对于第u个节点的形函数,Qv表示步骤(4)中提取的相对 于第v个节点的形函数,9 5表示位置矢量?在馈源坐标系下的俯仰角,下标s表示馈源坐标 系,Θ表示远场观察点俯仰角,表示单元e在口径面内的投影面积;
[0023] ⑶组集总体电性能一阶、二阶系数矩阵
[0024] 6a)通过下式组集总体电性能一阶系数矩阵: m
[0025] G=AM e=l
[0026] 其中,G表示总体电性能一阶系数矩阵,埒表示单元e的电性能一阶系数矩阵,上标 e表示从步骤(4)中提取的结构有限元模型中某一单元,m表示单元总数,A表示有限元组集 运算;
[0027] 6b)通过下式组集总体电性能二阶系数矩阵: m
[0028] M = A K e=l "
[0029] 其中,Η表示总体电性能二阶系数矩阵,Μ表示单元e的电性能二阶系数矩阵,上标 e表示从步骤(4)中提取的结构有限元模型中某一单元,m表示单元总数,A表示有限元组集 运算;
[0030] ⑵进行天线结构模态分析
[0031] 根据天线结构有限元模型,进行天线结构模态分析,确定天线结构无阻尼自由振 动的各阶振动频率与对应的阵型模态向量;
[0032] ⑶输出模态矩阵
[0033] 将天线结构的振型模态向量按列存放在模态矩阵中,并将其输出;
[0034] (9)计算电性能动力响应一阶导数列向量、二阶Hessian阵
[0035] 9a)通过下式计算电性能动力响应一阶导数列向量:
[0036] Gd=〇T * G
[0037] 其中,Gd为电性能动力响应一阶导数列向量,下标d表示动力响应,Φ为步骤⑶输 出的模态矩阵,上标T表示转置运算,G为步骤(6)得到的总体电性能一阶系数矩阵;
[0038] 9b)通过下式计算电性能动力响应二阶Hessian阵:
[0039] Η<ί=Φτ·Η·Φ
[0040] 其中,Hd为电性能动力响应二阶Hessian阵,下标d表示动力响应,Φ为步骤(8)输 出的模态矩阵,上标T表示转置运算,Η为步骤(6)得到的总体电性能二阶系数矩阵;
[0041] (10)施加动力载荷
[0042] 针对天线结构有限元模型施加随时间变化的动力载荷;
[0043] (11)计算模态坐标
[0044] 采用Duhamel积分公式,求解步骤(10)动力载荷作用下的天线结构动力响应的模 态坐标;
[0045] (12)近似计算天线远区电场变化量
[0046]结合步骤(11)获得的动力响应模态坐标与步骤⑶获得的电性能动力响应一阶导 数列向量、二阶Hessian阵,采用二阶近似计算公式计算天线远区电场变化量,通过下式计 算:
Figure CN106446469AD00081
[0048] 其中,At表示载荷作用下的远区电场变化量,j表示虚数单位,k表示自由空间波 数,η表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度, π表示圆周率,ϊ表示单位并矢,Μ表示单位矢量左的并矢,q (t)表示步骤(11)获得的动力 响应模态坐标,t为时间因子,Gd表示步骤(9)获得的电性能动力响应一阶导数列向量,上标 T表示转置运算,下标d表示动力响应,Hd表示步骤(9)获得的电性能动力响应二阶Hessian 阵;
[0049] (13)计算远区电场动力响应
[0050] 在步骤(2)和(12)的基础上,叠加理想天线远区电场和近似计算得到的远区电场 变化量,通过下式计算远区电场动力响应:
[0051] E = E0 + AE
[0052] 其中,启表示远区电场动力响应,At表示步骤(12)近似计算得到的远区电场变化 量,屋 &表示步骤(2)得到的理想天线远区电场。
[0053] (14)判断电性能是否满足要求
[0054] 判断步骤(13)得到的天线远区电场是否满足天线增益、波瓣宽度、副瓣电平、指向 精度在内的电性能要求,如果满足要求则转至步骤(15),否则转至步骤(16);
[0055] (15)输出天线结构设计方案
[0056] 当天线远区电场满足天线电性能要求时,输出天线结构设计数据;
[0057] (16)更新天线参数
[0058] 当天线远区电场不满足天线电性能要求时,更新天线参数,转至步骤(1)。
[0059]上述步骤(2)所述的物理光学法是一种基于面电流分布的高频近似方法,计算公 式如下:
Figure CN106446469AD00091
[0062] 其中,表示远区电场,f表示远场观察点位置矢量,j表示虚数单位,k表示自由 空间波数,η表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢 量幅度, π表示圆周率,_F表示单位并矢,Μ表示单位矢量及的并矢,Σ表示反射曲面, 表示反射面上位置矢量F处的面电流密度,i表示反射面位置矢量,及表示远场观察点的单 位矢量,σ表示投影口面,:/i表示单位法向矢量,//(7)表示反射面位置矢量/处的入射磁场。 [00 63] 上述步骤(11)所述的Duhamel积分公式是一种计算结构动力响应模态坐标的公 式,计算公式如下:
Figure CN106446469AD00092
L〇〇65j 其中,qi⑴表不糸统模态坐标列冋量中的弟i个模态坐标分量,t为时丨日」因于,c〇i 为第i阶自由振动的固有频率,Ci为第i阶模态振型阻尼比,Fi (τ)为第i阶模态动力载荷函 数,exp表示自然对数的指数运算,τ为积分变量,(1τ表示积分变量τ的微分。
[0066] 本发明的有益效果:本发明首先输入天线几何参数、材料参数与电参数信息,采用 物理光学法计算理想天线远区电场,与此同时,根据几何参数、材料参数信息建立天线结构 有限元模型;其次,根据天线结构有限元模型,提取节点、单元与形函数信息,计算单元电性 能一阶、二阶系数矩阵,并组集总体电性能一阶、二阶系数矩阵;在天线结构有限元模型基 础上,进行天线结构模态分析,输出模态矩阵;再次,结合总体电性能一阶、二阶系数矩阵与 模态矩阵,计算电性能动力响应一阶导数列向量、二阶Hessian阵;然后,对结构有限元模型 施加动力载荷,计算模态坐标;最后,将模态坐标代入电性能动力响应的二阶计算公式中, 通过叠加理想天线远区电场获得动力载荷作用下的天线电性能动力响应,以此进行天线结 构机电集成优化设计。
[0067] 与现有技术相比,本发明具有以下优点:
[0068] 1.本发明从天线结构有限元模型出发,进行天线结构模态分析,将模态分析得到 的模态矩阵与电性能一阶、二阶系数矩阵进行结合,获得了电性能动力响应的一阶导数列 向量与二阶Hessian阵,使电性能计算中包含了结构模态信息。
[0069] 2.本发明基于二阶近似计算公式获得空间网状天线电性能动力响应,在保证计算 精度的前提下,减少计算时间,提高计算效率,实现电性能动力响应的快速计算。
[0070] 以下将结合附图对本发明做进一步详细说明。
附图说明
[0071] 图1为本发明的流程图;
[0072]图2为网状天线结构示意图;
[0073]图3为网状天线投影示意图;
[0074]图4为网状天线施加动力载荷示意图;
[0075] 图5为动力载荷时间历程图;
[0076] 图6为电性能动力响应时间历程图。
具体实施方式
[0077]下面结合附图1,对本发明具体实施方式作进一步的详细描述:
[0078] 本发明提供了一种基于二阶近似计算公式的空间网状天线电性能动力响应分析 方法,包括如下步骤:
[0079] 步骤1,输入用户提供的空间网状天线的几何参数、材料参数以及电参数;其中几 何参数包括口径、焦距、偏置距离、前后网面最小距离;材料参数包括索结构、桁架结构和金 属丝网结构的材料密度、横截面积、杨氏弹性模量、泊松比;电参数包括工作波长、馈源参 数、馈源初级方向图以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求;
[0080] 步骤2,根据天线几何参数中的口径、焦距、偏置高度,电参数中的工作波长、馈源 参数、馈源初级方向图,采用下式物理光学法计算理想天线远区电场:
Figure CN106446469AD00101
[0083] 其中,氧表示远区电场,及表示远场观察点位置矢量,j表示虚数单位,k表示自由 空间波数,η表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢 量幅度J表示圆周率,F表示单位并矢,遍表示单位矢量龙的并矢,Σ表示反射曲面,J(?) 表示反射面上位置矢量f处的面电流密度,#表示反射面位置矢量,i表示远场观察点的单 位矢量,σ表示投影口面,3表示单位法向矢量,/7(/ :)表示反射面位置矢量F处的入射磁场;
[0084] 步骤3,根据用户提供的天线几何参数、材料参数建立天线结构有限元模型,其中 索结构采用只受拉的杆单元进行建模,桁架结构采用梁单元进行建模,金属丝网结构采用 壳单元进行建模;
[0085] 步骤4,以建立的结构有限元模型为基础,提取有限元模型中处于电磁波照射下的 反射面部分的节点、单元和形函数信息;
[0086] 步骤5,计算单元电性能一阶、二阶系数矩阵
[0087] 5a)通过下式计算单元电性能一阶系数矩阵:
Figure CN106446469AD00111
[0090] 其中,表示单元e的电性能一阶系数矩阵,上标e表示从步骤4中提取的结构有 限元模型中某一单元,下标i表示位于单元e上的节点编号,@丨表示单元e的电性能一阶系 数矩阵的第i个分量,符号e表示从属关系,NUM表示单元e上的节点总数,及表示单元e的法 向矢量,/7(/ 7)表示反射面位置矢量F处的入射磁场,?表示反射面位置矢量,exp表示自然对 数的指数运算,j表示虚数单位,k表示自由空间波数,及表示远场观察点的单位矢量,Qi表示 步骤4中提取的相对于第i个节点的形函数,0 3表示位置矢量?在馈源坐标系下的俯仰角,下 标s表示馈源坐标系,Θ表示远场观察点俯仰角,表示单元e在口径面内的投影面积;
[0091] 5b)通过下式计算单元电性能二阶系数矩阵:
Figure CN106446469AD00112
[0094] 其中,表示单元e的电性能二阶系数矩阵,上标e表示从步骤4中提取的结构有 限元模型中某一单元,u和v分别表示位于单元e上的节点编号,表示由节点u和v构成的 单元e的电性能二阶系数矩阵分量,符号e表示从属关系,NUM表示单元e上的节点总数,k表 示自由空间波数,及表示单元e的法向矢量,否(F)表示反射面位置矢量P处的入射磁场,?表 示反射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,身表示远场观察点的单 位矢量,Q u表示步骤4中提取的相对于第u个节点的形函数,Qv表示步骤4中提取的相对于第v 个节点的形函数,9S表示位置矢量F在馈源坐标系下的俯仰角,下标s表示馈源坐标系,Θ表 示远场观察点俯仰角,表示单元e在口径面内的投影面积;
[0095] 步骤6,组集总体电性能一阶、二阶系数矩阵
[0096] 6a)通过下式组集总体电性能一阶系数矩阵: m
[0097] G = AM B=1
[0098] 其中,G表示总体电性能一阶系数矩阵,¥表示单元e的电性能一阶系数矩阵,上标 e表示从步骤4中提取的结构有限元模型中某一单元,m表示单元总数,A表示有限元组集运 算;
[0099] 6b)通过下式组集总体电性能二阶系数矩阵: m
[0100] - \ e=l
[0101] 其中,Η表示总体电性能二阶系数矩阵,表示单元e的电性能二阶系数矩阵,上标 e表示从步骤4中提取的结构有限元模型中某一单元,m表示单元总数,A表示有限元组集运 算;
[0102] 步骤7,根据天线结构有限元模型,进行天线结构模态分析,确定天线
[0103] 结构无阻尼自由振动的各阶振动频率与对应的阵型模态向量;
[0104] 步骤8,将天线结构的振型模态向量按列存放在模态矩阵中,并将其输出;
[0Ί05] 步骤9,计算电性能动力响应一阶导数列向量、二阶Hessian阵
[0106] 9a)通过下式计算电性能动力响应一阶导数列向量:
[0107] Gd=〇T · G
[0108] 其中,Gd为电性能动力响应一阶导数列向量,下标d表示动力响应,Φ为步骤8输出 的模态矩阵,上标T表示转置运算,G为步骤6得到的总体电性能一阶系数矩阵;
[0109] 9b)通过下式计算电性能动力响应二阶Hessian阵:
[0110] Η<ί=Φτ·Η·Φ
[0111] 其中,Hd为电性能动力响应二阶Hessian阵,下标d表示动力响应,Φ为步骤8输出 的模态矩阵,上标T表示转置运算,Η为步骤6得到的总体电性能二阶系数矩阵;
[0112] 步骤10,针对天线结构有限元模型施加随时间变化的动力载荷;
[0113] 步骤11,采用下式Duhamel积分公式,求解步骤10动力载荷作用下的天线结构动力 响应的模态坐标:
Figure CN106446469AD00121
[0115] 其中,qi(t)表示系统模态坐标列向量中的第i个模态坐标分量,t为时间因子, 为第i阶自由振动的固有频率,Ci为第i阶模态振型阻尼比,Fi (τ)为第i阶模态动力载荷函 数,exp表示自然对数的指数运算,τ为积分变量,(1τ表示积分变量τ的微分;
[0116] 步骤12,近似计算天线远区电场变化量
[0117]结合步骤11获得的动力响应模态坐标与步骤9获得的电性能动力响应一阶导数列 向量、二阶Hessian阵,采用二阶近似计算公式计算天线远区电场变化量,通过下式计算:
Figure CN106446469AD00122
[0119] 其中,^表示载荷作用下的远区电场变化量,j表示虚数单位,k表示自由空间波 数,η表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度, π表示圆周率,f表示单位并矢,顧表示单位矢量及的并矢,q⑴表示步骤11获得的动力响 应模态坐标,t为时间因子,Gd表示步骤9获得的电性能动力响应一阶导数列向量,上标T表 示转置运算,下标d表示动力响应,Hd表示步骤9获得的电性能动力响应二阶Hessian阵; [0120] 步骤13,在步骤2和12的基础上,叠加理想天线远区电场和近似计算得到的远区电 场变化量,通过下式计算远区电场动力响应:
[0121] E = E{) + AE
[0122] 其中,J表示远区电场动力响应,At表示步骤12近似计算得到的远区电场变化量, 氧表示步骤2得到的理想天线远区电场。
[0123] 步骤14,判断步骤13得到的天线远区电场是否满足天线增益、波瓣宽度、副瓣电 平、指向精度在内的电性能要求,如果满足要求则转至步骤15,否则转至步骤16;
[0124] 步骤15,当天线远区电场满足天线电性能要求时,输出天线结构设计数据;
[0125] 步骤16,当天线远区电场不满足天线电性能要求时,更新天线参数,转至步骤1。
[0126] 本发明的优点可通过以下仿真实验进一步说明:
[0127] i.仿真条件:
[0128] 网状天线最大投影口径9.23m、焦距6m,偏置高度5m,前后网面最小间距0.2m,如图 2所示。工作频率2GHz,馈源采用Cosine-Q类型馈源,极化方式为右旋圆极化,馈源参数为Q x = Qy = 8.338,馈源倾斜角41.64度。网状反射面前网面在口径面内沿半径方向等分6段,如 图3所示,其中虚线代表周边桁架,实线代表索网。对天线结构上网面上端施加沿负z轴方向 的冲击载荷,冲击载荷施加位置如图4所;^。冲击载荷大小为25N,在0到0.02s内施加冲击载 荷,之后撤去该冲击载荷,计算天线电性能动力相应。冲击载荷作用的时间历程如图5所示。
[0129] 2.仿真结果:
[0130] 采用本发明的方法进行网状天线处于冲击载荷作用下的电性能动力响应分析,并 与传统方法进行比较。图6为采用传统方法与本发明方法得到的电性能动力响应时间历程 图。表1为最小主轴方向系数与计算时间表。从图6与表1中可以看出,在冲击载荷作用下,天 线主轴方向系数发生了振荡,振荡幅值逐渐减小;本发明方法得到的电性能动力响应时间 历程曲线与传统方法计算结果吻合,保证了计算精度;同时,本发明方法计算时间远远小于 传统方法计算时间,提高了计算效率。
[0131] 表1本方法与传统方法比较表
Figure CN106446469AD00131
[0134]综上所述,本发明首先输入天线几何参数、材料参数与电参数信息,采用物理光学 法计算理想天线远区电场,与此同时,根据几何参数、材料参数信息建立天线结构有限元模 型;其次,根据天线结构有限元模型,提取节点、单元与形函数信息,计算单元电性能一阶、 二阶系数矩阵,并组集总体电性能一阶、二阶系数矩阵;在天线结构有限元模型基础上,进 行天线结构模态分析,输出模态矩阵;再次,结合总体电性能一阶、二阶系数矩阵与模态矩 阵,计算电性能动力响应一阶导数列向量、二阶Hessian阵;然后,对结构有限元模型施加动 力载荷,计算模态坐标;最后,将模态坐标代入电性能动力响应的二阶计算公式中,通过叠 加理想天线远区电场获得动力载荷作用下的天线电性能动力响应,以此进行天线结构机电 集成优化设计。
[0135] 与现有技术相比,本发明具有以下优点:
[0136] 1.本发明从天线结构有限元模型出发,进行天线结构模态分析,将模态分析得到 的模态矩阵与电性能一阶、二阶系数矩阵进行结合,获得了电性能动力响应的一阶导数列 向量与二阶Hessian阵,使电性能计算中包含了结构模态信息。
[0137] 2.本发明基于二阶近似计算公式获得空间网状天线电性能动力响应,在保证计算 精度的前提下,减少计算时间,提高计算效率,实现电性能动力响应的快速计算。
[0138] 本实施方式中没有详细叙述的部分属本行业的公知的常用手段,这里不一一叙 述。以上例举仅仅是对本发明的举例说明,并不构成对本发明的保护范围的限制,凡是与本 发明相同或相似的设计均属于本发明的保护范围之内。

Claims (3)

1.基于二阶近似计算公式的空间网状天线电性能动力响应分析方法,其特征在于,包 括如下步骤: (1)输入天线几何参数、材料参数与电参数 输入用户提供的空间网状天线的几何参数、材料参数以及电参数;其中几何参数包括 口径、焦距、偏置距离、前后网面最小距离;材料参数包括索结构、桁架结构和金属丝网结构 的材料密度、横截面积、杨氏弹性模量、泊松比;电参数包括工作波长、馈源参数、馈源初级 方向图以及包括天线增益、波瓣宽度、副瓣电平、指向精度在内的电性能要求; ⑵计算理想天线远区电场 根据天线几何参数中的口径、焦距、偏置高度,电参数中的工作波长、馈源参数、馈源初 级方向图,采用物理光学法计算理想天线远区电场; ⑶建立天线结构有限元模型 根据用户提供的天线几何参数、材料参数建立天线结构有限元模型,其中索结构采用 只受拉的杆单元进行建模,桁架结构采用梁单元进行建模,金属丝网结构采用壳单元进行 建模; ⑷提取节点、单元与形函数信息 以建立的结构有限元模型为基础,提取有限元模型中处于电磁波照射下的反射面部分 的节点、单元和形函数信息; ⑶计算单元电性能一阶、二阶系数矩阵 .¾)诵]寸下宣里TT由丨牛能一阶畜翁拓陈.
Figure CN106446469AC00021
其中,\'表示单元e的电性能一阶系数矩阵,上标e表示从步骤⑷中提取的结构有限元 模型中某一单元,下标i表示位于单元e上的节点编号,ί;:表示单元e的电性能一阶系数矩 阵的第i个分量,符号e表示从属关系,NUM表示单元e上的节点总数,#表示单元e的法向矢 量,/?(/")表示反射面位置矢量:f处的入射磁场,P表示反射面位置矢量,exp表示自然对数的 指数运算,j表示虚数单位,k表示自由空间波数,表示远场观察点的单位矢量,Qi表示步骤 ⑷中提取的相对于第i个节点的形函数,9 S表示位置矢量F在馈源坐标系下的俯仰角,下标 s表示馈源坐标系,Θ表示远场观察点俯仰角,h表示单元e在口径面内的投影面积; 诵讨下忒i+笪里元电件能二阶系教钽陈,
Figure CN106446469AC00022
其中,表示单元e的电性能二阶系数矩阵,上标e表示从步骤⑷中提取的结构有限元 模型中某一单元,u和V分别表示位于单元e上的节点编号,表示由节点u和V构成的单元 e的电性能二阶系数矩阵分量,符号e表示从属关系,NUM表示单元e上的节点总数,k表示自 由空间波数,没表示单元e的法向矢量,/7(7)表示反射面位置矢量:?处的入射磁场,?表示反 射面位置矢量,exp表示自然对数的指数运算,j表示虚数单位,及表示远场观察点的单位矢 量,Qu表示步骤⑷中提取的相对于第u个节点的形函数,Q v表示步骤(4)中提取的相对于第V 个节点的形函数,9S表示位置矢量F在馈源坐标系下的俯仰角,下标s表示馈源坐标系,Θ表 示远场观察点俯仰角,O f3表示单元e在口径面内的投影面积; (6)组集总体电性能一阶、二阶系数矩阵 6a)通过下式组集总体电性能一阶系数矩阵:
Figure CN106446469AC00031
其中,G表示总体电性能一阶系数矩阵,表示单元e的电性能一阶系数矩阵,上标e表 示从步骤(4)中提取的结构有限元模型中某一单元,m表示单元总数,A表示有限元组集运 算; fih)诵ΐ寸下式拍隹益彼电性能二阶系数矩阵:
Figure CN106446469AC00032
共中,H表不思体电ΐ生能二阶系数矩阵,表示单元e的电性能二阶系数矩阵,上标e表 示从步骤(4)中提取的结构有限元模型中某一单元,m表示单元总数,A表示有限元组集运 算; ⑵进行天线结构模态分析 根据天线结构有限元模型,进行天线结构模态分析,确定天线结构无阻尼自由振动的 各阶振动频率与对应的阵型模态向量; ⑶输出模态矩阵 将天线结构的振型模态向量按列存放在模态矩阵中,并将其输出; (9) 计算电性能动力响应一阶导数列向量、二阶Hessian阵 9a)通过下式计算电性能动力响应一阶导数列向量: Gd= Φτ · G 其中,Gd为电性能动力响应一阶导数列向量,下标d表示动力响应,Φ为步骤⑶输出的 模态矩阵,上标T表示转置运算,G为步骤(6)得到的总体电性能一阶系数矩阵; 9b)通过下式计算电性能动力响应二阶Hessian阵: Hd= Φτ · H · Φ 其中,Hd为电性能动力响应二阶Hessian阵,下标d表示动力响应,Φ为步骤⑶输出的模 态矩阵,上标T表示转置运算,H为步骤(6)得到的总体电性能二阶系数矩阵; (10) 施加动力载荷 针对天线结构有限元模型施加随时间变化的动力载荷; (11) 计算模态坐标 采用Duhamel积分公式,求解步骤(10)动力载荷作用下的天线结构动力响应的模态坐 标; (12) 近似计算天线远区电场变化量 结合步骤(11)获得的动力响应模态坐标与步骤(9)获得的电性能动力响应一阶导数列 向量、二阶Hessian阵,采用二阶近似计算公式计算天线远区电场变化量,通过下式计算:
Figure CN106446469AC00041
其中,Δ及表示载荷作用下的远区电场变化量,j表示虚数单位,k表示自由空间波数,η 表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅度,表 示圆周率,J表示单位并矢,M表示单位矢量左的并矢,q CO表示步骤(11)获得的动力响应 模态坐标,t为时间因子,Gd表示步骤(9)获得的电性能动力响应一阶导数列向量,上标τ表 示转置运算,下标d表示动力响应,Hd表示步骤(9)获得的电性能动力响应二阶Hessian阵; (13) 计算远区电场动力响应 在步骤(2)和(12)的基础上,叠加理想天线远区电场和近似计算得到的远区电场变化 量,通过下式计算远区电场动力响应:
Figure CN106446469AC00042
其中,f表示远区电场动力响应,化表示步骤(12)近似计算得到的远区电场变化量, ^^表示步骤⑵得到的理想天线远区电场; (14) 判断电性能是否满足要求 判断步骤(13)得到的天线远区电场是否满足天线增益、波瓣宽度、副瓣电平、指向精度 在内的电性能要求,如果满足要求则转至步骤(15),否则转至步骤(16); (15) 输出天线结构设计方案 当天线远区电场满足天线电性能要求时,输出天线结构设计数据; (16) 更新天线参数 当天线远区电场不满足天线电性能要求时,更新天线参数,转至步骤(1)。
2. 根据权利要求1所述的基于二阶近似计算公式的空间网状天线电性能动力响应分析 方法,其特征在于:步骤(2)所述的物理光学法是一种基于面电流分布的高频近似方法,计 算公式如下:
Figure CN106446469AC00043
其中,爲表示远区电场,i表示远场观察点位置矢量,j表示虚数单位,k表示自由空间 波数,η表示自由空间波阻抗,exp表示自然对数的指数运算,R表示远场观察点位置矢量幅 度,表示圆周率,f表示单位并矢,M表示单位矢量1的并矢,Σ表示反射曲面,J(F)表示 反射面上位置矢量/处的面电流密度,及表示反射面位置矢量,左表示远场观察点的单位矢 量,σ表示投影口面,3表示单位法向矢量,/7(F)表示反射面位置矢量F处的入射磁场。
3. 根据权利要求1所述的基于二阶近似计算公式的空间网状天线电性能动力响应分析 方法,其特征在于:步骤(11)所述的Duhamel积分公式是一种计算结构动力响应模态坐标的 公式,计算公式如下:
Figure CN106446469AC00051
其中,qi (t)表不系统模态坐标列向量中的第i个模态坐标分量,t为时间因子,ω i为第i 阶自由振动的固有频率,Ii为第i阶模态振型阻尼比,Fi (τ)为第i阶模态动力载荷函数,exp 表示自然对数的指数运算,τ为积分变量,dT表示积分变量τ的微分。
CN201610997422.XA 2016-11-14 2016-11-14 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法 Active CN106446469B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610997422.XA CN106446469B (zh) 2016-11-14 2016-11-14 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610997422.XA CN106446469B (zh) 2016-11-14 2016-11-14 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法

Publications (2)

Publication Number Publication Date
CN106446469A true CN106446469A (zh) 2017-02-22
CN106446469B CN106446469B (zh) 2019-02-01

Family

ID=58207654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610997422.XA Active CN106446469B (zh) 2016-11-14 2016-11-14 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法

Country Status (1)

Country Link
CN (1) CN106446469B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842951A (zh) * 2017-03-08 2017-06-13 西安电子科技大学 面向电性能与控制的空间网状天线状态空间建模方法
CN107169174A (zh) * 2017-04-24 2017-09-15 西安电子科技大学 最优焦距下伞状天线机电集成优化设计方法
CN107590309A (zh) * 2017-08-03 2018-01-16 西安电子科技大学 基于近似计算公式的网状天线电性能时域特性分析方法
CN112526229A (zh) * 2020-11-30 2021-03-19 西北工业大学 展开臂关节间隙对可展开网状天线电性能影响的评估方法
CN113391138A (zh) * 2020-03-13 2021-09-14 中国人民解放军63756部队 基于跟踪轨迹拟合的天线副瓣识别及其自动转主瓣方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120050732A1 (en) * 2010-08-25 2012-03-01 Weixing Lu Sensor system with plasmonic nano-antenna array
CN104036093A (zh) * 2014-06-26 2014-09-10 西安电子科技大学 一种基于阵元互耦的大型变形阵列天线副瓣性能预测方法
CN104112051A (zh) * 2014-07-25 2014-10-22 西安电子科技大学 基于单元形函数的反射面天线机电集成设计方法
CN105160115A (zh) * 2015-09-14 2015-12-16 西安电子科技大学 基于近似与灵敏度分析的反射面天线机电集成优化设计方法
CN105302962A (zh) * 2015-10-30 2016-02-03 西安电子科技大学 基于结构-电磁混合单元的反射面天线机电集成优化设计方法
CN106021766A (zh) * 2016-05-30 2016-10-12 西安电子科技大学 基于机电集成优化的圆极化网状天线高指向精度设计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120050732A1 (en) * 2010-08-25 2012-03-01 Weixing Lu Sensor system with plasmonic nano-antenna array
CN104036093A (zh) * 2014-06-26 2014-09-10 西安电子科技大学 一种基于阵元互耦的大型变形阵列天线副瓣性能预测方法
CN104112051A (zh) * 2014-07-25 2014-10-22 西安电子科技大学 基于单元形函数的反射面天线机电集成设计方法
CN105160115A (zh) * 2015-09-14 2015-12-16 西安电子科技大学 基于近似与灵敏度分析的反射面天线机电集成优化设计方法
CN105302962A (zh) * 2015-10-30 2016-02-03 西安电子科技大学 基于结构-电磁混合单元的反射面天线机电集成优化设计方法
CN106021766A (zh) * 2016-05-30 2016-10-12 西安电子科技大学 基于机电集成优化的圆极化网状天线高指向精度设计方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842951A (zh) * 2017-03-08 2017-06-13 西安电子科技大学 面向电性能与控制的空间网状天线状态空间建模方法
CN106842951B (zh) * 2017-03-08 2020-02-07 西安电子科技大学 面向电性能与控制的空间网状天线状态空间建模方法
CN107169174A (zh) * 2017-04-24 2017-09-15 西安电子科技大学 最优焦距下伞状天线机电集成优化设计方法
CN107169174B (zh) * 2017-04-24 2019-11-05 西安电子科技大学 最优焦距下伞状天线机电集成优化设计方法
CN107590309A (zh) * 2017-08-03 2018-01-16 西安电子科技大学 基于近似计算公式的网状天线电性能时域特性分析方法
CN107590309B (zh) * 2017-08-03 2020-04-24 西安电子科技大学 基于近似计算公式的网状天线电性能时域特性分析方法
CN113391138A (zh) * 2020-03-13 2021-09-14 中国人民解放军63756部队 基于跟踪轨迹拟合的天线副瓣识别及其自动转主瓣方法
CN113391138B (zh) * 2020-03-13 2022-08-30 中国人民解放军63756部队 基于跟踪轨迹拟合的天线副瓣识别及其自动转主瓣方法
CN112526229A (zh) * 2020-11-30 2021-03-19 西北工业大学 展开臂关节间隙对可展开网状天线电性能影响的评估方法
CN112526229B (zh) * 2020-11-30 2021-09-14 西北工业大学 展开臂关节间隙对可展开网状天线电性能影响的评估方法

Also Published As

Publication number Publication date
CN106446469B (zh) 2019-02-01

Similar Documents

Publication Publication Date Title
CN106446469A (zh) 基于二阶近似计算公式的空间网状天线电性能动力响应分析方法
CN104112051B (zh) 基于单元形函数的反射面天线机电集成设计方法
US9399981B2 (en) Methods and systems for comparing vertical axis turbine arrays and providing configurations thereof
CN104732243A (zh) 基于cnn的sar目标识别方法
CN108738290B (zh) 一种超材料吸波器及设计方法
CN105160115B (zh) 基于近似与灵敏度分析的反射面天线机电集成优化设计方法
Grilli et al. Progress in fully nonlinear potential flow modeling of 3D extreme ocean waves
Tay et al. Hydrodynamic interactions of oscillating wave surge converters in an array under random sea state
CN106940739A (zh) 一种振动对机翼共形相控阵天线电性能影响的快速预测方法
CN106021764A (zh) 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法
CN102156764A (zh) 一种分析天线辐射和电磁散射的多分辨预条件方法
CN105302962A (zh) 基于结构-电磁混合单元的反射面天线机电集成优化设计方法
Zheng et al. An efficient method to reduce the numerical dispersion in the ADI-FDTD
CN108959772A (zh) 一种大规模有限周期阵列结构特征模式分析方法
CN106599504A (zh) 基于机电耦合模型的空间网状天线动力载荷分析方法
CN107240780A (zh) 基于面片积分公式的伞状天线结构优化设计方法
CN106021766A (zh) 基于机电集成优化的圆极化网状天线高指向精度设计方法
CN102298658B (zh) 一种行波管注波互作用的模拟方法
Ding et al. Effects of installation position of fin-shaped rods on wind-induced vibration and energy harvesting of aeroelastic energy converter
CN106934097A (zh) 面向电性能的空间网状天线关键动力学模态选取方法
CN106528954A (zh) 空间网状天线温度载荷作用下机电集成优化设计方法
Gardonio et al. Vibration energy harvesting from an array of flexible stalks exposed to airflow: A theoretical study
Tong et al. Multilevel fast multipole algorithm for acoustic wave scattering by truncated ground with trenches
CN113158527B (zh) 一种基于隐式fvfd计算频域电磁场的方法
CN106326579A (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