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

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

Dynamic analysis for axially moving viscoelastic panels

In this study, stability and dynamic behaviour of axially moving viscoelastic panels are investigated with the help of the classical modal analysis. We use the flat panel theory combined with the Kelvin–Voigt viscoelastic constitutive model, and we include the material derivative in the viscoelastic relations. Complex eigenvalues for the moving viscoelastic panel are studied with respect to the panel velocity, and the corresponding eigenfunctions are found using central finite differences. The governing equation for the transverse displacement of the panel is of fifth order in space, and thus five boundary conditions are set for the problem. The fifth condition is derived and set at the in-flow end for clamped–clamped and clamped-simply supported panels. The numerical results suggest that the moving viscoelastic panel undergoes divergence instability for low values of viscosity. They also show that the critical panel velocity increases when viscosity is increased and that the viscoelastic panel does not experience instability with a sufficiently high viscosity coefficient. For the cases with low viscosity, the modes and velocities corresponding to divergence instability are found numerically. We also report that the value of bending rigidity (bending stiffness) affects the distance between the divergence velocity and the flutter velocity: the higher the bending rigidity, the larger the distance.     

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.     

AEROELASTIC OPTIMIZATION DESIGN OF COMPOSITE WING WITH CONSIDERATION OF STRUCTURAL INSTABILITY

从工程数学求解和有限元分析角度对复合材料结构的稳定性分析方法进行研究,基于这两个方面分别建立了同时考虑壁板稳定性约束和气动弹性约 束的气动弹性优化技术,并以大展弦比复合材料机翼为对象,进行气动弹性综合优化设计。研究表明,机翼气动弹性优化中若不考虑稳定性约束条件,虽然可以获得较小结构重量,但往往不满足稳定性要求;相比从有限元角度考虑结构失稳特征的气动弹性综合优化设计方法,通过工程数学方法对机翼结构进行分区失稳分析优化可以更加精准地控制变量,在满足各项性能指标,特别是稳定性约束的同时,进一步减轻了结构重量,提高了结构失稳因子。     

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

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

Aeroelastic stability of a cantilevered plate in yawed subsonic flow

The aeroelastic stability of cantilevered plates with their clamped edge oriented both parallel and normal to subsonic flow is a classical fluid–structure interaction problem. When the clamped edge is parallel to the flow the system loses stability in a coupled bending and torsion motion known as wing flutter. When the clamped edge is normal to the flow the instability is exclusively bending and is referred to as flapping flag flutter. This paper explores the stability of plates during the transition between these classic aeroelastic configurations. The aeroelastic model couples a classical beam structural model to a three-dimensional vortex lattice aerodynamic model. The aeroelastic stability is evaluated in the frequency domain and the flutter boundary is presented as the plate is rotated from the flapping flag to the wing configuration. The transition between the flag-like and wing-like instability is often abrupt and the yaw angle of the flow for the transition is dependent on the relative spacing of the first torsion and second bending natural frequencies. This paper also includes ground vibration and aeroelastic experiments carried out in the Duke University Wind Tunnel that confirm the theoretical predictions.     

Instability of two streaming conducting and dielectric bounded fluids in porous medium under time-varying electric field

In this paper, we consider the instability of the interface between two superposed streaming conducting and dielectric fluids of finite depths through porous medium in a vertical electric field varying periodically with time. A damped Mathieu equation with complex coefficients is obtained. The method of multiple scales is used to obtain an approximate solution of this equation, and then to analyze the stability criteria of the system. We distinguish between the non-resonance case, and the resonance case, respectively. It is found, in the first case, that both the porosity of porous medium, and the kinematic viscosities have stabilizing effects, and the medium permeability has a destabilizing effect on the system. While in the second case, it is found that each of the frequency of the electric field, and the fluid velocities, as well as the medium permeability, has a stabilizing effect, and decreases the value of the resonance point, while each of the porosity of the porous medium, and the kinematic viscosities has a destabilizing effect, and increases the value of the resonance point. In the absence of both streaming velocities and porous medium, we obtain the canonical form of the Mathieu equation. It is found that the fluid depth and the surface tension have a destabilizing effect on the system. This instability sets in for any value of the fluid depth, and by increasing the depth, the instability holds for higher values of the electric potential; while the surface tension has no effect on the instability region for small wavenumber values. Finally, the case of a steady electric field in the presence of a porous medium is also investigated, and the stability conditions show that each of the fluid depths and the porosity of the porous medium ɛ has a destabilizing effect, while the fluid velocities have stabilizing effect. The stability conditions for two limiting cases of interest, the case of purely fluids), and the case of absence of streaming, are also obtained and discussed in detail.     

Stationary Fingering Instability in a Non-equilibrium Porous Medium with Coupled Molecular Diffusion

Finger type double diffusive convective instability in a fluid-saturated porous medium is studied in the presence of coupled heat-solute diffusion. A local thermal non-equilibrium (LTNE) condition is invoked to model the Darcian porous medium which takes into account the energy transfer between the fluid and solid phases. Linear stability theory is implemented to compute the critical thermal Rayleigh number and the corresponding wavenumber exactly for the onset of stationary convection. The effects of Soret and Dufour cross-diffusion parameters, inter-phase heat transfer coefficient and porosity modified conductivity ratio on the instability of the system are investigated. The analysis shows that positive Soret mass flux triggers instability and positive Dufour energy flux enhances stability whereas their combined influence depends on the product of solutal Rayleigh number and Lewis number. It also reveals that cell width at the convection threshold gets affected only in the presence of both the cross-diffusion fluxes. Besides, asymptotic solutions for both small and large values of the inter-phase heat transfer coefficient and porosity modified conductivity ratio are found. An excellent agreement is found between the exact and asymptotic solutions.     

超声速飞行器壁板非线性颤振响应分析的时域法与频域法对比研究

为考查基于假设模态法在时域中开展壁板非线性颤振分析的可行性,在相同的参数下,分别采用时域方法和频域方法研究了超声速飞行器壁板的非线性颤振响应,并从壁板的颤振幅值、颤振频率和颤振型态三个方面对时域和频域分析结果的一致性作了较详细的比较。首先,基于von Karman应变－位移关系和Mindlin板理论建立考虑几何非线性的壁板力学模型,应用一阶活塞理论分析壁板上单面承受的超声速准定常气动力,基于虚功原理和有限单元法推导壁板的运动微分方程。然后,用壁板的线性固有模态作为假设模态,减缩系统的自由度而得到降阶模型。采用四阶龙格-库塔法对降阶模型作时域数值积分,得到壁板的非线性颤振响应。另一方面,假设壁板的极限环颤振为简谐振荡,可对壁板的非线性刚度作等效线性化处理,进而在频域中直接在有限元（未降阶）模型的基础上分析壁板的颤振幅值、颤振频率和颤振型态。数值分析表明,当极限环颤振为简谐振荡时,时域方法和频域方法的计算结果符合一致。本文最后讨论了时域法和频域法应用在壁板非线性颤振分析中各自的优点和局限性。     

Aeroelastic analysis of an idealized landing gear door in subsonic potential flow

Landing gear doors on aircraft have experienced flutter during preliminary flight testing. While designs vary widely, landing gear doors are typically plate-like structures with a relatively rigid actuator attached to their inside surface. To better understand the aeroelasticity of landing gear doors, this study investigates the aeroelastic stability of an idealized model. The model consists of a hinged plate with an interior constraint approximating the actuator attachment. The plate is subject to uniform flow, and an unsteady vortex lattice model is coupled to the structural model to predict critical flow velocities. The location and footprint area of the internal constraint, along with plate aspect and mass ratios, are varied to investigate a large parameter space. Results reveal that the critical flow speed and instability mechanism are sensitive to the postulated actuator placement. In general, flutter is the dominant mode of instability when the actuator is postulated in the leading quarter of the plate. In other postulated locations, divergence dominates. However, the exact shape and location of the boundary between flutter and divergence is configuration dependent and found to be especially sensitive to changes in aspect ratio.     

Post-critical analysis of highly deformable Joined Wings: The concept of snap-divergence as a characterization of the instability

Theoretical modeling of the aeroelastic static instability with the possibility of tracing post-critical branches (via arc length technique) in the framework of Joined Wings has never been presented before. A complete formulation of the numerical iterative method of solution of the aeroelastic equations is presented.The true critical condition is compared with the divergence speed evaluated by solving an eigenvalue problem about a steady state equilibrium, showing how this last approach may be unreliable and even nonconservative.This work also explores the theoretical implications of using mechanical loads to mimic the real loading conditions.A physical interpretation based on the aeroelastic effects, overconstrained nature of the system, and the bending/torsion coupling, is provided to interpret the conditions that lead to the snap-divergence.     

An aeroelastic model for composite rotor blades with straight and swept tips. Part I: Aeroelastic stability in hover

An analytical model for predicting the aeroelastic behavior of composite rotor blades with straight and swept tips is presented. The blade is modeled by beam type finite elements along the elastic axis. A single finite element is used to model the swept tip. The non-linear equations of motion for the finite element model are derived using Hamilton's principle and based on a moderate deflection theory and accounts for: arbitrary cross-sectional shape, pretwist, generally anisotropic material behavior, transverse shears and out-of-plane warping. Numerical results illustrating the effects of tip sweep, anhedral and composite ply orientation on blade aeroelastic behavior are presented. It is shown that composite ply orientation has a substantial effect on blade stability. At low thrust conditions, certain ply orientations can cause instability in the lag mode. The flap-torsion coupling associated with tip sweep can also induce aeroelastic instability in the blade. This instability can be removed by appropriate ply orientation in the composite construction.     

直升机旋翼/机体动稳定性研究进展

首先对直升机旋翼/机体动不稳定性问题的种类进行了简要概述,包括旋翼挥舞/变距、变距/摆振、挥舞/摆振和挥舞/摆振/变距耦合等孤立旋翼动不稳定性问题,以及直升机地面共振和空中共振等旋翼/机体耦合动不稳定性问题,然后分别从气动力与结构的高精度数值模型、动稳定性的计算分析方法和实验模型测试3 个方面详细介绍了直升机旋翼/机体动不稳定性问题的研究现状,并着重讨论了直升机旋翼/机体动稳定性分析技术最近的主要研究方向:耦合CFD(computational fluid dynamics)/CSD(computational structuraldynamics) 的直升机旋翼气弹动稳定性分析、复合材料旋翼动稳定性分析及其材料不确定性影响、带减摆器的旋翼/机体动稳定性分析和先进直升机构型的旋翼/机体动稳定性分析,最后对直升机旋翼/机体动稳定性研究的发展趋势进行了展望.     

High-order shear theory for static analysis of functionally graded plates with porosities

The bending responses of porous functionally graded (FG) thick rectangular plates are investigated according to a high-order shear deformation theory. Both the effect of shear strain and normal deformation are included in the present theory and so it does not need any shear correction factor. The equilibrium equations according to the porous FG plates are derived. The solution to the problem is derived by using Navier's technique. Numerical results have been reported and compared with those available in the open literature for non-porous plates. The effects of the exponent graded and porosity factors are investigated.     

Aeroelastic analysis of helicopter rotor blade in hover using an efficient reduced-order aerodynamic model

This paper presents a coupled flap–lag–torsion aeroelastic stability analysis and response of a hingeless helicopter blade in the hovering flight condition. The boundary element method based on the wake eigenvalues is used for the prediction of unsteady airloads of the rotor blade. The aeroelastic equations of motion of the rotor blade are derived by Galerkin's method. To obtain the aeroelastic stability and response, the governing nonlinear equations of motion are linearized about the nonlinear steady equilibrium positions using small perturbation theory. The equilibrium deflections are calculated through the iterative Newton–Raphson method. Numerical results comprising steady equilibrium state deflections, aeroelastic eigenvalues and time history response about these states for a two-bladed rotor are presented, and some of them are compared with those obtained from a two-dimensional quasi-steady strip aerodynamic theory. Also, the effect of the number of aerodynamic eigenmodes is investigated. The results show that the three-dimensional aerodynamic formulation has considerable impact on the determination of both the equilibrium condition and lead-lag instability.     

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

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

Instability of Streaming Viscoelastic Fluids in the Presence of ��Effective Interfacial Tension�� Through Porous Medium

The ??effective interfacial tension?? effect on the instability of the plane interface between two uniform, superposed, and streaming Rivlin?CEricksen viscoelastic fluids through a porous medium is considered. The case of two uniform streaming Rivlin?CEricksen viscoelastic fluids separated by a horizontal boundary is studied. In the absence of ??effective interfacial tension??, stability/instability of the system as well as perturbations transverse to the direction of streaming are found to be unaffected by the presence of streaming if the perturbations in the direction of streaming are ignored, whereas for perturbations in all other directions, there exists instability for a certain wave number range. ??Effective interfacial tension?? is able to suppress this Kelvin?CHelmholtz instability for small wavelength perturbations, the medium porosity reduces the stability range given in terms of a difference in streaming velocities.     

Aeroelastic effect on modal interaction and dynamic behavior of acoustically excited metallic panels

This paper details the study of the aeroelastic effect on modal interaction and dynamic behavior of acoustically excited square metallic panels with fully clamped edges using finite element method. The first-order shear deformation plate theory and von Karman nonlinear strain–displacement relationships are employed to consider the structural geometric nonlinearity caused by large vibration deflections. Piston aerodynamic theory and Gaussian white noise are used to simulate the aerodynamic load and the acoustic load, respectively. Motion equations are derived by the principle of virtual work in the physical coordinates and then transformed into the truncated modal coordinates with reduced orders. Runge–Kutta method is employed to obtain the system response, and the modal interaction mechanism is quantitatively valued by the modal participation distribution. Results show that in the pre-/near-flutter regions, in addition to the dominant fundamental resonant mode, the first twin companion antisymmetric modes can be largely excited by the aeroelastic coupling mechanism; thus, aeroelastic modal participation distribution and the spectrum response can be altered, while the dynamic behavior still exhibits linear random vibrations. In the post-flutter region, the dominant flutter motion can be enriched by highly ordered odd order super-harmonic motion occurs due to 1:1 internal resonances. Correspondingly, the panel dynamic behavior changes from random vibration to highly ordered motions in the fashion of diffused limit-cycle oscillations (LCOs). However, this LCOs motion can be affected by the intensifying acoustic excitation through changing the aeroelastic modal interaction mechanism. Accompanied with these changes, the panel can experience various stochastic bifurcations.     

铝面板厚度对Nomex蜂窝夹层结构冲击后弯曲性能影响的试验研究

通过低速冲击试验和四点弯曲试验,研究了铝面板厚度对Nomex蜂窝夹层结构抗冲击能力和剩余强度的影响。结果表明:在冲击荷载作用下,面板发生变形的区域大小随面板厚度增加而变大,当面板厚度大于0.5mm时,变形区域直径趋于稳定;无论试件是否受到过冲击,在弯曲载荷作用下,0.2mm厚面板发生芯格内屈曲失稳,而其他厚度面板均发生格间失稳;对无冲击损伤的结构,0.2mm厚面板弯曲强度显著低于其他厚度面板;对含冲击损伤的结构,0.2mm厚面板的剩余强度百分比最高。     

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

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