CN113177269A - 一种多星分离安全距离参数优化方法 - Google Patents
一种多星分离安全距离参数优化方法 Download PDFInfo
- Publication number
- CN113177269A CN113177269A CN202110733652.6A CN202110733652A CN113177269A CN 113177269 A CN113177269 A CN 113177269A CN 202110733652 A CN202110733652 A CN 202110733652A CN 113177269 A CN113177269 A CN 113177269A
- Authority
- CN
- China
- Prior art keywords
- separation
- deviation
- satellite
- parameter
- point
- 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
Links
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种多星分离安全距离参数优化方法,以Adams软件为初始分离段分析工具,将各分离体真实外形下的质量特性、分离弹簧力、后效推力等输入参数的偏差工况,利用目标函数组合计算姿态干扰和两体最小间隙的单偏差叠加和偏差工况,以获取姿态偏差和易碰点组合,将Adams软件计算中获取的姿态偏差和易碰点组合作为六自由度运动学计算程序的输入参数,利用输入参数的偏差数据的概率分布建立随机数,开展随机打靶仿真计算,获取概率分布条件下的近场分离姿态偏差分布和星‑星、星‑箭安全距离分布,本发明克服了对各输入参数连续变化情况下的计算结果覆盖性不足以及对初始分离段真实外形下两体最小间隙精细化不足的问题。
Description
技术领域
本发明涉及卫星分离技术领域,具体涉及一种多星分离安全距离参数优化方法。
背景技术
随着星座计划的发展,卫星组网发射需求日益强烈,卫星设计也向着小型化、高集成度的方向发展。一箭多星发射提升了运载效率,但也对卫星布局、星箭分离装置和星箭分离方案设计提出了较高的要求。
卫星分离过程中,除了轨道精度外,关注的关键技术指标包括分离姿态角速度偏差和与其他分离体的安全距离。分离姿态角速度偏差的主要决定因素是分离弹簧力偏差与质心位置偏差,而不同批次释放卫星之间、卫星与末级箭体之间的安全距离受到发动机后效推力、分离前箭体质量特性偏差的影响,如固体发动机的后效推力存在较大的不确定性,预示结果难以精确模拟飞行状态;分离前箭体质量特性与称重精度、末修推进剂消耗量密切相关,质心偏差也存在较大的不确定性,这就造成了卫星多次释放后星-星、星-箭之间的安全距离预示结果受到诸多因素的影响。
通常的分离系统设计中由于分离体数量较少,采用偏差组合工况仿真分析实现了各输入参数在上、下偏差和标准状态下的特征量计算,往往能够覆盖所有恶劣情况,但是在现有技术中,星箭分离系统特别是多星、多批次分离系统设计中,分离体数量众多导致偏差组合工况的模式难以充分遍历,组合往往不能覆盖所有最恶劣情况。
发明内容
本发明的目的在于提供一种多星分离安全距离参数优化方法,以解决现有技术中分离体数量众多导致偏差组合工况的模式难以充分遍历,组合往往不能覆盖所有最恶劣情况的技术问题。
为解决上述技术问题,本发明具体提供下述技术方案:
一种多星分离安全距离参数优化方法,包括以下步骤:
步骤100、在Adams软件中建立分离坐标系,将分离卫星三维模型和末级箭体三维模型导入所述分离坐标系并在所述分离坐标系中将所述分离卫星三维模型和末级箭体三维模型按照布局位置进行装配;
步骤200、三维模拟每个分离卫星与末级箭体装配时的连接点,且在每个所述分离卫星与末级箭体的连接点模拟所述分离卫星与末级箭体分离时分别受到的推力作用参数,且在装配后的所述分离卫星和末级箭体上选择沿分离运动方向的多个边界点作为易碰点,并记录每个所述边界点在所述分离坐标系的三维坐标值;
步骤300、所述Adams软件根据分离时序和调姿时序模拟所述分离卫星和末级箭体在多星分离初始段的第一分离状态参数,分别设定所述Adams软件得到的推力作用参数和所述第一分离状态参数的偏差范围,生成所述推力作用参数对应偏差范围以及所述第一分离状态参数对应偏差范围的随机数序列;
步骤400、基于所述推力作用参数和第一分离状态参数的偏差范围选取间隙最小的易碰点组成易碰点组合;
步骤500、将所述随机数序列以及所述易碰点组合中每个所述易碰点的三维坐标值作为所述六自由度运动学计算程序的输入数据,从所述随机数序列内任意选择偏差数据组合作为多星分离试验的输入数据,所述六自由度运动学计算程序按照分离时序和调姿时序对多星分离工作进行分段模拟,计算所述分离卫星和末级箭体在下一次卫星分离时的运动学参数,且计算所述运动学参数在所述分离卫星自身建立的坐标系内对应的坐标值;
步骤600、定义坐标变换函数,将六自由度运动学计算程序得到的运动学参数通过所述坐标变换函数变换至所述分离坐标系,计算每个所述分离卫星的运动学参数在所述分离坐标系下对应的第二分离状态参数,并确定该换算的第二分离状态参数与所述第一分离状态参数的相对值;
步骤700、循环迭代所述随机数序列内不同的随机数,重复步骤500-步骤600以获得每一次循环中的相对运动参数计算结果,统计相对值的概率分布情况。
作为本发明的优选方案,在步骤200中,在每个所述分离卫星与末级箭体的连接点模拟所述分离卫星与末级箭体分离时分别受到的推力作用参数的实现方式为:
在每个所述分离卫星与末级箭体的连接面中心位置设置连接点;
在每个所述分离卫星与末级箭体上的连接点之间设置分离弹簧处,且在每个所述分离弹簧处设置与位移相关作用力来模拟弹簧力;
在发动机喷管摆动中心点设置具有轴向分力和横向分力的发动机后效推力以模拟发动机轴向推力和残余摆角引起的横向推力。
作为本发明的优选方案,在步骤300中,所述推力作用参数和所述第一分离状态参数包括力参数、质量特性参数以及空间位置参数,其中所述力参数表示弹簧力和后效推力、所述质量特性参数表示质量惯量以及转动惯量,所述空间位置参数表示质心位置、弹簧力位置以及后效推力位置。
作为本发明的优选方案,在步骤400中,所述易碰点选取为临近所述分离卫星的边缘点、或所述分离卫星与末级箭体之间的边缘点、或两个所述分离卫星之间的最小间隙点、或所述末级箭体的末级喷管的边缘点。
作为本发明的优选方案,所述第一分离状态参数中关于所述分离卫星以及所述末级箭体的数据偏差包括所述分离卫星的质量偏差、转动惯量偏差以及质心位置偏差;以及所述末级箭体的质量偏差、转动惯量偏差以及质心位置偏差。
作为本发明的优选方案,在步骤400中,所述随机数序列包括由所述推力作用参数下偏差和下偏差为边界形成的概率分布下的多组随机数,以及由所述分离卫星以及所述末级箭体的第一分离状态参数下偏差和下偏差为边界形成的概率分布下的多组随机数,所述Adams软件利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定易碰点组合的实现步骤为:
设置星-星分离时的弹簧力和力偏差、后效推力及力偏差;
设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
设置箭体质量偏差、转动惯量偏差以及质心位置偏差;
生成偏差工况组合,并确定组合内的元素数量;
按照分离时序和调姿时序,针对每个工况开展单次仿真计算;
获取偏差组合工况下的各分离体最小间隙、质心最小距离,选择最小间隙的点-点组合作为易碰点组合。
作为本发明的优选方案,所述Adams软件利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定分离姿态角速度的实现步骤为:
设置星-星分离时的弹簧力和力偏差;
设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
生成偏差组合工况,并确定组合内的元素数量;
针对每个偏差组合工况开展单次仿真计算,获取偏差组合工况下的卫星分离姿态角和姿态角速度。
作为本发明的优选方案,在步骤500中,按照分离时序、分离后调姿时序对星箭分离程序进行分段,每段结束时刻的末级速度、位置、姿态信息作为下一段星箭分离前的初始条件;
所述运动学参数包括加速度、线速度、位移、角速度、飞行角度。
本发明与现有技术相比较具有如下有益效果:
本发明通过Adams软件和六自由度运动学计算程序相结合,将Adams软件计算中获取的姿态偏差和易碰点组合作为六自由度运动学计算程序的输入条件和参考依据,利用分离卫星和末级箭体的输入参数(质量特性、分离弹簧力以及后效推力)偏差工况的概率分布建立随机数,实现多体之间运动学参数的随机打靶仿真,弥补偏差工况组合有限次计算可能带来的不覆盖性,提高多星分离安全性计算的工况覆盖性,为多星分离方案设计提供更加有效的仿真计算手段。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例提供安全具体参数优化方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种多星分离安全距离参数优化方法,不同批次释放卫星之间、卫星与末级箭体之间的安全距离受到发动机后效推力、分离前箭体质量特性偏差的影响,分离前箭体质量特性与称重精度、末修推进剂消耗量密切相关,质心偏差也存在较大的不确定性,采用偏差组合工况仿真分析实现了各输入参数在上、下偏差和标准状态下的特征量计算,往往能够覆盖所有恶劣情况,因此,使用参数概率分布下的随机变量作为输入,开展随机打靶分析成为越来越实用的分析手段。
Adams软件能够将真实的分离体外形直观显示,并能够实时获取两体真实外形下的最小间隙,但一套输入参数只能开展一次仿真,在输入参数偏差组合情况下,即使修改一个参数,也需要单独开展一次仿真,对多参数偏差组合工况仿真适应性不强,只能满足几十次、几百次手工修改参数下的仿真工况。
C#六自由度计算程序的缺点是不能直观表征真实分离体外形和两体最小间隙,在初始分离段的分离体间隙计算不够精细,但优点是能够通过将输入参数进行概率化分布,实现多体之间运动学参数的随机打靶仿真,弥补偏差工况组合有限次计算可能带来的不覆盖性。
因此为了结合Adams软件和C#六自由度计算程序的优点,并解决Adams软件和C#六自由度计算程序的缺点,本实施方式通过Adams软件和六自由度运动学计算程序相结合,将多星分离初始段的姿态偏差和两体最小间隙进行偏差工况计算,将计算结果、易碰点组合选取参考作为六自由度运动学计算程序的输入,开展输入参数概率分布下的随机打靶仿真,提高多星分离安全性计算的工况覆盖性,为多星分离方案设计提供更加有效的仿真计算手段。
具体包括以下步骤:
步骤100、在Adams软件中建立分离坐标系,将分离卫星三维模型和末级箭体三维模型导入分离坐标系并在分离坐标系中将分离卫星三维模型和末级箭体三维模型按照布局位置进行装配。
在Adams软件中将卫星和箭体按照布局位置进行装配,在每个卫星与箭体连接面中心位置设置Joint连接,在每个分离弹簧处设置点-点之间与位移相关的作用力模拟弹簧力;发动机后效推力作用于发动机喷管摆动中心点,共有轴向和横向两个分力,分别模拟发动机轴向推力和残余摆角引起的横向推力。
步骤200、三维模拟每个分离卫星与末级箭体装配时的连接点,且在每个分离卫星与末级箭体的连接点模拟分离卫星与末级箭体分离时分别受到的推力作用参数,且在装配后的分离卫星和末级箭体上选择沿分离运动方向的多个边界点作为易碰点,并记录每个边界点在分离坐标系的三维坐标值。
在步骤200中,推力作用参数包括表示弹簧力和后效推力的力参数,在每个分离卫星与末级箭体的连接点模拟分离卫星与末级箭体分离时分别受到的推力作用参数的实现方式为:
在每个分离卫星与末级箭体的连接面中心位置设置连接点;
在每个分离卫星与末级箭体上的连接点之间设置分离弹簧处,且在每个分离弹簧处设置与位移相关作用力来模拟弹簧力;
在发动机喷管摆动中心点设置具有轴向分力和横向分力的发动机后效推力以模拟发动机轴向推力和残余摆角引起的横向推力。
步骤300、Adams软件根据分离时序和调姿时序模拟分离卫星和末级箭体在多星分离初始段的第一分离状态参数,分别设定Adams软件得到的推力作用参数和第一分离状态参数的偏差范围,生成推力作用参数对应偏差范围以及第一分离状态参数对应偏差范围的随机数序列。
易碰点选取为临近分离卫星的边缘点、或分离卫星与末级箭体之间的边缘点、或两个分离卫星之间的最小间隙点、或末级箭体的末级喷管的边缘点。
第一分离状态参数包括表示质量惯量以及转动惯量的质量特性参数,以及表示质心位置、弹簧力位置以及后效推力位置的空间位置参数,推力作用参数包括表示力参数的弹簧力,第一分离状态参数中关于分离卫星以及末级箭体的数据偏差包括分离卫星的质量偏差、转动惯量偏差以及质心位置偏差;以及末级箭体的质量偏差、转动惯量偏差以及质心位置偏差,推力作用参数中关于分离卫星以及末级箭体的数据偏差包括分离卫星以及末级箭体受到的弹簧力偏差。
需要进一步补充说明的是,卫星分离过程中,除了轨道精度外,关注的关键技术指标包括分离姿态角速度以及与其他分离体的安全距离,分离姿态角速度的主要决定因素是分离弹簧力偏差、质量特性参数偏差以及质心位置偏差。
Adams软件利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定分离姿态角速度的实现步骤为:
设置星-星分离时的弹簧力和力偏差;
设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
生成偏差工况组合,并确定组合内的元素数量;
针对每个工况开展单次仿真计算,获取偏差组合工况下的卫星分离姿态角和姿态角速度。
步骤400、基于推力作用参数和第一分离状态参数的偏差范围选取间隙最小的易碰点组成易碰点组合。
由于弹簧力偏差、质量特性偏差、空间位置偏差会引起对末级箭体计算的姿态偏差,其偏差量级由作为输入条件的偏差量级所决定,易碰点组合的确定是在输入条件偏差选取的基础上,对存在最小间隙的点-点组合进行辨识,最终确定最小间隙的点-点组合作为易碰点组合,其中,易碰点间隙指的是分离路径上两体最可能发生碰撞的危险点之间的距离,静态安装下的间隙称为静态间隙,动态运动过程中的间隙称为动态间隙,一般来说,由于分离体存在一定的初始姿态角速度干扰,动态间隙会小于静态间隙。
Adams软件利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定易碰点组合的实现步骤为:
设置星-星分离时的弹簧力和力偏差、后效推力及力偏差;
设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
设置箭体质量偏差、转动惯量偏差以及质心位置偏差;
生成偏差工况组合,并确定组合内的元素数量;
按照分离时序和调姿时序,针对每个工况开展单次仿真计算;
获取偏差组合工况下的各分离体最小间隙、质心最小距离,选择最小间隙的点-点组合作为易碰点组合。
因此本实施方式以Adams软件为初始分离段分析工具,确定分离卫星以及末级箭体真实外形下的质量特性、分离弹簧力、后效推力的偏差工况组合,计算姿态干扰和两体最小间隙的单偏差叠加和偏差工况组合,以获取姿态偏差和易碰点组合。
将Adams软件计算中获取的姿态偏差和易碰点组合作为六自由度运动学计算程序的输入条件和参考依据,利用各分离体输入参数的概率分布建立随机数,开展随机打靶仿真计算,获取概率分布条件下的近场分离姿态偏差分布和星-星、星-箭安全距离分布。
步骤500、将随机数序列以及易碰点组合中每个易碰点的三维坐标值作为六自由度运动学计算程序的输入数据,从随机数序列内任意选择偏差数据组合作为多星分离试验的输入数据,六自由度运动学计算程序按照分离时序和调姿时序对多星分离工作进行分段模拟,计算分离卫星和末级箭体在下一次卫星分离时的运动学参数,且计算运动学参数在分离卫星自身建立的坐标系内对应的坐标值。
按照分离时序、分离后调姿时序对星箭分离程序进行分段,每段结束时刻的末级速度、位置、姿态信息作为下一段星箭分离前的初始条件。
随机数序列包括由推力作用参数下偏差和下偏差为边界形成的概率分布下的多组随机数,以及由分离卫星以及末级箭体的第一分离状态参数下偏差和下偏差为边界形成的概率分布下的多组随机数,随机数序列包括三个元素,分别为质量特性参数、力参数以及空间位置参数,偏差数据具体为星体质量偏差、转动惯量偏差、弹簧力偏差、后效推力偏差、质心位置偏差、弹簧力位置偏差、后效推力位置偏差等。
其中,运动学参数包括加速度、线速度、位移、角速度、飞行角度以及易碰点空间坐标,从而Adams软件与六自由度运动学程序之间完成易碰点空间坐标、分离速度、分离姿态、分离位移等实现数据交互。
步骤600、定义坐标变换函数,将六自由度运动学计算程序得到的运动学参数通过坐标变换函数变换至分离坐标系,计算每个分离卫星的运动学参数在分离坐标系下对应的转化数据,并确定该换算的转化数据与Adams软件获取的相同数据进行对比。
六自由度运动学计算程序按照分离时序和调姿时序对多星分离工作进行分段模拟,具体的实现方式为:
(1)选择从随机数序列内选择任一个弹簧力偏差数据、质量特性偏差数据,各分离体空间位置偏差数据及易碰点坐标作为输入参数;
(2)设置各星分离时序和调姿时序,设置分离点位置、速度、姿态信息传递方向,计算分离体在自身坐标系下的位置、速度和姿态参数
(3)将分离体运动学参数坐标变换至分离坐标系计算分离体在分离坐标系 下的位置、速度和姿态参数;
(4)计算运动学参数相对值,获取相对分离姿态、易碰点间隙、质心位置、速度等。
具体的,六自由度运动学计算程序利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定分离姿态角速度的实现步骤为:
1、设置星-星分离时的弹簧力和力偏差;
2、设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
3、建立各参数概率分布下的随机数,针对随机数开展随机打靶计算;
4、获取随机打靶工况下的卫星分离姿态角和姿态角速度。
将六自由度运动学计算程序计算得到的姿态参数、力参数、质性惯量参数以及空间位置参数与Adams软件计算得到的数据进行对比,选择差距较少的数据参数作为优化参数。
六自由度运动学计算程序利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定易碰点组合的实现步骤为:
设置分离弹簧力及力偏差后效推力及力偏差
设置卫星质量、转动惯量及偏差质心位置及偏差
设置箭体质量、转动惯量及偏差质心位置及偏差
建立各参数概率分布下的随机数,按照分离时序和调姿时序,针对随机数开展随机打靶计算
获取随机打靶工况下各分离体最小间隙、质心最小距离、相对分离速度等。
其中,第二分离状态参数同样包括运动学参数、力参数、质性惯量参数以及空间位置参数,力参数包括弹簧力以及后效推力;空间位置参数包括质心位置、弹簧力位置以及后效推力位置等,质性惯量参数包括质量和转动惯量。
根据上述六自由度运动学计算程序和Adams软件计算易碰点组合和姿态角速度的实现方式可知,Adams软件能够将真实的分离体外形直观显示,并能够实时获取两体真实外形下的最小间隙,但一套输入参数只能开展一次仿真,在输入参数偏差组合情况下,即使修改一个参数,也需要单独开展一次仿真,对多参数偏差组合工况仿真适应性不强,只能满足几十次、几百次手工修改参数下的仿真工况。
由于力参数偏差、质量特性偏差和空间位置偏差导致海量的仿真数据,因此Adams软件并不适合进行仿真工作,六自由度计算程序能够通过将力参数、质量特性参数和空间位置参数进行概率化分布,实现多体之间运动学参数的随机打靶仿真,弥补偏差工况组合有限次计算可能带来的不覆盖性。
因此本实施方式结合两种计算程序的优点,利用Adams软件能够将真实的分离体外形直观显示,并能够实时获取两体真实外形下的最小间隙,同时将力参数、质量特性参数和空间位置参数进行概率化分布,实现多体之间运动学参数的随机打靶仿真,弥补偏差工况组合有限次计算可能带来的不覆盖性。
步骤700、循环迭代随机数序列内不同的随机数,重复步骤500-步骤600以获得每一次循环中的相对运动参数计算结果,统计相对值的概率分布情况。
本实施方式实现了Adams软件与六自由度运动学计算程序的结合使用,使初始分离段和近场分离段姿态偏差和易碰点间隙、安全距离计算实现无缝衔接;克服了单纯使用Adams软件分离计算对各输入参数连续变化情况下的计算结果覆盖性问题;克服了单纯使用C#六自由度运动学计算程序对初始分离段真实外形下两体最小间隙精细化预示不足的问题。
因此本发明通过Adams软件和六自由度运动学计算程序相结合,将Adams软件计算中获取的姿态偏差和易碰点组合作为六自由度运动学计算程序的输入条件和参考依据,利用分离卫星和末级箭体的输入参数(质量特性、分离弹簧力以及后效推力)偏差工况的概率分布建立随机数,实现多体之间运动学参数的随机打靶仿真,弥补偏差工况组合有限次计算可能带来的不覆盖性,提高多星分离安全性计算的工况覆盖性,为多星分离方案设计提供更加有效的仿真计算手段。
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。
Claims (8)
1.一种多星分离安全距离参数优化方法,其特征在于,包括以下步骤:
步骤100、在Adams软件中建立分离坐标系,将分离卫星三维模型和末级箭体三维模型导入所述分离坐标系,并在所述分离坐标系中将所述分离卫星三维模型和末级箭体三维模型按照布局位置进行装配;
步骤200、三维模拟每个分离卫星与末级箭体装配时的连接点,且在每个所述分离卫星与末级箭体的连接点模拟所述分离卫星与末级箭体分离时分别受到的推力作用参数,且在装配后的所述分离卫星和末级箭体上选择沿分离运动方向的多个边界点作为易碰点,并记录每个所述边界点在所述分离坐标系的三维坐标值;
步骤300、所述Adams软件根据分离时序和调姿时序模拟所述分离卫星和末级箭体在多星分离初始段的第一分离状态参数,分别设定所述Adams软件得到的推力作用参数和所述第一分离状态参数的偏差范围,生成所述推力作用参数对应偏差范围以及所述第一分离状态参数对应偏差范围的随机数序列;
步骤400、所述Adams软件基于所述推力作用参数和第一分离状态参数的偏差范围选取间隙最小的易碰点组成易碰点组合;
步骤500、将所述随机数序列以及所述易碰点组合中每个所述易碰点的三维坐标值作为所述六自由度运动学计算程序的输入数据,所述六自由度运动学计算程序按照分离时序和调姿时序对多星分离工作进行分段模拟,计算所述分离卫星和末级箭体在下一次卫星分离时的运动学参数,且计算所述运动学参数在所述分离卫星自身建立的坐标系内对应的坐标值;
步骤600、定义坐标变换函数,将六自由度运动学计算程序得到的运动学参数通过所述坐标变换函数变换至所述分离坐标系,计算每个所述分离卫星的运动学参数在所述分离坐标系下对应的转化数据,并将该换算的转化数据与Adams软件获取的相同数据进行对比;
步骤700、循环迭代所述随机数序列内不同的随机数,重复步骤500-步骤600以获得每一次循环中的相对运动参数计算结果,统计相对值的概率分布情况。
2.根据权利要求1所述的一种多星分离安全距离参数优化方法,其特征在于,在步骤200中,在每个所述分离卫星与末级箭体的连接点模拟所述分离卫星与末级箭体分离时分别受到的推力作用参数的实现方式为:
在每个所述分离卫星与末级箭体的连接面中心位置设置连接点;
在每个所述分离卫星与末级箭体上的连接点之间设置分离弹簧处,且在每个所述分离弹簧处设置与位移相关作用力来模拟弹簧力;
在发动机喷管摆动中心点设置具有轴向分力和横向分力的发动机后效推力以模拟发动机轴向推力和残余摆角引起的横向推力。
3.根据权利要求2所述的一种多星分离安全距离参数优化方法,其特征在于:在步骤300中,所述推力作用参数和所述第一分离状态参数包括力参数、质量特性参数以及空间位置参数,其中所述力参数表示弹簧力和后效推力、所述质量特性参数表示质量惯量以及转动惯量,所述空间位置参数表示质心位置、弹簧力位置以及后效推力位置。
4.根据权利要求1所述的一种多星分离安全距离参数优化方法,其特征在于:在步骤400中,所述易碰点选取为临近所述分离卫星的边缘点、或所述分离卫星与末级箭体之间的边缘点、或两个所述分离卫星之间的最小间隙点、或所述末级箭体的末级喷管的边缘点。
5.根据权利要求1所述的一种多星分离安全距离参数优化方法,其特征在于:所述第一分离状态参数中关于所述分离卫星以及所述末级箭体的数据偏差包括所述分离卫星的质量偏差、转动惯量偏差以及质心位置偏差;以及所述末级箭体的质量偏差、转动惯量偏差以及质心位置偏差。
6.根据权利要求5所述的一种多星分离安全距离参数优化方法,其特征在于,在步骤400中,所述随机数序列包括由所述推力作用参数下偏差和下偏差为边界形成的概率分布下的多组随机数,以及由所述分离卫星以及所述末级箭体的第一分离状态参数下偏差和下偏差为边界形成的概率分布下的多组随机数,所述Adams软件利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定易碰点组合的实现步骤为:
设置星-星分离时的弹簧力和力偏差、后效推力及力偏差;
设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
设置箭体质量偏差、转动惯量偏差以及质心位置偏差;
生成偏差工况组合,并确定组合内的元素数量;
按照分离时序和调姿时序,针对每个工况开展单次仿真计算;
获取偏差组合工况下的各分离体最小间隙、质心最小距离,选择最小间隙的点-点组合作为易碰点组合。
7.根据权利要求6所述的一种多星分离安全距离参数优化方法,其特征在于,所述Adams软件利用弹簧力偏差、质量特性参数偏差以及空间位置偏差确定分离姿态角速度的实现步骤为:
设置星-星分离时的弹簧力和力偏差;
设置卫星质量偏差、转动惯量偏差以及质心位置偏差;
生成偏差组合工况,并确定组合内的元素数量;
针对每个偏差组合工况开展单次仿真计算,获取偏差组合工况下的卫星分离姿态角和姿态角速度。
8.根据权利要求1所述的一种多星分离安全距离参数优化方法,其特征在于:在步骤500中,按照分离时序、分离后调姿时序对星箭分离程序进行分段,每段结束时刻的末级速度、位置、姿态信息作为下一段星箭分离前的初始条件,所述运动学参数包括加速度、线速度、位移、角速度、飞行角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110733652.6A CN113177269B (zh) | 2021-06-30 | 2021-06-30 | 一种多星分离安全距离参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110733652.6A CN113177269B (zh) | 2021-06-30 | 2021-06-30 | 一种多星分离安全距离参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113177269A true CN113177269A (zh) | 2021-07-27 |
CN113177269B CN113177269B (zh) | 2021-09-21 |
Family
ID=76927925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110733652.6A Active CN113177269B (zh) | 2021-06-30 | 2021-06-30 | 一种多星分离安全距离参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113177269B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113919064A (zh) * | 2021-09-29 | 2022-01-11 | 中国航天空气动力技术研究院 | 并联级间分离安全性判定方法、电子设备及介质 |
CN113978764A (zh) * | 2021-11-30 | 2022-01-28 | 长光卫星技术有限公司 | 一种适用于批量组网卫星入轨的防碰撞星箭分离方法 |
CN115358098A (zh) * | 2022-10-20 | 2022-11-18 | 北京宇航系统工程研究所 | 一种远场安全性分析方法、系统、电子设备及存储介质 |
CN116244836A (zh) * | 2023-03-21 | 2023-06-09 | 蓝箭航天空间科技股份有限公司 | 基于adams的运载火箭起飞漂移安全间隙设计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060269130A1 (en) * | 2005-05-30 | 2006-11-30 | Renaud Maroy | Method of segmentation of a sequence of three-dimensional images, in particular in pharmaco-imagerie |
CA2769264A1 (en) * | 2009-07-30 | 2011-02-03 | Consulting Services Support Corporation | Core/satellite financial portfolio design methodology, system and computer readable medium |
CN105183927A (zh) * | 2015-05-11 | 2015-12-23 | 上海宇航系统工程研究所 | 一种多星分离参数优化方法 |
CN107589756A (zh) * | 2017-09-12 | 2018-01-16 | 北京理工大学 | 一种奔月卫星编队初始化方法 |
CN108516104A (zh) * | 2018-03-12 | 2018-09-11 | 中国人民解放军国防科技大学 | 多级运载航天器分离可靠性仿真平台 |
CN111538345A (zh) * | 2020-05-07 | 2020-08-14 | 上海宇航系统工程研究所 | 运载火箭星箭分离段的程序角生成方法 |
CN112445234A (zh) * | 2020-11-27 | 2021-03-05 | 航天科工火箭技术有限公司 | 一种航天器的姿态控制方法和装置 |
-
2021
- 2021-06-30 CN CN202110733652.6A patent/CN113177269B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060269130A1 (en) * | 2005-05-30 | 2006-11-30 | Renaud Maroy | Method of segmentation of a sequence of three-dimensional images, in particular in pharmaco-imagerie |
CA2769264A1 (en) * | 2009-07-30 | 2011-02-03 | Consulting Services Support Corporation | Core/satellite financial portfolio design methodology, system and computer readable medium |
CN105183927A (zh) * | 2015-05-11 | 2015-12-23 | 上海宇航系统工程研究所 | 一种多星分离参数优化方法 |
CN107589756A (zh) * | 2017-09-12 | 2018-01-16 | 北京理工大学 | 一种奔月卫星编队初始化方法 |
CN108516104A (zh) * | 2018-03-12 | 2018-09-11 | 中国人民解放军国防科技大学 | 多级运载航天器分离可靠性仿真平台 |
CN111538345A (zh) * | 2020-05-07 | 2020-08-14 | 上海宇航系统工程研究所 | 运载火箭星箭分离段的程序角生成方法 |
CN112445234A (zh) * | 2020-11-27 | 2021-03-05 | 航天科工火箭技术有限公司 | 一种航天器的姿态控制方法和装置 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113919064A (zh) * | 2021-09-29 | 2022-01-11 | 中国航天空气动力技术研究院 | 并联级间分离安全性判定方法、电子设备及介质 |
CN113978764A (zh) * | 2021-11-30 | 2022-01-28 | 长光卫星技术有限公司 | 一种适用于批量组网卫星入轨的防碰撞星箭分离方法 |
CN113978764B (zh) * | 2021-11-30 | 2023-07-21 | 长光卫星技术股份有限公司 | 一种适用于批量组网卫星入轨的防碰撞星箭分离方法 |
CN115358098A (zh) * | 2022-10-20 | 2022-11-18 | 北京宇航系统工程研究所 | 一种远场安全性分析方法、系统、电子设备及存储介质 |
CN116244836A (zh) * | 2023-03-21 | 2023-06-09 | 蓝箭航天空间科技股份有限公司 | 基于adams的运载火箭起飞漂移安全间隙设计方法 |
CN116244836B (zh) * | 2023-03-21 | 2024-05-07 | 蓝箭航天空间科技股份有限公司 | 基于adams的运载火箭起飞漂移安全间隙设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113177269B (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113177269B (zh) | 一种多星分离安全距离参数优化方法 | |
CN107679655B (zh) | 一种航天发射火箭落点预测系统 | |
CN107544067B (zh) | 一种基于高斯混合近似的高超声速再入飞行器跟踪方法 | |
CN108180910B (zh) | 一种基于气动参数不确定的飞行器快速高精度制导方法 | |
CN111325768B (zh) | 一种基于3d视觉和模仿学习的自由漂浮目标抓捕方法 | |
CN106019323A (zh) | 卫星导航接收机仿真测试系统 | |
CN110874055B (zh) | 两相流场作用下高超声速飞行器分离过程预示与控制方法 | |
CN110412868B (zh) | 一种使用星间光学图像的非合作航天器轨道确定方法 | |
CN111695193B (zh) | 一种全局相关三维气动力数学模型的建模方法及系统 | |
CN110044210B (zh) | 考虑任意阶地球非球型引力摄动的闭路制导在线补偿方法 | |
CN116382124B (zh) | 一种运载火箭的姿态控制仿真方法和系统 | |
CN106682361A (zh) | 一种基于gps仿真的无人机飞行轨迹模拟系统及方法 | |
CN114936471B (zh) | 一种基于并行计算的航天器碰撞预警分层快速筛选方法 | |
CN114239305B (zh) | 一种战场态势场景仿真激励系统 | |
CN114280934A (zh) | 一种可重复使用运载火箭全程轨迹规划方法 | |
CN112650076A (zh) | 一种星群协同控制地面仿真系统 | |
CN112964450B (zh) | 一种用于风洞多体分离试验模型轨迹预测方法 | |
CN115562314A (zh) | 运载火箭子级落区控制方法、系统、介质及计算机设备 | |
CN107883925B (zh) | 一种导航星座星间观测目标卫星图像模拟方法 | |
CN112287560A (zh) | 一种用于火箭在线轨迹规划的求解器设计方法 | |
CN116244836B (zh) | 基于adams的运载火箭起飞漂移安全间隙设计方法 | |
CN111026139A (zh) | 一种基于飞行轨迹的三维模型姿态调整控制方法 | |
CN112257172B (zh) | 一种基于离散点计算的卫星天线指向角正确性的验证方法 | |
CN113408158B (zh) | 一种适用于运载火箭级间冷分离的实现方法 | |
CN111880404B (zh) | 基于非奇异终端滑模控制的航天器集群控制方法 |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220802 Address after: 511400 room 1702, No. 1, Cuiying street, Nansha District, Guangzhou, Guangdong Province (office only) Patentee after: Guangzhou Zhongke Aerospace Exploration Technology Co.,Ltd. Address before: 100190, No. 15 West Fourth Ring Road, Beijing, Haidian District Patentee before: INSTITUTE OF MECHANICS, CHINESE ACADEMY OF SCIENCES |