地心轨道引力波探测无拖曳系统平动控制策略
doi: 10.11887/j.cn.202402004
郝立维1 , 张锦绣2 , 王继河2 , 张谕3 , 孙玥4,5
1. 哈尔滨工业大学 卫星技术研究所,黑龙江 哈尔滨 150001
2. 中山大学 航空航天学院,广东 深圳 518107
3. 中山大学 天琴中心, 广东 珠海 519082
4. 上海市空间智能控制技术重点实验室,上海 201109
5. 上海航天控制技术研究所,上海 201109
基金项目: 广东省基础与应用基础研究重大资助项目(2019B030302001)
Translational tracking strategy of drag-free system for gravitational wave detection in geocentric orbit
HAO Liwei1 , ZHANG Jinxiu2 , WANG Jihe2 , ZHANG Yu3 , SUN Yue4,5
1. Research Center of Satellite Technology, Harbin Institute of Technology, Harbin 150001 , China
2. School of Aeronautics and Astronautics, Sun Yat-Sen University, Shenzhen 518107 , China
3. TianQin Research Center, Sun Yat-Sen University, Zhuhai 519082 , China
4. Shanghai Key laboratory of Aerospace Intelligent Control Technology, Shanghai 201109 , China
5. Shanghai Institute of Spaceflight Control Technology, Shanghai 201109 , China
摘要
探讨了一种针对空间引力波探测任务的在轨无拖曳控制技术,基于未来可行的地心轨道探测任务背景进行分析设计,并对搭载两颗检验质量的在轨无拖曳系统进行航天器与质量块间相对运动动力学及耦合特性建模。同时,初步分析了任务中无拖曳系统指标和摄动,并设计了基于频域H最优控制理论的系统相对平动控制律。数值仿真结果表明,当双检验质量在轨无拖曳系统各检验质量按激光测距呼吸角排列时,采用无固定追踪点策略且在非敏感轴无悬浮控制输入的情况下,可以实现航天器对基准点的追踪,并满足系统频域性能指标的要求。同时,每颗检验质量的时域偏移量可以控制在微米级别,从而获得任务所需的纯引力基准。
Abstract
An on-orbit drag-free control technique for spaceborne gravitational wave detection missions was discussed. Based on the analysis and design of a possible future geocentric orbit detection mission, the relative motion dynamics and coupling characteristics between the spacecraft and mass blocks of an on-orbit drag-free system with with two test masses were modeled. At the same time, the performance index and perturbation of the drag-free system in the mission were preliminarily analyzed, and a relative translational control law based on frequency domain H optimal control theory was designed. Numerical simulation results show that when the test masses of the two-test-mass on-orbit drag-free system are arranged according to the breathing angle of laser rangefinder, without a fixed tracking point strategy and without suspension control input along the non-sensitive axis, the spacecraft can achieve tracking of the reference point while meeting the frequency domain performance index of the system. At the same time, the time domain displacement of each test mass can be controlled to the micron level, thus obtaining the pure gravitational reference required by the mission.
空间引力波探测计划以欧航局的LISA任务[1-2]为代表,该计划将通过三颗无拖曳航天器构成大型星间激光干涉仪,通过每两颗航天器间往返链路的延迟干涉测距,探测中低频段引力波。LISA任务以深空日心轨道为背景,而本文将基于未来可行的地心轨道引力波探测任务进行分析设计,例如天琴计划。天琴计划[3-6]是由中山大学牵头的国际空间引力波探测计划,该任务航天器编队运行在十万千米地心轨道,三颗航天器组成等边三角形构型,在科学模式下,每颗航天器通过无拖曳控制补偿其与引力基准的受摄差,提供可供引力波观测的窗口。
引力波探测航天器需要攻关的关键技术之一是多检验质量的在轨无拖曳控制技术[7],引力波探测编队中每颗航天器内含两个真空腔体,质量块处腔体中,两检验质量敏感轴所成夹角为测距呼吸角,航天器对两颗检验质量提供的惯性基准进行追踪,进而实现测距所需平台的超静超稳。航天器执行机构常用FEEP微推进器[8]、胶体或冷气微推进器[9],质量块执行机构为腔体内静电悬浮控制极板,航天器与质量块间敏感器为电容位移传感器,量测各质量块在腔体内的偏移量。
针对多检验质量在轨无拖曳控制方法与控制策略的研究主要有以下成果。Fichter等以质量块1的三轴位置、质量块2的两条非敏感轴向相对位置、质量块1沿X轴的姿态角六个自由度为无拖曳系统追踪对象[10],由于LISA Pathfinder双检验质量平行放置,使得这种追踪具有可行性。随后该团队在对系统参数量级分析的基础上,解耦LISA Pathfinder的无拖曳动力学模型,整个系统分为无拖曳系统、质量块悬浮系统和卫星姿态系统后,将动力学模型进一步简化为单输入单输出的12个回路[11],并基于H最优控制理论对各轴进行单独的控制器设计,随后将其嵌入非线性动力学进行仿真,验证了这种解耦模型处理方式的合理性和控制器的可行性,保证了无拖曳控制中的性能指标,并分析了系统在参数摄动下的稳定性[12],除此之外该团队对LISA任务基于同样的思路进行设计[13],该文追踪策略被描述为一个选择矩阵的形式,期望于平衡卫星对XOZ面内质量块2和质量块1的XZ轴的追踪;Wu等通过回路解耦的思路对LISA Pathfinder动力学进行处理,与文献[12]采用相同的控制策略,基于定量反馈理论对各轴控制律进行设计与仿真[14],通过对系统性能和鲁棒稳定性的权衡,在达成既定指标同时留有一定的稳定裕度。
当前各类文献针对无拖曳控制的研究多集中于单检验质量情况,对多检验质量无拖曳系统的研究和分析并没有一个成熟的方案,同时缺乏对各类平动控制策略的详解和仿真[10-15],本文也以此为出发点,重点考量双检验质量无拖曳系统追踪策略,特别是,目前的双检验质量研究文献,其背景大多针对LISA日心轨道[10-14],国内外均没有相关文献面向地心轨道多检验质量航天器进行平动控制策略研究分析,与日心轨道相比,地心轨道在非敏感轴向上承受更高的非保守力摄动,同时航天器与检验质量间的引力梯度摄动量级也将比日心轨道情况高,这都是本文在设计中需要注意的新问题。
引力波探测编队航天器姿态通过激光链路进行捕获和校准,在轨道面内达成指向姿态稳定,因此,对空间引力波探测编队航天器,其无拖曳追踪策略以平动自由度的追踪为主要任务,本文将对地心轨道引力波探测任务中单颗航天器内双检验质量无拖曳相对平动系统的动力学、控制算法、追踪策略等问题进行初步的研究和探讨。
1 无拖曳系统相对平动动力学
无拖曳航天器控制过程分为两个阶段,空间探测试验执行之前,将针对质量块进行第一阶段控制,用以补偿质量块初释放误差,这一阶段后,无拖曳科学模式运行,其中第二阶段的质量块控制信号会受到严格的指标约束,在后文会详细叙述。本节所研究追踪策略为后一阶段即无拖曳科学模式阶段控制策略,系统的边界条件为第一阶段控制结果。
无拖曳科学模式中,双质量块无拖曳系统所追踪惯性基准点与两个质量块相关联[1],因此有多种追踪策略,下面对三种典型平动追踪策略进行简述,无拖曳系统示意图见图1,图中TM1、TM2分别表示质量块1和质量块2,CoMsc表示质心,追踪策略示意图见图2
引力波探测航天器中首先考虑以单颗检验质量为基准的追踪方案。这种追踪策略中,航天器追踪其中一颗检验质量的平动,另一颗检验质量追踪其所在腔体运动,悬浮回路作用于后一颗检验质量。
LOS点作为两条激光链路的反向延长线交点,是一个很好确认、衡量、解算且切实有效的追踪点,这一策略中,无拖曳控制期望于航天器相对LOS点位置稳定,其对非敏感轴向的补偿达到最优。
1双检验质量无拖曳系统示意图
Fig.1Diagram of the two-test-mass drag-free system
2三种追踪策略示意图
Fig.2Schematic diagram of three tracking strategies
无固定追踪点策略期望于将问题转化为航天器体坐标系XOZ面内的三点追踪问题,航天器相对于两颗质量块所组成的系统保持稳定,将位移偏差解算到体坐标系下进行补偿。
上述追踪策略均可分为两种情况,即敏感轴向有无悬浮控制信号。引入敏感轴控制信号的优势在于额外的控制补偿会带来更小的时域偏差和更短的修正时间,但敏感轴无控制输入情况可减少检验质量残余噪声来源,简化测距轴向的噪声分析,而这一点正是引力波探测任务关注的重点。本文所探求的是无敏感轴悬浮控制情况下的可行追踪策略。
双检验质量无拖曳系统动力学模型包含航天器、质量块姿态动力学与无拖曳系统相对运动动力学。在无拖曳科学模式中,不采用轨道机动,本文研究中航天器编队每两颗之间组成激光链路,指向姿态角保持稳定,即视航天器轨道坐标系与体坐标系固连,主要探讨航天器与检验质量的平动控制策略。记航天器体坐标系三轴为XYZ,质量块体坐标系三轴为XiYiZii=1,2)。由于航天器XZ轴向的微推控制均会与质量块测距轴向产生耦合(即质量块的敏感轴向),因此在分析中,将两轴均视为航天器的敏感轴,Y轴视为航天器的非敏感轴。
1.1 检验质量与腔体间相对运动动力学模型
质量块、航天器及腔体间平动问题即为惯性空间三个质点间相对运动问题,假定检验质量理想悬浮点为腔体中心,以腔体中心为过渡点来推导航天器与质量块间动力学方程,方程的最终目的应为建立质量块在腔体内偏移量的微分方程。质量块与腔体间相对运动方程为:
ToSr¨OM=2ωS×ToSr˙MO+ω˙S×ToSrMO+rSO+ωS×ωS×TOSrMO+rSO+aMS-aS
(1)
式中:ToS为从腔体到航天器的坐标转换矩阵;ωS表示航天器角速度;rMO表示从质量块到航天器位置矢量;rSO表示航天器到腔体位置矢量;aMSaS分别表示质量块、航天器所受外力造成的加速度总和;(aMS-aS)项则用来描述质量块与航天器间受力差值,该项主要包含质量块与航天器的非保守力与执行机构控制力、质量块与航天器间耦合力以及质量块与航天器分属不同轨道引起的引力梯度。
地心引力波探测轨道背景下[5],航天器所受非保守力以太阳光压为主,记为dSC,航天器微推控制记为uSC,质量块非保守力记为dTM,质量块悬浮机构控制力记为uTM。质量块非保守力主要由四部分组成:行星际空间和内部磁场带来的磁干扰,宇宙射线、残余气体和光子辐射引入的环境摄动,温度波动引起的摄动,质量块与传感器间库仑力[16]
质量块与航天器间引力梯度记为ΔG,是由无拖曳系统航天器质心与质量块质心不重合导致,航天器与质量块分属不同轨道,受到星体引力不同。引力梯度包含中心天体引力梯度ΔGEarth,非球形摄动加速度差值ΔGJ2,其他星体造成引力梯度ΔGO
ΔG=ΔGEarth +ΔGJ2+ΔG0ΔGEarth =-GMEarth rTM3rTM+GMEarth rSC3rSCΔGJ2=r¨J2TM-r¨J2SCΔG0=GMOrmtrmt3-rmt-rmermt-rme3-rstrst3-rst-rserst-rse3
(2)
其中,MEarth表示地球质量,MO表示对应星体的质量, rTMrSC分别表示质量块、航天器相对地球的位置, rmtrme分别表示扰动星体分析中质量块到扰动星体和地球的位置矢量,rstrse分别表示扰动星体分析中航天器到扰动星体和地球的位置矢量,r¨J2TMr¨J2SC分别表示J2项对质量块与航天器的影响,G为万有引力常数。引力梯度由航天器与质量块间质心初始间距d和质量块在腔体内偏移量||rMO||决定,但由于||rMO||作为被控项,量级远低于图2d量级(0.1 m级),因此控制律设计将此项视为扰动项,并基于此进行模型线性化。
航天器与质量块间的相对运动耦合项主要包含航天器与质量块间的自引力、质量块与腔体间的电荷效应、磁效应等。这些摄动的作用随质量块在腔体中的偏移量发生改变,在无拖曳系统中将这一部分力统称为寄生刚度耦合[17],这一改变的系数称为刚度系数。定义平动自由度刚度耦合系数矩阵为[Ωi1  Ωi2]3×6,耦合力可表示为:
(3)
式中,rSii=1,2)表示质量块在航天器内的偏移量,下标i取1、2分别表示质量块1、2。则有质量块1与质量块2方程相对腔体方程如下:
r1S¨=ω˙s×-r1s+rOl+ωs×ωs×-r1s+rO1-2ωS×r˙1S+d1S+u1Sm1-dSC+uSCmSC+ΔG1+F1
(4)
r¨2S=ω˙S×-r2S+r02+ωS×ωS×-r2S+r02-2ωS×r˙2S+d2S+u2Sm2-dSC+uSCmSC+ΔG2+F2
(5)
式中, rOi表示腔体到质量块位移,dSiuSi、ΔGi分别表示质量块所受非保守力扰动、控制和引力梯度,dSCuSC分别表示航天器所受扰动和控制。
表示系统状态变量,Ai1Ai2表示式(4)~(5)中相对运动方程系数矩阵,Ωi1Ωi2表示刚度耦合系数矩阵,前文所分析系统总扰动标记为ωi,其中i取1,2用以表示不同质量块所处系统,以上矢量均表示在航天器体坐标系下。本文研究中,忽略航天器指向姿态误差,假定呼吸角为2α,电容位移传感器输出在腔体坐标系,转化到航天器坐标系,r1S=T1r1r2S=T2r2,其中转换矩阵T1=cosα0sinα010-sinα0cosαT2=cos-α0sin-α010-sin-α0cos-α,记η= T=diagT1T2T1T2,diag函数表示对角矩阵,则有:
Tη˙=03×303×3I3×3I3×3A11+Ω1103×3A12+Ω1203×303×3A21+Ω2103×3A22+Ω22Tη+06×1-I3×3mSCuSC-I3×3mSCuSC+06×1I3×3mTMu1I3×3mTM u2+06×6I6×6ω1ω26×1y=J3×12Tη
(6)
式(6)记为双检验质量平动系统动力学方程,为方便,将式(6)记为:
η˙=Adfη+BdfuuSC+BdfmuTM+Bdfdωy=J3×12Tη
(7)
1.2 选取标准检验质量追踪策略
传统的在轨无拖曳控制以一颗检验质量为追踪基准,通过微推补偿航天器所受非保守力,但在多检验质量系统中,各质量块和航天器分属于不同轨道,引入的引力梯度带来难被补偿的摄动,因此这种策略实现难度增大。双检验质量无拖曳系统中,记c=cosα,h=sinα,假定TM2为追踪基准,式(7)中航天器追踪系数矩阵为:
(8)
1.3 LOS点追踪策略
以两条激光链路交点作为追踪基准,策略的逻辑即为通过航天器追踪的形式保证敏感轴正交轴向上的自由,质量块的敏感轴上的自由移动对航天器的追踪控制无影响,其几何关系示意图见图3
LOS点的追踪策略,可以对两颗质量块的X轴外所有自由度进行补偿,航天器在XOZ面内自由度的追踪可以充分解耦到质量块非敏感轴即Z1Z2轴进行补偿。XOZ面内的LOS点几何关系方程如下:
3LOS点追踪策略示意图
Fig.3Schematic diagram of LOS point tracking strategy
d-z1+d+z=3dx1+3d-xdz2+d-z=3dx2+3d-x
(9)
其中,(xz)为LOS点的偏移坐标量,(x1z1)、(x2z2)为检验质量的坐标偏移量,2d为两检验质量中心初始距离。c3=3ch3=3h,解算结果如下:
J3×12T=c-h320h+c32c-h320-h-c3201200120 c-h3230h+c323c-h3230h+c323 03×6
(10)
1.4 无固定追踪点策略
无固定追踪点策略记为NDP(no defined drag-free point)策略,策略期望系统转化为敏感轴面即XOZ面内多点追踪问题,通过补偿航天器来平衡多检验质量与航天器间的相对运动,同时通过悬浮回路限制质量块在腔体的偏移量,式(11)为追踪矩阵。
J3×12T=c20h2c20-h20120012003×6-h20c2h20c2
(11)
这种追踪策略航天器在体坐标系三轴对纯引力基准点补偿,质量块悬浮回路控制信号对质量块在腔体内的非敏感轴偏移量进行辅助性的补偿。
1.5 检验质量平动动力学模型
ΩM1=λij23×6ΩM2=μij23×6为质量块1、2体坐标系下的耦合系数矩阵,刚度系数矩阵有如下基本结论[1318]。刚度耦合矩阵的主导因素为腔体内执行机构与质量块间的耦合,在腔体内沿各轴的刚度耦合摄动可忽略与相对速度相关系数,且与位置矢量相关系数对其他轴的串扰影响量级比所在主轴系数量级小两阶以上,即在分析中,耦合影响因素以矩阵前三列对角元素为主导。同时检验质量作为追踪基准,在理想情况下沿纯引力轨道运动,无拖曳系统中,由于系统耦合项的量级远高于深空轨道相对运动方程系数量级,质量块运动可简化为各轴上弹簧系统[13],刚度耦合系数体现为质量块受摄系数,式(7)中Adf项可简化为耦合系数矩阵的对角项。
以质量块1为例,检验质量动力学表示为:
(12)
ΩM1与速度关联项系数为0,与位移关联项各轴耦合函数为0,传递函数形式表示为:
PM1(s)=1mTMdiag1s2-λ112,1s2-λ222,1s2-λ332
(13)
式中,s为控制学科传递函数通用频域自变量。质量块的动力学被充分解耦,其传递函数形式是一个经典的二阶非最小相位系统。
1.6 系统指标分析
引力波探测任务科学模式中,位于航天器内的两颗检验质量提供测距基准,作为基准的检验质量理想情况下沿纯引力轨道运动,因此对于非引力项,都为质量块上的残余噪声,残余加速度噪声总体需求为3×10-15 m·s-2·Hz-1/2以下,引力波探测频段为1 mHz到1 Hz[313]。残余噪声通常分为质量块所受非保守力项、质量块与航天器间耦合力项、质量块执行机构输出,质量块内部非保守力项又会向下分级,包括残余气体、宇宙射线、残余电荷引发的电磁干扰等[16]。质量块内部非保守力为硬件指标,体现在对质量块的电荷管理、腔体内残余气体控制等方面,在本文后续仿真中,给出指标量级的白噪声来代替质量块内部非保守力,用以模拟最糟糕的内部条件。
式(3)中所分析耦合项噪声,耦合系数波动范围不超过10-7量级,因此设定相应偏移量指标,用以限制耦合噪声。敏感轴偏移量指标为3×10-9 m·Hz-1/2,非敏感轴指标为9×10-9 m·Hz-1/2,而在控制设计中,考虑耦合系数为10-6量级,保证设计具有裕度[13]
质量块悬浮回路控制输出指标,敏感轴为10-16 m·s-2·Hz-1/2。为了防止非敏感轴带来的串扰,对非敏感轴各项变量引入指标限制,非敏感轴质量块悬浮回路控制输出指标[12]为10-13 m·s-2·Hz-1/2
2 检验质量控制律设计
由于空间任务对相关变量的功率谱密度结果有着极为严格的指标约束,本文将基于频域设计方法对系统进行控制律计算。
依据前文所述刚度系数有关结论,分析耦合影响因素以矩阵前三列对角元素为主导。以十万千米近圆轨道为背景[3],取ωS=[0  1  0]Tω0ω0为航天器角速度标量,设计中忽略非对角元。在此简化基础上,以两质量块在腔体内位移偏移量[r1  r2]T作为系统输出,将系统解耦为航天器无拖曳控制回路和质量块悬浮控制回路。先对检验质量悬浮回路进行单轴控制律设计,针对不稳定悬浮系统,通过闭环反馈消除其不稳定极点;后对航天器进行三轴的平动控制律计算,航天器控制回路中,由于检验质量的控制信号受前述指标约束,量级远低于航天器控制信号,因此无拖曳控制回路设计中,将质量块控制信号视为微小扰动。
控制系统框图见图4,检验质量悬浮回路闭环构成稳定系统,此回路输出检验质量在腔体内的偏移量,解算两颗检验质量共6个自由度获得航天器追踪基准,再通过航天器控制律形成闭环,达成对基准的追踪,在设计中,要求无拖曳系统回路带宽高于0.1 Hz,悬浮系统回路带宽远低于无拖曳回路[12]
4控制系统框图
Fig.4Block diagram of the control system
2.1 质量块平动控制律
由于无拖曳控制需求,需静电悬浮控制在频域上满足极为苛刻的指标,因此系统需提供探测频段上的超低带宽控制,基于这一原则,参照回路成形理论的控制律结构[19],基于频域设计方法对质量块平动控制律设计。针对系统广义对象和前述的解耦处理,设计中三轴的名义对象均为:
1mTMGM(s)=1mTM1s2-λdiag 2
(14)
式中,λdiag表示耦合系数中位移关联项矩阵的对角元。为了保证探测频段内的性能,同时保证控制律输出在加速度指标之下,选择权值矩阵:
W1=ρsW2=s+0.0001s+0.01
(15)
为了简化控制律,系统输入端权值矩阵为常数,便于直接调整系统控制增益,输出端权值矩阵选取为一阶,且考虑质量块控制输出集中在低频段。
假定控制律usus =1s+ks,质量块位置动力学系统开环传递函数表示为:
L(s)=W11mTMGM(s)ususW2=ρs1mTM1s2-λdiag 21s+kss+1s+0.01
(16)
假定系统处于极度不稳定状况,取λ2diag=10-6mTM=1.45 kg,依据系统开环截止频率进行参数整定,配置ks=0.1,ρs=0.000 1,控制律表示为:
um=W2W1usus=0.001s+0.1s+1s+0.01
(17)
系统开环截止频率为0.005 Hz,闭环回路带宽为8.5×10-5 Hz,此时系统闭环传递函数表示为:
CM(s)=GM(s)mTM+um(s)GM(s)
(18)
式中,um即为质量块平动控制律。作为非最小相位系统,需配置闭环极点在复平面左半面,则可从闭环极点配置证明稳定性,此处系统闭环极点为-0.105 1,-0.002 4+0.020 8i,-0.002 4-0.020 8i,-1.467×10-5。易证系统参数摄动范围-1×10-6λdiag<0,闭环极点均在复平面左半面。即检验质量三轴控制律表示为:
uTM=diag0,-um,-um
(19)
2.2 航天器平动控制律
本节将结合前述所建立追踪策略,进行航天器平动控制律计算。检验质量悬浮回路动力学表示为:
PMi(s)=diagGM(s),CM(s),CM(s)
(20)
J3×12T前六列为Jr=[Jr1  Jr2],在图4中,针对腔体-检验质量系统,设计中期望降低航天器控制律构成的闭环系统对悬浮回路的依赖,因此腔体-检验质量系统名义对象选取各检验质量系统开环传递函数GM,即设计名义对象为:
Jrr1(s)r2(s)=Jr-T1-1mSCI3×3PM(s)-T2-1mSCI3×3PM(s)u(s)
(21)
2.2.1 非敏感轴控制律
基于前述刚度耦合系数基本结论、所建立方程(7)以及追踪策略式(8)、(10)~(11),对于非敏感轴向的追踪即选取质量块1、2偏移量的平均,同时由于Y轴受到的引力梯度较小,致使此轴相对运动量级较另外两轴小,因此将此轴单独设计,设计中期望于通过航天器控制律构成的闭环消除悬浮回路中的不稳定极点,降低系统对悬浮回路的依赖。考虑模型输入端摄动问题,基于频域H混合灵敏度PS/T问题[20]进行控制律求解,设计名义对象选取为GSCys=-1mSC.1s2-λdiag 2,优化目标minW1PSW2T
设计问题中权值矩阵W1W2选取期望于系统带宽后闭环性能衰减,同时考虑未建模动态及不确定性,确保鲁棒性,权值函数包含参数ρ1yρ2y用于调整增益,非敏感轴系统控制输出带宽要求较低,依此调整权值函数零极点,对于不同探测背景,可依据不同带宽选取权值函数进行整定,此处选取:
W1y=ρ1ys+20s+0.000001W2y=ρ2y-s2+0.002s+0.005(0.00001s+1)2
(22)
和日心轨道相异的是,地心轨道背景会带来较大的非敏感轴向太阳光压,需提高系统开环截止频率,取mSC=250 kg,计算过程中参数整定ρ1y=ρ2y=1,闭环极点为-63.257,-63.234,-2.377,-1.191+2.058i,-1.191-2.058i,-0.083+0.075i,-0.083-0.075i,闭环系统带宽为0.123 8,非敏感轴控制律为:
uy=0.15s4+3055s3+1.5×107s2+2.4×106s+167822s5+131.4s4+4635s3+21246s2+50404s+699.5
(23)
2.2.2 敏感轴控制律
依据追踪策略式(8)、(10)~式(11),结合式(14)可计算航天器敏感轴向设计名义对象。考虑最糟糕情况λ2diag=10-6,其传递函数形式将被表示为:
rx(s)rz(s)=-0.004(s+0.001)(s-0.001)1 1ux(s)uz(s)
(24)
与非敏感轴相比,需略抬高系统开环截止频率,选取权值函数:
W1xz=s+60s+0.000001ρ1xρ1zW2xz=-s2+0.002s+0.005(0.00001s+1)2ρ2xρ2z
(25)
选取参数ρ1x=ρ1z=2,ρ2x=ρ2z=1。配置闭环极点在复平面左半面,系统对角元的闭环极点为-63.258 4,-63.232 7,-2.636 9,-1.281 3+2.270 4i,-1.281 3-2.270 4i,-0.099 0+0.131 8i,-0.099 0-0.131 8i,闭环系统带宽为0.164 6,敏感轴向控制律结果为:
uxz=0.2s4+4111s3+2.1×107s2+3.9×106s+48691s5+132s4+4697s3+23458s2+61052s+1266
(26)
控制律可表示为:
diag0.2s4+4111s3+2.1×107s2+3.9×106s+48691s5+132s4+4697s3+23458s2+61052s+1266,0.15s4+3055s3+1.5×107s2+2.4×106s+167822s5+131.4s4+4635s3+21246s2+50404s+699.50.2s4+4111s3+2.1×107s2+3.9×106s+48691s5+132s4+4697s3+23458s2+61052s+1266
(27)
3 数值仿真及结果分析
航天器及检验质量仿真相关参数见表1,航天器初始轨道根数见表2。数值仿真中的初始条件设定为,质量块1、2在腔体内的三轴初始偏移量为100 nm,初始相对速度为1 nm/s。仿真中,质量块内部非保守力考虑前文所述最糟糕情况为仿真条件,即在非敏感轴为10-13 m/s2(全频段),在敏感轴为10-16 m/s2(全频段),利用该量级的白噪声来仿真模拟。本文分析中频域性能采用功率谱密度(power spectral density,PSD)指标。
1航天器及质量块相关参数
Tab.1 Parameters related of spacecraft and mass blocks
2航天器初始轨道根数[3]
Tab.2 Initial orbit elements of spacecraft
3.1 摄动仿真
本节对系统摄动项进行仿真,评估其频域性能。
对于十万千米高度的轨道背景,航天器所受非保守力以太阳光压为主,摄动结果如图5所示,引力梯度以中心天体引力梯度和日月星体引力梯度为主,引力梯度如图6所示。
5航天器所受太阳光压摄动
Fig.5Solar radiation pressure perturbation of spacecraft
6星体引力梯度
Fig.6Gravity gradient caused by celestial body
3.2 标准检验质量追踪策略
检验质量平行放置的双检验质量系统,引力梯度影响集中在单条敏感轴向,这类系统中,以一颗检验质量为追踪基准,可同时达成对两个检验质量的追踪,但在引力波探测航天器中很难实现。由于两质量块相对姿态稳定后,测距轴向X1X2轴相交成呼吸角,Z1Z2轴也无法沿同一方向平行,使质量块1、2的XZ轴向耦合问题变得复杂。以TM2为追踪对象,航天器补偿TM2三轴在腔体中偏移量(x2y2z2)。此时对于TM2,情况等效于单检验质量无拖曳系统的控制问题。
其控制结果见图7,TM2敏感轴的偏移量限制在3×10-9 m·Hz-1/2以下,非敏感轴偏移量限制在9×10-9 m·Hz-1/2。事实上,对TM2追踪控制结果,也是控制律对单检验质量无拖曳系统可行性的一种验证。仿真结果表明,这种追踪策略下,追踪基准检验质量的频域性能得到保证,但非基准检验质量很难在无额外补偿信号的情况下满足要求,其时域偏移不断发散,频域性能也无法满足科学任务技术指标需求,以TM1的Z轴向为例,偏移量结果如图8所示。
7质量块2三轴偏移量功率谱密度结果
Fig.7PSD results of offset along three axes on TM2
3.3 LOS点追踪策略
LOS点追踪策略,逻辑在于利用航天器对LOS点的平动追踪,对两颗检验质量非敏感轴向的所有自由度进行补偿。非敏感轴偏移结果见图9。从图9可以看出LOS点对非敏感轴追踪非常稳定,频域性能满足指标需求,偏移量均限制在9×10-9 m·Hz-1/2以下。通过航天器质心对LOS点的追踪,可以对两条敏感轴外所有自由度进行补偿,并保证系统非敏感轴向的频域性能达到所需技术指标。
8质量块1的Z轴偏移量
Fig.8Offset of TM1 along Z axis
9质量块1、2非敏感轴位移偏移量功率谱密度
Fig.9Displacement PSD results of offset along non-sensitive axes on TM1 and TM2
质量块1、2敏感轴偏移量如图10所示,这种策略中两检验质量敏感轴向偏移量将不断发散。后续通过引入两检验质量敏感轴额外信号进行补偿的方式可改善系统敏感轴时域偏移量和频域性能,情况相对于前述标准检验质量追踪策略有所简化。但在本文探究前提下,敏感轴不采取补偿信号时,此跟踪策略实现难度较大。
10质量块1、2敏感轴偏移量
Fig.10Offset of TM1 and TM2 along sensitive axes
3.4 无固定追踪点策略
无固定追踪点策略期望系统转化为敏感轴面的多点追踪问题,通过补偿航天器来平衡多检验质量与航天器间相对运动,同时通过悬浮回路限制质量块在腔体的偏移量。此节仿真时长高于航天器一个轨道周期,时间选为600 000 s。检验质量1、2仿真结果见图11
11质量块1、2的三轴位移偏移量
Fig.11Offset of TM1 and TM2 along all of the axes
这种追踪策略,航天器在体坐标系三轴对纯引力基准点进行补偿,质量块悬浮回路控制将对质量块在腔体内非敏感轴偏移量进行限制。图12为系统频域性能,包括检验质量1、2腔体内三轴偏移量功率谱和非敏感轴控制输出功率谱。结果表明非敏感轴控制保持在10-13 m·s-2·Hz-1/2前提下,检验质量腔体内偏移量可限制在敏感轴3 nm·Hz-1/2以下,非敏感轴9 nm·Hz-1/2以下。
对于这种追踪策略,可基于前述指标需求,引入敏感轴悬浮控制信号。为了探究可行的敏感轴信号对系统性能的影响,将信号保持在全频段10-15 m·s-2·Hz-1/2以下,输出信号见图13。以质量块1的X轴为例,引入敏感轴控制信号后时域偏移见图14(a),频域性能见图14(b)
仿真结果可见,由于敏感轴悬浮输出信号受指标约束,其对系统性能影响较小。敏感轴无补偿信号前提下,无固定基准点方案可以稳定提供较长的观测周期(在仿真的一个航天器周期时长内均可满足探测任务技术指标),质量块控制信号保持在指标需求(非敏感轴10-13 m·s-2 ·Hz-1/2)以下,质量块的运动在探测频段保持在敏感轴3 nm·Hz-1/2以下,非敏感轴9 nm·Hz-1/2以下。
12无固定追踪点策略系统频域性能
Fig.12Frequency performance of the NDP strategy system
13质量块1、2敏感轴静电执行机构输出功率谱
Fig.13PSD results of suspension actuators along sensitive axes on TM1 and TM2
14有无敏感轴悬浮对比图
Fig.14Suspension comparison diagram with and without sensitive axis
4 结论
本文针对引力波探测编队航天器内特有的双检验质量无拖曳系统进行了分析和设计,并得出了地心轨道背景下的初步结论。在无额外补偿信号前提下,标准检验质量追踪策略更适用于平行放置的双检验质量系统,如LISA Pathfinder卫星,LOS点追踪策略对非敏感轴的控制达到最优,但需在敏感轴向上施加额外控制信号来限制检验质量在测距轴的偏移量,无固定追踪点策略在激光链路成固定角度的呼吸角情况下,能够将质量块敏感轴非敏感轴悬浮控制信号保持在10-13 m·s-2·Hz-1/2以下的同时,使质量块的敏感轴平动量级保持在3 nm·Hz-1/2以下,达成任务所需技术指标。
随着太极一号[21-22]和天琴一号[23-28]的成功发射,无拖曳控制技术在我国成功地获得了在轨验证,相对于目前已有在轨验证结果[2123],本文所针对的深空背景,具有更高的实现精度,但当前多检验质量无拖曳系统研究和在轨试验在国内仍处于起步阶段,有许多领域需进一步探索,包括多检验质量的成功在轨验证、星间长距离激光干涉测量等。
1双检验质量无拖曳系统示意图
Fig.1Diagram of the two-test-mass drag-free system
2三种追踪策略示意图
Fig.2Schematic diagram of three tracking strategies
3LOS点追踪策略示意图
Fig.3Schematic diagram of LOS point tracking strategy
4控制系统框图
Fig.4Block diagram of the control system
5航天器所受太阳光压摄动
Fig.5Solar radiation pressure perturbation of spacecraft
6星体引力梯度
Fig.6Gravity gradient caused by celestial body
7质量块2三轴偏移量功率谱密度结果
Fig.7PSD results of offset along three axes on TM2
8质量块1的Z轴偏移量
Fig.8Offset of TM1 along Z axis
9质量块1、2非敏感轴位移偏移量功率谱密度
Fig.9Displacement PSD results of offset along non-sensitive axes on TM1 and TM2
10质量块1、2敏感轴偏移量
Fig.10Offset of TM1 and TM2 along sensitive axes
11质量块1、2的三轴位移偏移量
Fig.11Offset of TM1 and TM2 along all of the axes
12无固定追踪点策略系统频域性能
Fig.12Frequency performance of the NDP strategy system
13质量块1、2敏感轴静电执行机构输出功率谱
Fig.13PSD results of suspension actuators along sensitive axes on TM1 and TM2
14有无敏感轴悬浮对比图
Fig.14Suspension comparison diagram with and without sensitive axis
1航天器及质量块相关参数
2航天器初始轨道根数[3]
HECHLER F, FOLKNER W M. Mission analysis for the Laser Interferometer Space Antenna(LISA)mission[J]. Advances in Space Research,2003,32(7):1277-1282.
ARMANO M, AUDLEY H, BAIRD J,et al. Beyond the required LISA free-fall performance:new LISA Pathfinder results down to 20 μHz[J]. Physical Review Letters,2018,120(6):061101.
LUO J, CHEN L S, DUAN H Z,et al. TianQin:a space-borne gravitational wave detector[EB/OL].(2015-11-07)[2022-07-01].https://doi.org/10.48550/arXiv.1512.02076.
LUO J, MEI J W, SHAO C G,et al. TianQin mission concept[C]//Proceedings of the 12th Asia-Pacific International conference on Gravitation, Astrophysics,and Cosmology,2016.
YE B B, ZHANG X F, ZHOU M Y,et al. Optimizing orbits for TianQin[J]. International Journal of Modern Physics D,2019,28(9):1950121.
HU X C, LI X H, WANG Y,et al. Fundamentals of the orbit and response for TianQin[J]. Classical and Quantum Gravity,2018,35(9):095008.
LANGE B. The drag-free satellite[J]. AIAA Journal,1964,2(9):1590-1606.
SCHARLEMANN C, BULDRINI N, KILLINGER R,et al. Qualifciation test series of the indium needle FEEP micro-propulsion system for LISA Pathfinder[J]. Acta Astronautica,2011,69(9/10):822-832.
ZIEMER J K, MERKOWITZ S M. Microthrust propulsion for the LISA mission[C]//Proceedings of the 40th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit,2004.
FICHTER W, GATH P, VITALE S,et al. LISA Pathfinder drag-free control and system implications[J]. Classical and Quantum Gravity,2005,22(10): S139-S148.
FICHTER W, SCHLEICHER A, BENNANI S,et al. Closed loop performance and limitations of the LISA Pathfinder drag-free control system[C]//Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit,2007.
GATH P, FICHTER W, KERSTEN M,et al. Drag free and attitude control system design for the LISA Pathfinder mission[C]//Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit,2004.
GATH P, SCHULTE H R, WEISE D,et al. Drag free and attitude control system design for the LISA science mode[C]//Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit,2007.
WU S F, FERTIN D. Spacecraft drag-free attitude control system design with Quantitative Feedback Theory[J]. Acta Astronautica,2008,62(12):668-682.
LIAN X B, ZHANG J X, LU L,et al. Frequency separation control for drag-free satellite with frequency-domain constraints[J]. IEEE Transactions on Aerospace Electronic Systems,2021,57(6):4085-4096.
ZIEGLER T, FICHTER W. Test mass stiffness estimation for the LISA Pathfinder drag-free system[C]//Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit,2007.
SCHUMAKER B L. Disturbance reduction requirements for LISA[J]. Classical and Quantum Gravity,2003,20(10): S239-S253.
章楚, 贺建武, 段俐. 基于扩张状态观测器的无拖曳系统参数辨识[J]. 空间科学学报,2019,39(5):670-676. ZHANG C, HE J W, DUAN L. Parameter identification of drag-free system based on extended state observer[J]. Chinese Journal of Space Science,2019,39(5):670-676.(in Chinese)
FEYEL P. Loop-shaping robust control[M]. New York: John Wiley & Sons, Inc.,2013:33-41.
KWAKERNAAK H. Robust control and H-optimization[J]. Automatica,1993,29(2):255-273.
CAI Z M, DENG J F. Satellite architecture and preliminary in-orbit experiment of Taiji-1[J]. International Journal of Modern Physics A,2021,36(11n12):2140020.
The Taiji Scientific Collaboration. China′s first step towards probing the expanding universe and the nature of gravity using a space borne gravitational wave antenna[J]. Communications Physics,2021,4:34.
LUO J, BAI Y Z, CAI L,et al. The first round result from the TianQin-1 satellite[J]. Classical and Quantum Gravity,2020,37:185013.
LIANG Z C, HU Y M, JIANG Y,et al. Science with the TianQin observatory:preliminary results on stochastic gravitational-wave background[J]. Physical Review D,2022,105:022001.
XIAO C Y, BAI Y Z, LI H Y,et al. Drag-free control design and in-orbit validation of TianQin-1 satellite[J]. Classical and Quantum Gravity,2022,39(15):155001.
王继河, 张锦绣, 孟云鹤, 等. 空间引力波探测系统编队动力学与控制技术综述[J]. 中山大学学报(自然科学版),2021,60(1/2):156-161. WANG J H, ZHANG J X, MENG Y H,et al. Review of formation dynamics and control technology of space-borne gravitational wave detection system[J]. Acta Scientiarum Naturalium Universitatis SunYatseni,2021,60(1/2):156-161.(in Chinese)
张立华, 黎明, 高永新, 等. 空间引力波探测航天器系统及平台技术[J]. 中山大学学报(自然科学版),2021,60(1/2):129-137. ZHANG L H, LI M, GAO Y X,et al. The spacecraft system and platform technologies for gravitational wave detection in space[J]. Acta Scientiarum Naturalium Universitatis SunYatseni,2021,60(1/2):129-137.(in Chinese)
李洪银, 叶小容, 刘佳恒, 等. 天琴无拖曳控制研究的关键问题[J]. 中山大学学报(自然科学版),2021,60(1/2):213-224. LI H Y, YE X R, LIU J H,et al. Key issues in the research on drag-free control of TianQin[J]. Acta Scientiarum Naturalium Universitatis SunYatseni,2021,60(1/2):213-224.(in Chinese)