发明名称 一种赋形卡式天线主面精度和主副面调整的精确计算方法
摘要 本发明涉及一种赋形卡式天线主面精度和主副面调整的精确计算方法,属于天线技术领域,可用于精确计算赋形卡式天线的主面精度,同时可准确指导副反射面位置姿态的调整和主反射面面板的主动调整,尤其适用于大型射电望远镜天线领域。本发明的思路是基于等光程条件,通过分别沿采用x轴、y轴、z轴方向的平移和绕x<sub>1</sub>轴、y<sub>2</sub>轴的旋转及光程变化量6个自由度的最小二乘拟合计算最佳吻合反射面,从而可求得赋形卡式天线的主面精度。同时,本发明计算的结果可以指导副反射面位置姿态的调整和主反射面面板的主动调整,从而可减小主反射面的变形对天线性能的影响,提高天线效率。
申请公布号 CN103258075B 申请公布日期 2016.05.25
申请号 CN201310095282.3 申请日期 2013.03.25
申请人 中国电子科技集团公司第五十四研究所 发明人 闫丰;杜彪;刘国玺;周建寨
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 一种赋形卡式天线主面精度和主副面调整的精确计算方法,其特征在于包括以下步骤:(1)假设,赋形卡式天线理论设计主副反射面坐标系为xoyz;将坐标系xoyz分别沿着x轴、y轴、z轴方向进行平移,之后再分别绕x<sub>1</sub>轴、y<sub>2</sub>轴、z<sub>3</sub>轴方向进行旋转,得到描述移动后的副面和与移动后的副面相对应的最佳吻合反射面的坐标系x<sub>4</sub>o<sub>4</sub>y<sub>4</sub>z<sub>4</sub>;所述的x<sub>1</sub>轴指的是坐标系xoyz分别沿着x轴、y轴、z轴方向进行平移后形成的坐标系x<sub>1</sub>o<sub>1</sub>y<sub>1</sub>z<sub>1</sub>的横轴;y<sub>2</sub>轴指的是坐标系x<sub>1</sub>o<sub>1</sub>y<sub>1</sub>z<sub>1</sub>绕x<sub>1</sub>轴旋转后形成的坐标系x<sub>2</sub>o<sub>2</sub>y<sub>2</sub>z<sub>2</sub>的纵轴;z<sub>3</sub>轴指的是坐标系x<sub>2</sub>o<sub>2</sub>y<sub>2</sub>z<sub>2</sub>绕y<sub>2</sub>轴旋转后形成的坐标系x<sub>3</sub>o<sub>3</sub>y<sub>3</sub>z<sub>3</sub>的竖轴;x<sub>4</sub>o<sub>4</sub>y<sub>4</sub>z<sub>4</sub>坐标系的表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>x</mi><mn>4</mn><mo>&prime;</mo></msubsup><mo>=</mo><msup><mi>x</mi><mo>&prime;</mo></msup><mo>-</mo><msub><mi>u</mi><mi>x</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mi>z</mi></msub><msup><mi>y</mi><mo>&prime;</mo></msup><mo>-</mo><msub><mi>&alpha;</mi><mi>y</mi></msub><msup><mi>z</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>y</mi><mn>4</mn><mo>&prime;</mo></msubsup><mo>=</mo><msup><mi>y</mi><mo>&prime;</mo></msup><mo>-</mo><msub><mi>u</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>z</mi></msub><msup><mi>x</mi><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&alpha;</mi><mi>x</mi></msub><msup><mi>z</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>z</mi><mn>4</mn><mo>&prime;</mo></msubsup><mo>=</mo><msup><mi>z</mi><mo>&prime;</mo></msup><mo>-</mo><msub><mi>u</mi><mi>z</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mi>y</mi></msub><msup><mi>x</mi><mo>&prime;</mo></msup><mo>-</mo><msub><mi>&alpha;</mi><mi>x</mi></msub><msup><mi>y</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FSB0000150074520000011.GIF" wi="499" he="231" /></maths>式中:x′<sub>4</sub>、y′<sub>4</sub>、z′<sub>4</sub>分别为x<sub>4</sub>o<sub>4</sub>y<sub>4</sub>z<sub>4</sub>坐标系下的x<sub>4</sub>轴、y<sub>4</sub>轴、z<sub>4</sub>轴的坐标值;x′、y′、z′分别为xoyz坐标系下的x轴、y轴、z轴的坐标值;u<sub>x</sub>、u<sub>y</sub>、u<sub>z</sub>分别为沿x轴、y轴、z轴的平移量;α<sub>x</sub>、α<sub>y</sub>、α<sub>z</sub>分别为绕x<sub>1</sub>轴、y<sub>2</sub>轴、z<sub>3</sub>轴的旋转角度;(2)由于荷载因素导致理论设计主反射面产生了变形,通过测量或力学有限元软件分析得到变形后的主反射面表达式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><msub><mi>x</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&prime;</mo></msup><mo>=</mo><msub><mi>x</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>x</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msup><msub><mi>y</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&prime;</mo></msup><mo>=</mo><msub><mi>y</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>y</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msup><msub><mi>z</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&prime;</mo></msup><mo>=</mo><msub><mi>z</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>z</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FSB0000150074520000012.GIF" wi="313" he="288" /></maths>式中:x′<sub>m0</sub>、y′<sub>m0</sub>、z′<sub>m0</sub>分别为变形后的主反射面在xoyz坐标系下的x轴、y轴、z轴的坐标值;x<sub>m</sub>、y<sub>m</sub>、z<sub>m</sub>分别为理论设计主反射面在xoyz坐标系下的x轴、y轴、z轴的坐标值;ε<sub>y</sub>、ε<sub>y</sub>、ε<sub>z</sub>分别为理论设计主反射面在xoyz坐标系下沿x轴、y轴、z轴的变形量;(3)假设,最佳吻合反射面与变形后的主反射面之间的法向偏差为δ,结合步骤(2)中变形后的主反射面表达式,计算最佳吻合反射面表达式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><msub><mi>x</mi><mi>m</mi></msub><mo>&prime;</mo></msup><mo>=</mo><msup><msub><mi>x</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>n</mi><mi>x</mi></msub><mi>&delta;</mi><mo>=</mo><msub><mi>x</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>x</mi></msub><mo>+</mo><msub><mi>n</mi><mi>x</mi></msub><mi>&delta;</mi></mrow></mtd></mtr><mtr><mtd><mrow><msup><msub><mi>y</mi><mi>m</mi></msub><mo>&prime;</mo></msup><mo>=</mo><msup><msub><mi>y</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>n</mi><mi>y</mi></msub><mi>&delta;</mi><mo>=</mo><msub><mi>y</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>y</mi></msub><mo>+</mo><msub><mi>n</mi><mi>y</mi></msub><mi>&delta;</mi></mrow></mtd></mtr><mtr><mtd><mrow><msup><msub><mi>z</mi><mi>m</mi></msub><mo>&prime;</mo></msup><mo>=</mo><msup><msub><mi>z</mi><mrow><mi>m</mi><mn>0</mn></mrow></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>n</mi><mi>z</mi></msub><mi>&delta;</mi><mo>=</mo><msub><mi>z</mi><mi>m</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><mi>z</mi></msub><mo>+</mo><msub><mi>n</mi><mi>z</mi></msub><mi>&delta;</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FSB0000150074520000021.GIF" wi="655" he="289" /></maths>式中:x′<sub>m</sub>、y′<sub>m</sub>、z′<sub>m</sub>分别为最佳吻合反射面在xoyz坐标系下的x轴、y轴、z轴的坐标值;n<sub>x</sub>、n<sub>y</sub>、n<sub>z</sub>分别为理论设计主反射面的法向向量在x轴、y轴、z轴方向的分量;δ为变形后的主反射面与最佳吻合反射面之间的法向偏差;(4)计算移动后的副面表达式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>x</mi><mi>s</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>x</mi><mi>s</mi></msub><mo>+</mo><msub><mi>u</mi><mi>x</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>z</mi></msub><msub><mi>y</mi><mi>s</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mi>y</mi></msub><msub><mi>z</mi><mi>s</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>y</mi><mi>s</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>y</mi><mi>s</mi></msub><mo>+</mo><msub><mi>u</mi><mi>y</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mi>z</mi></msub><msub><mi>x</mi><mi>s</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>x</mi></msub><msub><mi>z</mi><mi>s</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>z</mi><mi>s</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>z</mi><mi>s</mi></msub><mo>+</mo><msub><mi>u</mi><mi>z</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>y</mi></msub><msub><mi>x</mi><mi>s</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mi>x</mi></msub><msub><mi>y</mi><mi>s</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FSB0000150074520000022.GIF" wi="552" he="248" /></maths>式中:x<sub>s</sub>、y<sub>s</sub>、z<sub>s</sub>分别为理论设计副面在xoyz坐标系下的x轴、y轴、z轴的坐标值;x′<sub>s</sub>、y′<sub>s</sub>、z′<sub>s</sub>分别为移动后的副面在xoyz坐标系下的x轴、y轴、z轴的坐标值;(5)依据步骤(3)、(4)中确定的最佳吻合反射面表达式和移动后的副面表达式,选取x<sub>4</sub>o<sub>4</sub>y<sub>4</sub>为等光程参考面,计算最佳吻合反射面与变形后的主反射面之间的法向偏差δ的表达式:<img file="FSB0000150074520000023.GIF" wi="1663" he="199" /><img file="FSB0000150074520000031.GIF" wi="1419" he="158" />式中:θ是赋形设计时馈源与副面相应点的连线与z轴的夹角;θ<sub>v</sub>是赋形设计时主副面相应点的连线与z轴的夹角;<img file="FSB0000150074520000032.GIF" wi="42" he="47" />为方位角;r<sub>s</sub>为原馈源和理论设计副面之间的距离,<img file="FSB0000150074520000033.GIF" wi="402" he="95" />光程差Δ<sub>ck</sub>=C′<sub>k</sub>‑C<sub>k</sub>。其中,C<sub>k</sub>=r<sub>s</sub>+S‑z<sub>m</sub>,C<sub>k</sub>为卡式天线赋形时理论设计的等光程常数,S为理论设计主副反射面之间的距离,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>m</mi></msub><mo>-</mo><msub><mi>x</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>m</mi></msub><mo>-</mo><msub><mi>y</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>m</mi></msub><mo>-</mo><msub><mi>z</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>;</mo></mrow>]]></math><img file="FSB0000150074520000034.GIF" wi="793" he="86" /></maths>(6)依据步骤(5)中确定的法向偏差δ的表达式,计算法向偏差δ的加权均方值即赋形卡式天线主面精度δ<sub>m</sub>的表达式;(7)依据步骤(6)中确定的赋形卡式天线主面精度δ<sub>m</sub>的表达式,应用最小二乘法计算得到求解最佳吻合反射面的6个自由度参数。所述的最佳吻合反射面的6个自由度参数具体为:分别沿x轴、y轴、z轴的平移量、分别绕x<sub>1</sub>轴、y<sub>2</sub>轴的旋转量和光程变化量;(8)将步骤(7)中得到的最佳吻合反射面的6个自由度参数代入步骤(5)中得到的法向偏差δ的表达式,得到法向偏差δ的数值;再将法向偏差δ的数值代入步骤(6)中得到的赋形卡式天线主面精度δ<sub>m</sub>的表达式,得到赋形卡式天线主面精度数值,完成赋形卡式天线主面精度的计算;(9)将步骤(7)中得到的最佳吻合反射面的5个自由度参数代入步骤(4)中得到的移动后的副面表达式,便可得到副面调整的精确位置,完成副面的调整。所述的最佳吻合反射面的5个自由度参数具体为分别沿x轴、y轴、z轴的平移量和分别绕x<sub>1</sub>轴、y<sub>2</sub>轴的旋转量;(10)将步骤(7)中得到的最佳吻合反射面的6个自由度参数代入步骤(5)中得到的法向偏差δ表达式,得到法向偏差δ的数值即将变形的主反射面调整为最佳吻合反射面的调整量,完成主面的调整。
地址 050081 河北省石家庄市中山西路589号第五十四所天伺部