CN110470916A - 日照温度场作用下反射面天线电性能快速评估与补偿方法 - Google Patents

日照温度场作用下反射面天线电性能快速评估与补偿方法 Download PDF

Info

Publication number
CN110470916A
CN110470916A CN201910759994.8A CN201910759994A CN110470916A CN 110470916 A CN110470916 A CN 110470916A CN 201910759994 A CN201910759994 A CN 201910759994A CN 110470916 A CN110470916 A CN 110470916A
Authority
CN
China
Prior art keywords
matrix
temperature
reflecting surface
reflector antenna
finite element
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.)
Pending
Application number
CN201910759994.8A
Other languages
English (en)
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.)
Xian University of Electronic Science and Technology
Original Assignee
Xian University of Electronic Science and Technology
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 Xian University of Electronic Science and Technology filed Critical Xian University of Electronic Science and Technology
Priority to CN201910759994.8A priority Critical patent/CN110470916A/zh
Publication of CN110470916A publication Critical patent/CN110470916A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/10Radiation diagrams of antennas

Abstract

本发明公开了一种日照温度场作用下反射面天线电性能快速评估与补偿方法,包括:依据反射面天线结构参数以及材料属性建立有限元模型;确定温度传感器在有限元模型中的位置;建立温度到热变形的映射矩阵、扰动矩阵和理想矩阵;建立扰动矩阵系数和理想矩阵系数;调用映射矩阵、扰动矩阵、理想矩阵、扰动矩阵系数和理想矩阵系数,依据热电矩阵耦合模型和实测温度快速评估日照温度场作用下的反射面天线的电性能;依据评估结果实时调整反射面天线的方位角和俯仰角,实现日照温度场作用下的反射面天线指向的快速补偿。本发明可快速准确地由测量温度评估相应的电性能,进而实时改善日照温度场作用下反射面天线的指向精度。

Description

日照温度场作用下反射面天线电性能快速评估与补偿方法
技术领域
本发明涉及天线技术领域,具体是一种日照温度场作用下的反射面天线电性能快速评估与补偿方法,应用于对工程中反射面天线在日照温度场作用下电性能的快速评估与补偿。
背景技术
大型反射面天线因其结构简单,具有高增益、高指向精度的特点,其在射电天文、深空探测等领域得到广泛应用。随着科学技术的发展,反射面天线向着更大口径、更高频段的方向发展,任何微小结构变形都将引起电性能的显著下降,温度作为大型反射面天线的外载荷之一,其对大型反射面天线电性能的影响越来越受到重视。
西班牙IRAM-30米射电望远镜工作频率高达230GHz,该天线整体配有温控设备,通过将背架结构整体封闭隔离,并安装风机空调设备,保持天线结构的温度均匀性,显著降低了温度载荷对电性能的影响。美国Haystack36.6米天线配备一个直径45.7米的球形天线罩,天线罩内部安装空调设备,保障温度的均匀性,可显著降低温度载荷的影响。然而,像上海佘山65米天线、德国Effelsberg-100米天线,美国GBT-100米天线,以及中国拟建造的新疆110米天线等,因其结构庞大,配置温控设备将导致成本极其高昂,故工程中大口径天线大多露天工作,故日照温度场对大型反射面天线电性能的影响不可忽视。
为了实现露天工作下反射面天线电性能快速评估与补偿,工程中多布置大量温度传感器,通过温度采集数据,重构天线整体温度场,借助有限元模型,仿真计算天线的热变形;进而,建立具有热变形的反射面电磁分析模型,利用商用电磁分析软件仿真评估热变形影响下反射面天线的电性能,根据反射面天线的电性能变化,通过调整天线的馈源或者副面改善天线的指向精度。然而,上述电性能评估过程较为繁琐,且耗时较长,难以满足工程上电性能调整的实时性要求。工程中亟需解决日照温度场作用下反射面天线电性能快速评估与补偿问题。
发明内容
本发明的目的在于克服上述现有方法的不足,将上述繁琐复杂的电性能计算过程,简化为一系列矩阵相乘过程,将电性能评估转换为一次矩阵运算,提出了一种日照温度场作用下反射面天线电性能快速评估与补偿方法,可快速准确地由测量温度评估相应的电性能,进而基于电性能快速评估结果实时改善日照温度场作用下反射面天线的指向精度。
本发明是通过下述技术方案来实现的。
一种日照温度场作用下反射面天线电性能快速评估与补偿方法,包括如下步骤:
(1)依据反射面天线结构参数以及材料属性建立有限元模型;
(2)依据实际天线中温度传感器布局确定其在有限元模型中的位置;
(3)依据温度传感器在有限元模型中的位置建立温度到热变形的映射矩阵;
(4)依据有限元模型和天线辐射场积分公式建立扰动矩阵和理想矩阵;
(5)依据映射矩阵和温度采集数据建立扰动矩阵系数和理想矩阵系数;
(6)依据实测温度和热电矩阵耦合模型快速评估日照温度场作用下的反射面天线的电性能;
(7)依据电性能评估结果实时调整反射面天线的方位角和俯仰角,实现日照温度场作用下的反射面天线指向的快速补偿。
进一步,所述步骤(1)中,反射面天线结构参数包括天线背架、中心体、主反射面和副反射面;所述材料属性包括密度、泊松比、弹性模量和热膨胀系数;建立有限元模型且将反射面划分为三角形网格单元。
进一步,所述步骤(2)中,依据实际天线中温度传感器布局确定其在有限元模型中的位置,即确定温度传感器在有限元模型中的节点坐标(us,vs,ws),其中s=1,2,…,Sensor,s表示第s个温度传感器,Sensor为温度传感器总数。
进一步,所述步骤(3)按如下过程进行:
(3a)假设第s个温度传感器测量数据为1℃,其余温度传感器测量数据为0℃,基于空间距离平方反比插值方法计算有限元模型所有节点的插值温度列向量为Ts,依次令s等于1,2,…,Sensor,可依次计算得到相应的插值温度列向量将其组成插值矩阵
(3b)假设第n个有限元节点温度为1℃,其余有限元节点温度为0℃,根据有限元模型,设置参考温度为0℃,计算所有反射面节点轴向热变形列向量为δn,依次令n等于1,2,…,N,其中N为天线结构节点总数,可依次计算得到相应的反射面节点轴向热变形列向量δ12,…,δN,将其组成影响矩阵Υ=[δ12,…,δN];
(3c)根据步骤(3a)和(3b)的插值矩阵和影响矩阵,建立映射矩阵Ω=Υ·Γ,将其存储为映射矩阵文本文件。
进一步,所述步骤(4)按如下过程进行:
(4a)根据有限元模型,提取反射面的节点集合{Node|m=1,2,…,M},其中m表示第m个反射面节点,M为反射面节点总数目,提取反射面的单元集合{Element|e=1,2,…,El},其中e表示第e个单元,El为单元总数目,提取反射面的单元中心点集合{Center|e=1,2,…,El};
(4b)依据有限元模型和天线辐射场积分公式建立理想矩阵;
(4c)依据有限元模型和天线辐射场积分公式建立扰动矩阵。
进一步,所述步骤(5)按如下过程进行:
(5a)将相位区间[-π,π]分成Q个子区间针对第q个子区间将相位指数在该区间中点位置进行一阶泰勒级数展开得到相位指数的展开形式依次令q等于1,2,…,Q,可依次得到所有子区间对应的展开系数;
(5b)依据步骤(3c)建立的温度到热变形的映射矩阵计算反射面节点热变形列向量δreflector=ΩTmeasure,其中Tmeasure为实时测量温度列向量,Ω为映射矩阵,依据热变形列向量δreflector和反射面单元连接关系,计算反射面单元中心点的热变形列向量其中 为第e个反射面单元中心点的热变形,上标“Tr”为转置符号,El为反射面单元总数目,分别为第e个反射面单元的三个顶点的热变形;
(5c)依据步骤(5b)得到的反射面三角形单元中心点热变形来计算其对应的相位误差ξcenter=kδcenter,进而根据相位误差ξcenter确定各个单元中心点相位误差所在的子区间及其相应的展开系数;
(5d)依据步骤(5c)所有反射面单元对应的展开系数,建立理想矩阵系数;
(5e)依据步骤(5b)得到的反射面三角形节点热变形列向量δreflector计算各个节点热变形对应的相位误差ξreflector=kδreflector,进而根据相位误差确定各个节点相位误差所在的子区间及其相应的展开系数;
(5f)依据步骤(5e)得到的所有反射面节点对应的展开系数,建立扰动矩阵系数。
进一步,所述步骤(6)按如下过程进行:
(6a)依据温度采集系统将实时测量的温度数据组成列向量Ttemperature
(6b)依据步骤(5b)快速计算反射面节点热变形列向量δreflector
(6c)依据步骤(5d)和(5f)快速建立理想矩阵系数R和扰动矩阵系数I;
(6d)从文本文件读取映射矩阵文本文件、理想矩阵文本文件、扰动矩阵文本文件,分别生成映射矩阵Ω、理想矩阵T和扰动矩阵H;
(6e)依据下述热电矩阵耦合模型快速评估反射面天线的远场电场;
(6f)依据远场电场的离散数据列向量E快速计算日照温度场作用下的反射面天线的增益损失、副瓣电平、指向误差指标,实现日照温度场作用下的反射面天线电性能快速评估。
进一步,所述步骤(7)按如下过程进行:
(7a)依据步骤(6f)评估的方位方向的指向误差Δθfw和俯仰方向的指向误差Δθfy,按如下公式计算天线方位角调整量Δφfw和俯仰角调整量Δφfy
(7b)将方位角调整量Δφfw和俯仰角调整量Δφfy实时发送到天线控制单元,实时修正反射面天线的方位角和俯仰角,实现日照温度场作用下反射面天线指向的快速补偿。
本发明与现有技术相比,具有以下优点:
1.本发明直接建立了测量温度与电性能之间的热电矩阵耦合模型,利用矩阵运算代替现有工程中一系列复杂的仿真过程,温度与电性能之间的函数关系简洁明了,方便工程中日照温度场作用下反射面天线电性能评估。
2.本发明将与温度无关的矩阵提前计算好,并存储为文本文件,工程中进行电性能评估时,直接调用文本文件即可,不必重新计算所有矩阵,计算量大幅度减小,使得电性能评估和补偿过程具有实时性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的不当限定,在附图中:
图1为本发明的实现流程图;
图2为映射矩阵文本文件生成流程图;
图3为反射面电性能快速评估流程图;
图4为7.3米反射面天线有限元模型;
图5为反射面温度测量云图;
图6为日照热环境下反射面指向误差时程曲线。
具体实施方式
下面将结合附图以及具体实施例来详细说明本发明,在此本发明的示意性实施例以及说明用来解释本发明,但并不作为对本发明的限定。
参照图1所示的本发明的实现流程图,本发明一种日照温度场作用下的反射面天线电性能快速评估与补偿方法,其实现步骤如下:
步骤一,建立反射面天线有限元模型。
依据反射面天线结构参数以及材料属性,通过结构分析软件ANSYS建立有限元模型,其中,反射面天线结构参数包括天线背架、中心体、主反射面和副反射面的几何参数信息,材料属性包括密度、泊松比、弹性模量和热膨胀系数;建立有限元模型且将反射面划分为三角形网格单元。
步骤二,确定温度传感器的布局位置。
依据实际天线中温度传感器布局确定其在有限元模型中的位置,即确定温度传感器在有限元模型中的节点坐标(us,vs,ws),其中s=1,2,…,Sensor,s表示第s个温度传感器,Sensor为温度传感器总数。
步骤三,建立温度到热变形的映射矩阵。
依据温度传感器在有限元模型中的位置建立温度到热变形的映射矩阵,参照图2所示的映射矩阵文本文件生成流程图,映射矩阵文本文件具体生成过程如下:
(3a)假设第s个温度传感器测量数据为1℃,其余温度传感器测量数据为0℃,基于空间距离平方反比插值方法计算有限元模型所有节点的插值温度列向量为Ts,依次令s等于1,2,…,Sensor,可依次计算得到相应的插值温度列向量将其组成插值矩阵
(3b)假设第n个有限元节点温度为1℃,其余有限元节点温度为0℃,根据有限元模型,设置参考温度为0℃,计算所有反射面节点轴向热变形列向量为δn,依次令n等于1,2,…,N,其中N为天线结构节点总数,可依次计算得到相应的反射面节点轴向热变形列向量δ12,…,δN,将其组成影响矩阵Υ=[δ12,…,δN];
(3c)根据步骤(3a)和(3b)的插值矩阵和影响矩阵,建立映射矩阵Ω=Υ·Γ,将其存储为映射矩阵文本文件。
步骤四,建立理想矩阵和扰动矩阵。
依据有限元模型和天线辐射场积分公式建立理想矩阵和扰动矩阵,理想矩阵和扰动矩阵文本文件生成过程如下:
(4a)根据有限元模型,提取反射面的节点集合{Node|m=1,2,…,M},其中m表示第m个反射面节点,M为反射面节点总数目,提取反射面的单元集合{Element|e=1,2,…,El},其中e表示第e个单元,El为单元总数目,提取反射面的单元中心点集合{Center|e=1,2,…,El};
(4b)依据有限元模型和天线辐射场积分公式建立理想矩阵
其中,元素其表示第e个反射面单元在第g个观测方向的辐射场,上标e表示第e个反射面单元,为全局坐标系下的球坐标分量,其表示第g个观测方向,G为远场观测点总数目,表示反射面节点向量r′位置的面电流,σe表示第e个单元的积分区域,j为虚数符号,k为波常数,为第g个观测方向的单位矢量,d为微分符号,将理想矩阵存储为理想矩阵文本文件;
(4c)依据有限元模型和天线辐射场积分公式建立扰动矩阵
其中,元素为扰动矩阵H的中间变量,表示第e个反射面单元对第g个观测方向辐射场造成的扰动量,其中i=1,2,3表示三角形单元的三个顶点序号,Ni为三角形单元的有限元插值函数,θs为面板节点向量r′与馈源坐标系z轴的夹角,表示反射面节点向量r′位置的面电流,σe表示第e个单元的积分区域,j为虚数符号,k为波常数,为第g个观测方向的单位矢量,d为微分符号,表示单元组集符号,即将 进行单元组集,将扰动矩阵存储为扰动矩阵文本文件。
步骤五,建立理想矩阵系数和扰动矩阵系数。
依据映射矩阵和温度采集数据建立扰动矩阵系数和理想矩阵系数,扰动矩阵系数和理想矩阵系数的具体生成过程如下:
(5a)将相位区间[-π,π]分成Q个子区间针对第q个子区间该区间的中点相位为将相位指数在该区间中点位置进行一阶泰勒级数展开得到其中ξq为第q个子区间内部的任意相位,aq、cq、bq和dq分别为该区间对应的展开系数,依次令q等于1,2,…,Q,可依次得到所有子区间对应的展开系数;
(5b)依据步骤(3c)建立的温度到热变形的映射矩阵计算反射面节点热变形列向量δreflector=ΩTmeasure,其中Tmeasure为实时测量温度列向量,Ω为映射矩阵,依据热变形列向量δreflector和反射面单元连接关系,计算反射面单元中心点的热变形列向量其中 为第e个反射面单元中心点的热变形,上标“Tr”为转置符号,El为反射面单元总数目,分别为第e个反射面单元的三个顶点的热变形;
(5c)依据步骤(5b)得到的反射面三角形单元中心点热变形来计算其对应的相位误差ξcenter=kδcenter,k为波常数;进而根据相位误差ξcenter确定各个单元中心点相位误差所在的子区间及其相应的展开系数,设第e个单元对应的展开系数为ae、ce、be和de
(5d)依据步骤(5c)所有反射面单元对应的展开系数,建立理想矩阵系数
(5e)依据步骤(5b)得到的反射面三角形节点热变形列向量δreflector计算各个节点热变形对应的相位误差ξreflector=kδreflector,进而根据相位误差确定各个节点相位误差所在的子区间及其相应的展开系数,设第m个节点对应的展开系数为am、cm、bm和dm
(5f)依据步骤(5e)得到的所有反射面节点对应的展开系数,建立扰动矩阵系数
步骤六,依据热电矩阵耦合模型和实测温度快速评估反射面天线的电性能。
参照图3所示的反射面电性能快速评估流程图,依据热电矩阵耦合模型快速评估反射面天线的电性能,具体实现过程如下:
(6a)依据温度采集系统将实时测量的温度数据组成列向量Ttemperature
(6b)依据步骤(5b)快速计算反射面节点热变形列向量δreflector
(6c)依据步骤(5d)和(5f)快速建立理想矩阵系数R和扰动矩阵系数I;
(6d)从文本文件读取映射矩阵文本文件、理想矩阵文本文件、扰动矩阵文本文件,分别生成映射矩阵Ω、理想矩阵T和扰动矩阵H;
(6e)依据下述热电矩阵耦合模型快速评估反射面天线的远场电场
E=Π[TR+HIΩTtemperature],
其中,η为波阻抗,r为远场观测点距离天线的距离,π为圆周率,
E为远场电场的离散数据列向量;
(6f)依据远场电场的离散数据列向量E快速计算日照温度场作用下的反射面天线的增益损失、副瓣电平、指向误差等指标,实现日照温度场作用下的反射面天线电性能快速评估。
步骤七,依据电性能评估结果调整方位角和俯仰角以改善天线指向精度。
依据电性能评估结果实时调整反射面天线的方位角和俯仰角,实现日照温度场作用下的反射面天线指向的快速补偿,具体实现过程如下:
(7a)依据步骤(6f)评估的方位方向的指向误差Δθfw和俯仰方向的指向误差Δθfy,按如下公式计算天线方位角调整量Δφfw和俯仰角调整量Δφfy
Δφfy=-Δθfy
其中,sign为符号函数,θfy为反射面天线当前的俯仰角;
(7b)将方位角调整量Δφfw和俯仰角调整量Δφfy实时发送到天线控制单元,
实时修正反射面天线的方位角和俯仰角,实现日照温度场作用下反射面天线指向的快速补偿。
本发明的效果可通过以下案例进一步说明:
1.实验对象
针对如图4所示的7.3米反射面天线,对本发明方法进行验证。面板材料为硬铝,其他结构材料为20号钢材,天线姿态为俯仰角90度(仰天),方位角90度(朝东),天线地理位置为西安地区(东经108度,北纬34度),面板安装21个温度传感器,背架安装48个温度传感器,温度测量日期为2012年11月23日。
2.计算结果
图5为反射面天线面板的温度测量云图,利用不同时刻温度测量数据,基于本发明的电性能快速计算方法,可快速计算不同时刻反射面天线的电性能。以反射面天线指向精度为例,图6给出反射面天线不同时刻温度场影响下E平面和H平面波束指向精度的计算结果。显然,随着太阳与天线之间相对转动,E平面和H平面波束指向精度也随之改变,基于指向精度的快速评估结果,可实现反射面天线指向精度的实时调整。
从上述仿真案例可以看出,应用本发明可以通过在天线结构上布置若干温度传感器,基于实时测量温度数据和热电矩阵耦合模型,快速评估反射面天线的电性能,基于评估结果,实现日照温度场作用下反射面天线电性能的实时补偿,使其满足电性能的指标要求,具有重要的工程应用价值。
本发明并不局限于上述实施例,在本发明公开的技术方案的基础上,本领域的技术人员根据所公开的技术内容,不需要创造性的劳动就可以对其中的一些技术特征作出一些替换和变形,这些替换和变形均在本发明的保护范围内。

Claims (8)

1.一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,包括如下步骤:
(1)依据反射面天线结构参数以及材料属性建立有限元模型;
(2)依据实际天线中温度传感器布局确定其在有限元模型中的位置;
(3)依据温度传感器在有限元模型中的位置建立温度到热变形的映射矩阵;
(4)依据有限元模型和天线辐射场积分公式建立扰动矩阵和理想矩阵;
(5)依据映射矩阵和温度采集数据建立扰动矩阵系数和理想矩阵系数;
(6)依据实测温度和热电矩阵耦合模型快速评估日照温度场作用下的反射面天线的电性能;
(7)依据电性能评估结果实时调整反射面天线的方位角和俯仰角,实现日照温度场作用下的反射面天线指向的快速补偿。
2.根据权利要求1所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(1)中,反射面天线结构参数包括天线背架、中心体、主反射面和副反射面;所述材料属性包括密度、泊松比、弹性模量和热膨胀系数;建立有限元模型且将反射面划分为三角形网格单元。
3.根据权利要求1所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(2)中,依据实际天线中温度传感器布局确定其在有限元模型中的位置,即确定温度传感器在有限元模型中的节点坐标(us,vs,ws),其中s=1,2,…,Sensor,s表示第s个温度传感器,Sensor为温度传感器总数。
4.根据权利要求1所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(3)中,依据温度传感器在有限元模型中的位置建立温度到热变形的映射矩阵,按如下过程进行:
(3a)假设第s个温度传感器测量数据为1℃,其余温度传感器测量数据为0℃,基于空间距离平方反比插值方法计算有限元模型所有节点的插值温度列向量为Ts,依次令s等于1,2,…,Sensor,Sensor为温度传感器总数,可依次计算得到相应的插值温度列向量将其组成插值矩阵
(3b)假设第n个有限元节点温度为1℃,其余有限元节点温度为0℃,根据有限元模型,设置参考温度为0℃,计算所有反射面节点轴向热变形列向量为δn,依次令n等于1,2,…,N,其中N为天线结构节点总数,可依次计算得到相应的反射面节点轴向热变形列向量δ12,…,δN,将其组成影响矩阵Υ=[δ12,…,δN];
(3c)根据步骤(3a)和(3b)的插值矩阵和影响矩阵,建立映射矩阵Ω=Υ·Γ,将其存储为映射矩阵文本文件。
5.根据权利要求1所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(4)中,依据有限元模型和天线辐射场积分公式建立扰动矩阵和理想矩阵,按如下过程进行:
(4a)根据有限元模型,提取反射面的节点集合{Node|m=1,2,…,M},其中m表示第m个反射面节点,M为反射面节点总数目,提取反射面的单元集合{Element|e=1,2,…,El},其中e表示第e个单元,El为单元总数目,提取反射面的单元中心点集合{Center|e=1,2,…,El};
(4b)依据有限元模型和天线辐射场积分公式建立理想矩阵:
其中,元素其表示第e个反射面单元在第g个观测方向的辐射场,上标e表示第e个反射面单元,为全局坐标系下的球坐标分量,其表示第g个观测方向,G为远场观测点总数目,表示反射面节点向量r′位置的面电流,σe表示第e个单元的积分区域,j为虚数符号,k为波常数,为第g个观测方向的单位矢量,d为微分符号,将理想矩阵存储为理想矩阵文本文件;
(4c)依据有限元模型和天线辐射场积分公式建立扰动矩阵:
其中,元素为扰动矩阵H的中间变量,表示第e个反射面单元对第g个观测方向辐射场造成的扰动量,其中i=1,2,3表示三角形单元的三个顶点序号,Ni为三角形单元的有限元插值函数,θs为面板节点向量r′与馈源坐标系z轴的夹角,表示反射面节点向量r′位置的面电流,σe表示第e个单元的积分区域,j为虚数符号,k为波常数,为第g个观测方向的单位矢量;表示单元组集符号,即将 进行单元组集,将扰动矩阵存储为扰动矩阵文本文件。
6.根据权利要求1所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(5)中,依据映射矩阵和温度采集数据建立扰动矩阵系数和理想矩阵系数,按如下过程进行:
(5a)将相位区间[-π,π]分成Q个子区间针对第q个子区间该区间的中点相位为将相位指数在该区间中点位置进行一阶泰勒级数展开得到其中ξq为第q个子区间内部的任意相位,j为虚数符号,aq、cq、bq和dq分别为该区间对应的展开系数,依次令q等于1,2,…,Q,可依次得到所有子区间对应的展开系数;
(5b)依据步骤(3c)建立的温度到热变形的映射矩阵计算反射面节点热变形列向量δreflector=ΩTmeasure,其中Tmeasure为实时测量温度列向量,Ω为映射矩阵,依据热变形列向量δreflector和反射面单元连接关系,计算反射面单元中心点的热变形列向量其中 为第e个反射面单元中心点的热变形,上标“Tr”为转置符号,El为反射面单元总数目,分别为第e个反射面单元的三个顶点的热变形;
(5c)依据步骤(5b)得到的反射面三角形单元中心点热变形来计算其对应的相位误差ξcenter=kδcenter,k为波常数;进而根据相位误差ξcenter确定各个单元中心点相位误差所在的子区间及其相应的展开系数,设第e个单元中心点对应的展开系数为ae、ce、be和de
(5d)依据步骤(5c)所有反射面单元对应的展开系数,建立理想矩阵系数:
式中,El为反射面单元总数目;
(5e)依据步骤(5b)得到的反射面三角形节点热变形列向量δreflector计算各个节点热变形对应的相位误差ξreflector=kδreflector,进而根据相位误差确定各个节点相位误差所在的子区间及其相应的展开系数,设第m个节点对应的展开系数为am、cm、bm和dm
(5f)依据步骤(5e)得到的所有反射面节点对应的展开系数,建立扰动矩阵系数
式中,M为反射面节点总数目。
7.根据权利要求1所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(6)中依据实测温度和热电矩阵耦合模型快速评估日照温度场作用下的反射面天线的电性能,按如下过程进行:
(6a)依据温度采集系统将实时测量的温度数据组成列向量Ttemperature
(6b)依据步骤(5b)快速计算反射面节点热变形列向量δreflector
(6c)依据步骤(5d)和(5f)快速建立理想矩阵系数R和扰动矩阵系数I;
(6d)从文本文件读取映射矩阵文本文件、理想矩阵文本文件、扰动矩阵文本文件,分别生成映射矩阵Ω、理想矩阵T和扰动矩阵H;
(6e)依据下述热电矩阵耦合模型快速评估反射面天线的远场电场
E=Π[TR+HIΩTtemperature],
其中,j为虚数符号,k为波常数,η为波阻抗,r为远场观测点距离天线的距离,E为远场电场的离散数据列向量;
(6f)依据远场电场的离散数据列向量E快速计算日照温度场作用下的反射面天线的增益损失、副瓣电平和指向误差指标,实现日照温度场作用下的反射面天线电性能快速评估。
8.根据权利要求7所述的一种日照温度场作用下反射面天线电性能快速评估与补偿方法,其特征在于,所述步骤(7)按如下过程进行:
(7a)依据步骤(6f)评估的方位方向的指向误差Δθfw和俯仰方向的指向误差Δθfy,按如下公式计算天线方位角调整量Δφfw和俯仰角调整量Δφfy
Δφfy=-Δθfy
其中,sign为符号函数,θfy为反射面天线当前的俯仰角;
(7b)将方位角调整量Δφfw和俯仰角调整量Δφfy实时发送到天线控制单元,实时修正反射面天线的方位角和俯仰角,实现日照温度场作用下反射面天线指向的快速补偿。
CN201910759994.8A 2019-08-16 2019-08-16 日照温度场作用下反射面天线电性能快速评估与补偿方法 Pending CN110470916A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910759994.8A CN110470916A (zh) 2019-08-16 2019-08-16 日照温度场作用下反射面天线电性能快速评估与补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910759994.8A CN110470916A (zh) 2019-08-16 2019-08-16 日照温度场作用下反射面天线电性能快速评估与补偿方法

Publications (1)

Publication Number Publication Date
CN110470916A true CN110470916A (zh) 2019-11-19

Family

ID=68511899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910759994.8A Pending CN110470916A (zh) 2019-08-16 2019-08-16 日照温度场作用下反射面天线电性能快速评估与补偿方法

Country Status (1)

Country Link
CN (1) CN110470916A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111274726A (zh) * 2020-01-16 2020-06-12 东南大学 一种考虑热效应的天线罩电磁性能分析方法
CN112347669A (zh) * 2020-10-09 2021-02-09 中国科学院国家天文台 一种大型天线背架温度测量与实时评估系统及方法
CN113297766A (zh) * 2021-05-26 2021-08-24 中国电子科技集团公司第五十四研究所 一种天线单元加工误差补偿方法
CN113625062A (zh) * 2021-07-29 2021-11-09 西安电子科技大学 基于泰勒展开法的天线罩电性能估计方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7139661B2 (en) * 2000-07-15 2006-11-21 Robert Bosch Gmbh Information and control system for vehicles
CN102243276A (zh) * 2011-03-28 2011-11-16 西安电子科技大学 大型地基面天线的温度载荷机电耦合分析方法
CN103488817A (zh) * 2013-09-02 2014-01-01 西安电子科技大学 一种变形大型单反射面天线的馈源位置补偿方法
CN104038295A (zh) * 2014-06-06 2014-09-10 西安电子科技大学 一种基于机电耦合的变形阵列天线散射性能分析方法
US20140334518A1 (en) * 2013-05-08 2014-11-13 New York University System, method and computer-accessible medium for determining specific absorption rate obtained based on magnetic resonance imaging and temperature property measurements
CN105760600A (zh) * 2016-02-17 2016-07-13 西安电子科技大学 基于机电耦合的星载有源相控阵天线组件热功耗确定方法
CN105787160A (zh) * 2016-02-17 2016-07-20 西安电子科技大学 基于机电耦合的星载有源相控阵天线结构轻量化设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7139661B2 (en) * 2000-07-15 2006-11-21 Robert Bosch Gmbh Information and control system for vehicles
CN102243276A (zh) * 2011-03-28 2011-11-16 西安电子科技大学 大型地基面天线的温度载荷机电耦合分析方法
US20140334518A1 (en) * 2013-05-08 2014-11-13 New York University System, method and computer-accessible medium for determining specific absorption rate obtained based on magnetic resonance imaging and temperature property measurements
CN103488817A (zh) * 2013-09-02 2014-01-01 西安电子科技大学 一种变形大型单反射面天线的馈源位置补偿方法
CN104038295A (zh) * 2014-06-06 2014-09-10 西安电子科技大学 一种基于机电耦合的变形阵列天线散射性能分析方法
CN105760600A (zh) * 2016-02-17 2016-07-13 西安电子科技大学 基于机电耦合的星载有源相控阵天线组件热功耗确定方法
CN105787160A (zh) * 2016-02-17 2016-07-20 西安电子科技大学 基于机电耦合的星载有源相控阵天线结构轻量化设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
连培园: "大型微波反射面天线机电耦合若干问题研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111274726A (zh) * 2020-01-16 2020-06-12 东南大学 一种考虑热效应的天线罩电磁性能分析方法
CN111274726B (zh) * 2020-01-16 2021-03-16 东南大学 一种考虑热效应的天线罩电磁性能分析方法
CN112347669A (zh) * 2020-10-09 2021-02-09 中国科学院国家天文台 一种大型天线背架温度测量与实时评估系统及方法
CN112347669B (zh) * 2020-10-09 2024-04-12 中国科学院国家天文台 一种大型天线背架温度测量与实时评估系统及方法
CN113297766A (zh) * 2021-05-26 2021-08-24 中国电子科技集团公司第五十四研究所 一种天线单元加工误差补偿方法
CN113625062A (zh) * 2021-07-29 2021-11-09 西安电子科技大学 基于泰勒展开法的天线罩电性能估计方法
CN113625062B (zh) * 2021-07-29 2022-05-17 西安电子科技大学 基于泰勒展开法的天线罩电性能估计方法

Similar Documents

Publication Publication Date Title
CN110470916A (zh) 日照温度场作用下反射面天线电性能快速评估与补偿方法
CN106405533B (zh) 基于约束加权最小二乘的雷达目标联合同步与定位方法
CN110132261A (zh) 一种基于数值拟合的高精度星上轨道预报方法
CN103675927B (zh) 固定翼航空电磁系统接收吊舱摆动角度的校正方法
CN109633526A (zh) 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法
Desroziers A coordinate change for data assimilation in spherical geometry of frontal structures
CN110220491A (zh) 一种无人机的光学吊舱安装误差角估算方法
CN104197839A (zh) 航天器装配精度受重力和温度影响的补偿方法
CN110081905A (zh) 一种基于单站光电经纬仪的光波大气折射误差计算方法
CN110244367B (zh) 一种基于地面多基站的ztem系统姿态补偿方法
CN108090306A (zh) 一种基于副面结构应变的变形天线副面形貌快速重构方法
CN104458123A (zh) 利用零部件的质量和质心信息计算卫星质心的方法
CN108153954B (zh) 一种基于结构热类比的反射面天线温度快速重构方法
CN110456304A (zh) 机载测向定位方法
Patra et al. Vertical ExB drifts from radar and C/NOFS observations in the Indian and Indonesian sectors: Consistency of observations and model
Guo et al. Evaluation of terrain effect on microwave radiometer measurement and its correction
Shen et al. Nonlinear magnetic gradients and complete magnetic geometry from multispacecraft measurements
Lian et al. Real‐time temperature estimation method for electromagnetic performance improvement of a large axisymmetric radio telescope under solar radiation
CN108470096A (zh) 一种基于卡尔曼滤波的有源相控阵天线形面重构方法
CN105868780A (zh) 基于指向调整的天线热变形补偿数据库匹配调用方法
CN103674000B (zh) 电子罗盘实时校准算法
Li et al. An improved velocity volume processing method
CN111025358B (zh) 基于导航卫星信号单短基线的定向方法
CN108828585B (zh) 一种空间平台目标无伺服跟踪测角方法
Lu et al. Motion model of floating weather sensing node for typhoon detection

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20191119

WD01 Invention patent application deemed withdrawn after publication