Neumann 随机有限元法计算涡轮叶片的可靠度

推导了具有随机性时Ｎｅｕｍａｎｎ随机有限元法的基本方程,用级数来拟合结构响应的概率密度函数,并通过数值积分计算结构的可靠度。对涡轮叶片的可靠度分析表明,该方法可利用较小的样本容量获得较高的计算精度,并具有较高的计算效率。     

利用遗传算法和神经网络响应面来实现复合材料结构优化设计

运用正交试验设计选择设计样本，建立神经网络响应面，以代替复合材料结构优化中的大量的有限元分析；将神经网络响应面作为目标函数或者约束条件，加上其他常规约束条件进行优化模型的建立，再应用遗传算法（GA）进行优化，这可以实现设计分析与设计优化的分离。以复合材料帽型加筋板的重量优化问题为例，建立了重量响应面目标函数、强度和翘曲稳定性响应面约束条件；并通过NASTRAN进行有限元计算，以获取用于响应面训练的样本点数据。研究表明，该方法能以较少的结构分析次数，取得高精度的响应面近似模型，从而使优化效率大为提高。神经网络响应面能够获得与传统响应面同等，甚至更好的精度。     

计算涡轮叶片结构可靠度的自适应重要抽样法

为直接准确地计算系统结构可靠度，提出了自适应重要抽样方法，由极限状态面失效区确定的最初抽样区开始，用与随机变量联合概率密度函数成比例的抽样密度，自适应递增调整抽样，使抽样区稍大于失效区，以减少安全区的追加抽样，提高抽样效率，并将其应用于航空发动机涡轮叶片的结构可靠性分析，在获得相同的准确计算结果时，所需计算量显著减少。     

伪随机谱载荷下的可靠性模型

齐次马尔可夫模型可以用于描述伪随机谱载荷下的裂纹扩展的随机过程。文中导出了用初始裂纹长度的分布及马尔可夫转移概率密度函数计算给定时刻裂纹长度分布的积分公式。建立了由实验数据估计马尔可夫转移概率密度函数的计算方法。在此基础上建立了计及裂纹检出概率的可检结构的可靠性模型。     

基于优化Kriging模型和重要抽样法的结构可靠度混合算法

 结构可靠度分析计算通常采用多项式响应面拟合隐式极限状态函数,但对于复杂航空航天机械结构产品极限状态方程往往表现出高度非线性,此时多项式响应面的模拟精度不够就会造成计算不收敛。为了提高结构可靠度计算的精度、效率和收敛性,提出了基于优化Kriging模型和重要抽样法的结构可靠度计算方法。首先,利用人工蜂群算法对Kriging模型的参数进行优化;再用优化后的模型模拟隐式极限状态函数,结合重要抽样法不断修正抽样重心,逐步提高模拟精度以达到给定要求;最后,结合一阶矩法(FORM)/二阶矩法(SORM)经典算法求解结构可靠度。该方法提高了高度非线性隐式极限状态方程可靠度计算的精度和收敛性,并且具有较高的计算效率。     

薄壁结构高温随机振动疲劳分析方法有效性验证

为了研究航空薄壁结构高温随机振动疲劳破坏机理,得到可靠的疲劳寿命分析方法,针对根部固支的GH188薄壁结构进行了数值仿真。重点研究了不同温度和不同振动量级组合下,薄壁结构危险点位置轴向动应力响应规律。采用改进的雨流计数法绘制出雨流循环矩阵和雨流损伤矩阵,结合疲劳累积损伤理论估算薄壁结构的疲劳寿命。通过高温随机振动疲劳试验对上述仿真结果进行验证。结果表明:数值仿真对结构破坏位置判断准确,不同温度和随机振动载荷作用下响应的计算值与试验值获得很好的一致性,响应峰值频率误差在1%~3%,结构疲劳寿命与试验结果处于同一量级,证明了高温随机振动疲劳分析方法的有效性和精度以及仿真结果的可靠性。     

结构可靠性分析中考虑高次项修正的组合响应面法

响应面法被广泛地应用于大型复杂结构隐式极限状态方程的可靠性分析,然而至今响应面法的精度仍是未完全得以解决的问题。分析并指出了影响响应面法精度的主要因素,针对这些影响响应面法精度的主要因素,提出了考虑高次项修正的组合响应面法,该方法在真实极限状态方程对结构失效概率贡献较大的区域以一组主次响应面来近似,并对每一个响应面采用高次项予以修正,从而使得响应法的精度得以大幅提高,算例结果证明了所提方法的可行性与精度。     

热声激励下高温合金薄壁结构振动响应试验验证与疲劳寿命预测

为了研究金属薄壁结构在热声循环激励下形成的内部残余应力(Residual Stress)对模态频率、响应规律及疲劳性能的影响，采用耦合的有限元/边界元法(FEM/BEM)对薄壁结构进行振动响应计算，重点研究带残余应力结构的动力学响应特性，并与试验结果进行对比分析。结果表明，不同温度下基频的计算值与试验值获得很好的一致性，计算的应变响应结果与试验测试结果量值相当，验证了热声响应计算方法与模型的有效性。最终结合改进雨流计数法(Improved Rain-flow Counting)、Morrow平均应力模型、Miner线性损伤累积理论计算随机疲劳寿命。结果表明，带残余应力结构在屈曲前后温度为50℃，150℃，250℃所对应的损伤量级分别为10-7，10-4，10-7，疲劳寿命也表现出先减小后增大趋势，且硬化现象使疲劳寿命整体降低，寿命谷值随温度变化发生偏移。     

对偶型最大熵概率密度函数模型及其优化解法

针对经典型最大熵概率密度函数模型及其计算目前存在的非线性程度高,优化不收敛,求解效率低等问题,提出了一种对偶型最大熵概率密度函数模型+逐次优化的方法.根据优化过程不稳定,重新推导了拉格朗日系数的线性变换公式.针对几种常见及一种复杂的概率密度函数,采用经典型与对偶型最大熵概率密度函数模型分别计算概率密度及可靠度的对比表明:与经典型最大熵概率密度函数模型相比,对偶型最大熵概率密度函数模型优化函数形式简单,非线性程度低.逐次优化法求解拉格朗日系数不仅克服了初始值敏感性问题,而且计算效率高.对偶型最大熵概率密度函数模型+逐次优化法与其他方法相比,计算精度最高,且能很好的应用于复杂概率分布及可靠性问题.      

热声载荷作用下薄壁结构非线性响应分析和试验验证

针对热声载荷作用下薄壁结构大挠度非线性响应问题，开展了固支约束金属薄壁板结构热声激励试验及数值模拟分析。通过计算结果与试验结果对比，表明两者结果存在一致性，进而验证了薄壁板在热声载荷作用下动态响应计算方法和数值模型的有效性。在此基础上，针对加筋板结构完成了多种热声载荷组合作用下的动力学响应计算，获得了时域位移响应。重点对该结构在后屈曲状态下的3种典型振动形式进行分析，总结出热载荷与声载荷之间的相对强度决定了板的不同跳变形式。采用统计分析方法建立了位移响应的概率谱密度函数（PDF）并绘图，清楚地显示了后屈曲板的PDF表现出双峰现象。使用功率谱密度（PSD）函数分析了响应频率和峰值随着温度升高的变化，并确定了板的软化和硬化区域。总结了屈曲前/后结构特定区域拉应力和压应力的变化规律，并阐述了造成这种变化的原因。本文工作可对热声载荷作用下薄壁结构响应分析和动强度设计提供参考依据。     

An Optimum First-Order Time Delay Tracker

The key to the general estimation problem is the computation ofthe posterior probability density function (PDF) of the desiredparameter, or PDF conditioned on knowing all availableobservations. Any optimum estimate can be computed from this PDF. The theory and implementation of a numerical solution ispresented for the computation of the posterior PDF in real time fortime delay tracking with a first-order Markov assumption. Benefitsinclude a lack of threshold effects usually associated with nonlineartrackers as well as theoretically optimum performance.     

测试性验证中抽样方案的精确算法及应用

通过对二项分布函数单调性的讨论,提出了求解二项分布函数联立不等式的精确算法,可迅速准确地获得测试性验证的抽样方案。并研制了测试性验证软件包,可以较好地满足测试性验证的要求。     

Reliability analysis based on a novel density estimation method for structures with correlations

Estimating the Probability Density Function (PDF) of the performance function is a direct way for structural reliability analysis,and the failure probability can be easily obtained by integration in the failure domain.However,efficiently estimating the PDF is still an urgent problem to be solved.The existing fractional moment based maximum entropy has provided a very advanced method for the PDF estimation,whereas the main shortcoming is that it limits the application of the reliability analysis method only to structures with independent inputs.While in fact,structures with correlated inputs always exist in engineering,thus this paper improves the maximum entropy method,and applies the Unscented Transformation (UT) technique to compute the frac tional moments of the performance function for structures with correlations,which is a very efficient moment estimation method for models with any inputs.The proposed method can precisely estimate the probability distributions of performance functions for structures with correlations.Besides,the number of function evaluations of the proposed method in reliability analysis,which is determined by UT,is really small.Several examples are employed to illustrate the accuracy and advantages of the proposed method.     

反应机理对湍流预混值班火焰稳燃影响的研究

针对反应机理对湍流预混值班火焰稳燃的影响,使用大涡模拟(LES)耦合输运概率密度(PDF)方法模拟了湍流射流火焰的燃烧流场。化学反应过程采用被广泛应用的SMOOKE (16组分,46步)和DRM22 (22组分,104步)机理,概率密度输运的确定采用欧拉随机场来实现。计算结果表明,LES/PDF模型可以准确地模拟该类火焰的速度分布,火焰褶皱等结构。反应机理的点火延迟时间的准确与否是模拟值班射流火焰的关键,预混射流出口处的小尺度掺混引起的点火和火焰传播过程是值班射流火焰稳燃的关键。     

Optimum time delay detection and tracking

An optimum solution to the time delay tracking problems has been previously achieved by recursively computing the time-delay posterior probability density function (PDF), or time-delay PDF conditioned on knowing all observed data, by numerical means. Any optimum estimate is determined from the posterior PDF. The same approach is extended to an optimum solution for the time delay detection problem with a randomly appearing and disappearing signal. The detector-tracker is suitable for a stand-alone detection and tracking system     

Non-Gaussian random vector identification using sphericallyinvariant random processes

With the modeling of non-Gaussian radar clutter in mind, elegant and tractable techniques are presented for characterizing the probability density function (PDF) of a correlated non-Gaussian radar vector. The need for a library of multivariable correlated non-Gaussian PDFs in order to characterize various clutter scenarios is discussed. Specifically,. the theory of spherically invariant random processes (SIRPs) is examined in detail. Approaches based on the marginal envelope PDF and the marginal characteristic function have been used to obtain several multivariate non-Gaussian PDFs. An important result providing the PDF of the quadratic form of a spherically invariant random vector (SIRV) is presented. This result enables the problem of distributed identification of a SIRV to be addressed     

可靠性全局灵敏度求解的Meta-IS-AK算法

可靠性全局灵敏度（GRS）可以衡量输入变量对结构系统失效概率的平均影响,但目前仍然缺乏具有广泛适应性的高效算法。针对此问题,本文将在元重要抽样和可靠性全局灵敏度的贝叶斯算法基础上建立一种新的高效算法。所提算法首先利用已有的贝叶斯算法,将可靠性全局灵敏度转换成由无条件失效概率及输入变量失效域条件下的概率密度函数（PDF）表达的形式,然后分3步来完成算法的组织。第1步是利用元重要抽样的迭代策略抽取失效域的重要抽样样本;第2步是在已有的元重要抽样法中嵌入自适应Kriging模型,高效计算出无条件失效概率;第3步是利用Metropolis-Hastings准则,将失效域的重要抽样样本转化成为原始密度函数在失效域的样本点,进而同时求得各个输入变量在失效域中的条件概率密度函数,并最终求得可靠性全局灵敏度。由于所提算法充分利用了已有的可靠性全局灵敏度贝叶斯算法的维度独立性、元重要抽样法对隐式多失效域的适应性以及元重要抽样法中嵌入式Kriging模型的高效性,因此所提算法具有广泛的适用范围和较高的效率,该结论得到了算例结果的充分验证。     

Combustion modeling in a model combustor

The flow-field of a propane-air diffusion flame combustor with interior and exterior conjugate heat transfers was numerically studied.Results obtained from four combustion models,combined with the re-normalization group(RNG) k-ε turbulence model,discrete ordinates radiation model and enhanced wall treatment are presented and discussed.The results are compared with a comprehensive database obtained from a series of experimental measurements.The flow patterns and the recirculation zone length in the combustion chamber are accurately predicted,and the mean axial velocities are in fairly good agreement with the experimental data,particularly at downstream sections for all four combustion models.The mean temperature profiles are captured fairly well by the eddy dissipation(EDS),probability density function(PDF),and laminar flamelet combustion models.However,the EDS-finite-rate combustion model fails to provide an acceptable temperature field.In general,the flamelet model illustrates little superiority over the PDF model,and to some extent the PDF model shows better performance than the EDS model.      

从探测概率的角度评价飞机的隐身性能

采用雷达散射截面(RCS)均值来衡量飞机的隐身性能并不能给出足够充分的信息,从信号检测概率的角度来衡量飞机的可探测性可以提供更完整的信息。本文从探测概率的角度详细分析了4种典型隐身飞机RCS起伏数据与虚警概率、探测概率、探测距离、信噪比(SNR)等参数间的关系,并给出了完整的推导过程,采用此方法可以对任意飞机目标在给定雷达参数下的可探测性进行准确计算。一般在雷达性能评估中需要用到RCS起伏模型,本文对飞机目标回波信号进行检测分析时采用数值计算方法,信号的概率密度函数(PDF)直接来源于原始回波起伏信号,避免了模型拟合带来的回波特性失真,且不会产生大的计算误差。通过对比分析计算结果发现:单脉冲检测下回波信号均值中值比越大,回波取值范围越小,则信号的检测概率越小,其中均值中值比一般相差3倍以上时可以得到明显不同的检测概率;在快起伏假设下,非相参积累检测的积累脉冲数Nin即便较小,也能得到较大的信噪比增益,此增益可能大于Nin。