含固定转速驱动电机的挠性卫星自旋稳定性分析

分析一类含驱动电机挠性卫星的自旋稳定性,应用不变集定理证明了系统轨迹收敛于系统平衡点集,并通过构造正定的单调递减的类能量函数来判断各平衡点的稳定性质。获得的结论包括：由于挠性振动阻尼特性,在无外扰力矩作用时,含固定转速驱动电机的挠性卫星存在稳定的自旋;不含驱动电机的挠性卫星,稳定的自旋轴为其最大主惯量轴;驱动电机可改变挠性卫星的自旋轴,给出了稳定的自旋轴与电机转速方向夹角公式,该夹角为锐角;当电机转动角动量幅值远小于系统角动量时,稳定自旋轴接近最大主惯量轴;电机转动角动量幅值越接近系统角动量,稳定自旋轴越接近电机转动方向。文中也给出了稳定自旋轴平行电机转轴的条件。     

On the stability of spinning satellites

We study the directional stability of rigid and deformable spinning satellites in terms of two attitude angles. The linearized attitude motion of a free system about an assumed uniform-spin reference solution leads to a generic MGK system when the satellite is rigid or deformable. In terms of Lyapunov’s stability theory, we investigate the stability with respect to a subset of the variables. For a rigid body, the MGK system is 6-dimensional, i.e., 3 rotational and 3 translational variables. When flexible parts are present the system can have any arbitrary dimension. The 2×2 McIntyre–Myiagi stability matrix gives sufficient conditions for the attitude stability. A further development of this method has led to the Equivalent Rigid Body method. We propose an alternative practical method to establish sufficiency conditions for directional stability by using the Frobenius–Schur reduction formula. As practical applications we discuss a spinning satellite augmented with a spring–mass system and a rigid body appended with two cables and tip masses. In practice, the attitude stability must also be investigated when the spinning satellite is subject to a constant axial thrust. The generic format becomes MGKN as the thrust is a follower force. For a perfectly aligned thrust along the spin axis, Lyapunov’s indirect method remains valid also when deformable parts are present. We illustrate this case with an apogee motor burn in the presence of slag. When the thrust is not on the spin axis or not pointing parallel to the spin axis, the uniform-spin reference motion does not exist and none of the previous methods is applicable. In this case, the linearization may be performed about the initial state. Even when the linearized system has bounded solutions, the non-linear system can be unstable in general. We illustrate this situation by an instability that actually happened in-flight during a station-keeping maneuver of ESA’s GEOS-I satellite in 1979.     

Variation of the Satellite Orbit Altitude by the Light Pressure Force

The possibility of the uncontrolled increase of the altitude of an almost circular satellite orbit by the force of the light pressure is investigated. The satellite is equipped with a damper and a system of mirrors (solar batteries can serve as such a system). The flight of the satellite takes place in the mode of a single-axis gravitational orientation, the axis of its minimum principal central moment of inertia makes a small angle with the local vertical and the motion of the satellite around this axis constitutes forced oscillations under the impact of the moment of force of the light pressure. The form of the oscillations and the initial orbit are chosen so that the transverse component of the force of the light pressure acting upon the satellite be positive and the semimajor axis of the orbit would continuously increase. As this takes place, the orbit remains almost circular. We investigate the evolution of the orbit over an extended time interval by the method which employs separate integration of the equations of the orbital and rotational motions of the satellite. The method includes outer and inner cycles. The outer cycle involves the numerical integration of the averaged equations of motion of the satellite center of mass. The inner cycle serves to calculate the right-hand sides of these equations. It amounts to constructing an asymptotically stable periodic motion of the satellite in the mode of a single-axis gravitational orientation for current values of the orbit elements and to averaging the equations of the orbital motion along it. It is demonstrated that the monotone increase of the semimajor axis takes place during the first 15 years of motion. In actuality, the semimajor axis oscillates with a period of about 60 years. The eccentricity and inclination of the orbit remain close to their initial values.     

地球同步卫星远地点变轨期间的三轴姿态稳定

本文提供了地球同步三轴稳定卫星在远地点变轨机动过程中的卫星动力学模型,其中考虑了刚体、液体晃动和帆板振动的耦合,这些模型经过实际工程验证,证明是正确的。本文还给出了控制器结构,它由PID和各种晃动滤波器及帆板挠性滤波器组成。工程设计证明,这种控制器结构可满足设计要求且简单可靠。本文还给出了系统稳定性分析结果及数字仿真结果。     

The use of refined model of aerodynamic moment in the problem of reconstruction of the <Emphasis Type="Italic">Foton</Emphasis> satellite rotational motion

A new mathematical model of the uncontrolled rotational motion of the Foton satellite is presented. The model is based on the Euler dynamic equations of rigid body motion and takes into account the action upon the satellite of four external mechanical moments: gravitational, restoring aerodynamic, moment with constant components in the satellite-fixed coordinate system, and moment arising due to interaction of the Earth’s magnetic field with the satellite’s proper magnetic moment. To calculate the aerodynamic moment a special geometrical model of the outer satellite shell is used. Detailed form of the formulas giving above-mentioned moments in the equations of satellite motion is agreed with the form of the considered motion. Model testing is performed by determining with its help the rotational motion of the Foton M-2 satellite (it was in orbit from May 31, 2005 to June 16, 2005) using the data of the onboard measurements of the Earth’s magnetic field strength. The use of the new model has led to a relatively small improvement in the accuracy of the motion determination, but allowed us to obtain physically real estimates of some parameters.     

柔性多体卫星自抗扰控制系统的研究

设计了一种新型的姿态控制器一姿态自抗扰控制器。采用简化了的刚体动力学模型，即把柔性太阳帆板与柔性天线等的振动都看作外干扰，对象的不确定性、未建模动态与外干扰的总量都由自抗扰控制器估计出，再用该估计总量实现动态反馈补偿。该姿态控制器，只需控制输入转换矩阵的估计值，无需精确的卫星动力学模型。改变星体转动惯量或设置干扰力矩情况下的数学仿真表明，所设计的自抗扰控制器具有较小的稳态误差、较快的动态响应及对干扰有较强的抑制能力。     

A heuristic design method for attitude stabilization of magnetic actuated satellites

In this paper a heuristic design strategy for stabilizing the satellite attitude has been proposed. It is assumed that the satellite is actuated by a set of mutually perpendicular magnetic coils. Using well-known Lyapunov direct stability method it is shown that the proposed controller causes to a global asymptotic stable system for all near polar orbits. The design procedure is based on analyzing of the conceptual effects of magnetic coils on the satellite attitude motion. Considering these effects lead to some intuitive results which determine the global stabilizing control law. The performance and robustness of the designed controller against actuators saturation and quantization error have been verified using a real-time-hardware–software in-loop (RTHSIL) simulation results. These results show that the global stability can be achieved although some disturbances and restrictions exist. This stabilizing controller can be simply combined with a linear explicit model predictive controller (EMPC) to achieve a full three-axis control law.     

Uncontrolled motion of the <Emphasis Type="Italic">Foton M-2</Emphasis> satellite and quasistatic microaccelerations on its board

The results of determining the uncontrolled rotational motion of the Foton M-2 satellite (in orbit from May 31 to June 16, 2005) are presented. The determination was made using the data of onboard measurements of the Earth’s magnetic field strength. Segments 270 min long (three orbits) were selected from these data covering the first two thirds of the flight. On each such segment the data were processed jointly by the least squares method using integration of the equations of attitude motion of the satellite. In processing, the initial conditions of motion and parameters of the used mathematical model were estimated. The thus obtained results gave a complete overview of the satellite motion. This motion, having started with a small angular velocity, gradually accelerated, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 09, 2005 (the last day of measurements) the angular velocity of the satellite relative to its lengthwise axis was about 1.1 deg/s, while the projection of the angular velocity onto a plane perpendicular to this axis had an absolute value of about 0.11 deg/s. Deviations of the lengthwise axis from a normal to the orbit plane did not exceed 60°. Based on the results of determination of the rotational motion of the satellite, calculations of quasi-static microaccelerations on its board are performed.     

A heuristic design method for attitude stabilization of magnetic actuated satellites

In this paper a heuristic design strategy for stabilizing the satellite attitude has been proposed. It is assumed that the satellite is actuated by a set of mutually perpendicular magnetic coils. Using well-known Lyapunov direct stability method it is shown that the proposed controller causes to a global asymptotic stable system for all near polar orbits. The design procedure is based on analyzing of the conceptual effects of magnetic coils on the satellite attitude motion. Considering these effects lead to some intuitive results which determine the global stabilizing control law. The performance and robustness of the designed controller against actuators saturation and quantization error have been verified using a real-time-hardware–software in-loop (RTHSIL) simulation results. These results show that the global stability can be achieved although some disturbances and restrictions exist. This stabilizing controller can be simply combined with a linear explicit model predictive controller (EMPC) to achieve a full three-axis control law.     

配合间隙对涡轮泵转子低速动平衡的影响研究

刚性涡轮泵转子低速动平衡是装配过程中的重要环节,残余不平衡量的大小是转子振动的关键影响因素之一,也是影响涡轮泵可靠性的重要因素。在某型发动机涡轮泵研制过程中,出现试车分解后不平衡量较装配时增加10~20倍的情况。对该涡轮泵转子建立了装配间隙对动平衡影响程度的理论模型,研究了不同零件的偏心距以及配合间隙对转子残余不平衡量的影响程度。根据影响程度对比,可以确定需要严格控制装配状态和单独动平衡的零件。通过反复拆装动平衡测试试验,对理论模型进行了验证。     

用于图像定位与配准的扫描辐射计扫描镜热变形模型研究

以某卫星扫描成像辐射计为研究对象,建立了扫描镜的热变形模型,并对模型进行了简化。模型应用结果表明,扫描镜镜面与扫描镜转轴需采用热特性接近的材料,通过增强仪器结构在Y向的刚度,能提高热变形模型的精度。用该模型获得的扫描镜热变形运动补偿(TMC)精度可满足卫星的使用要求。     

Analysis of Low-Frequency Microaccelerations onboard the Foton-11Satellite

The method and the results of investigating the low-frequency component of microaccelerations onboard the Foton-11satellite are presented. The investigation was based on the processing of data of the angular velocity measurements made by the German system QSAM, as well as the data of measurements of microaccelerations performed by the QSAM system and by the French accelerometer BETA. The processing was carried out in the following manner. A low-frequency (frequencies less than 0.01 Hz) component was selected from the data of measurements of each component of the angular velocity vector or of the microacceleration, and an approximation was constructed of the obtained vector function by a similar function that was calculated along the solutions to the differential equations of motion of the satellite with respect to its center of mass. The construction was carried out by the least squares method. The initial conditions of the satellite motion, its aerodynamic parameters, and constant biases in the measurement data were used as fitting parameters. The time intervals on which the approximation was constructed were from one to five hours long. The processing of the measurements performed with three different instruments produced sufficiently close results. It turned out to be that the rotational motion of the satellite during nearly the entire flight was close to the regular Eulerian precession of the axially symmetric rigid body. The angular velocity of the satellite with respect to its longitudinal axis was about 1 deg/s, while the projection of the angular velocity onto the plane perpendicular to this axis had an absolute value of about 0.2 deg/s. The magnitude of the quasistatic component of microaccelerations in the locations of the accelerometers QSAM and BETA did not exceed 5 × 10^–5–10^–4m/s²for the considered motion of the satellite.     

冗余空间机械臂抓捕自旋卫星后的消旋控制

旨在提出一种运动学冗余空间机器人抓捕自旋卫星后的消旋策略和协调控制方法。首先，给出运动学冗余空间机器人捕获目标后的动力学模型，作为协调控制器设计基础。然后，基于四阶Bézier曲线和满足特定约束的自适应微分进化(Differential Evolution, DE)算法提出抓捕后的最优消旋与路径规划策略，最优消旋策略中同时考虑对消旋时间和控制力矩进行优化。提出一种跟踪所设计参考轨迹的协调控制方法，调整基座的姿态达到期望值。所提方法有效地衰减了自旋卫星的初始角速度，同时实现对基座姿态的控制。文末给出利用7 DOF冗余空间机械臂消除目标自旋运动的仿真结果，表明所提方法的有效性。     

A hybrid attitude controller consisting of electromagnetic torque rods and an active fluid ring

In this paper, a novel hybrid actuation system for satellite attitude stabilization is proposed along with its feasibility analysis. The system considered consists of two magnetic torque rods and one fluid ring to produce the control torque required in the direction in which magnetic torque rods cannot produce torque. A mathematical model of the system dynamics is derived first. Then a controller is developed to stabilize the attitude angles of a satellite equipped with the abovementioned set of actuators. The effect of failure of the fluid ring or a magnetic torque rod is examined as well. It is noted that the case of failure of the magnetic torque rod whose torque is along the pitch axis is the most critical, since the coupling between the roll or yaw motion and the pitch motion is quite weak. The simulation results show that the control system proposed is quite fault tolerant.     

Estimating the accuracy of the technique of reconstructing the rotational motion of a satellite based on the measurements of its angular velocity and the magnetic field of the Earth

The paper has studied the accuracy of the technique that allows the rotational motion of the Earth artificial satellites (AES) to be reconstructed based on the data of onboard measurements of angular velocity vectors and the strength of the Earth magnetic field (EMF). The technique is based on kinematic equations of the rotational motion of a rigid body. Both types of measurement data collected over some time interval have been processed jointly. The angular velocity measurements have been approximated using convenient formulas, which are substituted into the kinematic differential equations for the quaternion that specifies the transition from the body-fixed coordinate system of a satellite to the inertial coordinate system. Thus obtained equations represent a kinematic model of the rotational motion of a satellite. The solution of these equations, which approximate real motion, has been found by the least-square method from the condition of best fitting between the data of measurements of the EMF strength vector and its calculated values. The accuracy of the technique has been estimated by processing the data obtained from the board of the service module of the International Space Station (ISS). The reconstruction of station motion using the aforementioned technique has been compared with the telemetry data on the actual motion of the station. The technique has allowed us to reconstruct the station motion in the orbital orientation mode with a maximum error less than 0.6° and the turns with a maximal error of less than 1.2°.     

基于T-H方程的卫星轨迹模型参考输出跟踪控制方法

当目标卫星沿椭圆轨道运行时,描述追踪星与目标星相对运动的线性化方程为T-H方程。将描述航天器相对运动的T-H方程变换为周期系统的状态空间形式,并给出卫星轨迹跟踪控制问题的数学描述。基于周期系统的参量Lyapunov方法和模型参考跟踪控制理论,提出了卫星轨迹跟踪控制器的设计方法。利用该方法设计了带有收敛速率保障的反馈镇定控制器和具有自由参数的前馈控制器。对追踪星相对目标星悬停任务进行了数值仿真,仿真结果表明提出的控制方案是有效的。     

航天器机动时DGMSCMG磁悬浮转子干扰补偿控制

双框架磁悬浮控制力矩陀螺（DGMSCMG）具有寿命长、综合效益好等突出优势,但航天器机动时,航天器及DGMSCMG内、外框架系统的转动均导致磁悬浮高速转子产生一定的耦合运动,影响磁悬浮转子系统的稳定性,同时使输出力矩精度下降,从而严重影响航天器姿态控制的精度。本文建立了基于DGMSCMG的航天器动力学模型,分析航天器、外框架、内框架、磁悬浮转子四者之间的动力学耦合关系。针对磁悬浮转子的非线性耦合干扰,提出一种基于复合控制的补偿方法,通过磁轴承产生相应的电磁力,对陀螺耦合力矩和惯性耦合力矩进行补偿控制。仿真结果表明,干扰补偿控制能有效抑制航天器及框架对磁悬浮转子的耦合干扰,也有效提高了磁悬浮转子系统的稳定性。     

卫星编队飞行相对位置自适应协同控制

针对卫星编队飞行相对位置协同控制问题,基于编队卫星相对运动非线性动力学方程和一致性理论设计了两种自适应协同控制器。首先,在卫星质量不确定和星间信息交互存在通信时延的条件下,设计了一种全状态反馈自适应协同控制器,并证明了该控制策略对空间摄动力的鲁棒性。其次,进一步考虑速度信息不可测的条件下,采用滤波器设计了一种无速度反馈的自适应协同控制器。最后,以编队构型重构为例对两种自适应协同控制器进行了仿真校验。仿真结果表明：两种自适应协同控制器均可有效应用于卫星编队飞行相对位置的协同控制,能够保证编队卫星对各自期望轨迹跟踪的同时暂态保持编队构型的稳定,具有较高的控制精度。     

A novel single thruster control strategy for spacecraft attitude stabilization

Feasibility of achieving three axis attitude stabilization using a single thruster is explored in this paper. Torques are generated using a thruster orientation mechanism with which the thrust vector can be tilted on a two axis gimbal. A robust nonlinear control scheme is developed based on the nonlinear kinematic and dynamic equations of motion of a rigid body spacecraft in the presence of gravity gradient torque and external disturbances. The spacecraft, controlled using the proposed concept, constitutes an underactuated system (a system with fewer independent control inputs than degrees of freedom) with nonlinear dynamics. Moreover, using thruster gimbal angles as control inputs make the system non-affine (control terms appear nonlinearly in the state equation). This necessitates the control algorithms to be developed based on nonlinear control theory since linear control methods are not directly applicable. The stability conditions for the spacecraft attitude motion for robustness against uncertainties and disturbances are derived to establish the regions of asymptotic 3-axis attitude stabilization. Several numerical simulations are presented to demonstrate the efficacy of the proposed controller and validate the theoretical results. The control algorithm is shown to compensate for time-varying external disturbances including solar radiation pressure, aerodynamic forces, and magnetic disturbances; and uncertainties in the spacecraft inertia parameters. The numerical results also establish the robustness of the proposed control scheme to negate disturbances caused by orbit eccentricity.