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

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

考虑气动弹性影响的机翼气动外形设计研究

采用三维Euler方程为控制方程,耦合静气动弹性平衡方程,进行机翼静气动弹性数值模拟;在机翼静气动弹性分析的基础上,结合Takanashi余量修正方法对三维大展弦比机翼进行气动外形反设计,以确定机翼的型架外形.以某型支线飞机的大展弦比机翼为算例,进行了静气动弹性数值模拟和机翼型架外形设计研究,设计结果表明发展的机翼静气动弹性数值模拟和型架外形设计方法是合理可行的.     

一种三维弹性机翼气动优化新方法

基于控制理论方法和有限元分析软件,通过将网格生成、流场计算、有限元分析、粘性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,研究与发展了计及静气动弹性影响的三维机翼气动优化设计方法,其中,雷诺平均Navier-Stokes方程为主控方程,气动载荷和结构静弹性变形量由气动/结构方程的耦合迭代求解得到,目标函数梯度信息由共轭方程数值求解得到。典型大展弦比机翼气动减阻设计结果表明：研究及发展的计及静气动弹性影响的三维机翼气动优化设计方法是有效的,能够有效考虑静气动弹性的影响。     

基于CFD/CSD方法的跨声速静气动弹性数值模拟应用研究

针对大型飞机跨声速静气动弹性问题,基于CFD/CSD流固耦合方法开展了数值模拟应用研究。具体方法是:采用RANS方程与静气动弹性平衡方程作为计算控制方程,CFD采用多块对接结构化网格进行分区并行计算,并利用多重网格技术加速收敛;使用RBF结合TFI方法进行网格变形;通过TPS插值实现多场数据交换。基于上述方法开发了计算程序,并与典型风洞试验结果开展了对比研究,验证了方法和程序的有效性。最后通过数值模拟结果,分析了跨声速时静气动弹性对典型大型飞机机翼的几何变形、表面压力及气动性能的影响特性。     

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

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

考虑气动弹性的风力机叶片性能分析

考虑气动弹性对风力机叶片的影响,采用叶素-动量理论计算气动力,采用盒形梁理论计算结构变形,耦合静气动弹性平衡方程,建立了风力机叶片静气动弹性分析程序。本文运用该程序进行了多种风速下叶片载荷及风轮性能的计算,分析了气动弹性对原设计的影响。结果表明,对于兆瓦级风力机,在大风速情况下,气动弹性对风轮性能有着明显影响,并会造成气动载荷的重新分布,影响结构设计的准确性。该方法可用于对叶片气动设计与载荷计算方法进行气动弹性修正。     

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

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

翼梢装置对某机翼气动弹性行为影响研究

以大客某方案机翼为基本翼,通过数值模拟的方法研究了翼梢装置对机翼气动弹性特性影响,包括静气动弹性及颤振特性。其中通过CFD/CSD弱耦合求解的方法研究其静气动弹性响应,气动力计算采用面元法,结构响应计算采用结构有限元法,通过插值实现翼面气动力与有限元节点力之间的传递,以及有限元模型与气动网格之间的变形传递。对基本翼及带翼梢装置机翼静力学有限元模型局部修改得到动力学模型,应用MSC NASTRAN进行颤振特性分析。研究发现翼梢装置使得机翼的气动弹性特性不同程度均有降低,而不同翼梢装置对其影响又有所不同,可见,翼梢装置的设计在追求气动特性改善的同时必须关注其带来的结构特性的损失。     

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

考虑静气动弹性影响的客机机翼气动/结构一体化设计研究

针对跨声速客机气动/结构一体化设计问题,建立了考虑静气动弹性影响的气动/结构一体化优化设计方法,并针对现代跨声速民用客机开展了气动/结构一体化设计研究。数值评估选择全速势方程加附面层修正,气弹分析采用基于RBF插值技术的松耦合分析方法,优化方法使用改进的微分进化算法。通过对CRM和DLR-F6标模进行计算并与实验数据对比,验证了采用的气动数值评估手段和静气动弹性分析方法可靠性。利用建立的优化设计方法对跨声速客机机翼进行了分别以扭转角分布和剖面翼型为设计变量的考虑静气动弹性影响的气动/结构一体化设计,航程分别提高了5.63%和3.05%。航程的提高主要得益于机翼的载荷分布和结构厚度分布的改变,以扭转角分布为设计变量的优化设计以2.56%的结构重量损失获得了6.53%的升阻比的提高,以剖面翼型外形为设计变量的优化重量减小了3.56%同时升阻比提高了1.53%。     

一种基于N—S方程的CFD／CSD耦合计算方法

对于气动弹性问题,在跨音速和大变形情况下,模态叠加法失效,为此需要探讨全场求解耦合方程的新方法。利用耦合计算流体动力学和计算结构动力学（CFD／CSD）方法,数值模拟了机翼的气弹响应问题。气动力的求解采用三维Navier—Stokes方程的有限体积法,实现了一种多块结构网格的动网格方法生成气动网格,采用有限元方法计算应...     

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

大攻角机翼的气动弹性计算方法

 本文利用大攻角带分离流机翼的完全非定常非线性气动力计算方法,通过与机翼运动方程的同时求解,在时间域内实现了大攻角机翼非线性气动弹性的数值模拟,根据弹性机翼各种状态下的运动过程,可以得到大攻角机翼的颤振速度等重要参数以及亚临界、超临界等飞行状态的运动规律。算例结果表明,大攻角下机翼的气动弹性问题需引起足够的重视。     

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

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

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

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

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

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

跨声速风扇叶片的静态气动弹性问题

使用时域的流固耦合数值计算方法,研究了跨声速风扇叶片在气动力和离心力共同作用下的静态气动弹性问题,分析了叶片在不同工况下的变形规律及叶片变形对整体气动性能的影响.NASA rotor 67的静态气动弹性计算说明气动力对叶片最大变形的贡献达13.07%, 而且叶片变形明显地改变了通道激波的位置和强度.宽弦空心跨声速风扇叶片的静态气动弹性计算说明叶片变形对总体气动效率的影响为0.15%~ 0.5%,其中气动力对变形贡献在叶片尖部的前缘可达41%,考虑气动力引起的变形使得该风扇的流量增大,气动特性线整体向右偏移.计算结果说明:气动力的非线性对跨声速风扇叶片静态变形问题有显著的影响,工程实践中从设计叶型到制造叶型的反扭过程应该采用流固耦合方法以得到更准确的叶型.      

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

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