间隙非线性对二元翼段颤振特性影响

为了研究间隙非线性对颤振特性的影响,本文采用基于非定常雷诺平均方程的非定常气动力求解方法,耦合结构运动方程建立了时域气动弹性分析系统,并运用该系统计算三自由度无间隙二元翼段构型的颤振速度。采用描述函数法对间隙问题进行处理,得到了间隙非线性所带来的极限环振荡现象,并分析亚跨音速阶段间隙大小对颤振特性的影响。通过研究预加载对颤振特性的影响,得出预加载能够减弱间隙非线性影响,有效提高系统颤振速度。     

采用非定常N—S方程的翼型颤振特性分析研究

以非定常N－S方程为主管方程,计算翼型振动的瞬时非定常气动力,并与颤振方程耦合求解,用时间推进的方法,计算了结构响应特性。经多个算例计算,研究了颤振的临界速度随马赫数的变化规律以及极限环振荡等非线性特性。计算结果与其它献计算结果吻合很好。     

考虑助力器动力学的舵系统结构非线性颤振特性分析

针对操纵系统的结构非线性问题,提出了2种考虑动刚度特性的舵系统(助力器系统和舵面耦合系统)颤振分析方法。利用分支模态法建立了以助力器系统动力学方程为基础的时域和频域非线性控制方程。使用描述函数法和一种改进的迭代方法得到系统控制命令为0时的预载间隙非线性颤振特性,比较了时域和频域结果。结果显示,在非线性条件下,控制命令的幅值和相位对颤振特性都有影响。     

超临界翼型的跨声速颤振特性研究

采用Navier-Stokes方程和二元弯-扭颤振运动方程耦合，用时间推进方法计算结构响应的时问历程，从而得到超临界翼型的跨声速颤振特性。研究了结构参数完全相同的3个超临界翼型(RAE2822，DFVLR-R2和NPU-3)的跨声速颤振特性。为了对比，同时计算了NACAOO系列的3个翼型(NACA0012，NACA0008，NACA0004)的跨声速颤振特性。研究结果表明，翼型形状对颤振特性有明显的影响，相对厚度较小的翼型颤振速度较高；厚度基本相同时，超临界翼型的颤振速度高于NACAOO系列翼型；跨声速范围内，由于气动力的非线性影响，普遍存在极限环振荡。     

基于参数不确定性的机翼颤振风险评定

随着飞机设计的发展,结构非线性颤振问题越来越突出,为了准确评价系统的颤振稳定性,不但要采用非线性颤振计算模型,同时还必须考虑模型中参数的不确定性,面对这些具有随机分布的参数,传统的确定性颤振分析方法已不能胜任;为了探索一套新的,且适合研究非线性和不确定性相互耦合的颤振计算分析方法,文中选取具有非线性扭转刚度的二元机翼颤振系统为研究对象,将概率统计学中的蒙特卡罗方法和核密度估计法结合确定性颤振分析,开展不确定性量化工作,并得到给定速度下系统发生颤振的概率,进而对系统的颤振风险进行评定。并给出算例中非线性扭转刚度表现为软弹簧特性时,不同速度下发生颤振的概率,从而说明系统在不同速度下的颤振风险。     

基于流固耦合方法的大展弦比机翼非线性颤振特性分析

针对大展弦比机翼的柔性大、变形大的特点,基于非定常涡格法求解机翼的非定常气动力,考虑了大展弦比机翼的几何非线性效应,提出了计算大展弦比机翼非线性颤振分析的新方法。以某平板机翼为例:计算了机翼静气动弹性变形、振动特性和颤振特性随着攻角增大的变化规律;比较了颤振结果线性解与非线性解的差别。相关的分析结果表明:大展弦比机翼颤振分析需同时考虑几何非线性效应和气动网格变形。     

几何大变形太阳能无人机非线性气动弹性稳定性研究

大柔性太阳能无人机在气动载荷的作用下产生较大的弯曲变形,机翼结构的刚度、质量分布等特性亦发生较大改变,线性理论无法满足这类飞机气动弹性稳定性分析的精度要求。基于Co-rotational(CR)理论,推导了结构变形后的切线刚度矩阵和质量矩阵,建立了大柔性机翼结构动力学模型;采用建立在局部气流坐标系下的片条非定常气动力模型,建立了考虑几何非线性效应的大柔性无人机气动弹性运动方程。引入准模态假设,采用P-k法研究了几何大变形对类"太阳神"布局太阳能无人机的气动弹性稳定性的影响。研究结果表明:随着弯曲变形的增加,非线性颤振速度可降低10%以上,非线性颤振频率可下降8%;合理的增加扭转刚度、前移弹性轴、前移剖面质心等,均可以有效改善几何大变形引起的不利影响。研究工作对大柔性飞机的气动弹性设计具有一定的参考意义。     

基于欧拉方程的一种机翼气动弹性计算方法

利用一种双时间方法求解三维非定常欧拉方程，采用无限插值理论生成O－H型代数网格，考虑了机翼变形时的网格生成问题，通过与气动力方程的联立求解，在时间域内用二阶龙－库塔方法求解机翼弹性运动方程。计算结果表明，本计算方法具有较高的计算效率，所计算的颤振临界速度与风洞实验一致。     

带间隙约束的二维悬臂亚音速壁板极限环颤振分析

研究了带间隙约束的二维悬臂壁板在亚音速气流作用下的极限环颤振。采用微分求积法对运动方程进行离散,并分析了结构参数对系统响应的影响。结果表明:增加间隙处的弹簧刚度能减小系统极限环颤振幅值;间隙大小不影响系统响应的拓扑结构。     

机翼-机身-尾翼结构的跨音速颤振分析研究

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

超音速、高超音速机翼的气动弹性计算方法

针对超音速和高超音速流动的特点，分析并检验了各种气动力工程算法(牛顿法，切楔／切锥法，活塞理论，激波膨胀波法等)，并将其推广运用于超音速和高超音速机翼的非定常气动力的计算中。通过与机翼结构运动方程的联立求解，在时间域内实现了超音速和高超音速机翼颤振的数值模拟。通过与实验结果的比较，证明该方法具有较高精度，误差能控制在10％左右。     

Navier-Stokes方程预处理方法及其对翼型绕流数值模拟的应用

在与定常方程相容的情况下,通过对N av ier-Stokes方程的时间导数项实施Cho i及Turkel矩阵预处理,发展了适用于低、亚、跨声速粘性流动的有效的预处理方法。对预处理后的控制方程采用了Jam eson格式进行有限体积空间离散和Runge-K utta显示时间推进求解,以及采用了FA S多重网格方法加速收敛。对RAE 2822、GAW-1等不同类型翼型进行了低速和跨声速流动的数值模拟。算例表明:Cho i预处理方法和Turkel预处理方法均有效改善了时间推进方法对低马赫数流动计算的收敛性,而且提高了计算精度;预处理方法同样适用于跨声速非线性流动的计算,且表现出了一定的加速收敛效果。     

基于最大李雅普诺夫指数的壁板热颤振特性分析

文章介绍了李雅普诺夫指数的定义,给出了壁板热颤振特性分析模型,计算得到超声速气流中受热壁板的非线性颤振运动形态以及相应的最大李雅普诺夫指数,画出了壁板稳定性区域图,结果显示不同参数下,超声速气流中受热壁板具有包括颤振在内的5种不同性质的振动形态:收敛于平壁板平衡点的衰减振动、收敛于屈曲平衡点的衰减振动、简谐型极限环颤振、非筒谐周期型极限环颤振和混沌型颤振,并给出了各种振动形态典型的时程响应图、相轨迹以及对应的最大李雅普诺夫指数.分析结果表明,来流速压和温升都严重影响着受扰壁板的非线性振动响应特性,壁板颤振容易发生在来流速压较大或者温度较高的参数区域,系统的最大李雅普诺夫指数能够定量表征壁板稳态响应时的运动形态,可以用该指数来判断超音速气流中受热壁板的稳定性并确定所研究参数范围内壁板的稳定性边界.     

高阶有限体积法研究

在有限体积法的基础上,提出了几种基于Caresian坐标系下的空间显式、隐式高阶离散格式和高阶解重构算法。采用高阶紧致滤波替代二阶、四阶人工耗散抑制Euler方程计算中的非线性高频振荡现象。以Sculley涡在二维均匀无粘流场中的非定常传输为算例,分析了不同格式的耗散误差和相误差。     

基于预处理方法的流动数值分析

以预处理方法为基础,采用了有限体积法,并将高精度高分辨率迎风格式与隐式时间推进法LU－SGS格式结合,对低速到超音速范围内的流动进行数值分析。为了提高AUSMPW^＋格式的计算精度,采用了三阶MUSCL格式,为了验证本方法对低速到超音速范围内的无粘和粘性流动的数值分析效果,对几个典型算例进行了数值试验。本的数值结果与献计算结果和实验数据基本一致,表明本方法可以有效地对低速到超音速范围内的流动进行数值分析。     

基于三阶矩的非最小相位MA模型参数估计

提出一种新的非高斯信号激励的非最小相位ＭＡ模型参数辨识方法。该方法基于三阶矩和自适应算法，并利用ＦＫＡ算法极小化输出信号的三阶矩估计方差。计算机模拟结果表明，木文方法的收敛速度比最大梯度下降法（ＮＰＧ）快，与ＮＰＯ算法同数量级，但每次递归的运算量比ＮＰＯ算法少，占内存量也小。     

T型橡胶减振器非线性特性分析与预压缩量设计

为研究T型橡胶减振器的结构型式和材料性能对其非线性动力学特性的影响,并提高其对飞行环境的适应能力,改善其在恶劣环境下的使用性能,提出了一种含分段线性刚度及阻尼和立方刚度的双层级非线性动力学模型。首先,对T型橡胶减振器进行了简化建模,并采用等价线性化方法,推导得到了分段线性系统的等效刚度和等效阻尼。其次,利用谐波平衡法分别数值模拟了等幅值外力载荷正弦扫频激励下分段线性刚度系统和立方刚度系统的非线性频响函数,验证了构建模型的有效性。最后,进行了T型橡胶减振系统在基础位移正弦扫频条件下的传递特性试验,并与数值计算结果进行了比对。研究结果表明,针对T型橡胶减振器,较低激励载荷即可导致结构型式的渐进软化,较高激励载荷才能引发材料性能的动态软化;该模型能够较好的模拟T型橡胶减振器具有的渐软刚度和渐小阻尼的非线性特性。此外,从航天工程应用的角度出发提出了针对大过载、强振动环境进行预压缩量设计的方法。     

整体式液体冲压发动机级间分离仿真研究

应用数值方法对整体式液体冲压发动机进气道整流罩打开及助推器分离过程进行了研究.动网格采用结构网格中的Chimera重叠网格技术,空间离散格式采用Roe格式,控制方程采用雷诺平均N-S方程,湍流模型为标准k-ε模型.计算结果发现,进气道整流罩打开后高速气流开始充填发动机内通道,发动机头部脱体激波迅速向唇口方向移动,进气道头部斜激波系建立,气流以当地声速向发动机出口方向移动.高压气流到达助推器头部时助推器开始分离,整个分离过程燃烧室压力呈大幅振荡,分离时间为0.0722 s.由于分离时间很短,冲压发动机可以在助推器推出后再点火起动.