基于不同天文标准计算地球引力对卫星轨道的影响

建立了一个考虑地球非球形引力摄动的卫星轨道摄动模型并用C＋＋语言编写程序进行轨道计算,其积分器采用Runge-Kutta-Fehlberg7（8）。建模过程中,对岁差、章动等量的计算分别考虑了美国海军天文台1981年和2005年的两种天文标准。通过与STK对比,发现无论采用哪一种标准,轨道传播两天后的位置误差都不会超过分米量级,速度误差不会超过毫米每秒的量级。说明了建模和编程计算的合理性,以及采用新旧两种天文标准计算卫星轨道的差别甚小。     

中尺度电离层扰动的统计色散特性

利用Millstone Hill(288.5E,42.6N)地区的数字测高仪漂移模观测数据,统计分析了中尺度电离层行扰(TID)的色散特性.分析结果表明,该地区的TID在扰动频率-相速度平面的分布中存在若干明显的扰动出现率峰值.比较相应的理论研究结果发现,观测到的TID大都可以用不同的受导重力波传播模式加以解释: 频率较高的TID可以用地面和热层底部维持的各阶完全制导模解释,频率相对较低的TID可以用传统理论给出的由热层底部陡峭的温度梯度维持的L₁模和我们新近提出的由热层背景风场所维持的各阶导管模解释.      

低重环境下圆柱贮箱内液体晃动特性研究的一种解析法

考虑低重环境下圆柱贮箱内由于表面张力影响而呈现弯曲自由液面的情况,利用Fourier Bessel级数对贮箱受横向激励时的自由液面处的边界条件进行展开,得到液体晃动系统的广义状态方程,给出了固有频率、晃动波高、晃动力和晃动力矩等晃动特性的计算公式。通过数值算例与文献结果对比,在验证文中方法正确性的同时,具体研究了各晃动特性随充液深度、外激励频率和Bond数等参数变化的规律。     

Harmonic analysis of the ionospheric electron densities retrieved from FORMOSAT-3/COSMIC radio occultation measurements

In order to investigate the regular variations of the ionosphere, the least-squares harmonic estimation is applied to the time series of ionospheric electron densities derived from about five years of Global Positioning System radio occultation observations by FORMOSAT-3/COSMIC satellites. The analysis is done for different latitudes and altitudes in the region of Iran. The least-squares harmonic estimation is found to be a powerful tool for the frequency analysis of the completely unevenly spaced time series of radio occultation measurements. Although the obtained results are slightly different from the exact expected cycles (i.e. annual and diurnal components with their Fourier decompositions, and the 27-day period) due to the low horizontal resolution of radio occultation measurements, high vertical resolution of the observations enables us to detect not only the total electron content variations but also periodic patterns of electron densities at different altitudes of the ionosphere. The dominant diurnal and annual signals together with their Fourier series decompositions are obtained, which are consistent with the previous analyses on the total electron content. In the equatorial anomaly band, the annual component is weaker than its Fourier decomposition periods. In particular, the semiannual period dominates the annual component, indicating the relationship between the semiannual variation of the electron densities and the ionospheric equatorial anomaly. From detection of the phases of the components, it is revealed that the annual signal generally has its maximum value in summers at high altitudes, and in the winters at low altitudes. This is probably due to the higher [O/N₂] ratios in winter than in the summer in the lower ionosphere. Furthermore, the semiannual component mostly peaks around solstices or about a month before/after them.     

非相干雷达探测低电离层的结果分析

位于波多黎各的Arecibo非相干雷达可以获得低电离层电子和离子密度, 利用此非相干雷达数据对中纬度低电离层的运动特征进行研究. 得到了电子密度随时间和高度的变化 情况, 结果显著呈现出周日变化特征, 并分析了电子密度随高度的变化规律. 进一步对数据进行频谱分析, 深入研究低电离层电子密度的周日变化效应. 得到电子密度的高度剖面, 发现从F层底部到E层有明显的等离子体沉降. 低电离层的层结构特征及电子密度变化表明, 在该区域还存在不同程度的等离子体扰动, 由此对低电离层的作用因素 进行分析, 认为大气潮汐或声重波可能对低电离层产生扰动, 即低电离层与大气存在一定程度的耦合作用.      

An experimental method to obtain the hard alpha anomaly distribution for titanium alloy aeroengine disks

A probability-based damage tolerance methodology has been proposed to improve the recognition of material anomalies, especially hard alpha (TiN) anomalies for aeroengine rotor disks. A key input to this method is hard alpha anomaly distribution, which reflects the occurrence rate and size of anomalies present in the finished part material of titanium rotors. Since anomalies rarely occur naturally, an experimental method is proposed to obtain the anomaly distribution for titanium alloy aeroengine disks to reflect and equivalently replace the manufacturing development in titanium industry. In general, the anomaly distribution information can be divided into two parts: the Probability of Detection (POD) curve and the detected anomaly distribution, which contains the size and frequency data of the detected anomalies. The distribution can be established based on several appropriate assumptions and derivation steps with different confidence levels of POD curves and detected anomaly distributions. In this case, the distribution can be obtained in a relatively short time as a key input to the probability-based damage tolerance methodology. Then the Probability of Failure (POF) can be calculated, and the value is found to vary with different confidence levels. On this basis, the conservative estimated POF can be obtained in conjunction with confidence levels.     

考虑重心位置不确定性的无尾无人机优化设计

飞机的飞行性能与重心(CG)位置密切相关,尤其是后掠式无尾飞机的重心位置对其飞行性能影响更甚,如果重心位置发生变化,升力分布随之改变,进而影响飞机航时。针对这个问题,从气动布局和设计方法两方面,设计了一种航时对重心位置不敏感的无尾无人机(UAV)。气动布局上,提出了利用螺旋桨动力配平纵向力矩的鸥翼(GW)布局,以减小重心位置变化对升阻特性的影响;设计方法上,采用稳健性优化设计(RDO)理论,分析重心位置不确定时的航时低敏感度问题。以一架小型电动无人机为研究对象,建立了无尾无人机稳健性优化设计环境,包括总体设计、代理模型构造以及稳健性优化。分析结果表明:利用螺旋桨动力配平的鸥翼布局使重心可用范围增加了5%;静安定裕度在5%~15%变化时,该布局可以有效提高航时稳健性。采用稳健性优化得到的无人机几何参数,大幅度降低了重心位置对航时的影响,显著提升了满足约束的概率。     

Thermal and TEC anomalies detection using an intelligent hybrid system around the time of the Saravan,Iran, (Mw = 7.7) earthquake of 16 April 2013

A powerful earthquake of M_w = 7.7 struck the Saravan region (28.107° N, 62.053° E) in Iran on 16 April 2013. Up to now nomination of an automated anomaly detection method in a non linear time series of earthquake precursor has been an attractive and challenging task. Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) have revealed strong potentials in accurate time series prediction. This paper presents the first study of an integration of ANN and PSO method in the research of earthquake precursors to detect the unusual variations of the thermal and total electron content (TEC) seismo-ionospheric anomalies induced by the strong earthquake of Saravan. In this study, to overcome the stagnation in local minimum during the ANN training, PSO as an optimization method is used instead of traditional algorithms for training the ANN method. The proposed hybrid method detected a considerable number of anomalies 4 and 8 days preceding the earthquake. Since, in this case study, ionospheric TEC anomalies induced by seismic activity is confused with background fluctuations due to solar activity, a multi-resolution time series processing technique based on wavelet transform has been applied on TEC signal variations. In view of the fact that the accordance in the final results deduced from some robust methods is a convincing indication for the efficiency of the method, therefore the detected thermal and TEC anomalies using the ANN + PSO method were compared to the results with regard to the observed anomalies by implementing the mean, median, Wavelet, Kalman filter, Auto-Regressive Integrated Moving Average (ARIMA), Support Vector Machine (SVM) and Genetic Algorithm (GA) methods. The results indicate that the ANN + PSO method is quite promising and deserves serious attention as a new tool for thermal and TEC seismo anomalies detection.     

Some ionospheric storm effects at equatorial and low latitudes

In this paper, the response of the equatorial and low latitude ionosphere to three intense geomagnetic storms occurred in 2002 and 2003 is reported. For that, critical frequency of F2-layer foF2 and the peak height hmF2 hmF2 for the stations Jicamarca (11.9°S), Ascension Is (7.92°S) and Tucuman (26.9°S) are used. The results show a “smoothing” of the Equatorial Anomaly structure during the development of the storms. Noticeable features are the increases in foF2 before the storm sudden commencement (SC) at equatorial latitudes and the southern crest of the Equatorial Anomaly. In some cases nearly simultaneous increases in foF2 are observed in response to the storm, which are attributed to the prompt electric field. Also, positive effects observed at equatorial and low latitudes during the development of the storm seem to be caused by the disturbance dynamo electric field due to the storm-time circulation. Increases in foF2 above the equator and simultaneous decreases in foF2 at the south crest near to the end of a long-duration main phase are attributed to equatorward-directed meridional winds. Decreases in foF2 observed during the recovery phase of storms are believed to be caused by composition changes. The results indicate that the prompt penetration electric field on the EA is important but their effect is of short lived. More significant ionospheric effects are the produced by the disturbance dynamo electric field. The role of storm-time winds is important because they modify the “fountain effect” and transport the composition changes toward low latitudes.     

Surface latent heat flux anomalies quasi-synchronous with ionospheric disturbances before the 2007 Pu’er earthquake in China

One of various mechanisms of pre-earthquake lithosphere–atmosphere–ionosphere coupling as possible explanation of the seismo-ionospheric effects is connected with the release of latent heat. Abnormal variations of ionospheric electromagnetic parameters possibly related to the 2007 Ms 6.4 Pu’er earthquake in China were reported. This paper attempts to examine whether there were abnormal changes of surface latent heat flux (SLHF) linked with this pre-earthquake ionospheric disturbances. The spatio-temporal statistical analyzes of multi-years SLHF data from USA NCEP/NCAR Reanalysis Project reveal that local SLHF enhancements appeared 11, 10 and 7 days before the Pu’er earthquake, respectively. As contrasted to the formerly reported local ionospheric Ne enhancement 9 days before the shock observed by DEMETER satellite, it is discovered that the SLHF and Ne anomalies are quasi-synchronous and have good spatial correspondence with the epicentre and the local active faults. This is valuable for understanding the seismogenic coupling processes and for recognizing earthquake anomaly with multiple parameters from integrated Earth observation system.