静压气体轴承气膜力及其与转子耦合动力学特性研究

针对静压气体轴承,进行三维实体建模,采用供气孔区域非结构化网格和非供气孔区域结构化网格相结合的网格划分方法;采用基于有限体积法的计算流体动力学(CFD)商业软件对三维稳态可压缩气体Navier-Stokes(N-S)方程进行求解;根据计算结果,通过数据拟合获得了考虑转子偏心和转速的静压气体轴承气膜支承力模型.基于有限元法建立了气体轴承-转子系统动力学模型,采用Newmark逐步积分法求解了系统的临界转速和不平衡响应.在此基础上进行实验测试,验证了数值仿真结果.研究结果表明:低速、小偏心下,气膜主支承力随偏心呈近似线性变化;高速、大偏心下,气膜主支承力急剧增大,气体轴承的动压效应显著增强;气膜x,y向耦合支承力随转速和偏心呈非线性变化;转子系统一、二阶临界转速对当前结构刚度变化比较敏感,而三阶临界转速对此不敏感.因此,气体轴承气膜支承力的非线性特性及其与转子耦合动力学特性较为复杂,在对气体轴承进行结构设计时,应充分考虑其与转子的耦合,合理设计工作转速范围.      

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

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

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

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

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

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

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

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

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

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

新型一体化热防护系统热力分析与试验研究

高速飞行器对结构效率的苛刻要求使得热防护系统不断趋于向轻质化、集成化方向发展,新型的力热耦合一体化热防护系统(ITPS)极具发展潜力.首先阐释了一种新型一体化热防护方案的概念与特点,总结了一体化结构设计的基本原则,数值分析了结构参数对背面温度响应、屈曲临界载荷的影响,结果表明腹板厚度对背面温度以及屈曲临界载荷的影响最大.然后设计并加工制备了ITPS的面板与单胞试验样件,分别展开了800℃的高温防隔热性能试验考核和屈曲性能的力学试验研究;试验表明腹板结构是引发热短路效应和屈曲的关键因素,屈曲试验与模拟结果吻合,高温屈曲分析表明温度梯度对屈曲特征有较大影响.     

高功率微波极化方式与入射方向对微带线耦合特性影响分析

针对高功率微波（HPM）耦合效应分析需求,采用时域多分辨分析法（MRTD）模型,给出了微带线端口电压和电流方程,验证了模型的适用性。基于MRTD,仿真分析了HPM作为入射波时其极化方式、入射方向对微带线电磁耦合特性的影响。结果显示：垂直线极化波和平行线极化波耦合系数峰值出现在不同的频率,且峰值相差约15 dB,采用圆极化波时耦合系数最大;入射方向变化时,耦合系数以水平入射最小,垂直入射时最大,高频段微带线边缘有绕射现象产生。该分析采用了分析计算电大尺寸HPM耦合效应的MRTD数学模型,提高了HPM与复杂结构的耦合计算效率。     

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

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

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

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