壁板颤振的分析模型、数值求解方法和研究进展

研究壁板颤振问题需要计及大挠度变形下结构的几何非线性效应,不仅涉及气动弹性稳定性,而且关心结构的非线性颤振响应.该文回顾了飞行器壁板颤振问题的国内外研究情况,评述了在壁板颤振研究中采用的分析模型、数值求解方法以及在理论分析和试验方面的研究成果,并提出了今后壁板颤振问题的4个研究方向.      

自适应非线性滤波及其在图像导航中的应用

本文首先介绍了飞行器末制导中采用图像导航的原理和方法，给出了系统模型。针对图像导航的实际情况，提出了自适应非线性Sage-Husa滤波，能在线估计系统的常值偏差及噪声方差。将滤波算法应用到图像导航中，进行仿真，滤波收敛，结果令人满意。     

考虑几何与材料及静风荷载的非线性因素的大跨径桥梁静风稳定分析法

针对现有的桥梁静风稳定分析方法中存在的问题，提出了增量与内外两重迭代相结合的新方法，并且考虑了结构几何、材料和静风荷载非线性。在上述方法的基础上，编制了桥梁非线性空气静力稳定分析程序BNAP，并进行了相应的算例分析，所得结果表明该方法具有计算稳定和速度快的优点。最后，以一座主跨1000米的斜拉桥为例，分析了结构几何非线性、材料非线性和静风荷载非线性对大跨径桥梁空气静力稳定性的影响。     

非线性被动隔振的若干进展

工程中航空航天、船舶与海洋结构物及其上装备和精密仪器易受极端环境干扰和破坏,使得非线性隔振理论在近十年来迅猛发展;针对日益严峻的隔振和抗冲击等要求,工程师和科学家们已发展出各种不同的非线性隔振系统,包括主动、半主动、被动和复合隔振.利用非线性改善的被动隔振兼具传统被动隔振的鲁棒性和主动隔振的高效性成为振动控制领域的先进技术.本文主要综述了非线性隔振理论和应用的近十年进展,包括非线性隔振设计、建模、分析、仿真和实验.在隔振系统的构建中,既考虑了刚度非线性又考虑了阻尼非线性;动力学响应的研究中,既有确定性分析又有随机分析.首先提出了适用于非线性隔振系统改进的评价方式;其次综述了高静态低动态刚度隔振及其加强形式非线性阻尼加强和双层非线性隔振,混沌反控制技术、内共振影响、非线性能量阱应用等振动机制利用型隔振和非线性隔振功能材料.最后,对非线性隔振研究发展的热点和关键性问题进行了分析和展望.     

以互信息为基础的广义相关系数及其应用

以往分析变量间相关性时，大多局限于仅反映线性相关程度的互相关函数或自相关函数上，但我们所研究的问题，却常常是非线性，因而对这种非线性结构上的广义相关程度进行度量显得更为重要。本文着重讨论了描述变量间这种广义相关程度的一种以互信息为基础的广义相关系数的概念、算法及操作中的技术问题，并实例其应用。     

一种微静电开关系统的非线性动力学

在微型电子机械（MEMS）中，存在着多种非线性因素，诸如：几何非线性、静电力的非线性以及空气阻力的非线性等，它们严重影响着微机械系统的行为。对于一种具有典型结构的微静电开关，考虑机械、静电、挤压气隙三部分的因素，导出了系统单自由度非线性运动方程（保留到三阶非线性），用数值方法分析了静电激励变化导致的分岔，发现倍周期分岔、逆倍周期分岔广泛存在于该类系统中。     

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

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

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

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

考虑气动力非线性影响的结构驰振分析

驰振现象对工程结构危害很大，因此开展结构驰振方面的理论研究十分重要，在一般的驰振分析中大多仅考虑气动力为线性时的情况，本文研究了气动力非线性对结构横向向动力失稳临界风速的影响，采用了以ｉ／ｖ的高次多项式形式表示的气动力模型，用等效线性法求角结构驰振的非线性微分运动方程，研究了非线性项的影响，为应用和研究提供参考。     

板壳结构中的非线性兰姆波

鉴于常规超声检测技术对分布式材料细微损伤和接触类结构损伤的检测效果不佳,近年来非线性超声技术逐渐引起广泛关注.超声波在板壳结构中通常以兰姆波的形式进行传播,然而由于兰姆波的频散及多模特性,使得非线性兰姆波的理论和实验研究进展缓慢.本文从经典非线性理论出发,总结了源于材料固有非线性诱发的非线性兰姆波的理论和实验两个方面的研究进展,井综述了兰姆波的二次谐波发生效应在材料损伤评价方面的若干应用;从接触声非线性理论出发,讨论了目前由于接触类结构损伤诱发的非线性兰姆波的研究现状.最后展望了非线性兰姆波的未来研究重点及发展趋势.     

Identification of Nonlinear Aeroelastic Systems Based on the Volterra Theory: Progress and Opportunities

The identification of nonlinear aeroelastic systems based on the Volterra theory of nonlinear systems is presented. Recent applications of the theory to problems in computational and experimental aeroelasticity are reviewed. Computational results include the development of computationally efficient reduced-order models (ROMs) using an Euler/Navier–Stokes flow solver and the analytical derivation of Volterra kernels for a nonlinear aeroelastic system. Experimental results include the identification of aerodynamic impulse responses, the application of higher-order spectra (HOS) to wind-tunnel flutter data, and the identification of nonlinear aeroelastic phenomena from flight flutter test data of the active aeroelastic wing (AAW) aircraft.     

飞机结构气动弹性分析与控制研究

随着主动控制技术的发展,飞机结构设计理念已由提高结构刚度的被动设计转变为随控布局的主动设计.主动设计理念不再刻意回避气动弹性问题,而是采用主动控制技术实时调节结构气动弹性,进而减轻结构重量、优化飞机性能. 在飞机随控布局主动设计中,必须深入分析结构与气流之间的耦合,才能更好发挥气动弹性主动控制技术的作用. 从20 世纪80 年代起,航空科技界对该问题进行了长期研究,对飞机结构-空气动力-主动控制相互耦合后的关键力学问题有了深入理解. 然而,已有研究多基于简化模型,导致研究结果难以直接应用于工程. 本文将针对气动弹性动态问题,综述空气动力非线性、控制面间隙非线性、时滞诱发失稳、颤振主动抑制、突风载荷减缓、风洞实验验证等方面的国内外研究进展,重点介绍近年来作者团队所提出的若干方法及相关算例和风洞实验. 最后,指出今后一个时期值得研究的若干气动弹性分析与控制问题.      

基于流形切空间插值的折叠翼参数化气动弹性建模

变体飞行器的气动弹性力学建模是当前先进飞行器设计的研究热点和难点.然而传统的气动弹性动力学建模方法对于具有结构参变特性的变体飞行器气动弹性力学研究存在建模效率低、计算复杂等问题.本研究提出了一种基于流形切空间插值的可折叠式变体机翼参数化气动弹性建模方法.首先,该方法建立若干个典型折叠角下的折叠翼结构有限元模型,通过流形...     

An efficient implementation of aeroelastic tailoring based on efficient computational fluid dynamics-based reduced order model

Current and future trends in the aerospace industry leverage on the potential benefits provided by lightweight materials that can be tailored to realize desired mechanical characteristics when loaded. For aircraft design, the deployment of aeroelastic tailoring is hindered by the need to re-compute, for any possible modification of the structure, the dependence of the aerodynamic field on the underlying structural properties. To make progress in this direction, the work presents a rapid computational fluid dynamics based aeroelastic tool which is built around a reduced order model for the aerodynamics that is updated for any modification of the structure by using the structural dynamics reanalysis method. The aeroelastic tailoring tool is demonstrated in transonic flow for the AGARD 445.6 wing, suitably modified with composite materials. It was found that the proposed method provides accurate engineering predictions for the aeroelastic response and stability when the structure is modified from the baseline model.     

An sequential optimization and aeroelastic constraint transformation method for strength-aeroelastic comprehensive design

With the development of modern aircraft technology, aeroelasticity plays a more and more crucial role in aircraft structural design. However, low efficiency of present aeroelastic analysis and optimization methods makes it difficult to apply in engineering practice. This paper presents a sequential optimization and aeroelastic constraint transformation method (SOACTM) for comprehensive design of airplane wings with strength and aeroelastic constraints. Optimization with structural strength constraint and aeroelastic constraint is transformed into a serial of cycles of decoupled structural strength sub-optimizations and aeroelastic sub-optimizations based on sequential optimization strategy. In structural strength sub-optimization, structural strength constraint is translated along its normal direction to make optimal design point satisfying aeroelastic constraint. And the goal of aeroelastic sub-optimization is to find the translational distance of structural strength constraint. Aeroelastic constraint is transformed to equivalent structural strength constraint via above approach. In this way, number of aeroelastic analyses in SOACTM is less than that in traditional optimization method and total computational time decreases. SOACTM is verified based on two examples. Traditional optimization method is applied for the sake of validation. The results demonstrated the accuracy and efficiency of SOACTM for wing comprehensive optimization considering both structural strength and aeroelastic constraints.     

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

Implications of cubic physical/aerodynamic non-linearities on the character of the flutter instability boundary

The present paper deals with a study of the benign and catastrophic characters of the flutter instability boundary of 2-D lifting surfaces in a supersonic flow field. The objectives of this work are: (i) to contribute to a better understanding of the implications of aerodynamic and physical non-linearities on the character of the flutter boundary and (ii), to outline the effects exerted in the same respect by some important parameters of the aeroelastic system. With the aim of addressing this problem, the method based on the First Liapunov Quantity is used to study the bifurcational behavior of the aeroelastic system in the vicinity of the flutter boundary. The expected outcomes of this study are: (a) to greatly enhance the scope and reliability of the aeroelastic analysis and design criteria of advanced aircraft and, (b) to provide a theoretical basis for the analysis of more complex non-linear aeroelastic systems.     

Effects of actuator nonlinearity on aeroelastic characteristics of a control fin

The influences of actuator nonlinearities on actuator dynamics and the aeroelastic characteristics of a control fin were investigated by using iterative V-g methods in subsonic flows; in addition, the doublet-hybrid method (DHM) was used to calculate unsteady aerodynamic forces. The changes of actuator dynamics induced by nonlinearities, such as backlash or freeplay, and the variations of flutter boundaries due to the changes of actuator dynamics were observed. Results show that the aeroelastic characteristics can be significantly dependent on actuator dynamics. Thus, the actuator nonlinearities may play an important role in the nonlinear aeroelastic characteristics of an aeroelastic system. The present results also indicate that it is necessary to seriously consider the influence of actuator dynamics on the flutter characteristics at the design stage of actuators to prevent aeroelastic instabilities of aircraft or missiles.     

Aircraft T-tail flutter predictions using computational fluid dynamics

The paper presents the application of computational aeroelasticity (CA) methods to the analysis of a T-tail stability in transonic regime. For this flow condition unsteady aerodynamics show a significant dependency from the aircraft equilibrium flight configuration, which rules both the position of shock waves in the flow field and the load distribution on the horizontal tail plane. Both these elements have an influence on the aerodynamic forces, and so on the aeroelastic stability of the system. The numerical procedure proposed allows to investigate flutter stability for a free-flying aircraft, iterating until convergence the following sequence of sub-problems: search for the trimmed condition for the deformable aircraft; linearize the system about the stated equilibrium point; predict the aeroelastic stability boundaries using the inferred linear model. An innovative approach based on sliding meshes allows to represent the changes of the computational fluid domain due to the motion of control surfaces used to trim the aircraft. To highlight the importance of keeping the linear model always aligned to the trim condition, and at the same time the capabilities of the computational fluid dynamics approach, the method is applied to a real aircraft with a T-tail configuration: the P180.     

Assessment of low-altitude atmospheric turbulence models for aircraft aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.