三维机翼跨音速颤振的分布式并行计算

以多台微机组成分布式计算系统,利用无限插值方法(TFI)生成三维多块贴体运动网格,以Navier Stokes方程为控制方程,求解三维机翼的跨音速非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义位移响应的时间历程,根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定机翼跨音速颤振临界条件.通过算例验证,计算结果与实验结果和理论分析相吻合.     

钝前缘梯形翼高超声速风洞颤振试验

为了研究钝前缘翼面的高超声速颤振特性,获得典型翼面高超声速颤振参数以校验非定常气动力和CFD计算,采用具有简单结构动力学特性的钝前缘梯形翼模型,在中国航天空气动力技术研究院FD-07高超声速风洞进行了高超声速风洞颤振试验研究.模型为9 mm厚钝前缘梯形平板翼,采用夹层设计:中间层为钢板,提供模型主要刚度和质量特性;两侧为泡沫,起维形作用.试验模型采用悬臂支撑安装于风洞试验段,试验Mach数分别为4.95和5.95.试验固定Mach数,通过缓慢增加动压以使模型达到颤振临界点,采用小波时频谱分析时域响应,结果显示试验模型发生了弯扭耦合经典颤振.试验采用直接观测法获得了颤振动压、颤振频率和对应的试验密度、总温等颤振相关参数.采用壳单元建立了结构有限元模型,并采用统一升力面理论对模型进行了颤振计算分析,研究了气流密度、结构阻尼、Mach数对颤振计算的影响,并对试验结果与理论计算的偏差进行了讨论.分析认为,计算气流密度、计算结构阻尼、结构建模偏差、试验结果散布特性等因素均会构成计算值和试验值之间的偏差,但即便在计算中考虑上述因素,计算结果与试验值仍存在较大偏差.      

一种基于连续核函数展开的跨声速气动力修正法

现有的跨声速修正面元目的 大多难以准确修正非定常压强系数的虚部.结合下洗修正目的 和连续核函数展开法(Successive Kernel Expansion Method),发展了一种新的修正目的.计算ONERA M6机翼的非定常气动力,与传统修正目的 比较,结果表明该目的 具有更高的精度.将改进的非定常气动力用于机翼的跨声速颤振分析,较好地预测了机翼在跨声速区域的非线性颤振特性.     

宽弦高速跨音风扇颤振特性研究

颤振是航空发动机、燃气轮机等运行安全的重要威胁,但颤振稳定性与流动结构之间的关系尚不清晰。本文使用行波法和影响系数法,对某宽弦复合掠型高速跨音风扇转子的一阶模态进行了颤振特性研究,计算了在100%转速下从堵塞点到近失速点的颤振表现。使用影响系数法时,分析了不同通道数的计算域对气动阻尼计算的影响,并与行波法得到的结果进行了对比。研究了流动结构与叶片表面气动阻尼之间的关系,旨在提高对流动致颤机理的认识。结果表明影响系数法和行波法均能对叶片的气动阻尼进行较好的预测;流动结构方面,激波、激波附面层分离、叶尖泄漏流以及吸力面前缘叶顶附近的非定常压力波动,对叶片的气动阻尼分布有较大的影响。     

振动叶栅叶片间振动相位差对非定常气动力的影响

一、前言 在叶轮机械非定常气动力的研究中,探索振动叶栅叶片间振动相位差对非定常气动力影响的规律,可以大大节省对叶片进行颤振预估的理论与实验方面的工作量,具有重要意义。在现已发表的文献中,虽然给出了几种相位差下某种振型的实验结果或理论分析,但只限于一种振型下的少数几个相位差,而且也没有得出定量的规律。由于以往计算非定常气动力的方法太烦杂且又受到各种限制,要进行系统分析是十分困难的。文献[3]提出的求解振动叶栅非定常气动力的“振荡中心流线法”,提供了一种适应性广、计算简单     

大展弦比气弹机翼柔性非线性变形的气动力效应分析与风洞试验

从飞行器刚弹耦合动力学模型出发，引入柔性机翼准定常假设，建立大柔性飞行器非线性静气动弹性气动力方程，利用非线性迭代求解思路模拟了柔性飞行器的静气动弹性响应行为，开展了大展弦比飞机静气动弹性风洞试验验证，采用气动力有限基本解与机翼的耦合计算，发现了大柔性飞机大变形状态下载荷及结构变形形式随风速的变化规律.传统基于小变形假设的线性分析方法和刚体分析由于无法考虑气动面随结构变形的曲面气动力因素和结构变形后的非线性刚度特性，均与风洞试验存在一定的误差.对于大展弦比柔性飞机的非线性静气动弹性分析十分必要.      

振荡机翼跨音速非定常粘性绕流的数值计算

以非定常Ｎ－Ｓ方程为主管方程，采用ＬＵ－ＮＮＤ混合差分格式，Ｃ和Ｃ－Ｈ型贴体运动网格，Ｂ－Ｌ双层代数紊流模型，求解绕振荡翼型和三维机翼的跨音速非定常粘性流场，分别计算了ＮＡＣＡ００１２翼型和Ｍ６机翼作俯仰振荡时跨音速非定常粘性绕流流场。研究了非定常绕流的气动特性，部分计算结果和风洞实验值作了比较。     

利用N-S方程模拟机翼气动弹性的一种计算方法

利用一种双时间方法求解三维非定常N-S方程,得到与任意非定常运动对应的气动力,在求解非定常气动力的同时,在时间域内用二阶龙格 库塔方法求解机翼弹性运动方程,从而模拟粘性流动中的气动弹性全过程.为保证网格生成效率,采用无限插值理论生成O-H型代数网格,考虑了机翼变形时的网格生成问题,并得到计算结果.     

叶轮机械叶片流体激振安全性的全功能分析

本文提出了对叶轮机械叶片流体激振安全性实行全功能数值分析的新概念和方法。实现气动力学和弹性力学的弱耦合计算,由弹性力学方法确定三维振型,用振荡流体力学方法确定三维条件下的三维振荡流场,确定在各种振幅条件下的非定常气动力作功与阻尼功,引入气动功等于阻尼功条件下的平衡振幅概念。根据平衡振幅确定叶片流体激振的动应力,并进行寿命评估。这种全功能分析将使叶轮机械流体激振的安全性分析更加符合工程实际的需要,也能够解释大量的流体激振破坏不属于瞬时突然破坏,而是一种疲劳破坏。     

飞行器大攻角非定常气动特性神经网络建模

大攻角气动特性预测与气动建模是新型飞行器提升飞行性能的重要内容.以轴对称导弹简化模型为研究对象，首先采用计算流体力学方法，对70°大攻角状态的非定常气动特性进行数值模拟，计算方法基于RANS的N-S方程，湍流模型采用SA模型，对流场采用有限体积法离散，无黏项采用Roe通量差分分裂格式，黏性项采用中心差分，时间推进采用LU-SGS格式的双时间步法.飞行器运动模式采用强迫振荡的方式，对5种不同振荡频率进行了非定常数值计算，并记录每一内迭代周期最终的气动力和力矩数值.其次，以CFD预测结果作为气动建模的样本，采用动导数模型、多项式模型等传统方法，进行气动建模，并分析其有效性和精度.最后采用神经网络方法对大攻角非定常气动力进行建模，并和动导数模型、多项式模型进行精度对比.结果表明，基于神经网络的人工智能气动建模方法具有较高的精度和适应性.该方法为飞行器大攻角非定常非线性气动建模，大攻角飞行稳定性分析与控制提供理论参考.      

Uncertainty analysis near bifurcation of an aeroelastic system

Variations in structural and aerodynamic nonlinearities on the dynamic behavior of an aeroelastic system are investigated. The aeroelastic system consists of a rigid airfoil that is supported by nonlinear springs in the pitch and plunge directions and subjected to nonlinear aerodynamic loads. We follow two approaches to determine the effects of variations in the linear and nonlinear plunge and pitch stiffness coefficients of this aeroelastic system on its stability near the bifurcation. The first approach is based on implementation of intrusive polynomial chaos expansion (PCE) on the governing equations, yielding a set of nonlinear coupled ordinary differential equations that are numerically solved. The results show that this approach is capable of determining sensitivity of the flutter speed to variations in the linear pitch stiffness coefficient. On the other hand, it fails to predict changes in the type of the instability associated with randomness in the cubic stiffness coefficient. In the second approach, the normal form is used to investigate the flutter (Hopf bifurcation) boundary that occurs as the freestream velocity is increased and to analytically predict the amplitude and frequency of the ensuing LCO. The results show that this mathematical approach provides detailed aspects of the effects of the different system nonlinearities on its dynamic behavior. Furthermore, this approach could be effectively used to perform sensitivity analysis of the system's response to variations in its parameters.     

超临界翼型的跨音速抖振特性

以二维非定常N-S方程为基本方程,计算跨音速翼型升力系数的时间历程.根据升力系数的脉动量急剧上升的起始点确定抖振起始边界,以超临界机翼DFVLR-R2和传统翼型NACA0012为研究对象,研究了两种翼型的抖振特性,计算结果表明,在超临界翼型的设计马赫数附近,超临界翼型具有良好的抖振特性.     

LIMIT CYCLE FLUTTER OF AIRFOILS IN STEADY AND UNSTEADY FLOWS

The limit cycle flutter of a two-dimensional wing with non-linear pitching stiffness is investigated. For modelling the aerodynamic forces of the wing steady linear and non-linear models as well as an unsteady model were used. The flutter speed was calculated using the harmonic balance method and by predicting Hopf bifurcation. Analytical solutions based on the centre manifold theory and normal forms were obtained as were results given by the harmonic balance method. The analytical solutions were compared with those obtained by numerical integration. The results show that the harmonic balance method can forecast flutter speed with a good accuracy while analytical solutions based on centre manifold theorem are accurate only in a small neighbourhood of the bifurcation point. The oscillation of the airfoil after flutter for two different models, linear and non-linear pitching stiffness were compared with each other and the flutter speeds for two linear steady and an unsteady aerodynamic model calculated. The obtained results show that flutter analysis based on the linear steady model is conservative only for the ratios of plunge frequency to pitch frequency lower than 1.     

Numerical study on the correlation of transonic single-degree-of-freedom flutter and buffet

Transonic single-degree-of-freedom(SDOF) flutter and transonic buffet are the typical and complex aeroelastic phenomena in the transonic flow. In this study, transonic aeroelastic issues of an elastic airfoil are investigated using Unsteady Reynolds-Averaged Navier-Stokes(URANS) equations. The airfoil is free to vibrate in SDOF of pitching. It is found that, the coupling system may be unstable and SDOF self-excited pitching oscillations occur in pre-buffet flow condition, where the free-stream angle of attack(AOA) is lower than the buffet onset of a stationary airfoil. In the theory of classical aeroelasticity, this unstable phenomenon is defined as flutter. However, this transonic SDOF flutter is closely related to transonic buffet(unstable aerodynamic models) due to the following reasons. Firstly, the SDOF flutter occurs only when the free-stream AOA of the spring suspended airfoil is slightly lower than that of buffet onset, and the ratio of the structural characteristic frequency to the buffet frequency is within a limited range. Secondly, the response characteristics show a high correlation between the SDOF flutter and buffet. A similar "lock-in" phenomenon exists, when the coupling frequency follows the structural characteristic frequency. Finally, there is no sudden change of the response characteristics in the vicinity of buffet onset, that is, the curve of response amplitude with the free-stream AOA is nearly smooth. Therefore, transonic SDOF flutter is often interwoven with transonic buffet and shows some complex characteristics of response, which is different from the traditional flutter.     

Bifurcation and chaos analysis for aeroelastic airfoil with freeplay structural nonlinearity in pitch

The dynamics character of a two degree-of-freedom aeroelastic airfoil with combined freeplay and cubic stiffness nonlinearities in pitch submitted to supersonic and hypersonic flow has been gaining significant attention. The Poincaré mapping method and Floquet theory are adopted to analyse the limit cycle oscillation flutter and chaotic motion of this system. The result shows that the limit cycle oscillation flutter can be accurately predicted by the Floquet multiplier. The phase trajectories of both the pitch and plunge motion are obtained and the results show that the plunge motion is much more complex than the pitch motion. It is also proved that initial conditions have important influences on the dynamics character of the airfoil system. In a certain range of airspeed and with the same system parameters, the stable limit cycle oscillation, chaotic and multi-periodic motions can be detected under different initial conditions. The figure of the Poincaré section also approves the previous conclusion.     

基于进化算法和复合参数化方法的低速翼型优化

针对设计工况为Re=1×106的翼型在2°和8°两种迎角下的气动性能进行综合优化, 提升其在巡航与起飞着陆时的性能。利用气动估算软件Xfoil结合NACA 4位数参数化方法进行翼型初步优化, 缩小优化空间, 减少计算量。通过求解Reynolds平均Navier-Stokes方程的方式进行精确优化, 最后利用Hicks-Henne形函数法对翼型表面施加扰动进行细节优化, 弥补NACA 4位数参数化方法细节表现不足的缺点。实现了2°和8°两种迎角下气动性能普遍提升20%以上的效果, 并发展了一种复合参数化方法下基于改进进化算法的多目标、多步骤气动优化方法。      

Modeling and identification of freeplay nonlinearity

A nonlinear analysis is performed for the purpose of identification of the pitch freeplay nonlinearity and its effect on the type of bifurcation of a two degree-of-freedom aeroelastic system. The databases for the identification are generated from experimental investigations of a pitch-plunge rigid airfoil supported by a nonlinear torsional spring. Experimental data and linear analysis are performed to validate the parameters of the linearized equations. Based on the periodic responses of the experimental data which included the flutter frequency and its third harmonics, the freeplay nonlinearity is approximated by a polynomial expansion up to the third order. This representation allows us to use the normal form of the Hopf bifurcation to characterize the type of instability. Based on numerical integrations, the coefficients of the polynomial expansion representing the freeplay nonlinearity are identified.     

气动/几何约束条件下翼型优化设计的最优控制理论方法

基于最优控制理论原理和Navier-Stokes方程,研究了气动/几何约束条件下多设计变量的翼型气动优化设计问题.根据给定的目标函数表达形式,在计算坐标下详细推导了相应的共轭方程及边界条件,以及梯度方程的具体数学形式.通过合理数学变换,得到了物理空间上适应于CFD数值求解的共轭方程直观表达形式,并发展了有效数值求解目的.通过将流动方程、共轭方程、目标函数敏感性导数和优化算法相结合,发展了一种新的气动优化设计目的.相关设计算例表明该目的在设计理论、适用性以及时间费用等方面具有显著特色和优点,且设计结果更为可靠.     

Aeroelastic dynamic response and control of an airfoil section with control surface nonlinearities

Nonlinearities in aircraft mechanisms are inevitable, especially in the control system. It is necessary to investigate the effects of them on the dynamic response and control performance of aeroelastic system. In this paper, based on the state-dependent Riccati equation method, a state feedback suboptimal control law is derived for aeroelastic response and flutter suppression of a three degree-of-freedom typical airfoil section. With the control law designed, nonlinear effects of freeplay in the control surface and time delay between the control input and actuator are investigated by numerical approach. A cubic nonlinearity in pitch degree is adopted to prevent the aeroelastic responses from divergence when the flow velocity exceeds the critical flutter speed. For the system with a freeplay, the responses of both open- and closed-loop systems are determined with Runge-Kutta algorithm in conjunction with Henon’s method. This method is used to locate the switching points accurately and efficiently as the system moves from one subdomain into another. The simulation results show that the freeplay leads to a forward phase response and a slight increase of flutter speed of the closed-loop system. The effect of freeplay on the aeroelastic response decreases as the flow velocity increases. The time delay between the control input and actuator may impair control performance and cause high-frequency motion and quasi-periodic vibration.