具有随机不确定性的机翼颤振优化

针对不确定性参数在颤振优化分析中的影响,提出了一种含颤振约束的、具有随机不确定性的机翼气动弹性优化设计方法.在MSC.NASTRAN平台上建立了结构的动力学模型、采用ZAERO程序进行了非定常气动力计算与颤振分析,用响应面方法计算了隐式失效功能函数,以重量最小为设计目标,对包含颤振概率约束的复合材料机翼进行了优化设计....     

可靠性灵敏度分析的一种新方法

基于极限状态函数矩估计的失效概率计算,提出一种新的可靠性灵敏度分析方法。推导极限状态函数的矩对基本变量分布参数的偏导数,并进而利用失效概率与极限状态方程矩的关系,推导失效概率对基本变量分布参数的偏导数,从而得到可靠性灵敏度。与改进一次二阶矩可靠性灵敏度分析方法相比,所提方法不用求极限状态方程的设计点,因而不需用到极限状态函数对基本变量的梯度函数,适用于隐式极限状态方程的可靠性灵敏度分析,算例结果也充分显示所提方法的合理性和精度。     

基于响应面法的可靠性稳健设计优化

针对极限状态函数和结构性能函数为隐式的结构可靠性稳健设计优化问题,将结构性能稳健性、可靠性灵敏度稳健性与6σ稳健设计思想相结合,建立了一种基于响应面法的可靠性稳健设计优化模型并给出了具体算法.采用多项式拟合结构性能函数和极限状态函数,推导了基于响应面函数的1阶可靠度指标以及性能函数均值和均方差的计算方程;将结构可靠度作为基本约束条件,采用序列二次规划法得到稳健设计优化结果.运用该方法对某涡轮叶片型线进行设计优化,经过优化后,静、动叶的动能效率均值分别提高1.5%和6.4%,动叶片最大应力均值下降7.9%,动叶叶身最大变形均值下降5.6%,动叶片结构强度可靠度由初始的91.3%提高到了99.9%.      

非线性壁板颤振分析

利用一种改进的计算流体力学与计算结构动力学(CFD/CSD)耦合方法研究了由气动和结构几何非线性引起的壁板颤振问题。在非定常气动力计算中,考虑了通量分裂格式、隐式时间推进方法和几何守恒律;二维和三维壁板的结构几何非线性建模则采用了有限元的协同旋转理论,并利用一种近似能量守恒算法求解结构的非线性响应。流场和结构求解器采用二阶松耦合方法联立求解,并将其应用于壁板在超声速、跨声速和亚声速的颤振计算中。当考虑结构几何非线性和气动非线性时,出现了典型的极限环振荡现象,并对颤振边界和极限环振荡幅度进行了对比分析。     

A reduced order model for coupled mode cascade flutter analysis

A Reduced Order Model (ROM) based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper. The unsteady aerodynamic model is established by a system identification technique combined with a set of Aerodynamic Influence Coefficients (AIC). Subsequently, the aerodynamic model is encoded into the state space and then coupled with the structural dynamic equations, resulting in a ROM of the cascade aeroelasticity. The cascade flutter can be determined by solving the eigenvalues of the ROM. Bending-torsional coupled mode flutter analysis for the Standard Configuration Eleven (SC11) cascade is used to validate the proposed method.     

操纵面对跨声速机翼气动弹性特性的影响

 运用基于非定常CFD的气动力辨识技术，得到跨声速非定常气动力降阶模型。耦合结构动力学方程，建立了基于状态空间的跨声速气动弹性分析模型。分析了典型三自由度二元机翼的颤振边界，分析结果与CFD/CSD直接耦合方法吻合。然后研究了操纵面结构参数（固有频率和重心位置）对跨声速气动弹性特性的影响。研究发现,一些传统的结构设计准则和颤振排除技术未必适用于跨声速状态;操纵面偏转模态常常成为诱发跨声速颤振的主要模态;经典的质量平衡技术可能会降低跨声速气动弹性系统的稳定性。     

两种跨声速气动弹性问题分析研究

在非结构运动网格基础上,采用中心有限体积法进行空间离散和双时间方法进行时间推进求解非定常欧拉方程.通过与气动力方程的联立求解,在时域内用四步龙格-库塔方法求解结构运动方程.分析和研究了二维嗡鸣和三维机翼颤振这两种跨声速非线性气动弹性问题.二维嗡鸣问题的研究考虑了翼面-舵面系统的缝隙间网格运动、缝隙对嗡鸣的影响和扰流片对嗡鸣的抑制.耦合多自由度Lagrange结构运动方程数值模拟了三维机翼的颤振问题.通过跨声速标模算例AGARD445.6机翼的颤振计算,计算的颤振临界速度与实验值有5%左右的误差,验证本方法的正确性.由于本方法是在对外形具有良好普适性的非结构动网格基础上完成的,具有良好的工程实践价值.     

Neural Networks Combined with Importance Sampling Techniques for Reliability Evaluation of Explosive Initiating Device

Concerning the issue of high-dimensions and low-failure probabilities including implicit and highly nonlinear limit state function, reliability analysis based on the directional importance sampling in combination with the radial basis function (RBF) neural network is used, and the RBF neural network based on first-order reliability method (FORM) is to approximate the unknown implicit limit state functions and calculate the most probable point (MPP) with iterative algorithm. For good efficiency, based on the ideas that directional sampling reduces dimensionality and importance sampling focuses on the domain contributing to failure probability, the joint probability density function of importance sampling is constructed, and the sampling center is moved to MPP to ensure that more random sample points draw belong to the failure domain and the simulation efficiency is improved. Then the numerical example of initiating explosive devices for rocket booster explosive bolts demonstrates the applicability, versatility and accuracy of the approach compared with other reliability simulation algorithm.     

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

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

基于谐波平衡法和 V-g 法的高效颤振预测分析

提出了一种高效的颤振预测分析方法,将之应用于跨声速颤振边界分析及结构参数影响研究中。本方法基于频域颤振分析 V-g 方法,为气弹系统提供一定的人工阻尼使之保持简谐运动状态,从而将结构动力学方程转换到频域内。然后通过高效的谐波平衡法得到一系列简谐运动频率下的气动力描述函数矩阵,结合频域结构方程,将气弹系统的稳定性问题转化为广义特征值求解。结果表明：本方法计算得到的颤振边界与高精度的时间推进方法非常吻合,分析效率有了明显的提升,而且当结构参数发生变化后,只需进行若干次广义特征值求解即可得到新的颤振边界,无需像时间推进方法一样开展大量的气动/结构耦合数值模拟。     

Flutter analysis of an airfoil with bounded random parameters in incompressible flow via Gegenbauer polynomial approximation

The effects of parameter uncertainty on the flutter characteristics of a two-dimensional airfoil in an incompressible flow were investigated through Gegenbauer polynomial approximation. The uncertain parameters, such as the linear and cubic pitch stiffness coefficients are modeled as bounded random variables with λ-PDFs (probability density functions). With the aid of Gegenbauer polynomial approximation, the two-dimensional stochastic airfoil system is transformed at first into its equivalent deterministic one. Then the Hopf-bifurcation point is determined through the equivalent deterministic system, and the onset of the flutter, together with the flutter frequency against the probability density distribution parameter and the intensity of the random variable is explored. In addition, the ranges of the peak pitch and plunge responses against the flight speed and the PDFs of the peak pitch response are obtained. Numerical results show how the probability density distribution parameters and the intensities of random parameters affect the flutter onset speed, flutter angular frequency, the upper and lower boundaries of peak responses, and the PDFs of peak responses.     

Control law design for transonic aeroservoelasticity

Existing computational transonic aeroservoelastic researches focus on directly coupling the structural dynamic equations, CFD solver and servo system in time domain, study the effect of the given feedback control laws on the responses of the aeroelastic system. These works have not involved the design of the flutter active control law. The non-linearity of transonic flow brings great difficulties to aeroservoelastic analysis and design. Recent research of the unsteady aerodynamic reduced order models (ROM) based on CFD provides a challenging approach for transonic aeroservoelastic analysis and design. Coupling the structural state equations with the aerodynamic state equations of the wing and the control surface based on the ROM, we construct a transonic aeroservoelastic model in state-space. Then the sub-optimal control method based on output feedback is used to design the flutter suppressing law. The study first demonstrates the open loop of the Benchmark Active Controls Technology (BACT) wing. The computational results of the CFD direct simulation method and the ROM analysis method are both agree well with the experimental data. Then both the closed loop time responses and the flutter results by ROM technique are compared with those of numerical aeroservoelastic simulation based on Euler codes to validate the correctness of the design method of the control law and aeroservoelastic analysis method. An increase of up to 20% of the speed index can be achieved by the control law designed by sub-optimal control method for this model.     

考虑状态与基本变量随机模糊性的广义响应面可靠性分析

针对复杂结构、机构可靠性分析中的隐式极限状态方程,建立了可以同时考虑基本变量和失效一安全状态随机模糊双重不确定性因素的广义响应面可靠性分析方法。该方法采用模糊基本变量的等效随机化变换,在不改变基本变量模糊分布的情况下,建立了基于有限元分析的复杂结构、机构多个随机一模糊基本变量情况下的广义响应面法,利用广义响应面函数和模糊随机事件概率的计算公式,定义了复杂结构、机构广义失效概率的计算公式,并给出了多个模式情况下体系同时考虑状态和基本变量随机模糊性的广义失效概率计算公式。在弹性连杆机构强度刚度多模式广义可靠性分析的应用算例说明了所提方法的合理性。     

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

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

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

跨声速颤振计算方法研究

采用改进的弹簧近似方法,以ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程为控制方程,耦合结构运动方程,计算了Isogai机翼模型的颤振边界。气动力求解采用双时间推进,外时间为物理时间,与结构运动方程同步;结构运动方程求解采用Rayleigh-Ritz方法,时间推进也采用双时间法,节省了大量计算时间。研究表明,计算结果与参考文献提供的结果吻合良好,证明了所提出的跨声速颤振计算方法的准确性。     

跨声速风扇的流固耦合的颤振分析(英文)

发展了时间推进的叶片颤振的数值方法,采用了时间推进求解流体和固体相互作用的过程。气动模型是基于求解三维雷诺平均的N-S方程,采用了多块结构化网格的有限体积格式,对流和耗散通量的计算使用了二阶迎风格式和中心格式。在变形的动网格上流体运动的守恒型方程的求解采用双时间步,隐式格式和多重网格方法。叶片振动采用了振型叠加的线性气动弹性模型。该方法在气动弹性标准算例4进行了验证,并用于求解跨音速风扇的颤振问题。     

Limit cycle oscillation control of wing with static output feedback control method

This paper presents a static output feedback controller (SOFC) for aeroelastic control of a cantilevered rectangular wing in low subsonic flow. For this purpose, an optimal formulation of this control method is developed, and a solution method is proposed for the related matrix equation. This optimal solution is obtained by solving combined Lyapunov and Riccati equations. At first these equations are transformed into a set of nonlinear algebraic equations and then are solved with iterative Newton–Raphson?s method. This solution method is applicable to the full state feedback case. The controller is designed to extend flutter boundary and suppress limit cycle oscillation (LCO) of a low aspect ratio rectangular nonlinear structural wing. This structural nonlinearity is given by Von Karman plate theory. Both full and reduced order aerodynamic models are examined based on the modified vortex lattice theory. Results show combination of SOFC with reduced order aerodynamic model would be an effective choice for aeroelastic stabilization, and this controller has a very comparable result with linear quadratic regulator (LQR).     

大后掠翼前缘涡对其颤振特性的影响

大迎角三角翼的前缘涡不仅可以改善其气动力特性,也会显著影响机翼的气动弹性特性。运用基于Euler方程的非定常气动力降阶模型（ROM）方法,耦合结构运动方程,在状态空间内建立了气动弹性分析模型,研究了70°削尖三角翼的大迎角颤振特性。研究结果显示前缘涡对该机翼颤振特性的影响不可忽略。颤振速度随迎角的增加而大幅降低,迎角α=20°时的颤振速度比α=0°时降低了22％。发现了颤振特性随迎角变化时出现的不连续现象,并揭示了该现象是由于系统颤振分支随着静态迎角的增加发生转移所致。     

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

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