考虑机翼几何非线性的气动弹性建模与分析

大展弦比柔性机翼结构重量轻、气动效率高,广泛应用于高空长航时无人机(UAVs)。飞行过程中,这类机翼在气动力作用下发生大变形,线性结构模型不再适用,需要建立考虑几何大变形的结构模型。采用牛顿力学方法推导了考虑结构几何非线性的机翼结构动力学模型,该方法推导过程简洁、物理意义明确,可以与Hodges基于哈密顿原理的推导方法相互补充,相互验证。为了能够更准确地求解大展弦比柔性机翼的非定常气动力,建立了能够考虑机翼三维效应且适用于机翼空间大变形的非定常气动力模型。基于建立的非线性结构模型和非定常气动力模型,采用松耦合方法建立了非线性气动弹性模型,并通过算例验证了气弹模型的准确性。研究结果表明,大展弦比柔性机翼颤振速度对来流迎角和机翼的展长均较为敏感;当来流速度大于颤振速度时,由于几何非线性,机翼振动并未发散而是形成稳定的极限环振荡(LCO);随着来流速度进一步增加,机翼再次穿过临界稳定点,由不稳定系统变为稳定系统,直到随着速度的增加系统再次达到临界稳定状态。     

基于CFD／CSD的大展弦比机翼气动弹性分析

大展弦比飞机在飞行过程中受气动载荷影响,其大展弦比机翼产生弯曲和扭转变形,这种弹性变形严重影响飞机的飞行性能和飞行安全,不能将此种飞机机翼当作传统的刚性机翼进行气动分析。针对一大展弦比机翼,采用气动／结构耦合分析方法,利用计算流体动力学（CFD）软件CFX和计算结构动力学（CSD）软件ANSYS联合求解,研究了在不同载荷情况下大展弦比机翼静气动弹性变形对机翼气动特性的影响。结果表明,大展弦比无人机机翼受载变形后升阻比降低,升力下降明显,阻力略有上升,机翼翼尖容易失速。     

侧向随动力作用下大展弦比柔性机翼的稳定性

 随动力能够诱发弹性结构发生颤振失稳。以侧向随动力和集中质量分别模拟发动机推力和外挂质量,考虑机翼垂直弯曲-扭转刚度比、集中质量大小、侧向随动力和集中质量的位置以及机翼后掠角和上反角的影响,研究了受侧向随动力作用的大展弦比柔性机翼的气动弹性稳定性。数值模拟所采用的大展弦比柔性机翼非线性气动弹性模型耦合了几何精确完全本征运动梁模型和ONERA动失速非定常气动力模型,该模型考虑了几何非线性、动失速和材料各向异性。模拟结果表明,侧向随动力对机翼颤振可以具有稳定作用,其具体表现依赖于若干变参数的影响,如：减小机翼垂直弯曲-扭转刚度比;发动机吊舱靠近翼根布置;使发动机推力作用点在法向上与机翼弹性轴靠近;单纯的集中质量避免布置在柔性机翼中部,且布置在机翼弹性轴之前或下方,这些设计或布置均有利于提高带发动机吊舱/有效载荷外挂的柔性机翼的气动弹性稳定性。     

基于CR理论的大柔性机翼几何非线性结构建模

大柔性机翼在气动载荷的作用下,将产生显著的弹性变形,常规线弹性理论的小变形假设不再成立,需要采用能够考虑几何非线性效应的结构模型进行求解。基于CR(Co-Rotational)共旋转有限元理论,把几何非线性大变形分解为刚体的旋转和平移及局部坐标系下的弹性变形,建立了适用于大柔性机翼几何非线性变形描述的结构模型。以大柔性悬臂梁为例,采用载荷增量法,研究了集中弯矩作用下的非线性变形特征,对静力学方法进行了验证,并讨论了耦合加载几何非线性变形特征;以类"太阳神"太阳能布局无人机(UAV)为例,研究了其几何非线性大变形特性。     

大展弦比柔性机翼气动弹性分析中的气动力方法研究进展

近20年来长航时飞行的需求强烈,大柔性飞行器的几何非线性气动弹性问题逐渐凸显,使得气动弹性力学面临新的挑战。本文针对大展弦比大变形的柔性飞行器,调研和分析了目前几何非线性气动弹性工程研究领域中主要使用的气动建模方法,着重介绍基于片条理论、面元法和计算流体力学(Computational Fluid Dynamics,CFD)技术等气动建模方法在静、动气动弹性分析中的主要特点、研究现状与应用状况,并对大展弦比大变形机翼的气动弹性分析中气动力方法的发展提出若干建议,供气动弹性基础研究和工程应用研究人员参考。     

柔性机翼阵风响应与被动减缓研究

高空长航时无人机的机翼展弦比大、柔性较强,飞行过程中极易受到阵风的影响。文章以几何精确本征理论建立结构模型,耦合 Pitt-Peters动力入流理论建立柔性机翼非线性气弹模型,研究了柔性机翼阵风响应以及翼尖被动阵风减缓效应。采用空间 -时间平行的有限元离散方法,将气弹方程转化为一阶微分代数方程,Newton-Raphson和 Generalized-α算法分别用于静态变形和动态响应的求解,通过算例研究了离散阵风载荷下柔性机翼的阵风响应,结果表明翼尖被动阵风减缓装置对机翼变形有明显的减缓效果。     

大变形状态机翼振动试验与气动弹性分析

本文针对结构大变形引起的几何非线性对于结构动力学特性的影响作了计算和试验研究,设计了一个试验,用于体现大展弦比平板机翼结构大变形对于动力学特性的改变。并计算分析了这种大变形对于该机翼气动弹性颤振特性的影响,并参考了一个双梁骨架结构机翼的计算分析结果。经分析发现几何非线性因素对于这两种机翼结构的影响体现在:第一,使得固有模态中的扭转模态频率大幅降低;第二,使得颤振特性发生变化。     

无人机大展弦比直机翼单双机身布局对结构刚度的影响

为了探讨单、双机身布局下的大展弦比直机冀之间的结构刚度问题以及双机身大展弦比直机翼布局中基于翼面弯曲刚度下机身与机翼的最佳布置问题，采用工程梁方法对这一问题进行简化与建模。理论分析与数值优化结果表明：就结构刚度而言，双机身布局优于单机身布局；针对双机身布局，通过优化给出了机身与机翼之间的相对位置，在此工作的基础上，研究了冀面扭转刚度，并由此得出：就结构的弯曲和扭转刚度而言，双机身布局具有一定的优势。     

考虑柔性气动力的机翼结构优化方法研究

分析研究了大展弦比弹性机翼结构静弹性变形对气动载荷的影响，在传统结构优化方法基础上进行了弹性机翼的结构优化设计方法研究。实现了在真实飞行条件下考虑静气动弹性变形影响的大展弦比弹性机翼的结构优化设计。结构优化设计算例验证了方法的可行性，优化结果表明结构重量收益明显。     

大展弦比机翼的气动弹性问题探讨

在考察国内外高空长航时飞机研究的基础上，对大展弦比机翼的气动弹性问题进行了探讨；分析了结构非线性对大展弦比机翼的气动弹性和飞行载荷的影响；提出了大展弦比机翼气动弹性问题的研究内容，并指出其技术难点。     

Gust response analysis and wind tunnel test for a high-aspect ratio wing

A theoretical nonlinear aeroelastic response analysis for a flexible high-aspect ratio wing excited by harmonic gust load is presented along with a companion wind tunnel test. A multidisci-plinary coupled numerical calculation is developed to simulate the flexible model wing undergoing gust load in the time domain via discrete nonlinear finite element structural dynamic analysis and nonplanar unsteady vortex lattice aerodynamic computation. A dynamic perturbation analysis about a nonlinear static equilibrium is also used to determine the small perturbation flutter bound-ary. A novel noncontact 3-D camera measurement analysis system is firstly used in the wind tunnel test to obtain the spatial large deformation and responses. The responses of the flexible wing under different static equilibrium states and frequency gust loads are discussed. The fair to good quanti-tative agreements between the theoretical and experimental results demonstrate that the presented analysis method is an acceptable way to predict the geometrically nonlinear gust response for flex-ible wings.     

Static aeroelastic analysis of very flexible wings based on non-planar vortex lattice method

A rapid and efficient method for static aeroelastic analysis of a flexible slender wing when considering the structural geometric nonlinearity has been developed in this paper. A non-planar vortex lattice method herein is used to compute the non-planar aerodynamics of flexible wings with large deformation. The finite element method is introduced for structural nonlinear statics analysis. The surface spline method is used for structure/aerodynamics coupling. The static aeroelastic characteristics of the wind tunnel model of a flexible wing are studied by the nonlinear method presented, and the nonlinear method is also evaluated by comparing the results with those obtained from two other methods and the wind tunnel test. The results indicate that the traditional linear method of static aeroelastic analysis is not applicable for cases with large deformation because it produces results that are not realistic. However, the nonlinear methodology, which involves combining the structure finite element method with the non-planar vortex lattice method, could be used to solve the aeroelastic deformation with considerable accuracy, which is in fair agreement with the test results. Moreover, the nonlinear finite element method could consider complex structures. The non-planar vortex lattice method has advantages in both the computational accuracy and efficiency. Consequently, the nonlinear method presented is suitable for the rapid and efficient analysis requirements of engineering practice. It could be used in the preliminary stage and also in the detailed stage of aircraft design.     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

基于本征方程的柔性机翼非线性颤振分析

柔性机翼在气动载荷作用下常常会产生较大的变形,颤振特性会随之发生变化,针对此问题线性理论常常难以进行合理的预测。以几何精确本征梁模型建立了机翼的运动方程,耦合ONERA-EDlin非线性气动模型,建立了柔性机翼的非线性气动弹性分析模型。利用Newton-Raphson和Backward-Differentiation-Formula(BDF)分别求解机翼的静态变形和动态响应,基于机翼平衡位置附近的线性化方程来判断系统的稳定性,进而确定颤振临界速度。通过算例验证了模型的准确性,并分析了不同刚度、后掠角、机翼安装角等参数对颤振速度的影响。     

压电纤维复合材料驱动的机翼动态形状控制

利用压电材料实现柔性机翼的主动形状控制,能够有效提高机翼结构效率和气动性能;要实现连续、光滑的高精确形状控制效果,机翼变形过程必须满足一定的动态要求。本文利用在上下翼面反对称铺设的新型压电纤维复合材料——宏纤维复合材料(MFC)提供驱动力矩,研究了机翼扭转变形的前馈轨迹跟踪控制。首先建立了机翼结构有限元模型和气动力载荷模型,采用载荷比拟法得到MFC作动器的控制载荷向量,给出了气动弹性控制方程及其状态空间表达形式。为跟踪预设的变形参考轨迹,以跟踪误差的时域积分为目标函数,对MFC作动器的电压加载历程进行了优化设计。结果表明,采用规划后的电压加载历程,机翼气动弹性响应很好地跟踪了预期参考轨迹,实现了连续、光滑的动态形状控制效果,提高了控制精度。     

Nonlinear Aeroelastic Response of High-aspect-ratio Flexible Wings

The aeroelastic analysis of high-altitude, long-endurance (HALE) aircraft that features high-aspect-ratio flexible wings needs take into account structural geometrical nonlinearities and dynamic stall. For a generic nonlinear aeroelastic system, besides the stability boundary, the characteristics of the limit-cycle oscillation (LCO) should also be accurately predicted. In order to conduct nonlinear aeroelastic analysis of high-aspect-ratio flexible wings, a first-order, state-space model is developed by combining a geometrically exact, nonlinear anisotropic beam model with nonlinear ONERA (Edlin) dynamic stall model. The present investigations focus on the initiation and sustaining mechanism of the LCO and the effects of flight speed and drag on aeroelastic behaviors. Numerical results indicate that structural geometrical nonlinearities could lead to the LCO without stall occurring. As flight speed increases, dynamic stall becomes dominant and the LCO increasingly complicated. Drag could be negligible for LCO type, but should be considered to exactly predict the onset speed of flutter or LCO of high-aspect-ratio flexible wings.     

基于高精度模态气动力的跨声速静弹高效分析方法

跨声速静弹分析一直是工程设计中的难点问题,以模态坐标系下的线性静弹方程为基础,提出了基于高精度模态气动力的跨声速静弹高效分析方法,该方法仍需求解线性静弹方程,但对于其中关键的模态变形引起的弹性气动力增量,采用由结构变形到气动力的单向计算流体动力学(Computational Fluid Dynamics,CFD)/计算结构动力学(Computational Structural Dynamics,CSD)耦合方法获得,实现了高效线性方法与高精度CFD/CSD耦合方法的有效融合。以某小展弦比机翼基本状态、舵偏状态以及某型战斗机跨声速副翼效率的静弹分析为例,对比分析了本文方法、经典线性方法、CFD/CSD耦合方法的计算结果以及某型机的试飞辨识结果。分析结果表明,所提出的方法在计算效率、精度和鲁棒性方面具备综合优势,具有较高的工程应用价值。     

Time-domain nonlinear aeroelastic analysis and wind tunnel test of a flexible wing using strain-based beam formulation

A theoretical formulation for time-domain nonlinear aeroelastic analysis of a flexible wing model is presented and validated by wind tunnel tests. A strain-based beam model for nonlinear structural analysis is combined with the Unsteady Vortex Lattice Method (UVLM) to form the complete framework for aeroelastic analysis. The nonlinear second-order differential equations are solved by an implicit time integration scheme that incorporates a Newton-Raphson sub-iteration technique. An advanced fiber optic sensing technique is firstly used in a wind tunnel for measuring large structural deformations. In the theoretical study, the nonlinear flutter boundary is determined by analyzing the transient response about the nonlinear static equilibrium with a series of flow velocities. The gust responses of the wing model at various gust frequencies are also studied. Comparisons of the theoretical and experimental results show that the proposed method is suitable for determining the nonlinear flutter boundary and simulating the gust response of flexible wings in the time domain.