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.     

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.     

基于降阶模型的气动弹性主动控制律设计

流体/结构耦合数值模拟是目前解决复杂气动弹性问题精度最高的方法。但由于计算效率比较低,模型阶数过高,不能直接用于气动弹性系统的主动控制律设计。为了对主动控制系统设计提供高效高精度状态空间模型,研究了气动弹性系统的时域正则正交分解(POD)/降阶模型(ROM)方法,并引入平衡截断(BT)技术进一步降低时域POD/ROM的阶数,从而有效克服了时域POD/ROM阶数过高的缺点。在此基础上建立了基于POD-BT/ROM的气动伺服弹性降阶方程。以AGARD445.6机翼为例,说明了时域POD/ROM建模的各个细节,并将其用于气动弹性主动控制律的设计。计算结果表明,POD/ROM具有接近计算流体力学(CFD)/计算结构动力学(CSD)耦合计算的精度,同时又大大提高了计算效率约1～2个量级,是一种高精度高效率的气动弹性主动控制系统设计工具。     

基于CFD/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

基于非定常气动力辨识技术的气动弹性数值模拟

选择离散型输入输出差分模型,运用最小二乘方法进行非定常气动力建模,并将辨识得到的降阶模型用于气动弹性的数值模拟。1个马赫数下的颤振临界点的计算仅需调用一次非定常流场求解器。计算精度保持与非定常欧拉方程计算方法相当的同时计算效率提高了1~2个量级。计算了跨声速具有S型颤振边界的气动弹性标准算例－Isogaiwing和三维气动弹性标模算例AGARD445.6,辨识模型计算边界与非定常Euler方程计算结果吻合。证明非定常气动力辨识技术可以提供高效的高精度的气动弹性分析。     

机翼/外挂系统的颤振主动抑制研究

 本文对颤振主动抑制控制律进行了研究。研究对象为一小展弦比带外侧导弹的机翼颤振模型,模型具有外侧后缘控制面。依据该模型的全部动力特性和刚度特性,以最优控制为基础,采用动态补偿器方法,设计了两阶控制律。对该控制律进行了风洞实验验证。实验结果表明:颤振临界速度提高了14％以上。理论计算结果与实验结果一致。     

数字式颤振主动抑制系统的研究

介绍了数字式颤振主动抑制系统的控制律综合、风洞实验及其工程特性。控制律的设计采用连续域离散化的方法。对４种不同控制律进行了风洞实验,验证了理论计算与实验设计的正确性。实验中控制律用微机和工业控制机联合实现;还根据工程应用的观点提出了颤振主动抑制系统应该具有的特性。对一个与实际飞机动力相似的机翼／外挂系统的研究表明,颤振速度提高了２０％左右,并且系统的工程特性也是基本可行的。     

飞行控制律对体自由度颤振特性影响试验

飞翼飞机易发生刚体短周期模态与机翼低阶弯曲模态耦合所致的体自由度颤振。飞行控制系统对飞机的短周期模态特性影响很大,因此考虑飞行控制系统的闭环体自由度颤振特性值得进一步研究。针对自主设计的颤振模型开发了相应的俯仰姿态保持控制律,综合运用风洞试验和仿真计算开展了相关研究,获得了不同刚体自由边界条件下的开环/闭环体自由度颤振特性,研究了闭环增益对体自由度颤振特性的影响规律,简要分析了影响机理。试验和仿真计算结果共同表明：俯仰姿态保持控制律明显地改变了俯仰模态阻尼的原有走势,闭环后的体自由度颤振特性变化明显。以开环颤振速度为基准,采用较小的比例回路增益K_P或较大的微分回路增益K_D,飞行控制律能增加飞行器俯仰阻尼,提高体自由度颤振速度,反之飞行控制律将导致颤振速度降低。就本文控制律而言,当K_P<0.07或K_D>0.2时俯仰姿态保持控制律能起到抑制体自由度颤振的作用。     

考虑机翼几何非线性的气动弹性建模与分析

大展弦比柔性机翼结构重量轻、气动效率高,广泛应用于高空长航时无人机(UAVs)。飞行过程中,这类机翼在气动力作用下发生大变形,线性结构模型不再适用,需要建立考虑几何大变形的结构模型。采用牛顿力学方法推导了考虑结构几何非线性的机翼结构动力学模型,该方法推导过程简洁、物理意义明确,可以与Hodges基于哈密顿原理的推导方法相互补充,相互验证。为了能够更准确地求解大展弦比柔性机翼的非定常气动力,建立了能够考虑机翼三维效应且适用于机翼空间大变形的非定常气动力模型。基于建立的非线性结构模型和非定常气动力模型,采用松耦合方法建立了非线性气动弹性模型,并通过算例验证了气弹模型的准确性。研究结果表明,大展弦比柔性机翼颤振速度对来流迎角和机翼的展长均较为敏感;当来流速度大于颤振速度时,由于几何非线性,机翼振动并未发散而是形成稳定的极限环振荡(LCO);随着来流速度进一步增加,机翼再次穿过临界稳定点,由不稳定系统变为稳定系统,直到随着速度的增加系统再次达到临界稳定状态。     

RESEARCH ON DIGITAL ACTIVEFLUTTER SUPPRESSION

ＲＥＳＥＡＲＣＨＯＮＤＩＧＩＴＡＬＡＣＴＩＶＥＦＬＵＴＴＥＲＳＵＰＰＲＥＳＳＩＯＮＺｏｕＣｏｎｇｑｉｎｇ，ＣｈｅｎＧｕｉｂｉｎ（５－ｔｈＤｅｐａｒｔｍｅｎｔ，ＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，...     

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

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

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).     

基于H_∞/H_2方法的弹性飞行器颤振抑制系统设计

在考虑机翼弹性的情况下,研究大展弦比高空长航时飞行器的主动控制技术,以达到颤振抑制的效果。为了解决飞行器模型阶数过高难以进行控制律设计的问题,采用改进的平衡截断方法对模型进行降阶处理。同时针对机翼弹性所带来的严重气动弹性问题,基于静态输出反馈H_∞/H_2控制方法对降阶后的模型设计了颤振抑制控制律。仿真结果表明,飞行器在海平面高度下的颤振临界马赫数提高了128%,说明H_∞/H_2控制器起到了很好的主动控制效果。     

翼面结构/颤振主动控制律一体化设计

在气动弹性领域内对结构 /控制一体化设计作了初步的探索与研究。针对采用加速度传感器的气动伺服弹性模型 ,建立了在拉普拉斯域中颤振速度对结构设计变量及控制设计变量的解析导数表达式。同时也建立了奇异值对结构设计变量的解析导数表达式。对一矩形机翼以输出反馈的形式进行了结构 /控制一体化设计。目标函数为结构重量最轻 ,满足闭环系统颤振速度和以奇异值为基础的鲁棒稳定性指标约束 ,取得了较好的计算结果     

基于CFD降阶模型的阵风减缓主动控制研究

飞行器飞行时会受到大气紊流的影响,降低飞行品质。阵风减缓控制是改善飞行器飞行性能的关键技术。现有的阵风响应分析多以离散阵风为研究对象,对更加真实描述大气紊流的连续型阵风时域分析关注较少。采用成形滤波器方法将频域形式给出的大气紊流信号转换为时域信号。在跨声速区域内,利用系统辨识技术,基于计算流体力学(CFD)方法建立阵风激励下的气动载荷状态空间降阶模型(ROM)。为方便控制器设计,借助平衡模态法进行模型的进一步降阶。使用模型预测控制(MPC)算法通过控制操纵面偏转实现阵风减缓主动控制。以AGARD445.6标模作为仿真算例,验证基于ROM设计的阵风减缓控制律的有效性。仿真结果表明,在跨声速飞行状态下,模型预测控制器能够在满足操纵面偏转范围的约束下,对连续阵风激励下的翼根弯矩输出进行有效抑制。     

Flutter control of a composite plate with piezoelectric multilayered actuators

Active flutter velocity enhancement scheme is presented for lifting surfaces, employing Linear Quadratic Gaussian based multi-input multi-output controller with multilayered piezoelectric actuators. To numerically test the developed concept, a composite plate wing, surface bonded with eight piezoelectric bender actuators and sensors has been considered. A modal flutter control model is formulated in state-space domain using coupled piezoelectric finite element procedures along with unsteady aerodynamics and optimal control theory. The bending – torsion flutter instability has been actively postponed from 44.13 to 55.5 m/s using the energy imparted by the multilayered piezoelectric actuators. As the power requirement by these actuators is comparatively very low with respect to stack actuators, they can be employed in an integrated form to develop active lifting surfaces for real time applications.     

气动弹性系统的阵风减缓与颤振主动抑制

针对Ｊ８飞机模型，研究了阵风减缓与颤振主动抑制的综合控制问题，应用现代控制理论设计控制系统，分别对机翼／外挂系统模型作开环和闭环分析。由数字式控制实现了阵风减缓与颤振主动抑制的风洞实验，风洞实验结果表明：设计控制律具有抑制颤振和减缓阵风响应的双重功能。     

Active flutter suppression of a multiple-actuated-wing wind tunnel model

In this study, a multi-input/multi-output（MIMO） time-delay feedback controller is designed to actively suppress the flutter instability of a multiple-actuated-wing（MAW） wind tunnel model in the low subsonic flow regime. The unsteady aerodynamic forces of the MAW model are computed based on the doublet-lattice method（DLM）. As the first attempt, the conventional linear quadratic-Gaussian（LQG） controller is designed to actively suppress the flutter of the MAW model. However, because of the time delay in the control loop, the wind tunnel tests illustrate that the LQG-controlled MAW model has no guaranteed stability margins. To compensate the time delay, hence, a time-delay filter, approximated via the first-order Pade approximation, is added to the LQG controller. Based on the time-delay feedback controller, a new digital control system is constructed by using a fixed-point and embedded digital signal processor（DSP） of high performance. Then, a number of wind tunnel tests are implemented based on the digital control system.The experimental results show that the present time-delay feedback controller can expand the flutter boundary of the MAW model and suppress the flutter instability of the open-loop aeroelastic system effectively.     

高超声速飞行器气动弹性力学研究综述

高超声速飞行器设计上的特点带来了一系列的气动弹性新问题。本文回顾高超声速飞行器气动弹性研究的历史与现状,着重介绍和分析了高超声速非定常气动力计算方法、热环境下的气动弹性问题、壁板颤振、推力影响下的气动弹性稳定性问题以及气动推进/气动弹性耦合的多学科交叉问题,相关的主动控制方法的研究进展亦有所介绍。在已有气动弹性问题研究发展的基础上,提出了高超声速飞行器在气动弹性领域需要解决和关注的若干问题,包括高超声速气动弹性试验、燃料消耗的质量变化对于飞行器气动弹性特性的影响以及气动弹性力学与飞行力学综合等方面。     

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 e...