基于Krylov子空间法的柔性航天器降阶研究

带大型柔性附件的复杂航天器,其弹性变形通常用有限单元法描述.由于有限元节点坐标数目庞大将会给动力学方程求解和控制策略的实现带来巨大负担,因此对柔性太阳电池阵的动力学模型降阶进行研究.利用单向递推组集推导带太阳电池阵的航天器的刚柔耦合动力学方程,分别采用模态截断法、模态价值分析法和Krylov子空间法对太阳电池阵进行模型降阶,比较不同降阶模型和全有限元模型的动力学仿真结果.仿真结果表明,无论是否考虑大范围运动的影响,采用Krylov方法只需要较低的自由度就可以得到和采用有限元方法完全一致的结果.说明Krylov方法能够有效地降低航天器柔性附件的自由度,提高动力学仿真的效率,便于驱动控制的实现.     

基于动力学递推算法的绳系卫星系统刚柔耦合多体模型

利用动力学递推算法建立了空间绳系卫星系统(TSS)的铰接柔性杆动力学模型.所研究的绳系系统包含两颗刚体绳端卫星和一段连接两颗卫星的柔绳,柔绳离散为一系列球铰连接的弹性杆,考虑杆不均匀的纵向一维形变,并利用球铰的运动模拟柔绳的弯曲和扭转;然后基于递推算法推导得出了该刚柔混合离散模型的动力学方程.数值仿真结果表明:所建立的绳系卫星系统模型能够很好地模拟空间系绳的各向摆振和扭转,此外,递推算法的引入能够有效地减少离散模型动力学方程的维数,从而大幅减少计算量.      

Comparison of ionosphere characteristic parameters obtained by ionosonde with IRI-2007 model over Southeast Asia

In this work, the foF2 and hmF2 parameters at the conjugate points near the magnetic equator of Southeast Asia are studied and compared with the International Reference Ionosphere (IRI) model. Three ionosondes are installed nearly along the magnetic meridian of 100°E; one at the magnetic equator, namely Chumphon (10.72°N, 99.37°E, dip angle 3.0°N), and the other two at the magnetic conjugate points, namely Chiang Mai (18.76°N, 98.93°E, dip angle 12.7°N) and Kototabang (0.2°S, 100.30°E, dip angle 10.1°S). The monthly hourly medians of the foF2 and hmF2 parameters are calculated and compared with the predictions obtained from the IRI-2007 model from January 2004 to February 2007. Our results show that: the variations of foF2 and hmF2 predicted by the IRI-2007 model generally show the similar feature to the observed data. Both parameters generally show better agreement with the IRI predictions during daytime than during nighttime. For foF2, most of the results show that the IRI model overestimates the observed foF2 at the magnetic equator (Chumphon), underestimates at the northern crest (Chiang Mai) and is close to the measured ones at the southern crest of the EIA (Kototabang). For hmF2, the predicted hmF2 values are close to the hmF2(M3000F2_OBS) during daytime. During nighttime, the IRI model gives the underestimation at the magnetic equator and the overestimation at both EIA crests. The results are important for the future improvements of the IRI model for foF2 and hmF2 over Southeast Asia region.     

Multi-Sphere Method for modeling spacecraft electrostatic forces and torques

The use of electrostatic (Coulomb) actuation for formation flying is attractive because non-renewable fuel reserves are not depleted and plume impingement issues are avoided. Prior analytical electrostatic force models used for Coulomb formations assume spherical spacecraft shapes, which include mutual capacitance and induced effects. However, this framework does not capture any orientation-dependent forces or torques on generic spacecraft geometries encountered during very close operations and docking scenarios. The Multi-Sphere Method (MSM) uses a collection of finite spheres to represent a complex shape and analytically approximate the Coulomb interaction with other charged bodies. Finite element analysis software is used as a truth model to determine the optimal sphere locations and radii. The model is robust to varying system parameters such as prescribed voltages and external shape size. Using the MSM, faster-than-realtime electrostatic simulation of six degree of freedom relative spacecraft motion is feasible, which is crucial for the development of robust relative position and orientation control algorithms in local space situational awareness applications. To demonstrate this ability, the rotation of a cylindrical craft in deep space is simulated, while charge control from a neighboring spacecraft is used to de-spin the object. Using a 1 m diameter craft separated by 10 m from a 3 by 1 m cylindrical craft in deep space, a 2 °/s initial rotation rate can be removed from the cylinder within 3 days, using electric potentials up to 30 kV.     

Finite element reconstruction approach for on-orbit spacecraft breakup dynamics simulation and fragment analysis

Breakup model is the key area of space debris environment modeling. NASA standard breakup model is currently the most widely used for general-purpose. It is a statistical model found based on space surveillance data and a few ground-based test data. NASA model takes the mass, impact velocity magnitude for input and provides the fragment size, area-to-mass ratio, velocity magnitude distributions for output. A more precise approach for spacecraft disintegration fragment analysis is presented in this paper. This approach is based on hypervelocity impact dynamics and takes the shape, material, internal structure and impact location etc. of spacecraft and impactor, which might greatly affect the fragment distribution, into consideration. The approach is a combination of finite element and particle methods, entitled finite element reconstruction (FER). By reconstructing elements from the particle debris cloud, reliable individual fragments are identified. Fragment distribution is generated with undirected graph conversion and connected component analysis. Ground-based test from literature is introduced for verification. In the simulation satellite targets and impactors are modeled in detail including the shape, material, internal structure and so on. FER output includes the total number of fragments and the mass, size and velocity vector of each fragment. The reported fragment distribution of FER shows good agreement with the test, and has good accuracy for small fragments.     

Quality variation of GPS satellite clocks on-orbit using IGS clock products

Over 60% clocks on board of the GPS satellites are working longer than their designed life. Therefore realizing their stabilities in a long time scales is essential to GPS navigation and positioning plus IGS time scale maintaining. IGS clock products from 2001 to 2010 are used to analyze the GPS satellite clock qualities such as frequency stabilities and clock noise level. We find out that for the clocks of Block IIA satellites the frequency stabilities and clock noise are 10 times worse than that of the Block IIR and IIR-M satellites. Moreover, the linear relationships between frequency stabilities and clock residuals have been deduced with an accuracy of better than 0.02 ns. Specially, it is noticed that the clock of the PRN27 is instable and the relationship between the frequency stability and residuals is at least a quadratic curve. Therefore, we suggested that GPS satellite clocks should be weighted by their quality levels in application, and the observations of the Block IIA should not be used for real-time positioning which required precision better than one meter.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

An improved approach to model regional ionosphere and accelerate convergence for precise point positioning

Given the severe effects of the ionosphere on global navigation satellite system (GNSS) signals, single-frequency (SF) precise point positioning (PPP) users can only achieve decimeter-level positioning results. Ionosphere-free combinations can eliminate the majority of ionospheric delay, but increase observation noise and slow down dual-frequency (DF) PPP convergence. In this paper, we develop a regional ionosphere modeling and rapid convergence approach to improve SF PPP (SFPPP) accuracy and accelerate DF PPP (DFPPP) convergence speed. Instead of area model, ionospheric delay is modeled for each satellite to be used as a priori correction. With the ionospheric, wide-lane uncalibrated phase delay (UPD) and residuals satellite DCBs product, the wide-lane observations for DF users change to be high-precision pseudorange observations. The validation of a continuously operating reference station (CORS) network was analyzed. The experimental results confirm that the approach considerably improves the accuracy of SFPPP. For DF users, convergence time is substantially reduced.     

被动型星载氢钟H型气体电离模型及参数优化

针对被动型星载氢钟气体电离装置低功耗、低电磁干扰的设计要求,设计了一种基于平面盘绕馈能天线的H型电离系统,并建立该电离方法的数学模型。首先,由Maxwell方程出发,推导出了电离泡中的电场和磁场的分布模型,得到气体电离机制;再由电磁场的分布方程推导出氢气击穿判据表达式,建立了气体电离条件的数学模型;最后,利用电离击穿条件判据方程优化了电离装置的设计参数。数值仿真与测试表明该模型是有效的,采用优化的设计参数可以使电离功耗和电磁干扰特性得到改善,且电离直流电源总功率可小于2W。     

基于SSI模型的加速应力试验定量评估方法

为了解决传统应力-强度干涉模型(SSI,Stress-Strength Interference)不能对加速应力试验(AST,Accelerated Stress Testing)进行可靠性指标随时间变化的定量评估问题,针对应力和强度均引入了随机过程进行描述.首先假设了产品强度退化过程和应力变化分别为漂移布朗运动和对称布朗运动.其次,为了解决Basu的布朗运动SSI模型在零时刻可靠度为1的问题,从工程角度,假设产品设计制造过程与试验或使用过程相互独立,提出了产品可靠性是设计生产决定的固有可靠性以及试验使用过程确定的使用可靠性的串联结果.通过将传统应力强度干涉模型结合Basu的布朗运动SSI模型,建立了与时间相关的应力强度干涉模型IMSDⅠ,并给出了相关评估方法.基于将加速应力试验与加速退化试验相结合,改进了已有的AST试验方法,扩大了试验数据量,最后,通过仿真算例验证了本文模型和评估方法的理论正确性和工程可实施性.