CN114970254B - 基于选取关键索准则的索网天线型面精度调整方法 - Google Patents

基于选取关键索准则的索网天线型面精度调整方法 Download PDF

Info

Publication number
CN114970254B
CN114970254B CN202210528465.9A CN202210528465A CN114970254B CN 114970254 B CN114970254 B CN 114970254B CN 202210528465 A CN202210528465 A CN 202210528465A CN 114970254 B CN114970254 B CN 114970254B
Authority
CN
China
Prior art keywords
cable
antenna
adjustment
network antenna
profile
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.)
Active
Application number
CN202210528465.9A
Other languages
English (en)
Other versions
CN114970254A (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 University
Original Assignee
Xidian University
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 University filed Critical Xidian University
Priority to CN202210528465.9A priority Critical patent/CN114970254B/zh
Publication of CN114970254A publication Critical patent/CN114970254A/zh
Application granted granted Critical
Publication of CN114970254B publication Critical patent/CN114970254B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明公开了一种基于选取关键索准则的索网天线型面精度调整方法,不同于采用索网调整方法使数学参考模型节点逼近实物模型的方法,本发明提出的数学模型修正方法基于索网天线平衡矩阵并测量部分独立索单元索力实现数学模型修正。不同于采用敏度矩阵方法计算的调整量直接调整索网天线,本发明提出的基于关键索选取准则的索网天线调整方法基于修正模型的敏度矩阵,引入质量份数衡量可调整索单元调整量对索网天线型面精度的影响因子,将质量份数作为关键索选取准则挑选部分可调整索单元作为关键索,并计算调整量对索网天线进行调整。

Description

基于选取关键索准则的索网天线型面精度调整方法
技术领域
本发明属于星载天线反射面调整技术领域,涉及一种基于选取关键索准则的索网天线型面精度调整方法。
背景技术
索网天线具有口径大、精度高、轻量化、收纳比高等优点,被广泛应用于通信、中继和侦察领域。
由于索网天线在设计阶段以及制造和装配的过程中不可避免地引入随机误差,将导致天线的反射面精度不能满足设计要求,影响天线的工作性能。为了尽量降低或者消除这些误差,提高索网天线型面精度,需要对其进行必要的索网调整工作。但随着天线口径不断增大,索网结构节点数量增多,调整索数目激增给天线型面精度调整效率带来很大难题。
索网天线是一种非线性结构,节点相互耦合,即其中一个节点偏移导致周围其他节点的移动。若竖向索的调整没有关联天线结构的耦合,导致即使对于每一根调整索精确计算调整长度并准确完成调整,在下一次依旧需要对该调整索进行调整。
为提高索网反射面调整效率,减少天线调整的工作量,需要研究新的高效的型面精度调整方法,使得仅用有限次数的测量和调整少量的索,就能使天线反射面的型面精度达到期望的要求。
发明内容
本发明的目的是提供一种基于选取关键索准则的索网天线型面精度调整方法,采用该方法能够提高索网天线型面精度调整效率。
本发明所采用的技术方案是,基于选取关键索准则的索网天线型面精度调整方法,具体包括如下步骤:
步骤1,根据设计实物天线时的理想天线建立有限元模型作为初始数学模型,测量实物索网天线所有节点的位置
Figure BDA0003645560830000021
组合所有节点平衡方程,得到索网天线自由节点的平衡矩阵B3(n-c)×m
步骤2,计算由步骤1得到的平衡矩阵B的秩RB,得到独立自由模态数s1=m-RB,计算挑选出s1根索单元编号作为待测量索力索单元编号序列N1,并将实物天线索力T分为待测量索力T1和待求解索力T2,待测量索力T1和待求解索力T2相对应的平衡矩阵分块为B1、B2
步骤3,测量实物索网天线索单元编号序列为N1的索力T1,通过分块平衡方程计算待求解索力T2
Figure BDA0003645560830000022
步骤4,令索网天线的调整次数k=0,记实物天线初始索力T(0)=[T1,T2];根据步骤1得到的实物索网天线位置
Figure BDA0003645560830000023
和实物天线索力T(0)修改索网天线数学模型中相应的节点位置为/>
Figure BDA0003645560830000024
修改数学模型的索力为T(0),完成索网天线数学模型参数的修正;
步骤5,令调整次数k=k+1,在步骤4修正的索网天线数学模型型面顶点建立直角坐标系Oxyz,根据数学模型全部节点坐标
Figure BDA0003645560830000025
得到数学模型的型面n1个标记点的坐标
Figure BDA0003645560830000026
计算n1个标记点的期望位移/>
Figure BDA0003645560830000031
计算索网天线型面上的标记点位置关于调整索索长的敏度矩阵/>
Figure BDA0003645560830000032
步骤6,挑选关键索的质量份数作为关键索选取准则,则第i根可调整索的质量份数Pi为:
Figure BDA0003645560830000033
式中,
Figure BDA0003645560830000034
为几何逼近误差δp的平方,Δz为型面标记点的期望z向位移,/>
Figure BDA0003645560830000035
为第i根调整索对应的敏度矩阵S(k)中第i列;
步骤7,将从全部的m1根可调整索质量份数进行排序,按照质量份数由高到低地挑选出s2根可调整索作为关键索,计算第k次调整关键索的调整量
Figure BDA0003645560830000036
Figure BDA0003645560830000037
步骤8,将步骤7所得的关键索的调整量
Figure BDA0003645560830000038
施加到实物模型上,测量得到调整后的实物反射面节点坐标/>
Figure BDA0003645560830000039
步骤9,对实物天线关键索的调整量进行修正,修正后的关键索调整量
Figure BDA00036455608300000310
为:
Figure BDA00036455608300000311
步骤10,将步骤9修正后的关键索调整量
Figure BDA00036455608300000312
施加到数学模型上,得到修正后的数学模型所有节点坐标/>
Figure BDA00036455608300000313
和标记点坐标/>
Figure BDA00036455608300000314
计算第k次调整之后的天线型面精度δ(k),至此完成第k次调整;
步骤11,重复步骤5~10,直到索网天线型面精度δ(k)≤[δ],完成索网天线关键索调整。
本发明的特点还在于:
步骤1中,根据设计实物天线的理想天线建立初始有限元数学模型,实物天线由n个节点、m根索单元组成,其中c个固定节点,索网天线型面的标记节点数为n1,可调整索单元m1根,初始型面误差为δ0,型面误差的调整目标为[δ],实物天线的平衡矩阵B3(n-c)×m为组合所有自由节点平衡方程中索力的系数矩阵,根据平衡索网天线中的任意自由节点应满足x、y、z方向的力平衡方程,其中自由节点i的x方向为:
Figure BDA0003645560830000041
式中,Ci为与自由节点i相连的节点编号的集合,索单元ij连接节点i、j,Tij为索单元ij的张力,Lij为索单元ij的长度,xi和xj分别为节点i、j的x方向坐标;
则节点i包含3个方向的力平衡方程:
Figure BDA0003645560830000042
组合所有自由节点平衡方程,记索网天线自由节点平衡方程中索力的系数矩阵为平衡矩阵B,则所有自由节点的平衡方程为:
B3(n-c)×mTm×1=03(n-c)×1 (7);
式中,B3(n-c)×m为索网天线由节点的索力平衡矩阵,Tm×1为索单元张力列向量,n为索网天线节点总数,c为索网天线固定节点数,(n-c)为索网天线自由节点总数,m为索网天线索单元总数。
步骤2中,由于平衡矩阵的秩RB的表达式如下:
RB=3(n-c)<m (8);
基于公式(8)的不等式,所以公式(7)存在不唯一解,为了对索网天线进行唯一修正,需要挑选独立不相关的索单元测量索力,然后对天线索力进行求解修正,计算平衡矩阵的秩RB得独立自应力模态数s1=m-RB,从全部的m根索单元中随机给定s1根索单元编号作为初始待测量索力编号序列N1,其余索力编号作为待求解索力编号序列N2,并依此将索力分为待测量索力T1和待求解索力T2,相应的平衡矩阵分块为B1、B2,则平衡矩阵表示为:
Figure BDA0003645560830000052
计算分块矩阵B2的秩,若B2不满秩则将编号序列N1最后一个索单元编号放到N2最前面,将N2最后一个编号放到N1最前面,然后再计算B2的秩,直到B2满秩,则完成挑选独立不相关待测量索力编号序列N1
步骤3中,实测实物模型的索力T1,则将T1代入公式(9)中,根据公式(1)得到待求解索力T2
步骤5中,索网型面上存在n1个标记节点,为确定第k次调整时的敏度矩阵S(k)和索网标记点期望的位移
Figure BDA0003645560830000051
需要计算索网天线的调整目标抛物面,步骤5的具体步骤为:
步骤5.1,在数学模型型面顶点建立直角坐标系Oxyz,其中,z轴沿抛物面轴向,拟合的最佳几何逼近抛物面为:
Figure BDA0003645560830000061
/>
其中,x、y、z为索网天线抛物面上点的坐标,f为抛物面焦距,Δh1为最佳几何逼近抛物面相对于坐标原点的z向偏移量,其中几何逼近误差由数值积分计算:
Figure BDA0003645560830000062
式中,s为上网面有效口径内划分三角形单元总数,
Figure BDA0003645560830000063
为第k1个三角形面片在口径面内的投影区域,/>
Figure BDA0003645560830000064
是该面片在天线孔径方向投影面积,(Xi,Yj)为三角形面上的数值积分高斯点,F(Xi,Yj)为型面高斯点的z向偏差,W(i),W(j)为X、Y方向上的权系数/>
Figure BDA0003645560830000065
步骤5.2,求解型面标记点的调整期望位移
Figure BDA0003645560830000066
需要计算标记点调整目标抛物面,即当型面标记点都调整到目标抛物面,索网型面几何逼近误差达到最小,则调整目标抛物面方程为:
Figure BDA0003645560830000067
其中,Δh2为天线在设计时的几何逼近抛物面相对于设计的理想抛物面的z向偏移量,Δh1和Δh2通过一维搜索进退法求解。
则将每个型面标记点坐标代入式(12),得到型面标记点投影到调整目标抛物面的投影坐标
Figure BDA0003645560830000068
由此可得型面标记点的期望位移为:
Figure BDA0003645560830000071
步骤5.3,根据第k次调整时的数学模型求解敏度矩阵,从第i=1根可调整索开始,改变索网天线第i根可调整索的索力ΔT,分析数学模型的平衡态得到型面节点的位移向量ΔXi和该可调整索的索长变化量ΔLi,敏度矩矩阵的第i列为
Figure BDA0003645560830000072
将该可调整索的索力恢复改变之前的状态,对下一根可调整索进行分析,直到改变了所有的m1根可调整索,组装敏度矩阵为:
Figure BDA0003645560830000073
步骤6中,质量份数Pi反映了第i根可调整索对于索网天线型面精度的调整能力,作为关键调整索的选取准则,质量份数Pi的值定义为几何逼近方差
Figure BDA0003645560830000074
的对索调整量ΔLi的偏导数的绝对值。
步骤7中,
Figure BDA0003645560830000075
为从所有的m1根可调整索中挑选s2根关键索对应的敏度矩阵中列的组合,/>
Figure BDA0003645560830000076
为/>
Figure BDA0003645560830000077
的广义逆,/>
Figure BDA0003645560830000078
为实物天线选取关键索的调整量,其它可调整索调整量为0。
本发明的有益效果如下:
1.考虑了数学模型和实际模型响应之间不可避免的存在一定的偏差,提出对索网天线的数学模型进行参数型修正,通过测量天线节点坐标和测量少量索力得到实物天线的节点位置和全部索力,依此建立修正参数的有限元数学模型。使数学模型的有限元分析参数和实物模型的测量值趋于一致,从而让数学模型能反映实物模型的真实状态。并且在数学模型指导实物模型调整过程中,不断对数学模型进行修正和调整,确保调整过程的收敛性,使得计算的调整量更准确,索网天线很快便能调整到位。
2.为唯一地修正索网天线数学模型,组合索网天线自由节点的平衡方程,挑选少量的待测量索力对实物索网天线进行测量,然后求解平衡方程得到其余的索力,完成索网天线数学模型的修正。这种测量必要参数修正的数学模型有明确的物理意义,可操作性高。
3.考虑了随着索网天线口径越大,调整索数目越来越多,型面调整工作量大且非常耗时。通过引入关键调整索的选取准则,采用选取部分关键调整索进行索网天线型面精度调整,减少了每次调整的调整索数量,并通过敏度矩阵计算的关键索的调整量,可显著减少索网天线调整时间。
附图说明
图1是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网桁架结构示意图;
图2是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线数学模型结构示意图;
图3是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线反射面示意图;
图4是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线最佳几何逼近抛物面和调整目标抛物面说明图;
图5是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中进退法挑选独立待测量索时待求解索对应分块矩阵秩的变化图;
图6是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线修正数学模型挑选的待测量索单元示意图;
图7是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线多种调整方式下计算的调整量对照图;
图8是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线多种调整方式下反射面精度调整过程图;
图9是本发明基于选取关键索准则的索网天线型面精度调整方法实施例中索网天线多种调整方式调整天线耗费时长说明图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明提出的基于关键索选取准则的索网天线型面精度调整方法。不同于采用索网调整方法使数学参考模型节点逼近实物模型的方法,本发明提出的数学模型修正方法基于索网天线平衡矩阵并测量部分独立索单元索力实现数学模型修正。不同于采用敏度矩阵方法计算的调整量直接调整索网天线,本发明提出的基于关键索选取准则的索网天线调整方法基于修正模型的敏度矩阵,引入质量份数衡量可调整索单元调整量对索网天线型面精度的影响因子,将质量份数作为关键索选取准则挑选部分可调整索单元作为关键索,并计算调整量对索网天线进行调整。
本发明基于选取关键索准则的索网天线型面精度调整方法,具体包括如下步骤:
步骤1,根据设计实物索网天线时的理想天线,建立初始的索网天线有限元数学模型,其中实物索网天线由n个节点、m根索单元组成,其中c个固定节点,索网天线型面的标记节点数为n1,可调整索单元m1根,索网天线调整之前的初始型面误差为δ0,型面误差的调整目标为[δ]。实物索网天线应满足节点平衡方程,测量实物索网天线所有节点坐标
Figure BDA0003645560830000101
建立总体索网天线自由节点平衡方程,索网天线的自由节点应满足的索力平衡方程为:
Figure BDA0003645560830000102
式中,Ci为与自由节点i相连的节点编号的集合,索单元ij连接节点i、j,Tij为索单元ij的张力,Lij为索单元ij的长度,xi和xj分别为节点i、j的x方向坐标;
根据索网天线节点的拓扑连接关系,组装总体索网自由节点的平衡矩阵B,则索网天线总体自由节点平衡方程应为:
B3(n-c)×mTm×1=03(n-c)×1 (2);
式中,B3(n-c)×m为索网天线由节点的索力平衡矩阵,Tm×1为索单元张力列向量,(n-c)为索网系统自由节点总数。
步骤2,根据索网天线平衡矩阵挑选测量部分索力,然后得到总体索力,修正索网天线数学模型,过程如下:
步骤2.1,计算平衡矩阵B的秩RB,可得独立自由模态数s1=m-RB,从所有实物天线的m根索单元中挑选出s1根独立索单元的初始编号序列N1,其余索力编号序列为N2,并按照这些编号将实物天线的索力T(0)分为待测量索力T1和待求解索力T2,将平衡矩阵相应地分块为B1、B2。计算分块矩阵B2的秩,若分块矩阵B2不满秩则将N2得最后一个编号放到N1最前面,将N1的最后一个编号放到N2最前面,再计算B2的秩,若仍不满秩则再次交换序号,直到B2满秩。
步骤2.2,根据平衡方程的分块形式:
Figure BDA0003645560830000119
测量实物索网天线索单元编号序列为N1的索力T1,代入式(3),计算待求解索力T2为:
Figure BDA0003645560830000111
令索网天线的调整次数k=0,记实物天线初始索力T(0)=[T1,T2]。根据上述的实物索网天线位置
Figure BDA0003645560830000112
和实物天线索力T(0)修改索网天线数学模型中相应的节点位置/>
Figure BDA0003645560830000113
修改数学模型的索力为T(0),完成索网天线数学模型参数的修正。
步骤3,令调整次数k=k+1,在数学模型型面顶点建立直角坐标系Oxyz,根据数学模型所有节点坐标
Figure BDA0003645560830000114
可得到调整型面上标记节点的坐标/>
Figure BDA0003645560830000115
计算索网天线标记点的期望调整位移/>
Figure BDA0003645560830000116
计算索网天线型面标记点关于可调整索索长的敏度矩阵/>
Figure BDA0003645560830000117
步骤3.1,采用进退法调整Δh1直到几何逼近误差δP最小,拟合计算得到最佳几何逼近抛物面方程为:
Figure BDA0003645560830000118
其中x、y、z为索网天线抛物面上点的坐标,f为抛物面焦距,Δh1为最佳几何逼近抛物面相对于坐标原点的z向偏移量。
其中索网天线型面几何逼近误差由数值积分计算:
Figure BDA0003645560830000121
式中,s为上网面有效口径内划分三角形单元总数,
Figure BDA0003645560830000122
为第k1个三角形面片在口径面内的投影区域,/>
Figure BDA0003645560830000123
是该面片在天线孔径方向投影面积,(Xi,Yj)为三角形面上的数值积分高斯点坐标,F(Xi,Yj)为型面高斯点的z向位移,W(i)、W(j)为X、Y方向上的权系数,其值均为/>
Figure BDA0003645560830000124
步骤3.2,则调整目标抛物面方程为:
Figure BDA0003645560830000125
其中,Δh2为天线在设计时的几何逼近抛物面相对于设计的理想抛物面z向偏移量。
将每个型面标记点的坐标代入式(7),得到型面标记点投影到调整目标抛物面的节点坐标
Figure BDA0003645560830000126
由此可得型面节点的期望节点位移:
Figure BDA0003645560830000127
步骤3.3,根据第k次调整时的数学模型求解敏度矩阵,从第i=1根可调整索开始,稍微改变索网天线第i根可调整索的索力ΔT,分析数学模型的平衡态得到型面节点的位移向量ΔXi和该可调整索的索长变化量ΔLi,于是敏度矩矩阵的第i列为
Figure BDA0003645560830000131
然后将该可调整索的索力恢复改变之前的状态,对下一根可调整索进行分析,直到改变了所有的m1根可调整索,组装敏度矩阵为:
Figure BDA0003645560830000132
步骤4,计算所有m1根可调整索的质量份数,挑选需要调整的关键索,并计算关键索的调整量
Figure BDA0003645560830000133
步骤4.1,定义第i个质量份数Pi为几何逼近误差δp对于第i根可调整索调整量ΔLi的偏导数的绝对值:
Figure BDA0003645560830000134
式中,
Figure BDA0003645560830000135
为几何逼近误差的平方,Δz为型面节点的z向期望位移,/>
Figure BDA0003645560830000136
为第i根调整索对应的敏度矩阵S(k)第i列。
步骤4.2,对全部m1根可调整索的质量份数进行排序,从所有的m1根可调整索中按照质量份数由高到低地挑选s2根作为关键索并计算调整量:
Figure BDA0003645560830000137
步骤5,将关键索调整量
Figure BDA0003645560830000138
施加到实物索网天线上,测量调整后的实物索网天线型面标记点位置/>
Figure BDA0003645560830000139
根据实物天线型面标记点的位移修正索网天线实际调整量:
Figure BDA00036455608300001310
步骤6,将修正后的索网天线实际调整量
Figure BDA00036455608300001311
施加到数学模型上,得到第k次调整之后数学模型的所有节点位置/>
Figure BDA0003645560830000141
和型面节点坐标/>
Figure BDA0003645560830000142
并计算索网天线的几何逼近误差/>
Figure BDA0003645560830000143
完成第k次调整。
步骤7,重复步骤3~6,直到索网天线的几何逼近误差达到调整目标值
Figure BDA0003645560830000144
完成索网天线型面的关键索调整。
下面结合仿真实例进一步说明本发明的调整效果:
采用AstroMesh索网天线作为调整实例进行仿真调整,图1是AstroMesh索网天线的索网桁架整体构型,图2为实例索网天线数学模型结构示意图。其中索网天线主要的参数为:口径5m,前后网面焦径比为0.8,天线高度为1m,分环数为7,自由节点数为242,固定节点数为36,竖向索单元总数为121,前后网面索单元总数为756,型面标记点为前网面自由节点。索单元为圆形横截面的芳纶绳,其直径为0.8mm,弹性模量为20GPa,密度为1685kg/m3,泊松比为0.3。图3为索网天线反射型面,其设计的理想几何逼近误差为0.64369mm。
索网天线调整过程中计算的最佳几何逼近抛物面和型面节点的关系如图4所示,型面上的点分布在最佳几何逼近抛物面两侧,调整目标理想抛物面位于最佳几何逼近抛物面下方。挑选独立待测量索单元的迭代曲线如图5所示,经过多次挑选待计算索力的索单元对应的分块矩阵最终达到满秩,挑选的独立待测量索单元在索网天线中的分布如图6所示。
如图7所示,为验证本发明的调整效果,对比分析了三种不同调整量计算方式的调整方法:基于修正模型敏度矩阵的调整方法,基于修正模型的关键索调整方法,基于修正模型敏度计算的最大调整量的调整方法:其中敏度矩阵调整方法基于修正模型的敏度矩阵,计算所有可调整索的调整量,对实物天线模型进行型面调整;最大调整量索的调整方法基于修正模型敏度矩阵,计算所有可调整索的调整量,从可调整索中选取调整量较大的百分之三十对索网天线进行型面调整;关键索调整方法基于修正模型的敏度矩阵计算质量份数,选取质量份数较大的百分之三十的可调整索作为关键索对索网天线进行型面调整。
由图8所示,索网天线型面初始几何逼近误差为3.49mm,采用关键索调整方法经过3次调整,型面几何逼近误差便达到调整目标值0.94mm;采用敏度矩阵方法调整全部可调整索,经过2次调整型面误差下降至调整目标值;采用最大调整量索的调整方法调整型面,调整5次后型面误差下降至目标值。
由图9所示,索网天线调整过程中每个步骤所花费时长分别为:单次测量索网天线节点坐标耗时15分钟;单次数学模型修正包含节点测量时长和修正数学模型的数值计算时长,耗时17.5分钟;单次计算调整量需要计算敏度矩阵,耗时2.5分钟;调整单根索耗费时长为40s。由此可得调整5m索网天线,采用敏度矩阵方法的总时长为213.833分钟,采用最大量调整量索调整方法耗费总时长为225分钟,采用关键索调整方法耗费总时长为142分钟。
采用关键索调整天线,虽然调整次数多导致测量节点位置耗时增加,但每次调整关键索耗费时长是敏度矩阵调整全部索的百分之三十,对于大口径天线耗费总时长远远少于采用敏度矩阵调整方法所耗费的时长。
综上,本发明基于修正模型的敏度矩阵定义了质量份数作为关键调整索的选取准则,挑选部分可调整索作为关键索调整索网天线型面。其优点在于通过修正数学模型保证计算敏度矩阵和计算调整量的准确性,而选取关键索调整方法减少了每次调整型面索的数量,减少了索网天线调整耗费的总时长,和现有的敏度矩阵调整方法相比,提高了提索网天线的调整速度。

Claims (5)

1.基于选取关键索准则的索网天线型面精度调整方法,其特征在于:具体包括如下步骤:
步骤1,根据设计实物天线时的理想天线建立有限元模型作为初始数学模型,测量实物索网天线所有节点的位置
Figure FDA0004172867130000011
组合所有节点平衡方程,得到索网天线自由节点的平衡矩阵B3(n-c)×m,n为实物索网天线节点个数、m为实物索网天线索单元个数;
步骤2,计算由步骤1得到的平衡矩阵B3(n-c)×m的秩RB,得到独立自由模态数s1=m-RB,挑选出s1根索单元编号作为待测量索力索单元编号序列N1,并将实物天线索力T分为待测量索力T1和待求解索力T2,待测量索力T1和待求解索力T2相对应的平衡矩阵分块为B1、B2
步骤3,测量实物索网天线索单元编号序列为N1的待测量索力T1,通过分块平衡方程计算待求解索力T2
Figure FDA0004172867130000012
步骤4,令索网天线的调整次数k=0,记实物天线初始索力T(0)=[T1,T2];根据步骤1得到的实物索网天线所有节点的位置
Figure FDA0004172867130000013
和实物天线索力T(0)修改索网天线数学模型中相应的节点位置为/>
Figure FDA0004172867130000014
修改数学模型的索力为T(0),完成索网天线数学模型参数的修正;
步骤5,令调整次数k=k+1,在步骤4修正的索网天线数学模型型面顶点建立直角坐标系Oxyz,根据数学模型全部节点坐标
Figure FDA0004172867130000015
得到数学模型的型面的n1个标记点的坐标/>
Figure FDA0004172867130000016
计算n1个标记点的期望位移/>
Figure FDA0004172867130000021
计算索网天线型面上的标记点位置关于调整索索长的敏度矩阵/>
Figure FDA0004172867130000022
步骤6,计算可调整索的质量份数Pi
Figure FDA0004172867130000023
式中,
Figure FDA0004172867130000024
为几何逼近误差δp的平方,Δz为型面标记点的期望z向位移,/>
Figure FDA0004172867130000025
为第i根调整索对应的敏度矩阵S(k)中第i列,ΔLi为可调整索的索长变化量;
步骤7,将全部的m1根可调整索质量份数进行排序,按照质量份数由高到低地挑选出s2根可调整索作为关键索,计算第k次调整中关键索的调整量:
Figure FDA0004172867130000026
所述步骤7中,
Figure FDA0004172867130000027
为从所有的m1根可调整索中挑选s2根关键索对应的敏度矩阵中列的组合,/>
Figure FDA0004172867130000028
为/>
Figure FDA0004172867130000029
的广义逆,在第k次调整中,关键索的调整量为/>
Figure FDA00041728671300000210
其余的可调整索的调整量为0;
步骤8,将步骤7所得的关键索的调整量
Figure FDA00041728671300000211
施加到实物模型上,测量得到调整后的实物反射面节点坐标/>
Figure FDA00041728671300000212
步骤9,对实物天线关键索的调整量进行修正,修正后的关键索调整量
Figure FDA00041728671300000213
为:
Figure FDA00041728671300000214
步骤10,将步骤9修正后的关键索调整量
Figure FDA00041728671300000215
施加到数学模型上,得到修正后的数学模型标记点坐标/>
Figure FDA00041728671300000216
计算第k次调整之后的天线型面精度δ(k),至此完成第k次调整;
步骤11,重复步骤5~10,直到索网天线型面精度δ(k)≤[δ],完成索网天线关键索调整,型面误差的调整目标为[δ]。
2.根据权利要求1所述的基于选取关键索准则的索网天线型面精度调整方法,其特征在于:所述步骤1中,根据设计实物天线的理想天线建立初始有限元数学模型,实物天线由n个节点、m根索单元组成,其中c个固定节点,索网天线型面的标记节点数为n1,可调整索单元m1根,初始型面误差为δ0,型面误差的调整目标为[δ],B3(n-c)×m为索网天线自由节点的平衡矩阵,根据平衡索网天线中的任意自由节点应满足x、y、z方向的力平衡方程,其中自由节点i的x方向为:
Figure FDA0004172867130000031
式中,Ci为与自由节点i相连的节点编号的集合,索单元ij连接节点i、j,Tij为索单元ij的张力,Lij为索单元ij的长度,xi和xj分别为节点i、j的x方向坐标;
则节点i包含3个方向的力平衡方程:
Figure FDA0004172867130000032
组合所有自由节点平衡方程,则所有自由节点的平衡方程为:
B3(n-c)×mTm×1=03(n-c)×1 (7);
式中,Tm×1为索单元张力向量,(n-c)为索网天线自由节点总数。
3.根据权利要求2所述的基于选取关键索准则的索网天线型面精度调整方法,其特征在于:所述步骤2中,由于平衡矩阵的秩RB的表达式如下:
RB=3(n-c)<m (8);
基于公式(8)的不等式,所以公式(7)存在不唯一解,为了对索网天线进行唯一修正,需要挑选独立不相关的索单元测量索力,然后对天线索力进行求解修正,计算平衡矩阵的秩RB得独立自应力模态数s1=m-RB,从全部的m根索单元中随机给定s1根索单元编号作为初始待测量索力编号序列N1,其余索力编号作为待求解索力编号序列N2,并依次将索力分为待测量索力T1和待求解索力T2,相应的平衡矩阵分块为B1、B2,则平衡矩阵表示为:
Figure FDA0004172867130000041
计算分块矩阵B2的秩,若B2不满秩则将编号序列N1最后一个索单元编号放到N2最前面,将N2最后一个编号放到N1最前面,然后再计算B2的秩,直到B2满秩,则挑选完成,独立不相关待测索力编号序列为N1
4.根据权利要求3所述的基于选取关键索准则的索网天线型面精度调整方法,其特征在于:所述步骤3中,实测实物模型的索力T1,则将T1代入公式(9)中,根据公式(1)得到待求解索力T2
5.根据权利要求4所述的基于选取关键索准则的索网天线型面精度调整方法,其特征在于:所述步骤5中,数学模型索网型面上存在n1个标记点,为确定第k次调整时索网标记点期望的位移
Figure FDA0004172867130000051
需要计算索网天线的调整目标抛物面,步骤5的具体步骤为:
步骤5.1,在数学模型型面顶点建立直角坐标系Oxyz,其中,z轴沿抛物面轴向,拟合的最佳几何逼近抛物面为:
Figure FDA0004172867130000052
其中,x、y、z为索网天线抛物面上点的坐标,f为抛物面焦距,Δh1为最佳几何逼近抛物面相对于坐标原点的z向偏移量,其中几何逼近误差由数值积分计算:
Figure FDA0004172867130000053
式中,s为上网面有效口径内划分三角形单元个数,
Figure FDA0004172867130000054
为第k1个三角形面片在口径面内的投影区域,/>
Figure FDA0004172867130000055
是该面片在天线孔径方向投影面积,(Xi,Yj)为三角形面上的数值积分高斯点,F(Xi,Yj)为型面高斯点的z向偏差,W(i),W(j)为X、Y方向上的权系数/>
Figure FDA0004172867130000056
步骤5.2,求解型面标记点的调整期望位移
Figure FDA0004172867130000057
需要计算标记点调整目标抛物面,索网型面几何逼近误差达到最小,则调整目标抛物面方程为:
Figure FDA0004172867130000058
其中,Δh2为天线在设计时的几何逼近抛物面相对于设计的理想抛物面的z向偏移量,Δh1和Δh2通过一维搜索进退法求解;
则将每个型面标记点坐标代入式(1),得到型面标记点投影到调整目标抛物面的投影坐标
Figure FDA0004172867130000061
由此能够得到型面标记点的期望位移为:
Figure FDA0004172867130000062
步骤5.3,根据第k次调整时的数学模型求解敏度矩阵,从第i=1根可调整索开始,改变索网天线第i根可调整索的索力ΔT,分析数学模型的平衡态得到型面节点的位移向量ΔXi和该可调整索的索长变化量ΔLi,敏度矩矩阵的第i列为
Figure FDA0004172867130000063
将该可调整索的索力恢复改变之前的状态,对下一根可调整索进行分析,直到改变了所有的m1根可调整索,组装敏度矩阵为:
Figure FDA0004172867130000064
/>
CN202210528465.9A 2022-05-16 2022-05-16 基于选取关键索准则的索网天线型面精度调整方法 Active CN114970254B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210528465.9A CN114970254B (zh) 2022-05-16 2022-05-16 基于选取关键索准则的索网天线型面精度调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210528465.9A CN114970254B (zh) 2022-05-16 2022-05-16 基于选取关键索准则的索网天线型面精度调整方法

Publications (2)

Publication Number Publication Date
CN114970254A CN114970254A (zh) 2022-08-30
CN114970254B true CN114970254B (zh) 2023-05-23

Family

ID=82983064

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210528465.9A Active CN114970254B (zh) 2022-05-16 2022-05-16 基于选取关键索准则的索网天线型面精度调整方法

Country Status (1)

Country Link
CN (1) CN114970254B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101587026A (zh) * 2009-06-19 2009-11-25 东南大学 基于索力监测来识别需调整索力的支承索的方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103037389B (zh) * 2012-11-22 2015-04-08 西安电子科技大学 基于敏度信息的索网反射面天线型面赋形方法
CN103761368B (zh) * 2014-01-02 2017-04-05 西安电子科技大学 一种索网反射面天线形面精度与张力状态的同时设计方法
CN104866666B (zh) * 2015-05-20 2017-09-26 西安电子科技大学 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法
CN105510881B (zh) * 2015-11-20 2018-02-06 西安电子科技大学 一种静电成形薄膜天线径向索长度调整装置及调整方法
CN108598716B (zh) * 2018-04-18 2020-07-14 大连理工大学 索网天线反射面索网生成的简易方法
CN110263429B (zh) * 2019-06-18 2021-11-19 西安电子科技大学 基于机器学习的索网反射面天线形面精度调整方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101587026A (zh) * 2009-06-19 2009-11-25 东南大学 基于索力监测来识别需调整索力的支承索的方法

Also Published As

Publication number Publication date
CN114970254A (zh) 2022-08-30

Similar Documents

Publication Publication Date Title
CN106934147B (zh) 一种空间桁架索网天线等几何找形方法
CN106020147A (zh) 一种复杂精密机械产品的装配精度系统分析方法
CN103712557B (zh) 面向特大型齿轮的激光跟踪多站位定位方法
CN110263429B (zh) 基于机器学习的索网反射面天线形面精度调整方法
CN108801218B (zh) 大尺寸动态摄影测量系统的高精度定向及定向精度评价方法
CN112539906B (zh) 挠性喷管型面的静态调试方法
CN108416083B (zh) 一种高耸电视塔结构二维动力模型分析方法及系统
CN110686655B (zh) 三维曲面工件的定位方法、装置及系统
CN109597962A (zh) 一种基于预测和迭代补偿的参数插值方法
CN114970254B (zh) 基于选取关键索准则的索网天线型面精度调整方法
CN111907729A (zh) 一种针对飞机大部件调姿的最佳位姿拟合方法
CN109902350B (zh) 对变截面梁的截面惯性矩进行模型修正中克服模态交换的方法
CN113221278A (zh) 一种车载光电平台照准架轻量化方法
CN108959702A (zh) 一种基于面偏差传递的火箭贮箱偏差分析方法
CN113932707B (zh) 基于弹性网络算法补偿三坐标测量机几何误差的方法
CN113954102B (zh) 一种基于离线编程的百叶轮抛磨叶片路径规划方法
CN110096798A (zh) 一种多状态有限元模型修正的方法
CN115618617A (zh) 一种基于地形高程匹配的落点预报方法
CN110110396B (zh) 有限元模型修正过程中克服模态交换的方法
CN113901690A (zh) 星载反射面天线在轨热变形性能评估方法
CN111474579B (zh) 一种基于覆盖次数均匀的自动变观方法
CN114707380B (zh) 一种亚毫米波望远镜主反射面促动器调整和精度计算方法
CN112307562A (zh) 综合热变形与重力变形的大型飞机上复杂部件的装配方法
KR101888794B1 (ko) 선체 설계방법
CN117428363B (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