航空发动机限寿件疲劳可靠度计算新方法

针对目前的航空发动机限寿件（ELLP）疲劳可靠性分析中的小失效概率事件以及其极限状态函数具有较强非线性的特点，提出了一种具有自更新机制的半径外自适应重要抽样（AUMCROAIS）疲劳可靠性分析方法。该方法首先利用蒙特卡罗自适应重要抽样（MCAIS）快速逼近真实设计验算点（MPP）附近，随后以近似设计验算点为中心进行极坐标抽样，并依次构造主动学习函数，对近极限状态函数和抽样半径进行最优选取，从而实现最优抽样半径的更新，通过不断的更新确定出最优抽样半径，加速失效概率计算的收敛。本方法提高了设计验算点的收敛速度同时保证了计算精度，解决了小失效概率事件以及强非线性极限状态函数可靠度计算难题，最后以某型发动机压气机轮盘为对象应用本方法，并与传统的蒙特卡罗仿真（MCS）方法、蒙特卡罗半径外自适应重要抽样法（MCROAIS）和一阶可靠性方法（FORM）进行了对比，验证了本方法的高效率、鲁棒性和仿真精度。     

相关变量模糊可靠性灵敏度分析的线抽样方法

 依据失效域具有模糊性时模糊失效概率的定义,提出了相关变量模糊可靠性灵敏度的分析方法。针对线性功能函数、正态基本变量和正态型隶属函数情况,推导了相关变量的模糊可靠性灵敏度计算的解析表达式。对于工程中的一般情况,给出了可靠性灵敏度分析的数字模拟方法。尽管数字模拟法适用范围广,但该方法的效率较低,尤其是针对高维和小失效概率问题。为了解决数字模拟方法效率低的问题,提出了相关变量模糊可靠性灵敏度分析的线抽样方法。通过离散模糊失效概率积分区域,建立了相关变量模糊可靠性灵敏度与离散区域随机可靠性灵敏度的关系,进而可以利用相关变量随机可靠性灵敏度分析的线抽样方法求得模糊可靠性灵敏度。相关变量模糊可靠性灵敏度分析线抽样方法的基本原理、计算公式及实现步骤被详细给出,文中算例充分验证了其精度高、收敛快及适用于高维和小失效概率等优点。     

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

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

两种基于方差的全局灵敏度分析W指标改进算法

在全局灵敏度分析(SA)中,基于方差的灵敏度分析指标(包括Sobol指标和W指标)应用广泛。其中Sobol指标是将输入随机变量固定于特定点时,求得其对输出响应量的平均影响;而W指标求解当输入随机变量在各自分布区间上缩减变化时,输入变量对输出响应量的影响程度。相比Sobol指标,W指标所反映的信息更加全面。但目前对W指标的求解方法还比较欠缺,双层重复抽样蒙特卡罗(DLRS MC)方法和双层一次抽样蒙特卡罗(DLSS MC)方法是两种传统的求解方法。针对W指标的求解问题,提出了两种新算法:改进的蒙特卡罗模拟(AMCS)和基于稀疏网格积分(SGI)的方法。AMCS只需抽取一组样本便可计算出所有变量的各阶W指标,由于该方法是通过筛选策略来计算条件区间上的方差,避免了DLSS MC法中由于小数取整带来的计数误差,从而提高了计算W指标的精度。基于SGI的方法则利用稀疏网格积分来计算三重矩进而得到W指标,由于该方法继承了稀疏网格积分的高效性,因而进一步提高了W指标的计算效率。最后,给出了两个数值算例和一个工程算例,用于验证所提方法求解W指标的准确性和高效性。     

一种改进的航空发动机结构概率风险评估方法

针对航空发动机适航条款FAR33.75中关于发动机限寿件(ELLP)结构失效概率要求,提出了一种基于Kriging和蒙特卡罗半径外重要抽样(MCROIS)混合的结构概率风险评估方法。该方法针对ELLP高维、小失效概率事件以及极限状态函数为隐式、高度非线性的特点,利用Kriging元模型模拟隐式极限状态函数,然后通过主动学习迭代算法,计算最优点(MPP,最接近设计验算点的样本点),更新实验设计(DOE)并提高Kriging元模型的模拟精度。在此基础上,利用Kriging元模型确定最优抽样半径,构造半径外重要抽样密度函数,在最优抽样半径确定区域进行抽样,通过构造主动学习函数,使样本点更多落在抽样半径确定的球区域附近,加速失效概率计算的收敛,并构建了ELLP风险概率模型,解决了高维、小失效概率事件以及隐式、非线性极限状态函数的发动机结构概率风险评估难题,以某型发动机低压压气机轮盘为应用示例,与传统的蒙特卡罗仿真(MCS)方法进行了对比,验证了该方法的高效率、鲁棒性和仿真精度。     

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

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

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

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

高维小失效概率下的改进线抽样方法

 针对隐式极限状态方程的可靠性分析，提出了一种基于失效域样本模拟的高效线抽样法。该法采用马尔可夫链来快速得到失效域中的条件样本，利用这些失效域中的样本可以得到线抽样的重要方向，并且这些样本还可以作为线抽样的样本，来得到隐式极限状态方程失效概率的估计值。由于该法可以较准确地得到线抽样方法的重要方向，因此线抽样的效率可以得到提高，另外失效域的样本又可以作为线抽样的样本，从而进一步降低所提方法的计算工作量。算例表明所提算法的计算精度和效率均高于传统线抽样方法。     

基于双重极值响应面法的叶盘联动可靠性分析

为了更准确地研究变形和应力对航空发动机叶盘性能的影响,在湍流模型理论和叶盘流场模型的基础上,提出了一种可靠性分析方法即双重极值响应面法(Dual Extremum Response Surface Method,DERSM)。在分析中考虑气体压力和离心力作用,合理选取随机输入变量,得到叶盘结构的总体最大变形和最大应力,并作为输出响应;对随机变量进行抽样,通过抽样点拟合各自的极值响应面方程;结合蒙特卡洛法对两个极值响应面模型进行联动抽样,计算出了航空发动机叶盘的可靠性概率。实例仿真表明:当叶盘最大允许变形量为4.2mm,最大许用应力为7.5×10~2MPa时,其可靠性概率为98.61%。与极值响应面法(Extremum Response Surface Method,ERSM)相比,DERSM在保证计算精度的前提下,计算效率提高23.52%,验证了DERSM在叶盘可靠性分析方面的有效性。     

基于蒙特卡罗模拟的飞机结构可靠度分析

飞机可靠度的分析要综合考虑各种因素的影响,由于各因素的影响较为复杂,很难通过解析的方法来计算飞机的可靠度.随着计算机速度的不断提升,通过蒙特卡罗模拟计算飞机结构的可靠性逐渐变得可行.本文基于蒙特卡罗法,将损伤、工作环境、检修周期等作为核心要素,提出了面向工程应用的飞机结构可靠度分析方法,最终得到了飞机全机结构失效概率的计算方法.本方法在计算机上得到了初步实现.     

基于MBD技术检验规程在飞机数字化装配中的应用

针对基于MBD技术的飞机数字化装配技术的大量应用,分析了传统检验模式在MBD技术环境中存在的问题和不足,提出基于MBD技术的检验规程,结合实际工作探讨了基于MBD技术编制检验规程的方法,并给出了解决其关键技术问题的方法。     

有限元的一种并行算法及实现

当前,由于有限元素法大都使用串行算法,处理时间较长。本文根据有限元素法的特点和采用的并行计算机结构,提出了一种并行有限元算法,能大大提高其处理速度,加快处理周期。且算法简单,实现方便。     

低反动度高负荷超声速轴流压气机气动设计方法

对超声速轴流压气机的发展历程进行了总结,分析了两类传统超声速压气机内部流动特点及其所存在的流动问题,并对未来超声速压气机的发展提出了展望.针对超声速压气机内部流动特点,提出了一种新的低反动度高负荷超声速轴流压气机气动设计原理.该原理可有效避免在动叶中进行流动控制,同时降低动叶出口绝对马赫数,并结合附面层抽吸控制静叶栅内部流动,最终实现级的高负荷气动设计.利用该原理进行设计验证,三维黏性数值模拟结果表明:在叶尖切线速度为360m/s的前提下,实现了2.3的级压比,级效率为86.5%.      

Structural optimization for multiple structure cases and multiple payload cases with a two-level multipoint approximation method

This paper is to address structural optimization problems where multiple structure cases or multiple payload cases can be considered simultaneously. Both types of optimization problems involve multiple finite element models at each iteration step, which draws high demands in opti-mization methods. Considering the common characteristic for these two types of problems, which is that the design domain keeps the same no matter what the structure cases or payload cases are, both problems can be formulated into the unified expressions. A two-level multipoint approxima-tion (TMA) method is firstly improved with the use of analytical sensitivity analysis for structural mass, and then this improved method is utilized to tackle these two types of problems. Based on the commercial finite element software MSC.Patran/Nastran, an optimization system for multiple structure cases and multiple payload cases is developed. Numerical examples are conducted to ver-ify its feasibility and efficiency, and the necessity for the simultaneous optimizations of multiple structure cases and multiple payload cases are illustrated as well.     

大型广义特征值问题的动力降阶迭代法

 本文在Guyan、Kuhar等人工作基础上,提出通过动力降阶和迭代,对大型广义特征值问题实施逼近求解的一个新方法,以能对大型、复杂动力系统实施动力分析。既解决计算机容量不足,又保证精度,亦不耗时过多。通过实际使用,证明该方法的有效性。     

Implicit discontinuous Galerkin method on agglomerated high-order grids for 3D simulations

High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is inves-tigated in combination with discontinuous Galerkin (DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry represen-tation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomer-ated high-order mesh works well with DG method in dealing with flow problems of curved geome-tries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant effi-ciency improvements compared to its explicit counterparts.     

混合场积分方程在开放结构散射分析中的应用

针对复合开放闭合结构电磁散射问题的矩量法求解,通过修改混合场积分方程的建立方式,将混合场积分方程扩展应用到这类开放结构,大大提高了求解效率,数值算例证明了该方法的精度和效率。     

Electrode design using revolving entity extraction for high-efficiency electric discharge machining of integral shrouded blisk

The integral shrouded blisk provides better performance with minimum weight, but its semi-open structure presents problems for its production. Manufacturing processes of these components require a removal of about 70%-90% of material from their blanks. Multi-axis electrical discharge machining (EDM) is commonly used for these processes, but its poor efficiency cannot meet the requirements of mass production. Strong flushing assisted high-current discharge is able to improve the machining efficiency. In this paper, this method was applied to manufacture the integral shrouded blisk. An electrode design method was proposed. Taking the largest revolving entity inside the flow channel as the base geometry, this laminated arc-shaped hollow electrode meets the requirements such as high pressure flushing, tool wear compensation, and easy to be made. A 4-axis linkage machining tool path with planetary motion was proposed. Taking an integral guide vane ring as an example, it has been experimentally verified that the time consumption for each channel using this method was reduced to a half of the ordinary EDM method, while the finished surface quality remains same.     

一种高效的隐式间断Galerkin方法研究

基于线性化处理,在时间方向上对间断Galerkin方程进行了隐式离散,从整体上对迭代过程进行了合理的优化,并以此求解了计算流体力学中的二维Euler方程。其中,LU-SGS方法得到了进一步的推广,被用来高效求解隐式格式对应的大型稀疏线性系统。数值实验表明,无论对于亚声速问题还是跨声速问题,该格式都是无条件稳定的;与显式的Runge-Kutta间断Galerkin格式相比,当残值下降到相同量级时,隐式格式所需的迭代步数和CPU时间均在很大程度上得到了减少。     

飞机主操纵系统疲劳试验载荷监控方法研究——应变数据分析法

针对某机全机主操纵系统疲劳试验中动态操纵情况下的应变测量以及试验载荷监控问题进行了分析和讨论,并对基于应变测量数据分析的操纵系统疲劳试验载荷监控方法的确定及其应用效果进行了阐述。在操纵系统疲劳试验过程中的实际监控结果表明,该方法对副翼/平尾/方向舵操纵系统多种试验状态下的试验载荷监控实用有效。