风洞实验中机翼振动对气动特性影响的计算分析

通过求解非定常N-S方程研究了弹性机翼在低速气流中振动时的气动特性,并与刚性机翼气动特性做了比较。通过结构有限元方法计算出模型机翼的结构模态和相应振动频率,假设机翼以最小刚度的弯曲模态振动,计算了不同振幅下不同迎角、尤其是失速迎角附近的升力特性。计算结果表明,在小振幅下机翼的平均升力特性和刚性机翼相差很小,但当振幅达到一定幅度后,振动机翼下平均升力特性与刚性机翼的结果有较明显的差别,主要表现在最大升力系数的降低和失速迎角的减小。此结果可以解释实际风洞实验中最大升力系数往往小于CFD结果、实验中最大升力系数对应的迎角也比CFD结果小的现象。     

机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

考虑动力影响的民用飞机静气动弹性分析方法

采用了一种基于多块网格的N-S方程和结构柔度影响系数法,考虑气动、结构非线性的基于RBF插值和RBFDelaunay动网格变形技术的静气动弹性分析方法对喷流对弹性机翼的气动力影响进行了研究。利用DLR F6翼身组合体构型对静气动弹性方法进行验证,保证了计算的可信性。采用该方法对比分析了某民用飞机无喷流/有喷流构型的静气动弹性特性,表明发动机喷流会给机翼带来一个正的扭转效应,抵消一部分机翼后掠效应的影响,使机翼前后缘挠度均会有所增大,弹性变形引起的多数剖面的附加扭转角有所减小。研究表明:喷流影响会使刚性机翼表面的压力分布发生变化,升力系数有所损失;考虑喷流的机翼静气动弹性变形是一个耦合效应,发动机喷流区主要受喷流影响,外翼段主要受弹性变形影响。数值模拟结果表明:无喷流影响时机翼的弹性变形使升力系数下降约16%,升阻比下降8.4%,考虑喷流影响时,升力系数下降达到18%,升阻比下降36%。因此,对于大展弦比机翼,考虑喷流影响的静气动弹性分析十分必要。     

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

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

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

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

基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

基于CFD的机翼突风响应计算

利用网格速度理论,计算机翼在锐边突风和1-cos突风下的响应,研究气动非线性和自由度耦合对翼尖加速度和升力系数的影响.采用中心格式有限体积法进行空间离散,并用双时间推进法求解非定常Euler方程.计算了刚性(沉浮)和弹性机翼在锐边突风下的加速度响应,以及刚性机翼(沉浮+俯仰)在1-cos突风下的升力响应过程,并分别与片...     

大展弦比柔性机翼的结构与气动耦合特性分析

针对人展弦比柔性机翼，本文提出了一种考虑结构与气动相互耦合的静气动弹性分析方法，通过以下迭代过程得到结构变形和气动载荷的稳定状态：采用亚音速面元法年和核函数法计算气动影响系数和翼面升力：采用有限元技术计算结构变形：采用B样条函数实现结构网格变形到气动网格变形的插值；采用载荷等效的二角形方法进行气动载筒和有限元节点载荷的转换。最后，通过一个后掠翼算例计算了柔性下的静气动弹性响应，并与刚性情况进行了对比。     

耦合六自由度运动弹性飞机阵风响应数值模拟

耦合飞机刚体六自由度运动,考虑飞机结构气动弹性变形,基于全湍流Navier—Stokes方程流场数值模拟方法．建立了离散阵风作用下弹性飞行器气动载荷与飞行姿态响应的数值模拟技术。其中模态空间中结构动力学方程以及六自由度运动方程分别采用四阶R—K进行时间推进求解,采用RadialBasicFunction（RBF）技术进行气动／结构数据耦合。对应飞机姿态以及弹性变形采用RBF与TFI组合模式的动网格方法进行网格更新;以飞翼布局弹性飞机为数值研究对象,研究了离散阵风作用下飞行器的气动载荷响应以及飞行姿态的变化,对比分析了刚性飞行器与弹性飞行器对离散阵风响应的区别,为进一步系统地分析真实飞机飞行品质提供数值平台。     

高精度的机翼静气动弹性数值模拟方法研究

气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。     

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

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

基于响应面法的机翼气动/结构一体化优化设计研究

基于响应面法进行了机翼气动/结构一体化优化设计研究,流动控制方程为三维欧拉方程,采用有限体积法进行数值求解,应力和结构变形采用有限元方法计算,静气动弹性分析采用强耦合迭代方式,响应面模型采用二次多项式来构造。以跨音速M6机翼为初始机翼,进行了多目标、多约束情形下的气动/结构综合优化设计,优化后所得到的新机翼具有更佳的气动/结构综合性能,升阻比增加了9.25%,而重量减轻了4.84%;响应面模型精度满足设计要求,拟合误差均不超过1%;这说明本文所发展气动/结构综合优化设计方法是成功且有效的,具有广泛的应用前景。     

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

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

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

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

一种基于有限元法和弹性接触理论的齿轮啮合刚度改进算法

提出了一种将有限元法和弹性接触理论相结合的齿轮啮合刚度计算方法.该方法利用子结构法提取齿面原始柔度矩阵并分离出接触点弯曲变形,根据线弹接触变形解析公式计算接触变形,通过求解非线性变形协调方程得到齿轮时变啮合刚度和齿面载荷分布.以一对齿轮副为例,计算的啮合刚度与航空标准计算结果相差在6%以内.该方法发挥了有限元法在预测物体整体变形方面的优势,同时结合弹性接触理论能够准确计算局部接触变形的优点,与常规有限元法相比,能够有效地提升计算效率.由于接触变形问题的非线性,啮合刚度随总啮合力增加呈现非线性增大的趋势.      

Coupling fluid (CFD) and structural (FE) models using finite interpolation elements

In the field of aeroelasticity, interactions between elastic structures and fluid flow are investigated. Recently, numerical aeroelastic models have been built composing those of the combining fluid dynamics (CFD) and the computational structural dynamics (CSD) domains. Since the fluid and the structural models differ in their formulation and discretization, an interface model has to be introduced that represents the connectivity and physical interaction between the two single domain models. In the following, a scheme for coupling fluid (CFD) and structural models (FE) in space is presented which is based on finite interpolation elements. It is applied to static aeroelastic problems, in order to predict the equilibrium of elastic wing models in transonic fluid flow. The structure is represented by finite elements and the related equations are solved using commercial FE analysis codes. The transonic fluid flow is described by the three-dimensional Euler equations, solved by an upwind scheme procedure. The resulting coupled field problem containing the fluid and the structural state equations, is solved by applying a partitioned solution procedure. In each solution step the interface and boundary conditions are exchanged and updated. Here, a fixed-point iteration scheme is used for the coupled aeroelastic equations.