结构首次穿越失效研究的新进展

在结构首次穿越失效分析中,初始条件、成群穿越、结构的非线性及抗力退化等问题影响着分析方法的复杂性和分析结果的精确性。结合近年来的研究进展,本文讨论和分析了解决上述问题的新思路、新技术和新方法,同时,归纳了几个新的研究热点,以便为结构首次穿越失效的进一步理论和应用研究提供有益的参考。     

利用奇异谱分析提高极移中长期预报精度

由于空间大地观测数据传输耗时及处理过程复杂, 导致极移测量值的获取存在时延, 无法满足对高精度的极移预报值有重大需求的应用领域. 针对极移复杂的时变特性, 提出一种基于奇异谱分析(singular spectrum analysis, SSA)的预报方法. 首先用SSA分离提取极移时序中的高频组分与低频组分; 其次建立最小二乘(least square, LS)外推与自回归(AutoreGressive, AR)模型对极移高频和低频组分进行组合预报. 结果表明, SSA方法能够准确地分离和提取极移低频和高频组分, 对低频和高频组分组合预报可以显著改善极移的中长期(30--365d)预报精度, 与国际地球自转和参考系服务局(International Earth Rotation and Reference Systems Service, IERS)提供的A公报中的极移预报值相比, SSA方法对极移PMX分量(本初子午线方向)和PMY分量(西90$^\circ$子午线方向)的中长期预报精度改进最高分别可达45.97%和62.44%. 研究结果验证了SSA方法对极移中长期预报改进的有效性.     

大型射电望远镜主反射面调整方法研究

随着毫米波天文学和空间通信的重要性日益提高, 对天线性能提出了越来越高的要求, 而天线性能往往受到其反射器表面精度的限制. 微波全息技术是一种快速有效的检测反射面天线表面轮廓的测量技术. 通过微波全息测量得到天线口径场, 计算天马65m射电望远镜反射面与理想抛物面的偏差. 天马65m射电望远镜的主反射面板是放射状的, 有14圈. 面板的每个角都固定在面板下方促动器的螺栓上进行上下移动, 且相邻面板交点处的拐角共用一个促动器. 采用平面拟合的方法可以计算各块面板拐角处的调整值, 但是同一个促动器会得到4个不同的调整量. 通过平面拟合, 同时以天线照明函数为权重的平差计算方法得到相邻面板拐角的一个平差值, 即天马65m射电望远镜1104个促动器的最佳调整值. 通过多次调整和新算法的应用, 天马65m射电望远镜反射面的面形精度逐渐提高到了0.24mm.     

高层大气密度的午夜极大值特征分析与建模

利用CHAMP (CHAllenging Minisatellite Payload)、GRACE-A (Gravity Recovery and Climate Experiment-A)、SWARM-C (The Earth''s Magnetic Field and Environment Explorers-C)等3颗极轨卫星的资料, 研究360—480km高层大气密度在低纬度区域的午夜极大值(Midnight Density Maximum, MDM)现象. MDM一般出现在23:00- 02:00 LT (Local Time)之间,峰值位置在低纬度15°以内,谷值位置在中纬度35°-45°附近,整体略偏向南半球,振幅约为平均密度的26%.随着高度增大以及太阳辐射水平的增强,MDM振幅呈减小趋势;冬至和夏至日附近的季节效应会减弱MDM振幅,在春秋分日的振幅最大.用3个主流大气模型DTM2000 (Drag Temperature Model 2000), NRLMSISE00 (US Naval Research Laboratory, Mass Spectrom-eter and Incoherent Scatter radar)和JB2008 (Jacchia- Bowman 2008 model)对MDM进行模拟,JB2008没有刻画出MDM现象;另两个模型低估了MDM效应,在360km和480km两个高度DTM2000模型的振幅仅为观测的46%和53%, NRLMSISE00模型仅为观测的33%和26%;模型没有准确刻画出MDM与高度、辐射水平和季节的关系.联合3颗卫星的资料,研究了-种基于地理纬度的6阶勒让德多项式,同时融合地方时和高度因素的经验函数,在振幅和相位上可以较好地刻画MDM特征,相关系数达到0.923,可为大气密度模型的修正提供借鉴,服务于低轨道航天器高精度轨道预报.     

Research on Parallel Algorithms for uv-faceting Imaging

The uv-faceting imaging is one of the widely used large field of view imaging technologies, and will be adopted for the data processing of the low-frequency array in the first stage of the Square Kilometre Array (SKA1). Due to the scale of the raw data of SKA1 is unprecedentedly large, the efficiency of data processing directly using the original uv-faceting imaging will be very low. Therefore, a uv-faceting imaging algorithm based on the MPI (Message Passing Interface)+OpenMP (Open Multi-Processing) and a uv-faceting imaging algorithm based on the MPI+CUDA (Compute Unified Device Architecture) are proposed. The most time-consuming data reading and gridding in the algorithm are optimized in parallel. The verification results show that the results of the proposed two algorithms are basically consistent with that obtained by the current mainstream data processing software CASA (Common Astronomy Software Applications), which indicates that the proposed two algorithms are basically correct. Further analysis of the accuracy and total running time shows that the MPI+CUDA method is better than the MPI+OpenMP method in both the correctness rate and running speed. The performance test results show that the proposed algorithms are effective and have certain extensibility.     

The circulation of the Indian Ocean at 32°S

Using inverse methods a circulation for a new section along 32°S in the Indian Ocean is derived with a maximum in the overturning stream function (or deep overturning) of 10.3 Sv at 3310 m. Shipboard and Lowered Acoustic Doppler Current Profiler (ADCP) data are used to inform the choice of reference level velocity for the initial geostrophic field. Our preferred solution includes a silicate constraint (−312 ± 380 kmol s⁻¹) consistent with an Indonesian throughflow of 12 Sv. The overturning changes from 12.3 Sv at 3270 m when the silicate constraint is omitted to 10.3 Sv when it is included. The deep overturning varies by only ±0.7 Sv as the silicate constraint varies from +68 to −692 kmol s⁻¹, and by ±0.3 Sv as the net flux across the section, driven by the Indonesian throughflow, varies from −7 to −17 Sv with an appropriately scaled silicate flux constraint. Thus, the overturning is insensitive to the size of the Indonesian throughflow and silicate constraint within their apriori uncertainties. We find that the use of the ADCP data adds significant detail to the horizontal circulation. These resolved circulations include the Agulhas Undercurrent, deep cyclonic gyres and deep fronts, features evidenced by long term integrators of the flow such as current meter and float measurements as well as water properties.     

Photometric Follow-up Transit (Primary Eclipse) Observations of WASP-43 b and TrES-3b and A Study on Their Transit Timing Variations

Two photometric follow-up transit (primary eclipse) observations on WASP-43 b and four observations on TrES-3 b are performed using the Xuyi Near-Earth Object Survey Telescope. After differential photometry and light curve analysis, the physical parameters of the two systems are obtained and are in good match with the literature. Combining with transit data from a lot of literature, the residuals (O ? C) of transit observations of both systems are fitted with the linear and quadratic functions. With the linear fitting, the periods and transit timing variations (TTVs) of the planets are obtained, and no obvious periodic TTV signal is found in both systems after an analysis. The maximum mass of a perturbing planet located at the 1:2 mean motion resonance (MMR) for WASP-43 b and TrES-3 b is estimated to be 1.826 and 1.504 Earth mass, respectively. By quadratic fitting, it is confirmed that WASP-43 b may have a long-term TTV which means an orbital decay. The decay rate is shown to be P? = (?0.005248 ± 0.001714) s·yr^?1, and compared with the previous results. Based on this, the lower limit of the stellar tidal quality parameter of WASP-43 is calculated to be ${Q^{\prime }}_{*}\ge 1.5\times {10}^5$, and the remaining lifetimes of the planets are presented for the different ${Q^{\prime }}_{*}$ values of the two systems, correspondingly.     

Low-energy Lunar Trajectories with Lunar Flybystwo

The low-energy lunar trajectories with lunar flybys are investigated based on the Sun-Earth-Moon bicircular problem (BCP). The characteristics of the distribution of trajectories in the phase space are summarized. Using the invariant manifolds in the BCP system, the low-energy lunar trajectories with lunar flybys are sought. Then, take time as an augmented dimension in the phase space of a nonautonomous system, we present the state space map and reveal the distribution of these lunar trajectories in the phase space. Consequently, we find that the low-energy lunar trajectories exist as families, and that the every moment in the Sun-Earth-Moon synodic period can be the departure date. Finally, we analyse the velocity increment, transfer duration, and system energy for the different trajectory families, and obtain the velocity-impulse optimal family and the transfer-duration optimal family, respectively.     

Quadrature-based approach for the efficient evaluation of surge hazard

The Joint Probability Method (JPM) has been used for hurricane surge frequency analysis for over three decades, and remains the method of choice owing to the limitations of more direct historical methods. However, use of the JPM approach in conjunction with the modern generation of complex high-resolution numerical models (used to describe winds, waves, and surge) has become highly inefficient, owing to the large number of costly storm simulations that are typically required. This paper describes a new approach to the selection of the storm simulation set that permits reduction of the JPM computational effort by about an order of magnitude (compared to a more conventional approach) while maintaining good accuracy. The method uses an integration scheme called Bayesian or Gaussian-process quadrature (together with conventional integration methods) to evaluate the multi-dimensional joint probability integral over the space of storm parameters (pressure, radius, speed, heading, and any others found to be important) as a weighted summation over a relatively small set of optimally selected nodes (synthetic storms). Examples of an application of the method are shown, drawn from the recent post-Katrina study of coastal Mississippi.     

高压超高压变质岩抬升—剥露的一种动力学机制分析

综述了高压超高压变质岩带的基本构造特征，认为运用体粘滞性流变学理论可以很好地解释高压超高压变质岩的发育过程及其抬升－剥露动力学机制。