CN107729705B - 一种面天线单块面板精度的测算方法 - Google Patents
一种面天线单块面板精度的测算方法 Download PDFInfo
- Publication number
- CN107729705B CN107729705B CN201711225543.3A CN201711225543A CN107729705B CN 107729705 B CN107729705 B CN 107729705B CN 201711225543 A CN201711225543 A CN 201711225543A CN 107729705 B CN107729705 B CN 107729705B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- panel
- axis
- local coordinate
- theoretical design
- 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
Links
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Nitrogen And Oxygen Or Sulfur-Condensed Heterocyclic Ring Systems (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种面天线单块面板精度的测算方法,属于天线技术领域。该方法包括建立面板局部坐标系、建立测量坐标系与面板局部坐标系的转换关系、建立面板局部坐标系与理论设计坐标系的转换关系,以及计算面板精度等步骤。本发明提出的方法可用于计算测量或仿真的面板精度,为理论设计工作提供重要指导,是对现有技术的一种重要改进。
Description
技术领域
本发明涉及天线技术领域,特别是指一种面天线单块面板精度的测算方法。
背景技术
反射面天线的面形精度是影响天线性能的重要因素,尤其对于大口径、高频段天线而言,实现主反射面面型的高精度成为其研制成功的关键技术之一。单块面板的精度对主反射面面型起着至关重要的影响,能否研制出高精度的面板,是天线研制的技术瓶颈之一。
随着卫星通信、射电天文和深空探测的发展需求,面天线正向高频段、大口径方向发展,为了实现单块面板的高精度,面板设计的尺寸不宜过大,然而小尺寸面板将增加天线反射器重量、座架重量和天线转动惯量,从而恶化天线动态性能,同时增加了制造成本,由此可见,大尺寸高精度面板的设计是现有技术中迫切需要实现的技术。特别地,对于面板精度的测算是高精度面板设计的一个重要组成部分,但是,现有技术中尚缺少面板精度测算的有关技术。
发明内容
有鉴于此,本发明提供一种面天线单块面板精度的测算方法,该方法具有测算过程简单、测算结果准确的特点。
为了实现上述目的,本发明所采用的技术方案如下:
一种面天线单块面板精度的测算方法,其包括以下步骤:
(1)建立面板局部坐标系x′o′y′z′:
在面板上选取3个不共线的点P1、P2、P3,以P1点为局部坐标系x′o′y′z′的原点o′,以P1、P2点的连线作为局部坐标系的x′轴,以P1、P2、P3三点所确定的平面作为局部坐标系x′o′y′z′的x′o′y′平面,以过P1点且与x′o′y′平面相垂直的直线为z′轴;
(2)建立测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换关系;
(3)建立面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换关系;
(4)计算面板精度δm:
根据测量坐标系与面板局部坐标系的转换关系以及面板局部坐标系与理论设计坐标系的转换关系,以面板局部坐标系为公共坐标系,将面板在测量坐标系下的测量数值转换到理论坐标系下,计算转换后的测量值与理论设计值之间的偏差δm,偏差δm即面板精度。
可选的,所述步骤(2)的具体步骤为:
(201)求解面板局部坐标系x′o′y′z′中x′轴的单位向量ex′:
(202)求解面板局部坐标系x′o′y′z′中z′轴的单位向量ez′:
(203)求解面板局部坐标系x′o′y′z′中y′轴的单位向量ey′:
(204)将测量坐标系x1o1y1z1的原点平移至点P1处,并使o1x1轴与o′x′轴重合、o1y1轴与o′y′轴重合、o1z1轴与o′z′轴重合,结合步骤(201)~(203),求得测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换关系为:
以上各式中,(x1p1,y1p1,z1p1)、(x1p2,y1p2,z1p2)和(x1p3,y1p3,z1p3)分别表示面板上的3个点P1、P2、P3在测量坐标系下的坐标,为从测量坐标系x1o1y1z1到面板局部坐标系x′o′y′z′的转换矩阵,(x1,y1,z1)和(x′,y′,z′)分别表示同一点在测量坐标系和面板局部坐标系下的坐标,a01、a02、a03为中间变量,分别表示测量坐标系x1o1y1z1的o1x1轴单位向量、o1y1轴单位向量和o1z1轴单位向量。
可选的,所述步骤(3)的具体步骤为:
(301)求解面板局部坐标系x′o′y′z′中x′轴的单位向量ex′:
(302)求解面板局部坐标系x′o′y′z′中z′轴的单位向量ez′:
(303)求解面板局部坐标系x′o′y′z′中y′轴的单位向量ey′:
(304)将面板局部坐标系x′o′y′z′的原点P1平移至理论设计坐标系xoyz的原点处,并使o′x′轴与ox轴重合、o′y′轴与oy轴重合、o′z′轴与oz轴重合,结合步骤(301)~(303),求得面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换关系为:
以上各式中,(xp1,yp1,zp1)、(xp2,yp2,zp2)、(xp3,yp3,zp3)分别为面板上的3个点P1、P2、P3在理论设计坐标系下的坐标,为从面板局部坐标系x′o′y′z′到理论设计坐标系xoyz的转换矩阵,(x,y,z)和(x′,y′,z′)分别表示同一点在理论设计坐标系和面板局部坐标系下的坐标,a′01、a′02、a′03为中间变量,分别表示理论设计坐标系xoyz的ox轴单位向量、oy轴单位向量和oz轴单位向量。
可选的,所述步骤(4)中δm的计算方式为:
(401)根据测量坐标系与面板局部坐标系的转换关系以及面板局部坐标系与理论设计坐标系的转换关系,建立测量坐标系x1o1y1z1与理论设计坐标系xoyz的转换关系:
其中,(xp1,yp1,zp1)和(x1p1,y1p1,z1p1)分别表示面板上的点P1在理论设计坐标系和测量坐标系下的坐标,为从测量坐标系x1o1y1z1到面板局部坐标系x′o′y′z′的转换矩阵y,为从面板局部坐标系x′o′y′z′到理论设计坐标系xoyz的转换矩阵,(x,y,z)和(x1,y1,z1)分别表示同一点在理论设计坐标系和测量坐标系下的坐标;
(402)根据步骤(401)的转换关系,将面板在测量坐标系下的测量值(xpn,ypn,zpn)转换为理论设计坐标系下的坐标值(x′m,y′m,z′m),转换后的径向值为r′m:
(403)选取理论设计径向值rm与转换后的径向值r′m相同的点,计算轴向偏差:Δz=zm-z′m;
其中,理论设计坐标系xoyz下面板的坐标为(xm,ym,zm),径向值为rm:
(404)根据法向偏差与轴向偏差的转换关系δ=Δz×cos(θv/2),得到面板精度δm为:
本发明与背景技术相比具有如下优点:
1、本发明测算方法可用于计算变形面板的精度,对理论设计提供重要指导。
2、本发明可用于计算力学仿真软件分析的面板形变、测量的面板形变精度。
附图说明
图1是本发明实施例中理论设计坐标系xoyz、面板局部坐标系x′o′y′z′及测量坐标系x1o1y1z1的对应关系原理图。
图2是本发明实施例中理论设计坐标系下理论面与测量面的对应关系原理图。
具体实施方式
如图1和2所示,一种面天线单块面板精度的测算方法,其包括如下步骤:
(1)建立面板局部坐标系x′o′y′z′
参照图1,选取面板上的3个点P1、P2、P3,其中以P1点为局部坐标系x′o′y′z′的原点o′,P1、P2点确定的直线为局部坐标系的为x′轴,P1、P2、P3三点确定的平面为局部坐标系x′o′y′z′的x′o′y′平面,通过P1点与x′o′y′平面垂直的为z′轴,确定x′轴、z′轴后通过右手法则可确定y′轴。
(2)建立测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换关系
参照图1,测量坐标系为x1o1y1z1,面板在测量坐标系下的坐标为(x1pn,y1pn,z1pn),选取的面板上3个点P1、P2、P3在测量坐标系下的坐标分别为(x1p1,y1p1,z1p1),(x1p2,y1p2,z1p2),(x1p3,y1p3,z1p3)。
(201)求解x′轴的单位向量ex′:
(202)求解z′轴的单位向量ez′:
ez′的求解公式为:
(203)求解y′轴的单位向量ey′:
(204)建立测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换矩阵
将测量坐标系x1o1y1z1的原点平移至P1,o1x1轴与o′x′轴重合,o1y1轴与o′y′轴重合,o1z1轴与o′z′轴重合,结合步骤(201)~(203),可求得测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换关系为:
(3)建立面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换关系
参照图1,理论设计坐标系为xoyz,面板在理论设计坐标系下坐标为(xpn,ypn,zpn),选取的面板上3个点P1、P2、P3在理论设计坐标系下的坐标点为(xp1,yp1,zp1),(xp2,yp2,zp2),(xp3,yp3,zp3)。
(301)求解x′轴的单位向量ex′:
(302)求解z′轴的单位向量ez′:
ez′的求解公式为:
(303)求解y′轴的单位向量ey′:
(304)建立面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换矩阵
将面板局部坐标系x′o′y′z′的原点P1平移至o,o′x′轴与ox轴重合,o′y′轴与oy轴重合,o′z′轴与oz轴重合,结合步骤(301)~(303),可求得面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换关系为:
(4)计算面板精度δm
参照图1,根据步骤(2)、(3),以面板局部坐标系为公共坐标系,依据公共坐标系可将面板在测量坐标系下的数值转换到理论坐标系下,转换的过程为面板在测量坐标系下的数值转换到面板局部坐标系,再由面板局部坐标系转换到理论坐标系,便可计算在理论坐标系下面板的精度值。
(401)建立测量坐标系x1o1y1z1与理论设计坐标系xoyz的转换关系
参照图1,以面板局部坐标系为公共坐标系,结合步骤(2)、(3),测量坐标系x1o1y1z1与理论设计坐标系xoyz的转换关系为:
(402)参照图2,结合步骤(401),面板在测量坐标系下的测量值(xpn,ypn,zpn)转换到理论设计坐标系xoyz的数值为(x′m,y′m,z′m),转换后的径向值为r′m:
(403)求解面板精度δm:
参照图2,理论设计坐标系xoyz下面板的坐标为(xm,ym,zm),径向值为rm:
参照图2,选取转换后的径向值r′m与理论设计径向值rm相同的点,则轴向偏差的计算公式为:
Δz=zm-z′m。
法向偏差与轴向偏差的转换关系为:
δ=Δz×cos(θv/2)
参照图2,式中:θv为赋形设计时主副面相应点的连线与z轴的夹角。
面板精度δm计算公式为:
式中:
N为测量的数据点数。
至此,完成单块面板精度的计算。
该方法具有测算过程简单、测算结果准确的特点,可用于计算变形面板的精度,并可用于计算力学仿真软件分析的面板形变以及测量的面板形变精度,为理论设计提供重要指导,是对现有技术的一种重要改进。
Claims (3)
1.一种面天线单块面板精度的测算方法,其特征在于,包括以下步骤:
(1)建立面板局部坐标系x′o′y′z′:
在面板上选取3个不共线的点P1、P2、P3,以P1点为局部坐标系x′o′y′z′的原点o′,以P1、P2点的连线作为局部坐标系的x′轴,以P1、P2、P3三点所确定的平面作为局部坐标系x′o′y′z′的x′o′y′平面,以过P1点且与x′o′y′平面相垂直的直线为z′轴;
(2)建立测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换关系;
(3)建立面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换关系;
(4)计算面板精度δm:
根据测量坐标系与面板局部坐标系的转换关系以及面板局部坐标系与理论设计坐标系的转换关系,以面板局部坐标系为公共坐标系,将面板在测量坐标系下的测量数值转换到理论坐标系下,计算转换后的测量值与理论设计值之间的偏差δm,偏差δm即面板精度;
所述步骤(4)的具体方式为:
(401)根据测量坐标系与面板局部坐标系的转换关系以及面板局部坐标系与理论设计坐标系的转换关系,建立测量坐标系x1o1y1z1与理论设计坐标系xoyz的转换关系:
其中,(xp1,yp1,zp1)和(x1p1,y1p1,z1p1)分别表示面板上的点P1在理论设计坐标系和测量坐标系下的坐标,为从测量坐标系x1o1y1z1到面板局部坐标系x′o′y′z′的转换矩阵,为从面板局部坐标系x′o′y′z′到理论设计坐标系xoyz的转换矩阵,(x,y,z)和(x1,y1,z1)分别表示同一点在理论设计坐标系和测量坐标系下的坐标;
(402)根据步骤(401)的转换关系,将面板在测量坐标系下的测量值(xpn,ypn,zpn)转换为理论设计坐标系下的坐标值(x′m,y′m,z′m),转换后的径向值为r′m:
(403)选取理论设计径向值rm与转换后的径向值r′m相同的点,计算轴向偏差:Δz=zm-z′m;
其中,理论设计坐标系xoyz下面板的坐标为(xm,ym,zm),径向值为rm:
(404)根据法向偏差与轴向偏差的转换关系δ=Δz×cos(θv/2),得到面板精度δm为:
2.根据权利要求1所述的面天线单块面板精度的测算方法,其特征在于,所述步骤(2)的具体步骤为:
(201)求解面板局部坐标系x′o′y′z′中x′轴的单位向量ex′:
(202)求解面板局部坐标系x′o′y′z′中z′轴的单位向量ez′:
(203)求解面板局部坐标系x′o′y′z′中y′轴的单位向量ey′:
(204)将测量坐标系x1o1y1z1的原点平移至点P1处,并使o1x1轴与o′x′轴重合、o1y1轴与o′y′轴重合、o1z1轴与o′z′轴重合,结合步骤(201)~(203),求得测量坐标系x1o1y1z1与面板局部坐标系x′o′y′z′的转换关系为:
3.根据权利要求1所述的面天线单块面板精度的测算方法,其特征在于,所述步骤(3)的具体步骤为:
(301)求解面板局部坐标系x′o′y′z′中x′轴的单位向量ex′:
(302)求解面板局部坐标系x′o′y′z′中z′轴的单位向量ez′:
(303)求解面板局部坐标系x′o′y′z′中y′轴的单位向量ey′:
(304)将面板局部坐标系x′o′y′z′的原点P1平移至理论设计坐标系xoyz的原点处,并使o′x′轴与ox轴重合、o′y′轴与oy轴重合、o′z′轴与oz轴重合,结合步骤(301)~(303),求得面板局部坐标系x′o′y′z′与理论设计坐标系xoyz的转换关系为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711225543.3A CN107729705B (zh) | 2017-11-29 | 2017-11-29 | 一种面天线单块面板精度的测算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711225543.3A CN107729705B (zh) | 2017-11-29 | 2017-11-29 | 一种面天线单块面板精度的测算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107729705A CN107729705A (zh) | 2018-02-23 |
CN107729705B true CN107729705B (zh) | 2020-10-16 |
Family
ID=61220066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711225543.3A Active CN107729705B (zh) | 2017-11-29 | 2017-11-29 | 一种面天线单块面板精度的测算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107729705B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110345923B (zh) * | 2018-04-08 | 2021-06-18 | 孟艳艳 | 天线主副反射面位姿测量系统与副反射面位姿调整方法 |
CN110686615B (zh) * | 2019-08-29 | 2022-01-04 | 西安空间无线电技术研究所 | 一种高精度伞状天线型面评价方法 |
CN111046463B (zh) * | 2019-11-28 | 2022-11-22 | 郑州大学 | 基于正交理论的桁架结构变形分解方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013775A (zh) * | 2006-12-22 | 2007-08-08 | 西安电子科技大学 | 基于多维无约束优化的天线表面精度调整方法 |
CN103258075A (zh) * | 2013-03-25 | 2013-08-21 | 中国电子科技集团公司第五十四研究所 | 一种赋形卡式天线主面精度和主副面调整的精确计算方法 |
CN104866666A (zh) * | 2015-05-20 | 2015-08-26 | 西安电子科技大学 | 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法 |
CN105740554A (zh) * | 2016-02-02 | 2016-07-06 | 西安电子科技大学 | 一种面向增益的大型变形抛物面天线面板精度调整方法 |
CN107240780A (zh) * | 2017-05-27 | 2017-10-10 | 西安电子科技大学 | 基于面片积分公式的伞状天线结构优化设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2378444B1 (de) * | 2010-04-13 | 2015-07-29 | CST-Computer Simulation Technology AG | Verfahren, Vorrichtung und Computerprogrammprodukt zur Bestimmung eines elektromagnetischen Nahfeldes einer Feldanregungsquelle eines elektrischen Systems |
-
2017
- 2017-11-29 CN CN201711225543.3A patent/CN107729705B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013775A (zh) * | 2006-12-22 | 2007-08-08 | 西安电子科技大学 | 基于多维无约束优化的天线表面精度调整方法 |
CN103258075A (zh) * | 2013-03-25 | 2013-08-21 | 中国电子科技集团公司第五十四研究所 | 一种赋形卡式天线主面精度和主副面调整的精确计算方法 |
CN104866666A (zh) * | 2015-05-20 | 2015-08-26 | 西安电子科技大学 | 一种基于有限元模型修正的索网反射面形面精度的鲁棒调整方法 |
CN105740554A (zh) * | 2016-02-02 | 2016-07-06 | 西安电子科技大学 | 一种面向增益的大型变形抛物面天线面板精度调整方法 |
CN107240780A (zh) * | 2017-05-27 | 2017-10-10 | 西安电子科技大学 | 基于面片积分公式的伞状天线结构优化设计方法 |
Non-Patent Citations (1)
Title |
---|
赋形卡式天线最佳吻合反射面的计算方法;闫丰 等;《电磁场与微波》;20111231;第38-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107729705A (zh) | 2018-02-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107729705B (zh) | 一种面天线单块面板精度的测算方法 | |
CN106020147B (zh) | 一种复杂精密机械产品的装配精度系统分析方法 | |
CN100437601C (zh) | 一种用于自动钻铆托架变形的修补数据处理方法 | |
CN105108215B (zh) | 一种自由曲面微细铣削让刀误差预测及补偿方法 | |
CN105426566B (zh) | 基于蒙特卡罗技术的大型高速回转装备形位公差分配方法 | |
CN101989682B (zh) | 天线反射面拟合方法 | |
CN105426565B (zh) | 基于概率密度技术的大型高速回转装备形位公差分配方法 | |
CN110052881B (zh) | 一种复杂零件多工位转换加工基准找正方法 | |
CN111780698B (zh) | 一种工件坐标系的标定方法及相关装置 | |
CN106405581B (zh) | 卫星结构变形引起的多类型载荷间耦合指向精度评估方法 | |
Zheng et al. | An efficient method for minimum zone cylindricity error evaluation using kinematic geometry optimization algorithm | |
CN104759964A (zh) | 光学非球面元件的变形加工方法 | |
CN111046584B (zh) | 星载仪器设备的精密调整方法 | |
CN101847262A (zh) | 一种快速三维点云搜索匹配方法 | |
CN106991210A (zh) | 一种基于机电耦合模型的赋形反射面天线电性能预测方法 | |
CN115046511A (zh) | 一种叶片周向圆弧锤足型榫头三坐标测量坐标系构建方法 | |
CN111891176A (zh) | 基于轨道几何状态测量仪的轨道测量数据整体平差方法 | |
CN112129322A (zh) | 一种捷联惯组与三轴转台的安装误差检测方法及校正方法 | |
CN106353033A (zh) | 一种航空发动机质心计算方法 | |
CN107066726A (zh) | 一种数控机床旋转轴垂直度误差建模方法 | |
CN103258075B (zh) | 一种赋形卡式天线主面精度和主副面调整的精确计算方法 | |
JP5001857B2 (ja) | アンテナアレーの向きを決定するための方法 | |
CN112631200A (zh) | 一种机床轴线测量方法以及装置 | |
JP2008528957A5 (zh) | ||
Zhong et al. | Modeling variation propagation in machining systems with different configurations |
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 |