基于自适应Kriging代理模型的交叉熵重要抽样法

对于复杂失效域和小失效概率耦合的可靠性分析问题,本文提出了一种交叉熵重要抽样（CE-IS）方法结合自适应Kriging （AK）代理模型的求解方法（CE-IS-AK）。所提方法基于交叉熵原理,用混合高斯模型逐步逼近最优重要抽样密度函数,并采用AK模型协助逼近过程中混合高斯模型的参数的更新,从而提高了CE-IS方法的计算效率。另外,本文还改进了CE-IS方法的收敛准则,避免了方法的冗余迭代,扩大了方法的适用范围。由于在CE-IS方法中引入了AK模型,因此,本文方法所构建的重要抽样函数在保证精度的基础上提高了效率。相较于AK-MCS方法,本文方法中引入了重要抽样的思想,因此在Kriging训练点数目基本相同的情况下,大幅缩减小失效概率计算时样本池规模,并且由于利用了混合高斯模型,因而对多失效域具有较好的适用性。算例分析也证明了本文所提方法的优越性。     

基于随机粒子追踪方法的环形燃烧室点火动态模拟

为了深入研究低排放燃烧室点火联焰规律,在全新的环形模型燃烧室中开展了点火模拟和试验研究。点火模拟采用随机粒子追踪方法,能够基于时均冷态流场的仿真结果快速模拟火焰传播过程。环形燃烧室包含16个中心分级旋流器,仅向预燃级通入丙烷,用于模拟航空发动机低排放燃烧室点火状态下的空气燃油分级。试验采用PIV技术测量3个头部区域流场,利用高速相机拍摄火焰CH*/C2*基团化学发光信号。对多个流量和当量比条件下的联焰过程、联焰时间和传焰速率进行了分析,试验和模拟的结果均表明：环形燃烧室内火焰双向传播,燃烧室内外环流速度差异导致了双向火焰传播速度差,传焰速率随燃烧室湍流速度和当量比的增加而增加。点火模型很好地捕捉了环形燃烧室点火动态,所得传焰速率也符合湍流火焰传播规律,表明该模型具有较强预测能力。     

多重应力下滚动轴承剩余寿命预测

为了实现多重应力下滚动轴承的剩余寿命预测,有效利用不用应力下的退化数据,提出了一种基于加速模型和贝叶斯（Bayesian）理论的滚动轴承剩余寿命预测方法。通过拟合优度检验和威布尔（Weibull）概率图检验法对滚动轴承试验中的数据进行有效性分析。利用switching Kalman filters(SKF)判断滚动轴承各时刻的退化状态。当滚动轴承进入加速退化时,用指数模型拟合轴承退化过程,利用广义线性对数模型表示退化模型参数与应力的关系,根据修正后的轴承实时退化数据利用贝叶斯算法更新模型参数,得到滚动轴承剩余寿命的概率密度函数,从而实现滚动轴承剩余寿命预测。采用XJTU-SY轴承数据集进行验证,预测结果的均方根误差在20 min以内,证明该方法能够有效预测滚动轴承的剩余寿命。     

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Global Navigation Satellite Systems (GNSS), in particular the Global Positioning System (GPS), have been widely used for high accuracy geodetic positioning. The Least Squares functional models related to the GNSS observables have been more extensively studied than the corresponding stochastic models, given that the development of the latter is significantly more complex. As a result, a simplified stochastic model is often used in GNSS positioning, which assumes that all the GNSS observables are statistically independent and of the same quality, i.e. a similar variance is assigned indiscriminately to all of the measurements. However, the definition of the stochastic model may be approached from a more detailed perspective, considering specific effects affecting each observable individually, as for example the effects of ionospheric scintillation. These effects relate to phase and amplitude fluctuations in the satellites signals that occur due to diffraction on electron density irregularities in the ionosphere and are particularly relevant at equatorial and high latitude regions, especially during periods of high solar activity. As a consequence, degraded measurement quality and poorer positioning accuracy may result.     

Comparison of ionospheric F2 peak parameters foF2 and hmF2 with IRI2001 at Hainan

Monthly median values of foF2, hmF2 and M(3000)F2 parameters, with quarter-hourly time interval resolution for the diurnal variation, obtained with DPS4 digisonde at Hainan (19.5°N, 109.1°E; Geomagnetic coordinates: 178.95°E, 8.1°N) are used to investigate the low-latitude ionospheric variations and comparisons with the International Reference Ionosphere (IRI) model predictions. The data used for the present study covers the period from February 2002 to April 2007, which is characterized by a wide range of solar activity, ranging from high solar activity (2002) to low solar activity (2007). The results show that (1) Generally, IRI predictions follow well the diurnal and seasonal variation patterns of the experimental values of foF2, especially in the summer of 2002. However, there are systematic deviation between experimental values and IRI predictions with either CCIR or URSI coefficients. Generally IRI model greatly underestimate the values of foF2 from about noon to sunrise of next day, especially in the afternoon, and slightly overestimate them from sunrise to about noon. It seems that there are bigger deviations between IRI Model predictions and the experimental observations for the moderate solar activity. (2) Generally the IRI-predicted hmF2 values using CCIR M(3000)F2 option shows a poor agreement with the experimental results, but there is a relatively good agreement in summer at low solar activity. The deviation between the IRI-predicted hmF2 using CCIR M(3000)F2 and observed hmF2 is bigger from noon to sunset and around sunrise especially at high solar activity. The occurrence time of hmF2 peak (about 1200 LT) of the IRI model predictions is earlier than that of observations (around 1500 LT). The agreement between the IRI hmF2 obtained with the measured M(3000)F2 and the observed hmF2 is very good except that IRI overestimates slightly hmF2 in the daytime in summer at high solar activity and underestimates it in the nighttime with lower values near sunrise at low solar activity.     

Upper limit of the total magnetic flux in an active region according to the photometric-magnetic dynamical model

The photometric-magnetic dynamical model handles the evolution of an individual sunspot as an autonomous nonlinear, though integrable, dynamical system. One of its consequences is the prediction of an upper limit of the sunspot areas. This upper limit is analytically expressed by the model parameters, while its calculated value is verified by the observational data. In addition, an upper limit for the magnetic strength inside the sunspot is also predicted, and then, we obtain the following significant result: The upper limit of the total magnetic flux in an active region is found to be of about 7.23 × 10²³ Mx, namely, phenomenologically equal to the magnetic flux concentrated in the totality of the granules of the quiet Sun, having a typical maximum magnetic strength of about 12G. Therefore, the magnetic flux concentrated in an active region cannot exceed the magnetic flux concentrated in the photosphere as a whole.     

Using the radial basis function neural network to predict ionospheric critical frequency of F2 layer over Wuhan

Neural networks (NNs) have been applied to ionospheric predictions recently. This paper uses radial basis function neural network (RBF-NN) to forecast hourly values of the ionospheric F₂ layer critical frequency(f_oF₂), over Wuhan (30.5N, 114.3E), China. The false nearest neighbor method is used to determine the embedding dimension, and the principal component analysis (PCA) is used to reduce noise and dimension. The whole study is based on a sample of about 26,000 observations of f_oF₂ with 1-h time resolution, derived during the period from January 1981 to December 1983. The performance of RBF-NN is estimated by calculating the normalized root-mean-squared (NRMSE) error, and its results show that short-term predictions of f_oF₂ are improved.     

液体火箭发动机涡轮泵用机械密封温度场及热载变形研究

基于ANSYS数值计算软件,建立了液体火箭发动机涡轮泵用机械密封的二维稳态传热模型,依靠经验公式确定了模型的对流换热系数。计算了密封环的温度场和热载变形,分析了密封端面比压、回流流量以及不同材质对密封温度场的影响规律。结果表明:密封端面最高温度发生在靠近密封环内径处,且密封端面比压越大密封环温度梯度越大;密封环热载变形呈收敛间隙,最大变形发生在动环端面的外径处,其值约为2.2μm;密封环端面最高温度随回流流量增加而减小,当回流流量从0.1~0.6 kg/s变化时,密封环端面最高温度可降低18%(从100℃降至82℃);当回流流量增大到0.3 kg/s时,继续提高对密封环端面温升的控制不再显著;采用高导热系数的摩擦副材料能够显著降低端面温升和温度梯度,提高密封工作可靠性。     

考虑侧窗约束的模型预测静态规划末制导方法

针对红外导引头侧窗探测模式下,非对称视场约束造成末制导阶段目标易丢失的问题,提出一种考虑侧窗约束的模型预测静态规划末制导方法。首先,基于体视线坐标系建立三维相对运动模型,得到不依赖于“小攻角”假设的准确模型。在处理无过程约束问题的模型预测静态规划方法基础上,引入松弛变量与虚拟控制量,设计出考虑侧窗视场约束的末制导算法。为了进一步降低末制导算法对初始猜测轨迹的依赖性,提高适应性与计算效率,提出逐步增加约束条件的计算策略。仿真结果表明,该方法在末制导过程中满足侧窗约束,相比于凸优化方法,优化变量减少,计算速度更快;相比于基于障碍李雅普诺夫函数的末制导律,能够满足侧窗约束,同时能适应不同的初始条件。     

复合材料桨叶结构损伤演化模型

与金属材料桨叶相比,复合材料桨叶因具有更加优良的抗疲劳性能而被广泛应用到直升机旋翼上。但由于复合材料破坏机理复杂,疲劳性能分散,影响因素众多,导致复合材料桨叶疲劳现象尚处于研究探索之中,在复合材料的微观失效机制与宏观结构的力学性能之间仍然缺少一座桥。鉴于此,文章利用典型复合材料试样的拉伸疲劳实验数据,建立了基体裂纹、纤维断裂和界面脱胶等损伤变量累积模型,从断裂能的角度出发构建了基体裂纹密度、纤维断裂面积与复合材料属性之间的函数关系,分析了基体裂纹密度、纤维断裂面积等损伤变量对复合材料工程性能参数的影响。利用复合材料宏观力学理论,研究了各物理损伤变量对桨叶刚度特性的影响,采用连续损伤变量的状态方程建立了复合材料桨叶的损伤演化模型,这种以有理多项式为状态转移函数微分模型能很好地体现复合材料桨叶在疲劳初期和疲劳末期刚度快速损伤的现象。