Bifurcation of the Curved Panel in Supersonic Air Flow

采用Galerkin方法建立了超音速气流中二维曲壁板的非线性热气动弹性运动方程。用von Karman大变形理论来考虑曲壁板的大变形。用准定常的一阶活塞理论模拟曲壁板上表面受到的气动力。在不同来流速压和温升条件下,基于分岔理论研究了具有不同初始几何曲率的曲壁板系统对应的定常状态方程（组）的解的个数、性态和动态稳定性,并对方程（组）进行了解曲线的跟踪分析。研究表明,不同条件下,方程组的解特性不同,并且随着初始几何曲率和温升条件的变化,系统的失稳机理发生变化。超音速气流中的二维曲壁板系统存在动态Hopf分岔和静态鞍－结点分岔两种失稳现象,但不会发生热屈曲失稳。     

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

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

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

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

机电耦联系统余维3动态分岔研究

以r_{sl}, r_{f}以及x_{c}为分岔参数，对具有串补电容的单 机无穷大电力系统的失稳振荡问题，运用动态分岔理论进行了研究. 对系统同时出现有3对 纯虚根特征值的一类多参数高余维分岔情况，运用中心流行方法降维后得到约化方程，对此 强非线性约化方程的求解难点，运用多参数稳定性理论、谐波平衡法、归一化技术和Normal Form方法，得到了系统的解析解. 由分析得知，系统会出现3种Hopf分岔情况、二维环面 情况，以及三维环面分岔解，甚至会出现四维环面，或者更高维的环面分岔. 详细讨论 了系统各种分岔解的稳定性条件和稳定区域，并作了详细的数值分析加以验证.     

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

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

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

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

受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用.将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度,是俘能器优化设计的重要措施之一.但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究.本文以受限亚音速气流中倒置的二维悬臂壁板为对象,以理论分析及风洞实验为手段,研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律.在理论分析中,首先应用镜像函数法来处理壁面约束条件,基于算子理论研究获得了以Possio积分方程为表征的壁板气动力,壁面效应实际表征为一包含移位Tricomi算子的复合算子;然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题;最后,依据Wererstrass定理并利用最小二乘法求解该最优函数,以获得系统的失稳临界参数.在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验.理论分析结果表明,壁板会发生发散(静气动弹性)失稳,临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况);通过理论与风洞实验结果的对比分析,验证了本文气动力及理论分析的适用性及准确性.针对倒置悬臂壁板结构的气动弹性失稳问题,本文提出的方法不涉及系统方程的离散及特征值求解问题,而是将其转化为了定区间上的函数逼近问题进行求解,这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.     

受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用. 将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度, 是俘能器优化设计的重要措施之一. 但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究. 本文以受限亚音速气流中倒置的二维悬臂壁板为对象, 以理论分析及风洞实验为手段, 研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律. 在理论分析中, 首先应用镜像函数法来处理壁面约束条件, 基于算子理论研究获得了以Possio积分方程为表征的壁板气动力, 壁面效应实际表征为一包含移位Tricomi算子的复合算子; 然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题; 最后, 依据Wererstrass定理并利用最小二乘法求解该最优函数, 以获得系统的失稳临界参数. 在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验. 理论分析结果表明, 壁板会发生发散(静气动弹性)失稳, 临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况); 通过理论与风洞实验结果的对比分析, 验证了本文气动力及理论分析的适用性及准确性. 针对倒置悬臂壁板结构的气动弹性失稳问题, 本文提出的方法不涉及系统方程的离散及特征值求解问题, 而是将其转化为了定区间上的函数逼近问题进行求解, 这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.      

热环境中功能梯度圆柱薄壳分岔屈曲的边界约束效应

基于前屈曲一致理论,研究了热环境中受轴压功能梯度圆柱薄壳分岔屈曲的边界约束效应问题.导出前屈曲变形的解析解,结合分离变量法与有限差分法求解分岔屈曲控制方程,由此导出确定临界轴压的非线性特征值问题.考虑材料热物性质与温度的相关性,分别就两端简支和两端固支边界条件,分析了温度梯度、初始几何缺陷、组分材料体积分数对分岔屈曲临...     

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

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

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.     

Nonlinear dynamics analysis of a two-dimensional thin panel with an external forcing in incompressible subsonic flow

Based on the potential theory of incompressible flow and the energy method, a two-dimensional simply supported thin panel subjected to external forcing and uniform incompressible subsonic flow is theoretically modeled. The nonlinear cubic stiffness and viscous damper in the middle of the panel is considered. Transformation of the governing partial differential equation to a set of ordinary differential equations is performed through the Galerkin method. The stability of the fixed points of the panel system is analyzed. The regions of different motion types of the panel system are investigated in different parameter spaces. The rich dynamic behaviors are presented as bifurcation diagrams, phase-plane portraits, Poincaré maps and maximum Lyapunov exponents based on carefully numerical simulations.     

Nonlinear dynamics of extensible fluid-conveying pipes,supported at both ends

In this paper, the post-divergence behaviour of extensible fluid-conveying pipes supported at both ends is studied using the weakly nonlinear equations of motion of Semler, Li and Païdoussis. The two coupled nonlinear partial differential equations are discretized via Galerkin's method and the resulting set of ordinary differential equations is solved either by Houbolt's finite difference method or via AUTO. Typically, the pipe is stable at its original static equilibrium position up to the flow velocity where it loses stability by static divergence via a supercritical pitchfork bifurcation. The amplitude of the resultant buckling increases with increasing flow, but no secondary instability occurs beyond the pitchfork bifurcation. The effects of the system parameters on pipe behaviour as well as the possibility of a subcritical pitchfork bifurcation have also been studied.     

A theoretical computational model of a plate in hypersonic flow

A theoretical model of an elastic panel in hypersonic flow is derived to be used for design and analysis. The nonlinear von Kármán plate equations are coupled with 1st order Piston Theory and linearized at the nonlinear steady-state deformation due to static pressure differential and thermal loads. Eigenvalue analysis is applied to determine the system’s stability, natural frequencies and mode shapes. Numerically time marching the equations provides transient response prediction which can be used to estimate limit cycle oscillation amplitude, frequency and time to onset. The model’s predictive capability is assessed by comparison to an experiment conducted at a free stream flow of Mach 6. Good agreement is shown between the theoretical and experimental natural frequencies and mode shapes of the fluid–structure system. Stability analysis is performed using linear and nonlinear methods to plot stability, flutter and buckling zones on a free stream static pressure vs temperature differential plane.     

Chaotic Analysis of Nonlinear Viscoelastic Panel Flutter in Supersonic Flow

In this paper chaotic behavior of nonlinear viscoelastic panels in asupersonic flow is investigated. The governing equations, based on vonKàarmàn's large deflection theory of isotropic flat plates, areconsidered with viscoelastic structural damping of Kelvin's modelincluded. Quasi-steady aerodynamic panel loadings are determined usingpiston theory. The effect of constant axial loading in the panel middlesurface and static pressure differential have also been included in thegoverning equation. The panel nonlinear partial differential equation istransformed into a set of nonlinear ordinary differential equationsthrough a Galerkin approach. The resulting system of equations is solvedthrough the fourth and fifth-order Runge–Kutta–Fehlberg (RKF-45)integration method. Static (divergence) and Hopf (flutter) bifurcationboundaries are presented for various levels of viscoelastic structuraldamping. Despite the deterministic nature of the system of equations,the dynamic panel response can become random-like. Chaotic analysis isperformed using several conventional criteria. Results are indicative ofthe important influence of structural damping on the domain of chaoticregion.     

Nonlinear dynamics of doubly curved shallow microshells

The nonlinear dynamical characteristics of a doubly curved shallow microshell are investigated thoroughly. A consistent nonlinear model for the microshell is developed on the basis of the modified couple stress theory (MCST) in an orthogonal curvilinear coordinate system. In particular, based on Donnell’s nonlinear theory, the expressions for the strain and the symmetric rotation gradient tensors are obtained in the framework of MCST, which are then used to derive the potential energy of the microshell. The analytical geometrically nonlinear equations of motion of the doubly microshell are obtained for in-plane displacements as well as the out-of-plane one. These equations of partial differential type are reduced to a large set of ordinary differential equations making use of a two-dimensional Galerkin scheme. Extensive numerical simulations are conducted to obtain the nonlinear resonant response of the system for various principal radii of curvature and to examine the effect of modal interactions and the length-scale parameter.     

Nonlinear vibration of slightly curved pipe with conveying pulsating fluid

The nonlinear governing motion equation of slightly curved pipe with conveying pulsating fluid is set up by Hamilton’s principle. The motion equation is discretized into a set of low dimensional system of nonlinear ordinary differential equations by the Galerkin method. Linear analysis of system is performed upon this set of equations. The effect of amplitude of initial deflection and flow velocity on linear dynamic of system is analyzed. Curves of the resonance responses about \(\varOmega \approx {\omega _\mathrm{{1}}}\) and \(\varOmega \approx \mathrm{{2}}{\omega _\mathrm{{1}}}\) are performed by means of the pseudo-arclength continuation technique. The global nonlinear dynamic of system is analyzed by establishing the bifurcation diagrams. The dynamical behaviors are identified by the phase diagram and Poincare maps. The periodic motion, chaotic motion and quasi-periodic motion are found in this system.