小卫星分离姿态干扰及抑制分析.pdf

2023 年 5 月第 19 卷 第 2 期系统仿真技术System Simulation TechnologyMay，2023Vol.19，No.2小卫星分离姿态干扰及抑制分析王金昌，毛鹏程，吕涛，徐奎，徐方舟（中国运载火箭技术研究院，北京 100076）摘要：本研究针对小卫星分离过程中的姿态干扰问题建立了动力学模型，采用MATLAB/Simulink对动力学模型进行求解，并与ADAMS环境下的分析结果进行比较。为评估导向机构在卫星分离过程中对姿态干扰的抑制作用，进一步比较了导向机构前后的卫星分离角速度，并评估不同参数的导向机构对卫星分离姿态干扰的改善作用。本研究所采用的分析方法适用于各类卫星分离过程，可为卫星分离系统的安装布局设计、弹簧参数设计、导向机构参数设计等提供技术参考。关键词:卫星分离；导向机构；虚拟样机；动力学；ADAMS仿真Analysis of Separation Attitude Disturbance and Inhibition to Small SatelliteWANG Jinchang，MAO Pengcheng，LV Tao，XU Kui，XU Fangzhou（China Academy of Launch Vehicle Technology，Beijing 100076，China）Abstract：In this paper，a dynamic model is established for the attitude interference problem in the separation process of small satellites.This model is solved by using MATLAB/Simulink and compared with the analysis results in the ADAMS environment.In order to evaluate the suppression effect of the guidance mechanism on the attitude disturbance during the satellite separation process，the satellite separation angular velocity before and after the guidance mechanism is further analyzed，and the improvement effect of the guidance mechanism with different parameters on the satellite separation attitude disturbance is evaluated.The analysis method in this paper is applicable to various satellite separation processes，and can provide reference for the installation layout design，spring parameter design，and guide mechanism parameter design of the satellite separation system.Key words：satellite separation；guidance mechanism；virtual prototype；dynamics simulation；ADAMS simulation近年来，商业航天快速发展，小卫星具有重要的价值和作用。以美国太空探索技术公司（SpaceX）为例，该公司将计划在太空部署12 000颗卫星，提供移动互联网服务。这或将引发未来移动互联网的技术革命1-2。受发射成本等因素限制，小卫星一般采用一箭多星发射，根据需要一次发射数量可达几十颗甚至上百颗3。一箭多星任务中，小卫星安装空间受限、布局方式复杂。因此，对小卫星的分离姿态、速度、角速度，以及分离安全性都提出了较高的要求。小卫星一般采用爆炸螺栓+弹簧的分离方式，受弹簧刚度与预紧力偏差、卫星与航天器质量特性偏差、安装偏差、初始运动参数等因素影响，小卫星的分离干扰存在较大的不确定性。因此，有必要对小卫星的分通信作者：王金昌，E-mail:中图分类号:V475.1 文献标志码:A DOI:10.16812/31-1945.2023.02.008王金昌，等：小卫星分离姿态干扰及抑制分析离干扰采取抑制措施，并对其分离安全性进行分析4。星箭分离动力学建模过程复杂，多进行简化处理，例如在文献 5 中，建立了卫星分离过程的动力学模型，通过模型可直观地获取各干扰因素对卫星分离过程的影响5。该模型较为复杂，求解难度大，适用于定性评估各类干扰因素的影响。为了能够求解卫星分离后的运动参数，本研究在此建模思想的基础上进行简化与改进，建立卫星分离过程的六自由度动力学模型，分析了卫星分离过程的速度和角速度，并通过ADAMS建立虚拟样机模型，通过分析结果的一致性来评估所建立的动力学模型的准确性。为评估导向机构对小卫星星箭分离干扰的抑制作用，分析了导向机构作用下的星箭分离过程，进一步比较了导向机构前后的小卫星分离速度、角速度，验证了导向机构对卫星分离偏差的改善作用，同时也对导向机构的参数进行了优化。1 卫星分离的动力学模型建立与求解 1.1动力学模型建立图1为某一颗搭载星分离过程的简化模型，假设该卫星质量均匀分布，质心位于直角坐标系O-XYZ的原点，其相对质心的转动惯量分别为IX、IY和IZ。该搭载星通过爆炸螺栓实现解锁，通过弹簧力提供分离速度，图1中A、B、C、D4个弹簧均匀分布，弹簧顶端为等效轻质小球，与卫星底部贴紧，相邻2个弹簧之间距离为L，弹簧刚度分别为KA、KB、KC、KD，弹簧预紧力分别为FA0、FB0、FC0、FD0，分离过程中弹簧力大小分别为FA、FB、FC、FD，初始压缩量分别为XA0、XB0、XC0、XD0，4个支撑点位移分别为XA、XB、XC、XD，X轴为卫星分离方向，卫星质心沿X轴位移为x，绕Y轴转动角度为Y，绕Z轴转动角度为Z6-7。为简化模型，进行如下假设：（1）由于火箭末级与卫星相比质量足够大（超过卫星质量的100倍），火箭运动时对卫星影响小到可以忽略不计，因此分析时认为火箭位置固定不变；（2）忽略卫星分离过程中的阻力与摩擦力；（3）4个弹簧的刚度系数相等。根据4个弹簧支点的几何关系，在一定转动角度范围内，4个弹簧支撑点位移等效为|XA=L2Z+L2Y+xXB=L2Z-L2Y+xXC=-L2Z-L2Y+xXD=-L2Z+L2Y+x（1）分离过程中弹簧力大小为|FA=KA()XA0-XAFB=KB()XB0-XBFC=KC()XC0-XCFD=KD()XD0-XD（2）在式（2）中，FA、FB、FC、FD的取值范围为非负值，当XAXA0、XBXB0、XCXC0、XDXD0时，FA、FB、FC、FD取值均为0。假设卫星Y轴及Z轴方向没有作用力，X轴方向没有作用力矩，根据牛顿第二定律8，有|FA+FB+FC+FD=mxL2FA+L2FB-L2FC-L2FD=IYYL2FA-L2FB-L2FC+L2FD=IZZ（3）为简化模型，假设4个弹簧的刚度系数均为K，由式（1）（3）可得|K()XA0+XB0+XC0+XD0-4x=mxL2K()XA0+XB0-XC0-XD0-2LZ=IZYL2K()XA0-XB0-XC0+XD0-2LY=IZZ（4）1.2动力学模型求解上述模型中，设卫星质量m为50 kg，火箭末级质量M为5 750 kg，相邻2个弹簧之间距离L为0.226 m，弹簧刚度系数K为8 750 N/m，转动惯量IX为1.7 kg m2，IY为1.9 kg m2，IZ为1.9 kg m2，弹簧初始压缩量XA0为39 mm，XB0为40 mm，XC0为41 mm，XD0为42 mm。根据上述公式推导，在MATLAB/Simulink环境下建立动力学模型并进行求解。需要特别说明的是，在图1分离仿真模型Fig.1Simulation model of separation157系统仿真技术第 19 卷 第 2 期Simulink环境下，须设置弹簧力的取值范围，模拟弹簧力在卫星分离之后取值为0。经仿真分析，卫星分离后绕Y轴转动的最大角速度Y为 0.0341（）/s，绕 Z 轴转动的最大角速度Z为6.117（）/s，沿 X轴方向分离速度x为 1.072 m/s。如图2所示。在上述分析过程中简化了动力学模型，而实际卫星分离过程受前文所述的多种因素影响，若将这些因素全部考虑进去，建模及求解难度极大。分析更复杂的、多干扰要素的卫星分离过程需借助虚拟样机等手段，下面采用ADAMS建立该卫星分离的虚拟样机模型，建模时考虑了多种因素干扰，并分析导向机构对卫星分离偏差的抑制效果。2 虚拟样机的模型建立 在ADAMS环境下建立如图3所示的虚拟样机模型。在运载火箭末级上共搭载8颗小卫星，其中4颗小卫星搭载于主卫星支架内部，其余4颗搭载于运载火箭末级的前锥段。以主卫星支架内部的其中1颗小卫星为研究对象，对小卫星的分离过程进行分析。其他小卫星的分离仿真可通过编制脚本程序，通过控制小卫星与火箭之间固定副的失效时刻，来实现卫星在确定时刻的分离。本研究所有卫星的参数设置及分析方法相同，其他几颗卫星的仿真过程分析在此不再赘述。如图 4所示，作为分析对象的小卫星通过 4组弹簧+爆炸螺栓和运载火箭末级连接，小卫星体表装有导向套，运载火箭末级装有导向杆，导向套与导向杆采用间隙配合，用于消除小卫星分离过程中的姿态干扰。在 ADAMS 环境下，导向套与导向杆之间添加CONTACT 接触力约束，通过控制 CONTACT 接触力是否失效，可比较有无导向装置小卫星分离姿态的变化。弹簧中的一端与火箭末级连接，另一端连接一轻质小球，模拟弹簧顶杆，小球与卫星底面之间添加CONTACT 接触力，实现了弹簧+爆炸螺栓功能的模拟。小卫星与运载火箭之间添加Fixed固定副，通过脚本控制固定副的失效时刻来实现卫星的分离9。3 卫星分离干扰分析 按照上文动力学模型参数设置虚拟样机模型参数，分析弹簧存在偏差和电分离插头异常分离情况下，对卫星分离姿态的干扰作用。3.1弹簧偏差对卫星的分离干扰分析按照上文Simulink环境下的弹簧参数对ADAMS模型进行设置，并进行仿真分析。卫星分离后绕X轴图2卫星分离后的运动参数Fig.2The motion parameters after satellite separation图3多星分离的虚拟样机模型Fig.3Virtual prototype of multi-satellite separation system图4星箭接口定义Fig.4Interface definition of satellite and launch vehicle158王金昌，等：小卫星分离姿态干扰及抑制分析转动无角速度，绕Y轴转动的角速度Y为0.0496（）/s，绕Z轴转动的角速度Z为5.946（）/s，沿X轴方向的分离速度x为 1.060 m/s。MATLAB/Simulink环境下的仿真结果与ADAMS环境下的分析结果相比差值分别为-0.016（）/s，0.171（）/s和0.012 m/s。见图5与Simulink环境分析结果相比，在ADAMS环境下还考虑了运载火箭末级质量、卫星在火箭上偏质心安装、接触阻尼与摩擦等因素的影响，所以2种方法分析结果有小量级差值，尤其是绕Y轴的角速度?Y，由于受各种干扰因素影响在初始分离阶段差异性明显，但整体来看，2种方法分析结果差异值不大，曲线一致性较好。3.2电分离插头拉力干扰对卫星的分离干扰分析为进一步分析电分离插头拉力干扰对卫星的分离干扰，在上一分析模型中加入电分离插头机械分离时产生的拉脱力。卫星与运载火箭末级一般采用电分离插头实现电气连接，在卫星分离过程中，电分离插头机械分离产生的拉力可能会给小卫星分离姿态带来较大干扰，甚至导致碰撞事故发生10。正常情况下，电分离插头通过电信号实现自动分离。若电分离功能失效，电分离插头则采用钢索拉脱方式实现机械分离。以某一型号电分离插头为例，拉脱力峰值约98 N，卫星分离60 mm时钢索起作用，拉力作用行程不大于5 mm，拉力方向与过卫星质心轴线距离为150 mm，分析导向装置对分离插头拉力干扰的消除效果。在ADAMS环境下，分离插头拉力的模拟可通过2个Step函数叠加实现11，拉力起作用的时刻由钢索长度决定，仿真分析后的电分离插头拉力作用曲线如图6所示，该拉力作用时间约4 ms。如图7所示，电分离插头机械分离条件下，卫星分离总角速度为4.648（）/s。而电分离插头正常电分离条件下，卫星分离总角速度为5.946（）/s，大于电分离插头机械分离条件下的角速度。从卫星分离后的总角速度来看，电分离插头机械分离反而减小了卫星分离后的姿态干扰。出现这种情况的原因是，电分离插头拉力与弹簧偏差形成的干扰相互抵消。在实际应用中，弹簧偏差干扰方向无法准确测量，2种干扰的合成具有不确定性。因此，任何非确定性干扰都不是人们所期望的。4 卫星分离干扰的抑制作用分析 4.1卫星分离干扰的抑制方法减小卫星分离过程中的各类干扰可提高卫星分离姿态精度，具体方法包括：提高运载器平台的姿态控制图5卫星分离后的运动参数Fig.5The motion parameters after satellite separation图6电分离插头作用力Fig.6The force of the electrical separation plug图7电分离插头机械分离条件下的总角速度Fig.7The total angular velocity of the satellite with the electrical seperation plug mechanical separated159系统仿真技术第 19 卷 第 2 期精度、提高弹簧制造精度及其安装精度、卫星分离方向尽量靠近运载器平台质心、避免电分离插头拉力干扰、增加导向装置等12。卫星分离过程中的姿态干扰因素不确定性较大，无法精确控制，而通过增加导向装置能够有效抑制分离姿态干扰，同时能够降低卫星分离时产生的碰撞风险。尤其在多星发射任务中，分离干扰因素多，并且卫星布局方式复杂，极易造成碰撞危险，因而增加导向装置尤为重要。增加导向装置，可以提高卫星分离姿态精度，同时能够降低碰撞危险。但增加导向装置的同时，也增加了结构重量，还带来潜在的卡滞风险。因此，需要对导向装置的材料、长度进行分析。4.2导向装置的参数设计及抑制干扰能力分析卫星姿态干扰抑制一般通过增加导向装置来实现，导向装置材料和导向距离对卫星分离姿态干扰抑制有重要的影响。为了分析导向装置参数对卫星分离姿态干扰的抑制作用，可选用不同材料的导向套和不同长度的导向杆进行模拟。导向杆与导向套材料的接触力参数设置如表1所示，具体分析工况如表2所示。按照表2所列的14个工况进行分析，分析结果曲线如图8-9所示，卫星分离后的总角速度如表2所示。从分析结果可得出如下结论。（1）导向距离越长，对卫星分离后的姿态干扰抑制作用越大；但同时导向距离越长，分离系统质量越大，导向杆随着长度增加而刚度减小，从而增加分离过程中导向装置卡滞的风险。金属导向套在导向装置作用期间，卫星角速度呈现出周期振荡趋势，并且振幅衰减不明显，在0.1 s后导向装置作用已经不明显，对于金属导向套来说，导向距离100 mm即可满足要求。（2）如表2所示，卫星分离后的角速度并非随导向距离的增加而减小，而是具有一定的不确定性。主要原因是，由于在导向期间卫星的角速度呈振荡衰减趋势，导向装置在不同时刻失去作用，导致分离角速度的表1不同材料的接触力参数Tab.1The contact force parameters for different materials序号12345678参数刚度/Nm-1力指数阻尼系数/Nsm-1渗透深度/m静摩擦系数动摩擦系数静摩擦转换速度/ms-1动摩擦转换速度/ms-1钢铜11081.55104110-40.080.05110-40.01钢塑料3.810621.52103110-40.080.05110-40.01表2不同材料及导向装置长度下的分离总角速度Tab.2The separation angular velocity for different materials and guide lengths工况1234567891011121314导向距离/mm50751001251501752005075100125150175200导向杆材料钢钢钢钢钢钢钢钢钢钢钢钢钢钢导向套材料铜铜铜铜铜铜铜塑料塑料塑料塑料塑料塑料塑料分离角速度/（）/s1.9531.4120.4121.5930.4701.3651.2171.1621.4491.0621.1500.8170.5510.996图8采用金属导向套时不同导向距离的分离角速度Fig.8The separation angular velocity in different length guiding mechanism with metal guiding ring图9采用非金属导向套时不同导向距离的分离角速度Fig.9The separation angular velocity in different length guiding mechanism with nonmetal guiding ring160王金昌，等：小卫星分离姿态干扰及抑制分析不确定性。（3）卫星分离后的角速度不能够彻底消除，一方面由于卫星在火箭末级上偏质心安装，与卫星和运载火箭相互作用产生相对角速度有关；另一方面，导向装置相互作用过程中，不能够完全吸收干扰能量。即使如此，设置导向装置后的卫星分离角速度已经得到极大改善，本研究所分析的14个工况中，采用导向装置，卫星分离后的角速度减小了58%91.1%。（4）如图9所示，采用金属导向套0.1 s时卫星角速度降到0.15（）/s，而采用非金属导向套0.2 s时卫星角速度降到0.813（）/s。可见，采用非金属材料的导向套，由于接触刚度较小，卫星姿态干扰消除较慢，需要较长的导向距离，增加了分离过程中的卡滞风险。因此，提高导向杆与导向套的接触刚度能够快速消除干扰，减小导向杆的长度。5 结 论 本研究建立了一箭多星分离过程的动力学模型，通过MATLAB/Simulink实现卫星空间六自由度的分离过程仿真。该模型只适用于分离弹簧存在偏差、卫星质量远小于火箭末级质量情况下的卫星分离过程仿真。若考虑卫星分离过程的全部干扰因素并评估干扰抑制效果，须借助虚拟样机手段进行仿真分析。本研究借助ADAMS虚拟样机软件，对卫星分离过程的干扰及干扰抑制效果进行仿真分析，该分析方法适用于所有卫星分离过程的仿真分析，通过分析得出如下结论。（1）MATLAB/Simulink 环 境 下 的 仿 真 结 果 与ADAMS环境下的分析结果相比，Y、Z、x之差分别为-0.016（）/s，0.171（）/s和0.012 m/s，一致性较好。采用建立数学模型、MATLAB求解的方法适合计算简单的卫星分离问题，包括分离弹簧参数设计、分离速度计算、分离角速度预估等分析。对于复杂卫星分离过程，业界通常采用ADAMS等虚拟样机手段。（2）弹簧安装偏差、电分离插头拉力、火箭姿态干扰是卫星分离过程的主要干扰因素，以本研究设置的干扰参数为例，多种因素共同作用可引起1.298（）/s的角速度偏差，通过增加导向装置可使角速度减小91.1%，有效抑制分离干扰。（3）导向杆与导向材料应选择接触刚度较大的材料，刚度增大有利于快速降低干扰，本研究采用接触刚度 较 大 的 金 属 导 向 套 0.1 s 时 卫 星 角 速 度 降 到0.150（）/s，而采用接触刚度较小的非金属导向套0.2 s卫星角速度才降到0.813（）/s。（4）导向杆的刚度随着长度增加而减小，增加导向杆长度也增加了分离卡滞风险，同时也增加了系统重量。在有效降低干扰前提下，导向距离越短越好。最优导向距离可通过仿真求解，以本研究中的金属导向套为例，导向距离100 mm为最优导向距离。（5）受火箭末级姿态干扰及卫星安装质心偏差等因素影响，即使设置导向装置也无法将姿态干扰消除，但至少能够有效降低卫星相对于火箭末级的相对角速度，从避免碰撞和安全分离角度来说具有一定的必要性。参考文献：1吴超，谢伟.“星链”计划未来发展分析 J.国际太空，2020，498（6）：13-17.WU Chao，XIE Wei.The future development analysis of Starlink Plan J.Space International，2020，498（6）：13-17.2蒋超，王兆魁，张育林.筒式偏心在轨分离角速度抑制方法 J.航空学报，2015，36（10）：3382-3392.JIANG Chao，WANG Zhaokui，ZHANG Yulin.Angular velocity depressing method of constrained and centroid biased onorbit separationJ.Acta Aeronautica Et Astronautica Sinica，2015，36（10）：3382-3392.3吴胜宝，胡冬生.国外“一箭多星”发射现状及关键技术分析 J.国际太空，2015（10）：18-22.WU Shengbao，HU Dongsheng.Current situation and key technology of multi-payload launch missionsJ.Space International，2015（10）：18-22.4胡星志.小卫星星箭分离系统设计、分析与优化研究 D.长沙：国防科学技术大学，2012.HU Xingzhi.Research on design，analysis，and optimization of separation systems for small satellitesD.Changsha：National University of Defense Technology，2012.5袁勇，陈辉，杜善亮，等.小卫星分离过程建模与分析 J.飞行力学，2019，37（1）：71-74.YUAN Yong，CHEN Hui，DU Shanliang，et al.Modeling and analysis on separation process for small satellite J.Flight Dynamics，2019，37（1）：71-74.6卢丽颖，孟宪红，邢依琳.卫星空间分离动力学研究 J.动力学与控制学报，2014，12（2）：165-166.LU Liying，MENG Xianhong，XING Yilin，et al.Dynamic study of satellite space separation J.Journal of Dynamics and Control，2014，12（2）：165-166.（下转第169页）161邓家胜，等：混输管网冬季常温集输工艺仿真与应用研究北京：中国石油大学（北京），2017.JIA Zhiyuan.Study on the boundary conditions of unheated gathering and transportation in high water cut oilfieldD.Beijing：China University of Petroleum（Beijing），2017.9刘晓燕，王德喜，韩国有，等.特高含水采油期安全混输 温 度 界 限 试 验 研 究J.石 油 学 报，2005（3）：102-105.LIU Xiaoyan，WANG Dexi，HAN Guoyou，et al.Temperature limit for oil-gas-water mixed transportation in safety during oil production with special high water-cut J.Acta Petrolei Sinica，2005（3）：102-105.10张莹.常温输送高含水稠油粘壁机理研究 D.北京：中国石油大学（北京），2018.ZHANG Ying.Study on the sticking wall mechanism of high water-cut heavy oil D.Beijing：China University of Petroleum Beijing，2018.邓家胜 男（1986-），湖北荆州人，硕士，工程师，主要研究方向是油气田开发。何旺达 男(1992-)，辽宁海城人，学士，工程师，主要研究方向是油气田开发。7范瑞祥，唐军刚，雷凯.内支撑式串联双星分离技术J.导弹与航天运载技术，2012（4）：1-4.FAN Ruixiang，TANG Jungang，LEI Kai.Separation technology of dual satellite for internal carrying layoutJ.Missiles and Space Vehicles，2012（4）：1-4.8TENG Lai，JIN Zhonghe.A composite optimization method for separation parameters of large-eccentricity pico-satellitesJ.Frontiers of Information Technology&Electronic Engineering，2018，19（5）：685-698.9张兵，岑拯.多星分离的ADAMS仿真 J.导弹与航天运载技术，2004，269（2）：1-6.ZHANG Bing，CEN Zheng，The ADAMS simulation of multi-satellite separation system J.Missiles and Space Vehicles，2004，269（2）：1-6.10秦远田，陈卫东.电分离连接器对小卫星星箭分离干扰 分 析J.航 天 器 环 境 工 程，2012，29（3）：292-296.QIN Yuantian，CHEN Weidong.Analysis of satellite-rocket separation interfered by the electromagnetic-separation-electric-connectorJ.Spacecraft Environment Engineering，2012，29（3）：292-296.11王金昌，闫波，张佳，等.基于虚拟样机的多星分离仿真分析 J.中国空间科学技术，2016，36（6）：70-76.WANG Jinchang，YAN Bo，ZHANG Jia，et al.Multi-satellite separation simulation based on virtual prototypeJ.Chinese Space Science and Technology，2016，36（6）：70-76.12吴昌聚，徐秀琴.基于导轨结构的新型笼式皮卫星分离机构 J.浙江大学学报（工学版），2014，48（3）：548-554.WU Changju，XU Xiuqin.New cage style pico-satellite deployer based on sliding guide structure J.Journal of Zhejiang University（Engineering Science），2014，48（3）：548-554.王金昌 男（1983-），黑龙江大庆人，硕士，高级工程师，主要研究方向为飞行器结构与机构设计。毛鹏程 男（1992-），新疆库尔勒人，硕士，工程师，主要研究方向为飞行器结构与机构设计。（上接第161页）169