CN113772127B - Space debris racemization control method - Google Patents
Space debris racemization control method Download PDFInfo
- Publication number
- CN113772127B CN113772127B CN202111066585.3A CN202111066585A CN113772127B CN 113772127 B CN113772127 B CN 113772127B CN 202111066585 A CN202111066585 A CN 202111066585A CN 113772127 B CN113772127 B CN 113772127B
- Authority
- CN
- China
- Prior art keywords
- star
- target
- parent
- tether
- speed
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000006340 racemization Effects 0.000 title claims abstract description 9
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000013016 damping Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 8
- 238000012544 monitoring process Methods 0.000 claims description 4
- 230000001960 triggered effect Effects 0.000 claims description 4
- 230000005486 microgravity Effects 0.000 claims description 3
- 239000000446 fuel Substances 0.000 abstract description 10
- 239000002699 waste material Substances 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000009471 action Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 210000000078 claw Anatomy 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/66—Arrangements or adaptations of apparatus or instruments, not otherwise provided for
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G4/00—Tools specially adapted for use in space
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
Description
技术领域Technical field
本发明涉及航空航天技术领域,具体涉及一种空间碎片消旋控制方法。The invention relates to the field of aerospace technology, and in particular to a method for controlling space debris racemation.
背景技术Background technique
根据美国卫星监测网络记录显示,目前地球轨道中大于10厘米的空间碎片数量已超过两万片,这无疑对在轨航天器形成巨大威胁。柔性连接捕获法是空间碎片主动移除最有效的方法之一。该方法是指捕获目标后主星与目标是通过柔性介质连接,通常为系绳。柔性连接捕获法具有轻质,发射成本低,作用距离远等特点。常用的柔性连接捕获法包括飞网捕获法,飞爪捕获法,鱼叉法等。由于空间碎片残留角动量的存在,其往往在空间做自旋运动。利用柔性连接捕获法成功捕获此类目标后,需通过系绳对目标进行消旋处理,否则容易导致目标缠绕系绳而与主星相撞。According to records from the U.S. satellite monitoring network, the number of space debris larger than 10 centimeters in the Earth's orbit currently exceeds 20,000, which undoubtedly poses a huge threat to spacecraft in orbit. The flexible connection capture method is one of the most effective methods for active space debris removal. This method means that after capturing the target, the main star and the target are connected through a flexible medium, usually a tether. The flexible connection capture method has the characteristics of light weight, low launch cost and long range. Commonly used flexible connection capture methods include flying net capture method, flying claw capture method, harpoon method, etc. Due to the existence of residual angular momentum of space debris, it often spins in space. After successfully capturing such a target using the flexible connection capture method, the target needs to be deracinated through the tether, otherwise it will easily cause the target to wrap around the tether and collide with the main star.
在现有的绳系消旋控制系统中,控制算法往往比较复杂且传统方法需对系绳状态不断检测,一旦系绳没有拉直就打开推进器调整母星位置使系绳处于拉直状态,需要推进器的不停作动,但是系绳没有拉直只是说明主星与目标的距离比系绳拉直的状态变近,如果速度依然同步也不会发生相撞,因此依据系绳是否拉直就做出调整,导致了主星燃料的极大浪费。In the existing tether race control system, the control algorithm is often relatively complex and the traditional method requires continuous detection of the tether status. Once the tether is not straightened, the thruster is turned on to adjust the position of the parent star to keep the tether in a straightened state. The propeller needs to be constantly moved, but the tether is not straightened. This only means that the distance between the main star and the target is closer than when the tether is straightened. If the speed is still synchronized, there will be no collision. Therefore, it depends on whether the tether is straightened. Adjustments were made, resulting in a huge waste of fuel for the main star.
发明内容Contents of the invention
有鉴于此,本发明提供了一种空间碎片消旋控制方法,能够通过有限的推进器数量进行有限次数的动作,将目标的角速度在短时间范围内降低至零附近,大大节省了主星上宝贵的燃料,且防止目标与主星距离过近相撞,保证主星与目标组成的系统安全离轨。In view of this, the present invention provides a space debris despin control method, which can perform a limited number of actions with a limited number of thrusters to reduce the angular velocity of the target to near zero in a short time range, greatly saving valuable space on the main star. fuel, and prevent the target from colliding too close to the main star, ensuring the safe deorbiting of the system composed of the main star and the target.
本发明的具体技术方案如下:The specific technical solutions of the present invention are as follows:
一种空间碎片消旋控制方法,用于绳系系统捕获空间碎片,所述绳系系统为由母星、空间碎片和连接于母星与空间碎片之间的系绳组成的组合体;在对目标即空间碎片消旋的过程中,监测目标与母星的速度,当目标的速度大于母星的速度时,通过推进器调整母星的位置,防止目标与母星相撞;当目标的速度小于或者等于母星的速度时,目标与母星不会相撞,无需调整母星的位置;A method for controlling space debris decycling, which is used for capturing space debris in a tether system. The tether system is a combination composed of a parent star, space debris, and a tether connected between the parent star and the space debris; During the process of decycling the target, that is, space debris, the speed of the target and the parent star is monitored. When the speed of the target is greater than the speed of the parent star, the position of the parent star is adjusted through the thruster to prevent the target from colliding with the parent star; when the speed of the target is greater than the speed of the parent star, When the speed is less than or equal to the speed of the parent star, the target and the parent star will not collide, and there is no need to adjust the position of the parent star;
目标的速度和母星的速度均包括x,y,z三个方向的速度分量,当目标只有其中一个方向的分量速度大于母星在对应方向的分量速度时,只打开对应速度分量方向的推进器,对母星的位置做出调整。The speed of the target and the speed of the parent star both include velocity components in the three directions of x, y, and z. When the component speed of the target in only one direction is greater than the component speed of the parent star in the corresponding direction, only the propulsion in the direction of the corresponding speed component is turned on. device to adjust the position of the parent star.
进一步地,通过监测目标与母星的速度,对母星的位置做出调整,包括以下步骤:Further, by monitoring the speed of the target and the parent star, the position of the parent star is adjusted, including the following steps:
步骤一、建立组合体动力模型,包括系绳动力学模型、母星动力学模型和目标动力学模型;Step 1: Establish a dynamic model of the assembly, including a tether dynamic model, a parent star dynamic model and a target dynamic model;
步骤二、根据母星动力学模型,利用PD控制算法求解出母星的动量轮输出控制力矩,对母星的姿态进行控制;Step 2: According to the parent star's dynamic model, use the PD control algorithm to solve the parent star's momentum wheel output control torque to control the parent star's attitude;
步骤三、根据监测的速度情况,通过推进器对母星的位置进行控制,其中,根据目标与母星的速度差ΔV和调整速度差所期望的时间Δt,决定推进器的推力Fthru的大小。Step 3: Control the position of the mother star through the thruster according to the monitored speed. The thrust F thru of the thruster is determined based on the speed difference ΔV between the target and the mother star and the time Δt expected to adjust the speed difference. .
进一步地,所述系绳动力学模型表示为:Further, the tether dynamic model is expressed as:
其中,m为系绳的集中质量,r为系绳的位置矢量,为系绳位置矢量r对时间的二阶导数;n为与P点相连的节点个数;Tj为对应的系绳对P点的拉力,G为P点所受空间微重力矢量,Fk为P点所受其他外力,p为P点所受外力的个数,P点是系绳中主绳和子绳的连接点;Among them, m is the concentrated mass of the tether, r is the position vector of the tether, is the second derivative of the tether position vector r with respect to time; n is the number of nodes connected to point P; T j is the pulling force of the corresponding tether on point P, G is the spatial microgravity vector of point P, F k are other external forces on point P, p is the number of external forces on point P, and point P is the connection point between the main rope and the sub-rope in the tether;
所述母星动力学模型表示为:The parent star dynamic model is expressed as:
mSaS=FS+Fthru,m S a S =F S +F thru ,
其中,mS为主星的质量,aS为主星的加速度,FS为主星所受到的系绳拉力,Fthru为推进器推力,JS为主星的惯量矩阵,ωS为主星旋转角速度,τS为系绳拉力及扭力对主星产生的力矩,τc为主星动量轮产生的控制力矩,为主星旋转角速度ωS对时间的一阶导数;Among them, m S is the mass of the main star, a S is the acceleration of the main star, F S is the tether tension experienced by the main star, F thru is the thrust of the propeller, J S is the inertia matrix of the main star, ω S is the rotation angular velocity of the main star, τ S is the torque produced by the tether tension and torsion on the main star, τ c is the control torque produced by the momentum wheel of the main star, is the first derivative of the main star’s rotation angular velocity ω S with respect to time;
所述目标动力学模型表示为:The target dynamics model is expressed as:
mTaT=FT,m T a T = F T ,
其中,mT为目标的质量,aT为目标的加速度,FT为目标所受的系绳拉力,JT为目标在其体坐标系下惯量矩阵,ωT为目标旋转角速度,τT为子绳拉力对目标产生的力矩,为目标旋转角速度ωT对时间的一阶导数。Among them, m T is the mass of the target, a T is the acceleration of the target, F T is the tether tension exerted by the target, J T is the inertia matrix of the target in its body coordinate system, ω T is the target rotation angular velocity, and τ T is The moment caused by the pulling force of the sub-rope on the target, is the first derivative of the target rotation angular velocity ω T with respect to time.
进一步地,所述母星姿态用欧拉角ΦS=[α,β,γ]T描述,其中,按z-x-z旋转方式对应的欧拉角分别为α,β,γ;Further, the attitude of the parent star is described by the Euler angle Φ S = [α, β, γ] T , where the Euler angles corresponding to the zxz rotation are α, β, and γ respectively;
所述动量轮输出控制力矩表示为:The momentum wheel output control torque is expressed as:
其中,αSd,βSd,γSd代表期望的母星欧拉角,表示母星欧拉角对时间的一阶导数,/>表示期望母星欧拉角对时间的一阶导数,/>表示期望母星欧拉角对时间的二阶导数,参数kp>0代表比例控制增益,kd>0代表微分控制增益,τcx,τcy,τcz表示动量轮控制力矩在x,y,z三个方向的分量,JSx,JSy,JSz表示主星惯量矩阵在x,y,z三个方向的分量,τSx,τSy,τSz为系绳对主星产生的力矩在x,y,z三个方向的分量。Among them, α Sd , β Sd , γ Sd represent the expected Euler angle of the parent star, Represents the first derivative of the Euler angle of the parent star with respect to time,/> Represents the first derivative of the Euler angle of the expected parent star with respect to time,/> Represents the second derivative of the expected parent star Euler angle with respect to time, the parameter k p > 0 represents the proportional control gain, k d > 0 represents the differential control gain, τ cx , τ cy , τ cz represents the momentum wheel control torque in x, y , the components in the three directions of z, J Sx , J Sy , J Sz represent the components of the main star inertia matrix in the three directions of x, y, and z, τ Sx , τ Sy , τ Sz are the torque generated by the tether on the main star in x , components in the three directions of y and z.
进一步地,所述根据目标与母星的速度差ΔV和调整速度差所期望的时间Δt,决定推进器的推力Fthru的大小,具体为:Further, the size of the thrust F thru of the thruster is determined based on the speed difference ΔV between the target and the parent star and the time Δt expected to adjust the speed difference, specifically as follows:
所述目标与母星的速度差表示为:ΔV=[vTx-vSx,vTy-vSy,vTz-vSz],其中vTx,vTy,vTz分别表示目标速度在x,y,z三个方向的速度分量,vSx,vSy,vSz分别表示母星速度在x,y,z三个方向的速度分量;The speed difference between the target and the parent star is expressed as: ΔV=[v Tx -v Sx , v Ty - v Sy , v Tz - v Sz ], where v Tx , v Ty and v Tz respectively represent the target speed at x, The velocity components in the three directions of y and z, v Sx , v Sy , and v Sz respectively represent the velocity components of the parent star’s velocity in the three directions of x, y, and z;
所述推进器的推力:The thrust of the propeller:
其中,mS为母星的质量。Among them, m S is the mass of the parent star.
进一步地,当只考虑目标与母星沿系绳方向即x轴方向的速度分量时,ΔV=[vTx-vSx,0,0],此时,只需打开与x轴方向对应位置的推进器即可。Furthermore, when only considering the velocity component of the target and the parent star along the tether direction, that is, the x-axis direction, ΔV = [v Tx - v Sx ,0,0], at this time, it is only necessary to turn on the position corresponding to the x-axis direction. Just propeller.
进一步地,所述绳系系统还包括系绳弹射装置和弹簧阻尼单元;所述系绳弹射装置安装于所述主星上;所述系绳和所述弹簧阻尼单元贮存在系绳弹射装置的内部,当需要捕获所述空间碎片时,触发系绳弹射装置,所述系绳和弹簧阻尼单元一并弹出;Further, the tether system also includes a tether ejection device and a spring damping unit; the tether ejection device is installed on the main satellite; the tether and the spring damping unit are stored inside the tether ejection device. , when it is necessary to capture the space debris, the tether ejection device is triggered, and the tether and the spring damping unit are ejected together;
所述系绳包括主绳和子绳;所述主绳与主星相连,所述子绳与空间碎片相连;The tether includes a main rope and a sub-tether; the main rope is connected to the main star, and the sub-tether is connected to the space debris;
所述弹簧阻尼单元设置于所述主绳和所述子绳上。The spring damping unit is provided on the main rope and the sub-rope.
进一步地,所述子绳的数量为n,n条所述子绳与所述主绳连接在同一点P点;其中n为大于等于2的正整数;Further, the number of the sub-ropes is n, and the n sub-ropes and the main rope are connected at the same point P; where n is a positive integer greater than or equal to 2;
所述主绳上设置有多个弹簧阻尼单元,每条所述子绳上均设置有一个弹簧阻尼单元。A plurality of spring damping units are provided on the main rope, and one spring damping unit is provided on each of the sub-ropes.
有益效果:Beneficial effects:
(1)本发明监测目标与母星的速度,当目标的速度大于母星的速度时,通过推进器调整母星的位置,与传统方法对系绳状态检测,一旦系绳没有拉直就打开推进器调整母星位置使其拉直的方法相比,避免了目标与主星距离拉近但速度相同不会相撞时调整主星位置的燃料浪费,节省了主星燃料的消耗;同时对速度在x,y,z三个方向的速度分量进行监测,可以根据单独方向分量的速度变化情况进行调整,无须每次都打开所有的推进器,进一步节省了主星燃料,同时也更加有利于维持系统的稳定性,保证目标与母星安全离轨。(1) The present invention monitors the speed of the target and the parent star. When the speed of the target is greater than the speed of the parent star, the position of the parent star is adjusted through the thruster, and the tether status is detected using the traditional method. Once the tether is not straightened, it is opened. Compared with the method of using the thruster to adjust the position of the main star to straighten it, it avoids the waste of fuel in adjusting the position of the main star when the distance between the target and the main star is close but the speed is the same and they will not collide, thus saving the fuel consumption of the main star; at the same time, when the speed is x , the velocity components in the three directions of y and z are monitored and can be adjusted according to the velocity changes of the individual direction components. There is no need to turn on all the thrusters every time, which further saves the main star fuel and is more conducive to maintaining the stability of the system. safety to ensure the target is safely de-orbited from the parent star.
(2)充分考虑系绳的柔性,针对系绳建立动力学模型,有利于更清晰掌握系统主体的状态,在监测速度的同时,依然考虑距离因素,避免速度一直相同,但距离越来越近直接相撞的现象。(2) Fully consider the flexibility of the tether and establish a dynamic model for the tether, which will help to understand the status of the system body more clearly. While monitoring the speed, the distance factor should still be considered to avoid that the speed is always the same but the distance is getting closer and closer. A direct collision.
(3)利用PD控制算法求解出母星的动量轮输出控制力矩,在消旋以及调整母星位置的过程中,对母星的姿态进行控制,进一步确保系统的稳定性,进一步降低目标与母星相撞的可能性。(3) Use the PD control algorithm to calculate the momentum wheel output control torque of the parent star. During the process of derotating and adjusting the position of the parent star, the attitude of the parent star is controlled to further ensure the stability of the system and further reduce the distance between the target and the parent star. The possibility of a star collision.
附图说明Description of the drawings
图1为本发明的组合体模型示意图。Figure 1 is a schematic diagram of the assembly model of the present invention.
图2为消旋过程目标摆动示意图。Figure 2 is a schematic diagram of the target swing during the racemization process.
图3为目标角速度仿真变化图。Figure 3 is a simulation change diagram of the target angular velocity.
具体实施方式Detailed ways
一种空间碎片消旋控制方法,基于绳系系统捕获空间碎片,绳系系统为由母星、空间碎片和连接于母星与空间碎片之间的系绳组成的组合体,在对目标即空间碎片消旋的过程中,监测目标与母星的速度,当目标的速度大于母星的速度时,通过推进器调整母星的位置,防止目标与母星相撞;当目标的速度小于或者等于母星的速度时,目标与母星不会相撞,无需调整母星的位置。同时对速度在x,y,z三个方向的速度分量进行监测,可以根据单独方向分量的速度变化情况进行调整,无须每次都打开所有的推进器,进一步节省了主星燃料。A space debris race control method that captures space debris based on a tether system. The tether system is a combination of a parent star, space debris, and a tether connected between the parent star and the space debris. When targeting the target, that is, space During the process of debris descycling, the speed of the target and the parent star is monitored. When the speed of the target is greater than the speed of the parent star, the position of the parent star is adjusted through the thruster to prevent the target from colliding with the parent star; when the speed of the target is less than or equal to At the speed of the parent star, the target and the parent star will not collide, and there is no need to adjust the position of the parent star. At the same time, the velocity components in the three directions of x, y, and z are monitored, and adjustments can be made according to the velocity changes of the individual direction components. There is no need to turn on all the thrusters every time, further saving the main star fuel.
下面结合附图并举实施例,对本发明进行详细描述。The present invention will be described in detail below with reference to the accompanying drawings and examples.
如图1所示,为本发明的组合体模型示意图,绳系系统为由母星、空间碎片和连接于母星与空间碎片之间的系绳组成的组合体。As shown in Figure 1, it is a schematic diagram of the assembly model of the present invention. The tether system is an assembly composed of a parent star, space debris, and a tether connected between the parent star and the space debris.
在对目标即空间碎片消旋的过程中,监测目标与母星的速度,当目标的速度大于母星的速度时,通过推进器调整母星的位置,防止目标与母星相撞;During the process of decycling the target, that is, space debris, the speed of the target and the parent star is monitored. When the speed of the target is greater than the speed of the parent star, the position of the parent star is adjusted through the thruster to prevent the target from colliding with the parent star;
目标的速度和母星的速度都包括x,y,z三个方向的速度分量,当目标只有其中一个方向的分量速度大于母星在对应方向的分量速度时,只打开对应速度分量方向的推进器,对母星的位置做出调整。The speed of the target and the speed of the parent star include velocity components in the three directions of x, y, and z. When the component speed of the target in only one direction is greater than the component speed of the parent star in the corresponding direction, only the propulsion in the direction of the corresponding speed component is turned on. device to adjust the position of the parent star.
监测目标与母星的速度,对母星的位置做出调整,首先需要建立动力学模型,在控制母星姿态的同时,再对母星的位置做出调整,保证目标不会与母星相撞,具体包括以下步骤:To monitor the speed of the target and the parent star and make adjustments to the parent star's position, you first need to establish a dynamic model. While controlling the attitude of the parent star, you can then adjust the parent star's position to ensure that the target does not phase with the parent star. Collision, including the following steps:
步骤一、建立组合体动力模型,包括系绳动力学模型、母星动力学模型和目标动力学模型。Step 1: Establish a dynamic model of the combination, including a tether dynamic model, a parent star dynamic model and a target dynamic model.
系绳动力学模型表示为:The tether dynamic model is expressed as:
其中,m为系绳的集中质量,r为系绳的位置矢量,为系绳位置矢量r对时间的二阶导数;n为与P点相连的节点个数;Tj为对应的系绳对P点的拉力,G为P点所受空间微重力矢量,Fk为P点所受其他外力,p为P点所受外力的个数,P点是系绳中主绳和子绳的连接点,如图1所示。Among them, m is the concentrated mass of the tether, r is the position vector of the tether, is the second derivative of the tether position vector r with respect to time; n is the number of nodes connected to point P; T j is the pulling force of the corresponding tether on point P, G is the spatial microgravity vector of point P, F k are other external forces on point P, p is the number of external forces on point P, and point P is the connection point between the main rope and the sub-rope in the tether, as shown in Figure 1.
由于绳索单元仅能承受拉力,不能承受压力,因此绳索力学模型为Since the rope unit can only bear tension and not pressure, the rope mechanics model is
式中,l0,l分别为系绳的原长和当前长度;k为系绳的拉伸刚度,c为系绳拉伸阻尼系数,rj为P点与各连接点间的距离,则为由P点指向各连接点的单位方向向量。In the formula, l 0 and l are the original length and current length of the tether respectively; k is the tensile stiffness of the tether, c is the tensile damping coefficient of the tether, r j is the distance between point P and each connection point, Then it is the unit direction vector from point P to each connection point.
将主星即母星与目标看作刚体,由于主星可以在动量轮和推进器作用下实现对自身位姿的控制,母星动力学模型表示为:The main star, that is, the parent star and the target are regarded as rigid bodies. Since the main star can control its own posture under the action of the momentum wheel and thruster, the parent star dynamic model is expressed as:
mSaS=FS+Fthru, (3)m S a S =F S +F thru , (3)
其中,mS为主星的质量,aS为主星的加速度,FS为主星所受到的系绳拉力,Fthru为推进器推力,JS为主星的惯量矩阵,ωS为主星旋转角速度,τS为系绳拉力及扭力对主星产生的力矩,τc为主星动量轮产生的控制力矩,为主星旋转角速度ωS对时间的一阶导数。Among them, m S is the mass of the main star, a S is the acceleration of the main star, F S is the tether tension experienced by the main star, F thru is the thrust of the propeller, J S is the inertia matrix of the main star, ω S is the rotation angular velocity of the main star, τ S is the torque produced by the tether tension and torsion on the main star, τ c is the control torque produced by the momentum wheel of the main star, is the first derivative of the main star's rotation angular velocity ω S with respect to time.
目标的动力学方程与主星具有相似形式,目标动力学模型表示为:The dynamic equation of the target has a similar form to that of the host star. The target dynamic model is expressed as:
mTaT=FT, (4)m T a T = F T , (4)
其中,mT为目标的质量,aT为目标的加速度,FT为目标所受的系绳拉力,JT为目标在其体坐标系下惯量矩阵,ωT为目标旋转角速度,τT为子绳拉力对目标产生的力矩,为目标旋转角速度ωT对时间的一阶导数。Among them, m T is the mass of the target, a T is the acceleration of the target, F T is the tether tension exerted by the target, J T is the inertia matrix of the target in its body coordinate system, ω T is the target rotation angular velocity, and τ T is The moment caused by the pulling force of the sub-rope on the target, is the first derivative of the target rotation angular velocity ω T with respect to time.
其中,in,
步骤二、根据母星动力学模型,利用PD控制算法求解出母星的动量轮输出控制力矩,对母星的姿态进行控制。Step 2: According to the dynamic model of the parent star, use the PD control algorithm to solve the output control torque of the momentum wheel of the parent star to control the attitude of the parent star.
为了使母星的太阳能帆板以及天线指向不受目标捕获过程影响,应控制母星姿态在捕获目标前后保持一致。母星姿态采用欧拉角ΦS=[α,β,γ]T描述,按z-x-z旋转方式对应的欧拉角分别为α,β,γ,则有欧拉角变化率与角速度的关系为In order to prevent the solar panel and antenna pointing of the parent star from being affected by the target acquisition process, the attitude of the parent star should be controlled to be consistent before and after the target is captured. The attitude of the parent star is described by the Euler angle Φ S = [α, β, γ] T. The corresponding Euler angles according to the zxz rotation method are α, β, and γ respectively. Then the relationship between the Euler angle change rate and the angular velocity is:
式中,ωx,ωy,ωz为母星角速度在体坐标系中绕坐标轴的各分量。In the formula, ω x , ω y , and ω z are the components of the parent star’s angular velocity around the coordinate axis in the body coordinate system.
假设卫星欧拉角初始值设为ΦS0=[90°,90°,90°]T,且卫星姿态以及角速度变化范围很小,对上式进行线性化处理得到欧拉角变化率与角速度的对应关系为,Assume that the initial value of the satellite Euler angle is set to Φ S0 = [90°, 90°, 90°] T , and the satellite attitude and angular velocity change range is very small. Linearize the above equation to obtain the Euler angle change rate and angular velocity. The corresponding relationship is,
将式(3)中母星姿态动力学方程展开写成分量形式为,The parent star attitude dynamics equation in equation (3) is expanded and written into component form as,
将式(7)代入式(8)中并忽略二次项,可得简化姿态动力学控制方程为,Substituting equation (7) into equation (8) and ignoring the quadratic term, the simplified attitude dynamics control equation can be obtained as,
采用传统PD控制算法,可得出动量轮输出控制力矩如下:Using the traditional PD control algorithm, the momentum wheel output control torque can be obtained as follows:
其中,αSd,βSd,γSd代表期望的母星欧拉角,表示母星欧拉角对时间的一阶导数,/>表示期望母星欧拉角对时间的一阶导数,/>表示期望母星欧拉角对时间的二阶导数,参数kp>0代表比例控制增益,kd>0代表微分控制增益,τcx,τcy,τcz表示动量轮控制力矩在x,y,z三个方向的分量,JSx,JSy,JSz表示主星惯量矩阵在x,y,z三个方向的分量,τSx,τSy,τSz为系绳对主星产生的力矩在x,y,z三个方向的分量。Among them, α Sd , β Sd , γ Sd represent the expected Euler angle of the parent star, Represents the first derivative of the Euler angle of the parent star with respect to time,/> Represents the first derivative of the Euler angle of the expected parent star with respect to time,/> Represents the second derivative of the expected parent star Euler angle with respect to time, the parameter k p > 0 represents the proportional control gain, k d > 0 represents the differential control gain, τ cx , τ cy , τ cz represents the momentum wheel control torque in x, y , the components in the three directions of z, J Sx , J Sy , J Sz represent the components of the main star inertia matrix in the three directions of x, y, and z, τ Sx , τ Sy , τ Sz are the torque generated by the tether on the main star in x , components in the three directions of y and z.
步骤三、根据监测的速度情况,通过推进器对母星的位置进行控制,其中,根据目标与母星的速度差ΔV和调整速度差所期望的时间Δt,决定推进器的推力Fthru的大小。Step 3: Control the position of the mother star through the thruster according to the monitored speed. The thrust F thru of the thruster is determined based on the speed difference ΔV between the target and the mother star and the time Δt expected to adjust the speed difference. .
假设主星在成功抓捕目标后,主星与目标沿系绳方向的相对速度为零。而目标在旋转过程中势必会牵拉系绳使系绳中产生拉力。该拉力将导致主星沿系绳方向的加速度减小,目标沿系绳方向加速度的增加,如此目标沿系绳方向的速度势必大于主星。务必对该构型下母星施加控制方能够防止目标与主星相撞。Assume that after the main star successfully captures the target, the relative speed between the main star and the target along the tether direction is zero. The target will inevitably pull the tether during the rotation process, causing tension in the tether. This pulling force will cause the main star's acceleration along the tether direction to decrease and the target's acceleration along the tether direction to increase, so the target's speed along the tether direction must be greater than the main star's. It is necessary to control the home star in this configuration to prevent the target from colliding with the host star.
设计如下控制条件和控制律:Design the following control conditions and control laws:
控制条件:vT>vS,即目标的速度大于母星的速度Control condition: v T > v S , that is, the speed of the target is greater than the speed of the parent star
控制律:Control law:
ΔV=[vTx-vSx,vTy-vSy,vTz-vSz]ΔV=[v Tx -v Sx ,v Ty -v Sy ,v Tz -v Sz ]
目标与母星的速度差表示为:ΔV=[vTx-vSx,vTy-vSy,vTz-vSz],其中vTx,vTy,vTz分别表示目标速度在x,y,z三个方向的速度分量,vSx,vSy,vSz分别表示母星速度在x,y,z三个方向的速度分量;The speed difference between the target and the parent star is expressed as: ΔV = [v Tx -v Sx , v Ty - v Sy , v Tz - v Sz ], where v Tx , v Ty , and v Tz respectively represent the target speed in x, y, The velocity components in the three directions of z, v Sx , v Sy , and v Sz respectively represent the velocity components of the parent star's velocity in the three directions of x, y, and z;
所述推进器的推力:The thrust of the propeller:
其中,mS为母星的质量。Among them, m S is the mass of the parent star.
在一具体实施例中,如图1所示,认为目标为对称立方体,受到的拉力也是对称的,系绳的子绳有4根,其对称连接于目标上,这时目标受到的来自系绳的拉力可以在z方向上自动抵消,由于目标的往复摆动,在y方向也仅需少量微调即可。因此,只考虑目标与母星沿系绳方向即x轴方向的速度分量时,ΔV=[vTx-vSx,0,0],此时,只需打开与x轴方向对应位置的推进器即可。In a specific embodiment, as shown in Figure 1, the target is considered to be a symmetrical cube, and the tensile force it receives is also symmetrical. There are four sub-ropes of the tether, which are symmetrically connected to the target. At this time, the target receives the force from the tether. The pulling force can be automatically offset in the z-direction, and due to the reciprocating swing of the target, only a small amount of fine-tuning is required in the y-direction. Therefore, when only considering the velocity component of the target and the parent star along the tether direction, that is, the x-axis direction, ΔV = [v Tx -v Sx ,0,0]. At this time, you only need to turn on the thruster corresponding to the x-axis direction That’s it.
此时,控制条件为vTx>vSx,At this time, the control condition is v Tx > v Sx ,
控制律为:The control law is:
ΔV=[vTx-vSx,0,0],ΔV=[v Tx -v Sx ,0,0],
此构态下系绳处于松弛状态,而目标角速度不为零,消旋过程尚未结束。此控制律通过计算母星与目标的速度差,期望母星在Δt时间内达到目标沿x轴的线速度,由此得出所需推进器推力。在该推力作用下系绳将被拉紧并能够对目标施加消旋力矩,进而对目标进行消旋。In this configuration, the tether is in a relaxed state, but the target angular velocity is not zero, and the racemization process has not yet ended. This control law calculates the speed difference between the parent star and the target, and expects the parent star to reach the target's linear velocity along the x-axis within Δt time, thereby obtaining the required propeller thrust. Under the action of this thrust, the tether will be tightened and can exert a despin torque on the target, thereby despinating the target.
如图2所示,为消旋过程目标的摆动状态,其中小立方体表示目标,大立方体表示母星,其之间是系绳。图2(a)表示初始状态,目标没有摆动,系绳处于拉紧状态。图2(b)表示目标沿着顺时针方向发生了摆动,系绳松弛,此时目标的速度大于母星的速度,需要触发推进器,产生一定推力拉直系绳,拉直的系绳产生拉力,目标则会在该拉力作用下使其转速降为零并反向逆时针旋转。图2(c)表示母星对目标沿逆时针旋转时进行消旋,当目标的转速经过上一步消旋处理后沿逆时针旋转,目标速度再次超越母星速度时,推进器则根据控制率触发并产生对应的推进力,对目标进行二次消旋处理,如此反复,直至目标的转速降低至一定范围且目标与主星的速度差距缩小至0或接近0的状态,系统恢复稳定。As shown in Figure 2, it is the swing state of the target during the racemization process, where the small cube represents the target, the large cube represents the parent star, and the tether is in between. Figure 2(a) shows the initial state, the target does not swing, and the tether is in a tight state. Figure 2(b) shows that the target swings in the clockwise direction and the tether is relaxed. At this time, the speed of the target is greater than the speed of the parent star. The thruster needs to be triggered to generate a certain thrust to straighten the tether. The straightened tether generates pulling force. , the target will reduce its rotation speed to zero and rotate counterclockwise under the action of this pulling force. Figure 2(c) shows that the parent star decelerates the target when it rotates counterclockwise. When the target's speed rotates counterclockwise after the previous deracination process, and the target speed exceeds the speed of the parent star again, the thruster rotates according to the control rate. It triggers and generates corresponding propulsion force, and performs a second racecycling process on the target. This is repeated until the target's rotational speed reduces to a certain range and the speed difference between the target and the main star is reduced to 0 or close to 0, and the system returns to stability.
根据表1的仿真参数对目标单轴旋转消旋进行仿真。The target single-axis rotation racemization is simulated according to the simulation parameters in Table 1.
表1仿真参数Table 1 Simulation parameters
目标的初始角速度为绕z轴0.3rad/s,通过消旋仿真得到图3所示的目标角速度变化过程图,其中x,y方向的角速度没有发生变化,两条线重合为0,z方向的角速度每次穿越直线x=0表示触发了推进器并使得目标角速度方向进行切换,由图3可知,在150秒的时间内进行了四次消旋动作,目标每经历一次消旋动作,其角速度降低一点,直至降至0附近,系统恢复至稳定状态。由此可见,系统仅需要4次推进器的触发动作即可以将目标从0.3rad/s降至0附近,极大的节省了燃料。The initial angular velocity of the target is 0.3rad/s around the z-axis. The target angular velocity change process diagram shown in Figure 3 is obtained through derotation simulation, in which the angular velocity in the x and y directions does not change, and the two lines coincide with 0, and the z direction Each time the angular velocity crosses the straight line x=0, it means that the thruster is triggered and the direction of the target's angular velocity is switched. As can be seen from Figure 3, four deceleration actions are performed within 150 seconds. Each time the target undergoes a deceleration action, its angular velocity Lower it a little until it drops to near 0, and the system returns to a stable state. It can be seen that the system only needs 4 thruster triggering actions to reduce the target from 0.3rad/s to near 0, which greatly saves fuel.
本发明的一种空间碎片消旋控制方法,是基于由母星、空间碎片和连接于母星与空间碎片之间的系绳组成的组合体即绳系系统进行的。如图1所示,绳系系统还包括系绳弹射装置和弹簧阻尼单元;系绳弹射装置安装于主星上,系绳和弹簧阻尼单元贮存在系绳弹射装置的内部,当需要捕获空间碎片时,触发系绳弹射装置,系绳和弹簧阻尼单元一并弹出。A space debris despin control method of the present invention is based on a combination of a parent star, space debris, and a tether connected between the parent star and the space debris, that is, a tether system. As shown in Figure 1, the tether system also includes a tether ejection device and a spring damping unit; the tether ejection device is installed on the main satellite, and the tether and spring damping unit are stored inside the tether ejection device when space debris needs to be captured. , triggering the tether ejection device, and the tether and spring damping unit eject together.
系绳包括主绳和子绳,主绳与主星相连,子绳与空间碎片相连;弹簧阻尼单元设置于主绳和子绳上。The tether includes a main rope and a sub-tether, the main rope is connected to the main star, and the sub-tether is connected to the space debris; the spring damping unit is provided on the main rope and the sub-tether.
子绳的数量为n,n条子绳与主绳连接在同一点P点;其中n为大于等于2的正整数;如图1所示,子绳的数量为4条,子绳的数量与空间碎片的形状有关,应该尽量保证空间碎片收到子绳的拉力均匀且对称。The number of sub-ropes is n, and n sub-ropes are connected to the same point P as the main rope; where n is a positive integer greater than or equal to 2; as shown in Figure 1, the number of sub-ropes is 4, and the number of sub-ropes is related to the space It depends on the shape of the debris. Try to ensure that the pulling force received by the space debris from the sub-rope is even and symmetrical.
主绳上设置有多个弹簧阻尼单元,每条子绳上均设置有一个弹簧阻尼单元。在具体实施过程中,如果需要也可以在子绳上设置多个弹簧阻尼单元。如图1所示,主绳上设置有3个弹簧阻尼单元,在具体实施过程中,也可以设置任意整数个。Multiple spring damping units are provided on the main rope, and one spring damping unit is provided on each sub-rope. In the specific implementation process, if necessary, multiple spring damping units can also be provided on the sub-ropes. As shown in Figure 1, three spring damping units are provided on the main rope. During the specific implementation process, any integer number can also be provided.
以上的具体实施例仅描述了本发明的设计原理,该描述中的部件形状,名称可以不同,不受限制。所以,本发明领域的技术人员可以对前述实施例记载的技术方案进行修改或等同替换;而这些修改和替换未脱离本发明创造宗旨和技术方案,均应属于本发明的保护范围。The above specific embodiments only describe the design principles of the present invention. The shapes and names of the components in this description can be different and are not limited. Therefore, those skilled in the field of the present invention can make modifications or equivalent substitutions to the technical solutions described in the foregoing embodiments; and these modifications and substitutions do not deviate from the creative purpose and technical solutions of the present invention, and should all fall within the protection scope of the present invention.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111066585.3A CN113772127B (en) | 2021-09-13 | 2021-09-13 | Space debris racemization control method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111066585.3A CN113772127B (en) | 2021-09-13 | 2021-09-13 | Space debris racemization control method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113772127A CN113772127A (en) | 2021-12-10 |
CN113772127B true CN113772127B (en) | 2023-12-08 |
Family
ID=78842695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111066585.3A Active CN113772127B (en) | 2021-09-13 | 2021-09-13 | Space debris racemization control method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113772127B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115416879B (en) * | 2022-09-14 | 2024-05-10 | 北京理工大学 | Device and method for speed reduction racemization of space nonmagnetic metal body |
CN118168561B (en) * | 2024-05-14 | 2024-08-09 | 北京航天驭星科技有限公司 | Satellite speed acquisition method and system, satellite track lifting method and control system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6116544A (en) * | 1997-09-12 | 2000-09-12 | Tethers Unlimited, Inc. | Electrodynamic tether and method of use |
CN106114919A (en) * | 2016-08-01 | 2016-11-16 | 北京理工大学 | A kind of space junk rope system pulls racemization and method for cleaning |
CN107364589A (en) * | 2017-07-04 | 2017-11-21 | 上海宇航系统工程研究所 | Racemization control method of being diversion based on more tether tie points to Spatial Instability target |
CN107643689A (en) * | 2017-10-19 | 2018-01-30 | 北京理工大学 | A kind of rope system towing stable control method of space junk |
CN109319171A (en) * | 2018-10-19 | 2019-02-12 | 北京航空航天大学 | A method for lateral angular velocity suppression and spin direction control of space debris |
CN109987258A (en) * | 2019-01-28 | 2019-07-09 | 西北工业大学深圳研究院 | A racemization method for space robots after capturing non-cooperative targets |
CN112597587A (en) * | 2020-12-23 | 2021-04-02 | 北京理工大学 | Swing suppression method for recovery of failure satellite tether |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8052092B2 (en) * | 2009-01-30 | 2011-11-08 | The Boeing Company | Method and apparatus for satellite orbital change using space debris |
-
2021
- 2021-09-13 CN CN202111066585.3A patent/CN113772127B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6116544A (en) * | 1997-09-12 | 2000-09-12 | Tethers Unlimited, Inc. | Electrodynamic tether and method of use |
CN106114919A (en) * | 2016-08-01 | 2016-11-16 | 北京理工大学 | A kind of space junk rope system pulls racemization and method for cleaning |
CN107364589A (en) * | 2017-07-04 | 2017-11-21 | 上海宇航系统工程研究所 | Racemization control method of being diversion based on more tether tie points to Spatial Instability target |
CN107643689A (en) * | 2017-10-19 | 2018-01-30 | 北京理工大学 | A kind of rope system towing stable control method of space junk |
CN109319171A (en) * | 2018-10-19 | 2019-02-12 | 北京航空航天大学 | A method for lateral angular velocity suppression and spin direction control of space debris |
CN109987258A (en) * | 2019-01-28 | 2019-07-09 | 西北工业大学深圳研究院 | A racemization method for space robots after capturing non-cooperative targets |
CN112597587A (en) * | 2020-12-23 | 2021-04-02 | 北京理工大学 | Swing suppression method for recovery of failure satellite tether |
Also Published As
Publication number | Publication date |
---|---|
CN113772127A (en) | 2021-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107364589B (en) | Flying around and despinning control method for space instability target based on multi-rope connecting points | |
CN112597587B (en) | Swing suppression method for recovery of failure satellite tether | |
CN113772127B (en) | Space debris racemization control method | |
CN109375648B (en) | Elliptical orbit satellite formation configuration initialization method under multi-constraint condition | |
CN106114919B (en) | A kind of space junk rope system pulls racemization and method for cleaning | |
CN107145071B (en) | Lagrange power model and controller of space tether system | |
CN108181807B (en) | A kind of satellite initial state stage self-adapted tolerance attitude control method | |
CN104570742B (en) | Feedforward PID (proportion, integration and differentiation) control based rapid high-precision relative pointing control method of noncoplanar rendezvous orbit | |
O’Connor et al. | Detumbling of space debris by a net and elastic tether | |
CN106774360B (en) | A target star attitude stabilization method using tethers and connecting rods in towed orbit change | |
CN106272380B (en) | A kind of attitude stabilization method for arresting mechanical arm assembly after high speed rotation target | |
CN111552326A (en) | Space debris rope system towing off-orbit target attitude nutation suppression control method and system | |
Safaei et al. | Lyapunov-based nonlinear controller for quadrotor position and attitude tracking with GA optimization | |
Ishijima et al. | The on-orbit maneuvering of large space flexible structures by free-flying robots | |
CN112362289A (en) | Aircraft split test method and device, computer equipment and readable storage medium | |
Wang et al. | Modeling and predefined-time anti-disturbance control for the aerial refueling phase of receiver aircraft | |
CN106020217A (en) | Reel-controlled towing orbital transfer anti-winding and anti-collision method | |
CN114706413B (en) | Attitude control method and system for low-earth orbit micro-nano satellites with variable center of mass | |
CN111392069A (en) | A tethered drag derotation and deorbit method for a failed liquid-filled spacecraft | |
CN109018442B (en) | Novel low-cost satellite three-axis attitude time-sharing decoupling high-multiplexing air injection control method | |
CN113607378B (en) | A method for simulating and suppressing forced free angular motion of a tethered aircraft model | |
CN110562500A (en) | non-cooperative target three-dimensional rolling motion spinning simulation air injection control method and system | |
CN111736459B (en) | Tether vibration rapid suppression control method without initial value dependence | |
Khoo et al. | Robust control of novel thrust vectored 3D printed multicopter | |
CN105799949A (en) | Pressure center design method, attitude control method and system of suborbital satellite |
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 |