控制面间隙对非线性二元机翼气动弹性响应的影响

采用Theodorsen非定常气动力建立同时含有俯仰立方非线性和控制面间隙非线性二元机翼的运动方程,并以状态空间形式描述。基于状态依赖Riccati方程控制方法,设计了非线性颤振控制律。综合运用Runge-Kutta数值方法与Henon方法,研究了控制面间隙对系统开环/闭环响应的影响。其中Henon方法用以准确快速地确定间隙非线性的转折点,从而可以避免间隙非线性引起的数值不稳定现象。仿真结果显示,俯仰立方非线性可以使间隙非线性系统产生振幅稳定的极限环振荡;控制面间隙对开环响应的影响随着来流速度的升高而减弱;在速度较高的情况下,控制面间隙导致闭环系统响应产生过阻尼现象。     

迟滞非线性二元机翼颤振特性分析

 采用多项式迟滞非线性模型建立二元机翼气动弹性运动方程，并用数值积分法进行求解。通过系统响应振幅随来流速度变化的分叉图和频谱分析发现，俯仰方向由于含有非线性因素，振动中的高阶分量随速度提高不断增加，并引起高次分叉。重点研究“机翼/空气”质量比以及“沉浮/俯仰”两个自由度的自然振动频率比对非线性颤振速度边界的影响，并提出可以通过提高自然振动频率比来减小迟滞非线性因素的不利影响。     

非线性气动弹性模型参考自适应控制

 首先将含有前/后缘双控制面二元机翼的动态方程以状态空间形式描述,然后考虑俯仰方向的迟滞非线性模型存在参数不确定性的情况下,利用Lyapunov稳定性理论进行了结构化模型参考自适应控制律设计。仿真结果显示：所设计的控制律能够使开环不稳定的气动弹性系统快速地达到稳定状态。由于最大控制面偏转的存在,来流速度较高时闭环系统仍会发生颤振;根据控制面最大偏转的不同取值,文中给出了闭环临界颤振速度的变化曲线。     

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

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

全动翼面间隙非线性气动弹性响应分析

战斗机全动翼面操纵系统中普遍存在间隙，其非线性特性不但影响系统地面振动试验数据的有效性，而且可能导致飞行器气动弹性稳定性边界产生偏差，因此有必要开展考虑间隙非线性的气动弹性响应分析研究。以全动翼面旋转方向间隙为研究对象，基于虚拟质量法线化模态振型，建立一组可以表达整个响应域变形的统一模态振型；采用有理函数拟合，将频域非定常气动力转换到时域；研究分析不同间隙参数对应的极限环响应特性。结果表明：当全动翼面结构系统中存在间隙时，在低于线性颤振边界的特定飞行区域，翼面会出现非线性极限环振荡现象，间隙参数会影响极限环振荡幅值、频率以及限环振荡的进入临界速压和发散临界速压。     

不同迎角的翼型气弹特性风洞实验研究

基于可在不同迎角下作沉浮、俯仰两自由度运动的翼段振动装置,在低速风洞中分别针对普通薄翼型NA-CA0012和风力机翼型NREL S809进行气动弹性测试,得到不同实验状态的气动弹性振动时域响应。分别观察到经典颤振和失速颤振现象,并证明了迎角改变对两种翼型颤振特性的影响。     

基于CFD/CSD耦合的跨声速非线性气动弹性研究

跨声速气动弹性数值数值模拟中,流场非线性特征明显,传统线化理论不再适用,需要采用时域方法,计算任意时刻的流场信息,进而得到正确的气动弹性现象。为AVICFD-X软件开发了气动弹性模块,该模块基于CFD/CSD(计算流体力学/计算结构力学)耦合模拟方法,包含流体、固体求解器以及交界面信息传输程序,在模块中分别求解流体控制方程和固体结果方程,通过界面进行流、固体之间的气动载荷以及运动状态等信息交换。AVICFD-X软件气动弹性模块数值模拟强迫振荡、静气弹、颤振等问题,验证气动弹性模块是否正确。数值模拟结果满足预期效果,与实验结果相吻合,充分说明了AVICFD-X软件气动弹性模块适用于跨声速非线性气动弹性研究。     

带结构刚度非线性的全动操纵面跨音速颤振特性研究

主要的研究工作是建立了以三维非定常N—S方程为主管方程的计算非定常气动力的高效计算方法，采用颤振方程和N—S方程耦合求解，时间推进的方法计算结构响应特性。发展了一种对既带有结构非线性又含有气动非线性的多重非线性的复杂颤振系统响应特性的分析方法。对全动操纵面转轴带有硬弹簧和间隙型非线性刚度的颤振系统进行了分析，计算结果表明这类系统具有复杂的极限环特性。     

非线性二元机翼气动弹性近似解析研究

 建立了不可压流动中多项式迟滞非线性二元机翼的气动弹性运动方程，然后利用谐波平衡法进行了求解。与数值积分结果比较分析表明，在系统发生二次分叉以前，谐波平衡法可以准确地预测极限环振荡的频率和振幅,通过频谱分析与时间响应历程讨论了谐波平衡法产生误差的原因。另外还研究了弹性轴位置对颤振特性的影响，随着弹性轴不断靠近翼弦中点，俯仰振幅不断增大，而沉浮振幅则存在一个极小值点。     

某型飞机操纵面间隙非线性颤振时域分析

操纵面间隙作为一种常见的结构非线性,是由飞机设计、制造、装配等众多环节产生,有可能引起极限环振荡（LCO）。极限环振荡通常表现为等幅振动,如果其振幅过大,也会影响机体结构完整性,引发结构失效。本文对操纵面间隙非线性颤振的时域分析方法进行研究,采用时域分析方法计算某型飞机的非线性颤振响应,并与频域描述函数法计算结果进行对比分析。结果表明：极限环振荡的临界速度和频率基本一致,时域分析方法能够准确计算全机操纵面间隙的非线性颤振临界速度,用来预测操纵面的极限环振荡是可行的,可以将其作为民用飞机适航符合性验证的理论分析方法之一。     

带有气动及结构非线性的二元机翼颤振分析

研究了翼型在低马赫数条件下的非定常气动特性,从翼型表面气流运动的角度对Leishman-Beddoes(L-B)模型进行了修正,并在此基础上建立了适合低马赫数颤振研究且带有气动及结构非线性的二元机翼气弹系统分析模型.对比低马赫数翼型气动载荷试验结果表明对L-B模型的修正是有效的,且机翼颤振试验结果亦验证了二元机翼气弹分析模型.研究结果表明:二元机翼气弹系统的失速颤振与初始变距角和来流速度密切相关,且耦合的三次非线性变距和浮沉刚度是造成系统呈现准周期运动的主要原因.      

Nonlinear aeroelastic analysis of the folding fin with freeplay under thermal environment

The nonlinear aeroelastic behavior of a folding fin in supersonic flow is investigated in this paper. The finite element model of the fin is established and the deployable hinges are represented by three torsion springs with the freeplay nonlinearity. The aerodynamic grid point is assumed to be at the center of each aerodynamic box for simplicity. The aerodynamic governing equation is given by using the infinite plate spline method and the modified linear piston theory. An improved fixed-interface modal synthesis method, which can reduce the rigid connections at the interface, is developed to save the problem size and computation time. The uniform temperature field is applied to create the thermal environment. For the linear flutter analyses, the flutter speed increases first and then decreases with the rise of the hinge stiffness due to the change of the flutter coupling mechanism. For the nonlinear analyses, a larger freeplay angle results in a higher vibration divergent speed. Two different types of limit cycle oscillations and a multiperiodic motion are observed in the wide range of airspeed under the linear flutter boundary. The linear flutter speed shows a slight descend in the thermal environment, but the effect of the temperature on the vibration divergent speed is different under different hinge stiffnesses when there exists freeplay.     

Nonlinear aeroelastic analysis of airfoils: bifurcation and chaos

Different types of structural and aerodynamic nonlinearities commonly encountered in aeronautical engineering are discussed. The equations of motion of a two-dimensional airfoil oscillating in pitch and plunge are derived for a structural nonlinearity using subsonic aerodynamics theory. Three classical nonlinearities, namely, cubic, freeplay and hysteresis are investigated in some detail. The governing equations are reduced to a set of ordinary differential equations suitable for numerical simulations and analytical investigation of the system stability. The onset of Hopf-bifurcation, and amplitudes and frequencies of limit cycle oscillations are investigated, with examples given for a cubic hardening spring. For various geometries of the freeplay, bifurcations and chaos are discussed via the phase plane, Poincaré maps, and Lyapunov spectrum. The route to chaos is investigated from bifurcation diagrams, and for the freeplay nonlinearity it is shown that frequency doubling is the most commonly observed route. Examples of aerodynamic nonlinearities arising from transonic flow and dynamic stall are discussed, and special attention is paid to numerical simulation results for dynamic stall using a time-synthesized method for the unsteady aerodynamics. The assumption of uniform flow is usually not met in practice since perturbations in velocities are encountered in flight. Longitudinal atmospheric turbulence is introduced to show its effect on both the flutter boundary and the onset of Hopf-bifurcation for a cubic restoring force.     

间隙非线性气动弹性系统的辨识

间隙非线性气动弹性系统有非常丰富的动力学特征。由于工程实际中间隙不易准确测量甚至无法测量,通常需要通过辨识来建立其数学模型,其中间隙开关点位置是辨识问题中的难点。本文利用间隙开关点将间隙非线性系统区分为3个线性子系统,用Hammerstein模型表示其非线性部分,构造了开关点的迭代序列;并采用非迭代和迭代相结合的方法进行求解,从而成功获得包含间隙开关点在内的该系统所有模型参数;以间隙非线性的二元翼段为例,验证了该辨识方法的有效性。     

Nonlinear flutter wind tunnel test and numerical analysis of folding fins with freeplay nonlinearities

The flutter characteristics of folding control fins with freeplay are investigated by numerical simulation and flutter wind tunnel tests. Based on the characteristics of the structures, fins with different freeplay angles are designed. For a 0° angle of attack, wind tunnel tests of these fins are conducted, and vibration is observed by accelerometers and a high-speed camera. By the expansion of the connected relationships, the governing equations of fit for the nonlinear aeroelastic analysis are established by the free-interface component mode synthesis method. Based on the results of the wind tunnel tests, the flutter characteristics of fins with different freeplay angles are analyzed. The results show that the vibration divergent speed is increased, and the divergent speed is higher than the flutter speed of the nominal linear system. The vibration divergent speed is increased along with an increase in the freeplay angle. The developed free-interface component mode synthesis method could be used to establish governing equations and to analyze the characteristics of nonlinear aeroelastic systems. The results of the numerical simulations and the wind tunnel tests indicate the same trends and critical velocities.     

Structural Nonlinear Flutter Characteristics Analysis for an Actuator-fin System with Dynamic Stiffness

The flutter characteristics of an actuator-fin system are investigated with structural nonlinearity and dynamic stiffness of the electric motor. The component mode substitution method is used to establish the nonlinear governing equations in time domain and frequency domain based on the fundamental dynamic equations of the electric motor and decelerator. The existing describing function method and a proposed iterative method are used to obtain the flutter characteristics containing preload freeplay nonlinea...     

Numerical studies of a non-linear aeroelastic system with plunging and pitching freeplays in supersonic/hypersonic regimes

The flutter and post flutter of a two-dimensional double-wedge lifting surface with combined freeplay and cubic stiffness nonlinearities in both plunging and pitching degrees-of-freedom operating in supersonic/hypersonic flight speed regimes have been analyzed. In addition to the structural nonlinearities, the third-order piston theory aerodynamics is used to evaluate the unsteady non-linear aerodynamic force and moment. Such model accounts for stiffness and damping contributions produced by the aerodynamic loads. Responses involving limit cycle oscillation and chaotic motion are observed over a large number of parameters that characterizes the aeroelastic system. Results of the present study show that the freeplay in the pitching degree-of-freedom and soft/hard cubic stiffness in the pitching and plunging degrees-of-freedom have significant effects on the LCOs exhibited by the aeroelastic system in the supersonic/hypersonic flight speed regimes. The simulations also show that the aeroelastic system behavior is greatly affected by physical structural parameters, such as the radius of gyration and the frequency ratio, especially in post-flutter regimes, when accounting for all system nonlinearities. It has been shown that at high Mach numbers the non-linear aerodynamic stiffness yields detrimental effects from the aeroelastic point of view, while the damping one do not.     

Flutter analysis of an airfoil with nonlinear damping using equivalent linearization

The equivalent linearization method （ELM） is modified to investigate the nonlinear flut- ter system of an airfoil with a cubic damping. After obtaining the linearization quantity of the cubic nonlinearity by the ELM, an equivalent system can be deduced and then investigated by linear flut- ter analysis methods. Different from the routine procedures of the ELM, the frequency rather than the amplitude of limit cycle oscillation （LCO） is chosen as an active increment to produce bifurca- tion charts. Numerical examples show that this modification makes the ELM much more efficient. Meanwhile, the LCOs obtained by the ELM are in good agreement with numerical solutions. The nonlinear damping can delay the occurrence of secondary bifurcation. On the other hand, it has marginal influence on bifurcation characteristics or LCOs.     

具有不对称间隙的二元机翼自激振动的混沌分析

利用数值方法分析了无预荷载情况下具有不对称间隙的气动弹性非线性与结构非线性耦合的二元机翼颤振问题的动力学响应.通过研究速度比U*/U*L从零逐渐增大到1时系统的全局演化规律,发现在无预荷载的不对称分段线性系统中存在混沌运动形式,并对其进行具体分析,得到产生混沌的原因,为进一步控制和减小颤振提供理论帮助,对飞机的飞行结构设计及飞机结构的日常维护具有指导意义.