近地轨道卫星的地影预报算法

针对传统地影预报算法计算量大,不适合星载计算机进行地影时刻自主预报的问题,提出一种地影预报星上算法。通过构造降维坐标系将卫星进出地影过程转换为“星-地-日”平面内的几何问题;基于常数变易法推导出阴影条件的隐式解析公式,其方程系数可由实时轨道参数确定。为简化计算,将星上算法设计为双层结构：构造迭代算法,依据实时轨道参数可精确预报卫星在每个节点上的地影时刻;而在任意两个节点之间采用解析算法进行近似预报。仿真结果表明,星上算法的预报精度高且计算量小,适合星载计算机进行地影时刻自主预报。     

垂直返回重复使用运载火箭技术分析

结合SpaceX公司近期进行的多次猎鹰火箭一子级垂直着陆返回技术试验,对比分析垂直返回重复使用运载火箭两种返回方式的工程应用价值。首先,建立运载火箭一子级动力返回段弹道设计动力学模型。随后,提出基于H-V飞行剖面分段返回弹道设计方法。然后针对“返回原场”（RTLS）和“不返回原场”(NRTLS)两种垂直返回方式,构建综合考虑上升段与返回段的推进剂耦合作用的一体化弹道优化设计策略。最后,通过数值仿真,对比分析了两种返回方式下的火箭运载能力。结果表明,采用“不返回原场”方式的运载能力损失仅占“返回原场”方式的一半,具有较好的工程应用价值。     

采用弱攻角补偿与脱敏设计的火星进入段制导

为提升火星进入段存在多种扰动时的制导末端精度,在现有三自由度脱敏设计的基础上,提出针对大气密度及升阻力系数波动的弱攻角补偿方法。通过攻角调整,协调升、阻力加速度测量值相对理论计算值的偏离程度,使加权形式的偏离程度指标趋近于1,从而降低气动参数波动对制导精度的影响。相平面分析和反证法证明了攻角调整过程的稳定性。蒙特卡洛仿真结果表明,该方法可达到较高的纵向末端状态精度。     

Direct trajectory optimization based on a mapped Chebyshev pseudospectral method

 In view of generating optimal trajectories of Bolza problems, standard Chebyshev pseudospectral (PS) method makes the points' accumulation near the extremities and rarefaction of nodes close to the center of interval, which causes an ill-condition of differentiation matrix and an oscillation of the optimal solution. For improvement upon the difficulties, a mapped Chebyshev pseudospectral method is proposed. A conformal map is applied to Chebyshev points to move the points closer to equidistant nodes. Condition number and spectral radius of differentiation matrices from both methods are presented to show the improvement. Furthermore, the modification keeps the Chebyshev pseudospectral method's advantage, the spectral convergence rate. Based on three numerical examples, a comparison of the execution time, convergence and accuracy is presented among the standard Chebyshev pseudospectral method, other collocation methods and the proposed one. In one example, the error of results from mapped Chebyshev pseudospectral method is reduced to 5% of that from standard Chebyshev pseudospectral method.     

Trajectory control strategy of cathodes in blisk electrochemical machining

A turbine blisk, which combines blades and a disk together, is one of the most important components of an aero engine. In the process of blisk electrochemical machining (ECM), the sheet cathode, which is usually used as a tool electrode, has a complicated structure. In addition to that, the channel between the adjacent blades is narrow and twisted, so interference is apt to happen when the sheet cathode feeds into the channel. Therefore, it is important to choose suitable trajectory control strategy. In this paper, a new trajectory control strategy of the sheet cathode is presented and corresponding simulation analysis is conducted on the basis of an actual blisk model. The simulation results demonstrate that the sheet cathode can feed into the channel by a spatial line trajectory without interference. Moreover, the verification experiments are carried out according to the simulation. The experimental results show that the cathode can move into the channel without interference. It is verified that the new trajectory control strategy is correct and can be used in the blisk ECM process successfully.     

飞机驾驶舱视觉舒适性研究

驾驶舱视觉舒适性会直接对驾驶员产生影响,可能引起驾驶员视觉疲劳和操作失误,对飞机驾驶的安全性和驾驶员的身体健康非常不利.采用日光眩光概率（DGP）、亮度对比度（LR）和人眼垂直照度（Ev）指标,对不同的季节、日照时间和天气条件下的飞机驾驶舱进行视觉舒适性的仿真研究,形成数据库.结果表明：该方法可在飞机设计阶段对飞机驾驶舱的视觉舒适性进行有效地评估,研究结果可为飞机驾驶舱防眩布置及照明设计提供一定的指导.     

Estimates of vertical land motion along the southwestern coasts of Turkey from coastal altimetry and tide gauge data

The differences between coastal altimetry and sea level time series of tide gauges in between March 1993 and December 2009 are used to estimate the rates of vertical land motion at three tide gauge locations along the southwestern coasts of Turkey. The CTOH/LEGOS along-track coastal altimetry retrieves altimetric sea level anomalies closer to the coast than the standard along-track altimetry products. However, the use of altimetry very close to the coast is not found to improve the results. On the contrary, the gridded and interpolated AVISO merged product exhibits the best agreement with tide gauge data as it provides the smoothest variability both in space and time compared with along track altimetry data. The Antalya gauge to the south (in the Mediterranean Sea) and the Mentes/Izmir gauge to the west (in the Aegean Sea) both show subsidence while the Bodrum tide gauge to the south (in the Aegean Sea) shows no significant vertical land motion. The results are compared and assessed with three independent geophysical vertical land motion estimates like from GPS. The GIA effect in the region is negligible. The VLM estimates from altimetry and tide gauge data are in good agreement both with GPS derived vertical velocity estimates and those inferred from geological and archaeological investigations.     

Equatorial F-region vertical ion drifts during quiet solar maximum

Median values of ionosonde h′F data acquired at Ibadan (Geographic:7.4°N, 3.9°E, Magnetic: dip 6°S, and magnetic declination, 3°W), Nigeria, West Africa, have been used to determine vertical ion drift (electric field) characteristics in the postsunset ionosphere in the African region during a time of high solar activity (average F_10.7 −208). The database spans from January and December 1958 during the era of International Geophysical Year (IGY) for geomagnetic quiet conditions. Bimonthly averaged diurnal variations patterns are very similar, but differ significantly in magnitude and in the evening reversal times. Also, monthly variations of F-region vertical ion drift reversal times inferred from the time of h′F maximum indicates early reversal during equinoxes and December solstice months except for the month of April. Late reversal is observed during the June solstice months. The equatorial evening prereversal enhancement in vertical ion drift (V_zp) occurs largely near 1900 LT with typical values 20–45 m/s. Comparison of Ibadan ionosonde V_zp with the values of prereversal peak velocity reported for Jicamarca (South America), Kodaikanal (India), and Scherliess and Fejer global model show considerable disparity. The changes of postsunset peak in virtual height of F-layer (h′F_P) with prereversal velocity peak V_zp are anti-correlated. Investigation of solar effects on monthly values of V_zp and h′F_P revealed that these parameters are independent of monthly averaged solar flux intensity during quiet-time sunspot maximum conditions.     

Application of support vector machine combined with K-nearest neighbors in solar flare and solar proton events forecasting

The support vector machine (SVM) combined with K-nearest neighbors (KNN), called the SVM-KNN method, is new classing algorithm that take the advantages of the SVM and KNN. This method is applied to the forecasting models for solar flares and proton events. For the solar flare forecasting model, the sunspot area, the sunspot magnetic class, and the McIntosh class of sunspot group and 10 cm solar radio flux are chosen as inputs; for the solar proton event forecasting model, the inputs include the longitude of active regions, the flux of soft X-ray, and those for the solar flare forecasting model. Detailed tests are implemented for both of the proposed forecasting models, in which the SVM-KNN and the SVM methods are compared. The testing results demonstrate that the SVM-KNN method provide a higher forecasting accuracy in contrast to the SVM. It also gives an increased rate of ‘Low’ prediction at the same time. The ‘Low’ prediction means occurrence of solar flares or proton events with predictions of non-occurrence. This method show promise for forecasting models of solar flare and proton events.     

Study of the Enhancement Events of Relativistic Electron at the Geosynchronous Orbit

The high flux of energetic electron on geostationary orbit can induce many kinds of malfunction of the satellite there, within which the bulk-charging is the most significant that several broadcast satellite failures were confirmed to be due to this effect. The electron flux on geostationary orbit varies in a large range even up to three orders accompanied the passage of interplanetary magnetic cloud and the following geomagnetic disturbances. Upon the investigation of electron flux enhancement events, two types of events were partitioned as recurrent events and random ones. Both of the two kinds of events relate to the interplanetary conditions such as solar wind parameters, IMF etc and their evolution characters as well. As for the recurrent events, we found that, (1) all of the events exhibits periodic recurrence about 27 days, (2) significant increase of electron flux relates to interplanetary index and characters of their distribution, (3) the electron flux also has relation to solar activity index. An artificial neural network was constructed to estimate the flux I day ahead. The random electron flux enhancement events are rare and present different distribution figures to the recurrent ones. The figure of the random events and the conditions of their occurrence is also discussed in this paper.