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

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

一种减小再入飞行器侧向气动非线性的布局优化方法

为改善升力式再入飞行器在跨声速段出现的侧向气动特性非线性问题,发展了一种基于Kriging代理模型的自适应迭代气动布局优化方法.设计了一种常规升力式再入飞行器布局,计算了该布局在跨声速段的侧向气动力,分析了可能影响侧向气动特性的机翼布局参数.根据气动布局优化流程,计算了气动布局样本气动特性,建立了布局参数到侧向力矩系数...     

基于Kriging自适应代理模型的气动优化方法

气动优化设计中,引入代理模型可以有效减少计算周期,而运用有效的插值和选样方法(自适应选样)可以大大减少建立代理模型的时间,因此提出了一种基于Kriging自适应代理模型的气动优化方法。使用Kriging方法建立代理模型,通过求解EI函数最大值得到添加样本点更新代理模型,提高了代理模型的拟合精度。针对Kriging自适应代理模型的精确性和有效性,分别进行典型函数测试分析和翼型算例验证。结果表明:基于Kriging自适应代理模型气动优化方法可以实现高效的翼型气动性能优化设计。     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

基于FLUENT的双翼末敏弹气动特性研究

为研究双翼末敏弹气动特性,基于计算流体力学方法,对S-C形尾翼组合的末敏弹气动外形流场进行计算.获得了末敏弹气动外形表面压力分布规律,以及阻力系数、升力系数、转动力矩系数和压力中心系数随迎角的变化规律.通过模型自由飞行试验获得了末敏弹动态阻力系数和转动力矩系数,结合末敏弹飞行姿态变化分析了双翼末敏弹动态气动特性.结果表明双翼末敏弹增阻效果明显,能够提供平衡弹体自身重量的阻力,达到减速要求.S-C形尾翼组合结构能够为末敏弹提供转动力矩,使末敏弹产生具有稳态扫描形式的运动.研究结果可以为改进末敏弹稳态扫描平台设计提供参考.     

基于 Kriging模型机翼平面外形气动优化设计

使用基于Kriging模型的优化设计方法,进行了非常规布局机翼的平面外形多目标优化设计。利用CFD技术进行机翼升力系数和阻力系数的气动计算,通过拉丁超立方试验设计生成样本点,建立了Kriging代理模型,结合多目标遗传算法对机翼平面外形进行多点多目标优化设计,最终得到了Pareto最优解集。根据设计需求,从Pa-ret0前沿选取一个非劣解作为优化结果。结果表明：陆ging模型与cFD计算误差很小,可信度高;在不问设计状态下,机翼气动性能都得到了提高,表明优化设计方法具有可行性和高效性。     

直升机飞行力学线性变参数建模方法研究

针对直升机飞行力学建模问题,分别基于模型缝合及模糊逼近提出了两种线性变参数(LPV)建模方法。在模型缝合方法中,基于经过滤波加权的飞行速度进行气动导数和操纵导数的插值,完成气动力的实时计算,结合非线性欧拉运动方程完成LPV模型的建立。在模型逼近方法中,建立了基于T-S模糊模型的非线性状态空间逼近模型,采用anfis结构完成了模糊模型训练。利用UH-60直升机非线性模型生成了不同飞行速度的线性状态空间模型,并分别通过两种方法建立了两套LPV模型。最后,基于频率响应和泰尔不等系数对典型LPV建模参数对模型精度的影响特性进行了分析。     

基于代理模型的气动外形平面参数多目标匹配设计

将Kriging代理模型和Pareto遗传算法引入气动外形平面参数匹配设计中,提出一种基于代理模型的多目标平面参数匹配设计方法。将拉丁超立方试验设计用于平面参数筛选,确定出参数匹配方案库;基于方案库的计算流体力学(CFD)分析结果构建Kriging气动代理模型;将Kriging模型替代CFD分析,用于气体布局参数匹配优化设计,提高了设计效率并保证了可信度;通过Pareto遗传算法优化解决多点设计要求下气动布局参数匹配问题,一次性给出参数匹配方案的最优解集,从Pareto前沿中根据设计偏向选择气动布局最佳匹配方案。以典型的双后掠布局平面参数多点匹配优化设计问题作为算例,研究结果表明:Kriging气动代理模型与实际CFD分析结果的误差均小于5%,满足精度要求;根据不同设计偏向,选择了3种参数匹配Pareto优化方案,与原样本方案相比超声速阻力减小6.0%~12.8%,跨声速升阻比增加0.01%~3.40%,证明了匹配设计方法的有效性;通过试验设计的Pareto分析与主、交互效应分析,获得了气动布局平面参数对气动性能影响的定量关系,能够为参数匹配设计提供依据。所提出的平面参数匹配设计方法可应用于其他常规与非常规气动布局型式。     

大展弦比前掠翼气动隐身多目标优化

前掠翼布局优越的气动性能为无人机气动布局设计提供了一条新的方向。采用CST方法对翼型几何外形进行参数化描述,实现前掠翼气动和隐身多学科优化设计模型的参数化描述。建立了基于N-S方程的计算流体力学方法的前掠翼气动分析模型和基于矩量法的计算电磁学方法的前掠翼隐身分析模型。提出了基于Kriging模型的前掠翼气动隐身多目标优化方法,采用拉丁超方试验设计方法获取样本点,建立前掠翼气动和隐身的Kriging代理模型。将Pareto多目标遗传算法与Kriging代理模型结合进行大展弦比前掠翼的气动隐身多目标优化设计。研究结果表明,所建立的分析模型是合理的,所提出的多目标优化设计方法是可行的,能够有效提高大展弦比前掠翼性能与优化效率。     

飞翼无人机进舰下滑纵向固有模态特性研究

针对舰载飞翼布局无人机与常规布局飞行器不同的气动外形,建立无人机飞行状态下附近流场的网格模型,并计算了该型无人机的气动参数。基于插值函数计算得出该型无人机的升力系数、阻力系数和俯仰力矩系数对于迎角变化的曲线。建立纵向小扰动方程进行仿真计算,求解出该型无人机在进舰下滑状态时的固有模态特性。计算结果表明,该模型在配平条件下可用于舰载飞翼布局无人机的着舰气动适配性研究。     

基于非线性试验气动力的飞机静气动弹性响应分析

分别使用线性和非线性静气动弹性分析方法,对某飞机纵向静气动弹性响应特性随飞行动压、攻角、平尾偏度和纵向过载变化的趋势进行了分析,并将分析结果和飞行试验进行了比较。线性方法的气动力计算使用亚音速偶极子格网法。非线性方法由风洞试验提供刚体气动力,并使用线性气动力影响系数矩阵对其进行弹性化处理。通过对两种方法的计算结果进行对比分析可以看出:①使用刚体试验气动力的非线性分析方法能够获得和飞行试验比较一致的结果;②在线性方法所提供的结果中,翼面的弯曲变形及剪力、弯矩和扭矩等部分参数能为初步设计提供大致参考,但气动力系数的弹性增量、翼面的扭转变形、翼面的剪力、弯矩和扭矩的弹性部分可能不能为型号设计提供正确的指导;③不论刚体气动力为线性还是非线性,气动力系数随动压变化的趋势均呈线性。     

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.     

基于CFD/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

基于代理模型的运载器头罩外形优化设计

引入Kriging代理模型,选取若干不同头罩外形的运载器,进行气动性能分析,利用运载器气动性能参数作为拟合样本建立代理模型.以具有足够精度的代理模型替代CFD分析,发展了一种基于代理模型的运载器头罩外形优化设计方法.在马赫数为3、飞行攻角为3°、飞行高度为8km条件下,利用该方法对运载器进行最小阻力系数、最大纵向压心系数的单目标优化和综合考虑上述2个目标的多目标优化.结果表明:2个目标存在冲突,为同时兼顾减小阻力和增强纵向稳定性,必须对运载器头罩外形进行多目标优化,得到的外形阻力系数减小了1.95%,纵向压心系数增大了5.93%,采用基于代理模型的优化设计方法能在保证精度条件下有效提高计算效率.      

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

Integrated optimization on aerodynamics-structure coupling and flight stability of a large airplane in preliminary design

The preliminary phase is significant during the whole design process of a large airplane because of its enormous potential in enhancing the overall performance. However, classical sequential designs can hardly adapt to modern airplanes, due to their repeated iterations, long periods, and massive computational burdens. Multidisciplinary analysis and optimization demonstrates the capability to tackle such complex design issues. In this paper, an integrated optimization method for the preliminary design of a large airplane is proposed, accounting for aerodynamics, structure, and stability. Aeroelastic responses are computed by a rapid three-dimensional flight load analysis method combining the high-order panel method and the structural elasticity correction. The flow field is determined by the viscous/inviscid iteration method, and the cruise stability is evaluated by the linear small-disturbance theory. Parametric optimization is carried out using genetic algorithm to seek the minimal weight of a simplified plate-beam wing structure in the cruise trim condition subject to aeroelastic, aerodynamic, and stability constraints, and the optimal wing geometry shape, front/rear spar positions, and structural sizes are obtained simultaneously. To reduce the computational burden of the static aeroelasticity analysis in the optimization process, the Kriging method is employed to predict aerodynamic influence coefficient matrices of different aerodynamic shapes. The multidisciplinary analyses guarantee computational accuracy and efficiency, and the integrated optimization considers the coupling effect sufficiently between different disciplines to improve the overall performance, avoiding the limitations of sequential approaches utilized currently.     

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

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

低速流场中柔性悬臂板的后颤振响应

建立了一个新的非线性气动弹性模型,对低速流场中柔性悬臂板的后颤振响应特性进行了分析。建模中考虑了结构几何非线性、气动力非线性以及两者之间的强耦合效应。通过实验数据对所建立的气动弹性模型进行了验证。发现在低速流场中柔性悬臂板可能会以周期加倍的方式进入混沌运动。结构几何非线性效应和翼尖涡引起的非定常气动力效应对柔性悬臂板的结构响应有显著影响,而尾涡变形引起的非定常气动力对结构运动的影响较小。还研究了不同耦合算法的差异,给出了小展弦比大柔性结构非线性气动弹性数值仿真时耦合策略的选择依据。     

基于能量法的跨声速风扇叶片气弹稳定性研究

基于能量法原理,采用弱耦合的方法对跨声速风扇叶片进行气弹稳定性分析;采用3维线性插值算法编程实现CSD/CFD数据交换。利用有限元法计算叶片的模态振型和固有振动频率;应用FLUENT动网格技术,对NASA R 67风扇叶片在不同模态振动下的非定常流场进行数值模拟,给出在前3阶模态振动下叶片表面的周期累积气动功和气动阻尼,并探讨了激波对跨声速风扇叶片气弹稳定性的影响。