天基激光移除空间碎片仿真平台研究与开发

为了从空间碎片和天基移除系统两方面分析激光驱动碎片变轨过程,优化移除任务规划和策略,基于真实的可观测空间碎片数据,设计并开发了一套天基激光移除空间碎片三维数值仿真平台。首先从总体设计出发,对三维数值仿真平台的需求分析、总体框架、模块功能进行了明确的描述。其次通过对激光驱动空间碎片变轨过程数学模型的分析,确定了各模块的具体实现方法。最后采用C++/Qt开发了三维数值仿真平台,通过仿真验证了设计平台的有效性。该仿真平台可用于不同天基平台和目标碎片的任务规划、碎片分布热点区域和航天器防护区域的方案设计及空间环境治理体系的优化设计。     

J2轨道模型对月地激光通信预瞄准精度影响研究

基于J_2轨道预测模型,设计了不同倾角和轨道高度圆轨道月球卫星,通过将J_2轨道预测模型预瞄准仿真结果与直接积分RKF7(8)法仿真结果相对比,研究了不同类型月球卫星轨道对-Y面卫星舱板激光通信终端瞄准精度的影响。仿真结果表明,J_2轨道预测模型可以满足月球极地卫星月地激光通信要求,当轨道高度为1 000 km和2 000 km的时候,10 min方位角偏差不超过40μrad,而俯仰角偏差仅为7μrad时,在一定程度上J2轨道预测模型可以满足月地激光通信预瞄准要求。     

基于压电纤维复合材料的航天器动力学建模与振动抑制

压电纤维复合材料(MFC)在柔性航天器的振动主动抑制中具有很好的应用前景。利用哈密顿原理和压电驱动的载荷比拟方法,建立了带MFC压电驱动的离散形式的刚柔耦合动力学方程,采用线性二次型最优控制(LQR)算法进行主动控制。结果表明：在航天器的柔性体受到脉冲载荷激励条件下,使用MFC驱动器可以实现航天器挠性振动的快速抑制,并且同时保持中心刚体姿态的稳定性,即能够实现挠性振动与姿态运动的协同控制。基于MFC的主动控制方法对于高频响应也具有较好的控制效果。对于柔性占优的航天器,采用MFC的主动控制优于被动控制。本文方法在处理具有复杂柔性体的航天器时更具优势,更适合于工程应用。     

SiC/CVD SiC复合涂层的抗氧化及抗热震性能研究

为了提高C/SiC复合材料耐高温性能,采用泥浆浸渍裂解与真空化学气相沉积(CVD)在材料表面制备了SiC/CVD SiC复合涂层,通过XRD、SEM分析了涂层组成与结构;研究了复合涂层的高温抗氧化(700～1 500℃)和抗热震性能。结果表明,泥浆浸渍法制备的SiC涂层具有一定的封孔效果,可使材料开孔率下降,但高温抗氧化效果并不佳,1 200℃氧化10 min后材料弯曲强度保留率下降明显仅有86%。CVD SiC涂层结构致密,与SiC封孔涂层结合较好,在700～1 500℃具有较好的抗氧化效果,随着氧化温度的升高,氧化后涂层完好,表面O元素逐渐增加,材料失重率缓慢增加但不大于0.5%,且材料性能并未下降。涂层材料在1 200℃-10 min短时热震5次后材料弯曲强度保留率仍有95%以上,且未出现开裂、剥落等热震损伤。在1 200℃-30 min长时热震10次后,涂层材料基本被完全氧化,材料失去保护作用,弯曲强度下降至90%左右。     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Cutting path-dependent machinability of SiCp/Al composite under multi-step ultra-precision diamond cutting

Particle-tool interactions, which govern the synergetic deformation of SiC particle reinforced Al matrix composites under mechanical machining, strongly depend on the geometry of particle position residing on cutting path. In the present work, we investigate the influence of cutting path on the machinability of a SiCp/Al composite in multi-step ultra-precision diamond cutting by combining finite element simulations with experimental observations and characterization. Be consistent with experimentally characterized microstructures, the simulated SiCp/Al composite is considered to be composed of randomly distributed polygonally-shaped SiC particles with a volume fraction of 25vol%. A multi-step cutting strategy with depths of cut ranging from 2 to 10 μm is adopted to achieve an ultimate depth of cut of 10 μm. Intrinsic material parameters and extrinsic cutting conditions utilized in finite element simulations of SiCp/Al cutting are consistent with those used in corresponding experiments. Simulation results reveal different particle-tool interactions and failure modes of SiC particles, as well as their correlations with machining force evolution, residual stress distribution and machined surface topography. A detailed comparison between numerical simulation results and experimental data of multi-step diamond cutting of SiCp/Al composite reveals a substantial impact of the number of cutting steps on particle-tool interactions and machined surface quality. These findings provide guidelines for achieving high surface finish of SiCp/Al composites by ultra-precision diamond cutting.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

Plasma transport process in the equatorial/low-latitude ionosphere

The behaviour of the equatorial/low-latitude ionosphere and the transport processes during magnetic disturbed and quiet periods of a high solar activity year, 2014, in the American sector are investigated. Parameters used include vertical drift (Vz), transport term (W), NmF2, hmF2 and scale-height (H). The F2 plasma variations followed the diurnal local solar pattern, being higher at daytime. The sunset maximum and sunrise minimum peaks of hmF2 were directly opposite to the scale height (H) pattern. The plasma distribution was basically controlled by combined actions of the electrodynamic convection/thermospheric composition, which is geomagnetic activity dependent. The annual, semi-annual and winter-anomalies of the F2 parameters were higher at the dip equator in comparison with the low-latitude. The Vz pre-reversal peak magnitude coincided with hmF2 peak and the effects are more pronounced during geomagnetic disturbed conditions. The transport term pattern was similar to that of the scale height and it is suggested as a proxy parameter for quantifying low-latitude plasma irregularities and distribution of thermospheric composition.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.