Space-Wise,Time-Wise,Torus and Rosborough Representations in Gravity Field Modelling

The decade of the geopotentials started July 2000 with the launch of the German high-low SST mission CHAMP. Together with the joint NASA-DLR low-low SST mission GRACE and the ESA gradiometry mission GOCE an unprecedented wealth of geopotential data becomes available over the next few years. Due to the sheer number of unknown gravity field parameters (up to 100 000) and of observations (millions), especially the latter two missions are highly demanding in terms of computational requirements. In this paper several modelling strategies are presented that are based on a semi-analytical approach. In this approach the set of normal equations becomes block-diagonal with maximum block-sizes smaller than the spherical harmonic degree of resolution. The block-diagonality leads to a rapid and powerful gravity field analysis tool. Beyond the more-or-less conventional space-wise and time-wise formulations, the torus approach and Rosborough's representation are discussed. A trade-off between pros and cons of each of the modelling strategies will be given. This revised version was published online in August 2006 with corrections to the Cover Date.     

离散事件动态系统的GSMP模型及扰动分析

对离散事件动态系统（ＤＥＤＳ）的仿真过程建立了包含中断事件集和驱动事件集的广义半马尔可夫过程（ＧＳＭＰ）模型，对参数的无穷小扰动所引起的系统状态转移时刻进行了扰动分析，最后，讨论了扰动分析所需要的可变换条件和样本性能函数的关系，并提出值得进一步研究的问题。     

飞机最快上升转弯轨迹的计算

70年代初,Miele A等人提出SGRA方法用于优化飞行轨迹计算。但方法存在时间加权系数难于确定、长时间的飞行轨迹计算积分精度不易控制及要计算数字微分等不足。对此提出一改进梯度方法,并计算了某型战斗机从亚音速三维转弯升高到超音速的最短时间飞行轨迹。结果表明,本方法稳定有效,为解决飞机最快机动轨迹的优,乜提供了一种行之有效的手段。     

制造误差对气体静压轴承涡流力矩影响分析方法研究

 采用有限元方法研究了制造误差对狭缝节流气体静压轴颈—止推轴承的涡流力矩的影响。对于轴颈—止推相连结构的气体轴承,通过相容变换进行统一编程计算;在离散化过程中,利用加权余量法将二阶偏微分方程降低一阶,放松了对插值函数连续度的要求,便于借助有限元技术分析狭缝节流气体静压轴承的流场参数。分析了狭缝气膜宽度误差和轴颈圆度误差对涡流力矩的影响,以及轴颈的不同安装角度、偏心等因素对涡流力矩的影响。经对比验证,有限元计算结果与实测结果基本一致,研究结果对于气体静压轴颈—止推轴承的设计、装配优化和性能预测有重要指导意义。     

钛合金板材激光弯曲成形的研究

 对常温下钛合金板材进行了激光弯曲成形研究。分析了激光功率、扫描速度、光斑大小对板材弯曲变形的影响,获得理想工艺参数范围为:功率1 .0～1 .2 k W、扫描速度2～3 m/min、光斑直径6～7mm。同时分析了板材厚度、宽度及扫描次数对材料弯曲变形的影响,为钛合金板材弯曲成形提供了一种新途径。     

承力操纵装置性能试验台系统方案设计

介绍了一个飞机承力操纵装置性能试验台的实现机理、采用的方案,有助于增加读者对液压系统、液压性能试验台及液压伺服控制系统的了解。     

DME Interference mitigation for L-DACS1 based on system identification and sparse representation

L-band digital aeronautical communication system 1 (L-DACS1) is a promising candi-date data-link for future air-ground communication, but it is severely interfered by the pulse pairs (PPs) generated by distance measure equipment. A novel PP mitigation approach is proposed in this paper. Firstly, a deformed PP detection (DPPD) method that combines a filter bank, correlation detection, and rescanning is proposed to detect the deformed PPs (DPPs) which are caused by mul-tiple filters in the receiver. Secondly, a finite impulse response (FIR) model is used to approximate the overall characteristic of filters, and then the waveform of DPP can be acquired by the original waveform of PP and the FIR model. Finally, sparse representation is used to estimate the position and amplitude of each DPP, and then reconstruct each DPP. The reconstructed DPPs will be sub-tracted from the contaminated signal to mitigate interference. Numerical experiments show that the bit error rate performance of our approach is about 5 dB better than that of recent works and is closer to interference-free environment.     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

低速大转矩永磁游标电机在发酵罐中的应用

采用秸秆生产丁醇时,为了达到高效率发酵,需要使发酵液中菌体分布均匀,所加物料需要及时搅拌。工业化生产丁醇的发酵罐尺寸较大,搅拌过程对驱动电机的输出转矩要求较高,因此高效率、小体积、低速大转矩的秸秆发酵搅拌电机对工业化发酵生产丁醇意义重大。设计了一种定转子永磁游标电机用于发酵液搅拌,并对电机的槽极数、电磁性能进行了具体的研究。制造样机并采用滞环控制实现了对液态发酵对象的快速响应,且实现了低能耗、高效率的发酵过程。     

Semi-stochastic modification of second-order radial derivative of Abel–Poisson’s formula for validating satellite gravity gradiometry data

The geoid can be used to validate the satellite gravity gradiometry data. Validation of such data is important prior to their downward continuation because of amplification of the data errors through this process. In this paper, the second-order radial derivative of Abel–Poisson’s formula is modified stochastically to reduce the effect of the far-zone geoid and generate the second-order radial derivative of geopotential at 250 km level. The numerical studies over Fennoscandia show that this method yields the gradients with an error of 10 mE and when the long wavelength of geoid is removed from the estimator and restored after the computations (remove–compute–restore) the error will be in 1 mE level. We name this method semi-stochastic modification. The best case scenario is found when the degree of modification of the integral formula is 200 and the long wavelength geoid to degree 100 is removed and restored. In this case the geoid should have a resolution of 15′ × 15′ and the integration should be performed over a cap size of 3°.