认知不确定性下可靠性灵敏度分析的新指标

传统的可靠性灵敏度指标只考虑了随机不确定性,当结构中还存在认知不确定性时,该指标难以反映出基本变量中认知不确定性对结构可靠性评估结果的影响程度.针对这一问题,引入概率包络来表达结构中的随机不确定性和认知不确定性,在此基础上建立了可靠性灵敏度分析的新指标——认知灵敏度系数,给出了其定义和代表的意义,并基于证据理论提出了认知灵敏度系数的计算方法.算例分析表明,认知灵敏度系数可以反映出基本变量中认知不确定性对结构可靠性评估结果的不确定性的影响程度,从而为有针对性地提高可靠性评估结果的准确性提供依据.另外,基于证据理论的认知灵敏度系数的计算方法适用于极限状态函数为单调、非单调的结构,以及变量服从正态、非正态的结构,具有普适性.     

概率盒全局灵敏度和活跃子空间跨层降维

航空航天仿真系统中的不确定性通常是多源的、混合的,并且系统参数的维数众多。针对高维混合不确定性量化问题,提出一种结合概率盒全局灵敏度和活跃子空间的跨层降维方法。在随机和认知不确定的概率盒表征基础上,使用不确定性缩减法分析参数的全局灵敏度继而进行参数筛选;基于输出梯度协方差矩阵的特征分解,使用活跃子空间法对参数进行降维;构造出一种概率盒表征下的参数筛选和跨层降维方法。最后以NASA多学科不确定性量化挑战问题为例,通过概率盒全局灵敏度分析进行第1层次的参数筛选,原有的21维输入参数减为13维;随后采用活跃子空间进行第2层次的参数降维,维数进一步降至一维。研究结果表明,所提出的方法能够对混合不确定性参数进行灵敏度排序,还能够有效降低模型输入参数的维度,为高维系统混合不确定性量化和进一步的优化工作奠定了基础。     

基于贝叶斯可靠性的设计优化在弹翼蒙皮设计中的应用

针对设计优化中存在的偶然不确定性和认知不确定性,采用贝叶斯理论,利用验前信息,进行可靠性分析,提出了基于贝叶斯可靠性的设计优化方法,并将其应用到弹翼蒙皮设计中。在仿真中,计算了两种不确定性下的可靠度分布,优化结果达到了可靠性要求,并且计算简单,易于工程实现。     

Epistemic uncertainty quantification in flutter analysis using evidence theory

Aimed at evaluating the structural stability and flutter risk of the system, this paper manages to quantify epistemic uncertainty in flutter analysis using evidence theory, including both parametric uncertainty and method selection uncertainty, on the basis of information from limited experimental data of uncertain parameters. Two uncertain variables of the actuator coupling system with unknown probability distributions, that is bending and torsional stiffness, which are both described with multiple intervals and the basic belief assignment(BBA) extricated from the modal test of actuator coupling systems, are taken into account. Considering the difference in dealing with experimental data by different persons and the reliability of various information sources, a new combination rule of evidence––the generalized lower triangular matrices method is formed to acquire the combined BBA. Finally the parametric uncertainty and the epistemic uncertainty of flutter analysis method selection are considered in the same system to realize quantification. A typical rudder of missile is selected to examine the present method, and the dangerous range of velocity as well as relevant belief and plausibility functions is obtained. The results suggest that the present method is effective in obtaining the lower and upper bounds of flutter probability and assessing flutter risk of structures with limited experimental data of uncertain parameters and the belief of different methods.     

基于交叉熵和空间分割的全局可靠性灵敏度分析

基于失效概率的全局灵敏度分析可以度量各个基本随机变量的不确定性对失效概率的影响程度,对如何降低结构的失效概率具有指导意义。基于交叉熵方法和空间分割提出一种新全局可靠性灵敏度分析方法。该方法采用交叉熵法自适应的确定重要抽样密度函数,有效地回避了传统重要抽样中设计点位置和个数求解的困难。基于评估失效概率所使用的样本,利用空间分割方法计算各个输入随机变量的全局可靠性灵敏度指标,能够提高样本的利用率和计算效率。文中利用一个数值算例和两个工程算例验证了所提方法的计算效率和精度。     

矩独立重要性分析的Kriging代理模型方法

结构重要性分析是结构不确定性分析的重要研究内容之一。作为一种矩独立的重要性测度指标,基于失效概率的重要性测度能够有效地反映输入变量的不确定性对结构失效概率的影响。然而,相比于其他形式的重要性测度,对基于失效概率的矩独立重要性测度进行高效准确的求解仍然存在一定的困难。基于此,结合Kriging代理模型提出了一种高效的矩独立重要性分析的新方法。所提方法首先通过较少的实验设计样本构建能够充分近似真实输入输出响应函数的Kriging代理模型,进而通过数值模拟策略对所构建的Kriging模型进行分析,最终求解得到各个输入变量基于失效概率矩独立重要性测度的主效应和总效应,给出各个输入变量对于失效概率的影响程度排序。相比于现有方法,所提方法由于引入Kriging代理模型极大地降低了对响应函数的计算次数,同时保证了重要性分析的计算精度。两个工程应用实例结果表明了所提方法在计算效率和计算精度方面的优势,体现了所提方法良好的适用性。     

结构系统可靠性及可靠性灵敏度分析的改进子集模拟法

在再生自适应子集模拟(RASS)法的基础上,提出了一种改进的再生自适应子集模拟(MRASS)法以用于结构系统的可靠性及可靠性灵敏度分析。MRASS法继承了RASS法中马尔可夫链再生、延迟拒绝、自适应马尔可夫过程及分量各自采样等优点,并改进了自适应采样过程。MRASS法通过对产生的样本点的接受率与最佳接受率进行比较,有针对性地寻找合适的建议分布方差,提高了抽样的效率。工程算例的数值结果表明:MRASS法相比于RASS法及传统的子集模拟(SS)法,在处理具有高维随机变量、小失效概率及高度非线性特点的结构系统时,有更好的适应性、稳健性及精度。     

Reliability and Sensitivity Analysis of Transonic Flutter Using Improved Line Sampling Technique

The improved line sampling (LS) technique, an effective numerical simulation method, is employed to analyze the probabilistic characteristics and reliability sensitivity of flutter with random structural parameter in transonic flow. The improved LS technique is a novel methodology for reliability and sensitivity analysis of high dimensionality and low probability problem with implicit limit state function, and it does not require any approximating surrogate of the implicit limit state equation. The improved LS is used to estimate the flutter reliability and the sensitivity of a two-dimensional wing, in which some structural properties, such as frequency, parameters of gravity center and mass ratio, are considered as random variables. Computational fluid dynamics (CFD) based unsteady aerodynamic reduced order model (ROM) method is used to construct the aerodynamic state equations. Coupling structural state equations with aerodynamic state equations, the safety margin of flutter is founded by using the critical velocity of flutter. The results show that the improved LS technique can effectively decrease the computational cost in the random uncertainty analysis of flutter. The reliability sensitivity, defined by the partial derivative of the failure probability with respect to the distribution parameter of random variable, can help to identify the important parameters and guide the structural optimization design.     

An efficient algorithm for calculating Profust failure probability

For efficiently estimating the Profust failure probability based on probability input variables and fuzzy-state assumption, a General Performance Function(GPF) expression is established under the strict mathematical derivation for the Profust reliability model. By constructing the GPF,the calculation of the Profust failure probability can be transformed into the calculation of the traditional failure probability. Then various existing methods for the traditional failure probability can be used to estimate the Profust failure probability. Due to the high efficiency of the Adaptive Kriging(AK) model and the universality of the Monte Carlo Simulation(MCS), AK inserted MCS(abbreviated as AK-MCS) has been proven to be an efficient method for estimating the failure probability. Therefore, the AK-MCS combined with the GPF(abbreviated as AK-MCS + GPF)is proposed for estimating Profust failure probability. The proposed method greatly reduces the computational cost while ensuring the accuracy. Finally, four examples are given to validate the proposed AK-MCS + GPF. The results of the examples show the rationality and the efficiency of the proposed AK-MCS + GPF.     

考虑几何分散性的涡轮盘寿命概率分析

针对考虑几何分散性的涡轮盘低循环疲劳（LCF）寿命概率分析中几何参数多、几何随机变量难确定、分布特征获取困难、模型需自动更新及计算成本高的问题,提出几何分散性的概率处理方法:采用试验设计方法对涡轮盘结构所有几何参数进行灵敏度分析,筛选出对应力影响较大的关键几何参数作为随机变量,使用K-S(Kolmogorov Smirnov)方法确定其分布类型和特征参数,最后建立代理模型进行Monte Carlo概率分析.基于此方法,开发出了涡轮盘概率分析系统,在该系统中筛选得到某发动机GH720Li涡轮盘内径、外径、盘缘厚度3个结构参数作为几何随机变量,完成对LCF寿命的概率分析工作得到寿命-可靠度分布曲线.分析结果表明涡轮盘外径对LCF寿命有较大影响.      

同时考虑基本变量和失效域模糊性的广义失效概率数字计算方法

针对基本变量和失效域均具有模糊性的广义可靠性分析问题,提出了一种基于模拟退火优化的高效数字计算方法。在数字模拟的过程中,由模拟退火优化的Metropolis准则逐渐优化提取样本的重要抽样密度函数。由于基本变量和失效状态均含有模糊不确定性,因此所提算法在构造重要抽样函数时考虑了两个因素的影响,其一是基本变量的等价联合概率密度函数,其二是状态变量对模糊失效域的隶属函数,从而使得对广义失效概率贡献大的样本出现的概率较大,提高了抽样效率和计算精度。所提算法在模拟退火逐渐寻优的过程中,充分利用了过程中的信息,进一步提高了计算的效率,算例的结果也表明本文所提方法是合理可行的。     

截断分布下广义失效概率的可靠性分析方法

提出基于随机变量为截断分布和失效域模糊下的结构可靠性分析方法。由于常规的可靠性分析方法不能直接用于截断分布情况下广义失效概率的计算,因此,先将截断分布情况下的可靠性分析模型转化为非截断分布情况下的多模式并联系统模型,在此基础上引入模糊失效域的隶属函数,将截断分布情况下的模糊可靠性分析问题转化为一般的随机可靠性问题,利用Monte Carlo法和重要抽样法求得广义失效概率。给出了方法的实现步骤和原理,并用不同算例验证了方法的合理性与可行性。     

Efficient adaptive Kriging for system reliability analysis with multiple failure modes under random and interval hybrid uncertainty

In the field of the system reliability analysis with multiple failure modes, the advances mainly involve only random uncertainty. The upper bound of the system failure probability with multiple failure modes is usually employed to quantify the safety level under Random and Interval Hybrid Uncertainty(RI-HU). At present, there is a lack of an efficient and accurate method for estimating the upper bound of the system failure probability. This paper proposed an efficient Kriging model based on nume...     

IMPROVED STRATIFIED SAMPLING MONTE CARLO METHOD TO ANALYZE RELIABILITY OF STRUCTURAL SYSTEMS

ＩＭＰＲＯＶＥＤＳＴＲＡＴＩＦＩＥＤＳＡＭＰＬＩＮＧＭＯＮＴＥＣＡＲＬＯＭＥＴＨＯＤＴＯＡＮＡＬＹＺＥＲＥＬＩＡＢＩＬＩＴＹＯＦＳＴＲＵＣＴＵＲＡＬＳＹＳＴＥＭＳＬｉＱｉａｎｇ；ＦｅｎｇＹｕａｎｓｈｅｎｇ（ＤｅｐａｒｔｍｅｎｔＯｆＡｉｒｃｒａｆｔＥｎ...     

基于失效概率的涡轮后机匣全局灵敏度分析

涡轮后机匣是航空发动机安全的关键部件,但是其具有工况复杂、不确定性因素多的缺点。为了探究输入随机变量的不确定性对涡轮后机匣结构失效概率的影响,建立参数化有限元模型进行确定性分析。考虑材料性能、几何参数及外部载荷的不确定性,对涡轮后机匣两种典型失效模式：强度失效以及刚度失效建立极限状态函数;通过构造自适应Kriging 代理模型并结合重要抽样方法评估涡轮后机匣结构失效概率,利用基于失效概率的全局灵敏度方法对涡轮后机匣结构可靠度的不确定性来源进行分析,对各输入随机变量重要性进行排序,构建一种涡轮后机匣全局灵敏度分析框架。结果表明：涡轮后机匣在两种失效模式以及系统失效模式下,发动机推力以及线性膨胀系数对结构失效概率影响最为显著,应对其重点考虑;内外机匣长度以及材料弹性模量对涡轮后机匣结构失效概率影响较小,可对其适当忽略。     

Sensitivity analysis for probabilistic anti-resonance design of aeronautical hydraulic pipelines

In an uncertainty scheme, reliability and global sensitivity analysis is studied in this work, to provide helpful information for probabilistic anti-resonance design of vibration systems.Discussions show that the resonance failure problem can be viewed as a series system, in which input uncertainties are modeled by random variables. In order to quantitatively measure the contributions of input variables to the system reliability, a global sensitivity index is proposed, the properties of which are also discussed. Then the proposed index is tested with an aeronautical hydraulic pipeline system, which is under the excitation of pump vibration and at a risk of resonance failure. Sensitivity results under different failure criteria and variation coefficients are obtained and studied, from which significant and insignificant input variables can be identified.The proposed method provides a relatively new insight for anti-resonance design of engineering structures.     

Research on the Plane Multiple Cracks Stress Intensity Factors Based on Stochastic Finite Element Method

Taylor stochastic finite element method (SFEM) is applied to analyze the uncertainty of plane multiple cracks stress intensity factors (SIFs) considering the uncertainties of material properties, crack length, and load. The stochastic finite element model of plane multiple cracks are presented. In this model, crack tips are meshed with six-node triangular quarter-point elements; and other area is meshed with six-node triangular elements. The partial derivatives of displacement and stiffness matrix with respect to all the random variables obtained by Taylor SFEM are derived. Meanwhile, the mean and variance expressions of SIF under uncertain factors are also derived. Parallel double-crack illustrative example of using the proposed method is given. The calculation results indicate that the uncertainty of SIF is influenced by the distance of the cracks, the smaller the distance of cracks is, the greater the SIF uncertainty is, and the uncertainties of material elastic modulus, load and crack length markedly affect the uncertainty of SIF; and the Poisson ratio of material has little influence. Among the variables, the elastic modulus has the greatest effect on SIF uncertainty. The next is external load. The crack length has lower effect on SIF uncertainty than both the elastic modulus and external load.     

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

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

基于密度权重的可靠性灵敏度分析方法简

 为了提高可靠性灵敏度求解数字模拟法的效率,提出了一种变量空间确定性低偏差均匀抽样与样本点处联合概率密度函数构造权重相结合的方法,来估计可靠性灵敏度。该方法通过均匀样本点处联合概率密度函数的权重保证了可靠性灵敏度的估计值收敛于真值,而由低偏差抽样代替原问题中的联合概率密度抽样则可以保证更低的误差阶以及在小失效概率条件下抽得的样本有更高的可能性落入失效域,从而保证了所提方法具有更高的收敛速度。另外,所提方法可以采用与独立变量相同的步骤来估计相关变量情况下的可靠性灵敏度,计算简便,适用范围广。算例充分证明了所提方法的优越性。     

基于密度权重的可靠性灵敏度分析方法

为了提高可靠性灵敏度求解数字模拟法的效率,提出了一种变量空间确定性低偏差均匀抽样与样本点处联合概率密度函数构造权重相结合的方法,来估计可靠性灵敏度。该方法通过均匀样本点处联合概率密度函数的权重保证了可靠性灵敏度的估计值收敛于真值,而由低偏差抽样代替原问题中的联合概率密度抽样则可以保证更低的误差阶以及在小失效概率条件下抽得的样本有更高的可能性落入失效域,从而保证了所提方法具有更高的收敛速度。另外,所提方法可以采用与独立变量相同的步骤来估计相关变量情况下的可靠性灵敏度,计算简便,适用范围广。算例充分证明了所提方法的优越性。