Orbit-attitude-vibration coupled dynamics of tethered solar power satellite

A new orbit-attitude-vibration coupled dynamic model of the tethered solar power satellite (Tethered SPS) is established based on absolute nodal coordinate formulation. The Hamilton’s equation of the system is derived by introducing generalized momentum through Legendre transformation. The correctness of the proposed model is verified by an example. The dynamic characteristics of the Tethered SPS are studied using the symplectic Runge-Kutta method. Simulation results show that the orbital radius and the total energy of the system are well preserved. The attitude of the system is unstable when the mass of the bus system is small. However, the attitude stability is dependent on some other parameters of the system, which requires further studies. It is also found that the average tether force/deformation can be roughly estimated by simplifying the solar panel as a particle. The proposed model can be used to study the orbit-attitude-vibration coupled dynamics and control problems.     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

Study of ionospheric D region changes during solar flares using MF radar measurements

During solar flares, the X-ray radiation suddenly increases, resulting in an increase in the electron density of the atmospheric D region and a strong absorption of short-wave radio waves. Based on Langfang medium frequency (MF) radar, this paper analyzed the variation characteristics of D region in the lower ionosphere from 62 km to 82 km. The analysis focused on multiple C-level and M-level solar flare events before and after the large-scale flare event at 11:53 (UT) on September 6, 2017. The results show that it is difficult to detect the electron density over 70 km in Langfang during solar flares, but the electron density value can be obtained as low as 62 km, and the stronger the flare intensity, the lower the detectable electron density height. Besides, the equal electron density height, the received power of X and O waves will also be significantly reduced during the flares, and the reduction of equal electron density height has a weak linear relationship with flare intensity.     

Thermosphere densities derived from Swarm GPS observations

After the detection of many anomalies in the Swarm accelerometer data, an alternative method has been developed to determine thermospheric densities for the three-satellite mission. Using a precise orbit determination approach, non-gravitational and aerodynamic-only accelerations are estimated from the high-quality Swarm GPS data. The GPS-derived non-gravitational accelerations serve as a baseline for the correction of the Swarm-C along-track accelerometer data. The aerodynamic accelerations are converted directly into thermospheric densities for all Swarm satellites, albeit at a much lower temporal resolution than the accelerometers would have been able to deliver. The resulting density and acceleration data sets are part of the European Space Agency Level 2 Swarm products.To improve the Swarm densities, two modifications have recently been added to our original processing scheme. They consist of a more refined handling of radiation pressure accelerations and the use of a high-fidelity satellite geometry and improved aerodynamic model. These modifications lead to a better agreement between estimated Swarm densities and NRLMSISE-00 model densities. The GPS-derived Swarm densities show variations due to solar and geomagnetic activity, as well as seasonal, latitudinal and diurnal variations. For low solar activity, however, the aerodynamic signal experienced by the Swarm satellites is very small, and therefore it is more difficult to accurately resolve latitudinal density variability using GPS data, especially for the higher-flying Swarm-B satellite. Therefore, mean orbit densities are also included in the Swarm density product.     

不同强度磁暴期间热层大气密度响应特征分析

利用CHAMP卫星数据，对2002-2008年12个不同强度磁暴事件期间的热层大气密度变化特征进行分析，并研究对应磁暴期间大气模式NRLMSISE-00分布特征.结果表明，大磁暴期间日侧大气密度峰值从高纬到低纬的时间延迟为2h，中小磁暴期间的延迟时间为3～4h；春秋季暴时大气密度分布基本呈南北对称分布，而夏冬季大气密度的分布是夏半球大于冬半球，春秋季暴时大气密度大于夏冬季；NRLMSISE-00大气模式得到的热层大气密度很好的体现了半球分布以及季节分布的特征，但模式模拟结果偏小；Dst指数峰值比ap指数峰值更能反应大气密度的变化情况.      

无刷双馈电机静止坐标系数学模型的研究

由于没有电刷和滑环,提高了运行的可靠性,降低了维护成本,因此无刷双馈电机(BDFM)在变速恒频发电领域有着广阔的应用前景。由于特殊的结构和工作原理,BDFM数学模型十分复杂。从BDFM多回路空间矢量数学模型出发,推导出了两个静止坐标系之间的坐标变换关系以及功率绕组静止坐标系数学模型和控制绕组静止坐标系数学模型,推导过程简洁明了。所得结果为BDFM理论分析和不同控制策略的研究提供了理论基础。MATLAB/Simulink环境下的仿真结果验证了数学模型及坐标变换关系的正确性。     

一种高动态低信噪比卫星导航信号捕获方法

为提升高动态低信噪比环境下卫星导航信号的捕获性能,提出了一种基于分数阶傅里叶变换(FrFT)及部分匹配滤波(PMF)的捕获方法。在该方法中,接收机首先利用PMF对接收信号做分段相干积分,随后借助快速傅里叶变换(FFT)对分段积分结果做离散快速FrFT,最后通过检测FrFT输出的峰值完成信号的捕获。由于具有多普勒频率变化率的卫星导航信号在FrFT后呈现能量聚焦特性,所提方法能够显著提高信号的长时间相干积分增益。同时对所提算法的捕获概率、平均捕获时间以及算法复杂度等性能指标进行了理论分析及计算机仿真验证。仿真表明,与传统的PMF-FFT方法相比,所提方法能够通过延长相干积分时间的方式有效提升高动态低信噪比卫星导航信号的捕获概率、降低捕获时间。     

考虑局部遮挡的太阳能无人机能源控制

针对太阳能无人机在飞行状态下可能出现的太阳能电池局部遮挡情况，开展相应的太阳能电池最大功率点追踪算法和能源控制研究。通过将发光亮度引入相对吸引力计算过程对萤火虫算法进行改进，实现了局部阴影情况下太阳能电池最大功率点的高效追踪。以此为基础，设计了考虑局部遮挡情况下太阳能无人机的太阳能电池/蓄电池混合能源状态机控制规则。以"蒲公英I"无人机为例，建立了太阳能电池阵列模型，开展了考虑局部遮挡情况下太阳能电池最大功率点追踪仿真实验；基于"蒲公英I"飞行剖面，开展了考虑局部遮挡情况的混合能源控制仿真试验。研究结果表明：改进的萤火虫算法可以实现在局部阴影情况下太阳能电池最大功率点的有效跟踪，与萤火虫算法相比收敛时间更短、且功率波动幅度更小；采用改进萤火虫算法和状态机能源管理策略，在考虑局部遮挡的飞行状态下可以实现太阳能电池/蓄电池之间的合理功率分配与控制。     

太阳能飞机循环飞行的高度剖面能量策略

为满足设定的太阳能飞机多日连续飞行条件，依据飞行过程中当前时刻的飞行高度、光伏输出功率、动力电池组余量等系统状态参数，研究如何分配动力电池组充放电和电推进系统输入等功率。所用策略立足于实时功率平衡，充分利用正午前后的光伏峰值功率用于飞机爬升及充电，在午后下滑过程中利用全部光伏输出，以最大化利用光伏资源；在光伏有效输出不足时则以一定的维持功率下滑，使能量的综合损失最小。方法能够提高以预定夜间飞行高度连续多日续航的成功率，提升飞行高度、纬度、季节范围或搭载能力，或者拓展这几种飞行条件的组合域，优化太阳能飞机的适用性。     

SSPS-OMEGA超大型球形聚光器结构拓扑构型与设计

为同时兼顾千米量级尺度空间太阳能电站聚光系统的在轨发射安装难度、光学收集效率、稳定性及结构工艺等多方面的因素，提出采用基于球面正多面体的经纬线划分法模块化构建球形聚光器。该方法不仅可实现基础划分单元分布均匀、规格种类少，降低在轨构建与运载发射难度。而且只要合理控制弧高等分数和基础划分单元的口径就可以同时保证聚光系统的光学收集效率。理论与仿真结果显示，当弧高等分数为24时，该划分法共有25种基础单元件，整个球面被分成5762块。并且此时模块化构建的实际球形聚光器的光学收集效率与理想球面相比基本保持不变，解决了太阳能电站聚光系统设计中存在的关键问题。