应用惯性释放方法的静气动弹性发散分析

翼面发散作为静气动弹性的重点分析对象,传统的分析方法通常将翼面根部固支约束,不考虑刚体模态的影响。根据咨询通告AC 25.629-1B 的要求,刚体模态或短周期模态也可能会造成发散。因此,在民用飞机静气动弹性分析及适航符合性验证中,要计及刚体模态对发散的影响。通过对惯性释放方法的研究,在发散分析中计及刚体模态,给出考虑刚体模态的模态法发散分析方法;以某型民用飞机为例,应用本文方法分别计算空机状态、不同燃油和商载的静气动弹性发散结果,并和不考虑刚体模态的模态法、考虑/不考虑刚体模态的颤振法的发散计算结果进行对比分析。结果表明：考虑刚体模态后的发散速度有所提高,临界发散模态不变,发散速度随燃油和商载的增加而降低。     

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

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

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

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

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

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

基于非结构网格的跨声速静气动弹性计算

基于非结构网格开发了静气动弹性计算程序,通过气动力与弹性变形的交替迭代计算实现流固耦合。气动力求解采用RANS方程,结构变形计算采用柔度系数法,通过面样条(IPS)方法进行柔度系数插值,使用改进弹簧近似法实现网格变形。通过对AGARD 445.6机翼的静气弹计算验证了方法的正确性,同时对机翼的静气动弹性特性进行了分析,计算表明,弹性变形使其升力线斜率下降、焦点前移。     

基于非线性气动力的飞机飞行载荷计算

为了研究气动力非线性对飞机静气动弹性特性及飞行载荷的影响,发展了一种基于非线性气动力的飞机静气弹及飞行载荷计算方法。通过引入CFD计算结果,实现对线性方法中的气动力影响系数矩阵的非线性修正,采用模态法求解静气弹配平方程得到飞机非线性的气弹变形及飞行载荷。以某翼身组合体构型为计算算例,分别针对刚性模型和柔性模型,基于非线性气动力完成了飞机静气动弹性分析及载荷计算,并与线性结果进行了对比分析。结果表明,非线性方法可对大迎角范围内飞机静气弹变形及飞行载荷做出较合理的预测。     

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

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

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

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

考虑物理参数摄动的静气动弹性鲁棒分析

为了研究静气动弹性系统在不确定性摄动下的稳定性和性能,建立了一种考虑物理参数摄动的静气动弹性鲁棒分析方法。从静气动弹性系统的物理方程出发,应用摄动理论和线性分式变换推导了由物理参数摄动引起的广义刚度、结构模态和广义定常气动力的摄动模型,然后得到了静气动弹性系统不确定性模型。将结构奇异值μ分析推广应用于静气动弹性鲁棒性分析。以一大展弦比长直机翼为例,分析了机翼前梁的材料弹性模态受到一定摄动时的鲁棒静发散稳定性和静气动弹性变形性能。数值结果验证了该方法是准确有效的。     

基于静气动弹性效应的飞机型架外形修正方法研究

民用客机主要针对巡航外形设计,这就需要通过静气动弹性修正得到型架外形,保证型架外形在设计巡航状态飞行时,气动性能能够恢复到巡航设计状态。本文基于混合网格的N-S方程和结构柔度矩阵方法的耦合求解,发展了大飞机静气动弹性计算方法。根据静气动弹性变形过程的可逆性,结合飞行器质量分布、发动机推力以及巡航飞行重量等设计参数,发展了考虑静气动弹性效应的飞机型架外形修正方法。利用所发展的修正方法,获得了某大型客机巡航飞行的型架修正外形,同时分析了其它飞行工况的弹性气动特性,并通过算例比较验证了本文发展方法的有效性。     

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

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

Static aeroelastic analysis including geometric nonlinearities based on reduced order model

This paper describes a method proposed for modeling large deflection of aircraft in non-linear aeroelastic analysis by developing reduced order model (ROM). The method is applied for solving the static aeroelastic and static aeroelastic trim problems of flexible aircraft containing geo-metric nonlinearities; meanwhile, the non-planar effects of aerodynamics and follower force effect have been considered. ROMs are computational inexpensive mathematical representations com-pared to traditional nonlinear finite element method (FEM) especially in aeroelastic solutions. The approach for structure modeling presented here is on the basis of combined modal/finite ele-ment (MFE) method that characterizes the stiffness nonlinearities and we apply that structure mod-eling method as ROM to aeroelastic analysis. Moreover, the non-planar aerodynamic force is computed by the non-planar vortex lattice method (VLM). Structure and aerodynamics can be cou-pled with the surface spline method. The results show that both of the static aeroelastic analysis and trim analysis of aircraft based on structure ROM can achieve a good agreement compared to anal-ysis based on the FEM and experimental result.     

基于响应面插值的非线性气动弹性计算

飞翼、大展弦比低雷诺数气动布局容易在小迎角的条件下出现气流分离,会带来明显的非线性气动力问题,同时气动弹性带来的影响亦不可忽略。针对此类布局提出一种建立基于径向基函数插值的非线性压力系数分布模型的方法,利用径向基函数插值建立面元上压力系数对迎角导数的响应面,将压力系数积分并通过无限板样条(IPS)方法进行气动结构多次迭代插值实现非线性气弹分析。结果验证了该方法对于静气动弹性分析的有效性,同时能准确地反映弹性带来气动效率的降低和变形对升力及阻力的影响。     

Aeroelastic prediction and analysis for a transonic fan rotor with the “hot” blade shape

This paper focuses on aeroelastic prediction and analysis for a transonic fan rotor with only its “hot” (running) blade shape available, which is often the case in practical engineering such as in the design stage. Based on an in-house and well-validated CFD solver and a hybrid structural finite element modeling/modal approach, three main aspects are considered with special emphasis on dealing with the “hot” blade shape. First, static aeroelastic analysis is presented for shape transformation between “cold” (manufacturing) and “hot” blades, and influence of the dynamic variation of “hot” shape on evaluated aerodynamic performance is investigated. Second, implementation of the energy method for flutter prediction is given and both a regularly used fixed “hot” shape and a variable “hot” shape are considered. Through comparison, influence of the dynamic variation of “hot” shape on evaluated aeroelastic stability is also investigated. Third, another common way to predict flutter, time-domain method, is used for the same concerned case, from which the predicted flutter characteristics are compared with those from the energy method. A well-publicized axial-flow transonic fan rotor, Rotor 67, is selected as a typical example, and the corresponding numerical results and discussions are presented in detail.     

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

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

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

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

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.     

静气弹中非线性气动力建模方法与分析

大展弦比低雷诺数气动布局容易较早出现气流分离,会带来明显的非线性气动力问题。针对此类布局提出了一种建立基于Kriging插值的非线性压力系数分布模型的方法。从Navier-Stokes方程计算的不同状态下飞机的压力系数中提取不同坐标的系数。利用Kriging插值函数建立CFD压力系数对迎角导数的响应面,将插值结果代入偶极子网格法(Double-Lattice Method,DLM)修正其线性方法。利用无限板样条(IPS)方法进行气动结构耦合,实现了有限元结构的非线性气弹响应分析。算例结果验证了方法对于静气弹分析的有效性,同时能准确地反映弹性带来的气动效率的降低和非线性力矩特征。     

压气机转子叶片颤振边界的预测方法

通过发展的基于能量法的颤振数值预测方法得到了压气机转子叶片的颤振边界.将计算结构动力学(CSD)分析得到的叶片表面节点位移插值到耦合面的流体网格点上,并将设计的多层动网格技术应用于计算流体动力学(CFD)方法,实现叶栅振荡作用下的非定常分析,得到叶片表面的非定常气动功以及模态气动阻尼比.以某第一级压气机转子叶片为例,对...     

篦齿封严结构气弹稳定性数值分析

针对航空发动机中篦齿封严装置的气弹稳定性问题,用所发展的气弹稳定性分析的能量法,考查了两种不同的篦齿封严结构模型的颤振特性,对比分析了在不同节径振动下的气动功分布和模态气动阻尼比。结果表明:能量法可以用于预测篦齿封严结构的气弹稳定性,通过数值计算给出了改进前后篦齿封严结构低阶振动模态的气动功和模态气动阻尼比,改进后的篦齿结构在2节径振型下不再发生颤振,说明改进后的结构比改进前的气弹稳定性更好。