气动伺服弹性系统不确定性建模与鲁棒稳定性

 气动伺服弹性系统在不确定性摄动下的鲁棒稳定性问题对于带有自动控制系统的弹性飞行器是非常关键的。利用线性分式变换形式,考虑参数和非参数不确定性的摄动,如广义刚度、广义质量、广义非定常气动力以及伺服系统的不确定性摄动,由各子系统到整个闭环系统依次建立气动伺服弹性的状态空间模型,并应用结构奇异值μ方法分析了系统的鲁棒稳定性。两个算例表明了该建模方法的方便适用以及μ方法在气动伺服弹性鲁棒稳定性分析中的应用前景。     

大展弦比机翼的静气动弹性分析

大型飞机普遍采用的大展弦比机翼,其气动弹性问题显得尤为突出。通过求解基于三维结构网格的Euler/边界层耦合方程组,采用结构模态分析法,对HIRENASD机翼的静气动弹性进行了数值模拟。计算得到了机翼静弹变形前和变形后的压力分布以及不同攻角下翼尖变形值。计算结果和实验值对比分析,验证了Euler/边界层耦合方程组求解器对静气动弹性数值模拟的准确性,表明所发展的数值模拟方法可以用于大展弦比机翼静气动特性的分析研究。     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

航空发动机结构参数和非结构参数不确定系统鲁棒H∞输出反馈控制

分析了航空发动机混合不确定性系统的鲁棒稳定性和鲁棒性能;引入结构奇异值μ综合法,研究了结构参数和非结构参数混合摄动下的鲁棒H∞输出反馈控制;对航空发动机不确定性系统设计的控制器在双转子涡喷发动机气动热力学非线性模型上。      

复合材料大展弦比机翼的弹性变形对飞机纵向气动导数的影响

针对复合材料大展弦比机翼进行了静气动弹性特性的分析研究。计算应用MSC／NASTRAN软件，以某飞机复合材料大展弦比机翼的静力分析有限元模型为基础，对模型进行合理的修改，建立起静气动弹性分析的有限元模型，使用结构柔度法和亚音速偶极子格网法进行静气动弹性计算，重点分析了纵向气动导数以及相应的弹性修正系数随马赫数、速压变化的特点，总结出复合材料大展弦比机翼的一些气动弹性特性。     

大展弦比机翼跨声速静气动弹性风洞试验

基于静气动弹性风洞试验研究了某翼身组合体的跨声速静气动弹性效应。试验结果表明:在设计巡航点,静气动弹性对大展弦比超临界机翼的气动特性影响明显,可使机翼的升力系数降低21%、升阻比增加8%、焦点前移约1%bA;在超过巡航马赫数后,静气动弹性效应使得机翼气动特性有恶化的趋势。跨声速时,马赫数和速压对机翼的静气动弹性效应具有较大影响,且影响规律呈复杂非线性,难以依据现有理论分析准确预计。     

静气动弹性模型高速风洞试验研究北大核心CSCD

对采用复合材料玻璃纤维、碳纤维加工的静气动弹性模型进行了高速风洞试验研究,测试了模型的柔度矩阵、气动力、表面压力、弯/扭应变信号及弯/扭变形,为静气动弹性模型刚度试验、弯/扭应变信号测量、模型变形视频测量(VMD)及风洞总压控制等静气动弹性风洞试验能力的提高积累了经验,为飞行器静气动弹性研究提供了良好的试验平台。研究表明:静气动弹性模型较刚性模型升力线斜率及襟副翼效率下降、气动焦点前移;静气动弹性模型与刚性模型表面压力差异明显;在小迎角范围内,静气动弹性机翼模型弯/扭应变信号随迎角增加基本呈线性变化;在正迎角时,大展弦比后掠机翼静气动弹性模型的剖面扭转变形使有效迎角减小,剖面越靠近翼尖弯/扭变形越大。     

静气动弹性模型高速风洞试验研究

对采用复合材料玻璃纤维、碳纤维加工的静气动弹性模型进行了高速风洞试验研究,测试了模型的柔度矩阵、气动力、表面压力、弯/扭应变信号及弯/扭变形,为静气动弹性模型刚度试验、弯/扭应变信号测量、模型变形视频测量(VMD)及风洞总压控制等静气动弹性风洞试验能力的提高积累了经验,为飞行器静气动弹性研究提供了良好的试验平台。研究表明:静气动弹性模型较刚性模型升力线斜率及襟副翼效率下降、气动焦点前移;静气动弹性模型与刚性模型表面压力差异明显;在小迎角范围内,静气动弹性机翼模型弯/扭应变信号随迎角增加基本呈线性变化;在正迎角时,大展弦比后掠机翼静气动弹性模型的剖面扭转变形使有效迎角减小,剖面越靠近翼尖弯/扭变形越大。     

TRIP软件的静气动弹性计算模块开发及精度验证

在结构网格流场求解软件TRIP基础上,发展了一套静气动弹性计算模块。首先简要介绍了静气动弹性计算模块总体架构、构成单元和耦合策略,然后详细介绍了构成单元中采用的一些数值计算方法。最后利用大展弦比机翼、DLR-F6翼身组合体模型和HIRENASD机翼模型对静气动弹性计算模块进行了测试和精度验证。大展弦比机翼计算结果表明柔度矩阵、模态叠加和有限元数值求解在模态选取合适的情况下能够得到一致的变形计算结果。DLR-F6和HIRENASD模型不同计算软件的结果一致性表明本文静气动弹性计算模块采用了正确的算法流程,而计算与试验的结果一致性表明静气动弹性计算模块具有很好的预测精度。     

NNW-FSI软件静气动弹性耦合加速策略设计与实现

在国家数值风洞（NNW）工程的资助下,依托NNW-FSI流固耦合模拟软件平台,从气动载荷作用下飞行器结构静变形大小与收敛过程无关的物理机制出发,基于变形增量叠加的方式,设计和实现了一种静气动弹性耦合加速策略,通过松弛因子对耦合迭代的收敛过程进行调整。结合超大展弦比无人机和CHN-T1模型两种不同外形,开展了不同松弛因子下的静气动弹性耦合数值模拟,对耦合加速策略的参数影响和加速效果进行了测试和评估。从计算误差控制角度对松弛因子加速耦合迭代收敛的作用机制进行了理论分析,弄清3种类型静气动弹性耦合模拟过程中松弛因子发挥的作用,并给出了松弛因子选取范围的建议。静气动弹性耦合模拟和理论分析结果表明,针对不同类型的静气动弹性耦合问题,选取合适的松弛因子,能够达到抑制振荡并加速收敛的效果。     

Match Point Solution for Robust Flutter Analysis in Constant-Mach Prediction

This paper presents a method for robust flutter computation which uses flight altitude as the perturbation variable in order to obtain a match point solution. The air density and sound speed of standard atmosphere model are approximated as the polynomial function of altitude, such that the flight altitude becomes the single perturbation variable that describes the aeroelastic system. The uncertainties of generalized stiffness and damping are considered and the uncertain aeroelastic system can be formulated as linear fractional transformation （LFT） representation which is suitable for/.t analysis framework. Finally, the match point solutions of robust flutter margins can be computed with structured singular value （SSV） theory. The robust flutter analysis method provided in this paper is suitable for constant-Mach flight flutter test and provides valuable reference for flight envelope expansion.     

智能结构不确定参数系统振动控制及其摄动分析

 结构建模中通常考虑不确定性等因素以确保系统及其控制系统具有良好的鲁棒性,由于参数不确定性引起的系统参数的变化将导致系统性能退化,甚至影响系统内部稳定性,所以不确定性概念在工程结构的分析与设计中起到重要的作用。研究了具有不确定参数系统鲁棒性理论,提出了抑制系统振动的控制规律;基于矩阵摄动法讨论了不确定参数对智能结构系统的影响,并利用不确定性凸模型理论分析了智能结构具有不确定参数系统稳定性的问题,提出了分析含不确定参数系统鲁棒性的方法。算例说明该方法的有效性。     

含不确定参数的复合材料壁板热颤振分析

研究了含有不确定结构参数的壁板颤振问题,利用vonKarman大变形应变-位移关系、气动力活塞理论和准定常热应力理论建立了复合材料壁板颤振的气动弹性力学模型,考虑在壁板颤振分析模型中存在的不确定参数,将其用区间向量定量化,基于区间扩张理论和Taylor级数展开,并结合有限元计算方法,提出了区间分析的方法来估计含有不确定参数的壁板结构颤振临界风速的区间,以及发生极限环振动时振幅的变化区间。通过数值算例,将本文提出的壁板颤振的区间有限元模型与随机有限元模型进行了比较,显示了本文方法的有效性和可行性。这种方法的优点是只需要知道不确定参数的所在范围界限,为解决含有不确定参数的壁板颤振这类复杂的气动弹性动力学问题提供了一个途径。     

一种基于光固化快速成型的飞机静弹性风洞试验模型设计方法

高速风洞静弹性模型设计和制造是静弹性风洞试验的一个关键。为解决模型设计周期长、制造费用高等问题,提出了一种基于立体光固化快速成型面向高速风洞大展弦比机翼静弹性模型研制方法。基于机翼刚度分布相似参数,采用机翼钢梁骨架和树脂蒙皮组合结构,通过优化结构尺寸完成静弹性模型结构设计;使用机械加工和快速成型技术完成模型制造,并通过地面刚度试验对加工模型进行了刚度分布验证。风洞试验结果表明：基于立体光固化成型技术设计和制造的静弹性风洞试验模型工程实用、可行,与传统静弹性模型研制过程相比,具有研制周期短、成本低而且不存在因填充物带来附加刚度的显著优势。     

复材机翼气弹特性工程化设计分析技术

本文在静气弹的基本原理分析基础上,发展了机翼气弹特性工程化设计方法：基于参数化的机翼气动弹性模型,研究了复合材料机翼整体蒙皮重要设计参数对舵效、静变形和静强度的影响,采用变参分析方法进行了灵敏度分析,获取影响舵效的关键参数变化趋势,用以指导舵效气动弹性剪裁优化设计。对优化设计结果进行了工程化的处理并进行校核,获得了工程可用的蒙皮铺层设计方案,在改善舵效使其满足飞机性能要求的同时,还获得了满意的减重105 kg（总质量的14%）效益。     

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

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

一种计及静气动弹性变形影响的跨声速机翼气动优化设计方法研究

耦合流场控制方程和结构静平衡方程求解得到大展弦比跨声速弹性机翼气动性能,分析研究了结构静弹性变形对气动载荷的影响.在此基础上以机翼典型剖面外形和机翼型架外形为设计变量进行了基于控制理论的跨声速弹性机翼气动优化设计方法研究,实现了在真实飞行条件下考虑静气动弹性变形影响的大展弦比跨声速弹性机翼一体化优化设计.优化设计算例结果表明所发展的方法是成功的,优化得到的弹性机翼在满足升力系数约束条件下提高了升阻比.     

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

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

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.