超音速气流中受热曲壁板的非线性颤振特性

基于von Karman 大变形理论及带有曲率修正的一阶活塞理论, 用Galerkin方法建立了超音速气流中受热二维曲壁板的非线性气动弹性运动方程; 采用牛顿迭代法计算得到由静气动载荷和热载荷引起的静气动弹性变形; 根据李雅谱诺夫间接法分析了壁板初始曲率与温升对颤振边界的影响; 对二维曲壁板的非线性气动弹性方程组进行数值积分求解,分析了动压参数对受热二维曲壁板分岔特性的影响, 给出了典型状态下曲壁板非线性颤振响应的时程图与相图. 分析结果表明对小初始曲率的曲壁板, 温升对其静气动弹性变形影响较大, 且随着温升的增加其颤振临界动压急剧减小; 对具有较大初始曲率的曲壁板, 温升对其静气动弹性变形的影响较弱, 且随着温升的增加颤振临界动压基本保持不变. 初始几何曲率与气动热效应使得曲壁板具有复杂的动力学特性, 不再像平壁板一样, 经过倍周期分岔进入混沌, 而会出现由静变形状态直接进入混沌运动的现象, 且在混沌运动区域中还会出现静态稳定点或谐波运动, 在大曲率情况下, 曲壁板不会产生混沌运动, 而是幅值在一定范围内的极限带振荡.      

采用压电材料提高超声速飞行器壁板结构的颤振特性

采用压电材料对结构进行振动主动控制已经进行了广泛研究,论文进一步采用压电材料改进超声速壁板结构的气动弹性颤振特性,研究中考虑压电材料力电耦合效应的影响.采用Hamilton原理和Rayleigh-Ritz方法建立壁板及压电材料整体结构的运动方程,采用超声速活塞理论模拟气动力,利用加速度反馈控制策略对压电材料施加外电压,获得结构的主动质量.求解运动方程的特征值问题获得固有频率,进而确定气动弹性颤振边界,分析了反馈控制增益对超声速飞行器壁板结构主动颤振特性的影响,研究表明,采用压电材料可以提高超声速壁板结构的气动弹性颤振特性.     

激波主导流动下壁板的热气动弹性稳定性理论分析

针对激波主导流动下弹性壁板的热气动弹性稳定性分析问题,建立了基于当地活塞流理论的分析模型,并用数值仿真方法来验证其正确性. 首先基于Hamilton原理和Von-Karman大变形理论,建立壁板的热气动弹性运动方程,其中假设壁板受热后温度均匀分布,激波前后区域的气动力模型采用当地一阶活塞流理论;利用Galerkin方法将具有连续参数系统的偏微分颤振方程离散为有限个自由度的常微分方程;基于李雅普诺夫间接法将非线性颤振方程组在平衡位置处进行线化,再用Routh-Hurwits判据来判断线性系统的稳定性,从而来推论出非线性颤振系统的气动弹性稳定性. 在时域中采用龙格--库塔法对非线性颤振方程进行数值积分,得到壁板非线性颤振响应的时间历程,与理论分析结果进行对比. 研究结果表明,壁板受到斜激波冲击时,更容易发生颤振失稳,并且激波强度越大,极限环幅值和频率越大;激波主导流场中的壁板失稳边界不同于传统单纯超声速气流中壁板颤振的失稳边界;只有在斜激波前后不同的动压值都满足颤振稳定性边界的条件下,壁板才可能保持其气动弹性稳定性.      

激波主导流动下壁板的热气动弹性稳定性理论分析

针对激波主导流动下弹性壁板的热气动弹性稳定性分析问题,建立了基于当地活塞流理论的分析模型,并用数值仿真方法来验证其正确性.首先基于Hamilton原理和Von-Karman大变形理论,建立壁板的热气动弹性运动方程,其中假设壁板受热后温度均匀分布,激波前后区域的气动力模型采用当地一阶活塞流理论;利用Galerkin方法将具有连续参数系统的偏微分颤振方程离散为有限个自由度的常微分方程;基于李雅普诺夫间接法将非线性颤振方程组在平衡位置处进行线化,再用Routh-Hurwits判据来判断线性系统的稳定性,从而来推论出非线性颤振系统的气动弹性稳定性.在时域中采用龙格-库塔法对非线性颤振方程进行数值积分,得到壁板非线性颤振响应的时间历程,与理论分析结果进行对比.研究结果表明,壁板受到斜激波冲击时,更容易发生颤振失稳,并且激波强度越大,极限环幅值和频率越大;激波主导流场中的壁板失稳边界不同于传统单纯超声速气流中壁板颤振的失稳边界;只有在斜激波前后不同的动压值都满足颤振稳定性边界的条件下,壁板才可能保持其气动弹性稳定性.     

转子-机匣局部碰摩的分叉与混沌行为分析

研究了转子-机匣系统发生碰摩时的分叉与混沌行为,分析了转子机匣频率比与刚度比、偏心质量等参数对系统分叉与混沌特性的影响.当转子机匣系统发生碰摩时除了通过倍周期、阵发性和拟周期分叉进入混沌外,还发现了孪生叉形分叉现象,呈现出非常丰富的动力学行为.     

热环境下复合材料壁板的振动特性分析

针对飞行器中常见的壁板结构,运用能量原理和变分方法,建立了定常温度场下复合材料壁板振动的控制方程以及相应的有限元分析模型。分析了热环境对壁板振动特性影响的机理,同时提出了一种针对热环境下复合材料壁板振动特性分析的线性化计算方法。采用这种方法,热环境的影响以一个热刚度项和一个热载荷项的形式出现在常温下的振动运动方程中,由此可以较准确地模拟热效应对结构振动特性的影响。通过对热环境下复合材料壁板振动固有特性数值分析结果的对比,验证了本文方法的可行性和计算精度。同时分析结果表明,热效应产生的诱导应力对结构刚度的影响是导致壁板固有振动频率降低的主要原因。     

高速飞行器壁板颤振的分析模型和分析方法

壁板颤振是壁板结构在高速气流中产生的一种自激振动,在超声速和高超声速飞行器上特别容易发生这种现象。壁板颤振引发的非线性振动将对高速飞行器结构的疲劳强度、飞行性能和飞行安全带来不利的影响。随着高速飞行器设计中各项研究工作的开展,壁板颤振问题受到了到越来越多的重视。本文阐述了目前国内外学者在高速飞行器壁板颤振分析领域的研究现状及壁板颤振研究中常用的六种分析模型,并根据壁板颤振分析中使用的结构理论和气动力理论,详述了这种分类的依据。文中还介绍了温度、气流偏角、壁板几何尺寸及边界条件对壁板颤振的影响规律和目前常用于分析壁板颤振问题的频域和时域方法,总结了各种分析方法的优缺点。最后归纳了目前在高速飞行器壁板颤振研究中得出的几个重要结论,提出了今后在高速飞行器壁板颤振研究中需要解决的若干问题。     

外激励作用下亚音速二维壁板的复杂响应研究

研究了亚音速流中二维壁板在外激励作用下的复杂响应问题。采用迦辽金方法将非线性运动控制方程离散为常微分方程组,采用数值方法进行计算,研究了壁板系统的复杂响应。应用最大李亚普诺夫指数和庞加莱截面方法对系统的运动性质进行了判定。结果表明,系统随着参数的变化呈现出复杂的响应,系统的周期运动与混沌运动会相间出现;系统由周期运动进...     

一类滞回系统混沌振动的数值分析

在不同参数下,Brouc-Wen滞回模型使系统具有软或硬式响应特性,导致系统有非线性振动特性。利用数值方法,本文给出单自由度滞回系统 稳态振动最大振幅与频率之间的关系曲线。分析了滞回参数对硬式响应特性滞回振动系统的分叉与混沌的影响,发现一些新的现象。     

资金项目:国家自然科学基金资助项目

采用能计及非线性结构刚度的颤振方程为控制方程,和非定常N-S方程耦合求解,运用龙格-库塔方法在时域内求解结构响应的时间历程,从而确定颤振临界条件.计算了带结构刚度非线性的跨音速颤振特性.计算结果表明,结构刚度非线性对颤振特性有明显的影响.由于同时具有结构和气动力非线性,导致了具有复杂振荡极限环的特性.     

热过屈曲功能梯度壁板的气动弹性颤振

功能梯度材料的宏观物理性能随空间位置连续变化,能充分减少不同组份材料结合部位界面性能的不匹配因素.功能梯度壁板用作高速飞行器的热防护结构,能有效消除气动加热带来的壁板内部热应力集中.本文考虑热过屈曲变形引入的结构几何非线性,分析功能梯度壁板的气动弹性颤振边界.基于幂函数材料分布假设,采用混合定律计算功能梯度材料的等效力学性能.根据一阶剪切变形板理论、冯·卡门应变-位移关系和一阶活塞理论,基于虚功原理建立超声速气流中受热功能梯度壁板的非线性气动弹性有限元方程.采用牛顿-拉弗森迭代法数值求解壁板的热屈曲变形,分析超声速气流对热屈曲变形的影响机理.在壁板热过屈曲的静力平衡位置分析动态稳定性,确定了壁板的颤振边界.研究表明,当陶瓷-金属功能梯度壁板的组份材料沿厚度方向梯度分布时,会破坏结构的对称性导致壁板在面内热应力作用下发生指向金属侧的热屈曲变形.超声速气流中壁板热屈曲变形最大的位置随气流速压增大向下游推移,并伴随屈曲变形量的减小.热过屈曲壁板的几何非线性效应会提高壁板的颤振边界,这种影响在高温、低无量纲速压且壁板发生大挠度热屈曲变形时表现显著.较高无量纲气流速压下由于壁板的热屈曲变形被气动力限定在小挠度范围,几何非线性效应不明显.      

On the aeroelastic stability and bifurcation structure of subsonic nonlinear thin panels subjected to external excitation

Dynamic behavior of panels exposed to subsonic flow subjected to external excitation is investigated in this paper. The von Karman’s large deflection equations of motion for a flexible panel and Kelvin’s model of structural damping is considered to derive the governing equation. The panel under study is two-dimensional and simply supported. A Galerkin-type solution is introduced to derive the unsteady aerodynamic pressure from the linearized potential equation of uniform incompressible flow. The governing partial differential equation is transformed to a series of ordinary differential equations by using Galerkin method. The aeroelastic stability of the linear panel system is presented in a qualitative analysis and numerical study. The fourth-order Runge-Kutta numerical algorithm is used to conduct the numerical simulations to investigate the bifurcation structure of the nonlinear panel system and the distributions of chaotic regions are shown in the different parameter spaces. The results shows that the panel loses its stability by divergence not flutter in subsonic flow; the number of the fixed points and their stabilities change after the dynamic pressure exceeds the critical value; the chaotic regions and periodic regions appear alternately in parameter spaces; the single period motion trajectories change rhythmically in different periodic regions; the route from periodic motion to chaos is via doubling-period bifurcation.     

To a formulation of the flutter problem for a conical shell with a small apex angle

The well-known piston theory formula for the excess aerodynamic pressure is used in the majority of works devoted to the panel flutter of shells. In this paper a refined expression for the excess pressure is proposed to take into account the irregularity of undisturbed flow parameters. The case of moderate supersonic velocities is studied in detail. The critical velocity problem is reduced to a new eigenproblem in the panel flutter theory.     

Nonlinear dynamics of an aeroelastic airfoil with free-play in transonic flow

Nonlinear dynamic behaviors of an aeroelastic airfoil with free-play in transonic air flow are studied. The aeroelastic response is obtained by using time-marching approach with computational fluid dynamics (CFD) and reduced order model (ROM) techniques. Several standardized tests of transonic flutter are presented to validate numerical approaches. It is found that in time-marching approach with CFD technique, the time-step size has a significant effect on the calculated aeroelastic response, especially for cases considering both structural and aerodynamic nonlinearities. The nonlinear dynamic behavior for the present model in transonic air flow is greatly different from that in subsonic regime where only simple harmonic oscillations are observed. Major features of the responses in transonic air flow at different flow speeds can be summarized as follows. The aeroelastic responses with the amplitude near the free-play are dominated by single degree of freedom flutter mechanism, and snap-though phenomenon can be observed when the air speed is low. The bifurcation diagram can be captured by using ROM technique, and it is observed that the route to chaos for the present model is via period-doubling, which is essentially caused by the free-play nonlinearity. When the flow speed approaches the linear flutter speed, the aeroelastic system vibrates with large amplitude, which is dominated by the aerodynamic nonlinearity. Effects of boundary layer and airfoil profile on the nonlinear responses of the aeroelastic system are also discussed.     

THE DYNAMIC BEHAVIORS OF VISCOELASTIC PIPE CONVEYING FLUID WITH THE KELVIN MODEL

Based on the differential constitutive relationship of linear viscoelastic, material, a solid-liquid coupling vibration equation for viscoelastic pipe conveying fluid is derived by the D'Alembert's principle. The critical flow velocities and natural frequencies of the cantilever pipe conveying fluid with the Kelvin model (flutter instability) are calculated with the modified finite difference method in the form of the recurrence formula. The curves between the complex frequencies of the first, second and third mode and flow velocity of the pipe are plotted. On the basis of the numerical, calculation results, the dynamic behaviors and stability of the pipe are discussed. It should be pointed out that the delay time of viscoelastic material with the Kelvin model has a remarkable effect on the dynamic characteristics and stability behaviors of the cantilevered pipe conveying fluid, which is a gyroscopic non-conservative system.     

Influence of boundary conditions relaxation on panel flutter with compressive in-plane loads

The influence of boundary conditions relaxation on two-dimensional panel flutter is studied in the presence of in-plane loading. The boundary value problem of the panel involves time-dependent boundary conditions that are converted into autonomous form using a special coordinate transformation. Galerkin's method is used to discretize the panel partial differential equation of motion into six nonlinear ordinary differential equations. The influence of boundary conditions relaxation on the panel modal frequencies and LCO amplitudes in the time and frequency domains is examined using the windowed short time Fourier transform and wavelet transform. The relaxation and system nonlinearity are found to have opposite effects on the time evolution of the panel frequency. Depending on the system damping and dynamic pressure, the panel frequency can increase or decrease with time as the boundary conditions approach the state of simple supports. Bifurcation diagrams are generated by taking the relaxation parameter, dynamic pressure, and in-plane load as control parameters. The corresponding largest Lyapunov exponent is also determined. They reveal complex dynamic characteristics of the panel, including regions of periodic, quasi-periodic, and chaotic motions.     

超音速气流中受热壁板的稳定性分析

采用Galerkin方法建立二维壁板的非线性气动弹性运动方程，用一阶活塞理论模拟壁板 受到的气动力. 基于李雅普诺夫间接法分析了平壁板的稳定性，得到了壁板失稳的边界 曲线；采用牛顿迭代法分析了壁板的屈曲变形，进而分析了后屈曲状态下壁板的稳定性； 在时域中分析了后屈曲状态下壁板的颤振边界. 分析结果表明，为了保证计算精度， 在二维壁板的静态失稳及过屈曲变形分析中，至少要取二阶谐波模态；在平壁板的超音速颤 振(动态失稳)边界分析中至少应取四阶模态. 还对壁板的温升，壁板长厚比、壁板密 度和气流马赫数作了无量纲变参分析，研究了这些参数的变化对壁板稳定性的影响规律. 研 究中发现，当气流速压较低时壁板一般会稳定在低阶谐波模态的屈曲变形位置，但是如果系 统出现多个渐近稳定的不动点，即使作用在壁板上的气流速压很低，壁板也有可能在较低速 压下发生二次失稳型颤振.     

Viscous and inviscid interactions of an oblique shock with a flexible panel

The complex self-sustained oscillations arising from the interaction of an oblique shock with a flexible panel in both the inviscid and viscous regimes have been investigated numerically. The aeroelastic interactions are simulated using either the Euler or the full compressible Navier–Stokes equations coupled to the nonlinear von Karman plate equations. Results demonstrate that for a sufficiently strong shock limit-cycle oscillations emerge from either subcritical or supercritical bifurcations even in the absence of viscous separated flow effects. The critical dynamic pressure diminishes with increasing shock strength and can be much lower than that corresponding to standard panel flutter. Significant changes in panel dynamics were also found as a function of the shock impingement point and cavity pressure. For viscous laminar flow above the panel without a shock, high-frequency periodic oscillations appear due to the coupling of boundary-layer instabilities with high-mode flexural deflections. For a separated shock laminar boundary layer interaction, non-periodic self-excited oscillations arise which can result in a significant reduction in the extent of the time-averaged separation region. This finding suggests the potential use of an aeroelastically tailored flexible panel as a means of passive flow control. Forced panel oscillations, induced by a specified variable cavity pressure underneath the panel, were also found to be effective in reducing separation. For both inviscid and viscous interactions, the significant unsteadiness generated by the fluttering panel propagates along the complex reflected expansion/recompression wave system.