基于成对推力器连续工作的卫星质心位置测量方法技术

本发明专利技术提供一种基于成对推力器连续工作的卫星质心位置测量方法，对推力器安装布局进行设计；选择成对推力器进行连续喷气工作；地面处理遥测下传的陀螺测量数据，对卫星的质心位置进行计算。本发明专利技术通过选定推力器进行固定时间长度连续工作，可避免推力器的输出特性变化带来的测量误差；在数据处理中通过对推力器产生的力矩取比值的方式，避免了推力大小偏差引起的测量误差，提高了质心测量精度。

【技术实现步骤摘要】

本专利技术涉及卫星姿态与轨道控制技术，尤其是对于储箱平铺安装结构的地球静止轨道卫星在转移轨道段期间的姿轨控技术，利用本专利技术对卫星的质心位置进行测量，为储箱的均衡排放、混合比调节提供参考。
技术介绍
受运载能力所限，地球静止轨道卫星需要星上的变轨发动机进行远地点变轨，将卫星由大椭圆轨道送入地球静止轨道。变轨过程中随着燃料消耗，由于储箱的不均衡排放、混合比偏差等因素的影响，卫星的质心位置会发生较大变化，可能会导致发动机推力偏心产生较大的干扰，增加姿控消耗的燃料。由于不均衡排放、混合比偏差等不确定因素，地面难以准确计算出卫星的质心位置。为减小变轨时的姿态扰动力矩，实现对卫星质心位置在轨调节，有必要对卫星的质心位置进行在轨测量。当前大多数方案均基于推力器随机离散喷气的方式，推力器工作在小脉冲宽度条件下，利用递推算法对质心位置进行在轨计算，星上处理数据量大，计算复杂；另外，在小脉冲宽度条件下，推力器实际产生的推力大小变化较大，并且与标称值之间存在较大偏差，导致质心测量误差大。
技术实现思路
针对现有技术存在的缺点，本专利技术提供一种基于成对推力器连续工作的方式，利用陀螺测量数据实现对卫星质心位置测量。该专利技术通过选择相应的推力器在固定时间内连续工作，将陀螺测量数据遥测下传，结合推力器的安装位置及角度，通过对陀螺测量数据进行处理计算得到卫星的质心位置。为了达到上述目的，本专利技术的技术方案是提供一种基于成对推力器连续工作的卫星质心位置测量方法：布置至少一组推力器，每一组中的两个推力器产生相反控制力矩；选取被布置的其中任意一组推力器，控制这一组的两个推力器在固定时间内同时连续喷气；将这一组的两个推力器同时连续喷气期间的陀螺测量数据遥测下传进行处理，得到卫星三轴惯性角速度；根据卫星三轴惯性角速度与卫星转动惯量，基于姿态动力学方程计算这一组中两个推力器在同时连续喷气期间各自产生的力矩及力矩的比值；被布置的其中一组推力器是第一组推力器(2A、3A)，所述第一组推力器(2A、3A)安装位置的y向坐标，相对于卫星y向平均质心位置对称；计算第一组推力器(2A、3A)同时连续喷气产生的俯仰力矩与偏航力矩，根据俯仰力矩与偏航力矩的比值以及第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标，计算卫星x向质心位置。优选地，被布置的其中一组推力器是第二组推力器(4A、5A)，所述第二组推力器(4A、5A)安装位置的x向坐标，相对于卫星x向平均质心位置对称；计算第二组推力器(4A、5A)同时连续喷气产生的滚动力矩与偏航力矩，根据滚动力矩与偏航力矩的比值以及第二组推力器(4A、5A)在卫星布局坐标系下的y向位置坐标，计算卫星y向质心位置。优选地，被布置的其中一组推力器是第三组推力器(6A、7A)，所述第三组推力器(6A、7A)安装位置的z向坐标，相对于卫星z向平均质心位置对称；计算第三组推力器(6A、7A)同时连续喷气产生的滚动力矩与俯仰力矩，根据滚动力矩与俯仰力矩的比值以及第三组推力器(6A、7A)在卫星布局坐标系下的z向位置坐标，计算卫星z向质心位置。优选地，每一组中的各个推力器均倾斜安装，保证每个推力器产生的推力与卫星本体的轴不平行。优选地，第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标，第二组推力器(4A、5A)在卫星布局坐标系下的y向位置坐标，第三组推力器(6A、7A)在卫星布局坐标系下的z向位置坐标，是在地面安装相应推力器时分别通过测量得到。本专利技术提出了一种基于成对推力器连续工作的卫星质心位置测量方法。依据此方法可实现对卫星的质心位置进行测量。本专利技术中对推力器安装布局进行设计；选择推力器进行连续喷气工作；地面处理遥测下传的陀螺测量数据，对卫星的质心位置进行计算。本专利技术采用的方法与现有技术相比，其优点在于：本专利技术采用产生相反控制力矩的两个推力器同时连续喷气的方式，保证推力器产生的推力近似为恒值，从而使得喷气过程中卫星的角加速度近似为常值，避免了推力器的推力大小变化所引起的测量误差；此外，通过对力矩取比值的方式，可消除推力器稳态推力大小偏差带来的误差。附图说明图1为星上姿控推力器的安装布局示意图；图2为质心测量过程流程图。具体实施方式本专利技术提供一种基于成对推力器连续工作的卫星质心位置测量方法，以下将结合附图和实施例对本专利技术作进一步说明。1.推力器安装布局为了保证能够利用该专利技术进行质心位置计算，对推力器的安装布局要求如下：(1)推力器倾斜安装，保证每个推力器产生的推力与卫星本体轴不平行；(2)滚动推力器的推力与星体的YOZ平面相互平行；俯仰推力器的推力与星体的XOZ平面相互平行；偏航推力器的推力与星体的XOY平面相互平行；(3)产生相反控制力矩的两个推力器按照平均质心位置对称安装。本专利技术中典型的推力器安装布局如图1所示，推力器2A、3A安装位置的y向坐标相对于卫星y向平均质心位置对称；推力器4A、5A安装位置的x向坐标相对于卫星x向平均质心位置对称；推力器6A、7A安装位置的z向坐标相对于卫星z向平均质心位置对称。2.选择推力器及工作时间长度在质心测量过程中，选择产生方向相反控制力矩的两个推力器在固定时间内连续喷气的方式，可保证测量过程中推力器产生的推力近似为常值，同时可保证卫星的角速度不会过大；此外，延长推力器的工作时间，可以获取足够多的测量信息，减小陀螺测量噪声对测量精度的影响。3.对遥测下传的陀螺测量的角速度信息进行处理将推力器喷气期间陀螺测量得到的三轴惯性角速度遥测下传，通过拟合的方式可估算出推力器喷气引起的三轴惯性角加速度；由姿态动力学方程可计算出推力器产生的作用力矩然后根据推力器的安装位置及角度(地面安装时可测量得到)，由力矩公式可计算出卫星的质心位置。按照图1中推力器的安装布局图，选择成对工作的推力器2A与3A、4A与5A、6A与7A分别“喷气”，实现对卫星X、Y、Z向质心位置的测量。如图2所示，本专利技术的方法包含以下步骤：①t0时刻注数使系统进入到质心测量模式后，星上计算机停止发送控制指令，通过地面注数控制推力器2A、3A在t0时刻开始同时工作，喷气时间长度为10s(暂定，可进行修改)，将喷气期间卫星三轴惯性角速度遥测下传。②在t0+100s时刻，地面注数使系统进入惯性保持模式。③当卫星三轴角速度均小于0.01°/s时，通过地面注数使系统再次进入到质心测量模式。④t1时刻注数使系统进入到质心测量模式后，通过地面注数控制推力器4A、5A在t1时刻开始同时工作，喷气时间长度为10s，将喷气期间卫星三轴惯性角速度遥测下传。⑤在t1+100s时刻，地面注数使系统进入惯性保持模式。⑥t2时刻注数使系统进入到质心测量模式后，通过地面注数控制推力器6A、7A在t2时刻开始同时工作，喷气时间长度为10s，将喷气期间卫星三轴惯性角速度遥测下传。⑦在t2+100s时刻，地面注数使系统退出质心测量模式。⑧对遥测下传的陀螺测量数据进行处理，计算卫星质心位置。对上述①、④、⑥中遥测下传的数据进行处理的算法如下：假设每个控制周期内的惯性角加速度不变，将遥测下传的陀螺测量数据按照形式进行拟合，其中,将为卫星三轴惯性角速度，由陀螺组合测量遥测下传得到，t为遥测下传的时间，t＝t0+k*T(k＝1,2,…),t0为第一拍数据对应的时间本文档来自技高网...

【技术保护点】
一种基于成对推力器连续工作的卫星质心位置测量方法，其特征在于，布置至少一组推力器，每一组中的两个推力器产生相反控制力矩；选取被布置的其中任意一组推力器，控制这一组的两个推力器在固定时间内同时连续喷气；将这一组的两个推力器同时连续喷气期间的陀螺测量数据遥测下传进行处理，得到卫星三轴惯性角速度；根据卫星三轴惯性角速度与卫星转动惯量，基于姿态动力学方程计算这一组中两个推力器在同时连续喷气期间各自产生的力矩及力矩的比值；被布置的其中一组推力器是第一组推力器(2A、3A)，所述第一组推力器(2A、3A)安装位置的y向坐标，相对于卫星y向平均质心位置对称；计算第一组推力器(2A、3A)同时连续喷气产生的俯仰力矩与偏航力矩，根据俯仰力矩与偏航力矩的比值以及第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标，计算卫星x向质心位置。

【技术特征摘要】
1.一种基于成对推力器连续工作的卫星质心位置测量方法，其特征在于，布置至少一组推力器，每一组中的两个推力器产生相反控制力矩；选取被布置的其中任意一组推力器，控制这一组的两个推力器在固定时间内同时连续喷气；将这一组的两个推力器同时连续喷气期间的陀螺测量数据遥测下传进行处理，得到卫星三轴惯性角速度；根据卫星三轴惯性角速度与卫星转动惯量，基于姿态动力学方程计算这一组中两个推力器在同时连续喷气期间各自产生的力矩及力矩的比值；被布置的其中一组推力器是第一组推力器(2A、3A)，所述第一组推力器(2A、3A)安装位置的y向坐标，相对于卫星y向平均质心位置对称；计算第一组推力器(2A、3A)同时连续喷气产生的俯仰力矩与偏航力矩，根据俯仰力矩与偏航力矩的比值以及第一组推力器(2A、3A)在卫星布局坐标系下的x向位置坐标，计算卫星x向质心位置。2.如权利要求1所述的卫星质心位置测量方法，其特征在于，被布置的其中一组推力器是第二组推力器(4A、5A)，所述第二组推力器(4A、5A)安装位置的x向坐标，相对于卫星x向平均质心位置对称；计算第二组推力器(4A...

【专利技术属性】
技术研发人员：俞洁，郭正勇，何益康，赵永德，马雪阳，余维，
申请(专利权)人：上海航天控制技术研究所，
类型：发明
国别省市：上海;31