CN107863599A - 基于椭球模型的可展开天线展开臂容差设计方法 - Google Patents

基于椭球模型的可展开天线展开臂容差设计方法 Download PDF

Info

Publication number
CN107863599A
CN107863599A CN201711036984.9A CN201711036984A CN107863599A CN 107863599 A CN107863599 A CN 107863599A CN 201711036984 A CN201711036984 A CN 201711036984A CN 107863599 A CN107863599 A CN 107863599A
Authority
CN
China
Prior art keywords
mtd
mrow
msub
mtr
error
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
CN201711036984.9A
Other languages
English (en)
Other versions
CN107863599B (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
Publication of CN107863599A publication Critical patent/CN107863599A/zh
Application granted granted Critical
Publication of CN107863599B publication Critical patent/CN107863599B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q1/00—Details of, or arrangements associated with, antennas
    • H01Q1/08—Means for collapsing antennas or parts thereof
    • 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
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q1/00—Details of, or arrangements associated with, antennas
    • H01Q1/08—Means for collapsing antennas or parts thereof
    • H01Q1/084—Pivotable antennas

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Aerials With Secondary Devices (AREA)

Abstract

本发明涉及一种基于相关椭球模型的可展开天线展开臂关节的容差设计方法。该方法首先建立从展开臂端部反射器到卫星坐标系的坐标转换矩阵；同时通过调研分析得到容差‑成本关系。针对展开臂整体结构形式，将多级展开臂中各种误差按照相关性进行分类；同时针对内部相关的每一类误差，通过试验得到表征其相关性的椭球数学模型；继而在该椭球中进行蒙特卡洛随机抽样，基于上述坐标转换矩阵得到最终反射器的角度误差。最终，以总造价极小化为目标，以容差为设计变量，以指向误差满足设计要求为约束，进行容差优化设计。仿真结构表明，本方法可以有效的进行展开臂关节的容差设计，同时，避免了传统方法未考虑相关性时所带来的设计困难。

Description

基于椭球模型的可展开天线展开臂容差设计方法
技术领域
本发明属于可展开天线的关节容差设计技术领域，具体涉及一种基于椭球模型的可展开天线展开臂容差设计方法，用于指导天线展开臂关节的选型。
背景技术
近年来，空间可展开天线在资源探测、气象监测、空间通信等领域得到了越来越广泛的应用。可展开天线在发射时处于收拢状态并被固定在火箭整流罩内，到达预定轨道后，在地面控制指令下自动展开至工作状态。而此展开过程主要指的就是展开机构的展开。考虑到各级展开机构都存在着展开角度误差，导致展开后的连接臂的位置与理想位置存在一定的偏差，而连接臂和反射器是固连在一起，进而会对反射器指向精度产生不利的影响。考虑到展开天线的服役性能的好坏在最终很大程度上取决于反射器的指向精度，所以，分析展开臂展开角度误差对反射器指向精度的影响机理是很有必要的。
同时，考虑到铰链的存在，展开臂关节误差在各个方向的数值是相关的。一般意义上来讲，对于一个铰链三个方向的误差，当其中一个取极大值时，另外两个取极大值的可能性将减小；而针对统一制造商、同类驱动器、同种工作环境下的不同展开铰链，其取值则往往统一地偏向某一个方向。而在传统的设计方法思路中，这些不同方向误差之间的相关性并没有考虑到反射器指向误差的分析流程中去。一般的处理方法是将各个方向上误差可能出现的最大值带入到系统中去，分析此时反射器出现的最大指向偏差。而最终只要保证此最大偏差在系统设计要求之内即可。严格来说，考虑到不同方向误差之间相关性的存在，传统方法中误差的最大值集合，即为最差情况组合，并不会出现。所以，这种方法所带来的分析结果将过于保守。
因而，从误差相关性的角度出发，分析展开臂展开角度误差对反射器指向精度的影响机理是很有工程意义的。同时，在此基础之上，寻求在系统设计指标要求(反射器指向精度指标)下，展开臂合理的容差设计也是很有必要的。因为过大的容差会导致设计指标的不满足，而过小的容差会导致制造、装配环节的费用增加和加工难度。针对这个研究课题，目前在展开天线领域，现有的文献和资料均没有给出响应的解决方案。
发明内容
本发明的目的是提供一种基于椭球模型的可展开天线展开臂容差设计方法，避免上述现有技术的不足，提供一种考虑误差相关性下的可展开天线展开臂关节的容差设计方法，指导此类天线展开臂的设计，提高天线的展开精度。
实现本发明目的的技术方案是，基于椭球模型的可展开天线展开臂容差设计方法，其特征是至少包括如下步骤：
1)获取所涉及展开臂类型的容差-成本关系Cost＝C(Δ)，其中，Δ＝|ΔXmax|为各种角度误差的容差；Cost为成本；C为关系函数；
2)建立从展开臂根部到端部的不包含误差的坐标转换矩阵T，得到各个关节的展开角度误差时的转换矩阵反射器理想指向矢量在卫星坐标系中的矢量P，以及考虑误差时指向矢量在卫星坐标系中的矢量两矢量之差即为天线支撑臂展开锁定后反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差；
3)对多级展开臂中存在的各种误差X按照相关性分为N大类；
4)针对相关的每一大类误差，通过试验数据，得到描述其相关性的多维椭球；
5)利用当前容差设计值，更新椭球；
6)针对每一大类，在由Wi,i＝1,2,...,N确定的高维椭球中进行蒙特卡洛抽样M组，利用上述已经得到的坐标传递矩阵针对每一组样本，得到对应的反射器在卫星坐标系下的滚动、俯仰以及偏航角度误差向量ΔP；
7)对得到的一系列反射器的指向误差值进行统计，从而得到指向误差的最大值；
8)利用上述步骤1)中，计算系统造价；
9)以总造价极小化为目标，以容差为设计变量，以指向误差满足设计要求为约束，进行容差优化设计。
所述步骤2)按如下步骤进行：
步骤201，针对展开臂整体结构中的每一节，建立一个与其固结
的三维直角坐标系；
步骤202，按照空间几何关系，以上述的固结在每一节上的坐标系为基础，建立从展开臂端部反射器到展开臂根部卫星的坐标转换矩阵T和
步骤203，取反射器指向矢量上两个点P1，P2，他们在与反射器固结的转接臂坐标系中的坐标分别为(P1x,P1y,P1z)和(P2x,P2y,P2z)；则P＝T·((P1x,P1y,P1z)-(P2x,P2y,P2z))T；考虑展开误差时，
步骤204，两矢量之差即为天线支撑臂展开锁定后反射器在卫星坐标系下的滚动ΔPx、俯仰ΔPy以及偏航角度偏差ΔPz。
所述的步骤202包括：
a)转接臂上坐标系O4X4Y4Z4，由于反射器固定安装在转接臂上，反射器在伸展臂展开过程中相对于坐标系O4X4Y4Z4坐标位置是固定的；
b)臂间关节坐标系O′3X′3Y′3Z′3，坐标系固定在展开臂上，原点O′3与转接臂坐标系原点O3重合，当臂间机构不存在展开角度误差时，两坐标系Y轴重合，展开角度为θ3；当臂间机构存在展开角度误差时，角度误差为绕X′3偏角为α3,绕Y′3夹角为θ3+β3,绕Z′3偏角为γ3；
c)展开臂1坐标系O3X3Y3Z3也固定在展开臂上，为坐标系O′3X′3Y′3Z′3的平移坐标系，坐标系O′3X′3Y′3Z′3在坐标系O3X3Y3Z3X轴向的距离为t2；
d)回转关节坐标系O'2X'2Y2'Z'2固定在回转关节和根部展开关节的连接部件上，与展开臂坐标系O3X3Y3Z3原点重合，当回转机构不存在展开角度误差时，两坐标系X轴重合，展开角度为θ2；当臂间机构存在展开角度误差时，坐标系O3X3Y3Z3与回转关节坐标系O'2X'2Y′2Z'2绕X'2夹角为θ2+α2,绕Y′2夹角为β2,绕Z'2夹角为γ2；
e)根部展开关节坐标系O2X2Y2Z2与坐标系O'2X'2Y′2Z'2固定在同一部件上，为坐标系O'2X'2Y2'Z'2的平移坐标系，坐标系O'2X'2Y′2Z'2在坐标系O2X2Y2Z2X轴向的距离为t1；
f)根部支座坐标系O1X1Y1Z1固定在根部支座上，与根部展开关节坐标系O2X2Y2Z2原点重合，当根部机构不存在展开角度误差时，两坐标系Z轴重合，展开角度为θ1；当根部机构存在展开角度误差时，坐标系O2X2Y2Z2与坐标系O1X1Y1Z1，绕X1夹角为α1,绕Y1夹角为β1,绕Z1夹角为θ1+γ1；
g)坐标系O0X0Y0Z0为卫星坐标系，通过转换矩阵Q可以将点在坐标系O1X1Y1Z1下的坐标转换为在卫星坐标系下的坐标；
h)得到T和如下，
所述步骤3)按如下过程进行：
步骤301，将展开臂关节分为回转关节和展开关节两大类；
步骤302，考虑到回转关节和展开关节的工作形式的不同，回转关节铰链内部三个方向的误差互相相关；同时，针对同一制造商、同类驱动器、同种工作环境下的不同展开关节的铰链，其各个方向的误差互相相关，两大类之间互相没有相关性。
所述步骤4)按如下过程进行：
步骤401，将展开臂展开n次，分别测量各个关节在各个方向上的误差；
步骤402，在每一大类内部，取出一对误差变量Xi,Xj,i≠j，将对应的二维误差离散点绘制在平面坐标系内，通过最小体积法获得将所有离散点包含在内的最小椭圆，获取每一位变量分别取值区间XI＝[XL,XR]以及椭圆长轴与坐标轴的夹角θ；
步骤403，由椭圆方程获得这两个误差的相关矩阵
其中，
步骤404，回到步骤402，继续抽取Xi,Xj,i≠j，重复步骤403；直到每一大类的任意两组误差变量均被抽取过，继而得到，相关性矩阵Wij,k,k＝1,2,...,N；
步骤405，在每一大类误差中，通过对二维Wij,k的组集，得到每一大类的相关性矩阵Wk,k＝1,2,...,N
步骤406最终，确定每一大类误差的椭球方程为
XTWkX≤1。
所述步骤5)按如下过程进行：
步骤501，在容差优化过程中，不同误差的容差值将会被调整为新的XI＝[XL,XR]；
步骤502，利用θ值和不同维变量的取值上下界，利用以及步骤4，得到更新后的Wij。
所述步骤6)按如下过程进行：
步骤601，在步骤4)所构造的多维椭球的基础上，构造其外接超长方体；
步骤602，在其外接长方体内，假设每一维均为独立的服从均匀分布的随机变量，进行蒙特卡洛抽样，从而构成样本集S；
步骤603，循环遍历样本集S中的每一个样本，判断其是否属于(3f)所定义椭球内，所有落在椭球内的样本点构成样本集S'；
步骤604，将S'各组样本值带入步骤2)，得到反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差。
所述步骤9)按如下过程进行：
步骤901，建立如下可展开天线展开臂关节容差设计优化模型，
FindΔ
Min Cost(Δ)
s.t.max{ΔPx,ΔPy,ΔPz}-ΔP*≤0
式中，Δ为容差向量，ΔPx,ΔPy,ΔPz分别是步骤2)得到的反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差；ΔP*设计人员给定的角度偏差的上限值；
步骤902，利用序列二次规划法来求解上述非线性规划问题，其中，目标函数相对于设计变量的敏度可以由解析表达式求导获得；约束函数的敏度有有限差分法获得；需要判断得到的反射器指向误差是否满足要求，如果满足要求则容差设计合理；否则，修改设计变量值，并重复步骤5)到8)，利用数学规划算法直到得到满足设计要求的容差值。
本发明与现有技术相比，具有如下优点：
1、本方法利用相关性椭球模型，从而能够表征不同方向角度误差之间的相关性，更贴合工程实际。
2、利用优化设计的思路，以系统造价最小化为目标，以各个关节不同角度误差的容差为设计变量，实现了展开臂的容差设计。
3、通过相关性模型以及容差设计，能够有效降低设计的保守性，从而得到既可靠稳健又经济的设计方案。
4、针对某型给定的展开臂形式，利用相关试验，得到任意两维展开臂关节误差的相关椭圆参数，在此基础上，构造得到多维误差的相关椭球。在此相关椭球内进行蒙特卡洛抽样，针对每一组样本，进行天线反射器指向误差的分析。最终得到在椭球模型下的指向误差分布情况。最终，以上述分析流程为基础，以展开臂关节的容差为设计变量，以极小化天线的造价为目标函数，天线指向误差为约束函数，建立容差优化模型，对可展开天线的展开臂关节误差的容差进行设计。
下面结合实施例附图对本发明作进一步说明：
附图说明
图1是本发明考虑误差相关性的可展开天线展开臂容差设计方法流程图；
图2是某类型展开臂结构示意图；
图3是由二维离散试验点构造二维椭球示意图；
图4是二维椭球倾角θ示意图；
图5是卫星-天线整体结构示意图；
图6是仿真实验中采用的展开臂结构示意图。
具体实施方式
如图1所示，基于椭球模型的可展开天线展开臂容差设计方法，其特征是至少包括如下步骤：
1.获取所涉及展开臂类型的容差-成本关系Cost＝C(Δ)，其中，Δ＝|ΔXmax|为各种角度误差的容差；Cost为成本；C为关系函数。
2.建立从展开臂根部到端部的不包含误差的坐标转换矩阵T，得到各个关节的展开角度误差时的转换矩阵反射器理想指向矢量在卫星坐标系中的矢量P，以及考虑误差时指向矢量在卫星坐标系中的矢量两矢量之差即为天线支撑臂展开锁定后反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差；
3.对多级展开臂中存在的各种误差X按照相关性分为N大类；
4.针对相关的每一大类误差，通过试验数据，得到描述其相关性的多维椭球；
5.利用当前容差设计值，更新更新相关椭球的相关性矩阵Wk,k＝1,2,...,N；
6.针对每一大类，在由Wi,i＝1,2,...,N确定的高维椭球中进行蒙特卡洛抽样M组，利用上述已经得到的坐标传递矩阵针对每一组样本，得到对应的反射器在卫星坐标系下的滚动、俯仰以及偏航角度误差向量ΔP；
7.对得到的一系列反射器的指向误差值进行统计，从而得到指向误差的最大值；
8.利用上述步骤(1)中，计算系统造价；
9.以总造价极小化为目标，以容差为设计变量，以指向误差满足设计要求为约束，进行容差优化设计。
如图2所示，所述步骤2按如下过程进行
步骤201，针对展开臂整体结构中的每一节，建立一个与其固结的三维直角坐标系；
步骤202，按照空间几何关系，以上述的固结在每一节上的坐标系为基础，建立从展开臂端部反射器到展开臂根部卫星的坐标转换矩阵T和
如图2所示，所述的步骤202包括：
a)转接臂上坐标系O4X4Y4Z4，由于反射器固定安装在转接臂上，反射器在伸展臂展开过程中相对于坐标系O4X4Y4Z4坐标位置是固定的；
b)臂间关节坐标系O′3X′3Y′3Z′3，坐标系固定在展开臂上，原点O′3与转接臂坐标系原点O3重合，当臂间机构不存在展开角度误差时，两坐标系Y轴重合，展开角度为θ3；当臂间机构存在展开角度误差时，角度误差为绕X′3偏角为α3,绕Y′3夹角为θ3+β3,绕Z′3偏角为γ3；
c)展开臂1坐标系O3X3Y3Z3也固定在展开臂上，为坐标系O′3X′3Y′3Z′3的平移坐标系，坐标系O′3X′3Y′3Z′3在坐标系O3X3Y3Z3X轴向的距离为t2；
d)回转关节坐标系O'2X'2Y′2Z'2固定在回转关节和根部展开关节的连接部件上，与展开臂坐标系O3X3Y3Z3原点重合，当回转机构不存在展开角度误差时，两坐标系X轴重合，展开角度为θ2；当臂间机构存在展开角度误差时，坐标系O3X3Y3Z3与回转关节坐标系O'2X'2Y′2Z'2绕X'2夹角为θ2+α2,绕Y2'夹角为β2,绕Z'2夹角为γ2；
e)根部展开关节坐标系O2X2Y2Z2与坐标系O'2X'2Y′2Z'2固定在同一部件上，为坐标系O'2X'2Y′2Z'2的平移坐标系，坐标系O'2X'2Y′2Z'2在坐标系O2X2Y2Z2X轴向的距离为t1；
f)根部支座坐标系O1X1Y1Z1固定在根部支座上，与根部展开关节坐标系O2X2Y2Z2原点重合，当根部机构不存在展开角度误差时，两坐标系Z轴重合，展开角度为θ1；当根部机构存在展开角度误差时，坐标系O2X2Y2Z2与坐标系O1X1Y1Z1，绕X1夹角为α1,绕Y1夹角为β1,绕Z1夹角为θ1+γ1；
g)坐标系O0X0Y0Z0为卫星坐标系，通过转换矩阵Q可以将点在坐标系O1X1Y1Z1下的坐标转换为在卫星坐标系下的坐标；
h)得到T和如下，
步骤203，取反射器指向矢量上两个点P1，P2，他们在与反射器固结的转接臂坐标系中的坐标分别为(P1x,P1y,P1z)和(P2x,P2y,P2z)。则P＝T·((P1x,P1y,P1z)-(P2x,P2y,P2z))T。考虑展开误差时，
步骤204，两矢量之差即为天线支撑臂展开锁定后反射器在卫星坐标系下的滚动ΔPx、俯仰ΔPy以及偏航角度偏差ΔPz。
所述步骤3按如下过程进行：
步骤301，如图2所示，将展开臂关节分为回转关节和展开关节两大类；
步骤302，考虑到回转关节和展开关节的工作形式的不同，回转关节铰链内部三个方向的误差互相相关；同时，针对同一制造商、同类驱动器、同种工作环境下的不同展开关节的铰链，其各个方向的误差互相相关。两大类之间互相没有相关性。
所述步骤4按如下过程进行：
步骤401，将展开臂展开n次，分别测量各个关节在各个方向上的误差；
步骤402，在每一大类内部，取出一对误差变量Xi,Xj,i≠j，将对应的二维误差离散点绘制在平面坐标系内，通过最小体积法获得将所有离散点包含在内的最小椭圆，获取每一位变量分别取值区间XI＝[XL,XR]以及椭圆长轴与坐标轴的夹角θ；其中，XL,XR分别为XI的下界和上界。如图3，4所示。
步骤403，由椭圆方程获得这两个误差的相关矩阵
其中，
步骤404，回到步骤402，继续抽取Xi,Xj,i≠j，重复步骤403。直到每一大类的任意两组误差变量均被抽取过。继而得到，相关性矩阵Wij,k,k＝1,2,...,N。
步骤405，在每一大类误差中，通过对二维Wij,k的组集，得到每一大类的相关性矩阵Wk,k＝1,2,...,N
步骤406，最终，确定每一大类误差的椭球方程为
XTWkX≤1。
所述步骤5)按如下过程进行：
步骤501，在容差优化过程中，不同误差的容差值将会被调整。即为新的XI＝[XL,XR]，考虑到同一厂家，同类元件的类似性，我们认为θ值保持不变；
步骤502，利用θ值和不同维变量的取值上下界，利用以及步骤4，可以得到更新后的Wij。
所述按如下过程进行：
步骤601，在步骤4所构造的多维椭球的基础上，构造其外接超长方体；
步骤602，在其外接长方体内，假设每一维均为独立的服从均匀分布的随机变量，进行蒙特卡洛抽样，从而构成样本集S；
步骤603，循环遍历样本集S中的每一个样本，判断其是否属于(3f)所定义椭球内。所有落在椭球内的样本点构成样本集S'；
步骤604，将S'各组样本值带入步骤2，得到反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差。
所述步骤9)按如下过程进行：
步骤901，建立如下可展开天线展开臂关节容差设计优化模型，
FindΔ
Min Cost(Δ)
s.t.max{ΔPx,ΔPy,ΔPz}-ΔP*≤0
式中，Δ为容差向量，ΔPx,ΔPy,ΔPz分别是步骤2)得到的反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差；ΔP*设计人员给定的角度偏差的上限值。
步骤902，利用序列二次规划法来求解上述非线性规划问题，其中，目标函数相对于设计变量的敏度可以由解析表达式求导获得；约束函数的敏度有有限差分法获得；需要判断得到的反射器指向误差是否满足要求，如果满足要求则容差设计合理；否则，修改设计变量值，并重复步骤5)到8)，利用数学规划算法直到得到满足设计要求的容差值。
本发明的优点可通过以下仿真实验进一步说明：
1.仿真条件
将本发明的考虑误差相关性的可展开天线展开臂容差设计方法应用于某卫星天线展开臂关节容差设计上。天线的整体结构如图5所示。图6为展开臂的详细结构示意，此结构包括1个回转关节，4个展开关节。各关节初始给定的容差如表1所示。设计要求为反射器各个方向角度误差均小于0.06°。
表1展开臂初始给定关节容差
X方向(°) Y方向(°) Z方向(°)
回转关节 ±0.02 ±0.02 ±0.025
展开关节 ±0.03 ±0.03 ±0.03
2.仿真结果
1)天线反射器指向误差分析
表2由展开臂角度误差引起的天线指向误差
X方向误差(°) Y方向误差(°) Z方向误差(°)
传统分析方法 0.138 0.139 0.166
新分析方法 0.077 0.077 0.045
在初始给定的容差基础上，分别利用传统分析方法和本发明所提出的分析方法进行天线指向误差的分析。由表2可以看出，因为传统分析方法忽视了各个展开角度误差之间的相关性，所以导致的天线三个角度的误差远远大于设计指标要求。而本发明提出的新方法，可以有效降低天线的指向误差。
2)容差设计
考虑到表2中，在初始容差的基础上，天线指向误差依然不能满足0.06°。利用本发明的容差设计方法，以回转关节和展开关节展开误差的容差为设计变量，最优结果对比如下，
表3容差设计的最优值以及对应的反射器指向误差值
可见，通过容差设计，当回转关节和展开关节的容差如上表所示时，天线三个方向的误差能满足设计要求。上述仿真数值试验证明，采用本发明可有效。

Claims (8)

1.基于椭球模型的可展开天线展开臂容差设计方法，其特征是至少包括如下步骤：
1)获取所涉及展开臂类型的容差-成本关系Cost＝C(Δ)，其中，Δ＝|ΔXmax|为各种角度误差的容差；Cost为成本；C为关系函数；
2)建立从展开臂根部到端部的不包含误差的坐标转换矩阵T，得到各个关节的展开角度误差时的转换矩阵反射器理想指向矢量在卫星坐标系中的矢量P，以及考虑误差时指向矢量在卫星坐标系中的矢量两矢量之差即为天线支撑臂展开锁定后反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差；
3)对多级展开臂中存在的各种误差X按照相关性分为N大类；
4)针对相关的每一大类误差，通过试验数据，得到描述其相关性的多维椭球；
5)利用当前容差设计值，更新椭球；
6)针对每一大类，在由Wi,i＝1,2,...,N确定的高维椭球中进行蒙特卡洛抽样M组，利用上述已经得到的坐标传递矩阵针对每一组样本，得到对应的反射器在卫星坐标系下的滚动、俯仰以及偏航角度误差向量ΔP；
7)对得到的一系列反射器的指向误差值进行统计，从而得到指向误差的最大值；
8)利用上述步骤1)中，计算系统造价；
9)以总造价极小化为目标，以容差为设计变量，以指向误差满足设计要求为约束，进行容差优化设计。
2.根据权利要求1所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述步骤2)按如下步骤进行：
步骤201，针对展开臂整体结构中的每一节，建立一个与其固结的三维直角坐标系；
步骤202，按照空间几何关系，以上述的固结在每一节上的坐标系为基础，建立从展开臂端部反射器到展开臂根部卫星的坐标转换矩阵T和
步骤203，取反射器指向矢量上两个点P1，P2，他们在与反射器固结的转接臂坐标系中的坐标分别为(P1x,P1y,P1z)和(P2x,P2y,P2z)；则P＝T·((P1x,P1y,P1z)-(P2x,P2y,P2z))T；考虑展开误差时，
步骤204，两矢量之差即为天线支撑臂展开锁定后反射器在卫星坐标系下的滚动ΔPx、俯仰ΔPy以及偏航角度偏差ΔPz。
3.根据权利要求2所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述的步骤202包括：
a)转接臂上坐标系O4X4Y4Z4，由于反射器固定安装在转接臂上，反射器在伸展臂展开过程中相对于坐标系O4X4Y4Z4坐标位置是固定的；
b)臂间关节坐标系O3'X3'Y3'Z3'，坐标系固定在展开臂上，原点O3'与转接臂坐标系原点O3重合，当臂间机构不存在展开角度误差时，两坐标系Y轴重合，展开角度为θ3；当臂间机构存在展开角度误差时，角度误差为绕X3'偏角为α3,绕Y3'夹角为θ3+β3,绕Z3'偏角为γ3；
c)展开臂1坐标系O3X3Y3Z3也固定在展开臂上，为坐标系O3'X3'Y3'Z3'的平移坐标系，坐标系O3'X3'Y3'Z3'在坐标系O3X3Y3Z3X轴向的距离为t2；
d)回转关节坐标系O'2X'2Y2'Z'2固定在回转关节和根部展开关节的连接部件上，与展开臂坐标系O3X3Y3Z3原点重合，当回转机构不存在展开角度误差时，两坐标系X轴重合，展开角度为θ2；当臂间机构存在展开角度误差时，坐标系O3X3Y3Z3与回转关节坐标系O'2X'2Y2'Z'2绕X'2夹角为θ2+α2,绕Y2'夹角为β2,绕Z'2夹角为γ2；
e)根部展开关节坐标系O2X2Y2Z2与坐标系O'2X'2Y2'Z'2固定在同一部件上，为坐标系O'2X'2Y2'Z'2的平移坐标系，坐标系O'2X'2Y2'Z'2在坐标系O2X2Y2Z2X轴向的距离为t1；
f)根部支座坐标系O1X1Y1Z1固定在根部支座上，与根部展开关节坐标系O2X2Y2Z2原点重合，当根部机构不存在展开角度误差时，两坐标系Z轴重合，展开角度为θ1；当根部机构存在展开角度误差时，坐标系O2X2Y2Z2与坐标系O1X1Y1Z1，绕X1夹角为α1,绕Y1夹角为β1,绕Z1夹角为θ1+γ1；
g)坐标系O0X0Y0Z0为卫星坐标系，通过转换矩阵Q可以将点在坐标系O1X1Y1Z1下的坐标转换为在卫星坐标系下的坐标；
h)得到T和如下，
<mrow> <mi>T</mi> <mo>=</mo> <mi>Q</mi> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>t</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mi>t</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <mover> <mi>T</mi> <mo>^</mo> </mover> <mo>=</mo> <mi>Q</mi> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mi>t</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>t</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <msub> <mi>t</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>t</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>sin</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>.</mo> </mrow>
4.根据权利要求1所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述步骤3)按如下过程进行：
步骤301，将展开臂关节分为回转关节和展开关节两大类；
步骤302，考虑到回转关节和展开关节的工作形式的不同，回转关节铰链内部三个方向的误差互相相关；同时，针对同一制造商、同类驱动器、同种工作环境下的不同展开关节的铰链，其各个方向的误差互相相关，两大类之间互相没有相关性。
5.根据权利要求1所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述步骤4)按如下过程进行：
步骤401，将展开臂展开n次，分别测量各个关节在各个方向上的误差；
步骤402，在每一大类内部，取出一对误差变量Xi,Xj,i≠j，将对应的二维误差离散点绘制在平面坐标系内，通过最小体积法获得将所有离散点包含在内的最小椭圆，获取每一位变量分别取值区间XI＝[XL,XR]以及椭圆长轴与坐标轴的夹角θ；
步骤403，由椭圆方程获得这两个误差的相关矩阵
<mrow> <msub> <mi>W</mi> <mrow> <mi>i</mi> <mi>j</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>X</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>j</mi> </msub> <mo>,</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>,</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> </mrow>
其中， D(Xi)≠D(Xj)；
步骤404，回到步骤402，继续抽取Xi,Xj,i≠j，重复步骤403；直到每一大类的任意两组误差变量均被抽取过，继而得到，相关性矩阵Wij,k,k＝1,2,...,N；
步骤405，在每一大类误差中，通过对二维Wij,k的组集，得到每一大类的相关性矩阵Wk,k＝1,2,...,N
<mrow> <msub> <mi>W</mi> <mi>k</mi> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>X</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>n</mi> </msub> <mo>,</mo> <msub> <mi>X</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>n</mi> </msub> <mo>,</mo> <msub> <mi>X</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow>
步骤406，最终，确定每一大类误差的椭球方程为
XTWkX≤1。
6.根据权利要求1所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述步骤5)按如下过程进行：
步骤501，在容差优化过程中，不同误差的容差值将会被调整为新的XI＝[XL,XR]；
步骤502，利用θ值和不同维变量的取值上下界，利用以及步骤4，得到更新后的Wij。
7.根据权利要求1所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述步骤6)按如下过程进行：
步骤601，在步骤4)所构造的多维椭球的基础上，构造其外接超长方体；
步骤602，在其外接长方体内，假设每一维均为独立的服从均匀分布的随机变量，进行蒙特卡洛抽样，从而构成样本集S；
步骤603，循环遍历样本集S中的每一个样本，判断其是否属于(3f)所定义椭球内，所有落在椭球内的样本点构成样本集S'；
步骤604，将S'各组样本值带入步骤2)，得到反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差。
8.根据权利要求1所述的基于椭球模型的可展开天线展开臂容差设计方法，其特征是，所述步骤9)按如下过程进行：
步骤901，建立如下可展开天线展开臂关节容差设计优化模型，FindΔ
Min Cost(Δ)
s.t.max{ΔPx,ΔPy,ΔPz}-ΔP*≤0
式中，Δ为容差向量，ΔPx,ΔPy,ΔPz分别是步骤2)得到的反射器在卫星坐标系下的滚动、俯仰以及偏航角度偏差；ΔP*设计人员给定的角度偏差的上限值；
步骤902，利用序列二次规划法来求解上述非线性规划问题，其中，目标函数相对于设计变量的敏度可以由解析表达式求导获得；约束函数的敏度有有限差分法获得；需要判断得到的反射器指向误差是否满足要求，如果满足要求则容差设计合理；否则，修改设计变量值，并重复步骤5)到8)，利用数学规划算法直到得到满足设计要求的容差值。
CN201711036984.9A 2017-07-21 2017-10-30 基于椭球模型的可展开天线展开臂容差设计方法 Active CN107863599B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201710599557 2017-07-21
CN2017105995575 2017-07-21

Publications (2)

Publication Number Publication Date
CN107863599A true CN107863599A (zh) 2018-03-30
CN107863599B CN107863599B (zh) 2019-11-22

Family

ID=61697556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711036984.9A Active CN107863599B (zh) 2017-07-21 2017-10-30 基于椭球模型的可展开天线展开臂容差设计方法

Country Status (1)

Country Link
CN (1) CN107863599B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113110548A (zh) * 2021-04-21 2021-07-13 北京控制工程研究所 一种航天器椭球集合演化的设计方法
CN114139299A (zh) * 2021-10-25 2022-03-04 西北工业大学 降低展开臂关节间隙对网状天线电性能影响的设计方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040017332A1 (en) * 2002-07-26 2004-01-29 Rao Sudhakar K. Multiple beam antenna using reflective and partially reflective surfaces
CN102073754A (zh) * 2010-09-20 2011-05-25 西安电子科技大学 基于误差因素的反射面天线机电综合分析方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040017332A1 (en) * 2002-07-26 2004-01-29 Rao Sudhakar K. Multiple beam antenna using reflective and partially reflective surfaces
CN102073754A (zh) * 2010-09-20 2011-05-25 西安电子科技大学 基于误差因素的反射面天线机电综合分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡乃岗等: "大型相控阵天线结构与调整机构一体化设计", 《机械工程学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113110548A (zh) * 2021-04-21 2021-07-13 北京控制工程研究所 一种航天器椭球集合演化的设计方法
CN114139299A (zh) * 2021-10-25 2022-03-04 西北工业大学 降低展开臂关节间隙对网状天线电性能影响的设计方法
CN114139299B (zh) * 2021-10-25 2022-09-23 西北工业大学 降低展开臂关节间隙对网状天线电性能影响的设计方法

Also Published As

Publication number Publication date
CN107863599B (zh) 2019-11-22

Similar Documents

Publication Publication Date Title
Cavagna et al. NeoCASS: an integrated tool for structural sizing, aeroelastic analysis and MDO at conceptual design level
CN103991555B (zh) 一种用于飞机数字化装配的自动化测量方法
CN109470268B (zh) 一种提高卫星姿态确定精度和效率的方法
Aftosmis et al. Adjoint-based low-boom design with Cart3D
CN104121992A (zh) 空间目标红外辐射全任务空域动态特性分析方法
CN104899365A (zh) 一种可减小气动干扰不利影响的直升机气动布局优化方法
CN110059422B (zh) 频率选择表面曲面天线罩的电磁散射特性仿真方法
CN103106305A (zh) 一种基于实测模态的空间网格结构模型分步修正方法
Wei et al. Drogue dynamic model under bow wave in probe-and-drogue refueling
CN107863599B (zh) 基于椭球模型的可展开天线展开臂容差设计方法
CN104156523A (zh) 索网反射面天线的测地线索网生成方法
CN106971043B (zh) 一种嵌套网格的建立方法及装置
CN105843074A (zh) 一种获取天线在轨振动影响的动力学建模方法
Tonelli et al. Stability of statics aware voronoi grid-shells
CN107491594B (zh) 一种索网天线形面精度计算方法
Dalle et al. Adjustments and Uncertainty Quantification for SLS Aerodynamic Sectional Loads
Xue et al. Structural optimization of a nose landing gear considering its fatigue life
CN107344630A (zh) 一种多视场星敏感器星上构形布局确定方法
CN112347561A (zh) 飞行器的静气动弹性分析方法、装置、设备和存储介质
Ramadoss et al. Design, construction, and control of curves and surfaces via deployable mechanisms
Brezillon et al. Aerodynamic inverse design framework using discrete adjoint method
Peng et al. Overset structured grids assembly method for numerical simulations of multi-bodies and moving objects
Lora-Clavijo et al. PBH mass growth through radial accretion during the radiation dominated era
CN117094071A (zh) 一种卫星sar天线动力学快速建模方法及系统
Sitaraman et al. Prediction of helicopter maneuver loads using a fluid-structure analysis

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