缺口件疲劳寿命预测新方法

采用有限元方法针对缺口件多轴疲劳实验结果进行模拟。结果表明,相同路径条件下,随着缺口半径减小缺口根部附近的应力梯度显著增加。基于缺口根部应变值,采用等效应变法进行疲劳寿命预测,预测结果随缺口半径的减小而偏于保守。采用应力梯度法确定有效距离,相同路径下,随着缺口半径的减小有效距离减小;依据该有效距离处的等效应变进行疲劳寿命预测,总体预测结果较为分散且偏于不安全。基于实验及有限元模拟结果,提出了基于应变梯度的有效距离确定的新方法,大部分疲劳寿命预测结果位于2倍分散带内。     

电磁测力系统研究

介绍了电磁测力系统中静磁场分析与优化、通电线圈位置控制方法及提高线性度的自校准方法。运用Maxwell 3D进行仿真计算，得到了相对优化的磁密及磁均匀性，通过控制通电线圈在磁场中位置及分段线性微分技术，提高了电磁测力准确度。构建了试验装置，对比试验数据表明，该研究措施对于提高电磁测力准确度取得了较好效果。     

基于Kriging模型梯度解析解的改进一次二阶矩方法

改进一次二阶矩（AFOSM）法是一种基于功能函数梯度的结构可靠性分析方法，鉴于其对隐式函数的梯度较难求解，提出了一种基于Kriging模型梯度解析解的AFOSM方法，利用Kriging代理模型的解析表达式推导得到功能函数对输入变量的梯度解析解，为AFOSM中设计点的确定提供高精度的梯度信息。通过Kriging与AFOSM的结合，很好地解决了基于有限元模型的隐式情况下梯度计算量相当大、可靠性分析难的问题。数值与工程算例验证了所提Kriging梯度解析解的较高精确性，同时验证了所提基于Kriging解析解的AFOSM结构可靠性分析方法的正确性与较高精度。     

高主流湍流度下密度比对涡轮导叶全气膜冷却特性的影响

为了研究高主流湍流度下二次流密度比对涡轮导叶全气膜冷却特性的影响,使用热色液晶测量了在主流湍流度为15%,二次流密度比为1.0和1.5下三维涡轮导叶的气膜冷却效率和换热系数。二次流与主流质量流量比为7.0%和12.5%。结果表明:二次流密度比增大可以降低冷气射流的动量,小流量比工况下,在叶片前缘和压力面前半段,动量较低的二次流在高主流湍流度的影响下更易耗散,增大二次流密度比使冷却效率明显降低;大流量比工况下,二次流动量降低使气膜孔后区域冷气贴附性增强,气膜冷却效率和冷气覆盖效果均得到提升。小流量比工况下,二次流密度比增大对叶片表面换热的影响较小;大流量比工况下,二次流密度比增大使吸力面中弦区域和压力面后半段的平均换热系数比分别降低15%和25%。     

中国民航旅客量季节性问题研究

中国民航旅客量的季节性波动是经济学研究的重要命题。选取中国民航旅客量时间序列数据作为分析对象，采用年度季节强度指数（R）与月季指数（Zt），揭示中国民航旅客量的年度季节集中性与月度季节变化规律。 R值整体上呈现逐年递增的趋势，指标随季节变化明显，且按一定规律分布。 Zt 值指示4月、7月、8月、9月、10月为旺季，其余月份为淡季。民航旅客量淡旺季与游客出游行为、节假日长短、气候等因素有直接联系。研究结论为解决中国民航旅客量的淡旺季问题提供基础数据，也为调整航空公司在淡旺季时的人员、航班及其他工作安排，为缓解航空客流量供需不平衡的矛盾提供决策参考。     

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.     

含多部位损伤搭接结构应力强度因子三维有限元分析

应力强度因子(Stress intensity factor,SIF)分析是含多部位损伤(Multiple site damage,MSD)结构剩余强度和裂纹扩展寿命预测的基础和关键。考虑接触与摩擦,建立了含MSD搭接结构的三维有限元模型,研究了不同裂纹长度、铆钉类型以及损伤模式下裂纹尖端SIF分布情况和变化规律。结果表明,搭接件孔边裂纹Ⅰ型SIF起主导作用,Ⅱ型和Ⅲ型SIF可忽略不计。由于次弯曲、铆钉变形和板厚度等因素,SIF在外表面最小,接触面一侧较大,最大值多位于蒙皮内部。MSD会使裂纹间的干涉作用增强,SIF增大,且裂纹间距离越近干涉作用越强。裂纹长度相同时,埋头铆钉的孔边裂纹SIF积分均值大于平头铆钉,且接触面的SIF埋头铆钉大于平头铆钉,外表面则相反。     

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.     

An Artificial Neural Network based approach for estimation of rain intensity from spectral moments of a Doppler Weather Radar

By using a Doppler Weather Radar (DWR) at Shriharikota (13.66°N & 80.23°E), an Artificial Neural Network (ANN) based technique is proposed to improve the accuracy of rain intensity estimation. Three spectral moments of a Doppler spectra are utilized as an input data to an ANN. Rain intensity, as measured by the tipping bucket rain gauges around the DWR station, are considered as a target values for the given inputs. Rain intensity as estimated by the developed ANN model is validated by the rain gauges measurements. With the help of a developed technique, reasonable improvement in the estimation of rain intensity is observed. By using the developed technique, root mean square error and bias are reduced in the range of 34–18% and 17–3% respectively, compared to Z–R approach.     

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°.