基于序列近似优化算法的固体运载火箭弹道设计

采用打靶法研究固体运载火箭弹道优化设计问题。针对打靶法中智能优化算法效率、优化问题中约束条件提法等方面的不足,引入序列近似优化算法,提出算法中径向基函数高斯核宽度的高效确定方法。数值仿真实验表明,提出的方法可在基本不带来计算量增加的前提下,显著提高代理模型精度,使序列近似优化算法得到改进,提高优化效率。设计固体运载火箭飞行程序。建立弹道优化设计问题数学模型,并将模型中相关等式约束合理转化为不等式约束,降低优化问题求解难度;基于改进后的序列近似优化算法完成某固体运载火箭弹道优化,原始计算模型调用308次之后便搜索到最优解,较传统智能优化算法显著提高了优化效率,优化方案末助推级液体推进剂消耗比原方案减少25.7%。     

无筒空射运载火箭重力出舱机箭耦合动力学

研究了以运输机为平台的内装式空射运载火箭重力出舱载机-火箭两体动力学。根据火箭受力条件和相对运动自由度,将其出舱运动顺次分为5个阶段。给出了各阶段过渡的力学条件,基于牛顿-欧拉法分别建立了前4个阶段的机-箭两体动力学模型。定义了火箭出舱过程可能发生危险的多种异常情况,给出了发生异常的力学或几何条件。然后对火箭出舱全程作了数值仿真。若发射初始条件合理,仿真将顺次经历前4个或3个正常阶段,直至火箭离舱。如果异常情况的条件满足,火箭出舱将发生危险,仿真失去意义,故中止,需要调整发射初始条件后进一步仿真验证。本文给出的系统动力学模型和仿真数据可为工程部门设计空射型火箭的本体参数,以及设置空射初始条件提供参考。     

基于间接法在线能力评估和自主规划技术研究

针对运载火箭推力故障条件下的自主救援任务,研究了一种基于间接法的在线能力评估和自主规划方法.构建了以入轨点能量最大为最优性能指标的最优控制问题.通过数值方法探讨了在满足轨道倾角和入轨点路径角约束条件下,入轨点速度、远地点和近地点跟入轨高度之间的关系,设计了一种与入轨高度相关的特征量来对入轨能力进行评估.以近地点最高为救...     

Mathematical modeling of plasma drifts over equatorial low latitude regions

This paper presents a mathematical model to simulate ionospheric plasma drifts at equatorial low latitude regions by coupling of E- and F-regions. The governing non-linear differential equations (of elliptic and parabolic nature) are solved numerically through finite-difference schemes and obtained neutral winds and electric fields. The temperature and electron density profiles are generated utilizing MSIS-86 atmospheric model. The continuity equation is employed to obtain night-time E-region density profile using measured ionograms at Trivandrum (India). The computed vertical and zonal plasma drifts are comparable with measured Jacamarca plasma drifts with little variations during noon and evening times. The plasma drifts at Trivandrum (8.5° N, 76.5° E, dip 0.5° N) are compared with those of Jicamarca (12° S, 76.9° W, dip 2° N). Neutral wind simulations of present model agree well with those of horizontal wind model (HWM-93). The post-sunset enhancement and its reversal are also discussed.     

Vertical motions in Thailand after the 2004 Sumatra–Andaman Earthquake from GPS observations and its geophysical modelling

Following previous findings from ongoing GPS research in Thailand since 2004 we continue to exploit the GPS technique to monitor and model land motions induced by the Sumatra–Andaman Earthquake. Our latest results show that up to the end of 2010, Thailand has been co-seismically displaced and is subsequently undergoing a post-seismic horizontal deformation with total displacements (co-seismic plus post-seismic) ranging from 10.5 to 74.7 cm. We observed the largest horizontal displacements in the southern part of Thailand and moderate and small displacements in the central and northern parts. In addition to horizontal displacements throughout Thailand, continuous GPS measurements show that large parts of Thailand are subsiding at rates up to 1 cm/yr. It is the first time that such vertical post-seismic deformations at large distances (650–1500 km away from the Earthquake’s epicentre) have been recorded. We have investigated the physical processes leading to the observed subsidence. While after-slip on the subduction interface induces negligible or even slightly positive vertical motions, relaxation in the asthenosphere is associated with a sizable subsidence. Predictions from a 3D finite element model feature an asthenosphere with an effective viscosity of the order of 3 * 10¹⁸ Pas, fit the horizontal post-seismic data and the observed subsidence well. This model is then used to predict the subsidence over the whole seismic cycle. The subsidence should go on with a diminishing rate through the next two decades and its final magnitude should not exceed 10 cm in the Bangkok area.     

牵制缓释放过程中火箭动力响应特性分析

为研究运载火箭在牵制缓释放过程中的结构动力响应特性,将牵制缓释放发射过程依次分为静态竖立、点火牵制和缓释放三个阶段进行计算。运用MSC.Patran有限元程序的场功能来实现上一阶段的全部计算结果场向下一阶段初始条件场的传递,运用MSC.Nastran有限元程序的非线性弹簧单元的非线性位移载荷函数功能来模拟运载火箭缓释放过程中的缓释力,实现运载火箭牵制缓释放过程结构动力响应的数值计算,并对比分析运载火箭几种常用发射释放方式的结构动力响应特性。结果表明:采用牵制缓释放系统可有效减小运载火箭释放时所受冲击载荷,减小运载火箭全箭结构动力响应。     

Influence of the gravitational vertical on geometric visual illusions

The occurrence of geometric orientation illusions and the perception of ambiguous figures were analyzed in 24 subjects during static body tilt relative to gravity on Earth. Results showed that illusions such as the Rock's diamond/square, the Ponzo illusion, and orientation contrast illusions occurred less frequently, and that depth reversal of ambiguous figures took more time when subjects were lying on their side or supine compared to upright, thus suggesting that the gravitational reference plays a significant role in these “visual” illusions. The structure of images, our representation of the environment, and orientation relative to gravity are all integral parts in interpreting visual images. In a weightless environment where no gravitational reference can be used, it is expected that similar alterations in visual perception will occur.     

多等级区域侦察弹性星座设计方法

针对传统侦察星座目标单一、弹性低的问题，提出了多等级区域侦察弹性星座的设计方法。该方法将星座设计过程按区域等级信息分为多个子星座逐步设计，直到整体星座对所有的区域性能满足设计要求。以区域被划分为3个等级为例，首先对星座设计需求、设计指标及设计步骤进行了分析。其次推导了地面最低分辨率和轨道高度的关系并确定了不同子星座的轨道高度。最后考虑轨道倾角、一箭多星发射、光照和升交点漂移同步约束，构建基础星座、子星座1和子星座2的优化模型。最终设计星座为3层混合星座，共8个轨道面和70颗卫星，星座对各等级区域的最大重访时间分别为10937s，12241s和17437s，弹性指数为2213%，2420%和6361%。结果表明该方法设计的星座可实现对区域覆盖和弹性分级的设计要求，证明了方法的有效性。对比Walker星座设计方法，在同等设计要求下，Walker星座所需卫星数为156颗，多等级区域侦察弹性星座所需卫星数远低于Walker星座，结果进一步证明了该星座设计方法的优越性。     

Numerical study of traveling ionosphere disturbances with vertical incidence data

High frequency ionosphere vertical sounding, as an important and representative application for detecting the ionosphere and studying the characteristics of radio propagation, can be utilized to monitor the ionosphere continuously variation and to acquire the ionosphere asymmetrical features of diverse scale above the ionosphere vertical sounding stations. This is a first article on real time application of numerical methods to obtain the parameters of traveling ionosphere disturbances (TIDs) using vertical incident ionograms. In this paper, the distribution of ionosphere electron density with TIDs is constructed using a background ionosphere model superimposed a perturbation theory model. The background ionosphere electron density is modelled by the inversion of vertical incident ionograms which are observed before the appearance of the disturbance. Based on the fourth order Adams-Bashforth-Moulton (the so-called ABM) predictor corrector method, instead of Runge-Kutta method, the fast digital ray tracing method is established. According to process of the disturbed trace simulation and parameters inversion, the characteristic parameters of ionosphere disturbance at different detection time can be obtained in real time. The numerical analysis of TIDs is then captured completely.     

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.