大后掠翼前缘涡对其颤振特性的影响

大迎角三角翼的前缘涡不仅可以改善其气动力特性,也会显著影响机翼的气动弹性特性.运用基于Euler方程的非定常气动力降阶模型(ROM)方法,耦合结构运动方程,在状态空间内建立了气动弹性分析模型,研究了70°削尖三角翼的大迎角颤振特性.研究结果显示前缘涡对该机翼颤振特性的影响不可忽略.颤振速度随迎角的增加而大幅降低,迎角α=20°时的颤振速度比α=0°时降低了22%.发现了颤振特性随迎角变化时出现的不连续现象,并揭示了该现象是由于系统颤振分支随着静态迎角的增加发生转移所致.     

大后掠翼前缘涡对其颤振特性的影响

大迎角三角翼的前缘涡不仅可以改善其气动力特性,也会显著影响机翼的气动弹性特性。运用基于Euler方程的非定常气动力降阶模型（ROM）方法,耦合结构运动方程,在状态空间内建立了气动弹性分析模型,研究了70°削尖三角翼的大迎角颤振特性。研究结果显示前缘涡对该机翼颤振特性的影响不可忽略。颤振速度随迎角的增加而大幅降低,迎角α=20°时的颤振速度比α=0°时降低了22％。发现了颤振特性随迎角变化时出现的不连续现象,并揭示了该现象是由于系统颤振分支随着静态迎角的增加发生转移所致。     

涡破裂诱导的垂尾抖振数值模拟

 采用数值模拟方法研究了大后掠三角翼前缘涡破裂诱导的垂尾抖振问题,分析了大迎角条件下的垂尾抖振特性。采用Navier-Stokes方程求解非定常气动力、耦合结构动力学方程,建立了气动弹性方程,在时域内采用松耦合方式推进以得到垂尾结构响应。研究结果表明:涡破裂流的脉动频带覆盖了垂尾扭转模态的固有频率,诱发了垂尾抖振现象;与传统的颤振频域响应特性不同,垂尾抖振响应的各阶位移与加速度响应主频均位于各阶结构模态固有频率附近。此外,弯曲与扭转响应存在耦合效应,且耦合作用的频率与提取的垂尾表面气动载荷脉动频率一致。垂尾的位移响应由一阶弯曲模态主导,振幅不大;加速度响应主要由扭转模态产生,量级较大,使结构持续遭受严重的附加惯性载荷作用。     

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

边条翼前缘涡非定常涡场特性研究中PIV技术的应用

描述了应用PIV技术在水槽中对边条机翼上旋涡及破裂旋涡流场进行的测量和分析。实验是在北航水槽中进行的。通过PIV技术的测量，揭示了旋涡及破裂旋涡中的非定常特性，这种非定常特性同飞机上机翼、尾翼的抖振密切相关。实验结果表明，对于未破裂的边条涡，存在着两种非定常特性，其一是剪切层中不断地有小涡沿剪切层输运和合并。其二是由一次涡诱导的二次涡与剪切层中的小涡互相诱导引起的非定常现象。对于破裂涡，则发现与未破裂的涡相比，截面上涡量分布的区域突然扩大很多，最大涡量的绝对值也比上游未破裂区截面上的涡量最大值小。此外还发现在涡量分布区域出现反涡量，这同涡破裂后出现涡核螺旋变形有关。对于同一截面处涡量分布是非定常的。     

三角翼跨声速动态失速与涡破裂特性研究

通过数值方法研究了高速气流中细长三角翼作匀速上仰机动飞行时的背风面分离流与涡破裂特性和气动力特性。结果表明,在三角翼匀速上仰非定常绕流中,涡破裂起始攻角与上仰速度关系比较复杂,在本计算速度范围内,涡破裂起始攻角变化很小;动态时涡破裂点发展速度随攻角变化规律与静态结果差别甚大;三角翼大攻角非定常绕流前缘涡涡轴附近的流动以不稳定涡和涡破裂为主要特征;在一定的上仰速度以后即使在中等攻角前非定常升力与定常升力仍有较大差别,与低速情况有所不同;上仰运动使涡破裂点发展速度被延缓,从而提高了失速攻角和最大升力,并使失速攻角与涡破裂起始攻角的间隔进一步拉大,同时不改变升阻比     

涡襟翼振动对三角翼涡的影响

 <正> 1.引言 许多实验表明,用主动干扰的方式是控制分离的一个有效方法,而且可以把分离区从一死水区(或紊乱的区域)变成一个有序流动区域。二维流动实验还表明非定常机动可以较大地改变流动结构。三维流态显示表明强迫振动对集中涡的形成过程等有显著影响,并影响涡和物面之间的距离。目前对三维非定常流动特性,以及非定常干扰对机翼绕流中涡破裂等的作用尚未充分了解。 影响集中涡破裂的关键因素是沿涡轴方向的压力梯度。涡环量对破裂的影响是双向的,涡强过大易使涡破裂;涡强过小时涡结构松散,也易于破裂,甚至迅速耗散掉。从二维结果看,非定常强迫扰动可使涡的结构更紧凑、清晰,而且扰动可以改     

三角翼匀速上仰跨声速非定常流场数值模拟

用有限差分法求解Ｎ－Ｓ方程，数值模拟了细长三角翼作匀速上仰运动时的距声速绕流流场。给出了非定常流场结构，涡破裂特性和非定常升力，阻力特性，研究了上仰运动对动态失速的影响。     

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

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

三角翼涡破裂非定常特性实验研究

本文依据染色液流动显示结果,通过子波和频谱分析,探讨了70°三角翼前缘涡涡核轴向速度的变化规律及其子波特性、涡破裂位置的非定常特性,指出涡破裂点位置的变化属于低频高幅振荡,这主要是左右涡之间的相互作用造成的,当两个涡的时间平均涡破裂点位置彼此靠得更近时,相应的振荡就更大一些,此外本实验得到涡破裂位置振荡的折合频率在St=0.2以内.     

涡破裂诱导的垂尾抖振气动弹性分析

通过试验方法分析了三角翼前缘分离涡与垂尾抖振之间的关系,深入研究了尾迹流动对垂尾抖振各阶模态的激励作用。计算得到了垂尾模型固有频率及各阶模态。在风洞试验中,应用激光片光烟流场显示技术,得到了三角翼模型在风速为30 m/s下,各迎角的涡结构;使用加速度传感器测量了垂尾翼根和翼梢的抖振响应;使用热线风速仪测量了垂尾翼根和翼梢位置的脉动速度分量。结果表明:前缘涡破裂后产生的高湍流度的尾迹是垂尾抖振的直接原因,抖振边界与涡破裂的强度和位置有关;涡破裂后尾迹与垂尾产生共振,使得抖振加速度响应频率与垂尾固有频率一致;涡破裂后,在较小迎角下,尾迹对垂尾的高频振动模态的激励较为明显,在较大迎角下,涡破裂流动对垂尾低频振动模态的激励加强了。     

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

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

Unsteady flow phenomena associated with leading-edge vortices

This paper presents selected results from extensive experimental investigations on turbulent flow fields and unsteady surface pressures caused by leading-edge vortices, in particular, for vortex breakdown flow. Such turbulent flows may cause severe dynamic aeroelastic problems like wing and/or fin buffeting on fighter-type aircraft. The wind tunnel models used include a generic delta wing as well as a detailed aircraft configuration of canard-delta wing type. The turbulent flow structures are analyzed by root-mean-square and spectral distributions of velocity and pressure fluctuations. Downstream of bursting local maxima of velocity fluctuations occur in a limited radial range around the vortex center. The corresponding spectra exhibit significant peaks indicating that turbulent kinetic energy is channeled into a narrow band. These quasi-periodic velocity oscillations arise from a helical mode instability of the breakdown flow. Due to vortex bursting there is a characteristic increase in surface pressure fluctuations with increasing angle of attack, especially when the burst location moves closer to the apex. The pressure fluctuations also show dominant frequencies corresponding to those of the velocity fluctuations. Using the measured flow field data, scaling parameters are derived for design purposes. It is shown that a frequency parameter based on the local semi-span and the sinus of angle of attack can be used to estimate the frequencies of dynamic loads evoked by vortex bursting.     

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.     

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

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

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。     

Transient aeroelastic responses and flutter analysis of a variable-span wing during the morphing process

To investigate the transient aeroelastic responses and flutter characteristics of a variablespan wing during the morphing process,a novel frst-order state-space aeroelastic model is proposed.The time-varying structural model of the morphing wing is established based on the Euler-Bernoulli beam theory with time-dependent boundary conditions.A nondimensionalization method is used to translate the time-dependent boundary conditions to be time-independent.The time-domain aerodynamic forces are calculated by the reduced-order unsteady vortex lattice method.The morphing parameters,i.e.,wing span length and morphing speed,are of particular interest for understanding the fundamental aeroelastic behavior of variable-span wings.A test case is proposed and numerical results indicate that the flutter characteristics are sensitive to both of the two morphing parameters.It could be noticed that the aeroelastic characteristics during the wing extracting process are more serious than those during the extending process at the same morphing speed by transient aeroelastic response analysis.In addition,a faster morphing process can get better aeroelastic performance while the mechanism comlexity will arise.