幂律流体液膜破裂的线性稳定性分析

为对凝胶推进剂撞击式喷嘴产生的液膜的破裂行为进行研究,使用线性稳定性分析方法并结合幂定律本构方程对幂律流体液膜的稳定性进行了研究,并对流体物性参数、流变参数等对液膜流动稳定性的影响规律进行了分析.结果发现:稠度系数、流动指数和表面张力越大,液膜表面扰动波的增长率越小,波长越大.而较高的液膜速度、气液密度比和液膜厚度都使扰动波增长率变大,但波长随液膜速度和气液密度比的增大而减小,而液膜厚度的情况却相反.理论计算与实验值的对比表明,线性稳定性分析能较准确地预测液膜的破碎长度,但预测的表面波长却大于实验值.      

基于泵式流体回路的小卫星姿控/热控一体化执行机构设计

 主要研究了小卫星姿控/热控一体化执行机构的设计问题.首先,根据流体回路中液体流速变化对小卫星产生力矩实现姿态控制、液体流动吸/散废热实现热控制的原理,提出一种姿控/热控一体化执行机构设计方案.然后针对该设计方案,利用以电机转速为变量的流体回路内压强和电磁力矩方程,推导了一体化执行机构姿控力矩模型;利用散热量随流体回路流速的变化,建立了一体化执行机构热控模型.最后,针对某小卫星设计了基于姿控/热控一体化执行机构的闭环控制系统,并针对该一体化执行机构设计了一种姿控/热控解耦算法,对其姿控/热控能力进行数学仿真验证,仿真结果证明了该一体化执行机构的有效性.     

低声爆静音锥设计方法研究

 声爆精确预测及低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。基于计算流体力学(CFD)、波形参数法以及MARK-VII方法构建了高精度声爆预测方法,利用该方法对低声爆静音锥的设计展开研究。研究结果表明,静音锥的设计存在临界长度,静音锥长度小于临界长度时,静音锥产生的激波在传播过程中与机头弓形激波合并,静音锥无法起到降低声爆的作用;静音锥长度大于临界长度时,声爆水平也会略有上升。静音锥临界长度随飞行高度和飞行马赫数的变化而变化,可以根据实际飞行状态采用可伸缩设计,达到最佳的低声爆效果。多级静音锥利用多道弱激波取代机头强弓形激波,其声爆水平较单级静音锥也更低,同样,多级静音锥每一级的长度也要达到临界长度。不同静音锥头部形状产生的脱体激波形状不同,脱体距离也不同,导致阻力系数以及静音锥壁面温度有所不同,但静音锥头部形状对远场声爆信号的影响并不明显。采用静音锥的低声爆方案与原始方案比较,声爆水平得到大幅降低,阻力系数略有上升。     

高超声速飞行器飞行动力学特性不确定分析

 高超声速飞行器分析模型存在较大不确定性,其给出的动力学特性与真实值之间存在偏差,因此研究飞行动力学特性分析结果的可信度水平对控制系统设计具有重要意义。针对典型的乘波体构型高超声速飞行器,在建立气动/结构/推进相互耦合的动力学模型基础上,利用非概率区间来描述模型中的不确定参数,并将复特征根不确定范围求解问题转化为频率和阻尼比两个实数的不确定区间分别求解,给出了动力学模态特征根、频率及阻尼比的不确定边界。分别采用直接蒙特卡罗(DMC)模拟方法、基于泰勒展开的区间分析方法(TIAM)和基于多项式逼近的区间分析方法(CIAM)对高超声速飞行器飞行动力学不确定性进行了研究。结果表明:CIAM计算时间适中,且给出的边界更为准确、安全,适合在控制系统设计和验证过程中使用。     

几何非线性机翼本征梁元素模型的高效化改进

 采用Hodges等提出的时间-空间离散化的几何精确非线性本征梁通用模型处理柔性机翼结构动力学问题时,当离散化的节点数增大时,该方法的未知数数量成倍地增长,而且方程组是严重病态的,因此数值模拟计算的速度非常缓慢。针对机翼中最常见的悬臂梁结构,根据空间离散化的边界条件,提出了空间缩聚法把空间离散差分方程缩聚为常系数矩阵格式,得到了只与时间相关的微分方程组,进一步推导得到了该方程组的雅可比矩阵,因而大大减少了方程组的数量以及求解过程的循环和迭代步数。采用Gear方法分别求解了原始的本征梁元素模型和本文提出的缩聚模型,结果表明空间缩聚模型在相同条件下可提高运算速度约5.1倍,而且对不同类型的外载荷都具有较好的通用性、稳定性和高效性。     

弹性飞行器飞行动力学建模研究

随着飞行器结构模态频率与刚性模态频率愈加接近,弹性效应对弹性飞行器飞行动力学特性的影响变得愈加明显,特别是操稳特性将变得更加的复杂和严峻,已经不能用“刚性飞行器”的分析方法进行研究,因此迫切需要为弹性飞行器建立能够包含多学科耦合的飞行动力学模型.本文对弹性飞行器飞行动力学建模的相关研究进行了总结与发展:首先简要阐述了建立多学科耦合弹性飞行器飞行动力学模型的必要性;然后对传统的弹性飞行器飞行动力学模型进行了研究,分析总结了这些方法的优缺点;最后在此基础上提出了一种新的体轴系(瞬态坐标系),并利用拉格朗日方程和有限元思想推导了该坐标系下的动力学模型,该模型克服了传统模型的缺点,并准确自然地耦合了结构动力学、飞行动力学、空气动力学与控制等学科,较现有的模型而言,该模型能更充分更全面的描述弹性飞行器飞行过程中流场、结构、控制和飞行力学之间的交叉耦合特性.本文的研究成果可为弹性飞行器的动态特性分析提供必要的理论基础.     

基于NURBS曲线的涡控蛇形进气道设计

利用NURBS(non-uniform rational B-spline)曲线成功实现了涡控蛇形进气道参数化描述,并运用数值仿真方法对其中两个关键设计参数进行参数化研究.仿真结果表明:①第二S弯上壁面两侧后掠状凸起型面诱导的受控旋涡能够将低能流牵引至出口两侧,从而抑制大范围的气流分离,但凸起角取值需权衡选取,否则将不利于涡控蛇形进气道综合性能的改善.②通过抬高第二S弯下壁面能够减缓上壁面沿程逆压力梯度,进而影响第二S弯上壁面的流态,恰当的取值能够以微小的总压损失换取大幅度的畸变改善.③当设计参数选取恰当时,涡控蛇形进气道在设计状态下总压恢复系数为0.9667,畸变指数为0.2451.进气道性能较传统方案有显著改善,使得蛇形进气道迈向工程实用成为可能.      

剪刀式尾桨悬停状态气动力及噪声特性计算研究

建立了一个基于计算流体力学(CFD) /FW-H(Ffowcs Williams-Hawkings)方程的预测剪刀式尾桨悬停状态气动性能和噪声特性的分析方法.该方法首先采用CFD方法对尾桨流场进行求解,并应用嵌套网格技术对流场空间进行离散.控制方程采用非惯性坐标系下的Navier-Stokes方程,空间方向采用二阶迎风格式(Roe格式)进行求解,时间方向采用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式进行推进.在此基础上,采用FW-H方程将尾桨噪声声压扰动传播至远场,以获得尾桨的噪声特性.应用该方法对两种剪刀式尾桨构型(“L”构型和“U”构型)进行了计算研究,对比分析了剪刀式尾桨在气动力和噪声方面与常规尾桨的差别,以及两个重要构型参数(剪刀角和轴向间距)对剪刀式尾桨气动力和噪声的影响规律.计算结果表明,构型参数对剪刀式尾桨气动力和噪声特性影响很大,合理地选择构型参数可以降低尾桨噪声水平.     

高速车辆侧窗风振噪声仿真分析与控制

采用大涡模拟的计算方法,对某轿车在不同侧窗开启工况下风振噪声特性进行了数值仿真.揭示了侧窗风振噪声的产生机理,分析了不同侧窗开启方式下驾驶员耳旁声压级存在差异的原因.分析结果表明:后窗附近湍流强度弱,能量耗散少,开启时风振噪声最大;开启2个侧窗时,由于气流导出效应,可以有效地降低风振噪声.提出了使用扰流装置抑制风振噪声的措施,探讨了不同形状扰流装置对风振噪声影响.仿真结果表明:方形扰流装置的降噪效果最明显,最大降幅可达15dB.      

基于解析解的人字齿轮齿廓修形动态分析

根据人字齿轮传动的结构,考虑时变啮合刚度、齿侧间隙,并把齿廓修形作为一种时变齿侧间隙计入,建立了单级人字齿轮传动的弯-扭-轴耦合动力学模型和相应的非线性动力学方程.对方程的数值解进行分析,发现合理的修形使人字齿轮传动不发生齿面的完全分离.基于多尺度法对不发生齿面完全分离时的人字齿轮传动系统进行了摄动分析,得到了系统在主共振频率附近幅频响应的近似解析解.相比数值方法,解析法具有很高的求解效率.最后,根据近似解析解研究了修形参数对动载系数幅值的影响.研究发现:螺旋角越大使动载系数最小所需的修形量越大,而所需的修形长度越小.另外,使动载系数最小的最优修形参数解域呈新月形,其随着主动轮支承刚度的增大而先增大后减小;随着被动轮支承刚度的增大而持续减小.